数字信号处理课程设计指导书90172VIP免费

数字信号处理上机指导设计一正余弦信号的谱分析【一】设计目的1.用DFT实现对正余弦信号的谱分析；2.观察DFT长度和窗函数长度对频谱的影响；3.对DFT进行谱分析中的误差现象获得感性认识。【二】设计原理一、谱分析原理数字信号处理方法的一个重要用途是在离散时间域中确定一个连续时间信号的频谱，通常称为频谱分析，更具体地说，它也包括确定能量谱和功率谱。数字频谱分析可以应用在很广阔领域，频谱分析方法是基于以下的观测：如果连续时间信号是带限的，那么它的离散时间等效信号的DFT进行谱分析。然而，在大多数情况下，是在范围内定义的，因此也就定义在的无线范围内，要估计一个无限长信号的频谱是不可能的。实用的方法是：先让模拟连续信号通过一个抗混叠的模拟滤波器，然后把它采样成一个离散序列。假定反混叠滤波器的设计是正确的，则混叠效应可以忽略，又假设A/D变换器的字长足够长，则A/D变换中的量化噪声也可忽略。假定表征正余弦信号的基本参数，如振幅、频率和相位不随时间改变，则此信号的傅立叶变换可以用计算它的DTFT得到（1.1）实际上无限长序列首先乘以一个长度为M的窗函数，使它变成一个长为M的有限长序列，，对求出的DTFT应该可以作为原连续模拟信号的频谱估计，然后求出在区间等分为N点的离散傅立叶变换DFT。为保证足够的分辨率，DFT的长度N选的比窗长度M大，其方法是在截断了的序列后面补上N－M个零。计算采用FFT算法。二、MATLAB函数介绍1.输入函数input()格式：R=input(string)功能：在屏幕上显示input括号后的’string’内容，提示用户从键盘输入某值，并将输入的值赋给R。例如，在命令窗口输入R=input(‘Howmanyapples’)会显示Howmanyapples从键盘输入3会显示R=32.一维快速傅里叶变换函数fft()格式：y=fft(x)y=fft(x,n)说明：fft函数用于计算矢量或矩阵的傅里叶变换。３．幅度函数abs()格式：abs（X）功能：对X取绝对值，当X是复数时，得到X的复模值。例1-1对余弦序列进行8点谱分析，绘出序列和幅频特性曲线。MATLAB命令如下：N=input('输入谱分析的长度')n=0:N-1;x=cos(pi*n/4);subplot(1,2,1),stem(n,x),xlabel('n');ylabel('x1(n)');%绘制序列的波形title('余弦序列');X=abs(fft(x,N));%求余弦序列的N点FFTsubplot(1,2,2)k=0:N-1;stem(k,X)%绘制序列的幅频特性曲线xlabel('k');ylabel('X(k)');string=[num2str(N),'点FFT幅频曲线'];title(string);程序运行结果如图1.1所示。图1.1例1-1的运行结果【三】设计内容1.对一个频率为10Hz的连续余弦信号，以采样频率为64Hz得到离散余弦序列，对此余弦序列进行32点进行谱分析，画出其频谱图；若将信号频率改为11Hz，其他参数不变，重新画出该序列的频谱图，观察频谱泄漏现象，分析原因。2.考察DFT的长度对双频率信号频谱分析的影响。设待分析的信号为（1.2）令两个长度为16的正余弦序列的数字频率为及。取N为四个不同值16，32，64，128。画出四个DFT幅频图，分析DFT长度对频谱分辨率的影响。3.在上题中若把两个正弦波的频率取得较近，令，，试问怎样选择FFT参数才能在频谱分析中分辨出这两个分量？【四】课程设计报告要求完成课程设计任务后，应按要求提交课程设计报告。设计报告应包含如下几个方面的内容。1.课程设计目的2.课程设计要求。3.详细设计过程。4.调试分析。5.结果分析与体会。6.附录或参考资料。设计二数字滤波器的设计及实现【一】设计目的1.熟悉IIR数字滤波器和FIR数字滤波器的设计原理和方法；2.学会调用MATLAB信号处理工具箱中的滤波器设计函数设计各种IIR和FIR数字滤波器，学会根据滤波要求确定滤波器指标参数；3.掌握用IIR和FIR数字滤波器的MATLAB实现方法，并能绘制滤波器的幅频特性、相频特性；4.通过观察滤波器的输入、输出信号的时域波形及其频谱，建立数字滤波的概念。【二】设计原理抑制载波单频调幅信号的数学表达式为（2.1）其中，称为载波，为载波频率，称为单频调制信号，为调制正弦波信号频率，且满足>。由（2.1）式可见，所谓抑制载波单频调制信号，就是两个正弦信号相乘，它有2个频率成分：和频+，差频-，这两个频率成分关于载波频率对称。所以，1路抑制载波单频调幅信号的频...