t=0:0.001:0.255 y=rand(size(t)) x=2*sin(2*3.14*50*t)+2.5*sin(2*3.14*320*t)+0.5*y subplot(331) plot(t,x) title('E(t)') N=length(t) f=(0:N-1)*(1000/N) X=fft(x) Pxx=X.*conj(X) subplot(332) plot(f,Pxx) title('Sortie frequentiel de sys ') axis([0 500 0 110000]) grid on T=1/(200*pi) h=(1/T)*(exp(-t/T)) h2=(1/T*T).*t.*(exp(-t/T)) subplot(333) plot(t,h) axis([0 0.015 0 700]) H=fft(h) ph=H.*conj(H) subplot(334) plot(t,h) semilogx(f,20*log(ph)) axis([0 500 240 290]) grid on s=conv(x,h) S=fft(s,256) Pxx1=S.*conj(S) subplot(336) plot(f,Pxx1) title('Sortie de sys par application de filtre 1er ordre') subplot(337) plot(t,h2) H2=fft(h2) ph2=H2.*conj(H2) subplot(338) plot(t,h2) semilogx(f,20*log(ph2)) grid on s2=conv(x,h2) S2=fft(s2,256) Pxx12=S2.*conj(S2) subplot(339) plot(f,Pxx12) title('Sortie de sys par application de filtre 2er ordre')