ms18r22.m

《数字信号处理实验(MATLAB版)》这

%切比雪夫Ⅰ型模拟滤波器
fp=6000;Omgp=2*pi*fp;    %输入滤波器的通带截止频率
fs=15000;Omgs=2*pi*fs;    %输入滤波器的阻带截止频率
Rp=1;As=30;              %输入滤波器的通阻带衰减指标
%计算滤波器的阶数和通带截止频率
[n,Omgn]=cheb1ord(Omgp,Omgs,Rp,As,'s')
%计算n阶模拟低通原型，得到左半平面零极点
[z0,p0,k0]=cheb1ap(n,Rp)
b0=k0*real(poly(z0))     %求滤波器系数b0
a0=real(poly(p0))        %求滤波器系数a0
[H,Omg]=freqs(b0,a0);    %求系统的频率特性
dbH=20*log10(abs(H)/max(abs(H))); %化为分贝值

subplot(2,2,1),plot(Omg*Omgc/(2*pi),abs(H)),grid 
axis([0,20000,0,1.1]);
ylabel('幅度');xlabel('f(Hz)');
subplot(2,2,2),plot(Omg*Omgc/(2*pi),angle(H)),grid
axis([0,20000,-4,4]);
ylabel('相位');xlabel('f(Hz)');

subplot(2,2,3),plot(Omg*Omgc/(2*pi),dbH),grid 
axis([0,20000,-50,2]);
ylabel('幅度(dB');xlabel('f(Hz)');
subplot(2,2,4),pzmap(b0,a0);
axis square,axis equal,grid on

Omgx0=[Omgp,Omgs]/Omgc; %设置频率向量
Hx=freqs(b0,a0,Omgx0);  %计算该两点的频率特性
dbHx=-20*log10(abs(Hx)/max(abs(H))) %化为分贝值