信号实验报告(三).docx

课程：测试信号处理技术实验名称：信号分析处理综合实验图片目录图片 1 加入两正弦信号后的频谱4图片 2 原型及滤波后频谱5图片 3 局部放大信号5图片 4 设计的IIR低通滤波器6图片 5 设计的FIR低通滤波器6图片 6 实验程序截图7图片 7 信号频谱对比7实验三：信号分析处理综合实验（信号叠加、分析及信号滤波频谱）实验目的综合应用信号频谱分析和数字滤波器设计的知识，实现信号的滤波。比较分析干扰信号加入前、后的频谱有何区别，分析噪声的影响。加深理解信号时域和频域分析的物理概念，认识不同滤波器的特性和适用范围。掌握设计IIR、FIR 数字滤波器的原理和特性，通过比较对两种滤波器的滤波作用有直观认识。实验原理及方法（1） FIR 数字滤波器总是稳定的系统，且可以设计成具有线性相位。M 阶FIR 数字滤波器的系统函数为系统的单位脉冲响应h[k]是长度为M+1的有限长因果序列。当满足h[k] = ±h[M ? k]的对称条件时，该FIR 数字滤波器具有线性相位。FIR 数字滤波器设计方法主要有窗口法、频率取样法、优化设计法。MATLAB 中常用FIR 数字滤波器设计函数：fir1 窗函数法FIR 数字滤波器设计：低通、高通、带通、带阻、多频带滤波器。fir2 频率抽样法FIR 数字滤波器设计：任意频率响应。firls FIR 数字滤波器设计：指定频率响应。（2） IIR 数字滤波器的设计原理：IIR 数字滤波器一般为线性时不变的因果离散系统，N 阶IIR 数字滤波器的系统函数可以表达为z?1 的有理多项式，即系数中至少有一个非零。对于因果IIR 数字滤波器，满足M ≤ N 。IIR 数字滤波器的设计主要通过成熟的模拟滤波器设计方法来实现。首先在频域将数字滤波器设计指标转化为模拟低通滤波器设计指标，将任意模拟滤波器转换为原型模拟低通滤波器设计指标，根据模拟低通滤波器设计指标设计出模拟低通滤波器，由经过相应的复频域变换得到H(s)，由H(s)经过脉冲响应不变法或双线性变换法得到所需的IIR 数字滤波器H(z)，由此可见，IIR 数字滤波器的设计的重要环节是模拟低通滤波器的设计，涉及低通模拟滤波器的方法哟Butterworth、Chebyshev 和椭圆等滤波器设计方法。用 MATLAB 设计IIR 滤波器的常用函数：IIR 滤波器阶数选择：Buttord-------巴特沃斯（Butterworth）滤波器阶数选择。Cheblord-----切比雪夫（Chebyshev）I 型滤波器阶数选择。Cheb2ord----切比雪夫（Chebyshev）II 型滤波器阶数选择。Ellpord-------椭圆（Elliptic）滤波器阶数选择。IIR 滤波器设计：Butter-------巴特沃斯（Butterworth）滤波器设计。cheby1------切比雪夫（Chebyshev）I 型滤波器设计。cheby2------切比雪夫（Chebyshev）II 型滤波器设计。maxflat------通用的巴特沃斯（Butterworth）低通滤波器设计。运行结果由matlab语句y(k)=10*x1(k)+sin(f1*2*pi*k*dt)+sin(f2*2*pi*k*dt);即为如下图所示的蓝色信号；s=sin(500*2*pi*t)+sin(1000*2*pi*t)即为如下图所示的红色干扰信号；图片 1 加入两正弦信号后的频谱将加入两个不同频率的正弦干扰后的信号进行快速傅里叶变换（fft），频谱如下图中最上面的一幅图像；设计IIR滤波器，原型为切比雪夫1型，设计matlab函数为语句为：[b1,a1]=cheby1(3,0.5,wn1);使第一个信号通过该设计的滤波器，然后进行fft，得到频谱图下图中间一幅图像，可以看到：在f=500hz处幅度被减弱，更高频率的幅度被一直，频谱保留了低频部分;对应下图中的最下面一幅图像，此时为使用窗函数法设计的FIR低通滤波器，经滤波及fft处理后的原信号频率即如该图所示；图片 2 原型及滤波后频谱（3） 将两正弦干扰信号及经处理后的原信号放大后，如下图所示，其中蓝色为原信号，红色为正弦干扰信号；图片 3 局部放大信号（4） 将设计的IIR滤波器幅频及相频特性输出，结果如下图所示：图片 4 设计的IIR低通滤波器（5）将设计的FIR滤波器幅频及相频特性输出，结果如下图所示：图片 5 设计的FIR低通滤波器图片6实验程序截图分析及说明通过在程序最末段加入文件输出语句：wavwrite(y1,fs,sole);wavwrite(y2,fs,sole1);可以获得经滤波处理后的音频文件，通过试听可以发现，经过滤波的音频大致上可以听到源文件的节奏，但声音明显嘈杂，细节部分（如高音）几乎