DSP研究性学习报频谱计算.docVIP

《近代数字信号处理》课程研究性学习报告 姓名 舒晓萌 学号 同组成员 指导教师 魏杰 时间 2010-4-10 Matlab习题补充： M2-1 利用DFT的性质，编写一MATLAB程序，计算下列序列的循环卷积。 g[k]={1,-3,4,2,0,-2}, h[k]={3,0,1,-1,2,1}; x[k]=cos(pi*k/2), y=,k=0,1,2,3,4,5. （1） N=6; g=[1,-3,4,2,0,-2]; k=[3,0,1,-1,2,1]; x1=[g zeros(1,N-length(g))]; x2=[k zeros(1,N-length(k))]; m=0:N-1; x1=x1(mod(-m,N)+1); H=zeros(N,N); for n=1:N; H(1,:)=x1; H(n+1,:)=circshift(H(n,:),[0,1]); end y=H*x2; y=y’ y = 6 -3 17 -2 7 -13 6 自主学习方面：利用到mod求余函数，还利用到circshift求循环位移函数。 b = circshift(a,[x,y]) 其中a为待移位的矩阵，x表示上下移位位数（正数表示向下移位），y表示左右移位位数（正数表示向右移位）N=6; k=0:N-1; x=cos(pi*k/2); y=3.^k; x1=[g zeros(1,N-length(x))]; x2=[k zeros(1,N-length(y))]; x1=x1(mod(-k,N)+1); H=zeros(N,N); for n=1:N; H(1,:)=x1; H(n+1,:)=circshift(H(n,:),[0,1]); end y=H*x2; y=y y = 5 25 3 -7 -5 9 5 M2-2 已知序列 计算该序列DTFT的表达式，并画出N=10时，的曲线。 编写一MATLAB程序，利用fft函数，计算N=10时，序列x[k]的DTFT在的抽样值。利用hold函数，将抽样点画在曲线上。 (1) N=10; k=0:N-1; L=512; x=cos(k*pi/2/N); X=fftshift(fft(x,L)); plot(2*pi*(0:L-1)/L,abs(X)) grid on; (2) N=10; k=0:N-1; L=512; x=cos(k*pi/2/N); X=fftshift(fft(x,L)); w=2*pi*(-1/2+(0:L-1)/L); plot(w,abs(X),r);grid on;hold on; m=N*(-1/2+(0:N-1)/N); stem(m,abs(X),b); hold off; 怎么规定横坐标？ DFT近似计算信号频谱专题研讨 【目的】 (1) 掌握利用DFT近似计算不同类型信号频谱的原理和方法。 (2) 理解误差产生的原因及减小误差的方法。 (3) 培养学生自主学习能力，以及发现问题、分析问题和解决问题的能力。 【研讨题目】 基本题 M2-0 已知一离散序列为 (1)用L=32点DFT计算该序列的频谱，求出频谱中谱峰的频率； (2)对序列进行补零，然后分别用L=64、128、256、512点DFT计算该序列的频谱，求出频谱中谱峰的频率； (3)讨论所获得的结果，给出你的结论。该结论对序列的频谱计算有何指导意义？ 【题目分析】 本题讨论补零对离散序列频谱计算的影响。 【温磬提示】 在计算离散非周期序列频谱时常用(/(作为横坐标，称(/(为归一化频率? 【仿真结果】 (1) peak value 14.5561 frequency of peak value 3 (2) peak value