还剩3页未读,继续阅读
文本内容:
数字信号处理实验五用FFT做频谱分析实验目的1通过本实验,加深对DTFT和IDFT以及DFT和FFT的理解,熟悉FFT子程序2熟悉应用FFT对典型信号进行频谱分析的方法3了解应用FFT进行信号频谱分析过程中可能出现的问题,以便在实际中正确应用FFT4熟悉应用FFT实现两个序列线性卷积的方法5初步了解用周期图法做随机信号频谱分析的方法实验内容1已知有限长序列xn=[
10.
500.
5110.50]要求:用FFT求该序列的DFT、IDFT的图形程序如下:xn=[
10.
500.
5110.50];n=lengthxn;k=0:n-1;subplot221;stemkxnk.;titlexn;Xk=fftxnn;subplot212;stemkabsXk;titleXk=DFTxn;xn1=ifftXkn;subplot222;stemkxn1;titlexn=IDFTXk;波形如下假设采样频率Fs=20Hz序列长度N分别取
8、32和64,用FFT计算幅度谱和相位谱程序如下:clear;closeallfs=20;T=1/fs;N=
[83264];form=1:3x=[
10.
500.
5110.50];x1=fftxNm;x2=ifftxNm;subplot322*m-1;stem[0:Nm-1]absx1o;title幅度谱;subplot322*m;stem[0:Nm-1]absx2o;title相位谱;end波形如下:2用FFT计算下面连续信号的频谱,并观察选择不同的采样周期Ts和序列长度N值对频谱特性的影响:程序如下:clear;closeallfs=4;T=1/fs;Tp=4;N=Tp*fs;N1=[N4*N8*N];T1=[T2*T4*T];form=1:3n=1:N1m;x1=exp-
0.01*T;x2=sin2*n*T+sin
2.1*n*T+sin
2.2*n*T;x=x1*x2;xk=fftx51fk=[0:50]/51/Tsubplot322*m-1;stemfkabsxk/maxabsxkk.;xlabelT1m;ylabelN1m;ifm==1title序列长度endendform=1:3n=1:14;x1=exp-
0.01*T1m;x2=sin2*n*T1m+sin
2.1*n*T1m+sin
2.2*n*T1m;x=x1*x2;xk=fftx51fk=[0:50]/51/Tsubplot322*m;stemfkabsxk/maxabsxkk.;xlabelT1m;ylabelN1m;ifm==1title采样周期endend波形如下:。