第三章 离散傅里叶.ppt

第三章 离散傅立叶变换 主讲 温川雪 3.2 离散傅里叶变换的基本性质 3.2.3 循环卷积定理 循环卷积 频域循环卷积定理 任何有限长序列都可以表示成其共轭对称分量和共轭反对称分量之和。 DFT的共轭对称性 3.3 频率域采样 MATLAB求解程序ep331.m如下： %《数字信号处理(第三版）》第3章例3.3.1程序ep331. % 频域采样理论验证 M=26; N=32; n=0:M;  xa=0:M/2; xb= ceil(M/2)-1:-1:0; xn=［xa, xb］; 　　　　　　　　 %产生M长三角波序列x(n) Xk=fft(xn, 512); 　 %512点FFT［x(n)］ X32k=fft(xn, 32); 　%32点FFT［x(n)］ x32n=ifft(X32k); 　%32点IFFT［X32(k)］得到x32(n) X16k=X32k(1:2:N); %隔点抽取X32k得到X16(k) x16n=ifft(X16k, N/2); %16点IFFT［X16(k)］得到x16(n) 以下绘图部分省略。 重叠相加法计算线性卷积 适用条件：两序列长度相差很大 设序列h(n)长度为N， x(n)为无限长序列。 将x(n)均匀分段， 每段长度取M， 则 于是， h(n)与x(n)的线性卷积可表示为 　【例3.4.1】 假设h(n)=R5(n)，x(n)=［cos(πn/10)　+cos(2πn/5)］u(n)，用重叠相加法计算y(n)=h(n)*x(n)，并画出h(n)、x(n)和y(n)的波形。 解 h(n)的长度为N=5, 对x(n)进行分段，每段长度为M=10。计算h(n)和x(n)的线性卷积的MATLAB程序如下： Lx=41; N=5; M=10; %Lx为信号序列x(n)长度 hn=ones(1, N); hn1=［hn zeros(1, Lx-N)］; 　　　　　%产生h(n)，其后补零是为了绘图好看 n=0:L-1;  xn=cos(pi*n/10)+cos(2*pi*n/5); %产生x(n)的Lx个样值 yn=fftfilt(hn, xn, M); %调用fftfilt用重叠相加法计算卷积 3.4.2 用DFT对信号进行谱分析 用DFT对连续信号进行谱分析 所谓信号的谱分析就是计算信号的傅里叶变换。 用DFT对序列进行谱分析 周期序列，用整数倍周期的长度作FFT 非周期序列，尽量选择信号的观察时间长一些 　27． 选择合适的变换区间长度N， 用DFT对下列信号进行谱分析， 画出幅频特性和相频特性曲线。 （2） x2(n)=sin(0.45πn)sin(0.55πn) 解： x2(n)=sin(0.45πn) sin(0.55πn) =0.5［cos(0.1πn)－cos(πn)］， 前半部分周期为20，后半部分周期是2，所以x2(n)的周期为20， 所以取N=20； 用DFT进行谱分析的误差问题 截断效应 用DFT对模拟信号进行谱分析 按采样定理变为时域离散信号 按周期或非周期序列，选择适当长度作FFT ex327.m 时域循环卷积定理 3.2.4 复共轭序列的DFT 　设x*(n)是x(n)的复共轭序列，长度为N，X(k)=DFT［x(n)］N，则  且X(N)=X(0)。 3.2.5 DFT的共轭对称性 【例3.2.2】 利用DFT的共轭对称性，设计一种高效算法，通过计算一个N点DFT，就可以计算出两个实序列x1(n)和x2(n)的N点DFT。 解：1、直接求解，需作两次N点的DFT. 2、将两个实数序列合并为一个复数序列 对x(n)求DFT得到： 所以，由X(k)可以求得两个实序列x1(n)和x2(n)的N点DFT： 　　12． 已知f(n)=x(n)+jy(n)， x(n)与y(n)均为长度为N的实序列。 设 　　　　　F(k)=DFT［f(n)］N 0≤k≤N－1 (1) 　　 试求X(k)=DFT［x(n)］N， Y(k)=DFT［y(n)］N　　 解: 由DFT的共轭对称性可知 　　　　　　x(n)　　X(k)=Fep(k) 　　　　　　jy(n)　　jY(k)=Fop(k) 　　7． 证明: 若x(n)为实序列， X(k)=DFT［x(n)］N， 则X(k)为共轭对称序列，