DFT_利用DFT分析模拟信号频谱实验报告[实用论文].docVIP

第一篇 实验三 利用DFT分析模拟信号频谱 实验目的： 应用傅里叶变换DFT，分析各种模拟信号x(t)的频谱。 实验原理： 应用离散傅里叶变换DFT，通过MATLAB实现对连续信号频谱分析，以满足理论分析和工程实际的需要。 实验内容： 利用FFT计算信号x(t)=exp(-2x)[u(t)-u(t-2)]的离散谱X(m) DFT计算参数的确定 Tp=2s；Fs=100HZ；由公式F=1/Tp=0.5HZ;N=200;取L=512; 实验程序及频谱图如下： N=200; L=512; fs=100; T=1/fs; ws=2*pi*fs; t=(0:N-1)*T; f=exp((-2)*t); F=fftshift(fft(f,L)); w=(-ws/2+(0:L-1)*ws/L)/2*pi; plot(w,abs(F)); xlabel(w); ylabel(exp((-2*t))*[u(t)-u(t-2)]的幅度谱) 2.利用DFT计算信号x(t)=exp(-2t)u(t)的离散频谱X(m) DFT计算参数的确定 Tp=16s；Fs=16HZ；由公式F=1/Tp=0.06HZ;N=256;取L=256; 实验程序及频谱图如下： N=256; L=256; fs=16; T=1/fs; ws=2*pi*fs; t=(0:N-1)*T; f=exp((-2)*t); F=fftshift(fft(f,L)); w=(-ws/2+(0:L-1)*ws/L)/2*pi; plot(w,abs(F)); xlabel(w); ylabel(exp((-2*t))的幅度谱) 从图中可见，理论频谱与由DFT近似计算的频谱之间存在一定的误差，特别在高频段近似计算的频谱误差更大，这是由于信号x(t)不是限带信号，在时域抽样时产生混频。由于信号x(t)也不是时限信号，由DFT分析频谱时，也存在时域加窗截短造成频谱泄漏。同时，对连续信号抽样时，抽样频率也直接影响频谱的精度。 析淹没在噪声中的信号x(t),例如由50HZ和120HZ的正弦信号以及一个零均值的随机噪声叠加而成的频谱和功率谱。 确定N=50,L=512,Fs=300HZ,其频谱图如下： N=50; L=512; f1=50; f2=120; fs=300; T=1/fs; ws=2*pi*fs; t=(0:N-1)*T; f=sin(2*pi*f1*t)+sin(2*pi*f2*t)+randn(1,N); wh=(hamming(N)); f=f.*wh; F=fftshift(fft(f,L)); w=(-ws/2+(0:L-1)*ws/L)/(2*pi); plot(w,abs(F)); xlabel(w); ylabel(x(t)的幅度谱) 利用Welch法作其功率谱如下： N=512; Nfft=1024; Fs=2*pi; f1=50; f2=120; fs=300; f1=f1/fs; f2=f2/fs; n=0:N-1; xn=sin(2*pi*f1*n)+sin(2*pi*f2*n)+randn(size(n)); L=input(L=); window=boxcar(L); noverlap=L/2; [pxx2,f]=psd(xn,Nfft,Fs,window,noverlap); plot(f,10*log(pxx2));grid axis([0 pi -60 60]) xlabel(Frequency) ylabel(Power Spectral(db)) title([Welch PSD Estimate]) 检测信号的原理及步骤如下： 利用滤波的方法，在频域将有用信号提出，然后计算出时域信号x1(t),用x(t)-x1(t)即可得到噪声信号。 利用Hamming窗和Kaiser窗对长度为0.4s信号 ,其表达式为： x(t)=cos(2*pi*f1*t)+0.75cos(2*pi*f2*t) 其中f1=100HZ,f2=110HZ进行频谱分析。 Hamming窗的频谱如下： 取Fs=250HZ, N=100; L=512; f1=100; f2=110; fs=250; T=1/fs; ws=2*pi*fs; t=(0:N-1)*T; f=cos(2*pi*f1*t)+0.75*cos(2*pi*f2*t); wh=(hamming(N)); f=f.*wh; F=fftshift(fft(f,L)); w=(-ws/2+(0:L-1)*ws/L)/(2*pi); plot(w,abs(F)); xlabel(w); ylabel(x(t)的幅度谱) title([Hamming窗函数]) 通过实验分析可以确定TL=1/△F。 K