pul=[1 zeros(1,63)];%產生單位取樣序列
b=[0.2 0.1];
a=[1 -0.4 -0.5];%b和a 都是系統參數
h=filter(b,a,pul);
hl=impz(b,a,64);
figure(1) %作圖
stem(h)
title('FILTER function')
figure(2)
stem(hl)
title('IMPZ function')