batpq.m

我做了很久的

function batpq(n,wc) %绘制巴特沃斯滤波器零极点图，n:阶数，wc:截止频率
a=[1./((j*wc)^(2*n)) zeros(1,2*n-1) 1]; %定义系统函数分母多项式系数
b=[1];
p=roots(a)  
ljd(a,b)
u=0:pi/200:2*pi;
r=wc*exp(j*u);
plot(r,':')%作以wc半径的单位圆
set(gcf,'color',[1 1 1])
set(gca,'box','on')
title('巴特沃斯滤波器零极点分布图')
xlabel('s平面实轴')
ylabel('s平面虚轴')