ex31.m

matlab exprience,M文件,实验代码

%ex3(1)
N=100;
implse=[1 zeros(1,99)];
heaviside=ones(1,N);
t=1:100;
b=[1,-3,11,27,18];
a=[1,16,12,2,-4,-1];

y=impz(b,a);
y1=filter(b,a,heaviside);
subplot(411)
stem(t,implse);
title('implse')
subplot(412)
stem(heaviside)
title('u(t)')
subplot(413)
plot(y,'k-');
title('implse responce')
grid on;
ylabel('y(n)')
xlabel('n')
subplot(414)
plot(y1)
title('u(t) response')