数字信号处理课程设计指导书4.docVIP

数字信号处理上机指导 设计一 正余弦信号的谱分析 设计目的 用DFT实现对正余弦信号的谱分析； 观察DFT长度和窗函数长度对频谱的影响； 对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(‘How many apples’) 会显示How many apples 从键盘输入 3 会显示 R=3 2. 一维快速傅里叶变换函数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点FFT subplot(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的运行结果 设计内容 对一个频率为10Hz的连续余弦信号，以采样频率为64Hz得到离散余弦序列，对此余弦序列进行32点进行谱分析，画出其频谱图；若将信号频率改为11Hz，其他参数不变，重新画出该序列的频谱图，观察频谱泄漏现象，分析原因。 考察DFT的长度对双频率信号频谱分析的影响。设待分析的信号为 （1.2） 令两个长度为16的正余弦序列的数字频率为及。取N为四个不同值16，32，64，128。画出四个DFT幅频图，分析DFT长度对频谱分辨率的影响。 3. 在上题中若把两个正弦波的频率取得较近，令，，试问怎样选择FFT参数才能在频谱分析中分辨出这两个分量？ 【四】课程设计报告要求 完成课程设计任务后，应按要求提交课程设计报告。设计报告应包含如下几个方面的内容。 课程设计目的课程设计要求。详细设计过程。调试分析。结果分析与体会。 附录或参考资料。 （2.1） 其中，称为载波，为载波频率，称为单频调制信号，为调制正弦波信号频率，且满足。由（2.1）式可见，所谓抑制载波单频调制信号，就是两个正弦信号相乘，它有2个频率成分：和频+，差频-，这两个频率成分关于载波频率对称。所以，1路抑制载波单频调幅信号的频谱图是关于载波频率对称的两根谱线。 复合信号st产生函数mstg清单： function st=mstg %产生信号序列st，并显示st的时域波形和频谱