samp5_14.m

MATLAB7.x数字信号处理 光盘内容

%Samp5_14
wp=[0.35 0.65]*pi;   %给出滤波器的通带边界频率，用弧度表示
ws=[0.2 0.8]*pi;   %给出滤波器的阻带边界频率，，用弧度表示
Rp=1;Rs=60;   %给出滤波器的通带波纹和阻带衰减
[n,wn]=cheb1ord(wp,ws,Rp,Rs,'s');  %根据滤波器的要求求滤波器的最小阶数
[z,p,k]=cheb1ap(n,Rp);  %用求得的滤波器的最小阶数设计Chebyshev I型滤波器
[b,a]=zp2tf(z,p,k); %将零点极点增益形式转换为传递函数形式
Wo=sqrt(wn(1)*wn(2)); %计算中心点频率
Bw=wn(2)-wn(1);  %计算频带宽度
[bt,at]=lp2bp(b,a,Wo,Bw);  %模拟原型低通滤波器转换为带通滤波器
[h,w]=freqs(bt,at,128);  %计算滤波器的复数频率响应
plot(w/pi,20*log10(abs(h)));   %绘制滤波器的振幅响应
xlabel('角频率/pi');
ylabel('|H(jw)|/dB')
grid on;hold on;
plot([0.2 0.2],ylim, 'r'); plot([0.8 0.8],ylim, 'r');%绘出阻带界限以显示达到的阻带衰减
%plot([0.35 0.35],ylim); plot([0.65 0.65],ylim);%可绘出通带界限以显示达到的通带衰减