Fe = 22050;
f = 440; T = 1; dt = 1/Fe; t = 0:dt:T;
s = sin(2*pi*f*t);
sound(s,Fe);
disp('son pur');
disp('appuyez sur une touche pour continuer');
disp(' ');
pause
%-------------------------------------------------
a = 1; p = 0;
s = synthad(a,f,p,T,Fe);
sound(s,Fe);
disp('le meme via synthad');
disp('appuyez sur une touche pour continuer');
disp(' ');
pause
%-------------------------------------------------
a = [100 60 20 10 5 1];
f = 440*[1 2 3 4 5 6];
p = 0*a;
s = synthad(a,f,p,T,Fe);
sound(s,Fe);
disp('harmoniques 1 a 6');
disp('appuyez sur une touche pour continuer');
disp(' ');
pause
%-------------------------------------------------
a = [100 60 20 10 5 1];
f = 440*[1 3 5 7 9 11];
p = 0*a;
s = synthad(a,f,p,T,Fe);
sound(s,Fe);
disp('harmoniques impairs');
disp('appuyez sur une touche pour continuer');
disp(' ');
pause
%-------------------------------------------------
m = 10; M = 1:m;
a = 1./M; f = 440*M; p = 0*M;
s = synthad(a,f,p,T,Fe);
sound(s,Fe);
disp('son plus riche');
disp('appuyez sur une touche pour continuer');
disp(' ');
pause
%-------------------------------------------------
a = [350,950,500,150,700,100,250,370,1000,180,300,100,...
150,300,100,100,50,20,10,35,5,15];
f = 220*[0.5,1,1.188,1.530,2.000,2.470,2.607,2.650,2.991,3.367,...
4.137,4.487,4.829,5.385,5.863,6.709,8.077,8.547,...
9.017,9.530,11.026,12.393];
p = 0*a;
s = synthad(a,f,p,T,Fe);
sound(s,Fe);
disp('partiels de la cloche');
disp('appuyez sur une touche pour continuer');
disp(' ');
pause
%-------------------------------------------------
s = bruit(T,Fe);
sound(s,Fe);
disp('bruit blanc');
disp('appuyez sur une touche pour continuer');
disp(' ');
pause
%-------------------------------------------------
s = sirene(800,2,200,3,Fe);
sound(s,Fe);
disp('sirene');
disp('appuyez sur une touche pour continuer');
disp(' ');
pause
%-------------------------------------------------
a = [10, 5, 1];
f = [440,880,1320];
p = [0, 0, 0]+pi;
T = 1; Fe = 11025;
s1 = synthad(a,f,p,T,Fe);
s2 = synthad(a,1.5*f,p,T,Fe);
s3 = synthad(a,2*f,p,T,Fe);
s = [s1,s2,s3];
sound(s,Fe);
disp('trois notes consecutives');
disp('appuyez sur une touche pour continuer');
disp(' ');
pause
plot(s);
disp('analyse temporelle');
disp('cliquez dans la figure sur + pour zoomer');
disp('et delimitez avec la souris la zone a agrandir');
disp('indication valable aussi pour les figures suivantes');
disp('appuyez sur une touche pour continuer');
disp(' ');
pause

z = fft(s);
T = length(s)/Fe;
fr = 0:1/T:4000;
nf = length(fr);
plot(fr,abs(z(1:nf)));
disp('analyse frequentielle');
disp('appuyez sur une touche pour continuer');
disp(' ');
pause

specgram(s,512,Fe);
disp('analyse temps frequence');
disp('appuyez sur une touche pour continuer');
disp(' ');
pause
%-------------------------------------------------
%[s,Fe] = wavread('siffleur');
[s,Fe] = wavread('steinwayE');
sound(s,Fe);
plot(s);
disp('son et analyse temporelle d''une note de piano');
disp('appuyez sur une touche pour continuer');
disp(' ');
pause

z = fft(s);
T = length(s)/Fe;
fr = 0:1/T:4000;
nf = length(fr);
plot(fr,abs(z(1:nf)));
disp('analyse frequentielle de cette note');
disp('appuyez sur une touche pour continuer');
disp(' ');
pause

specgram(s,8*1024,Fe);
axis([0 T 0 3000]);
disp('analyse temps frequence de cette note');
disp('appuyez sur une touche pour continuer');
disp(' ');
pause

%-------------------------------------------------
b = fir1(100,.006);
s1 = filter(b,1,s);
soundsc(s1,Fe);
disp('filtrage passe-bas de cette note');
disp('appuyez sur une touche pour continuer');
disp(' ');
pause

disp('filtrage passe-haut de cette note');
b = fir1(100,.05,'high');
s1 = filter(b,1,s);
soundsc(s1,Fe);