线性相位FIR滤波器hn
- 格式:ppt
- 大小:4.99 MB
- 文档页数:20
滤波器的线性相位和非线性相位设计滤波器在信号处理中扮演着重要的角色,它能够去除频域中不感兴趣的频率成分,同时保留我们所关注的信号频率。
然而,在滤波器设计中,存在线性相位和非线性相位两种设计选择。
本文将讨论滤波器的线性相位和非线性相位设计以及其相应的特点和应用。
一、线性相位设计线性相位滤波器是指滤波器的相位响应与信号频率成正比的情况。
具体而言,线性相位滤波器具有以下特点:1.1 来源与时域中因果系统线性相位滤波器的设计通常基于时域中因果系统,这意味着滤波器对信号的响应仅依赖于当前时刻以及之前的时刻。
这种设计能够保持信号的波形特性,不会引入额外的时间延迟。
1.2 相位响应与频率成正比线性相位滤波器的相位响应与信号频率成正比,从而在频域中引入了一个线性增加或减少的相位特性。
这种相位特性在某些应用中十分重要,例如音频信号处理中的均衡器。
1.3 启发性和直观性由于线性相位滤波器的相位响应与频率成正比,因此它具有较强的启发性和直观性。
在频域中,我们可以直观地分析滤波器对信号的相位响应,并进行相应的调整。
二、非线性相位设计非线性相位滤波器是指滤波器的相位响应与信号频率不成正比的情况。
具体而言,非线性相位滤波器具有以下特点:2.1 来源与时域中非因果系统非线性相位滤波器的设计通常基于时域中非因果系统,这意味着滤波器对信号的响应不仅取决于当前时刻以及之前的时刻,还可能涉及到之后的时刻。
这种设计在某些应用中能够引入更复杂的滤波效果。
2.2 非线性的相位变化非线性相位滤波器的相位响应与信号频率不成正比,因此在频域中引入了非线性的相位变化。
这种非线性相位特性在某些应用中能够提供更丰富的音效处理,例如混响效果。
2.3 指数函数的特殊相位特性非线性相位滤波器中,一种常见的设计是基于指数函数的相位特性。
这种设计能够产生一种特殊的相位延迟,通常用于音频合成和特殊效果处理中。
三、线性相位和非线性相位的应用线性相位滤波器和非线性相位滤波器在不同的应用中具有不同的优劣势。
西南科技大学课程设计报告课程名称:数字通信课程设计设计名称:线性相位FIR数字滤波器的设计姓名:学号:班级:指导教师:起止日期:2011.6.21-2011.7.3西南科技大学信息工程学院制课 程 设 计 任 务 书学生班级: 学生姓名: 学号:设计名称: 线性相位FIR 数字滤波器的设计 起止日期: 2011.6.21-2011.7.3 指导教师:设计要求:1、用窗函数法设计一个线性相位FIR 高通数字滤波器。
要求:FIR 高通数字滤波器指标为:)阻带衰减()通带衰减(度)数字阻带截止频率(弧度)数字通带截止频率(弧dB dBs A dB dB p R s p 4013.05.0====πωπω根据技术指标选择合适的窗形状,并绘制FIR 高通数字滤波器的幅度响应曲线和相位响应曲线;2、用窗函数法设计一个线性相位FIR 低通数字滤波器。
要求:FIR 低通数字滤波器指标为:)阻带衰减()通带衰减(度)数字阻带截止频率(弧度)数字通带截止频率(弧dB dBs A dB dB p R s p 5014.02.0====πωπω根据技术指标选择合适的窗形状,并绘制FIR 低通数字滤波器的幅度响应曲线和相位响应曲线;FIR 数字滤波器的设计可以使用matlab 工具箱中的函数课程设计学生日志时间设计内容2011.6.21-6.24 查阅资料,确定方案2011.6.25-6.26 设计总体方案2011.6.27-6.30 编写程序2011.7.1-7.3 撰写设计报告2011.7.4 答辩课程设计考勤表周星期一星期二星期三星期四星期五课程设计评语表指导教师评语:成绩:指导教师:年月日线性相位FIR 数字滤波器的设计一、 设计目的和意义 1.目的(1)掌握用窗函数法设计FIR 数字滤波器的原理与方法。
(2)熟悉线性相位FIR 数字滤波器的特性。
(3)了解各种窗函数对滤波特性的影响。
2.意义通过做这个设计,我们可以加深对线性相位FIR 数字滤波器原理以及设计方法的了解,有助于夯实进一步学习的基础。
设计题目及要求设计一个线性相位FIR 低通滤波器,技术指标如下:通带截止频率fp=1500Hz ,阻带起始频率fst=2250H ,通带允许的最大衰减为Rp=0.25dB ,阻带应达到的最小衰减为As=50dB 。
滤波器的采样频率为fs=15000Hz 。
设计原理 (包括滤波器工作原理、涉及到的matlab 函数的说明) 滤波器,顾名思义,其作是对输入信号起到滤波的作用数字滤波器滤波器结构x[k]为输入,h[k]为单位脉冲序列MATLAB 信号处理中提供的窗函数 (1)矩形窗 W=boxcar(N) (2)汉宁窗 W=hanning(N) (3)Bartlerr 窗 W=Bartlett(N) (4)Backman 窗 W=Backman (N ) (5)三角窗 W=triang (N ) (6)Kaiser 窗∑∑==-+-=M k k N k k k n x b k n y a n y 01][][][∑∞-∞=-=k k n x k h n y ][][][W=kaiser(n,beta)其中,beta是kaiser窗参数,影响窗旁瓣幅值的衰减率Kaiser窗用于滤波器设计时,选择性大,使用方便Nf=512;Nwin=20;%窗函数数据长度figure(1)for ii=1:4switch iicase 1w=boxcar(Nwin);stext='矩形窗';case 2w=hanning(Nwin);stext='汉宁窗';case 3w=hamming(Nwin);stext='哈明窗';case 4w=bartlett(Nwin);stext='Bartlett窗';end[y,f]=freqz(w,1,Nf);%求解窗函数特性,窗函数相当于一个数字滤波器mag=abs(y);%求得窗函数幅频特性posplot=['2,2,',int2str(ii)];subplot(posplot);plot(f/pi,20*log10(mag/max(mag)));%绘制窗函数的幅频特性xlabel('归一化频率');ylabel('振幅/dB');title(stext);grid on;FIR滤波器设计的主要方法函数设计方法说明工具函数窗函数法理想滤波器加窗处理Fir1(单频带)fir2(多频带)kaiserord最优化设计平方误差最小化逼近理想幅频响应或Park-McClellan算法产生等波纹滤波器Firls,remez,remezord约束最在满足最大误差限制条件下使整Fircls firclsl小二乘逼近个频带平方误差最小化升余弦函数具有光滑,正弦过渡带的低通滤波器设计Fircos各种窗函数的特点窗函数主瓣宽第一旁瓣相对主瓣衰减dB 矩形窗4pi/N -13汉宁窗8pi/N -31哈明窗8pi/N -41Bartlett 8pi/N -25Backman 12pi/N -57三角窗8pi/N -25Kaiser窗可调整可调整Chebyshew窗可调整可调整主旁瓣频率宽度与窗函数长度N有关。
滤波器的线性相位和非线性相位设计方法滤波器是一种常用的信号处理器件,它可以通过选择特定频率范围内的信号,对信号进行滤波和处理。
滤波器的设计涉及到很多方面,其中一个重要的考虑因素是相位特性。
本文将介绍滤波器的线性相位和非线性相位设计方法。
一、线性相位设计方法线性相位滤波器是指滤波器的相位响应与频率成线性关系。
线性相位滤波器一般使用FIR (Finite Impulse Response) 滤波器来实现,其特点是具有稳定的相移特性,适用于实时信号处理应用。
线性相位滤波器的设计方法有两种常用的方式:窗函数法和最小相位反演法。
1.1 窗函数法窗函数法是一种常用的设计线性相位滤波器的方法。
该方法的基本思想是将滤波器的频率响应与理想滤波器的频率响应进行近似拟合。
常见的窗函数有矩形窗、汉宁窗、黑曼窗等。
在窗函数法中,首先确定滤波器所需的通带、阻带和过渡带的频率范围,然后选择合适的窗函数进行设计。
通过对窗函数进行傅立叶变换,可以得到滤波器的冲激响应。
最后,通过将冲激响应作为滤波器的系数,即可实现线性相位的滤波器设计。
1.2 最小相位反演法最小相位反演法是另一种常用的设计线性相位滤波器的方法。
该方法的基本原理是通过对滤波器的幅度响应进行傅立叶变换,并计算其对数幅度谱,然后将对数幅度谱反变换得到滤波器的冲激响应。
最小相位反演法的优点是可以设计出更短的线性相位滤波器,适用于信号处理时延较为严格的应用场景。
然而,该方法的计算复杂度较高,需要进行频域的计算和反变换,因此在实际应用中需要根据具体情况进行权衡和选择。
二、非线性相位设计方法非线性相位滤波器是指滤波器的相位响应与频率不成线性关系。
非线性相位滤波器常用于对信号的组成部分进行时间或相位延迟的处理。
非线性相位滤波器的设计方法有FIR型和IIR型两种。
2.1 FIR型非线性相位滤波器FIR型非线性相位滤波器是通过设计多通的滤波器来实现的。
其基本思想是在滤波器的频域响应上引入不同频率的群延迟,从而实现非线性相位特性。
DN0409: 四种情况线性相位FIR 滤波器的幅度特性:通信与系统处理专业:王修启(031120507)情况(1):)1()(n N h n h --=偶对称,N 为奇数为奇数åå-+=----=-+-+=1212123)()21()()(N N n nj N j N n nj j e n h e N h en h e H w w w w在上式右端对第一、三项进行变量代换得到在上式右端对第一、三项进行变量代换得到åå-=-+----=÷øöçèæ----++-+--=211)21(2121121)21()21()21()(N m N m j N j N m m N j j e N m h e N h em N h e H w w w w由于÷øöçèæ-+=úûùêëé÷øöçèæ----=÷øöçèæ--2121121N mh m N N h m N h ,所以可以合并上面第一项和第三项,从而得到下面上面第一项和第三项,从而得到下面{}úúúûùêêêëé-++--=å-=---)21()21()(211)21(N h e e n N h ee H Nn n j n j N j j w w w wïþïýüïîïíì-+--=å-=--211)21()21(cos )()21(2Nn N jN h n n h n N h ew w 令21,2,1,212)(,21)0(-=úûùêëé--=÷øöçèæ-=N n n N h n a N ha 则 úúúûùêêêëé=å-=--21021)cos()()(N n N j jn n a ee H w ww 其频率响应幅度特性为å-==21)cos()()(N n n n a H w w ,关于ω=0,π,2π偶对称。
截止频率可变的线性相位FIR滤波器设计线性相位FIR滤波器是一种数字滤波器,能够在频域上实现线性相位特性。
在实际应用中,经常需要设计一种截止频率可变的线性相位FIR滤波器。
这种滤波器可以根据应用要求,实现对输入信号的不同频率成分进行滤波处理。
在开始设计之前,我们需要明确设计的目标和要求。
首先,我们要确定滤波器的截止频率范围,以及滤波器的通带和阻带要求。
其次,我们需要选择一种滤波器设计方法,例如频域设计方法或时域设计方法。
最后,我们需要确定滤波器的阶数,以及滤波器参数的数量和范围。
设计一种截止频率可变的线性相位FIR滤波器的一种常用方法是使用窗函数设计法。
具体步骤如下:1.确定滤波器的截止频率范围和通带/阻带要求。
根据应用需求设定截止频率范围,并确定通带/阻带的频率响应要求。
2.计算滤波器的阶数。
根据通带/阻带的频率响应要求,使用滤波器设计公式计算出滤波器的阶数。
3.选择滤波器的窗函数。
常用的窗函数包括矩形窗、汉宁窗、汉明窗、布莱克曼窗等。
根据设计要求选择合适的窗函数。
4.计算滤波器的频率响应。
根据窗函数的性质,计算出滤波器的频率响应。
5.选择滤波器的截止频率。
根据应用需求选择滤波器的截止频率。
6.通过滤波器设计公式计算出滤波器的权重系数。
7.对滤波器的权重系数进行归一化处理,保证其幅度范围在0到1之间。
8.实现滤波器的频率响应。
将归一化的权重系数作为滤波器的频率响应。
9.对滤波器频率响应进行离散化处理,得到滤波器的数字滤波器系数。
10.实现滤波器的时间域响应。
使用滤波器的数字滤波器系数,通过离散时间的卷积操作,实现滤波器的时间域响应。
以上是一种设计截止频率可变的线性相位FIR滤波器的一般步骤。
具体的设计方法和参数选择要根据具体的应用场景来确定。
在设计过程中,需要注意滤波器的性能要求和设计的复杂程度之间的平衡。
同时,还需要注意计算和处理的精度,以及实现滤波器的计算和存储的效率等因素。
通过截止频率可变的线性相位FIR滤波器的设计,我们可以实现对输入信号的不同频率成分进行滤波处理,从而满足不同应用场景的需求。
设计题目及要求设计一个线性相位FIR 低通滤波器,技术指标如下:通带截止频率fp=1500Hz ,阻带起始频率fst=2250H ,通带允许的最大衰减为Rp=0.25dB ,阻带应达到的最小衰减为As=50dB 。
滤波器的采样频率为fs=15000Hz 。
设计原理 (包括滤波器工作原理、涉及到的matlab 函数的说明) 滤波器,顾名思义,其作是对输入信号起到滤波的作用数字滤波器滤波器结构x[k]为输入,h[k]为单位脉冲序列MATLAB 信号处理中提供的窗函数 (1)矩形窗 W=boxcar(N) (2)汉宁窗 W=hanning(N) (3)Bartlerr 窗 W=Bartlett(N) (4)Backman 窗 W=Backman (N ) (5)三角窗 W=triang (N ) (6)Kaiser 窗∑∑==-+-=M k k N k k k n x b k n y a n y 01][][][∑∞-∞=-=k k n x k h n y ][][][W=kaiser(n,beta)其中,beta是kaiser窗参数,影响窗旁瓣幅值的衰减率Kaiser窗用于滤波器设计时,选择性大,使用方便Nf=512;Nwin=20;%窗函数数据长度figure(1)for ii=1:4switch iicase 1w=boxcar(Nwin);stext='矩形窗';case 2w=hanning(Nwin);stext='汉宁窗';case 3w=hamming(Nwin);stext='哈明窗';case 4w=bartlett(Nwin);stext='Bartlett窗';end[y,f]=freqz(w,1,Nf);%求解窗函数特性,窗函数相当于一个数字滤波器mag=abs(y);%求得窗函数幅频特性posplot=['2,2,',int2str(ii)];subplot(posplot);plot(f/pi,20*log10(mag/max(mag)));%绘制窗函数的幅频特性xlabel('归一化频率');ylabel('振幅/dB');title(stext);grid on;FIR滤波器设计的主要方法各种窗函数的特点主旁瓣频率宽度与窗函数长度N有关。