资源描述：
《DFT 变换的性质及应用》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。
1、实验卡和实验报告 信息科学与工程学院课程名称数字信号处理实验课时4学时实验项目名称和编号DFT变换的性质及应用同组者姓名实验目的1、实现信号的DFT变换2、了解DFT应用：(1)用DFT计算卷积(2)用DFT对序列进行谱分析设计型实验，综合型实验实验环境MATLAB实验内容和原理任务1、用三种不同的DFT程序计算x(n)=(0.9)n(n=0,1,2,…,7)的傅立叶变换X(k)，并比较三种程序的计算机运行时间。步骤：a.用for循环语句编制函数文件，实现循环计算X(k)；
2、b.编写矩阵运算的函数文件，实现矩阵计算X(k)；c.调用FFT函数直接计算X(K)任务2、给定x(n)=nR16(n)，h(n)=R8(n)利用DFT实现两序列的线性卷积运算，并研究DFT的点数与混叠的关系，并用stem(n,y)画出相应的图形。任务3、讨论序列补零及增加数据长度对信号频谱的影响(1)求出序列x(n)=cos(0.48n)+cos(0.52n)基于有限个样点n=10的频谱；(2)求n=100时,取x(n)的前10个,后90个设为零，得到x(n)的频谱;(3)增加x(n)有效的样点数，取100个样点得到x(n)的频谱实验卡和实验报告
3、 信息科学与工程学院实验步骤方法关键代码一、任务1a.用for循环语句编制函数文件，实现循环计算X(k)；function[am,pha]=dft1(x)N=length(x);w=exp(-1i*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);endb.编写矩阵运算的函数文件，实现矩阵计算X(k)；function[am,pha]=
4、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);c.调用FFT函数直接计算X(K)function[amfft,phafft]=dft3(x)N=length(x);Xk=fft(x);amfft=abs(Xk);phafft=angle(Xk);三种程序的计算机运行时间x=ones(1,8);figure(1)[am,pha]=det1(x);t1=cputime实验卡和实验报
5、告 信息科学与工程学院subplot(3,1,1);stem(x);subplot(3,1,2);stem(am);subplot(3,1,3);stem(pha);figure(2)[am,pha]=dft2(x)t2=cputimesubplot(3,1,1);stem(x);subplot(3,1,2);stem(am);subplot(3,1,3);stem(pha);figure(3)[amfft,phafft]=dft3(x)t3=cputimesubplot(
6、3,1,1);stem(x);subplot(3,1,2);stem(am);subplot(3,1,3);stem(pha);二、任务2%N1+N2-1=23<32N=32;x=[0:15];xx=[x,zeros(1,16)];h=[ones(1,8),zeros(1,24)];Xk=fft(xx,N);Hk=fft(h,N);Yk=Xk.*Hk;y=ifft(Yk,N);n=0:N-1;实验卡和实验报告 信息科学与工程学院stem(n,y);holdon%N=N1=1
7、6N1=16;x1=[0:15];h1=[ones(1,8),zeros(1,8)];Xk1=fft(x1,N1);Hk1=fft(h1,N1);Yk1=Xk1.*Hk1;y1=ifft(Yk1,N1);n1=0:N1-1;stem(n1,y1,'.','m');(接上页)实验步骤方法关键代码三、任务3%x(n)基于10个样点的频谱figure(1)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(2,1,1);stem(n1,y1);title('
8、signalx(n),0<=n<=9');xlabel('n')axis([0,10,-2.5,2.5])Y1=fft(y
