重点1.离散频谱（DFT）概念2.DFT性质与计算3.DFT应用.pptVIP

5.Parseval定理 揭示信号时域能量和频域能量之间的关系 6. DFT和z变换的关系 单位圆上的z变换 频域等间隔抽样 因此 DFT是z平面单位圆上的等间隔抽样 6. DFT和z变换的关系 DFT是z平面单位圆上的等间隔抽样 ] [ m X ? ? ? ? IDFT ] [ k x ? ? ? ? ? 变换 z ) ( z X (内插公式) 由序列DFT表示序列z变换 4 利用DFT计算线性卷积 — DFT 应用之一 求解系统响应： 但计算机可以DFT计算循环卷积(根据卷积定理)： 线性卷积和循环卷积之间存在什么联系？ 线性卷积： 卷积后频谱： 频谱的L点抽样(即L点DFT)： L点IDFT（还原原序列） 推导过程 满足 ,则有： 利用DFT计算循环卷积框图( ) f=[1,2,2,1,0,0]; h=[3,2,1,0,0,0]; F=fft(f); H=fft(h); y=ifft(F.*H) % Linear Convolution through DFT p113 x = [1 2 0 1]; h = [2 2 1 1]; % Determine the length of the result of convolution L = length(x)+length(h)-1; % Compute the DFTs by zero-padding XE = fft(x,L); HE = fft(h,L); % Determine the IDFT of the product y1 = ifft(XE.*HE); % Plot the sequence generated by circular convolution % and the error from direct linear convolution k = 0:L-1; subplot(2,1,1); stem(k,real(y1));axis([0 6 0 7]); title(Result of Linear Convolution); xlabel(Time index k);ylabel(Amplitude);grid; y2 = conv(x,h); error = y1-y2; subplot(2,1,2); stem(k,abs(error)); xlabel(Time index k);ylabel(Amplitude); title(Error Magnitude);grid; 利用DFT计算线性卷积 5 利用DFT分析连续非周期信号的频谱 1)连续信号Fourier变换和DFT的关系 DTFT和Fourier变换的关系 DTFT和DFT的关系 --DFT 应用之二 5 利用DFT分析连续非周期信号的频谱 1. 无限长，其频带有限 加窗 抽样 DFT 5 利用DFT分析连续非周期信号的频谱 2. 有限长，其频带无限 抽样 DFT 5 利用DFT分析连续非周期信号的频谱 3. 无限长，其频带无限 加窗 出现三种现象：混叠、泄漏、栅栏 抽样 DFT ] [ m X m 2 / N 1 - N ) ( w j X w 2 s w 2 s w - 在 2 / N m ￡ , ] [ m X 对应于 ) ( w j X N m s w w = 的抽样值 在 1 2 / - ￡ ￡ N m N , ] [ m X 对应于 频率点m和信号模拟频率之间的对应关系 1 混叠现象(aliasing) 避免混叠：（使满足抽样定理） (1)对于带限信号，提高抽样率。 (2)对于非带限信号或高频处无用的信号， 通过预滤波去除高频分量 （由于不满足抽样定理引起） 1 混叠现象(续) 解决方法：减小抽样间隔T，抗混滤波 抗混滤波 抽样间隔T 抽样 DFT 例： 用 DFT 近似计算信号 t e t f - = ) ( 的频谱。 ) ( t f 2 / p t - t 2 / p t ] [ k f k 返回到时域也发生混叠 N=8点抽样 %compute the spectrum of f(t)=exp(-|t|)through DFT %draw the theoretical spectrum w=linspace(0,10,1024); plot(w,2./(w.^2+1)); % sampling points and frequency(rad/s) N=input(抽样点数=); Ws=input(抽样角频率=); Ts=2*pi/Ws; %compute the sampling points