实验三 离散傅里叶变换.doc

实验三 离散傅里叶变换 23320078104194 07通信工程 戴文 一·实验目的 1·掌握离散傅里叶级数 2·掌握DFT变换 3·掌握DFT特性 4·掌握利用DFT计算线性卷积 5·掌握快速傅里叶变换（FFT） 二·实验原理 1·离散傅里叶级数 2·离散傅里叶变换 3·DFT特性 4·利用DFT计算线性卷积 5·快速傅里叶变换 三·实验内容 1（2）. clear all; xn=[1,2,3,0] Xk=dfs(xn,4) 运行结果Xk = 6.0000 -2.0000 - 2.0000i 2.0000 + 0.0000i -2.0000 + 2.0000i 2_2. clear all; x=[1,1,1,1,1,1]; a=dft(x,6) figure(1); subplot(2,1,1); stem(abs(a)) subplot(2,1,2); stem(angle(a)) figure(2); b=idft(a,6) subplot(2,1,1); stem(abs(b)) subplot(2,1,2); stem(angle(b)) 2（3）_a. n=0:10; x=cos(0.47*pi.*n)+cos(0.53*pi*n); a=dft(x,11); figure(1); subplot(2,1,1); stem(abs(a)) subplot(2,1,2); stem(angle(a)) 2（3）_b. n=0:10; x=cos(0.47*pi.*n)+cos(0.53*pi*n); a=[x,zeros(1,90)]; figure(1); stem(a) y=dft(a,101); figure(2); subplot(2,1,1); stem(abs(y)) subplot(2,1,2); stem(angle(y)) figure(3); yy=fft(a,101); subplot(2,1,1); stem(abs(yy)) subplot(2,1,2); stem(angle(yy)) 2（3）_c. n=0:500; x=cos(0.47*pi.*n)+cos(0.53*pi*n); a=dft(x,501); figure(1); subplot(2,1,1); stem(abs(a)) subplot(2,1,2); stem(angle(a)) 3_1. n=0:10; x=10*(0.8).^n; figure(1); stem(x) figure(2); a=cirshftt(x,15,11); stem(a) 3(2). x1=[1,2,2]; x2=[1,2,5,4]; y1 = circonvt(x1,x2,5) x1=[1,2,2]; x2=[1,2,5,4]; y2 = circonvt(x1,x2,6) 4. x=[1,2,3,4,5,6,7,8,9,10]; h=[1,0,-1]; figure(1); subplot(2,1,1); stem(x) subplot(2,1,2); stem(h); b=ovrlpsav(x,h,4); figure(2); subplot(2,1,1); stem(abs(b)) subplot(2,1,2); stem(angle(b)) 5.1 n=0:44; x=2*sin(4*pi*0.1.*n)+5*cos(8*pi*0.1.*n); a=fft(x,45); stem(abs(a)) 5.2 w=randn(1,64); x=2*sin(4*pi*0.1*n)+5*cos(8*pi*0.1.*n)+0.8*w; a=fft(x,64); stem(abs(a)) 图形： 2_2 2_3a 2_3b 2_3c 3 4 5 四·思考题 1·clear all;clc; [x,fs,bits]=wavread(ding); Sound(x,fs,bits); n=0:999; t=0.1*n; y=fft(x); a=size(x) %数据量 subplot(211); stem(abs(x)) subplot(212); stem(abs(y)) 2·clear all;clc; N=64;n=0:N-1;t=0.001.*n; noise=randn(1,N); x=cos(100*pi*t)+cos(240*pi*t); x1=cos(100*pi*t)+cos(240*pi*t)+noise; y=fft(x); y1=fft(x1); subplot(2,1,1);stem(n,abs(y)