资源描述：
《八、DFT的原理及实现.doc》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。
1、实验题目DFT原理及实现 实验时间 实验台号 实验目的DFT算法原理；序列的DFT运算； 实验仪器PC台式机，软件MATLAB实验预习实验预习：1、DTFT和IDTFT公式；2、DFT变换公式； 3、DFT性质；实验原理：1、DTFT：可改写成：2、逆变换IDTFT：3、DFT：实验步骤应用函数for();首先用定义法定义DFT运算函数：1、二次循环法：functionX=DFTfor(xn)N=length(xn);X=zeros(1,N);fork=0:N-1;forn=0:N-1;X(k+1)=X(k+1)+xn(n+1)*exp(-
2、j*2*pi*n*k/N);endend2、矩阵乘法计算DFT：functionXk=DFTmat(xn);N=length(xn);n=0:N-1;k=n;nk=n'*k;WN=exp(-j*2*pi/N);Wnk=WN.^nk;Xk=xn*Wnk;题目8-1：设，计算在11点的DFT？n=-5:5;xn=(-0.9).^n;k1=0:1000;w=(pi/500)*k1;X=xn*(exp(-j*pi/500)).^(n'*k1);%求连续频谱magX=abs(X);angX=angle(X);Xk=DFTmat(xn);%调用矩阵相乘
3、法N=length(xn);k=0:N-1;Xk1=Xk.*exp(j*2*pi*5*k/N);%利用xn实际取值范围求解magXk=abs(Xk1);angXk=angle(Xk1);subplot(2,1,1);plot(w/pi,magX,'--');holdon;stem(2*k/N,magXk);holdoff;axis([0,2,0,15]);gridon;xlabel('omega(xpi)');ylabel('幅度');subplot(2,1,2);plot(w/pi,angX/pi,'--');holdon;stem(
4、2*k/N,angXk/pi);holdoff;axis([0,2,-1,1]);gridon;xlabel('omega(xpi)');ylabel('相位（弧度/pi)');%结果8-1题目8-2：二、设用直接法计算其4点DFT？n=0:3;xn=[1,1,1,1];k1=0:1000;w=(pi/500)*k1;X=xn*(exp(-j*pi/500)).^(n'*k1);magX=abs(X);angX=angle(X);N=length(xn);n1=0:N-1;k=n1;nk=n1'*k;WN=exp(-j*2*pi/N)
5、;Wnk=WN.^nk;Xk=xn*Wnk;magXk=abs(Xk);angXk=angle(Xk);subplot(2,1,1);plot(w/pi,magX,'k--');holdon;stem(2*k/N,magXk);holdoff;axis([0,2,0,5]);gridon;xlabel('omega(xpi)');ylabel('幅度
6、X(k)
7、');subplot(2,1,2);plot(w/pi,angX/pi,'k--');holdon;stem(2*k/N,angXk/pi);holdoff;axis([0,2,
8、-1,1]);gridon;xlabel('omega(xpi)');ylabel('相位/pi');%图像8-2实验数据结果及分析8-18-2结果分析：1、从函数定义上可以看出矩阵法比二次循环法计算要快。思考题教师评语 成绩: 日期:
