eg10_7.m

来自「matlab7.0从入门到精通(光盘内容)」· M 代码 · 共 29 行

%由系统的差分方程求出系统的冲激响应，画出其幅度、相位响应，并指出系统的零极点
A=[1,-1,0.6];
B=[1,-2];
%指出系统的零极点
[R,p,C]=residuez(B,A)
%求系统的冲激响应
n=[-10:50];
x=[(n)==0];
h=filter(B,A,x);
figure(1);
plot(n,h);
xlabel('n');
ylabel('h(n)');
title('系统的冲激响应');
%求其幅度相位响应
[H,W]=freqz(B,A,100);
magH=abs(H);
phaH=angle(H);
figure(2);
plot(W,magH);
xlabel('频率');
ylabel('幅度');
title('幅度响应');
figure(3);
plot(W,phaH);
xlabel('频率');
ylabel('幅度');
title('相位响应');