上海交通大学电子信息与电气工程学院DSP离散时间信号实验报告三(总10页).docVIP

应用傅里叶变换对信号进行频谱分析实验报告 刘运鹏 5080309226 实验目的： １、加深对DFT的理解 ２、学会应用DFT对信号进行频谱分析的方法 实验内容： １、复习DFT和FFT有关内容，掌握调用方法 ２、对所得结果加以讨论 题目一：31 为了说明高密度频谱与高分辨率频谱之间的区别，考虑序列： 求出它基于有限个样本的频谱。 当时，确定并画出x（n）的离散时间傅利叶变换。 当时，确定并画出x（n）的离散时间傅利叶变换。 参数修改：将0.48改为0.47，0.52改为0.53。 程序代码： N=11; n=0:N-1; x=cos(0.47*pi*n)+cos(0.53*pi*n); k=0:N-1; nk=n*k; Wn=exp(j*2*pi*-1/N); Wnnk=Wn*nk; X=x*Wnnk; Hx=abs(X); stem(k,Hx); N=101; n=0:N-1; x=cos(0.47*pi*n)+cos(0.53*pi*n); k=0:N-1; nk=n*k; Wn=exp(j*2*pi*-1/N); Wnnk=Wn*nk; X=x*Wnnk; Hx=abs(X); stem(k,Hx); 实验结果和分析： a的DFT图 b的DFT图 题目二：34 设的基本周期N=50，它的一个周期由下式给出： 的基本周期N=100，它的一个周期由下式给出： 这两个周期序列的区别在于它们的周期，但它们具有相同的非零样本。 a.求出的DFS，用stem函数画出它的幅度和角度样本。 b.求出的DFS，画出它的幅度和角度样本。 c.上面两个DFS图有什么不同。 提示：本题说明时域补零，则频域取样更密 参数修改： 将指数系数-0.3改为-0.2。 程序代码： x1=zeros(1,50); for i=1:26 x1(1,i)=(i-1)*exp(-0.2*(i-1)); end; X1=fft(x1); X1=[X1 X1]; k1=0:99; h1=abs(X1); pha1=angle(X1); stem(k1,h1); stem(k1,pha1); x2=zeros(1,100); for i=1:26 x2(1,i)=(i-1)*exp(-0.2*(i-1)); end; X2=fft(x2); X2=[X2 X2]; k2=0:199; h2=abs(X2); pha2=angle(X2); stem(k2,h2); stem(k2,pha2); 实验结果和分析： 后者比前者的采样更密集 X1幅度图 X1相角图 X2幅度图 X2相角图 X幅度图 X相角图 第三题：36 12点序列定义为： a．求出的DFT，画出它的幅度和相位曲线。 b．MATLAB画出的DTFT的幅度和相位曲线。 c．验证a中的DFT是的采样，采用hold函数把两图放在一幅图里。 d．有无可能从DFT重构DTFT？如果可能，给出重构所需的内插公 式；不可能则说明理由。 提示：参照题13（c），a调用fft（），b调用fplot（） 参数修改：没有修改参数。 程序代码： x=[1 2 3 4 5 6 6 5 4 3 2 1]; X=fft(x); k=0:11; h=abs(X); pha=angle(X); stem(k*2*pi/12,h); stem(k*2*pi/12,pha); x=[1 2 3 4 5 6 6 5 4 3 2 1]; n=0:11; w=linspace(0,12,1000000); X=x*exp(-j*n*w); h=abs(X); pha=angle(X); plot(w,h); plot(w,pha); 实验结果和分析： 前者是后者的频域采样 不能重构，频域采样频率不符合乃奎斯特采样定律 X DFT幅度图 X DFT相角图 X DTFT幅度图 X DTFT相角图 幅度比较 相角比较