数字信号处理离散傅里叶变换性质及应用.docx

数字信号处理实验题目:离散傅里叶变换的性质及应用学院:专业:学生姓名:班级/学号指导老师:一、、实验仪器及材料计算机,MATLAB软件三、,并比较三种程序计算机运行时间。(1),用循环变量逐点计算;(2),完成下列矩阵运算:(3)调用fft库函数,直接计算;(4)分别调用上述三种不同方式编写的DFT程序计算序列的离散傅里叶变换,并画出相应的幅频和相频特性,再比较各个程序的计算机运行时间。M函数文件如下::function[Am,pha]=dft1(x)N=length(x);w=exp(-j*2*pi/N);fork=1:N  sum=0;  forn=1:N     sum=sum+x(n)*w^((k-1)*(n-1));  end  Am(k)=abs(sum);  pha(k)=angle(sum);end :function[Am,pha]=dft2(x)N=length(x);n=[0:N-1];k=[0:N-1];w=exp(-j*2*pi/N);nk=n'*k;wnk=w.^(nk);Xk=x*wnk;Am=abs(Xk);pha=angle(Xk);:function[Am,pha]=dft3(x)Xk=fft(x);Am=abs(Xk);pha=angle(Xk);源程序、运行结果及结论1、源程序<1>function[Am,pha]=dft1(x)N=length(x);w=exp(-j*2*pi/N);fork=1:Nsum=0;forn=1:Nsum=sum+x(n)*w^((k-1)*(n-1));endAm(k)=abs(sum);pha(k)=angle(sum);endend<2>function[Am,pha]=dft2(x)N=length(x);n=[0:N-1];k=[0:N-1];w=exp(-j*2*pi/N);nk=n'*k;wnk=w.^(nk);Xk=x*wnk;Am=abs(Xk);pha=angle(Xk);end<3>function[Am,pha]=dft3(x)Xk=fft(x);Am=abs(Xk);pha=angle(Xk);end<4>clearall;clc;x=[ones(1,8),zeros(1,256-8)];t=cputime;[Am1,pha1]=dft1(x);t1=cputime-t;t=cputime;[Am2,pha2]=dft2(x);t2=cputime-t;t=cputime;[Am3,pha3]=dft3(x);t3=cputime-t;subplot(6,1,1);stem(Am1);title('幅频特性1');subplot(6,1,2);stem(pha1);title('相频特性1');subplot(6,1,3);stem(Am2);title('幅频特性2');subplot(6,1,4);stem(pha2);title('相频特性2');subplot(6,1,5);stem