数字信号处理实验指导书(学生版).docVIP

“数字信号处理”实验指导书（一） 一、实验课程编码：105003 二、实验课程名称：数字信号处理 三、实验项目名称： 应用MATLAB分析离散信号频谱 四、实验目的 掌握应用MATLAB分析离散信号频谱的方法，即熟悉应用MATLAB分析离散信号的函数。 五、主要设备 安装有MATLAB软件的电脑 六、实验内容 编写MATLAB程序，实现下面题目： 1. 用快速卷积法计算下面两个序列的线性卷积。 ， 2．已知序列 （1）计算该序列DTFT的表达式，并画出N=10时的曲线； （2）编写MATLAB程序，利用FFT函数，计算N＝10时，序列x[k]的DTFT在的抽样值。利用hold函数，将抽样点画在的曲线上。 3.理解高密度频谱和高分辨率频谱的概念。 设 取0≤n≤9，求 将（1）中的补零加长到0≤n≤99，求 增加取样值的个数，取0≤n≤99，求 4. 用DFT对连续信号做谱分析。 设，用DFT分析的频谱结构，选择不同的截取长度Tp，观察截断效应，试用加窗的方法减少谱间干扰。 选取的参数： 频率 采样信号序列，是窗函数。选取两种窗函数：矩形窗函数和Hamming窗，后者在程序中调用函数Hamming产生宽度为N的Hamming窗函数向量。 对做2048点DFT，作为的近似连续频谱。其中N为采样点数，，为截取时间长度，取三种长度0.04s、2×0.04s、4×0.04s、8×0.04s。 5. 已知一连续信号为，试利用DFT近似分析其频谱。要求频率分辨率为1Hz，确定抽样频率、抽样点数N以及持续时间Tp。 说明：连续信号x(t)的频谱X(jΩ)可以由其离散信号x(n)的DFT近似求得： X(jΩ) ≈ T．FFT[x(n)]，其中T为采样周期。 计算连续信号x(t)的频谱，要考虑几个问题：1）频率混叠：如果x(t)不是有限带宽的，则采样频率fs选择的依据是，使由时域采样所造成的频率混叠小到可以忽略；2）频率分辨率：频率分辨率和信号在时域的持续时间成反比。3）截断效应：如果x(t)是无限长的，就必须进行截断，截断会引起吉布斯效应（波动），也会把窗函数的频谱引入信号频谱，造成混叠。 分析连续信号频谱的一般步骤为：1）先初步选择时间记录长度T0,使得其包括了大部分非零的数据，然后用逐次增大采样频率fs的方法来选择采样频率。最终使得时域采样造成的频率混叠可以近似的忽略不计。2）在确定了采样频率后，进一步选择时间记录长度T0，即逐次将T0加倍，因为此时采样频率相同，两个频谱之间的差别是由截断效应不同引起的，如果两个频谱的误差不大，则这个记录长度T0可以接受。 （本题直接给出了频率分辨率的要求，也就是记录长度T0可以直接确定（不必进行步骤2），只需进行步骤1逐次选择采样频率。） 6．验证频域采样定理。 产生一个三角波序列x(n)，长度为M=40； 计算N=64时的X(k)=DFT[x(n)]，并图示x(n)和X(k) 对X(k)在上进行32点抽样，得到X1(k)＝X(2k)，k＝0,1,…,31。 求X1(k)的32点IDFT，即x1(n)＝IDFT[X1(k)]。 绘制x1((n))32的波形图，观察x1((n))32和x(n)的关系，并加以说明。 七、实验步骤 1、熟悉与离散信号频谱分析相关的MATLAB函数（参考附录1）； 2、通过运行附录2中提供的例题，熟悉用MATLAB分析离散信号频谱的基本方法； 3、根据“六、实验内容”中各个题目的要求，编写MATLAB程序代码，调试程序，分析并保存结果。 八、实验结果 对实验练习题编写MATLAB程序并运行，在计算机上输出仿真结果。 附录1 主要的相关MATLAB函数 1.fft.m和ifft.m 调用格式：〔X〕＝fft（x） 〔x〕＝ifft（X） 〔X〕＝fft（x,N） 〔x〕＝ifft（X,N） 2.czt.m 调用格式：〔y〕＝czt（x,m,w,s） 3.fftshift.m 调用格式：〔y〕＝fftshift（x） 附录2 例题 例1 利用DFT的性质，编写MATLAB程序，计算下列序列的6点圆周卷积。 （1）x [n]= {1,-3,4,2,0,-2}, h[n]= {3,0,1,-1,2,1} （2）x[n]=cos(πn/2), h[n]=3n, n=0,1,2,3,4,5 [MATLAB程序]： N=6; xn=[1,-3,4,2,0,-2]; hn=[3,0,1,-1,2,1]; Xk=fft(xn,N); %计算N点的DFT[x(n)] Hk=fft(hn,N); %计算N点的DFT[h(n)] Yk=Xk.*Hk; %DFT[x(n