离散傅里叶变换(DFT).ppt

*Xk=fft(xn,L); %L点FFT[x(n)]Hk=fft(hn,L); %L点FFT[h(n)]Yk=Xk.*Hk;%频域相乘得Y(k)yn=ifft(Yk,L);%L点IFFT得到卷积结果y(n)toc%快速卷积计时结束subplot(2,2,1),stem(nx,xn,.);ylabel(x(n))subplot(2,2,2),stem(nh,hn,.);ylabel(h(n))subplot(2,1,2);ny=1:L;stem(ny,real(yn),.);ylabel(y(n))tic,yn=conv(xn,hn); %直接调用函数conv计算卷积与快速卷积比较toc**01DFT(实际中用FFT计算)可用来对连续信号和数字信号进行谱分析。用DFT进行谱分析的误差问题02栅栏效应。截断效应。根据傅里叶变换的频域卷积定理有混叠现象。*其中幅度谱RN(ω)~ω曲线如图所示(RN（ω）以2π为周期，只画低频部分)。图中，|ω|2π/N的部分称为主瓣，其余部分称为旁瓣。例如，x(n)=cos(ω0n),ω0=π/4其频谱为*图3.4.11矩形窗函数的幅度谱*单击此处添加正文，文字是您思想的提炼，为了演示发布的良好效果，请言简意赅地阐述您的观点。泄露图3.4.12加矩形窗前后的频谱谱间干扰*在计算机上用DFT对模拟信号进行谱分析时，只能以有限大的采样频率fs对模拟信号采样有限点样本序列（等价于截取模拟信号一段进行采样）作DFT变换，得到模拟信号的近似频谱。其误差主要来自以下因素：截断效应（频谱泄露和谱间干扰）频谱混叠失真因素①使谱分辨率（能分辨开的两根谱线间的最小间距）降低，并产生谱间干扰；因素②使折叠频率（fs/2）附近的频谱产生较大失真。例7.19用DFT求连续信号频谱*加大截取长度Tp可提高频率分辨率；选择合适的窗函数可降低谱间干扰；而频谱混叠失真要通过提高采样频率fs和（或）预滤波（在采样之前滤除折叠频率以外的频率成分）来改善。编写程序q719.m验证截断效应及加窗的改善作用，先选取以下参数：?采样频率fs=400Hz，T=1/fs?采样信号序列?对x(n)作4096点DFT作为xa(t)的近似频谱Xa(jf)。例7.19用DFT求连续信号频谱*%第七章例7.19程序q719%用DFT作谱分析clear;closeallfs=400;T=1/fs; %采样频率为400HzTp=0.04;N=Tp*fs; %采样点数NN1=[N,4*N,8*N]; %设定三种截取长度供调用st=[|X1(jf)|;|X4(jf)|;|X8(jf)|]; %设定三种标注语句供调用%矩形窗截断form=1:3n=1:N1(m);xn=cos(200*pi*n*T)+sin(100*pi*n*T)+cos(50*pi*n*T);%产生采样序列x(n)Xk=fft(xn,4096);%4096点DFT，用FFT实现fk=[0:4095]/4096/T;subplot(3,2,2*m-1)plot(fk,abs(Xk)/max(abs(Xk)));ylabel(st(m,:))ifm==1title(矩形窗截取);endend*%加hamming窗改善谱间干扰form=1:3n=1:N1(m);wn=hamming(N1(m));%调用工具箱函数hamming产生N长hamming窗序列wnxn=(cos(200*pi*n*T)+sin(100*pi*n*T)+cos(50*pi*n*T)).*wn;Xk=fft(xn,4096);%4096点DFT，用FFT实现fk=[0:4095]/4096/T;subplot(3,2,2*m)plot(fk,abs(Xk)/max(abs(Xk)));ylabel(st(m,:))ifm==1title(Hamming窗截取);endend*图7.19*例7.19用DFT求连续信号频谱如图7.19所示。图中X1(jf)，X4(jf)和X8(jf)分别表示Tp=0.04s，0.04*4s和0.04*8s时的谱分析结果。由图可见，由于截断使原频谱中的单频谱线展宽（又称之为泄漏），Tp越大泄漏越小，频率分辨率越高。Tp=0.04s时，25Hz与50Hz两根谱线已分辨不清了。所以实际谱分析的截取时间