1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27
| clear all ;
f_min=input('valeur de la frequence minimun : ')
f_max=input('valeur de la frequence maximun : ')
p=input('nombre d echantillon : ')
tau=input('valeur du retard : ')
for j=1:p
f=f_min:1/1000:f_max ;
s(j)=sin(2*pi*f*(2*j+1))./(2*j+1);
end
h=s(j)*4/pi;
M=max(h)
subplot(1,2,1);
plot(f,10*log10(h/M))
title('module du spectre ')
xlabel('fréquence/GHz')
ylabel('module')
subplot(1,2,2);
plot(1:p,h)
title('réponse impulsionnelle bruitée')
xlabel('numéro du coefficient ')
ylabel('valeur du coefficient') |