当前位置：文档之家 › 数字信号处理课程设计-巴特沃斯带阻滤波器的设计

数字信号处理课程设计-巴特沃斯带阻滤波器的设计

  • 格式：doc
  • 大小：190.00 KB
  • 文档页数：11

郑州轻工业学院课程设计任务书题目模拟巴特沃斯带阻滤波器的设计专业、班级信息08-1 学号 200807070129 姓名秦伟伟主要内容、基本要求、主要参考资料等：1、主要内容1）设计巴特沃斯模拟低通原型滤波器；2）按频率变换设计巴特沃斯模拟带阻滤波器；3）对比模拟低通原型滤波器和模拟带阻滤波器并加以分析2、基本要求（1）编制MATLAB下的m文件实现主要内容。

（2）书写课程设计报告。

3、主要参考资料杨永双等编.数字信号处理实验指导书.郑州：郑州轻工业学院，2007 丁玉美等编著.数字信号处理第三版.西安电子科技大学出版社，2006 完成期限：指导教师签名：课程负责人签名：年月日第一章、理论简单介绍 (1)1、1 MATLAB概述 (1)1、2 滤波器设计 (2)1、3 基本定理 (4)第二章、设计目的、要求、指标 (5)2、1 设计目的 (5)2、2 设计要求 (5)2、3 设计指标 (6)第三章、程序代码和结果分析 (6)3、1 程序代码 (6)3、2结果与分析 (7)心得体会 (8)参考文献: (10)第一章、理论简单介绍1、1 MATLAB概述MATLAB 是一个可视化的计算程序，被广泛地应用在科学运算领域里。

它具有功能强大、使用简单等特点，内容包括：数值计算、符号计算、数据拟合、图形图像处理、系统模拟和仿真分析等功能。

此外，用Matlab 还可以进行动画设计、有限元分析等。

MATLAB系统包括五个主要部分：1）开发环境：这是一组帮助你使用MATLAB的函数和文件的工具和设备。

这些工具大部分是图形用户界面。

它包括MATLAB桌面和命令窗口，命令历史，和用于查看帮助的浏览器，工作空间，文件和查找路径。

2）MATLAB数学函数库：这里汇集了大量计算的算法，范围从初等函数如：求和，正弦，余弦和复数的算术运算，到复杂的高等函数如：矩阵求逆，矩阵特征值，贝塞尔(Bessel)函数和快速傅立叶变换等。

3）MATLAB语言：这是一种高水平的矩阵/数组语言，含有控制流语句，函数，数据结构，输入/输出，和面向对象编程特征。

它允许“小型编程”以迅速创立快速抛弃型程序，以及“大型编程”以创立完整的大型复杂应用程序。

4）句柄制图：这是MATLAB制图系统。

它包括高级别的二维、三维数据可视化，图像处理，动画，以及表现图形的命令。

它还包括低级别的命令，这使你不但能在MATLAB的应用中建立完整的图形用户界面，而且还能完全定制图形的外观。

5）MATLAB应用程序界面(API)：这是使你编写与MATLAB相合的C或Fortran程序的程序库。

它包括从MATLAB中调用程序（动态链接），调用MATLAB为计算引擎，和读写MAT-文件的设备。

1、2 滤波器设计在数字信号处理中，滤波占有极其重要的作用，滤波器是谱分析、雷达信号处理、通信信号处理应用中的基本处理算法。

目前常用的滤波器设计方法普遍采用Matlab 仿真。

1）巴特沃斯滤波器巴特沃斯滤波器模拟低通滤波器的平方幅频响应函数为：()()()N c A j H 22211ωωωω+== (1.1) 式中，c ω为低通滤波器的截止频率，N 为滤波器的阶数。

巴特沃斯滤波器的特点：通带内具有最大平坦的频率特性，且随着频率增大平滑单调下降；阶数愈高，特性愈接近矩形，过渡带愈窄，传递函数无零点。

这里的特性接近矩形，是指通带频率响应段与过渡带频率响应段的夹角接近直角。

通常该角为钝角，如果该角为直角，则为理想滤波器。

2）带阻滤波器的设计规范本设计中分别用h a (t)、H a (s)、Ha(j Ω)表示模拟滤波器的单位脉冲相应、系统函数、频率响应函数，三者的关系如下：Ha(s) = LT[ha(t)] =⎰∞∞-)(thae st-dtHa(jΩ)=FT[ha(t)] =⎰∞∞-)(thae t jΩ-dt可以用h a(t)、H a(s)、Ha(jΩ)中任一个描述模拟滤波器，也可以用线性常系数微分方程描述模拟滤波器。

但是设计模拟滤波器时，设计指标一般由幅频相应函数|Ha(jΩ)|给出，而模拟滤波器设计就是根据设计指标，求系统函数H a(s)。

工程实际中通常用所谓的损耗函数即衰减函数A（Ω）来描述滤波器的幅频响应特性，对归一化幅频响应函数，A（Ω）定义如下：A（Ω）=-20lg|Ha(jΩ)|=-10lg|Ha(jΩ)|2 dB模拟带阻滤波器的设计指标参数有Ω1p、Ω1s、Ω0、Ωsu、Ωpu。

Ω1p和Ωpu分别表示带阻滤波器的通带下边界频率和通带上边界频率；Ω1s和Ωsu分别表示带阻滤波器的阻带下边界频率和阻带上边界频率。

δ1和δ2分别称为通带和阻带波纹幅度。

带阻滤波器1、3 基本原理1）无失真传输所谓信号无失真传输是指输入信号通过系统后，输出信号的幅值和输入信号的幅值成正比。

允许有一定的延时，但没有波形上的畸变。

即系统的幅频响应|()ωj H |应为常数，相频响应()ωj H ∠应与频率ω成比例。

或者说，滤波器应具有无限宽的定值幅频与线性相频。

通常定义群延迟为信号系统的延迟时间为td(()ωωϕd d t d -=)，用函数表示为：|()ωj H |=C (常数)且t d =C(常数)2）理想滤波器理想滤波器应能无失真地传输有用信号，而又能完全抑制无用信号。

有用信号和无用信号往往占有不同的频带。

信号能通过滤波器的频带称为通带,信号被抑制的频带称为阻带。

由以前所学知识可以知道可能实现的,在具体实现的方面,我们只能想办法让实际滤波器的频率特性只能“逼近”理想滤波器。

滤波器的幅频响应在通带内不是完全平直的，而是呈波纹变化；在阻带内，幅频特性也不为零，而是衰减至某个值；在通带和阻带之间存在一个过渡带，而不是突然下降。

3）传递函数模拟滤波器的设计的理论基础通常在Laplace 域内进行讨论，模拟滤波器的技术指标可由平方幅值响应函数()()22ωωj H A =的形式给出，而()2ωj H 和传递函数H(s)存在下面关系： ()()()()ωωωj s s H s H j H A =-==|22 当给定模拟滤波器的技术指标后，由22)()(ωωj H A =求出A(-s 2)，再适当地选择分配零极点可求出H(s)。

为了使滤波器稳定，H(s)的极点必须落在s平面左半平面。

滤波器的零点选择可任取A(-s2)的一半零点，这是因为滤波器对Laplace域表示的传递函数并无特殊要求，但如果要求H(s)具有最小相位，零点也必须选择在s左半平面。

第二章、设计目的、要求、指标2、1 设计目的本次数字电路课程设计的主要目的有以下几点：1.熟练掌握MATLAB 软件系统；2.使学生的数字信号知识系统化、完整化；3.掌握在MATLAB环境下进行信号产生、信号变换、滤波器设计的方法；4. 初步了解和掌握工程实践中的设计思想和设计方法，具有开发信号分析、仿真系统的基本能力；5. 为后继课程的教学打下扎实基础；2、2 设计要求基于巴特沃斯滤波器设计中心频率为200Hz，带宽为150Hz的模拟带阻滤波器；（1）给定模拟滤波器的性能指标，选择合适的滤波器类型等；（2）确定滤波器阶数；（3）设计模拟低通原型滤波器；（4）按频率变换设计模拟带阻滤波器；2、3 设计指标模拟低通滤波器的设计指标有Wp, Wp, αp 和αs 。

Wp ：通带截止频率；Ws ：阻带截止频率；αp ：通带中最大衰减系数；αs ：阻带最小衰减系数。

αp 和αs 一般用dB 数表示。

对于单调下降的幅度特性，可表示成：( 如果Ω=0处幅度已归一化到1，即|Ha(j0)|=1,αp 和αs 表示为：第三章、程序代码和结果分析3、1 程序代码%设计中心频率为200Hz ，带宽为150Hz 的模拟带阻滤波器for n=1:6; %设计阶数为1~6的Butterworth 模拟原型滤波器[z,p,k]=buttap(n); %[z0，p0，k0]=buttap(N) 用于计算N 阶巴特沃斯归一化(3dB 截止频率Ωc=1)模拟低通原型滤波器系统函数的零、极点和增益因子[b,a]=zp2tf(z,p,k); %将零点极点增益形式转化传递函数形式[H,w]=freqs(b,a,0:0.01:2); %计算指定频率点的复数频率响应 subplot(1,1,1),plot(w,abs(H)); %绘制幅度响应xlabel('w/wc');ylabel('Butterworth ︳H(jw) ︳');title('Butterworth 低通原型滤波器(wc=1)')text(0.4,0.85,'n=1') %作必要的标记text(0.7,0.8,'n=2')text(0.89,0.7,'n=4')text(1.18,0.1,'n=6')hold on;figure(1) grid on;2222(0)(0)10lg 10lg ()()a a p s a s a p H j H j H j H j αα==ΩΩ2210lg ()10lg ()p a p s a s H j H j αα=-Ω=-Ωendfigure(2)hold on;for n=1:6; %设计阶数为1~6的Butterworth模拟带阻滤波器[z,p,k]=buttap(n); %[z0，p0，k0]=buttap(N) 用于计算N阶巴特沃斯归一化(3dB截止频率Ωc=1)模拟低通原型滤波器系统函数的零、极点和增益因子[b,a]=zp2tf(z,p,k); %将零点极点增益形式转化传递函数形式f0=200;bw=150; %中心点频率和阻带带宽[bt,at]=lp2bs(b,a,f0,bw); %频率转换（低通变带阻）[Ht,ft]=freqs(bt,at,0:2:400); %计算带阻滤波器指定频率点的复数频率响应plot(ft,abs(Ht));%绘制幅度响应xlabel('f/hz');ylabel('Butterworth ︳H(jw) ︳');title('Butterworth 带阻滤波器(f0=200,bw=150)');text(68,0.9,'n=1') %作必要的标记text(115,0.8,'n=2')text(305,0.8,'n=4')text(306,0.9,'n=6')grid on;end3、2 结果与分析1）仿真结果图5.1 Butterworth型模拟原型低通滤波器图5.2 Butterworth型模拟带阻滤波器2）结果分析从图5.2可以看到，中心频率为200hz，带宽为150hz，满足设计要求。

合集下载

数字信号处理实验数字巴特沃思滤波器的设计

数字信号处理实验数字巴特沃思滤波器的设计

数字信号处理实验数字巴特沃思滤波器的设计数字信号处理技术是现代通信、音频、图像等领域中不可或缺的一门技术。

数字信号处理的核心是数字滤波器设计，本文将介绍一种常用的数字滤波器——数字巴特沃斯滤波器的设计方法。

一、数字滤波器简介数字滤波器是将连续时间信号转换成离散时间信号，实现对离散时间信号的滤波处理，具有实时性好、精度高、可重复性强等优点。

数字滤波器有两种类型：有限冲激响应（FIR）滤波器和无限冲激响应（IIR）滤波器。

二、数字巴特沃斯滤波器数字巴特沃斯滤波器是一种常用的IIR滤波器，其主要特点是具有平坦的通/阻带，通/阻带边缘陡峭。

因此在实际应用中，数字巴特沃斯滤波器应用较为广泛。

数字巴特沃斯滤波器的设计方法一般包括以下步骤：确定滤波器类型、确定通/阻带的截止频率、确定滤波器的阶数、计算滤波器的系数。

1、确定滤波器类型在实际应用中，数字巴特沃斯滤波器有四种类型：低通、高通、带通和带阻滤波器，应根据实际需求选择。

2、确定通/阻带的截止频率通常情况下，固定本例中采用的是低通滤波器，需要确定的就是通带和阻带的截止频率。

对于低通滤波器，通带截止频率ωc应该比信号频率fs的一半小，阻带截止频率ωs 应该比ωc大一些，通常ωs/ωc取0.5~0.7比较好。

滤波器的阶数一般是与滤波器的性能相关的。

阶数越高，性能越好，但同时计算量也会更大。

在实际应用中，一般取4~8的阶数即可。

4、计算滤波器的系数根据上述参数计算滤波器的系数，这里介绍两种常用的方法：一种是脉冲响应不变法（Impulse Invariant Method），另一种是双线性变换法（Bilinear Transformation）。

脉冲响应不变法是一种较为简单的设计方法，但由于其数字滤波器与连续时间滤波器之间的不同，可能会引入一定程度的失真。

双线性变换法可以使二阶系统和一阶系统的增益分别为1和0dB，这是一种比较理想的设计方法。

四、实验步骤本实验采用Matlab软件进行数字滤波器的设计，具体步骤如下：1、打开Matlab软件，新建一个.m文件；2、输入需要滤波的数字信号，此处可以使用Matlab自带的signal工具箱中的一些模拟信号；4、使用filter函数实现数字滤波器对信号的滤波过程；5、通过比较信号的频谱图，评估滤波器的性能。

设计巴特沃斯数字带通滤波器

设计巴特沃斯数字带通滤波器

设计巴特沃斯数字带通滤波器，要求通带范围为：0.25π rad ≤ω≤0.45π rad,通带最大衰减为3dB ，阻带范围为0≤ω≤0.15π rad 和0.55π rad ≤ω≤πrad ，阻带最小衰减为40dB 。

利用双线性变换设计，写出设计过程，并用MATLAB 绘出幅频和相频特性曲线。

设计思路及计算：（1）确定技术指标，求得数字边缘频率：10.25Pp rad ωπ=，20.45Pp rad ωπ=，3p a dB =10.15Ps rad ωπ=，20.55Ps rad ωπ=，40s a dB =（2）将数字带通滤波器的技术指标转换为模拟带通滤波器技术指标： 用双线性变换法，则2tan 2T ωΩ=，可得 112tan tan 0.1250.4142/2Pp Pp rad s T ωπΩ===222tan tan 0.2250.8541/2Pp Pp rad s T ωπΩ=== 112tan tan 0.0750.2401/2Ps Ps rad s T ωπΩ=== 222tan tan 0.275 1.1708/2Ps Ps rad s T ωπΩ=== （3）将带通滤波器的指标转换为模拟低通指标。

模拟低通归一化边界频率为：1Lp Ω=，()()2212221 1.9748Ps Pp Pp Ls Ps Pp Pp Ω-ΩΩΩ==ΩΩ-Ω（4）确定低通滤波器阶数N4020100.01s δ-==，()2211lg 1lg 10.01 6.76812lg 1.97482lg ss p N δ⎛⎫⎛⎫-- ⎪ ⎪⎝⎭⎝⎭≥==⎛⎫Ω ⎪ ⎪Ω⎝⎭取N ＝7。

（5）c c ΩΩ=Ω=1c Ω≈巴特沃兹模拟滤波器：(217)14711H (),()j K a k kk s p es p π++===-∏再由双线性变换即可得到所求。

代码实现：>> [N,Wn]=buttord([.25 .45],[.15 .55],3,40) N = 7 Wn = 0.2482 0.4525 >> [b,a]=butter(7,[.2482 .4525]) b = Columns 1 through 10 0.0001 0 -0.00070 0.0022 0 -0.0036 0 0.0036 0 Columns 11 through 15 -0.0022 0 0.00070 -0.0001 a = Columns 1 through 10 1.0000 -5.3094 16.2918 -34.7303 56.9401 -74.5112 80.0108 -71.1129 52.6364 -32.2233 Columns 11 through 15 16.1673 -6.4607 1.9827 -0.4217 0.0523>> [h,w]=freqz(b,a,100); >>subplot(211) >>h1=20*log10(abs(h)); >>plot(w/pi,h1); >>axis([0 1 -50 10]); >>subplot(212)>>plot(w/pi,angle(h))则滤波器传递函数为：24681012141234567890.00010.00070.00220.00360.00360.00220.00070.0001()1 5.309416.291834.730356.940174.511280.010871.112952.636432.223316.1673z z z z z z z H z z z z z z z z z z z ----------------⎛⎫-+- ⎪ ⎪+-+-⎝⎭=-+-+-+-+-+10111213146.4607 1.98270.42170.0523z z z z -----⎛⎫ ⎪ ⎪-+-+⎝⎭幅频和相频曲线：由上图可知，已满足设计要求0.10.20.30.40.50.60.70.80.91-50-40-30-20-1001000.10.20.30.40.50.60.70.80.91-4-2024。

数字信号处理巴特沃斯滤波器的汉明窗设计

数字信号处理巴特沃斯滤波器的汉明窗设计

1 绪论1.1课程设计背景《数字信号处理》课程是一门理论和技术发展十分迅速、应用非常广泛的前沿性学科，他的理论性和实践性都很强，他的特点是：要求的数学知识多，包括高等代数、数值分析、概率统计、随机过程等。

要求掌握的基础知识强，网络理论、信号与系统是本课程的理论基础与其他学科密切相关，即与通信理论、计算机、微电子技术不可分，又是人工智能、模式识别、神经网络等新兴学科的理论基础之一。

学生在学习这门课程时，普遍感到数字信号处理的概念抽象，对其中的分析方法与基本理论不能很好地理解与掌握。

因此，如何帮助学生理解与掌握课程中的基本概念、基本原理、基本分析方法以及综合应用所学知识解决实际问题的能力，是本课程教学中所要解决的关键问题。

为了配合《数字信号处理》专业基础课的理论教学，我们在电子信息工程专业教学计划中安排了二周的《数字信号处理》课程设计，他是针对《数字信号处理》的基础理论和算法进行实践环节的一个综合训练，以便学习巩固所学的知识，加强理论和实际结合的能力，培养学生的综合设计能力与实际工作能力。

Matlab语言是一种广泛应用于工程计算及数值分析领域的新型高级语言，Matlab功能强大、简单易学、编程效率高，深受广大科技工作者的欢迎。

特别是Matlab还具有信号分析工具箱，不需具备很强的编程能力，就可以很方便地进行信号分析、处理和设计。

因此，选择用Matlab进行课程设计。

1.2 课程设计目的1.掌握数字信号处理的基本概念，基本理论和基本方法。

2.熟悉离散信号和系统的时域特性。

3.掌握序列快速傅里叶变换方法。

4.学会MATLAB的使用，掌握MATLAB的程序设计方法。

5.掌握利用MATLAB对语音信号进行频谱分析。

6.掌握滤波器的网络结构。

2 课程设计原理2.1用窗函数法设计FIR 滤波器根据过渡带宽及阻带衰减要求，选择窗函数的类型并估计窗口长度N （或阶数M=N-1），窗函数类型可根据最小阻带衰减As 独立选择，因为窗口长度N 对最小阻带衰减As 没有影响，在确定窗函数类型以后，可根据过渡带宽小于给定指标确定所拟用的窗函数的窗口长度N ，设待求滤波器的过渡带宽为Δw ，它与窗口长度N 近似成反比，窗函数类型确定后，其计算公式也确定了，不过这些公式是近似的，得出的窗口长度还要在计算中逐步修正，原则是在保证阻带衰减满足要求的情况下，尽量选择较小的N ，在N 和窗函数类型确定后，即可调用MATLAB 中的窗函数求出窗函数wd （n ）。

巴特沃斯带阻滤波器设计方法

巴特沃斯带阻滤波器设计方法

巴特沃斯带阻滤波器设计方法巴特沃斯滤波器是一种常用的电子滤波器，用于滤除特定频率范围内的信号。

带阻滤波器是巴特沃斯滤波器的一种特殊类型，用于在特定频率范围内阻止信号通过。

设计带阻滤波器需要考虑一些关键因素，包括滤波器的阶数、截止频率、通带和阻带的衰减等。

首先，确定滤波器的阶数。

巴特沃斯滤波器的阶数决定了滤波器的斜率和衰减率。

一般来说，阶数越高，滤波器的性能越好，但相应的设计和实现也更复杂。

根据实际需求和可实现的复杂度，选择适当的阶数。

其次，确定滤波器的截止频率。

带阻滤波器需要指定两个截止频率，即通带截止频率和阻带截止频率。

通带是允许信号通过的频率范围，而阻带则是需要被滤除的频率范围。

根据应用需求和信号特性，确定这两个截止频率的数值。

然后，计算滤波器的元件数值。

根据巴特沃斯滤波器的设计公式，计算出滤波器的元件数值，包括电容和电感的数值。

这些数值决定了滤波器的实际工作特性，需要精确计算和选择。

接下来，进行滤波器的电路设计和实现。

根据计算得到的元件数值，设计滤波器的电路图并选择合适的元件进行实现。

在设计过程中需要注意元件的精度、稳定性和可获得性，确保设计的可实现性和稳定性。

最后，进行滤波器的调试和性能验证。

实现滤波器电路后，需要进行调试和性能验证。

通过实际测试和测量，验证滤波器的性能是否符合设计要求，如通带衰减、阻带衰减、相位特性等。

总的来说，设计巴特沃斯带阻滤波器需要考虑阶数、截止频率、元件数值、电路设计和性能验证等多个方面。

综合考虑这些因素，可以设计出满足特定要求的带阻滤波器电路。