exa4_14.m

来自「matlab应用实例。如：FIR、IIR数字滤波器设计等。」· M 代码 · 共 20 行

% exa4-14_zeropole.m , for example 4-14
% using zero pole method to design highpass DF

clear all;
syms rad z hz hrad angle;		
Wh=0.8*pi;
a=cos(Wh)/(1+sin(Wh))
hz=(2/(1-a))*((1-1/z)/(1+a/z));
hrad=subs(hz,exp(i*rad));
subplot(211);
%ezplot(abs(hrad),[-2*pi,2*pi]);
ezplot(abs(hrad));
title('幅度特性曲线')
grid on;
subplot(212);
angle=atan(imag(hrad)/real(hrad));
ezplot(angle);
title('相位特性曲线')
grid on;