模拟低通原型滤波器的MATLAB设计

一、实验目的
1．掌握巴特沃兹、切比雪夫Ⅰ、Ⅱ型和椭圆型模拟滤波器的特性和技术指标。
2．掌握用MATLAB 设计巴特沃兹、切比雪夫Ⅰ、Ⅱ型和椭圆型低通滤波器的方法。
二、实验原理
1．巴特沃兹低通数字滤波器幅度平方函数定义为

四、实验报告要求
1.简述实验目的及原理
2.整理好经过运行并证明是正确的实验程序并加上注释。绘出相应的图形。
3.比较各种常用的模拟原型滤波器的优缺点。选取原型滤波器的标准是什么？

op=0.2*pi;
rp=1;
os=0.3*pi;
as=16;
figure;
[b,a]=afd_butt(op,os,rp,as);
[H,w]=freqs(b,a);
h=impulse(b,a);
t=0:length(h)-1;
subplot(221),plot(w/(2*pi),abs(H));title('巴特幅频特性');xlabel('Hz');
subplot(222),plot(w/(2*pi),angle(H));title('巴特相频特性');xlabel('Hz');
subplot(223),plot(w/(2*pi),20*log10(abs(H)/max(abs(H))));title('巴特幅度分贝表示');xlabel('Hz');
subplot(224),plot(t,h);title('巴特单位冲激响应h(t)');xlabel('sec');
figure;
[b,a]=afd_chb1(op,os,rp,as);
[H,w]=freqs(b,a);
h=impulse(b,a);
t=0:length(h)-1;
subplot(221),plot(w/(2*pi),abs(H));title('切比1幅频特性');xlabel('Hz');
subplot(222),plot(w/(2*pi),angle(H));title('切比1相频特性');xlabel('Hz');
subplot(223),plot(w/(2*pi),20*log10(abs(H)/max(abs(H))));title('切比1幅度分贝表示');xlabel('Hz');
subplot(224),plot(t,h);title('切比1单位冲激响应h(t)');xlabel('sec');
figure;
[b,a]=afd_chb2(op,os,rp,as);
[H,w]=freqs(b,a);
h=impulse(b,a);
t=0:length(h)-1;
subplot(221),plot(w/(2*pi),abs(H));title('切比2幅频特性');xlabel('Hz');
subplot(222),plot(w/(2*pi),angle(H));title('切比2相频特性');xlabel('Hz');
subplot(223),plot(w/(2*pi),20*log10(abs(H)/max(abs(H))));title('切比2幅度分贝表示');xlabel('Hz');
subplot(224),plot(t,h);title('切比2单位冲激响应h(t)');xlabel('sec');
figure;
[b,a]=afd_elip(op,os,rp,as);
[H,w]=freqs(b,a);
h=impulse(b,a);
t=0:length(h)-1;
subplot(221),plot(w/(2*pi),abs(H));title('椭圆幅频特性');xlabel('Hz');
subplot(222),plot(w/(2*pi),angle(H));title('椭圆相频特性');xlabel('Hz');
subplot(223),plot(w/(2*pi),20*log10(abs(H)/max(abs(H))));title('椭圆幅度分贝表示');xlabel('Hz');
subplot(224),plot(t,h);title('椭圆单位冲激响应h(t)');xlabel('sec');
调用的函数：
(1)afd_butt.m
function [b,a]=afd_butt(Wp,Ws,Rp,As);
%Analog Lowpass Filter Design :Butterworth
%--------------------------------------------
%[b,a]=afd_butt(Wp,Ws,Rp,As);
%b=Numberator coefficients of Ha(s)
%a=Denominator coefficients of Ha(s)
%Wp=Passband edge frequency in rad/sec;Wp>0
%Ws=Stopband edge frequency in rad/sec;Ws>Wp>0
%Rp=Passband ripple in +dB;(Rp>0)
%As=Stopband attenuation in +dB;(As>0)
if Wp<=0
error('Passband edge must be larger than 0')
end
if Ws<=Wp
error('Stopband edge must be larger than Passband edge')
end
if (Rp<=0)|(As<0)
error('PB ripple and/or SB attenuation must be larger than 0')
end
N=ceil(log10((10^(Rp/10)-1)/(10^(As/10)))/(2*log10(Wp/Ws)));
fprintf('\n***Butterworth Filter Order=%2.0f\n',N)
OmegaC=Wp/((10^(Rp/10)-1)^(1/(2*N)));
[b,a]=u_buttap(N,OmegaC);
(2) afd_chb1.m
% Chebyshev I 型模拟低通滤波器原型设计;
% afd_chb1.m;
function [b a]=afd_chb1(Wp,Ws,Rp,As);
% Anolog Lowpass Filter Design:chebyshev-1
% ^^^^^^^^^^^^^^^^^^^
% [b a]=afd_chb1(Wp,Ws,Rp,As);
% b=numerator polynomial coefficients of Ha(s);
% a=denominator polynomial coefficients of Ha(s);
% Wp=passband edge frequency in rad/sec;Wp>0;
% Ws=stopband edge frequency in rad/sec;Ws>Wp>0;
% Rp=passband ripple in +dB;(Rp>0);
% As=stopband attenuation in + dB;(As>0);
% ^^^^^^^^^^^^^^^^^^^^
if Wp<=0
error('passband edge must be larger than 0')
end
if Ws<=Wp
error('stopband dege must be larger than passband edge')
end
if(Rp<=0)|(As<0)
error('PB ripple and/or SB attenuation must be larger than 0')
end
ep=sqrt(10^(Rp/10)-1);
A=10^(As/20);
OmegaC=Wp;
OmegaR=Ws/Wp;
g=sqrt(A*A-1)/ep;
N=ceil(log10(g+sqrt(g*g-1))/log10(OmegaR+sqrt(OmegaR*OmegaR-1)));
fprintf('\n***Chebyshev-1 filter order=%2.0f\n',N)
[b a]=u_chblap(N,Rp,OmegaC);
(3) afd_chb2.m
% Chebyshev II 型模拟低通滤波器原型设计;
% afd_chb2.m;
function [b a]=afd_chb2(Wp,Ws,Rp,As);
% Anolog Lowpass Filter Design:chebyshev-2
% ^^^^^^^^^^^^^^^^^^^
% [b a]=afd_chb1(Wp,Ws,Rp,As);
% b=numerator polynomial coefficients of Ha(s);
% a=denominator polynomial coefficients of Ha(s);
% Wp=passband edge frequency in rad/sec;Wp>0;
% Ws=stopband edge frequency in rad/sec;Ws>Wp>0;
% Rp=passband ripple in +dB;(Rp>0);
% As=stopband attenuation in + dB;(As>0);
% ^^^^^^^^^^^^^^^^^^^^
if Wp<=0
error('passband edge must be larger than 0')
end
if Ws<=Wp
error('stopband dege must be larger than passband edge')
end
if(Rp<=0)|(As<0)
error('PB ripple and/or SB attenuation must be larger than 0')
end
ep=sqrt(10^(Rp/10)-1);
A=10^(As/20);
OmegaC=Wp;
OmegaR=Ws/Wp;
g=sqrt(A*A-1)/ep;
N=ceil(log10(g+sqrt(g*g-1))/log10(OmegaR+sqrt(OmegaR*OmegaR-1)));
fprintf('\n***Chebyshev-1 filter order=%2.0f\n',N)
[b,a]=u_chb2ap(N,As,OmegaC);
(4) afd_elip.m
%椭圆模拟低通滤波器原型设计
%afd_elip.m
function [b a]=afd_elip(Wp,Ws,Rp,As);
%Anolog lowpass filter design:Elliptic
%^^^^^^^^^^^^^^^^^^^^^^^^^^^^
%[b a]=afd_elip(Wp,Ws,Rp,As);
%b=Numberator coefficients of Ha(s)
% a=denominator polynomial coefficients of Ha(s);
% Wp=passband edge frequency in rad/sec;Wp>0;
% Ws=stopband edge frequency in rad/sec;Ws>Wp>0;
% Rp=passband ripple in +dB;(Rp>0);
% As=stopband attenuation in + dB;(As>0);
%%%%%%%%%%%%%%%%%%%%%%%%%
if Wp<=0
error('passband edge must be larger than 0')
end
if Ws<=Wp
error('stopband dege must be larger than passband edge')
end
if(Rp<=0)|(As<0)
error('PB ripple and/or SB attenuation must be larger than 0')
end
ep=sqrt(10^(Rp/10)-1);
A=10^(As/20);
OmegaC=Wp;
k=Wp/Ws;
k1=ep/sqrt(A*A-1);
capk=ellipke([k.^2 1-k.^2]);
capk1=ellipke([k1.^2 1-k1.^2]);
N=ceil(capk(1)*capk1(2)/(capk(2)*capk1(1)));
fprintf('\n***Elliptic Filter Order=%2.0f\n',N)
[b a]=u_elipap(N,Rp,As,OmegaC)

模拟低通原型滤波器的MATLAB设计相关推荐

[就酱的新征程]数字处理仿真与应用实验二：模拟低通 Butterworth 滤波器 IIR 滤波器（Matlab）
@数字处理仿真与应用实验二:模拟低通 Butterworth 滤波器 IIR 滤波器(Matlab) 这学期数字处理仿真与应用实验的记录不多说,进入正题吧~ 实验要求利用模拟滤波器原型设计 IIR ...
低通采样的matlab实现,基于matlab的FIR滤波器设计(低通,频率取样法)
基于matlab的FIR滤波器设计(低通,频率取样法)一.参考程序 M=63;%所需频率采样点个数 Wp=0.5*pi;%通带截止频率 m=0:(M+1)/2;%通频带上的采样点 Wm=2*pi*m. ...
matlab使用矩形窗设计一个具有线性相位的低通数字滤波器,用矩形窗设计一个FIR线性相位低通数字滤波器...
1 用矩形窗设计一个FIR线性相位低通数字滤波器.已知215.0Nc.求出nh并画出log20jeH曲线.分析此题给定的是理想线性相位低通滤波器故.-- 0- ccccjjdeeH 解deeHnhnj ...
低通采样定理 matlab,基于matlab的低通抽样定理仿真
基于matlab的低通抽样定理仿真 DSP 课程设计专业: 电子信息技术工程年级: 2011 级不姓名: 陈兰兰学号: 20113015 指导教师: 刘德春阿坝师专电子信息工程系 ...
低通采样定理 matlab,基于Matlab的低通抽样定理仿真.docx
DSP课程设计专业: 电子信息技术工程年级: 2011级不姓名: 陈兰兰学号: 指导教师: 刘德春阿坝师专电子信息工程系DSP课程设计专业: 电子信息技术工程年级: 2011级不姓名: ...
低通采样定理 matlab,通信原理MATLAB验证低通抽样定理实验报告
通信原理MATLAB验证低通抽样定理实验报告 p通信原理实验报告/pp一.实验名称/ppMATLAB验证低通采样定理/pp二.实验目的/pp1.掌握取样定理的工作原理./pp2.通过MATLAB编程实 ...
高阶低通无源滤波器的设计
题目:设某设备的工作频率为1MHZ,输入阻抗为72欧姆,干扰频率为5MHZ,根据原型滤波器参数计算,设计六阶低通无源滤波器. 图一原型滤波器参数表解:6阶原型滤波器的参数为:c1=0.5176F, ...
低通采样定理 matlab,低通采样定理
引入低通滤波器或提高低通滤波器的参数;该低通滤波器通常称为抗混叠滤波器抗混叠滤波器可限制信号的带宽,使之满足采样定理的条件.从理论上来说,这是可行的, ...... 3.理想低通滤波器的响应 3.8 ...
matlab设计高通系统,用matlab设计高通滤波器雪比切夫、fir两种方法课程设计HPF.doc...
课程设计通信电子电路课程设计通信电子电路课程设计 --数字滤波器的设计张静设计题目张静设计题生姓名生姓名光信息08-3班学光信息08-3班学号指导教师专业班级张静胡磊艾永春赵亚龙张腾 ...

模拟低通原型滤波器的MATLAB设计

模拟低通原型滤波器的MATLAB设计相关推荐

最新文章

热门文章