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

实验三 FFT在卷积积分和谱分析中的应用 实验目的 在理论学习的基础上，通过本实验，加深对FFT的理解，熟悉FFT子程序。 熟悉掌握FFT实现两个序列的线性卷积的方法 熟悉应用FFT进行信号谱分析过程中可能出现的问题以便在实际中正确应用FFT. 学习用FFT对连续信号进行谱分析的方法和参数选择原则。 掌握用FFT对连续信号进行谱分析时，由模拟信号采样得到序列的X(k)=DFT[xn]中k值与模拟信号实际频率f或w的对应关系。 实验原理 三．实验内容 1.已知序列xn={1/2,1,1,1/2}，利用FFT计算线性卷积xn*xn,画出图形。 clc,clear; N=4; M=4; L=N+M; n=0:L-1; xn=[1/2,1,1,1/2]; xn=[xn,zeros(1,(L-M))]; Xk=fft(xn,L); Yk=Xk.*Xk; subplot(211) stem(n,xn) title(原序列xn); xlabel(n) ylabel(xn); grid on; subplot(212) yn=ifft(Yk,L); stem(n,yn); title(); xlabel(k); ylabel(线性卷积Yk); grid on; 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; N1=16; N2=32; N3=64; F1=Fs/N1;%F1=4 F2=Fs/N2;%F2=2 F3=Fs/N3;%F3=1 %fk=k*F=k/(N*Ts)=k/Tp; k=0:N-1; f1=0:F1:(N1-1)*F1; f2=0:F2:(N2-1)*F2; f3=0:F3:(N3-1)*F3; Ts=1/Fs;%Ts=1/64 n1=0:N1-1; n2=0:N2-1; n3=0:N3-1; Tp1=1/F1; Tp2=1/F2; Tp3=1/F3; t1=0:Ts:Tp1; t2=0:Ts:Tp2; t3=0:Ts:Tp3; xn1=cos(8*pi*t1)+cos(16*pi*t1)+cos(20*pi*t1);%周期= 1/4 f1=1/4,f2=1/8 ,f3=1/10hz Xk1=fft(xn1,N1); 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); XaKF1=Ts*Xk1; XaKF2=Ts*Xk2; XaKF3=Ts*Xk3; subplot(311); stem(f1,abs(XaKF1)/max(abs(XaKF1))); title(N=16,Tp=0.5个周期? F=4); xlabel(f/hz); ylabel(幅度); subplot(312); stem(f2,abs(XaKF2)/max(abs(XaKF2))); title(N=32,Tp=1个周期 F=2) xlabel(f/hz); ylabel(幅度); subplot(313); stem(f3,abs(XaKF3)/max(abs(XaKF3))); title(N=64,Tp=2个周期 F=1) xlabel(f/hz); ylabel(幅度); 绘制三种情况的功率谱 Fs=64; N1=16; N2=32; N3=64; F1=Fs/N1; F2=Fs/N2; F3=Fs/N3; f1=0:F1:(N1-1)*F1; f2=0:F2:(N2-1)*F2; f3=0:F3:(N3-1)*F3; Ts=1/Fs; n1=0:N1-1; n2=0:N2-1; n3=0:N3-1; Tp1=1/F1; Tp2=1/F2; Tp3=1/F3; t1=0:Ts:Tp1; t2=0:Ts:Tp2; t3=0:Ts:Tp3; xn1=cos(8*pi*t1)+cos(16*pi*t1)+cos(20*pi*t1); Xk1=fft(xn1,N1); 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); subplot(311); plot(f1,Xk1.*conj(