数字信号处理实验指导书（学生版）-广电工.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)]上进行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)].*DFT[h(n)] y=ifft(Yk,N) %计算N点的IDFT[Y(k)],即为x(n)和h(n)的圆周卷积 [运行结果]： y = 6.0000 -3.0000 17.0000 -2.0000 7.0000 -13.0000 [MATLAB程序]： N=6; n=0:N-1; xn=cos(pi*n/2); hn=3*n; Xk=fft(xn,N); %计算N点的DFT[x(n)] Hk=fft(hn,N); %计算N点的DFT[h(n)] Yk=Xk.*Hk; %DFT[x(n)].*DFT[h(n)] y=ifft(Yk,N) %计算N点的IDFT[Y(k)],即为x(n)和h(n)的圆周卷积 [运行结果]： y = -6.0000 -3.0000 18.0000 21.0000 6.0000 9.0000 例2、基本序列的离散傅立叶变换计算 复正弦序列：，余弦序列： 分别对以上序列求当N=16和N=8时的DFT，并绘出幅频特性曲线，对其结果进行分析。 [MATLAB程序]： %基本序列的离散傅立叶变换计算 N=16;N1