%用脉冲响应不变法设计一个椭圆数字滤波器
wp1=0.3*pi;wp2=0.7*pi;
ws1=0.1*pi;ws2=0.9*pi;
Rp=1;As=15;
Fs=2000;T=1/Fs;
Omgp1=wp1*Fs;Omgp2=wp2*Fs;
Omgp=[Omgp1,Omgp2];
Omgs1=ws1*Fs;Omgs2=ws2*Fs;
Omgs=[Omgs1,Omgs2];
bw=Omgp2-Omgp1; w0=sqrt(Omgp1*Omgp2);
[N,Omgn]=ellipord(Omgp,Omgs,Rp,As,'s')
[z0,p0,k0]=ellipap(N,Rp,As);
ba1=k0*real(poly(z0));
aa1=real(poly(p0));
[ba,aa]=lp2bp(ba1,aa1,w0,bw);
[bd,ad]=impinvar(ba,aa,Fs)
[H,w]=freqz(bd,ad);
dbH=20*log10((abs(H)+eps)/max(abs(H)));
subplot(2,2,1),plot(w/pi,abs(H));
subplot(2,2,2),plot(w/pi,angle(H)/pi);
subplot(2,2,3),plot(w/pi,dbH); subplot(2,2,4),zplane(bd,ad);
该
系统的传递函数是z z z z z z z z z z z z H 654321543212917.0162.09128.00937.09247.03254.012409.06059.00677.02048.01777
.0592.0)(-----------+-+-+--+-+-=
%用脉冲响应不变法设计一个切比雪夫Ⅱ型数字带阻滤波器
wp1=0.1*pi;wp2=0.9*pi;
ws1=0.3*pi;ws2=0.7*pi;
Rp=1;As=40;
Fs=2000;T=1/Fs;
Omgp1=wp1*Fs;Omgp2=wp2*Fs;
Omgp=[Omgp1,Omgp2];
Omgs1=ws1*Fs;Omgs2=ws2*Fs;
Omgs=[Omgs1,Omgs2];
bw=Omgp2-Omgp1; w0=sqrt(Omgp1*Omgp2);
[N,Omgn]=cheb2ord(Omgp,Omgs,Rp,As,'s');
[z0,p0,k0]=cheb2ap(N,As);
ba1=k0*real(poly(z0));
aa1=real(poly(p0));
[ba,aa]=lp2bs(ba1,aa1,w0,bw);
[bd,ad]=impinvar(ba,aa,Fs);
%模拟滤波器与数字滤波器的冲击响应
t=0:T:(30*T); nt=length(t);
ha=impulse(ba,aa,t);
h=impz(bd,ad,nt);
subplot(2,1,1),plot(t,ha*T,'r');hold on;
stem(t,h,'k');
title('模拟滤波器与数字滤波器的冲击响应');
%模拟滤波器与数字滤波器的幅频响应
wb=[0:Fs]*2*pi;
Ha=freqs(ba,aa,wb);
H=freqz(bd,ad,wb/Fs);
subplot(2,1,2),plot(wb/(2*pi),abs(Ha)/max(abs(Ha)),'r');hold on; plot(wb/(2*pi),abs(H)/max(abs(H)),'k');
title('模拟滤波器与数字滤波器的幅频响应');
以下两个滤波器都是切比雪夫I型数字滤波器，不是巴特沃尔滤波器，请使用者注意！ 1.带通滤波器 function y=bandp(x,f1,f3,fsl,fsh,rp,rs,Fs) %带通滤波 %使用注意事项：通带或阻带的截止频率与采样率的选取范围是不能超过采样率的一半%即，f1,f3,fs1,fsh,的值小于Fs/2 %x:需要带通滤波的序列 % f 1：通带左边界 % f 3：通带右边界 % fs1：衰减截止左边界 % fsh：衰变截止右边界 %rp：边带区衰减DB数设置 %rs：截止区衰减DB数设置 %FS：序列x的采样频率 % f1=300;f3=500;%通带截止频率上下限 % fsl=200;fsh=600;%阻带截止频率上下限 % rp=0.1;rs=30;%通带边衰减DB值和阻带边衰减DB值 % Fs=2000;%采样率 % wp1=2*pi*f1/Fs; wp3=2*pi*f3/Fs; wsl=2*pi*fsl/Fs; wsh=2*pi*fsh/Fs; wp=[wp1 wp3]; ws=[wsl wsh]; % % 设计切比雪夫滤波器； [n,wn]=cheb1ord(ws/pi,wp/pi,rp,rs);
[bz1,az1]=cheby1(n,rp,wp/pi); %查看设计滤波器的曲线 [h,w]=freqz(bz1,az1,256,Fs); h=20*log10(abs(h)); figure;plot(w,h);title('所设计滤波器的通带曲线');grid on; y=filter(bz1,az1,x); end 带通滤波器使用例子 %-------------- %带通滤波器测试程序 fs=2000; t=(1:fs)/fs; ff1=100; ff2=400; ff3=700; x=sin(2*pi*ff1*t)+sin(2*pi*ff2*t)+sin(2*pi*ff3*t); figure; subplot(211);plot(t,x); subplot(212);hua_fft(x,fs,1); % y=filter(bz1,az1,x); y=bandp(x,300,500,200,600,0.1,30,fs); figure; subplot(211);plot(t,y); subplot(212);hua_fft(y,fs,1); %调用到的hua_fft()函数代码如下 function hua_fft(y,fs,style,varargin) %当style=1,画幅值谱；当style=2,画功率谱;当style=其他的，那么花幅值谱和功率谱
1 绪论 (1) 1.1 数字滤波器的优点 (2) 1.2 数字滤波器的发展概况 (3) 1.3 开发工具 (4) 1.3.1 MATLAB软件简介 (4) 1.3.2 MATLAB特点 (5) 2 数字滤波器理论研究 (6) 2.1 滤波器的设计 (6) 2.2 数字滤波器的定义 (7) 2.3 滤波器的设计步骤 (8) 2.4 数字滤波器的类型 (8) 2.5 滤波器的选择 (9) 2.5.1 FIR和IIR数字滤波器的比较 (9) 2.5.2 FIR或IIR滤波器的选取原则 (10) 2.6 数字滤波器的实现方法 (10) 3 IIR滤波器的设计 (11) 3.1 数字滤波器 (11) 3.2 IIR数字滤波器设计方法 (12) 3.2.1用脉冲相应不变法设计IIR数字滤波器 (13) 3.2．2 用双线性变换法设计IIR数字滤波器 (15) 4 IIR数字带通滤波器设计过程 (18) 4.1 设计步骤： (18) 4.2 程序流程框图： (19) 4.3 MATLAB程序： (19) 结论 (21) 参考文献 (22) 致谢 (23)
基于MATLAB的IIR数字带通滤波器设计 1 绪论 凡是有能力进行信号处理的装置都可以称为滤波器。在近代电信设备和各类控制系统中，滤波器应用极为广泛；在所有的电子部件中，使用最多，技术最为复杂的要算滤波器了。滤波器的优劣直接决定产品的优劣，所以，对滤波器的研究和生产历来为各国所重视。 1917年美国和德国科学家分别发明了LC滤波器，次年导致了美国第一个多路复用系统的出现。20世纪50年代无源滤波器日趋成熟。自60年代起由于计算机技术、集成工艺和材料工业的发展，滤波器发展上了一个新台阶，并且朝着低功耗、高精度、小体积、多功能、稳定可靠和价廉方向努力，其中小体积、多功能、高精度、稳定可靠成为70年代以后的主攻方向。导致RC有源滤波器、数字滤波器、开关电容滤波器和电荷转移器等各种滤波器的飞速发展，到70年代后期，上述几种滤波器的单片集成已被研制出来并得到应用。80年代，致力于各类新型滤波器的研究，努力提高性能并逐渐扩大应用范围。90年代至现在主要致力于把各类滤波器应用于各类产品的开发和研制。当然，对滤波器本身的研究仍在不断进行。 我国广泛使用滤波器是50年代后期的事，当时主要用于话路滤波和报路滤波。经过半个世纪的发展，我国滤波器在研制、生产和应用等方面已纳入国际发展步伐，但由于缺少专门研制机构，集成工艺和材料工业跟不上来，使得我国许多新型滤波器的研制应用与国际发展有一段距离。 目前数字滤波器的设计有许多现成的高级语言设计程序，但他们都存在设计效率较低，不具有可视图形，不便于修改参数等缺点，而Matlab为数字滤波的研究和应用提供了一个直观、高效、便捷的利器。他以矩阵运算为基础，把计算、可视化、程序设计融合到了一个交互式的工作环境中。尤其是Matlab工具箱使各个领域的研究人员可以直观方便地进行科学研究与工程应用。其中的信号处理工具箱、图像处理工具箱、小波工具箱等更是为数字滤波研究的蓬勃发展提供了可能。 1
m a t l a b程序之——滤波器(带通-带阻)
matlab程序之——滤波器（带通，带阻） 以下两个滤波器都是切比雪夫I型数字滤波器，不是巴特沃尔滤波器，请使用者注意！ 1.带通滤波器 function y=bandp(x,f1,f3,fsl,fsh,rp,rs,Fs) %带通滤波 %使用注意事项：通带或阻带的截止频率与采样率的选取范围是不能超过采样率的一半 %即，f1,f3,fs1,fsh,的值小于 Fs/2 %x:需要带通滤波的序列 % f 1：通带左边界 % f 3：通带右边界 % fs1：衰减截止左边界 % fsh：衰变截止右边界 %rp：边带区衰减DB数设置 %rs：截止区衰减DB数设置 %FS：序列x的采样频率 % f1=300;f3=500;%通带截止频率上下限 % fsl=200;fsh=600;%阻带截止频率上下限 % rp=0.1;rs=30;%通带边衰减DB值和阻带边衰减DB值 % Fs=2000;%采样率 % wp1=2*pi*f1/Fs; wp3=2*pi*f3/Fs; wsl=2*pi*fsl/Fs; wsh=2*pi*fsh/Fs; wp=[wp1 wp3]; ws=[wsl wsh]; % % 设计切比雪夫滤波器； [n,wn]=cheb1ord(ws/pi,wp/pi,rp,rs); [bz1,az1]=cheby1(n,rp,wp/pi); %查看设计滤波器的曲线 [h,w]=freqz(bz1,az1,256,Fs); h=20*log10(abs(h));
figure;plot(w,h);title('所设计滤波器的通带曲线');grid on; y=filter(bz1,az1,x); end 带通滤波器使用例子 %-------------- %带通滤波器测试程序 fs=2000; t=(1:fs)/fs; ff1=100; ff2=400; ff3=700; x=sin(2*pi*ff1*t)+sin(2*pi*ff2*t)+sin(2*pi*ff3*t); figure; subplot(211);plot(t,x); subplot(212);hua_fft(x,fs,1); % y=filter(bz1,az1,x); y=bandp(x,300,500,200,600,0.1,30,fs); figure; subplot(211);plot(t,y); subplot(212);hua_fft(y,fs,1); %调用到的hua_fft()函数代码如下 function hua_fft(y,fs,style,varargin) %当style=1,画幅值谱；当style=2,画功率谱;当style=其他的，那么花幅值谱和功率谱 %当style=1时，还可以多输入2个可选参数 %可选输入参数是用来控制需要查看的频率段的 %第一个是需要查看的频率段起点 %第二个是需要查看的频率段的终点 %其他style不具备可选输入参数，如果输入发生位置错误 nfft= 2^nextpow2(length(y));%找出大于y的个数的最大的2的指数值（自动进算最佳FFT步长nfft） %nfft=1024;%人为设置FFT的步长nfft y=y-mean(y);%去除直流分量 y_ft=fft(y,nfft);%对y信号进行DFT，得到频率的幅值分布 y_p=y_ft.*conj(y_ft)/nfft;%conj()函数是求y函数的共轭复数，实数的共轭复数是他本身。
目录 一．设计概述 二．设计任务及要求 2.1 设计任务 2.2 设计要求 三．设计方案 3.1设计结构 3.2元件参数的理论推导 3.3仿真电路构建 3.4仿真电路分析四．所用器件 五．实验结果 5.1 实验数据记录 5.2 实验数据分析六．实验总结 6.1 遇到的主要问题 6.2 解决问题的措施 6.3 实验反思与收获 附图 参考文献
一．设计概述 根据允许的通过的频率范围，可以将滤波器分为低通滤波器，高通滤波器，带通滤波器和带阻滤波器4种。其中，带通滤波器是指允许某一频率范围内的频率分量通过，其他范围的频率分量衰减到极低水平的滤波器。 在滤波器中，信号能够通过的范围成为通频带或通带，信号受到很大衰减或完全被抑制的频率范围成为阻带，通带和阻带之间的界限称为截止频率。对于一个理想的带通滤波器，通带范围内则完全平坦，对传输信号基本没有增益的衰减作用，其次，通带之外的所有频率均能被完全衰减掉，通带和阻带之间存在一定的过渡带。 在带通滤波器的实际设计过程中，主要参数包括中心频率f0,频带宽度BW，上限截止频率fH和下限截止频率fL。一般情况下，为使滤波器在任意频段都具有良好的频率分辨能力，可采用固定带宽带通滤波器（如收音机的选频）。所选带宽越窄，则频率选择能力越高。但为了覆盖所要检测的整个频率范围，所需要的滤波器数量就很大。因此，在很多场合，固定带宽带通滤波器不一定做成固定中心频率的，而是利用一个参考信号，使滤波器中心频率跟随参考信号的频率而变化，其中，参考信号是由信号发生器提供的。上述可便中心频率的固定带宽带通滤波器，经常用于滤波和扫描跟踪滤波应用中。 二．设计任务及要求 1）设计任务 带通滤波器的设计方案有很多，本实验将采用高通滤波器和低通滤波器级联的设计方案实现一个带通滤波器，通过多级反馈，减少干扰信号对滤波器的影响。为了检测滤波电路的通带特性，设计一个带宽检测电路，通过发光二极管的亮灭近似检测电路的带宽范围。 设计要求 2）设计要求 （1）性能指标要求 1.输入信号：有效值为1V的电压信号。 2.输出信号中心频率f0通过开关切换，分别为500Hz 1.5KHz 3KHz 10KHz 误差10%。 3.带通滤波器带宽BW
数字信号处理综合设计 实验报告 一、实验目的: （1） 深刻理解滤波器的设计指标及根据指标进行数字滤波器设计的过程 （2） 了解滤波器在通信系统中的使用 二、实验步骤: 1.通过SYSTEMVIEW 规划整个系统，确定系统的采样频率、观测时间、细化并设计整个系统，仿真调整并不断改进达到正确调制、正确滤波、正确解调的目的。（参考文件zhan3.svu ） （1） 检查滤波器的波特图，看是否达到预定要求； （2） 检查幅度调制的波形以及相加后的信号的波形和频谱是否正常； （3） 检查解调后的的基带信号是否正常，分析波形变形的原因和解决措施； （4） 实验中必须体现带通滤波器的物理意义和在实际中的使用价值。 2.熟悉matlab 中的仿真系统； 3.将1.中设计的SYSTEMVIEW （如zhan3.svu ）系统移植到matlab 中的仿真环境中，使其达到相同的效果； 4.或者不用仿真环境，编写程序实现该系统，并验证调制解调前后的信号是否一致。 实验总共提供三个单元的时间（6节课）给学生，由学生自行学习和自行设 sin ω2 基带信号2
计和移植 三、实验内容： 1.使用MATLAB软件中的图形化工具按照zhan3连接带通滤波器、低通滤波器等如下图： 2.其中各参数如下： （1）Plus Generator设置如下： 脉冲类型为Sample based，幅值1，周期100，脉冲宽度50，采样时间0.001s
（2）载波设置如下： 100HZ的载波：幅度为5，采样时间为0.001s 300HZ的载波：幅度为5，采样时间为0.001s
带通滤波器1：
实验八 有源滤波器的设计 一．实验目的 1． 学习有源滤波器的设计方法。 2． 掌握有源滤波器的安装与调试方法。 3． 了解电阻、电容和Q 值对滤波器性能的影响。 二．预习要求 1． 根据滤波器的技术指标要求，选用滤波器电路，计算电路中各元件的数值。设计出 满足技术指标要求的滤波器。 2. 根据设计与计算的结果，写出设计报告。 3. 制定出实验方案，选择实验用的仪器设备。 三．设计方法 有源滤波器的形式有好几种，下面只介绍具有巴特沃斯响应的二阶滤波器的设计。 巴特沃斯低通滤波器的幅频特性为： n c uo u A j A 21)(??? ? ??+= ωωω ， n=1，2，3，. . . （1） 写成： n c uo u A j A 211) (??? ? ??+=ωωω （2） )(ωj A u 其中A uo 为通带内的电压放大倍数，ωC A uo 为截止角频率，n 称为滤波器的阶。从（2） 式中可知，当ω=0时，（2）式有最大值1； 0.707A uo ω=ωC 时，（2）式等于0.707，即A u 衰减了3dB ；n 取得越大，随着ω的增加，滤波器的输出电压衰减越快，滤波器的幅频特性越接近于理想特性。如图1所示。ω 当 ω>>ωC 时， n c uo u A j A ??? ? ??≈ωωω1 )( （3） 图1低通滤波器的幅频特性曲线
两边取对数，得： lg 20c uo u n A j A ωω ωlg 20)(-≈ （4） 此时阻带衰减速率为： -20ndB/十倍频或-6ndB/倍频，该式称为衰减估算式。 表1列出了归一化的、n 为1 ~ 8阶的巴特沃斯低通滤波器传递函数的分母多项式。 在表1的归一化巴特沃斯低通滤波器传递函数的分母多项式中，S L = c s ω，ωC 是低通 滤波器的截止频率。 对于一阶低通滤波器，其传递函数： c c uo u s A s A ωω+= )( （5） 归一化的传递函数： 1 )(+= L uo L u s A s A （6） 对于二阶低通滤波器，其传递函数：2 22)(c c c uo u s Q s A s A ωωω++ = （7） 归一化后的传递函数： 1 1)(2 ++= L L uo L u s Q s A s A （8） 由表1可以看出，任何高阶滤波器都可由一阶和二阶滤波器级联而成。对于n 为偶数的高阶滤波器，可以由2n 节二阶滤波器级联而成；而n 为奇数的高阶滤波器可以由2 1-n 节二
基于matlab的带通、带阻滤波器设计实例 以下两个滤波器都是切比雪夫I型数字滤波器，不是巴特沃尔滤波器，请使用者注意！ 1.带通滤波器 function y=bandp(x,f1,f3,fsl,fsh,rp,rs,Fs) %带通滤波 %使用注意事项：通带或阻带的截止频率与采样率的选取范围是不能超过采样率的一半%即，f1,f3,fs1,fsh,的值小于 Fs/2 %x:需要带通滤波的序列 % f 1：通带左边界 % f 3：通带右边界 % fs1：衰减截止左边界 % fsh：衰变截止右边界 %rp：边带区衰减DB数设置 %rs：截止区衰减DB数设置 %FS：序列x的采样频率 % f1=300;f3=500;%通带截止频率上下限 % fsl=200;fsh=600;%阻带截止频率上下限 % rp=0.1;rs=30;%通带边衰减DB值和阻带边衰减DB值 % Fs=2000;%采样率 % wp1=2*pi*f1/Fs; wp3=2*pi*f3/Fs; wsl=2*pi*fsl/Fs; wsh=2*pi*fsh/Fs; wp=[wp1 wp3]; ws=[wsl wsh]; % % 设计切比雪夫滤波器； [n,wn]=cheb1ord(ws/pi,wp/pi,rp,rs); [bz1,az1]=cheby1(n,rp,wp/pi); %查看设计滤波器的曲线 [h,w]=freqz(bz1,az1,256,Fs); h=20*log10(abs(h)); figure;plot(w,h);title('所设计滤波器的通带曲线');grid on; y=filter(bz1,az1,x); end 带通滤波器使用例子 %-------------- %带通滤波器测试程序 fs=2000; t=(1:fs)/fs; ff1=100;
课程设计报告 专业班级 课程 题目 学号 学生姓名 指导教师 年月
一、设计题目：IIR 数字带通滤波器设计 二、设计目的 1、巩固所学理论知识。 2、提高综合运用所学理论知识独立分析和解决问题的能力。 3、更好地将理论与实践相结合。 4、掌握信号分析与处理的基本方法与实现。 5、熟练使用MATLAB 语言进行编程实现。 三、设计要求 采用适当方法基于MATLAB 设计一个IIR 带通滤波器，其中带通的中心频率为ωp0=0.5π，;通带截止频率ωp1=0.4π,ωp2=0.6π;通带最大衰减αp =3dB;阻带最小衰减αs =15dB;阻带截止频率ωs2=0.7π. 四、设计原理 1.用脉冲相应不变法设计IIR 数字滤波器 利用模拟滤波器来设计数字滤波器，也就是使数字滤波器能模仿模拟滤波器的特性，这种模仿可以从不同的角度出发。脉冲响应不变法是从滤波器的脉冲响应出发，使数字滤波器的单位脉冲响应序列h (n )模仿模拟滤波器的冲激响应 h a (t )，即将h a (t )进行等间隔采样，使h (n )正好等于h a (t )的采样值，满足 h (n )=h a (nT ) 式中,T 是采样周期。 如果令H a (s )是h a(t )的拉普拉斯变换，H (z )为h (n )的Z 变换，利用采样序列的 Z 变换与模拟信号的拉普拉斯变换的关系得 (1-1) 则可看出，脉冲响应不变法将模拟滤波器的S 平面变换成数字滤波器的Z 平面，这个从s 到z 的变换z =e sT 是从S 平面变换到Z 平面的标准变换关系式。 ??? ?? -= Ω-= ∑∑ ∞ -∞=∞ -∞ ==k T j s X T jk s X T z X k a s k a e z sT π21 )(1) (
二阶有源模拟带通滤波器设计 摘要 滤波器是一种具有频率选择功能的电路，它能使有用的频率信号通过。而同时抑制(或衰减)不需要传送频率范围内的信号。实际工程上常用它来进行信号处理、数据传送和抑制干扰等，目前在通讯、声纳、测控、仪器仪表等领域中有着广泛的应用。 以往这种滤波电路主要采用无源元件R、L和C组成，60年代以来，集成运放获得迅速发展，由它和R、C组成的有源滤波电路，具有不用电感、体积小、重量轻等优点。此外，由于集成运放的开环电压增益和输入阻抗都很高，输出阻抗比较低，构成有源滤波电路后还具有一定的电压放大和缓冲作用。 通常用频率响应来描述滤波器的特性。对于滤波器的幅频响应，常把能够通过信号的频率范围定义为通带，而把受阻或衰减信号的频率范围称为阻带，通带和阻带的界限频率叫做截止频率。 滤波器在通带内应具有零衰减的幅频响应和线性的相位响应，而在阻带内应具有无限大的幅度衰减。按照通带和阻带的位置分布，滤波器通常分为低通滤波器、高通滤波器、带通滤波器和带阻滤波器。文中结合实例，介绍了设计一个二阶有源模拟带通滤波器。 设计中用RC网络和集成运放组成，组成电路选用LM324不仅可以滤波，还可以进行放大。 关键字：带通滤波器 LM324 RC网络
目录 目录 (2) 第一章设计要求 (3) 1.1基本要求 (3) 第二章方案选择及原理分析 (4) 2.1.方案选择 (4) 2.2 原理分析 (5) 第三章电路设计 (7) 3.1 实现电路 (7) 3.2参数设计 (7) 3.3电路仿真 (9) 1.仿真步骤及结果 (9) 2.结果分析 (11) 第四章电路安装与调试 (12) 4.1实验安装过程 (12) 4.2 调试过程及结果 ..................................................................................................... 错误!未定义书签。 4.2.1 遇到的问题 .................................................................................................. 错误!未定义书签。 4.2.2 解决方法 ...................................................................................................... 错误!未定义书签。 4.2.3 调试结果与分析 (12) 结论 (13) 参考文献 (14)
燕山大学 课程设计说明书题目：椭圆带通滤波器的设计 学院（系）：电气工程学院 年级专业： 10级精仪二班 学号： 学生姓名： 指导教师：刘永红 教师职称：副教授
燕山大学课程设计（论文）任务书 课程名称：数字信号处理课程设计 基层教学单位：指导教师： 学号学生姓名专业（班级）10级精仪二班设计题目2椭圆带通滤波器设计 设 计 技术参数采样频率为100Hz,采样点数100，低频、中频、高频信号频率分别为5Hz、 15Hz 、30Hz 设计要求产生一个连续信号，包含低频率，中频，高频分量，对其进行采样，进行频谱分析。设计带通滤波器对信号进行滤波处理，观察滤波后信号的频谱。 （熟悉函数freqz,butter,filter,fft） 参 考资料数字信号处理方面资料MATLAB方面资料 周次前半周后半周工 作计划收集消化资料、学习MATLAB软件，进行 相关参数计算。 编写仿真程序、调试。 指导教师签字基层教学单位主任签字 说明：此表一式四份，学生、指导教师、基层教学单位、系部各一份。 年月日
目录 第1章摘要．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．4 第2章引言．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．4 第3章基本原理．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．5 3.1 模拟滤波器的基本原理．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．5 3.2 椭圆滤波器的特点．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．5 第4章设计过程．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．6 4.1 椭圆滤波器设计结构图．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．6 4.2 设计椭圆模拟滤波器．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．7 4.3 模拟滤波器的MATLAB实现和滤波器分析．．．．．．．．．．．．．．．．．．．．．．．．7 第5章仿真程序和仿真图．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．． 10 5.1、%连续信号的产生及采样．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．． 10 5.2、%椭圆带通滤波器的设计．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．11 5.3、%信号通过椭圆带通滤波器的波形图．．．．．．．．．．．．．．．．．．．．．．．．．12 5.4、信号通过椭圆带通滤波器的仿真图．．．．．．．．．．．．．．．．．．．．．．．．．．12 第6章分析及总结．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．． 14 心得体会．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．15 参考文献．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．15
LC椭圆函数带通滤波器设计 要求带通滤波器，在15kHz～ZOkHz的频率范围内，衰减最大变化1dB，低于14.06kHz和高于23kHz频率范围，最小衰减为50dB,Rs=RL=10kΩ。 ③运行Filter Solutions程序。点击“阻带频率”输人框，在“通带波纹（dB）”内输人0.18，在“通带频率”内输人1，在“阻带频率”内输人1.456，选中“频率单位－弧度”逻辑框。在“源阻抗”和“负载阻抗”内输人1。 ④点击“确定阶数”控制钮打开第二个面板。在“阻带衰减（dB）”内输人50，点击“设置最小阶数”按钮并点击“关闭”，主控制面板上形式出“6阶”，选中“偶次阶模式”逻辑框。 ⑤点击“电路”按钮。Filter s。lutions提供了两个电路图。选择“无源滤波器1”，如图1（a）所示。 ⑥这个滤波器必须变换为中心频率ω0=1的归一化带通滤波器。带通滤波器的Q 值为： 把所有的电感量和电容值都乘以Qbp°然后用电感并联每一个电容、用电容串联每一个电感使其谐振频率为ω0=1，该网络被变换为带通滤波器。使用的谐振元仵是原元件值的倒数，如图1（b）所示。 ⑦按照图1的方式转换Ⅱ型支路。
变换后的滤波器见图1（c）。在原理图下标出了以rad/s为单位的谐振频率。 ⑧用中心频率fo=17.32kHz和阻抗10kΩ对滤波器进行去归一化以完成设计。将所有的电感乘以Z/FSF，所有的电容除以z×FSF，其中z=104， FSF=2πfe=1.0882×105。最终的滤波器见图1（d）。图1（c）中的归一化谐振频率直接乘以几何中心频率fo=17.32kHz即可得到谐振频率。频率响应见图1（e）。
带通滤波电路设计一．设计要求 （1）信号通过频率范围 f 在100 Hz至10 kHz之间； （2）滤波电路在 1 kHz 电路的幅频衰减应当在 的幅频响应必须在± 1 kHz 时值的± 3 dB 1 dB 范围内，而在 范围内； 100 Hz至10 kHz滤波 （3）在10 Hz时幅频衰减应为26 dB ，而在100 kHz时幅频衰减应至少为16 dB 。 二．电路组成原理 由图（ 1）所示带通滤波电路的幅频响应与高通、低通滤波电路的幅频响应进行比较， 不难发现低通与高通滤波电路相串联如图（2），可以构成带通滤波电路，条件是低通滤波电路的截止角频率 W H大于高通电路的截止角频率 W L，两者覆盖的通带就提供了一个带通响应。 V I V O 低通高通 图（ 1） 1 W H低通截止角频率 R1C1 1 W L高通截止角频率 R2C2 必须满足W L │A│ O │A│ O │A│ O 低通 W w H 高通 W w L 带通 W W w L H 图（ 2） 三．电路方案的选择 参照教材 10.3.3 有源带通滤波电路的设计。这是一个通带频率范围为100HZ-10KHZ的带通滤波电路，在通带内我们设计为单位增益。根据题意，在频率低端f=10HZ 时，幅频响应至少衰减 26dB。在频率高端 f=100KHZ 时，幅频响应要求衰减不小于16dB。因此可以选择一个二阶高通滤波电路的截止频率fH=10KHZ,一个二阶低通滤波电路的fL=100HZ，有源器件仍选择运放 LF142，将这两个滤波电路串联如图所示，就构成了所要求的带通滤波电路。 由教材巴特沃斯低通、高通电路阶数n 与增益的关系知 A vf1 =1.586 ，因此，由两级串联的带通滤波电路的通带电压增益（Avf1 ） 2=（ 1.586 ）2=2.515, 由于所需要的通带增益为0dB, 因此在低通滤波器输入部分加了一个由电阻R1、 R2组成的分压器。 学院毕业设计(论文) 摘要 在模拟信号的处理中，滤波器起着重要的作用，一个好的滤波电路对一个产品的整个性能有着重要的影响。然而用传统的方法设计无源的高阶的滤波电路，不仅耗时，计算也比较复杂。因此本文在总结传统方法的基础上，运用计算机进行辅助设计，不仅省时，而且效率高。 本文对无源低通椭圆滤波器进行了研究。首先阐述了滤波器设计的理论，在此基础上，了解了常用的滤波器的原理和结构。接着根据椭圆滤波器的结构和原理，得出了理论上的椭圆滤波器。 然后分别利用滤波器设计软件Filter Solutions和MATLAB软件对椭圆滤波器进行设计和仿真。在比较两者优点的基础上，对设计出的电路进行了改进。 最后，对仿真后的电路，利用PROTEL 99SE画出电路的原理图，生成PCB板，并对元器件的值进行了修正，以便在实际中选择元器件。对焊接好的电路，进行了测试，结果表明设计出的椭圆滤波器符合规定的设计要求。 关键词：滤波器，传递函数，椭圆函数 无源低通椭圆滤波器的设计和实现 Abstract The filters play the vital role in processing analog signals. A good filter circuit has important influence on the whole performance of a product. However, in traditional methods, designing the passive and higher-order filter circuits wastes a great deal time and the computation is quite complex. Therefore, on the basis of summarizing traditional methods, this paper proposes the novel method which uses computers to aid design. The method makes time saved and efficiency improved. In the paper, the passive low pass elliptic filter is studied. Firstly, the theory about filter design is explained, and usual filters are introduced. On the basis of understanding structure and principle of the elliptic filter, the elliptic filter is theoretically designed. And then separately in the light of filter design software-Filter Solution and MATLAB, elliptic filter is realized and simulated. According to their respective advantage, designed circuit is improved. Finally, Schematic and PCB of designed circuit are charted by using Protel 99. In order to choose proper components, the values of components are amended. And sealed circuit is tested. The result has indicated that proposed elliptic filter satisfies design demands. Keywords: filters; transfer function; elliptic function 基于matlab的FIR低通-高通-带通-带阻滤波器设计 ————————————————————————————————作者：————————————————————————————————日期： 北京师范大学 课程设计报告 课程名称： DSP 设计名称：FIR 低通、高通带通和带阻数字滤波器的设计姓名： 学号: 班级： 指导教师： 起止日期： 课程设计任务书 学生班级： 学生姓名： 学号： 设计名称： FIR 低通、高通带通和带阻数字滤波器的设计 起止日期： 指导教师： 设计目标： 1、采用Kaiser 窗设计一个低通FIR 滤波器 要求： 采样频率为8kHz ； 通带：0Hz~1kHz ，带内波动小于5%； 阻带：1.5kHz ，带内最小衰减：Rs=40dB 。 2、采用hamming 窗设计一个高通FIR 滤波器 要求： 通带截至频率wp=rad π6.0, 阻带截止频率ws=rad π4.0， 通带最大衰减dB p 25.0=α,阻带最小衰减dB s 50=α 3、采用hamming 设计一个带通滤波器 低端阻带截止频率 wls = 0.2*pi ； 低端通带截止频率 wlp = 0.35*pi ； 高端通带截止频率 whp = 0.65*pi ； 高端阻带截止频率 whs = 0.8*pi ； 4、采用Hamming 窗设计一个带阻FIR 滤波器 要求： 通带：0.35pi~0.65pi ，带内最小衰减Rs=50dB ； 阻带：0~0.2pi 和0.8pi~pi ，带内最大衰减：Rp=1dB 。 FIR 低通、高通带通和带阻数字滤波器的设计 一、 设计目的和意义 1、熟练掌握使用窗函数的设计滤波器的方法，学会设计低通、带通、带阻滤波器。 2、通过对滤波器的设计，了解几种窗函数的性能，学会针对不同的指标选择不同的窗函数。 二、 设计原理 一般，设计线性相位FIR 数字滤波器采用窗函数法或频率抽样法，本设计采用窗函数法，分别采用海明窗和凯泽窗设计带通、带阻和低通。 如果所希望的滤波器的理想频率响应函数为)(jw d e H ，如理想的低通，由信号系统的知识知道，在时域系统的冲击响应h d (n)将是无限长的，如图2、图3所示。 H d (w) -w c w c 图2 图3 若时域响应是无限长的，则不可能实现，因此需要对其截断，即设计一个FIR 滤波器频率响应∑-=-=1 0)()(N n jwn jw e n h e H 来逼近)(jw d e H ，即用一个窗函数w(n)来 截断h d (n)，如式3所示： )()()(n w n h n h d = （式1）。 最简单的截断方法是矩形窗，实际操作中，直接取h d (n)的主要数据即可。 )(n h 作为实际设计的FIR 数字滤波器的单位脉冲响应序列，其频率响应函数为： ∑-=-=1 0)()(N n jwn jw e n h e H （式2） 令jw e z =，则 ∑-=-=1 0)()(N n n z n h z H （式3）， 式中，N 为所选窗函数)(n w 的长度。 带通滤波器设计步骤 1、根据需求选择合适的低通滤波器原型 2、把带通滤波器带宽作为低通滤波器的截止频率，根据抑制点的频率距离带通滤波器中心频点距离的两倍作为需要抑制的频率，换算抑制频率与截止频率的比值，得出m 的值，然后根据m 值选择低通滤波器的原型参数值。 滤波器的时域特性 任何信号通过滤波器都会产生时延。Bessel filter 是特殊的滤波器在于对于通带内的所有频率而言，引入的时延都是恒定的。这就意味着相对于输入，输出信号的相位变化与工作的频率是成比例的。而其他类型的滤波器（如Butterworth, Chebyshev,inverse Chebyshev,and Causer ）在输出信号中引入的相位变化与频率不成比例。相位随频率变化的速率称之为群延迟（group delay ）。群延迟随滤波器级数的增加而增加。 模拟滤波器的归一化 归一化的滤波器是通带截止频率为w=1radian/s, 也就是1/2πHz 或约0.159Hz 。这主要是因为电抗元件在1弧度的时候，描述比较简单，XL=L, XC=1/C ，计算也可以大大简化。归一化的无源滤波器的特征阻抗为1欧姆。归一化的理由就是简化计算。 Bessel filter 特征：通带平坦，阻带具有微小的起伏。阻带的衰减相对缓慢，直到原理截止频率高次谐波点的地方。原理截止频率点的衰减具有的经验公式为n*6dB/octave ，其中，n 表示滤波器的阶数，octave 表示是频率的加倍。例如，3阶滤波器，将有18dB/octave 的衰减变化。正是由于在截止频率的缓慢变化，使得它有较好的时域响应。 Bessel 响应的本质截止频率是在与能够给出1s 延迟的点，这个点依赖于滤波器的阶数。 逆切比雪夫LPF 原型参数计算公式（Inverse Chebyshev filter parameters calculate equiations ） ) (cosh )(cosh 11Ω=--Cn n 其中 1101.0-=A Cn , A 为抑制频率点的衰减值，以dB 为单位；Ω为抑制频率与截止频率的比值 例：假设LPF 的3dB 截止频率为10Hz,在15Hz 的频点需要抑制20dB,则有： 95.91020*1.0==Cn ；Ω=15/10=1.5 1.39624.0988.2) 5.1(cosh )95.9(cosh 11===--n ,因此，滤波器的阶数至少应该为4 基于MA TLAB的IIR数字带通滤波器设计 摘要 窗函数法在IIR 数字滤波器的设计中有着广泛的应用, 但这不是最优化的设计。介绍了一种基于等波纹切比雪夫逼近准则的IIR 数字滤波器的最优化设计方法,通过MA TLAB 的仿真实现, 证明了该方法是一种最优化的设计。传统的数字滤波器设计方法繁琐且结果不直观,本文利用MA TLAB具有强大的科学计算和图形显示这一优点,与窗函数法设计理论相结合共同设计IIR数字滤波器,不但使设计结果更加直观,而且提高了滤波器的设计精度,从而更好地达到预期效果。 关键词：IIR数字滤波器;窗函数，等波纹切比雪夫逼近，MA TLAB 仿真 ABSTRACT Window function method in the design of IIR digital filter has a wide range of applications, but this is not the most optimal design. Such as corrugated paper, a Chebyshev approximation for IIR digital filter criteria for the optimization design method to achieve through the MA TLAB simulation proved that the method is one of the most optimized design. Conventional digital filter design method is cumbersome and results are not intuitive, this paper, MA TLAB has a powerful scientific computing and graphics display the advantages, with the window function method combines design theory to design IIR digital filter design results not only more intuitive, but also improve the accuracy of the filter design to better achieve the desired results. KEY WORDS: IIR digital filters，window function，such as ripple Chebyshev approximation，MA TLAB simulation 目录 引言.............................................第页 第1 章数字滤波器................................第页 第2 章IIR数字滤波器设计方法......................第页 2.1用脉冲相应不变法设计IIR数字滤波器..........第页 2.2 脉冲响应不变法优缺点........................第页 2.3用双线性变换法设计IIR数字滤波器............第页 2.4双线性变换法优缺点..........................第页 第3章IIR数字带通滤波器设计过程...................第页 3.1设计步骤.....................................第页 3.2程序流程框图.................................第页 3.3 MA TLAB程序..................................第页 第四章结果及分析.................................第页 第五章总结.......................................第页 参考文献..........................................第页 致谢..............................................第页 附录..............................................第页椭圆滤波器的实现
基于matlab的FIR低通高通带通带阻滤波器设计
带通滤波器设计步骤
基于MATLAB的IIR数字带通滤波器设计
相关主题
文本预览