数字信号处理三级项目报告分析.docxVIP

数字信号处理三级项目报告Report of Discussion on Digital Signal Processing摘要本三级项目要求使用C或C++利用数字信号处理知识设计并实现一个具有信号采集、信号时域及频域分析、FIR滤波器设计、数字信号滤波等功能的多频窄带数字信号处理软件仿真系统。通过课程研究项目的实施，加深对数字信号处理基础知识理解的基础上，初步具备运用所学知识进行信号处理的能力，将学生的计算机操作能力、分析能力、软件设计能力与应用实践结合起来，引导学生由浅入深地掌握数字信号处理理论与开发工具，初步具备实际应用的信号处理软件开发与制作基础。目录摘要I1系统整体设计思想12连续信号时域采集模拟实现模块12.1将模拟多频窄带信号理想采样转换为数字信号12.2采样过程13输入信号的频域分析模块13.1DFT对信号进行谱分析13.2信号的时域分析23.3信号的频域分析23.4DFT的高分辨率频谱与高密度频谱之间的区别23.5 谱分析的误差及改进方法34FIR滤波器设计模块4设计方法如下：44.1线性相位法44.1.1第一类线性相位：44.1.2线性相位FIR数字滤器的零点分布特点64.2理想低通滤波器的逼近64.3窗函数设计法原理74.3.1窗函数法设计FIR滤波器的步骤84.3.2加窗处理对矩形频率响应的影响85输入信号的滤波处理及结果分析模块9系统整体设计思想a)信号采集：将模拟单频或多频窄带 信号（可由多个正弦信号叠加）模拟信号理想采样转换为数字信号 b)信号的时域波形：观察输入信号的时域波型。c) 信号频域分析：采用DFT或FFT对信号进行频谱分析，观察信号的频谱分布特性；观察DFT的高分辨率频谱与高密度频谱之间的区别。d) 设计FIR滤波器：利用窗函数方法设计FIR滤波器，观察不同窗函数的选取对滤波器性能的影响。e) 利用FIR滤波器对输入信号进行滤波：观察采用不同窗函数的滤波器对输入信号进行滤波处理后的时域、频域波形。连续信号时域采集模拟实现模块将模拟多频窄带信号理想采样转换为数字信号模拟信号：Xa(t)=cos(Ω1t)+cos(Ω2t)+cos(Ω3t) 其中f1=201Hz , f2=208Hz , f3=214Hz采样频率fs=1000Hz , 观测时间Tp=0.1s 采样过程Xa(t)=xa(t)p(t),其中p(t)为周期冲激脉冲，则有：X(n)即为连续xa(t)以采样频率fs=1/T时得到的离散信号，当fs≥2fmax时，采样信号可以无失真还原原模拟信号，其中fmax为模拟信号的最高频率。输入信号的频域分析模块DFT对信号进行谱分析DFT变换要求时域和频域均离散化且为有限长，所以对连续信号用离散傅里叶进行谱分析时，需要有一个逼近变换过程。DFT进行谱分析时的几个重要参数分析：要增加信号的最高频率，则必须增大，N一定时，F必须增大，那么分辨率就会下降。要提高分辨率则F减小，观测时间就会增大，当N一定时，则减小，且不能产生混叠，则也得减小。同时提高信号的最高频率和频率分辨率需要增加采样点N。那么N增大后工作量大大增加。信号的时域分析信号采集f=201HZ,208HZ,214HZ的100点时域图信号的频域分析采样64点，64点FFTDFT的高分辨率频谱与高密度频谱之间的区别采样64点，64点FFT采样64，补零到128,128点FFT采样128,128点FFTDFT的频谱分辨率：指对信号中两个靠的较近的频谱分量的识别能力。它仅取决于截取连续信号的长度，在采样频率不变时，通过改变采样点数N可以改变DFT的分辨率。DFT的高密度频谱：指当信号的时间域长度不变时，在频域内对它的频谱进行提高采样频率，而得到高密度谱，它只可以更细化当前分辨率下的频谱，但不能改变DFT的分辨率。3.5 谱分析的误差及改进方法1）频谱混叠在对模拟信号 xa(t) 进行采样时，必须满足采样定理，即采样频率 fs≥2fc，而通过上面分析时域有限的信号不可能是锐截止的，并且信号中不可避免地有一些高频杂散信号，因此在采样之前，一般都要对模拟信号进行滤波，滤除高频杂散信号。2）截断效应对无限长的模拟信号，用DFT进行谱分析时，用DFT对连续信号谱分析的误差问题先进行截断，通过采样才能得到有限点的序列，这样必然产生误差。截断可以理解为加窗，即：y(n)=x(n)w(n) 式中：x(n)为模拟信号经采样得到的时域离散信号；w(n)为窗函数序列。根据频域卷积定理，加窗后的信号频谱为：Y(ejω)=1/2πX(ejω)*W(ejω)显然与原序列的频谱是不同的。3） 栅栏效应栅栏效应是指用DFT对连续信号进行谱分析时，由于DFT的离散特性，使离散点之间的频谱无法得到，相当于透过栅栏观察频谱，只看到等间隔的离散点的频线，其他的频谱都被栅栏挡住