g.m

来自「信号与系统课程设计全部内容」· M 代码 · 共 12 行

num=[1]
den=[2 4 4 2]
sys=tf(num,den)
poles=roots(den)
figure(1);pzmap(sys)
t=0:0.02:10
h=impulse(num,den,t)
figure(2);plot(t,h)
[H,w]=freqs(num,den)
figure(3);plot(w,abs(H))