untitled3.m

来自「对数字图象进行分步处理,第一步滤波, 第二步去噪,第三步显示」· M 代码 · 共 11 行

%2.11(1)
b=[0,0.2];a=[1,0,0.8];n=1:100;
[H,w]=freqz(b,a,256,'whole',1);
Hr=abs(H);
Hphase=angle(H);
Hphase=unwrap(Hphase);
subplot(211)
plot(w,Hr);
subplot(212);
plot(w,Hphase);