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

开始

↓

读入数字滤波器技术指标

↓

将指标转换成归一化模拟低通滤波器的指标

↓

设计归一化的模拟低通滤波器阶数N和截止频率

↓

模拟域频率变换,将H(P变换成模拟带通滤波器H(s

↓

用双线性变换法将H(s转换成数字带通滤波器H(z

↓

输入信号后显示相关结果

求相应的幅频响应与相频响应

↓

50

100

150

-202

t

x1(t

x1的波形

50

100

150

-202

t

x2(t

x2的波形

50

100

150

-2

02

tx(t

输入信号x的波形

-0.01

-0.005

00.005

0.01

t

y

滤波器输出y的波形

clc;

clearall;

结束

%数字滤波器的技术指标

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([01