7电子0501郭标.docVIP

西安院It is desired to extract a constant signal x (n)=s(n)+ν(n)=s+ ν (n),where ν(n) is zero-mean white Caussian noise of varianceбv2 . (1) First-order IIR smoother; To this end, the following IIR lowpass filter is used:, where the parameter a is restrict to the range 0a1.This can be achieve in two ways. (2) FIR averaging filters: Consider ,for example, the third-order filter H(z)=h0+h1z-1+h2z-2+h3z-3 2 实验目的： （1）.运用MATLAB的平台编辑和设计程序，熟悉MATLAB语言及其调试方法。 （2）进一步加深对卷积定理和采样定理的理解。 （3）了解滤波器的工作原理和设计方法。 （4）学习用滤波器的方法把信号从噪声中滤出来。 （5）巩固数字信号处理的基本知识，能够运用MATLAB知识解决数字信号中的问题。 （6）让确定信号分别通过有限长脉冲响应（FIR）数字滤波器和无限长脉冲响应（IIR）数字滤器，通过分析输出图像，掌握滤波器各参数对性能的影响。 3 实验原理： 由于待提取的分量为直流分量，而干扰信号的频谱占据整个频带，本题将从零点、极点分布特性出发，利用所谓的零、极点累试方法进行滤波器的设计。滤波器的特性取决于零、极点的情况。极点位置主要影响系统的幅度特性峰值位置及尖锐程度，零点位置主要影响系统的幅度特性谷值位置及其凹下去的程度；且通过零、极点的反复调整，并将得到的幅度特性与期望的幅度特性惊醒比较，最终得到所期望系统的幅度特性。时域中y(n)=x(n)*h(n)，而频域中Y(Z)=X(Z)·H(Z)，而对于随机信号而言，其频谱复杂，因此利用时域卷积来计算较为简单。 4 实验步骤：: 生成直流信号 生成噪声信号 确定系数 对h和x进行卷积。 输出图形。 5 程序框图：: 6 源程序： close all; clear all; n=1:200; N=200; s=5*ones(1,N); v=randn(1,N); x=v+s; a1=0.97; a2=0.03; b1=1-a1; b2=1-a2; u(n)=ones(1,N); h1(n)=b1*(a1.^n).*u(n); y1=conv(x,h1(n)); h2(n)=b2*(a2.^n).*u(n); y2=conv(x,h2(n)); h3=[0.25,0.25,0.25,0.25]; h4=[0.2,0.02,0.5,0.01]; y3=conv(x,h3); y4=conv(x,h4); figure(1); subplot(2,1,1),plot(y1,r-);xlabel( a趋近于1) title(通过IIR滤波器) subplot(2,1,2),plot(y2,r-);xlabel(a趋近于0) figure(2); subplot(3,1,1),plot(y3,b-);xlabel(四个系数取值相等) title(通过FIR滤波器) subplot(3,1,2),plot(y4,b-);xlabel(四个系数取值不等) subplot(3,1,3),plot(x,b-);xlabel(叠加噪声后的波形) 实验结果与实验结果分析： 上图为.a=0.97与a=0.03时的滤波效果图时。通过比较可知，极点的位置越靠近单位圆（即参数a越大），在w=0处其频率特性越陡峭，则对噪声的滤波效果越明显，故其滤波效果越好。 上图为h0=h1=h2=h3=0.25四个系数相等时与h0=0.2，h1=0.02，h2=0.5，h3=0.01四个系数不等时的滤波效果图。通过比较可得当h0=h1=h2=h3=0.25四个系数相等时，其噪声抑制效果最佳。 第五题 题目五：DTMF信号的生成是通过将两个有限长数字序列相加而实现；而对DTMF信号的检测是通过计算DTMF信号的DFT；然后测量在给定8个频率上的能量而实现。已知采样频率为8KHz，DTMF信号点数N＝205，对DTMF信号进行205点DFT。右表为DTMF数字。 1209Hz 1336Hz 1447Hz 1663Hz 1 697Hz 1 2 3 A 2 770Hz 4 5 6 B 3 852H