当前位置：文档之家 › FIR滤波器设计与实现

FIR滤波器设计与实现

  • 格式：doc
  • 大小：672.00 KB
  • 文档页数：19

实验二 FIR 滤波器设计与实现10通信： 1028401031 沈欢健 1028401030 龚俊 1028401025 鲍也 一、实验目的（1） 通过实验巩固FIR 滤波器的认识和理解。

（2） 熟练掌握FIR 低通滤波器的窗函数设计方法。

（3） 理解FIR 的具体应用。

二、实验内容 在通信、信息处理以及信号检测等应用领域广泛使用滤波器进行去噪和信号的增强。

FIR 滤波器由于可实现线性相位特性以及固有的稳定特征而等到广泛应用，其典型的设计方法是窗函数设计法。

设计流程如下：（1）设定指标：截止频率fc ，过渡带宽度△f ，阻带衰减A 。

（2）求理想低通滤波器（LPF ）的时域响应hd （n ）。

（3）选择窗函数w （n ），确定窗长N 。

（4）将hd （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 ）=hd （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)。

四、实验结果（1）MATLAB 程序如下：function FIR = fir(fc,N,A,a,B,b,C,c) %窗函数设计n = 0 : N-1; %窗函数、LPF 的横坐标 fs = 10000; %采样频率10KHZ wc = 2 * pi * fc/fs; hd = n;for k = 1 : (N-1)/2hd(k) = sin(wc * (k-1-(N-1)/2)) ./ ( pi *(k-1-(N-1)/2)); %LPF 纵坐标 endfor k =(N-1)/2 + 2 : Nhd(k) = sin(wc * (k-1-(N-1)/2)) ./ ( pi *(k-1-(N-1)/2)); %LPF 纵坐标 endhd((N-1)/2+1) = wc/pi;wn = 0.54 - 0.46*cos(2*pi*n/(N-1)); %哈明窗窗函数 hn = hd .* wn; %FIR 滤波器时域表达式%窗函数频谱w = 0 : 0.001*pi :0.999* pi; W = 0:999;for k = 1 : 1000;W(k) = 0; n = 1;while(n < N+1)W(k) = hd(n) * exp(-j*k*2*pi/1000*n)+ W(k); n = n + 1;endend%输入频谱n = 0 : 999;x = A * sin ( a * pi * n ) + B * cos ( b * pi * n) + C* sin ( c * pi * n );X = x;s = 0;for k = 1 : 1000X(k) = 0;n = 1;while( n < 1001)X(k) = X(k) + x(n) * exp(-j * 2* pi *(k-1)*(n-1)/1000);n = n + 1;endif k < 500if s < abs(X(k))& (s ~= abs(X(k)))s = abs(X(k));m = k - 1;endendendY = X .* W;t = 0 : N-1; %过渡变量k = 0 : 0.002 * pi : 1.998 *pi;subplot(2,2,1); %图1plot(t,hn); %窗函数时域图subplot(2,2,2); %图2plot(k,20*log10(abs(W))); %窗函数频谱图l=0:999;subplot(2,2,3); %图3plot(l/1000*fs,abs(X),'.-'); %输入函数频谱图subplot(2,2,4);plot(k/1000*fs,abs(Y),'.-');f = m /1000 * fssendπ）+cos（0.8nπ），截止频率fc=2000Hz，(3)当输入信号为 x（n）=3.0sin（0.16nN=65时理论分析：根据公式 w=2*π*f/fs得到此时的输入信号的频率为800Hz和4000Hz，由滤波器滤波后得到的频率为800Hz，把4000Hz的频率滤除了。

实验仿真结果：滤波器单位脉冲响h（n）应如下图（a）所示，滤波器频谱响应如图（b）所示，输入信号的频谱如图（c）所示，输出信号的频谱如图（d）所示图（a）图(b)如图（b）所示哈明窗的阻带衰减为-53dB与理论-54dB相近，此时的截止角频率为0.4π如图（c）所示输入信号的频率为800Hz和4000Hz图（d）如图（d）所示输出信号的频率为800Hz，滤波器把4000Hz的频率成分滤除了，实验得到的结果与理论相一致。

π)-cos(0.4nπ)+1.2sin(0.9nπ)截止频率fc=1100Hz， (4) x（n）=1.5sin(0.2nN=65时理论分析：根据公式 w=2*π*f/fs得到此时的输入信号的频率为1000Hz，2000Hz 和4500Hz，由滤波器滤波后得到的频率为1000Hz，把4500Hz和2000Hz的频率滤除了。

实验仿真结果：滤波器单位脉冲响h（n）应如下图（a）所示，滤波器频谱响应如图（b）所示，输入信号的频谱如图（c）所示，输出信号的频谱如图（d）所示图（a）图（b）如图（b）所示哈明窗的阻带衰减为-54.39dB与理论-54dB相近，此时的截止角频率为0.22π图（c）如图（c）所示输入信号的频率为1000Hz，4500Hz和2000Hz图（d）如图（d）所示输出信号的频率为1000Hz，滤波器把4500Hz和2000Hz的频率成分滤除了，实验得到的结果与理论相一致。

π)-cos(0.4nπ)+1.2sin(0.9nπ)截止频率fc=2100Hz， (5) x（n）=1.5sin(0.2nN=65时理论分析：根据公式 w=2*π*f/fs得到此时的输入信号的频率为1000Hz，2000Hz 和4500Hz，由滤波器滤波后得到的频率为1000Hz和2100，把4500Hz的频率滤除了。

实验仿真结果：滤波器单位脉冲响h（n）应如下图（a）所示，滤波器频谱响应如图（b）所示，输入信号的频谱如图（c）所示，输出信号的频谱如图（d）所示图（a）图（b）如图（b）所示哈明窗的阻带衰减为-55dB与理论-54dB相近，此时的截止角频率为0.42π如图（c）所示输入信号的频率为1000Hz，4500Hz和2000Hz滤除了，实验得到的结果与理论相一致。

π)-cos(0.4nπ)+5.0sin(0.9nπ)截止频率fc=1100Hz，(6) x（n）=1.5sin(0.2nN=65时理论分析：根据公式 w=2*π*f/fs得到此时的输入信号的频率为1000Hz，2000Hz 和4500Hz，由滤波器滤波后得到的频率为1000Hz，把4500Hz和2100Hz的频率滤除了。

实验仿真结果：滤波器单位脉冲响h（n）应如下图（a）所示，滤波器频谱响应如图（b）所示，输入信号的频谱如图（c）所示，输出信号的频谱如图（d）所示图（a）图（b）如图（b）所示哈明窗的阻带衰减为-55dB与理论-54dB相近，此时的截止角频率为0.22π图（c）如图（c）所示输入信号的频率为1000Hz，4500Hz和2000Hz图（d）如图（d）所示输出信号的频率为1000Hz，滤波器把4500Hz和2000Hz的频率成分滤除了，实验得到的结果与理论相一致。

π)-cos(0.4nπ)+1.2sin(0.9nπ)截止频率fc=1990Hz，(7) x（n）=1.5sin(0.2nN=65时理论分析：根据公式 w=2*π*f/fs得到此时的输入信号的频率为1000Hz，2000Hz 和4500Hz，由滤波器滤波后得到的频率为1000Hz，把4500Hz和2100Hz的频率滤除了。

合集下载

数字信号处理实验报告-FIR滤波器的设计与实现

数字信号处理实验报告-FIR滤波器的设计与实现

数字信号处理实验报告-FIR滤波器的设计与实现在数字信号处理中，滤波技术被广泛应用于时域处理和频率域处理中，其作用是将设计信号减弱或抑制被一些不需要的信号。

根据滤波器的非线性抑制特性，基于FIR（Finite Impulse Response）滤波器的优点是稳定，易设计，可以得到较强的抑制滤波效果。

本实验分别通过MATLAB编程设计、实现、仿真以及分析了一阶低通滤波器和平坦通带滤波器。

实验步骤：第一步：设计一阶低通滤波器，通过此滤波器对波型进行滤波处理，分析其对各种频率成分的抑制效果。

为此，采用零极点线性相关算法设计滤波器，根据低通滤波器的特性，设计的低通滤波器的阶次为n=10，截止频率为0.2π，可以使设计的滤波器被称为一阶低通滤波器。

第二步：设计平坦通带滤波器。

仿真证明，采用兩個FIR濾波器組合而成的阻礙-提升系統可以實現自定義的總三值響應的設計，得到了自定義的總三值響應函數。

实验结果：1、通过MATLAB编程，设计完成了一阶低通滤波器，并通过实验仿真得到了一阶低通滤波器的频率响应曲线，证明了设计的滤波器具有良好的低通性能，截止频率为0.2π。

在该频率以下，可以有效抑制波形上的噪声。

2、设计完成平坦通带滤波器，同样分析其频率响应曲线。

从实验结果可以看出，此滤波器在此频率段内的通带性能良好，通带范围内的信号透过滤波器后，损耗较小，滞后较小，可以满足各种实际要求。

结论：本实验经过实验操作，设计的一阶低通滤波器和平坦通带滤波器具有良好的滤波特性，均已达到预期的设计目标，证明了利用非线性抑制特性实现FIR滤波处理具有较强的抑制滤波效果。

本实验既有助于深入理解FIR滤波器的设计原理，也为其他应用系统的设计和开发提供了指导，进而提高信号的处理水平和质量。

fir数字滤波器设计与软件实现数字信号处理实验原理

fir数字滤波器设计与软件实现数字信号处理实验原理

fir数字滤波器设计与软件实现数字信号处理实验原理FIR数字滤波器设计的基本原理是从理想滤波器的频率响应出发，寻找一个系统函数，使其频率响应尽可能逼近滤波器要求的理想频率响应。

为了实现这一目标，通常会采用窗函数法进行设计。

这种方法的基本思想是，将理想滤波器的无限长单位脉冲响应截断为有限长因果序列，并用合适的窗函数进行加权，从而得到FIR滤波器的单位脉冲响应。

在选择窗函数时，需要考虑其频率响应和幅度响应。

常见的窗函数包括矩形窗、三角形窗、汉宁窗、汉明窗、布莱克曼窗和凯泽窗等。

每种窗函数都有其特定的特性，如主瓣宽度、旁瓣衰减等。

根据实际需求，可以选择合适的窗函数以优化滤波器的性能。

在软件实现上，可以使用各种编程语言和信号处理库进行FIR滤波器的设计和实现。

例如，在MATLAB中，可以使用内置的`fir1`函数来设计FIR滤波器。

该函数可以根据指定的滤波器长度N和采样频率Fs，自动选择合适的窗函数并计算滤波器的系数。

然后，可以使用快速卷积函数`fftfilt`对输入信号进行滤波处理。

此外，还可以使用等波纹最佳逼近法来设计FIR数字滤波器。

这种方法的目标是找到一个最接近理想滤波器频率响应的实数序列，使得在所有可能的实
数序列中，该序列的误差平方和最小。

通过优化算法，可以找到这个最优序列，从而得到性能更优的FIR滤波器。

总的来说，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滤波器的MATLAB设计与实现

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滤波器对音频信号进行去噪处理，或者对图像进行锐化处理等。