MATLAB离散信号的产生和频谱分析实验报告.docxVIP

实验一 离散信号的产生和频谱分析 一、实验目的 仿真掌握采样定理。 学会用FFT进行数字谱分析。 掌握FFT进行数字谱分析的计算机编程实现方法。 培养学生综合分析、解决问题的能力，加深对课堂内容的理解。 二、实验要求 掌握采样定理和数字谱分析方法；编制FFT程序；完成正弦信号、线性调频信号等模拟水声信号的数字谱分析； 三、实验内容 单频脉冲（CWP）为。式中，是矩形包络，是脉冲持续时间，是中心频率。 矩形包络线性调频脉冲信号（LFM）为。式中，是线性调频指数。瞬时频率是时间的线性函数，频率调制宽度为。 设参数为，，，采样频率。 1．编程产生单频脉冲、矩形包络线性调频脉冲。 2．编程实现这些信号的谱分析。 3．编程实现快速傅立叶变换的逆变换。 四、实验原理 1、采样定理 所谓抽样，就是对连续信号隔一段时间T抽取一个瞬时幅度值。在进行模拟/数字信号的转换过程中，当采样频率fs大于信号中最高频率f的2倍时(fs=2f)，采样之后的数字信号完整地保留了原始信号中的信息，一般实际应用中保证采样频率为信号最高频率的5～10倍；采样定理又称奈奎斯特定理。 2、离散傅里叶变换（FFT） 长度为N的序列的离散傅立叶变换为： N点的DFT可以分解为两个N/2点的DFT，每个N/2点的DFT又可以分解为两个N/4点的DFT。依此类推，当N为2的整数次幂时（），由于每分解一次降低一阶幂次，所以通过M次的分解，最后全部成为一系列2点DFT运算。以上就是按时间抽取的快速傅立叶变换(FFT)算法。当需要进行变换的序列的长度不是2的整数次方的时候，为了使用以2为基的FFT，可以用末尾补零的方法，使其长度延长至2的整数次方。 3、离散傅里叶反变换 序列的离散傅立叶反变换为 五、实验步骤 1．编程产生单频脉冲、矩形包络线性调频脉冲。 2．应用快速傅立叶变换（FFT）求这两种信号的频谱，分析离散谱位置、归一化频率、实际频率的关系。 调用函数Y=fft(x) or Y=fft(x,N) or Y=fft(x,N,dim)。 3．对于步骤2的结果，应用快速傅立叶变换的逆变换（IFFT）求两种信号的时域波形，并与已给的单频脉冲、矩形包络线性调频脉冲和伪随机脉冲信号波形进行对照。 调用函数x=ifft(Y) or x=ifft(Y,N) or x=ifft(Y,N,dim)。 4．对于步骤2的结果，进行频谱移位调整。将FFT变换的结果Y（频谱数据）进行移位调整，使其符合频谱常观表示方法，调整后，频谱的直流成分（即频率为0处的值）移到频谱的中间位置。分析离散谱位置、归一化频率、实际频率的关系。 移位调整调用函数Z=fftshift(Y)。频率间隔为Fs/N，频率范围为Fs/N*[-N/2：N/2-1]。 六、实验结果及分析 图1 CW信号及LFM信号时域波形图 如图1，为产生的单频脉冲、矩形包络线性调频脉冲信号的抽样信号。 图2 CW信号及LFM信号频谱图 如图2，为两信号的傅里叶变换波形，为使运算简单，取采样点数为2的整数次方，本体中总点数为5000，故采样点数可去4096个；在频谱移位调整之前，FFT所得的频谱图并不能够反映真实信号的频率。如CW信号真实频率为20KHZ，而图形中显示的为30KHZ。 图3 频谱移位前后CW信号及LFM信号频谱图对比 如图3，为两信号的FFT频谱移位调整后频谱图，由图可知，调整后频谱为真实信号的频谱。 图4 傅里叶反变换前后CW信号及LFM信号时域波形对比 如图4，IFFT后的图形与原信号一致，当改变FFT取样点数时（比如取样点过少），傅里叶反变换后的图形与原图形就会发生变化。 图5 归一化后CW信号及LFM信号频谱图 六、讨论、思考题 1．回忆课堂上所讲的DFT、FFT的概念和数字谱分析方法。 设x(n)是一个长度为M的有限长序列，则定义x(n)的N点离散傅里叶变换为X(k)=DFT[x(n)]=n=0N-1x(n)WNkn k=0，1，…，N-1。然而直接计算DFT的计算量与变换区间长度N的平方成正比，当N较大时，计算量太大，直接用DFT算法进行谱分析和信号的实时处理是不切实际的。而FFT算法是将N电DFT分解为几个较短的DFT，可使乘法次数大大减少。常用的FFT算法是DIT-FFT和DIF-FFT。 而所谓的数字信号谱分析，就是计算信号的离散傅里叶变换，就是对连续信号和系统，可以通过时域采样，应用DFT进行近似谱分析。 2．总结用MATLAB工具箱函数实现谱分析的方法。 （1）在MATLAB中使用工具箱函数y=fft(x)得出连续信号的离散谱。 （2）使用y=fftshift(y)进行频谱搬移。 3．分析离散谱位置、归一化频率、实际频率的关系。 在经过fftshift进行频谱搬移将所分析信号的频谱放在正常频谱位置上后