Matlab课程设计利用MATLAB结合双线性变换法设计数字切比雪夫带通IIR滤波器代码.docxVIP

开始 ↓ 读入数字滤波器技术指标 ↓ 将指标转换成归一化模拟低通滤波器的指标 ↓ 设计归一化的模拟低通滤波器阶数 N 和截止频率 ↓ 模拟域频率变换，将 H(P)变换成模拟带通滤波器 H(s) ↓ 用双线性变换法将H(s)转换成数字带通滤波器H(z) ↓ 输入信号后显示相关结果 求相应的幅频响应与相频响应 ↓ 结束 x1的波形 x2的波形 2 2 (1t) 0 ( x -2 2t) 0 (x ( -2 )t(x 0 50 100 150 t 输入信号x的波形 2 0 y -2  0.01 0.005 0 -0.005 0 50 100 150 t 滤波器输出y的波形 clc; clear all;  0 50 -0.01 100 150 0 t  10 20 t  30 40 %数字滤波器的技术指标 Rp = 1; % 通带最大衰减Rs = 40;% 阻带最小衰减 OmegaS1_1=350; % 通带截止频率OmegaS1_2=550;% 通带截止频率OmegaP1_1=400; % 阻带截止频率OmegaP1_2=500;% 阻 带 截 止 频 率 Fp=2000; % 抽样频率Wp1=2*pi*OmegaP1_1/Fp; % 模数频率变换 Wp2=2*pi*OmegaP1_2/Fp; Ws1=2*pi*OmegaS1_1/Fp; Ws2=2*pi*OmegaS1_2/Fp; OmegaP1=2*Fp*tan(Wp1/2); % 非线性变换 OmegaP2=2*Fp*tan(Wp2/2); % 非线性变换 OmegaS1=2*Fp*tan(Ws1/2); % 非线性变换OmegaS2=2*Fp*tan(Ws2/2); % 非线性变换OmegaP0=sqrt(OmegaP1*OmegaP2);% 等效中心频率Bw=OmegaP2-OmegaP1; % 带通滤波器的通带宽度 Eta_P0=OmegaP0/Bw; % 归一化处理 Eta_P1=OmegaP1/Bw; % 归一化处理 Eta_P2=OmegaP2/Bw; % 归一化处理 Eta_S1=OmegaS1/Bw; % 归一化处理 Eta_S2=OmegaS2/Bw; % 归一化处理 Lemta_P_EquivalentLowPass=Eta_P2/(Eta_P2^2-Eta_P0^2); % 转换成低通参数 Lemta_S1_EquivalentLowPass=-Eta_S1/(Eta_S1^2-Eta_P0^2); % 转换成低通参数 Lemta_S2_EquivalentLowPass=Eta_S2/(Eta_S2^2-Eta_P0^2); % 转换成低通参数Lemta_S_EquivalentLowPass=min(Lemta_S1_EquivalentLowPass,Lemta_S2 _EquivalentLowPass); % 取最小值 % E求滤波器阶数 [N, Wn]=cheb2ord(Lemta_P_EquivalentLowPass, Lemta_S_EquivalentLowPass, Rp, Rs,s); % 滤波器设计 [num1,den1]=cheby2(N,Rs,Wn,s); [num2,den2]=lp2bp(num1,den1,OmegaP0,Bw); [num,den]=bilinear(num2,den2,Fp); [Z,P,K]=cheb1ap(N,Rp); w=linspace(1,1000,100)*2*pi; [M1,N1]=zp2tf(Z,P,K); %将零极点形式转换为传输函数形式 [M,N]=lp2bp(M1,N1,OmegaP0,Bw); %对低通滤波器进行频率变换转换为带通滤波器 % 计算增益响应 w = 0:pi/255:pi; h = freqz(num,den,w); g = 20*log10(abs(h)); %绘制切比雪夫带通滤波器幅频特性 figure; plot(w/pi,g);grid axis([0 1 -60 5]); xlabel(\频率/\pi); ylabel(增益/dB); title(切比雪夫II型带通滤波器幅频响应); %Plot the poles and zeros [z,p,k]=tf2zp(num,den); figure; zplane(z,p); %绘制传输函数零极点 title(?传输函数的零极点) f1=450;f2=600; t=0:0.0001:1 x1=sin(2*pi*f1*t); x2=sin(2*pi*f2*t); x=x1+x2; figure; subplot(2,2,1)%绘制x1的波形plot(x1);grid on; axis