Modern Digital Signal Processing curriculm project数字信号处理课程设计.docVIP

目 录 课程设计目的………………………………………….1页 设计任务以及内容………………………………….1-3页 设计方案……………………………………………….4页 设计结果分析讨论…………………………………4-18页 代码……………………………………………….18-25页 心得体会………………………………………….25-26页 参考资料………………………………………………26页 附录一（主界面图）…………………………………….27页 课程设计目的 现代信号处理是一门以算法为核心的，理论和实践性较强的学科。本次课程设计目的主要有：全面复习课程所学知识提高综合运用所学知识独立分析和解决问题的能力；1、单位抽样序列 2）矩形序列 4）反三角波序列 5）Gaussian序列 6） 正弦序列 ?? ? 进行采样，可得到测试序列 令A=50,采样周期T=1ms,即fs=1000Hz,f0=62.5,a=100。 2. 对上述信号完成下列信号分析 1）对三角波序列和反三角波序列，作N=8点的FFT，观察比较它们的幅频特性，说明它们有什么异同？绘出两序列及其它们的幅频特性曲线。 在和的尾部补零，作N=16点的FFT，观察它们的幅频特性发生了什么变化？分析说明原因。 ??2)、观察高斯序列，固定信号中的参数p=8，令q分别等于2，4，8，观察它们的时域和幅频特性，了解当q取不同值时，对信号序列的时域幅频特性的影响；固定q=8，令p分别等于8，13，14，观察参数p变化对信号序列的时域及幅频特性的影响，观察p等于多少时，会发生明显的泄漏现象，混叠是否也随之出现？记录实验中观察到的现象，绘出相应的时域序列和幅频特性曲线。 3）对于正弦序列 4)、观察衰减正弦序列的时域和幅频特性，绘出幅频特性曲线，改变采样频率fs,使fs=300Hz， 观察此时的频谱的形状和谱峰出现位置？说明产生现象的原因。 3．设有一连续时间信号s(t),其由20Hz、220Hz和750Hz的正弦信号叠加而成，分析确定采样频率及数据分析长度，计算并绘出信号的频谱，指出各个频率份量。 （二）、数字滤波器设计 由图可以看出FFT变换具有对称性，幅频关于N/2偶对称。三角波序列和反三角波序列8点的FFT相同，而16点的FFT不同。这是因为三角波序列[1,2,3,4,4,3,2,1]以及反三角波序列[4,3,2,1,1,2,3,4]在不补零的情况下进行周期延拓得到的序列是一样的，所以得到的DFS是一样的，所以DFT也是一样的；但是在补零之后再周期延拓得到的序列就不一样，所以DFT也不一样了。 从图中还可以看到：同一序列16点的FFT是在8点的FFT基础上在其相邻两点的中间内插了一个点。我们知道x(n)的X(K)就是它的Z变换X(Z)在Z平面单位圆上N个等分点上的采样值。同一序列的Z变换是一样的。 高斯序列 上图是高斯序列中改变参数p,q所对应的时域以及幅频特性图。 由图可以看出在p=8的情况下，改变q,随着q的增加，时域波形图中对称轴p=8两边的点所对应的幅度均同等比例增加，零点个数逐渐减少。在q=8的情况下，改变p,随着p的增加，峰值右移，而且左边零点数逐渐增加，右边零点数减少。这种现象可以从高斯序列的数学表达式中得到解释。 再看幅频图，p=8的情况下，在附近幅度为0，呈现出带阻特性，改变q,随着q的增加，阻带范围逐渐变宽。在q=8的情况下，改变p,随着p的增加，幅值逐渐减小。如图在q=8,p=25的情况下, 在[0,2]中幅度几乎为0，频谱发生明显的泄露，同时混叠也随着出现。 1.3正弦序列FFT 下图是正弦序列时域以及FFT图： 序列周期N=2*/(/4)=8。在对序列进行M点FFT时，若取M为N的整数倍，则不会发生频谱泄露。由于FFT具有对称性，而且正弦序列16点的时域图是8点的时域图的一个周期延拓，所以在进行FFT时得到的序列形状是一样的，只是增加了点数。而20不是8的整数倍，所以在做20点FFT时，会发生频谱泄露现象。 1.4正弦衰减序列频率改变 A=50, f0=62.5,a=100,T=1/fs 下图是正弦衰减序列在频率为fs1=1000Hz和fs2=300Hz时的幅频图。分析如下： 此正弦衰减序列是离散的非周期的，所以其FFT后频谱是连续的周期的。由图可以看出fs=1000Hz时幅频图要陡峻些，谱峰位置在0.12，以及1.88附近，而fs=300Hz时，幅频图要平缓些，谱峰位置在0.38以及1.62附近。而且fs=1000Hz的谱峰要高fs为300Hz时的谱峰。 fs的改变会引起x(n)表达式中幅值的改变，这是产生谱峰高度不一的原因。同样fs的改变会引起sin项的改变，使谱峰位置发生变化，并且由表达式容易得出f