实验三FFT在卷积积分和谱分析中的应用.docxVIP

实验三FFT在卷积积分和谱分析中的应用 实验目的 在理论学习的基础上，通过本实验，加深对FFT的理 解，熟悉FFT 了程序。 熟悉掌握FFT实现两个序列的线性卷积的方法 熟悉应用FFT进行信号谱分析过程中可能出现的问题 以便在实际屮正确应用FFT. 学习用FFT对连续信号进行谱分析的方法和参数选择 原则。 常握用FFT对连续信号进行谱分析时，由模拟信号采 样得到序列的X(k)二DFT[xn]屮k值与模拟信号实际频 率f或w的对应关系。 实验原理 gridon; gridon; subplot(212) yn=ifft(Ykz L); stem(nz yn); title (11); xlabel ( k); ylabel『线性卷积YkJ ; gridon; 实验内容 1?已知序列xn={l/2,l,l,l/2},利用FFT计算线性卷积xn*xn,画出图形 clc,clear; ylabel(* xn *); N=4; M=4; L=N+M; n=0:L-l; xn=[l/2zlz1,1/2]; xn=[xn,zeros(1,(L-M))]; Xk=f ft(xn,L); 原序列°0 \ ) 1 1 1 1~IlliillIlli■IIill■II iiii Illi iiii7 7 1 I rI 原序列 °0 \ ) 1 1 1 1~ Illi ill Illi ■II ill ■II iiii Illi iiii 7 7 1 I r I 1 1 i i 1 i i i i a i i i i i a i i i i i a i i 4 4 4 6 3 4 ￡ 0. 2 7 ? 、 1 》 ： 人 1 1 1 1 1 d r 1 1 1 1 S 弋 才 1 7 r $ 1 V 1 1 1 1 序列xn的线性卷积 2 7 0 2 4 6 兰1 2?对如下模拟信号进行谱分析 Xt=cos(8*pi*t)+cos(16*pi*t)+cos(20*pi*t) 选择Fs=64hz,对变换区间N=16,32,64三种情况进行谱分析。 (1)用stem绘制三种情况下的幅频特性曲线，显示三种情况的 谱分辨率F各是多少。 Fs=64; Nl=16; N2=32; N3=64; Fl=Fs/Nl;%Fl=4 F2 二 Fs/N2;%F2=2 F3=Fs/N3;%F3=l %fk=k*F=k/(N*Ts)=k/Tp; k=O:N-l; fl=O:Fl:(Nl-l)*Fl; f2=0:F2:(N2-l)*F2; f3=0:F3:(N3-l)*F3; Ts=l/Fs;%Ts=l/64 nl=O:Nl-l; n2=0:N2-l; n3=0:N3-l; Tpl=l/Fl; Tp2=l/F2; Tp3=l/F3; tl=O:Ts:Tpl; t2=0:Ts:Tp2; t3=0:Ts:Tp3; xnl=cos(8*pi*tl)+cos(16*pi*tl) +cos(20*pi*tl);% 周期二 1/4 fl=l/4,f2=l/8 J3=l/10hz Xkl=fft(xnl/Nl); xn2=cos(8*pi*t2)+cos(16*pi*t2) +cos(20*pi*t2); Xk2=fft(xn2,N2); xn3=cos(8*pi*t3)+cos(16*pi*t3) +cos(20*pi*t3); Xk3=fft(xn3,N3); XaKFl=Ts*Xkl; XaKF2 二 Ts*Xk2; XaKF3=Ts*Xk3; subplot(311); stem(fl,abs(XaKFl]/max(abs(Xa KF1))); titleCN 二 16,Tp 二 0.5 个周期 ZF=4); xlabel(f/hz); ylabef 幅度); subplot(312); stem(f2,abs(XaKF2)/max(abs(Xa KF2))); title(N 二 32,Tp 二 1 个周期 F=2) xlabelCf/hz*); ylabelf 幅度 J; subplot(313); stem(f3,abs(XaKF3]/max(abs(Xa KF3))); title(N 二 64,Tp 二 2 个周期 F=l) xlabelCf/hzJ; ylabelC 幅度)； N=16.Tp=05个周期，F=4幅频特性曲线 1 log 60 70 {~* 5 1 1 1 1~~5 -Q- -G N=32,Tp=1 个周期，F=2 0.5 10 20 50 log 30 40 f/hz N=64Jp=2个周期，F=1 f/hz log 绘制三种情况的功率谱 Fs=64; Nl=16; t2=0:Ts:Tp2; N2=3