batpq.m

来自「信号与系统分析及MATLAB实现一书中的所有源代码」· M 代码 · 共 18 行

function batpq(n,wc)
%绘制巴特沃兹滤波器零极点图程序
%n：巴特沃兹滤波器的阶数
%wc：巴特沃兹滤波器的截止频率
hold off
a=[1./((i*wc)^(2*n)) zeros(1,2*n-1) 1];	%定义系统函数分母多项式系数向量
b=[1];												%定义系统函数分子多项式系数向量
p=roots(a)
sjdt(a,b)											%绘制系统零极点分布图
u=0:pi/200:2*pi;
r=wc*exp(i*u);
plot(r,':')											%画半径为wc的圆
set(gcf,'color',[1 1 1])
set(gca,'box','on')
title('巴特沃兹滤波器极点分布图')
xlabel('S平面实轴')
ylabel('S平面虚轴')