ex_7.m

来自「信号与系统分析以matlab实现」· M 代码 · 共 12 行

%P289 7
clear all;
A=input('please input A=');
B=input('please input B=');
b=[1 -2*A*cos(B) A^2];
a=[1 -2/A*cos(B) 1/A^2];
[z,p]=tf2zp(b,a);
figure(1)
subplot(2,1,1),zplane(z,p)
subplot(2,1,2),impz(b,a)
figure(2)
freqz(b,a)