华东理工大学MATLAB分析.ppt

b=[1 1/0.9];a=[1 0.9]; [Fh,w] = freqz(b,a); [Gd,w] = grpdelay(b,a); subplot(311) plot(w/pi,abs(Fh));ylabel(|H(w)|); grid on; axis([0 max(w/pi) 0 1.5]) subplot(312) plot(w/pi,angle(Fh)); ylabel(ang[H(w)]); grid on; subplot(313) plot(w/pi,Gd);ylabel(grd[H(w)]); grid on; b=[1 0 -1/0.81];a=[1 0 -0.81]; [Fh,w] = freqz(b,a); [Gd,w] = grpdelay(b,a); subplot(311) plot(w/pi,abs(Fh));ylabel(|H(w)|); grid on; axis([0 max(w/pi) 0 1.5]) subplot(312) plot(w/pi,angle(Fh)); ylabel(ang[H(w)]); grid on; subplot(313) plot(w/pi,Gd);ylabel(grd[H(w)]); grid on; 6、最小相位系统