q616.m

来自「这是matlab在电子信息课程中的应用」· M 代码 · 共 20 行

%《MATLAB在电子信息课程中的应用》第六章例6.16程序q616
% 二阶巴特沃斯低通滤波器的频率响应
% 电子工业出版社出版  陈怀琛 吴大正 高西全合著 2001年10月

% 二阶巴特沃斯低通滤波器的离散系统函数为：
% H(z)=(z^2+2z+1)/((2+sqrt(2))z^2+(2+sqrt(2));
% 求其频率响应可将z用exp(i*w)代入

b=[1,2,1]; a=[2+sqrt(2),0,2-sqrt(2)];	% 给出滤波器分子分母系数
N=input('取频率数组的点数N= ');		
w=[0:N-1]*pi/N;							% 给出0到pi之间的频率数组
H=polyval(b,exp(i*w))./polyval(a,exp(i*w));	% 求频率响应
figure(1)							% 在线性坐标内画频率特性
subplot(2,1,1),plot(w,abs(H)),grid
subplot(2,1,2),plot(w,unwrap(angle(H))),grid
figure(2)							% 在对数坐标内画频率特性
subplot(2,1,1),semilogx(w,20*log10(abs(H))),grid
subplot(2,1,2),semilogx(w,unwrap(angle(H))),grid