m0714.m

来自「matlab的经典讲义」· M 代码 · 共 14 行

N=128;
x=[1 zeros(1,N-1)];                 %产生单位冲激函数
num=[0.01 -0.03 0.06 -0.03 0.01];
den=[1 2.5 3 1.8 0.5];
y=filter(num,den,x);                %计算单位冲激响应
figure(1);
n=1:N;
stem(n,y,'fill');grid on; 
title('单位冲激响应');
figure(2)
Fs=1024;
freqz(num,den,N,Fs);           %绘制幅频和相频响应曲线
grid on;