FIR滤波器设计与实现
- 格式:doc
- 大小:1.26 MB
- 文档页数:15
数字信号处理实验报告-FIR滤波器的设计与实现在数字信号处理中,滤波技术被广泛应用于时域处理和频率域处理中,其作用是将设计信号减弱或抑制被一些不需要的信号。
根据滤波器的非线性抑制特性,基于FIR(Finite Impulse Response)滤波器的优点是稳定,易设计,可以得到较强的抑制滤波效果。
本实验分别通过MATLAB编程设计、实现、仿真以及分析了一阶低通滤波器和平坦通带滤波器。
实验步骤:第一步:设计一阶低通滤波器,通过此滤波器对波型进行滤波处理,分析其对各种频率成分的抑制效果。
为此,采用零极点线性相关算法设计滤波器,根据低通滤波器的特性,设计的低通滤波器的阶次为n=10,截止频率为0.2π,可以使设计的滤波器被称为一阶低通滤波器。
第二步:设计平坦通带滤波器。
仿真证明,采用兩個FIR濾波器組合而成的阻礙-提升系統可以實現自定義的總三值響應的設計,得到了自定義的總三值響應函數。
实验结果:1、通过MATLAB编程,设计完成了一阶低通滤波器,并通过实验仿真得到了一阶低通滤波器的频率响应曲线,证明了设计的滤波器具有良好的低通性能,截止频率为0.2π。
在该频率以下,可以有效抑制波形上的噪声。
2、设计完成平坦通带滤波器,同样分析其频率响应曲线。
从实验结果可以看出,此滤波器在此频率段内的通带性能良好,通带范围内的信号透过滤波器后,损耗较小,滞后较小,可以满足各种实际要求。
结论:本实验经过实验操作,设计的一阶低通滤波器和平坦通带滤波器具有良好的滤波特性,均已达到预期的设计目标,证明了利用非线性抑制特性实现FIR滤波处理具有较强的抑制滤波效果。
本实验既有助于深入理解FIR滤波器的设计原理,也为其他应用系统的设计和开发提供了指导,进而提高信号的处理水平和质量。
fir数字滤波器设计与软件实现数字信号处理实验原理FIR数字滤波器设计的基本原理是从理想滤波器的频率响应出发,寻找一个系统函数,使其频率响应尽可能逼近滤波器要求的理想频率响应。
为了实现这一目标,通常会采用窗函数法进行设计。
这种方法的基本思想是,将理想滤波器的无限长单位脉冲响应截断为有限长因果序列,并用合适的窗函数进行加权,从而得到FIR滤波器的单位脉冲响应。
在选择窗函数时,需要考虑其频率响应和幅度响应。
常见的窗函数包括矩形窗、三角形窗、汉宁窗、汉明窗、布莱克曼窗和凯泽窗等。
每种窗函数都有其特定的特性,如主瓣宽度、旁瓣衰减等。
根据实际需求,可以选择合适的窗函数以优化滤波器的性能。
在软件实现上,可以使用各种编程语言和信号处理库进行FIR滤波器的设计和实现。
例如,在MATLAB中,可以使用内置的`fir1`函数来设计FIR滤波器。
该函数可以根据指定的滤波器长度N和采样频率Fs,自动选择合适的窗函数并计算滤波器的系数。
然后,可以使用快速卷积函数`fftfilt`对输入信号进行滤波处理。
此外,还可以使用等波纹最佳逼近法来设计FIR数字滤波器。
这种方法的目标是找到一个最接近理想滤波器频率响应的实数序列,使得在所有可能的实
数序列中,该序列的误差平方和最小。
通过优化算法,可以找到这个最优序列,从而得到性能更优的FIR滤波器。
总的来说,FIR数字滤波器设计与软件实现数字信号处理实验原理是基于对理想滤波器频率响应的逼近和优化,通过选择合适的窗函数和算法,实现信号的滤波处理。
实验四FIR数字滤波器设计与软件实现
实验目的:
掌握FIR数字滤波器的设计与软件实现方法,了解滤波器的概念与基
本原理。
实验原理:
FIR数字滤波器全称为有限脉冲响应数字滤波器,其特点是具有有限
长度的脉冲响应。
滤波器通过一系列加权系数乘以输入信号的延迟值,并
将这些值相加得到输出信号。
FIR滤波器的频率响应由滤波器系数所决定。
实验步骤:
1.确定所需的滤波器的设计规格,包括截止频率、通带波纹、阻带衰
减等。
2.选择适当的滤波器设计方法,如窗函数、最佳近似法、最小二乘法等。
3.根据所选方法,计算滤波器的系数。
4.在MATLAB环境下,使用滤波器的系数实现滤波器。
5.输入所需滤波的信号,经过滤波器进行滤波处理。
6.分析输出的滤波信号,观察滤波效果是否符合设计要求。
实验要求:
1.完成FIR数字滤波器的设计和软件实现。
2.对比不同设计方法得到的滤波器性能差异。
3.分析滤波结果,判断滤波器是否满足设计要求。
实验器材与软件:
1.个人电脑;
2.MATLAB软件。
实验结果:
根据滤波器设计规格和所选的设计方法,得到一组滤波器系数。
通过
将滤波器系数应用于输入信号,得到输出滤波信号。
根据输出信号的频率
响应、通带波纹、阻带衰减等指标,评估滤波器的性能。
实验注意事项:
1.在选择设计方法时,需要根据滤波器要求和实际情况进行合理选择。
2.在滤波器实现过程中,需要注意滤波器系数的计算和应用。
3.在实验过程中,注意信号的选择和滤波结果的评估方法。
FIR滤波器的MATLAB设计与实现FIR滤波器(Finite Impulse Response Filter)是一种数字滤波器,其特点是其响应仅由有限长度的序列决定。
在MATLAB中,我们可以使用信号处理工具箱中的函数来设计和实现FIR滤波器。
首先,需要明确FIR滤波器的设计目标,包括滤波器类型(低通、高通、带通、带阻)、通带和阻带的频率范围、通带和阻带的增益等。
这些目标将决定滤波器的系数及其顺序。
在MATLAB中,我们可以使用`fir1`函数来设计FIR滤波器。
该函数的使用方式如下:```matlabh = fir1(N, Wn, type);```其中,`N`是滤波器长度,`Wn`是通带边缘频率(0到0.5之间),`type`是滤波器的类型('low'低通、'high'高通、'bandpass'带通、'stop'带阻)。
该函数会返回一个长度为`N+1`的滤波器系数向量`h`。
例如,如果要设计一个采样频率为10kHz的低通滤波器,通带截止频率为2kHz,阻带频率为3kHz,可以使用以下代码:```matlabfc = 2000; % 通带截止频率h = fir1(50, fc/(fs/2), 'low');```上述代码中,`50`表示滤波器的长度。
注意,滤波器的长度越大,滤波器的频率响应越陡峭,但计算成本也更高。
在设计完成后,可以使用`freqz`函数来分析滤波器的频率响应。
例如,可以绘制滤波器的幅度响应和相位响应曲线:```matlabfreqz(h);```除了使用`fir1`函数外,MATLAB还提供了其他函数来设计FIR滤波器,如`fir2`、`firpm`、`firls`等,具体使用方式可以参考MATLAB的文档。
在实际应用中,我们可以将FIR滤波器应用于音频处理、图像处理、信号降噪等方面。
例如,可以使用FIR滤波器对音频信号进行去噪处理,或者对图像进行锐化处理等。
fir滤波器设计实验报告一、实验目的本次实验的目的是设计FIR滤波器,从而实现信号的滤波处理。
二、实验原理FIR滤波器是一种数字滤波器,它采用有限长的冲激响应滤波器来实现频率选择性的滤波处理。
在FIR滤波器中,系统的输出只与输入和滤波器的系数有关,不存在反馈环路,因此具有稳定性和线性相位的特性。
FIR滤波器的设计最常采用Window法和最小二乘法。
Window法是指先对理想滤波器的频率特性进行窗函数的处理,再通过离散傅里叶变换来得到滤波器的时域响应。
最小二乘法则是指采用最小二乘法来拟合理想滤波器的频率特性。
本次实验采用的是Window法。
三、实验步骤1.设计滤波器的频率响应特性:根据实际需要设计出需要的滤波器的频率响应特性,通常采用理想滤波器的底通、高通、带通、带阻等特性。
2.选择窗函数:根据设计的滤波器的频率响应特性选择相应的窗函数,常用的窗函数有矩形窗、汉宁窗、汉明窗等。
3.计算滤波器的时域响应:采用离散傅里叶变换将设计的滤波器的频率响应特性转化为时域响应,得到滤波器的冲激响应h(n)。
4.归一化:将得到的滤波器的冲激响应h(n)进行归一化处理,得到单位加权的滤波器系数h(n)。
5.实现滤波器的应用:将得到的滤波器系数h(n)应用于需要滤波的信号中,通过卷积的方式得到滤波后的信号。
四、实验结果以矩形窗为例,设计一阶低通滤波器,截止频率为300Hz,采样频率为8000Hz,得到的滤波器系数为:h(0)=0.0025h(1)=0.0025滤波效果良好,经过滤波后的信号频率响应相对于滤波前有较明显的截止效应。
五、实验总结通过本次实验,我们掌握了FIR滤波器的设计方法,窗函数的选择和离散傅里叶变换的应用,使我们能够更好地处理信号,实现更有效的信号滤波。
在日常工作和学习中,能够更好地应用到FIR滤波器的设计和应用,提高信号处理的精度和效率。
2北京邮电大学信息与通信工程学院概述:IIR 和FIR 比较IIR 与FIR 性能比较IIR 数字滤波器:幅频特性较好;但相频特性较差; 有稳定性问题;FIR 数字滤波器:可以严格线性相位,又可任意幅度特性因果稳定系统可用FFT 计算(计算两个有限长序列的线性卷积但阶次比IIR 滤波器要高得多3北京邮电大学信息与通信工程学院概述:IIR 和FIR 比较IIR 与FIR 设计方法比较IIR DF :无限冲激响应,H(Z 是z -1的有理分式,借助于模拟滤波器设计方法,阶数低(同样性能要求。
其优异的幅频特性是以非线性相位为代价的。
缺点:只能设计特定类型的滤波器,不能逼近任意的频响。
FIR DF :有限冲激响应,系统函数H(Z 是z -1的多项式,采用直接逼近要求的频率响应。
设计灵活性强缺点:①设计方法复杂;②延迟大;③阶数高。
(运算量比较大,因而在实现上需要比较多的运算单元和存储单元FIR DF 的技术要求:通带频率ωp ,阻带频率ωs 及最大衰减αp ,最小衰减αs 很重要的一条是保证H(z 具有线性相位。
4北京邮电大学信息与通信工程学院概述:FIR DF 设计方法FIR 数字滤波器设计FIR 滤波器的任务:给定要求的频率特性,按一定的最佳逼近准则,选定h(n 及阶数N 。
三种设计方法:n 窗函数加权法o 频率采样法p FIR DF 的CAD --切比雪夫等波纹逼近法5北京邮电大学信息与通信工程学院概述:FIR DF 零极点FIR 滤波器的I/O 关系:10N r y(nh(rx(n r−==−∑0121(, ,,,...,=−h n n N FIR 滤波器的系统传递函数:1211011N N N rN r h(z h(z .....h(N H(zh(rzz −−−−−=++−==∑⇒在Z 平面上有N-1 个零点;在原点处有一个(N-1阶极点,永远稳定。
FIR 系统定义:一个数字滤波器DF 的输出y(n,如果仅取决于有限个过去的输入和现在的输入x(n, x(n-1,. ......, x(n-N+1,则称之为FIR DF 。
dsp实验报告 fir实验报告DSP实验报告:FIR实验报告引言:数字信号处理(Digital Signal Processing,DSP)是一门研究如何对数字信号进行处理和处理的学科。
其中,滤波器是数字信号处理中最常用的技术之一。
本实验报告旨在介绍FIR(Finite Impulse Response)滤波器的原理、设计和实现过程,并通过实验验证其性能。
一、FIR滤波器的原理FIR滤波器是一种线性时不变系统,其输出信号仅由输入信号的有限个历史样本决定。
其基本原理是将输入信号与滤波器的冲激响应进行卷积运算,以实现对输入信号的滤波处理。
二、FIR滤波器的设计方法1. 理想低通滤波器设计方法理想低通滤波器的频率响应在截止频率之前为1,在截止频率之后为0。
通过对理想低通滤波器的频率响应进行采样和离散化,可以得到FIR滤波器的系数序列。
2. 窗函数法设计FIR滤波器窗函数法是一种常用的FIR滤波器设计方法。
其基本思想是将理想低通滤波器的频率响应与一个窗函数进行乘积,从而得到实际可实现的FIR滤波器的系数序列。
常用的窗函数有矩形窗、汉宁窗、汉明窗等。
三、FIR滤波器的实现FIR滤波器可以通过直接形式和间接形式两种方式实现。
直接形式是按照滤波器的差分方程进行计算,而间接形式则是利用FFT(Fast Fourier Transform)算法将滤波器的系数序列转换为频域进行计算。
四、FIR滤波器的性能评估1. 幅频响应幅频响应是评估FIR滤波器性能的重要指标之一。
通过绘制滤波器的幅频响应曲线,可以直观地观察滤波器在不同频率下的衰减情况。
2. 相频响应相频响应是评估FIR滤波器性能的另一个重要指标。
相频响应描述了滤波器对输入信号的相位延迟情况,对于某些应用场景,相频响应的稳定性和线性性非常重要。
3. 稳态误差稳态误差是指FIR滤波器在达到稳态后输出信号与理想输出信号之间的差异。
通过对滤波器的输入信号进行模拟或实际测试,可以计算出滤波器的稳态误差,并评估其性能。
实验四:FIR数字滤波器设计与软件实现1.实验目的(1)掌握用窗函数法设计FIR数字滤波器的原理和方法。
(2)学会调用MA TLAB函数设计FIR滤波器。
(3)通过观察频谱的相位特性曲线,建立线性相位概念。
(4)掌握FIR数字滤波器的MATLAB软件实现方法。
2.实验原理设计FIR数字滤波器一般采用直接法,如窗函数法和频率采样法。
本实验采用窗函数法设计FIR滤波器,要求能根据滤波需求确定滤波器指标参数,并按设计原理编程设计符合要求的FIR数字滤波器。
本实验软件实现是调用MATLAB提供的fftfilt函数对给定输入信号x(n)进行滤波,得到滤波后的输出信号y(n)。
3. 实验内容及步骤(1) FIR数字滤波器设计根据窗函数法设计FIR数字滤波器的原理和步骤,设计一个线性数字低通滤波器,要求通带临界频率fp=120Hz,阻带临界频率fs=150Hz,通带内的最大衰减Ap=0.1dB,阻带内的最小衰减As=60db,采样频率Fs=1000Hz。
观察设计的滤波器频率特性曲线,建立线性相位概念。
(2) FIR数字滤波器软件实现利用第(1)步设计的数字滤波器,调用fftfilt函数对信号进行滤波,观察滤波前后的信号波形变化。
4.思考题(1)简述窗函数法设计FIR数字滤波器的设计步骤。
(2)简述信号在传输过程中失真的可能原因。
5.实验报告要求(1)结合实验内容打印程序清单和信号波形。
(2)对实验结果进行简单分析和解释。
(3)简要回答思考题。
常用窗函数技术参数及性能比较一览表窗类型最小阻带衰减主瓣宽度精确过渡带宽窗函数矩形窗21dB 4π/M 1.8π/M boxcar三角窗25dB 8π/M 6.1π/M bartlett汉宁窗44dB 8π/M 6.2π/M hanning哈明窗53dB 8π/M 6.6π/M hamming 布莱克曼窗74dB 12π/M 11π/M blackman 取凯塞窗时用kaiserord函数来得到长度M和βkaiser附录:(1)FIR数字滤波器设计clear;clc;close all;format compactfp=120, Ap=0.1, fs=150, As=60 ,Fs=1000,wp=2*pi*fp/Fs,ws=2*pi*fs/Fs ,Bt=ws-wp; M=ceil(11*pi/Bt);if mod(M,2)==0; N=M+1, else N=M, end;wc=(wp+ws)/2,n=0:N-1;r=(N-1)/2;hdn=sin(wc*((n-r)+eps))./(pi*((n-r)+eps));win=blackman(N); hn=hdn.*win',figure(1);freqz(hn,1,512,Fs);grid on;图(一)FIR数字滤波器(2)FIR数字滤波器软件实现n=[0:190];xn=sin((2*pi*120/1000)*n)+sin((2*pi*150/1000)*n);yn=fftfilt(hn,xn);figure(2)subplot(2,1,1);plot(xn);title('滤波前信号') ;subplot(2,1,2);plot(yn);title('滤波后信号');图(2)FIR数字滤波器软件实现思考题:(1) 用升余弦窗设计一线性相位低通FIR数字滤波器,并读入窗口长度。
FIR滤波器设计与实现实验报告实验报告:FIR滤波器设计与实现一、实验目的本实验旨在通过设计和实现FIR滤波器来理解数字滤波器的原理和设计过程,并且掌握FIR滤波器的设计方法和实现技巧。
二、实验原理1.选择滤波器的类型和阶数根据滤波器的类型和阶数的不同,可以实现不同的滤波效果。
常见的滤波器类型有低通滤波器、高通滤波器、带通滤波器和带阻滤波器。
选择适当的滤波器类型和阶数可以实现对不同频率分量的滤波。
2.确定滤波器的系数在设计FIR滤波器时,系数的选择对滤波器的性能有重要影响。
通常可以使用窗函数法、最小二乘法、频率采样法等方法来确定系数的值。
常见的窗函数有矩形窗、汉明窗和布莱克曼窗等。
三、实验步骤1.确定滤波器的类型和阶数根据实际需求和信号特点,选择合适的滤波器类型和阶数。
例如,如果需要设计一个低通滤波器,可以选择实验中使用的巴特沃斯低通滤波器。
2.确定滤波器的频率响应根据滤波器的类型和阶数,确定滤波器的频率响应。
可以通过matlab等软件来计算和绘制滤波器的频率响应曲线。
3.确定滤波器的系数根据频率响应的要求,选择合适的窗函数和窗长度来确定滤波器的系数。
可以使用matlab等软件来计算和绘制窗函数的形状和频率响应曲线。
4.实现滤波器的功能将滤波器的系数应用于输入信号,通过加权求和得到输出信号的采样点。
可以使用matlab等软件来模拟和验证滤波器的功能。
四、实验结果在实际实验中,我们选择了一个4阶低通滤波器进行设计和实现。
通过计算和绘制滤波器的频率响应曲线,确定了窗函数的形状和窗长度。
在实际实验中,我们通过实现一个滤波器功能的matlab程序来验证滤波器的性能。
通过输入不同频率和幅度的信号,观察滤波器对信号的影响,验证了设计的滤波器的功能有效性。
五、实验总结通过本实验,我们深入了解了FIR滤波器的设计原理和实现方法。
通过设计和实现一个具体的滤波器,我们掌握了滤波器类型和阶数的选择方法,以及系数的确定方法。
FIR实验报告范文实验报告:FIR滤波器设计与实现一、实验目的1.了解和掌握FIR滤波器的基本概念和设计原理;2. 学会使用Matlab软件设计和实现FIR滤波器;3.分析和评价不同阶数(N)和窗函数类型对FIR滤波器性能的影响。
二、实验设备和原材料1.计算机;2. Matlab 2024a 软件;3.信号发生器。
三、实验原理1.FIR滤波器概念:FIR(Finite Impulse Response)即有限脉冲响应滤波器,是一种非递归滤波器。
其输出是当前输入和过去若干个输入的加权和。
2.FIR滤波器设计原理:FIR滤波器的设计一般通过确定滤波器系数实现。
常用的设计方法有直接设计法、频率采样法、窗函数法等。
其中,窗函数法是最常用的设计方法之一窗函数法步骤:a.选择适当的滤波器阶数N;b.选择合适的窗函数;c.根据滤波器的频率响应要求,计算出窗函数长度L;d.确定滤波器的理想频率响应;e.执行窗函数和理想频率响应的点乘运算得到滤波器的系数。
四、实验步骤step1. 设计滤波器频率响应:a. 确定采样频率fs和截止频率fc;b. 根据fs和fc确定滤波器通带、阻带宽度;c.设计理想频率响应,生成指定的通带增益和阻带衰减。
step2. 确定滤波器阶数和窗函数类型:a.根据滤波器的过渡带宽度、通带纹波和阻带衰减要求,选择适当的滤波器阶数N;b. 对于给定的滤波器长度L和滤波器阶数N,结合通带和阻带宽度,借助Matlab软件选择合适的窗函数类型。
step3. 确定窗函数长度L:根据滤波器的过渡带宽度,计算滤波器的窗函数长度L。
一般,L=N+1step4. 生成滤波器系数:通过窗函数与理想频率响应的点乘运算,生成滤波器的系数。
step5. 滤波器模拟与实现:a. 在Matlab软件中,使用fir1函数生成滤波器系数h;b. 使用filter函数实现滤波器的模拟和实现。
五、实验结果与分析本实验选择了截止频率fc = 1000Hz,采样频率fs = 8000Hz。
FIR 滤波器设计与实现姓名__孔亮__ 学号__0928401116__ 专业__电子信息工程__一、实验目的(1) 通过实验巩固FIR 滤波器的认识和理解。
(2) 熟练掌握FIR 低通滤波器的窗函数设计方法。
(3) 理解FIR 的具体应用。
二、实验内容在通信、信息处理以及信号检测等应用领域广泛使用滤波器进行去噪和信号的增强。
FIR 滤波器由于可实现线性相位特性以及固有的稳定特征而等到广泛应用,其典型的设计方法是窗函数设计法。
设计流程如下:(1)设定指标:截止频率fc ,过渡带宽度△f ,阻带衰减A 。
(2)求理想低通滤波器(LPF )的时域响应h d (n )。
(3)选择窗函数w (n ),确定窗长N 。
(4)将h d (n )右移(N-1)/2点并加窗获取线性相位FIR 滤波器的单位脉冲响应h (n )。
(5)求FIR 的频域响应H (e jw ),分析是否满足指标。
如不满足,转(3)重新选择,否则继续。
(6)求FIR 的系统函数H (z )。
(7)依据差分方程由软件实现FIR 滤波器或依据系统函数由硬件实现。
实验要求采用哈明窗设计一个FIR 低通滤波器并由软件实现。
哈明窗函数如下:w (n )=0.54-0.46cos (12-N nπ),0≤n ≤N-1;设采样频率为fs=10kHz 。
实验中,窗长度N 和截止频率fc 应该都能调节。
具体实验内容如下:(1)设计FIR 低通滤波器(FIR_LPF )(书面进行)。
(2)依据差分方程编程实现FIR 低通滤波器。
(3)输入信号x (n )=3.0sin (0.16n π)+cos (0.8n π)到fc=2000Hz ,N=65的FIR_LPF ,求输出信号y (n ),理论计算并画出0≤f ≤fs 范围输入信号x (n )和输出信号y (n )的幅度谱,标出峰值频率,观察滤波器的实际输出结果,分析其正确性。
(4)输入信号x (n )=1.5sin (0.2n π)-cos (0.4n π)+1.2sin (0.9n π)到fc=1100Hz ,N=65的FIR_LPF ,求输出信号y (n ),理论计算并画出0≤f ≤fs 范围输入信号x (n )和输出信号y (n )的幅度谱,标出峰值频率,观察滤波器的实际输出结果,分析其正确性。
(5)输入信号x (n )=1.5sin (0.2n π)-cos (0.4n π)+1.2sin (0.9n π)到fc=2100Hz ,N=65的FIR_LPF ,求输出信号y (n ),理论计算并画出0≤f ≤fs 范围输入信号x (n )和输出信号y (n )的幅度谱,标出峰值频率,观察滤波器的实际输出结果,分析其正确性。
(6)输入信号x (n )=1.5sin (0.2n π)-cos (0.4n π)+5.0sin (0.9n π)到fc=1100Hz ,N=65的FIR_LPF ,求输出信号y (n ),理论计算并画出0≤f ≤fs 范围输入信号x (n )和输出信号y (n )的幅度谱,标出峰值频率,观察滤波器的实际输出结果,分析其正确性。
(7)输入信号x (n )=1.5sin (0.2n π)-cos (0.4n π)+1.2sin (0.9n π)到fc=1990Hz ,N=65的FIR_LPF ,求输出信号y (n ),理论计算并画出0≤f ≤fs 范围输入信号x (n )和输出信号y (n )的幅度谱,标出峰值频率,观察滤波器的实际输出结果,分析其正确性。
三、实验分析如果采用哈明窗函数设计因果线性相位低通滤波器FIR ,所设计的FIR 低通滤波器单位脉冲响应h (n )如式。
h (n )=h d (n )w (n )=)]12cos(46.054.0[)21()]21(sin[------N n N n N n wc ππ,0≤n ≤N-1如果输入函数为x (n ),则y (n )=x (n )*h(n)。
这里要注意N 的取值,若取值为奇数,则会导致分母为零,频谱出错,故若N 为奇数,则取N+1。
因为ω=2πf/fs ,所以有2πf=ωfs 。
以此为依据对输入信号作如下理论分析:(1) 输入信号x(n)=3.0sin(0.16πn)+cos(0.8πn)到fc=2000Hz,N=65的FIR_LPF 时,由幅度值值峰值频率由正弦分量决定f=ωfs/2π=800Hz,输入信号x(n)的频率分量有:f1=800Hz ,f2=4000Hz ,进过低通滤波器后变成输出信号y(n),峰值频率应该为fm=800Hz ;(2)输入信号x(n)=1.5sin(0.2πn)-cos(0.4πn)+ 1.2sin(0.9πn)到fc=1100Hz,N=65的FIR_LPF 时,由幅度值值峰值频率由正弦分量决定f=ωfs/2π=1000Hz,输入信号x(n) 的频率分量有:f1=1000Hz ,f2=2000Hz ,f3=4500Hz ,进过低通滤波器后变成输出信号y(n),所以峰值频率应该为fm=1000Hz ;(3)输入信号x(n)=1.5sin(0.2πn)-cos(0.4πn)+ 1.2sin(0.9πn)到fc=2100Hz,N=65的FIR_LPF 时,由幅度值值峰值频率由正弦分量决定f=ωfs/2π=1000Hz,输入信号x(n) 的频率分量有:f1=1000Hz ,f2=2000Hz ,f3=4500Hz ,进过低通滤波器后变成输出信号y(n),所以峰值频率应该为fm=1000Hz ;(4)输入信号x(n)=1.5sin(0.2πn)-cos(0.4πn)+ 5.0sin(0.9πn)到fc=1100Hz,N=65的FIR_LPF 时,由幅度值值峰值频率由正弦分量决定f=ωfs/2π=4500Hz,输入信号x(n) 的频率分量有:f1=1000Hz ,f2=2000Hz ,f3=4500Hz ,进过低通滤波器后变成输出信号y(n),所以峰值频率还是应该在1000Hz ; (5)输入信号x(n)=1.5sin(0.2πn)-cos(0.4πn)+ 1.2sin(0.9πn)到fc=1990Hz,N=65的FIR_LPF 时,由幅度值值峰值频率由正弦分量决定f=ωfs/2π=1000Hz,输入信号x(n) 的频率分量有:f1=1000Hz ,f2=2000Hz ,f3=4500Hz ,进过低通滤波器后变成输出信号y(n),所以峰值频率应该1000Hz ;四、实验结果(1)MATLAB 程序如下:function fir(A,a,B,b,C,c,fc)N=66;wc=fc*2*3.14/10000;for i=1:Nx(i)=A*sin(a*3.14*i)+B*cos(b*3.14*i)+C*sin(c*3.14*i);h(i)=sin(wc*(i-(N-1)/2))/(3.14*(i-(N-1)/2))*(0.54-0.46*cos(2*3.14*i/( N-1)));y =conv(x,h);endsubplot(4,1,1);plot((abs(fft(x))),'-*')subplot(4,1,2);plot((abs(fft(y))),'-*')subplot(4,1,3);plot((abs(fft(h))),'-*')disp(max((abs(fft(x)))));disp(max((abs(fft(y)))));Y=(abs(fft(y)));X=(abs(fft(x)));maxim=0;maxi=0;maxk=0;maxl=0;for l=1:Nif X(l)>maximaxi=X(l);maxl=l;endfor k=1:Nif Y(k)>maximmaxim=Y(k);maxk=k;endendenddisp((maxl-1)/N*10000)disp(maxl)disp((maxk-1)/N/2*10000)disp(maxk)subplot(4,1,4);plot(y,'-*')π)+cos(0.8nπ),fc=2000Hz(3) x(n)=3.0sin(0.16n第一行为X(e jw),第二行为Y(e jw),第三行为H(e jw),第四行为y(n)当f=757.5758Hz时,峰值为90.9039。
与理论值相接近,可认为系统满足要求。
第一行为X(e jw),第二行为Y(e jw),第三行为H(e jw),第四行为y(n)当f=984.8485Hz时,峰值为42.3833。
与理论值相接近,可认为系统满足要求。
第一行为X(e jw),第二行为Y(e jw),第三行为H(e jw),第四行为y(n)当f=984.8485Hz时,峰值为50.9443。
与理论值相接近,可认为系统满足要求。
第一行为X(e jw),第二行为Y(e jw),第三行为H(e jw),第四行为y(n)当f=984.8485Hz时,峰值为43.5255。
与理论值相接近,可认为系统满足要求。
第一行为X(e jw),第二行为Y(e jw),第三行为H(e jw),第四行为y(n)当f=984.8485Hz时,峰值为50.9352。
与理论值相接近,可认为系统满足要求。
五、思考题1、当哈明窗长度N比65小(32)或大(129)的话,实验结果如何变化?答:(1)当N=32时A:x(n)=3.0sin(0.16nπ)+cos(0.8nπ),fc=2000HzB:x(n)=1.5sin(0.2nπ)-cos(0.4nπ)+1.2sin(0.9nπ),fc=1100HzC:x(n)=1.5sin(0.2nπ)-cos(0.4nπ)+1.2sin(0.9nπ),fc=2100D:x(n)=1.5sin(0.2nπ)-cos(0.4nπ)+5.0sin(0.9nπ),fc=1100Hz(2)当N=129时A:x(n)=3.0sin(0.16nπ)+cos(0.8nπ),fc=2000HzC:x(n)=1.5sin(0.2nπ)-cos(0.4nπ)+1.2sin(0.9nπ),fc=2100HzE、x(n)=1.5sin(0.2nπ)-cos(0.4nπ)+1.2sin(0.9nπ),fc=1990Hz2、当采用矩形窗的话,实验(3)(4)的结果是怎样的?答:矩形窗设计的FIR 低通滤波器的单位脉冲响应为:h (n )=hd (n )w (n )=)21()]21(sin[----N n N n wc π ,0≤n ≤N-1 (A )x (n )=3.0sin (0.16n π)+cos (0.8n π)(B )x (n )=1.5sin(0.2n π)-cos(0.4n π)+1.2sin(0.9n π)3、(6)的实验结果说明什么?答:(6)的实验结果表明该低通滤波器系统能够有效的滤除高于截止频率的分量。