通信10-3李国庆 324.doc

《数字信号处理实验》指导书 赵红怡 姓名：李国庆 学号：10101060324 2013年5月实验一 离散傅里叶变换（DFT）对确定信号进行谱分析 实验目的 加深对DFT算法原理和基本性质的理解。 熟悉DFT算法和原理的编程方法。 学习用DFT对信号进行谱分析的方法，了解可能出现的误差及其原因，以便在实际中正确利用。 实验原理 一个连续信号的频谱可以用其傅里叶变换表示，即 若对进行理想采样可得采样序列 对进行DTFT，可得其频谱为： 其中数字频率与模拟频率的关系为： 的DFT为 若是限带信号，且在满足采样定理的条件下，是的周期延拓，是在单位圆上的等间隔采样值，即。 为在计算机上分析计算方便，常用来近似，这样对于长度为N的有限长序列（无限长序列也可用有限长序列来逼近），便可通过DFT求其离散频谱。 实验步骤 复习DFT变换及其性质。 编写DFT变换程序，分析DFT结果。 3．按实验内容要求，上机实验，并写出实验报告。 实验内容 1．用DFT对下列序列进行谱分析。 （1） （2） %a08.m N=100;n=0:1:N-1; xn=3*sin(0.04*pi*n); subplot(3,2,1);stem(n,xn);title(连续信号序列); Xk=dft(xn,N); %离散傅里叶变换 magXk=abs(Xk) %求幅度谱 phaXk=angle(Xk); %求相位谱 k=0:N-1; subplot(3,2,3);stem(k,magXk);title(幅度谱); subplot(3,2,5);stem(k,phaXk);title(相位谱); xn1=[1,1,1,1,0,0,0,0,0,0,0,0]; N1=length(xn1);n1=0:1:N1-1; subplot(3,2,2);stem(n1,xn1); title(离散信号序列); %信号序列 Xk1=dft(xn1,N1); %离散傅里叶变换 magXk1=abs(Xk1) phaXk1=angle(Xk1); k1=0:N1-1; subplot(3,2,4);stem(k1,magXk1);title(幅度谱); subplot(3,2,6);stem(k1,phaXk1);title(相位谱); 2．为了说明高密度频谱和高分辨率频谱之间的区别，考察序列 （1）当0≤n≤10时，确定并画出x(n)的离散傅里叶变换。 （2）当0≤n≤100时，确定并画出x(n)的离散傅里叶变换。 %a09.m n=[0:1:99];x=cos(0.48*pi*n)+cos(0.52*pi*n); n1=[0:1:9];y1=x(1:1:10); subplot(3,3,1);stem(n1,y1);title(信号x(n),0=n=9); axis([0,10,-2.5,2.5]);text(10.2,-2.5,n) Y1=fft(y1,10);magY1=abs(Y1(1:1:6)); k1=0:1:5;w1=2*pi/10*k1; subplot(3,3,4);stem(w1/pi,magY1);title(FFT幅度); axis([0,1,0,10]);text(1.05,0,w) xlabel(以pi为单位的频率) subplot(3,3,7);stem(w1/pi,angle(Y1(1:1:6))),title(FFT相位谱); n2=[0:1:99];y2=[x(1:1:10),zeros(1,90)]; subplot(3,3,2);stem(n2,y2); title(信号x(n),0=n=9+90zeros); axis([0,100,-2.5,2.5]);text(100.5,-2.5,n) Y2=fft(y2,100);magY2=abs(Y2(1:1:51)); k2=0:1:50;w2=2*pi/100*k2; subplot(3,3,5);stem(w2/pi,magY2);title(FFT幅度); axis([0,1,0,10]);text(1.05,0,w) xlabel(以pi为单位的频率) subplot(3,3,8);stem(w2/pi,angle(Y2(1:1:51))),title(FFT相位谱); axis([0,1,-5,5]); subplot