lsc滤波器设计.docx

名 称：信号分析与处理 题 目: ?基于窗函数法的 FIR数字高通滤波器设计院 系：自动化系 班 级：测控1303 班姓 名：廖思成 指导教师：金秀章 成 绩： 日期： 2016 年 6 月18日滤波器设计一、课程设计的目的?1、通过课程设计把自己在大学中所学的知识应用到实践当中。?2、深入了解利用Matlab设计FIR数字滤波器的基本方法。?3、在课程设计的过程中掌握程序编译及软件设计的基本方法。4、提高自己对于新知识的学习能力及进行实际操作的能力。?5、锻炼自己通过网络及各种资料解决实际问题的能力。二、主要设计内容 利用窗函数法、频率取样法及优化设计方法设计FIR滤波器，绘制出滤波器的特性图。利用所设计的滤波器对多个频带叠加的正弦信号进行处理，对比 滤波前后的信号时域和频域图，验证滤波器的效果三、设计原理?FIR?滤波器具有严格的相位特性，对于信号处理和数据传输是很重要的。目前?FIR滤波器的设计方法主要有三种：窗函数法、频率取样法和切比雪夫等波纹逼近的最优化设计方法。常用的是窗函数法和切比雪夫等波纹逼近的最优化设计方法。本实验中的窗函数法比较简单，可应用现成的窗函数公式，在技术指标要求高的时候是比较灵活方便的。四、用窗函数设FIR 滤波器的基本方法基本思路：从时域出发设计 h(n)逼近理想 hd(n)。设理想滤波器的单位响应在时域表达为hd(n),则Hd(n) 一般是无限长的，且是非因果的，不能直接作为FIR 滤波器的单位脉冲响应。要想得到一个因果的有限长的滤波器单位抽样响应 h(n)，最直接的方法是先将hd(n)往右平移，再迕行截断，即截取为有限长因果序列：h(n)=hd(n)w(n)，并用合适的窗函数迕行加权作为 FIR 滤波器的单位脉冲响应。按照线性相位滤波器的要求，线性相位FIR数字低通滤波器的单位抽样响应h(n)必须是偶对称的。对称中心必须等于滤波器的延时常数，即用矩形窗设计的FIR 低通滤波器，所设计滤波器的幅度函数在通带和阻带都呈现出振荡现象，且最大波纹大约为幅度的9%，返个现象称为吉布斯（Gibbs）效应。为了消除吉布斯效应，一般采用其他类型的窗函数。MATLAB 设计 FIR 滤波器有多种方法和对应的函数。窗函数设计法不仅在数字滤波器的设计中占有重要的地位，同时可以用于功率谱的估计，从根本上讲，使用窗函数的目的就是消除由无限序列的截短而引起的Gibbs现象所带来的影响。高通滤波器高通滤波器定义：高通滤波器的英文名称为high-pass filter ，其被全国科学技术名词审定委员会定义为让某一频率以上的信号分量通过，而对该频率以下的信号分量大大抑制的电容、电感与电阻等器件的组合装置。高通滤波器是一个使高频率比较容易通过而阻止低频率通过的系统。它去掉了信号中不必要的低频成分或者说去掉了低频干扰。其特性在时域及频域中可分别用冲激响应及频率响应描述。后者是用以频率为自变量的函数表示，一般情况下它是一个以复变量jω为自变量的的复变函数，以H（jω）表示。它的模H（ω）和幅角φ（ω）为角频率ω的函数，分别称为系统的“幅频响应”和“相频响应”，它分别代表激励源中不同频率的信号成分通过该系统时所遇到的幅度变化和相位变化。可以证明，系统的“频率响应”就是该系统“冲激响应”的傅里叶变换。当线性无源系统可以用一个N阶线性微分方程表示时，频率响应H（jω）为一个有理分式，它的分子和分母分别与微分方程的右边和左边相对应。最简单的高通滤波器称之为一阶高通滤波器，它的的特性一般用一阶线性微分方程表示，它的左边与一阶低通滤波器完全相同，仅右边是激励源的导数而不是激励源本身。如果激励源通过一个电阻给电容器构成一个充电回路，并以电阻两端的电压作为响应，就构成了一个以一阶微分方程描述的“一阶高通滤波器”，它的幅频响应在零频率处及其附近等于或接近于零，随着频率的增加，这个系统的幅频响应逐渐平滑地增加到趋近于1。也就是说，较低的频率通过该系统时，没有或几乎没有什么输出，而当较高的频率通过该系统时，将会受到较小的衰减。实际上，对于极高的频率而言，电容器相当于“短路”一样，这些频率，基本上都可以在电阻两端获得输出。换言之，这个系统适宜于通过高频率而对低频率有较大的阻碍作用，是一个最简单的“高通滤波器”。设计一个高通滤波器:频率特性为由设计任务可得：20lg0.05=-26dB，阻带衰减应该小于-26dB，因此可以选择汉宁窗。2、根据题目要求知：因为过渡带宽=0.5π-0.3π=0.2π，所以可以选择汉宁窗，由于N必须为奇数，故N=41。3、理想的频率响应为:4、对其进行傅里叶反变换得：且=0.40求得：将其向右平移得：加窗之后得：5、利用计算机编程进行处理：以下是程序代码：clear;N=41;f