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 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
| close all;
clear all;
K=[0.9 0.2 1 -0.9 1.1];
N=2048;
n=1:N;
g=zeros(5,N);
g(:,1)=1;
for p=1:5
for n=2:2048
g(p,n)=K(p).*g(p,n-1);
end
end
n=1:N;
subplot(5,1,1)
stem(n,g(1,:));
title(['Reponse impulsionnelle pour K=',num2str(K(1))]);
ylabel('amplitude');
xlabel('n');
axis([0 10]);
grid on;
subplot(5,1,2)
stem(n,g(2,:));
title(['Reponse impulsionnelle pour K=',num2str(K(2))]);
ylabel('amplitude');
xlabel('n');
axis([0 10]);
grid on;
subplot(5,1,3)
stem(n,g(3,:));
title(['Reponse impulsionnelle pour K=',num2str(K(3))]);
ylabel('amplitude');
xlabel('n');
axis([0 10]);
grid on;
subplot(5,1,4)
stem(n,g(4,:));
title(['Reponse impulsionnelle pour K=',num2str(K(4))]);
ylabel('amplitude');
xlabel('n');
axis([0 10]);
grid on;
subplot(5,1,5)
stem(n,g(5,:));
title(['Reponse impulsionnelle pour K=',num2str(K(5))]);
ylabel('amplitude');
xlabel('n');
axis([0 10]);
grid on;
print('QP1.png');
pause; |
Partager