数字信号处理上机实验滤波器讲义.doc

数字信号处理上机练习 7.12 （1）双线性变换法 程序： %双线性变换法设计切比雪夫1型带通滤波器 %w计算式：w=(f/ADfs)*2*pi 单位：rad clc;clear;close all; ADfs=2000;%ADC的采样频率，单位：Hz ws1=(100/ADfs)*2*pi;%下阻带截止频率，单位：rad wp1=(200/ADfs)*2*pi;%通带下截止频率，单位：rad wp2=(400/ADfs)*2*pi;%通带上截止频率，单位：rad ws2=(600/ADfs)*2*pi;%上阻带截止频率，单位：rad display(数字滤波器的边界频率 单位：rad);ws1, wp1,wp2,ws2 Rp=2;%通带最大衰减，单位：dB As=20;%阻带最小衰减，单位：dB T=1;%模拟滤波器数字化时的采样间隔，单位：s Ws1=(2/T)*tan(ws1/2);Wp1=(2/T)*tan(wp1/2); Wp2=(2/T)*tan(wp2/2);Ws2=(2/T)*tan(ws2/2);%转换为模拟滤波器指标 Ws=[Ws1,Ws2];Wp=[Wp1,Wp2]; display(转换后滤波器的边界频率 单位：rad/s);Ws1, Wp1,Wp2,Ws2 [n,wn] = cheb1ord(Wp,Ws,Rp,As,s);%计算切比