f.m

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

w=linspace(0,5,200)
b=[1]
a=[2 4 4 2]
H=freqs(b,a,w)
subplot(2,1,1)
plot(w,abs(H))
set(gca,'xtick',[0 1 2 3 4 5]) 
set(gca,'ytick',[0 0.4 0.707  1]);grid
subplot(2,1,2)
plot(w,angle(H))
set(gca,'xtick',[0 1 2 3 4 5]);grid