绘制Chebshev 1型模拟低通滤波器的平方幅频响应曲线,阶数为2、4、6、8:
n=0:0.01:2;
for i=1:4
switch i
case 1 ,N=2;
case 2; N=4;
case 3;N=6;
case 4; N=8;
end
Rp = 1; %设置通滤波纹为1db
[z,p,k] =cheb1ap(N,Rp); %设计chebyshev 1 滤波器
[b,a] = zp2tf(z,p,k); %将零点、级点、增益形式转换为传递函数的形式
[H,w] =freqs(b,a,n); %按n指定的频率点给出频率响应
magH2 = (abs(H)).^2; %给出传递函数的幅度平方
posplot=['2,2',num2str(i)]; %将数字i转换为字符转,与‘2,2’合并并赋给posplot
subplot(posplot);
plot(w,magH2);
xlabel ('w/wc');
ylabel('Chebyshev 1 | H(jw)|^2');
grid on;
end