快速傅立叶变换(FFT)算法及其应用.docVIP

实验八 快速傅立叶变换（FFT）算法及其应用 1. 忽略 2. 利用MATLAB编程完成计算实验内容1、2、及3的要求，绘出相应图形。并与理论计算相比较，说明实验结果的原因。 1、观察三角波序列和反三角波序列的时域和幅频特性 ①用8点FFT分析信号和的幅频特性，观察两者的序列形状和频谱曲线有什么异同？（注意：这时候的可以看作是经过圆周移位后得到的）绘制两者的序列和幅频特性曲线。 ②在的和末尾补零，用32点 FFT分析这两个信号的幅频特性，观察幅频特性发生了什么变化?两个信号之间的FFT频谱还有没有相同之处？这些变化说明了什么？ 程序如下： N=8; n=0:N-1; xc=0:7; xd=0:7; for m=0:7; if(0=m m=3) xc(m+1)=m; xd(m+1)=4-m; elseif(4=m m=7) xc(m+1)=8-m;xd(m+1)=m-4; end end figure(1); subplot(2,1,1);stem(n,xc); title(三角波序列); subplot(2,1,2); stem(n,xd); title(反三角波序列); xck=fft(xc,N); xdk=fft(xd,N); figure(2); subplot(2,1,1); stem(n,abs(xck)); title(三角波序列8点fft); subplot(2,1,2); stem(n,abs(xdk)); title(反三角波序列8点fft); figure(3); xck1=fft(xc,32); xdk1=fft(xd,32); subplot(2,1,1); stem(0:4*N-1,abs(xck1)); title(三角波序列32点fft); subplot(2,1,2); stem(0:4*N-1,abs(xdk1)); title(反三角波序列32点fft); 2.2N点实数序列 N=64。用一个64点的复数FFT程序，一次算出，并绘出。 解 ： 原理：①把x[n] (0n2N-1)中的奇数位取出来组成o[n] (0nN-1)，把x[n] (0n2N-1)中的偶数位取出来组成e[n] (0nN-1) ②设v[n]=e[n]+j*o[n]，对其做64点DSP得到V[k] ③则e[n]和e[n]的DSP分别为E[k]=(V[k]+V* [N-k])/2，O[k]=(V[k]-V* [N-k])/2 ④则x[n]的DSP为 ??? ?? ??? 程序如下： clear;clf; N=64; n=0:2*N-1; x=cos(2*pi*7*n/N)+cos(2*pi*19*n/N)/2; subplot(3,2,1); X=fft(x); stem(real(x)); xlabel(原信号); subplot(3,2,3); stem(real(X)); xlabel(x的fft的实部（对照）); subplot(3,2,4); stem(imag(X)); xlabel(x的fft的虚部（对照）); for i=1:N e(i)=x(2*i); o(i)=x(2*i-1); end; v=o+j*e; V=fft(v); Vt=real(V)-j*imag(V); for i=1:N O(i)=(V(i)+Vt(mod((N+1-i),N)+1))/2; E(i)=(V(i)-Vt(mod((N+1-i),N)+1))/(2*j); end; for i=1:2*N X0(i)=O(mod(i-1,N)+1)+exp(-j*2*pi*(i-1)/(2*N))*E(mod(i-1,N)+1); end; subplot(3,2,5); stem(real(X0)); xlabel(通过64点FFT算法出的x的实部); subplot(3,2,6); stem(imag(X0)); xlabel(通过64点FFT算法出的x的虚部); 运行结果如图所示： 3. 已知某序列在单位圆上的N=64等分样点的Z变换为 用N点IFFT程序计算，绘出和。 解 ： 程序如下：clear;clf; N=64; k=0:N-1; X=1./(1-0.8*exp(-j*2*pi*k/N)); subplot(3,1,1); stem(abs(X)); xlabel(|X[k]|); x=ifft(X); subplot(3,1,2); stem(abs(x)); xlabel(abs(x[n])=a