实验3-离散信号的离散频域分析.doc

数字信号处理A 实验项目名称: 3离散信号的离散频域分析 学 院:______计算机与通信工程____ 专 业:______ _通信工程 _________ 学 号:______201xxxxxxxxxx_______ 班 级:______ 通信1x-x_________ 报 告 人:_____ ___xx__________ 指 导 老 师:___ 胡双红 _ _______ 实 验 时 间:_______201x-1x-2x________ 实验目的： 1、掌握离散时间系统的DFT的MATLAB实现； 2、熟悉DTFT和DFT之间的关系。 3、了解信号不同变形的DFT与原信号DFT之间的关系 实验内容： ： 实验： 2、在信号每两个相邻样本之间插入一个零值，扩充为16点序列，作DFT，画出幅度谱和相位谱，并与原序列的DFT进行比较。 3、将信号以8为周期扩展，得到长为16的两个周期，作DFT，画出幅度谱和相位谱，并与原序列的DFT进行比较。 实验： n=0:7; x=(0.9*exp(j*pi/3)).^n; x1=[x zeros(1,8)];x2=[x1 zeros(1,16)]; Xk=dft(x,8);Xk1=dft(x1,16);Xk2=dft(x2,32); w=0:pi/200:2*pi;X=dtft(x,n,w); magX=abs(X);angX=angle(X); magXk=abs(Xk);angXk=angle(Xk);magXk1=abs(Xk1);angXk1=angle(Xk1); magXk2=abs(Xk2);angXk2=angle(Xk2); subplot(4,2,1);plot(w/pi,magX);xlabel(w/pi);ylabel(幅度|X|);grid on; subplot(4,2,2);plot(w/pi,angX);xlabel(w/pi);ylabel(相位(rad/π)); subplot(4,2,3);stem(0:7,magXk);xlabel(k);ylabel(幅度|X(k)|);axis([0,8,0,6]); subplot(4,2,4);stem(0:7,angXk);xlabel(k);ylabel(相位(rad/π));axis([0,8,-2,2]); subplot(4,2,5);stem(0:15,magXk1);xlabel(k);ylabel(幅度|X1(k)|);axis([0,16,0,6]); subplot(4,2,6);stem(0:15,angXk1);xlabel(k);ylabel(相位(rad/π));axis([0,16,-2,2]); subplot(4,2,7);stem(0:31,magXk2);xlabel(k);ylabel(幅度|X2(k)|);axis([0,32,0,6]); subplot(4,2,8);stem(0:31,angXk2);xlabel(k);ylabel(相位(rad/π));axis([0,32,-2,2]); n=0:7; x=(0.9*exp(j*pi/3)).^n; for i=1:8 %将原序列每两个相邻样本之间插入一个零值，扩展为长为16的序列 y1(2*i-1:2*i)=[x(i) 0]; end Yk1=dft(y1,16); Warning: Function call dft invokes inexact match D:\MATLAB7\work\DFT.M. magYk1=abs(Yk1);angYk1=angle(Yk1); y2=[x x]; %将信号以8为周期扩展，得到长为16的两个周期 Yk2=dft(y2,16); magYk2=abs(Yk2);angYk2=angle(Yk2); Xk=dft(x,8); magXk=abs(Xk);angXk=angle(Xk); subplot(3,2,1);stem(0:7,magXk);xlabel(k);ylabel(幅度|X(k)|); subplot(3,2,2);stem(0:7,angXk);xlabel(k);ylabel(相位(rad/π)); subplot(3,2,3);stem(0:15,magYk1);xlabel(k);ylabel(幅度|Y1(k)|); subplot(3,2,4);stem(0:15,angYk1);xlabel(k);ylabel(相位(rad/π)); subplot(3,2,5);stem(