DFT变换的性质与应用.docxVIP

.专业整理 .专业整理. . .学习帮手. 课 程 名 称 数字信号处理 实验课时 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)，并比较三种程序的计算机运行时间。 步骤： 用for循环语句编制函数文件，实现循环计算 X(k); 编写矩阵运算的函数文件，实现矩阵计算 X(k); 调用FFT函数直接计算X(K) 任务2、给定x(n) = nR16 (n) , h(n) = R8 (n) 利用DFT实现两序列的 线性卷积运算，并研究 DFT的点数与混叠的关系，并用 stem(n,y)画出相应的图形。 任务3、讨论序列补零及增加数据长度对信号频谱的影响 (1)求出序列x(n)=cos(0.48 n)+cos(0.52 n) 基于有限个样点 n=10的频谱； ⑵求n=100时，取x(n)的前10个，后90个设为零，得到 x(n)的频谱； (3)增加x(n)有效的样点数，取100个样点得到x(n)的频谱  实 验 步 骤 方 法 关 键 代 码 一、任务1 用for循环语句编制函数文件，实现循环计算 X(k); function [am,pha]=dft1(x) N=le ngth(x); w=exp(-1i*2*pi/N); for k=1:N sum=0; for n=1:N sum=sum+x( n)*w^( (k_1)*( n-1)); end am(k)=abs(sum); pha(k)=a ngle(sum); end 编写矩阵运算的函数文件，实现矩阵计算 X(k)； function [am,pha]=dft2(x) N=le ngth(x); n=[0:N-1]; k=[0:N-1]; w=exp(-j*2*pi/N); n k= n *k; wn k=w.A( nk); Xk=x*w nk; am=abs(Xk); pha=a ngle(Xk); 调用FFT函数直接计算X(K) fun ctio n [amfft,phafft]=dft3(x) N=le ngth(x); Xk=fft(x); amfft=abs(Xk); phafft=a ngle(Xk); 三种程序的计算机运行时间 x=on es(1,8); figure(1) [am,pha]=det1(x); t1=cputime  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=cputime subplot(3,1,1); stem(x); subplot(3,1,2); stem(am); subplot(3,1,3); stem(pha); figure(3) [amfft,phafft]=dft3(x) t3=cputime subplot(3,1,1); stem(x); subplot(3,1,2); stem(am); subplot(3,1,3); stem(pha); 二、任务2 %N1+N2-仁2332 N=32; x=[0:15]; xx=[x,zeros(1,16)]; h=[o nes(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); hold on %N=N1=16 N1=16; x1= [0:15]; h1=[o nes(1,8),zeros(1,8)]; Xk1=fft(x1,N1); Hk1=fft(h1,N1); Yk仁Xk1.*Hk1; y1= ifft(Yk1,N1); n1=0:N1-1; stem( n1,y1,.,m); 三、任务3 %x(n)基于10个样点的频谱 figure⑴ 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(signal x(n),0=n=9); 页 xlabel( n) ) axis([0,10,-2.5,2.5]) 丫仁 fft(y1);magY1=abs(Y1(1