![数字信号处理及MATLAB实现](https://wfqqreader-1252317822.image.myqcloud.com/cover/677/41398677/b_41398677.jpg)
第三节 DFT的性质
本章第一节介绍了时间离散信号傅里叶变换的很多性质,这一节将介绍DFT的一些重要性质。读者可以发现DFT的性质与时间离散信号的傅里叶变换性质非常相似,这些性质归纳在表3-3中。
表3-3 DFT的性质
![](https://epubservercos.yuewen.com/688291/21570843601309506/epubprivate/OEBPS/Images/61_02.jpg?sign=1739337248-5ZfjsbJCWnVlL96TPgCTU388TSWvX0GD-0-3061a5502720575791abe2f65d27c288)
一、周期性
由于,k、m为整数,N为自然数,所以式(3-33)满足:
![](https://epubservercos.yuewen.com/688291/21570843601309506/epubprivate/OEBPS/Images/62_02.jpg?sign=1739337248-C39BJhhJOYPcDEWeMkAKsvriRPGBzTkQ-0-1d99f83ff73ad39f21e966ea7245e2ca)
即X(k)是周期为N的周期函数。通常研究其中的一个周期0≤k≤N-1,并把它称为主值区间,于是有公式X(k)=X(k+mN)RN(k)。
二、线性性质
若序列x1(n)和x2(n)是两个有限长的序列,长度分别为N1和N2,它们的离散傅里叶变换为X1(k)=DFT[x1(n)]和X2(k)=DFT[x2(n)],那么序列y(n)=ax1(n)+bx2(n)(a,b是任意常数)的离散傅里叶变换Y(k)可表示为
![](https://epubservercos.yuewen.com/688291/21570843601309506/epubprivate/OEBPS/Images/62_03.jpg?sign=1739337248-GdJ4aAgH451AQdyQRtSAaR0E6Lz7kHkR-0-1ea23a9c271eeafb68286eaefdf1ae3d)
三、循环移位
一个长度为N的有限长序列x(n),在区间[0,N-1]内有非零值。将该序列以N为周期作延拓,得到的周期序列用x((n))N或表示,它与原序列x(n)的关系可表示为
![](https://epubservercos.yuewen.com/688291/21570843601309506/epubprivate/OEBPS/Images/62_05.jpg?sign=1739337248-CMp32npQh3BxnOwAqDqF8jYgBeULcENY-0-5b1e4a4ba1190c2d61c8e9e4f00ca376)
将向左平移m个单位,然后取主值区间[0,N-1]的值,可表示成
![](https://epubservercos.yuewen.com/688291/21570843601309506/epubprivate/OEBPS/Images/62_07.jpg?sign=1739337248-vcSegsX2waE3QRgZVnea4cwal4gTh1x5-0-e108028ea27d9f3a4da71d302ee3688a)
y(n)的产生过程可理解为将原序列x(n)向左平移m个单位,移出区间[0,N-1]的序列从坐标轴右边循环回来补充到[N-M,N-1]区间内,这种移位称作循环移位。y(n)的长度仍然为N,且区间范围与原序列一样。图3-4表示一个长度为5的序列循环移位的过程。图3-4a表示一个长度为5的有限长序列,图3-4b表示将该序列以5为周期进行周期延拓,图3-4c表示将图3-4b向左平移两个单位,图3-4d是取主值区间的序列。对照图3-4c发现,向左移出的两个点又从右边补充进来。
下面用MATLAB程序验证循环移位性质。
例3-10 已知一个8点的序列x(n)=10(0.5)n,0≤n≤7,绘出x((n-3))10。
解:MATLAB脚本如下:
![](https://epubservercos.yuewen.com/688291/21570843601309506/epubprivate/OEBPS/Images/62_08.jpg?sign=1739337248-2FV1TkmSmDMwwZx0XGx8uyZqYJ8M9dkq-0-7b57fdaef93e6a0917a3da38e617a8ac)
![](https://epubservercos.yuewen.com/688291/21570843601309506/epubprivate/OEBPS/Images/63_01.jpg?sign=1739337248-BNYHAW5iDVwXdesFeb6089CfqRcQUD9A-0-6521635ae2620ca1f1a166c8321ebc86)
![](https://epubservercos.yuewen.com/688291/21570843601309506/epubprivate/OEBPS/Images/63_02.jpg?sign=1739337248-h3SvEUIkMVWPYI9nePWmWkATKff5TyJg-0-3ef9dcfa4f127ad104c5ad5b2cd77271)
图3-4 循环移位示意图
其中函数cirshftt的程序如下:
![](https://epubservercos.yuewen.com/688291/21570843601309506/epubprivate/OEBPS/Images/63_03.jpg?sign=1739337248-xqoCJNEMMvSqG6XtfbKUaE7OJ23bVZfx-0-9ad7e9f00f66fb4dc785c04b12a2fc85)
运行结果如图3-5所示。
![](https://epubservercos.yuewen.com/688291/21570843601309506/epubprivate/OEBPS/Images/64_01.jpg?sign=1739337248-78fjz1r9ASxyrAC5eey1WBXgxVYUi30v-0-d4d221e369e17c3c7adc822db651ab7b)
图3-5 例3-10运行结果
四、时移性质
已知x(n)是长度为N的有限长序列,y(n)是x(n)的循环移位,表示成
y(n)=x((n+m))NRN(n)
y(n)的离散傅里叶变换Y(k)为
![](https://epubservercos.yuewen.com/688291/21570843601309506/epubprivate/OEBPS/Images/64_02.jpg?sign=1739337248-VVr9q7sHbfmucGwxJPskF4ANEg3KmpI9-0-a71fd118fe16a79ed6c751bb2a703176)
其中
X(k)=DFT[x(n)]N,0≤k≤N-1
证明:
![](https://epubservercos.yuewen.com/688291/21570843601309506/epubprivate/OEBPS/Images/64_03.jpg?sign=1739337248-g67cKjg0Km98oEeS6zer2bna7hLYBHAl-0-23311ee2e197f54a46c47a7a0018ba33)
令n+m=n′,有
![](https://epubservercos.yuewen.com/688291/21570843601309506/epubprivate/OEBPS/Images/64_04.jpg?sign=1739337248-4GPIu3TxqhRLjFHsQLxpxhpz7wd1b4xe-0-d2d568501b57d7ffbb965d29cef49a95)
由于上式中求和项是以N为周期,对其在任一周期上的求和结果相同。将上式的求和区间改在主值区间内,则
![](https://epubservercos.yuewen.com/688291/21570843601309506/epubprivate/OEBPS/Images/64_06.jpg?sign=1739337248-LpvsH8CPMwclT1IK00DJZnUmUDAZnZzr-0-e9b30edfb68b7b71625c5bd8720deb83)
五、频移性质
如果
![](https://epubservercos.yuewen.com/688291/21570843601309506/epubprivate/OEBPS/Images/64_07.jpg?sign=1739337248-uB56xla16a3Wu1n4L2b3rGPjs7oEZDoL-0-9fcee448213fc13c7a8a3f1db6bd74ed)
则
![](https://epubservercos.yuewen.com/688291/21570843601309506/epubprivate/OEBPS/Images/65_01.jpg?sign=1739337248-w0h8f8VhRuHuO4iGKR40sy4w7z55CsZ9-0-3edaaa1f44f8f04813ff6ef2deea0684)
式(3-49)的证明方法与时移性质类似。
六、复共轭序列的DFT
已知x*(n)是x(n)的复共轭序列,x(n)的N点傅里叶变换为X(k),那么x*(n)的DFT可表示为
![](https://epubservercos.yuewen.com/688291/21570843601309506/epubprivate/OEBPS/Images/65_02.jpg?sign=1739337248-AUcDAr43yyUc48T4BDiQDwrLuholAUSu-0-fa64316a49058d4c33564b7c78add0df)
证明:
![](https://epubservercos.yuewen.com/688291/21570843601309506/epubprivate/OEBPS/Images/65_03.jpg?sign=1739337248-fTIzUU9AP1rTDw5J3kpZsGKFszlPX2Fd-0-0ee90582f1b1483cd7d9b7df9ec8efbb)
由于X(k)的隐含周期性,有X(N)=X(0)。
用同样的方法可以证明:
![](https://epubservercos.yuewen.com/688291/21570843601309506/epubprivate/OEBPS/Images/65_04.jpg?sign=1739337248-DAjPR6FlCLxeoKkJbhGOCt5PHugPIXGh-0-83bbc95b1c6cd2a731fd644ddc6a8142)
七、DFT的共轭对称性质
本章第一节介绍了无限长时间离散信号在时域和频域的对称性质,它们都是关于坐标原点共轭对称或反对称。类似地,有限长序列在时域和频域也具有对称性质,但它们是关于n=N/2对称。
(一)有限长序列的共轭对称性和共轭反对称性
为了区别本章第一节里的傅里叶变换的对称性,用xep(n)和xop(n)分别表示有限长共轭对称序列和有限长共轭反对称序列。它们满足下面的关系式:
![](https://epubservercos.yuewen.com/688291/21570843601309506/epubprivate/OEBPS/Images/65_05.jpg?sign=1739337248-p7Lrce40UZemFYFIwmHtBsghJaTJlhLI-0-7859f3e1e8da9a952c62cf9099101147)
当N为偶数时,用N/2-n代替上面式中的n,可得到
![](https://epubservercos.yuewen.com/688291/21570843601309506/epubprivate/OEBPS/Images/65_06.jpg?sign=1739337248-2AyI0pCfUPYJr3hsdqOmsQTb6Lxb9DCn-0-3aba47084044e59684c6de75b6512061)
以上两个式子说明有限长共轭对称序列关于n=N/2对称。类似于前面介绍的任何一个无限长序列可以写成共轭对称序列和反对称序列相加的形式,有限长序列也可以写成共轭对称分量和共轭反对称分量之和的形式,即
![](https://epubservercos.yuewen.com/688291/21570843601309506/epubprivate/OEBPS/Images/65_07.jpg?sign=1739337248-2AUnvzqp6ghv0LFi3PWzNulkcIsPytLo-0-9faca77092ba6ea67610e848d321ccc0)
将上式中的n用N-n代替,并取共轭,可以得到下面的式子
![](https://epubservercos.yuewen.com/688291/21570843601309506/epubprivate/OEBPS/Images/65_08.jpg?sign=1739337248-bipokCh1GMvD7Qg6gl7dqLIhxYqLvjbe-0-4e8c60d1cc93442c13fb8b0bade94ef7)
由式(3-54)和式(3-55)相加或相减可得
![](https://epubservercos.yuewen.com/688291/21570843601309506/epubprivate/OEBPS/Images/66_01.jpg?sign=1739337248-vGzJVcyztlb5Nv12HbseUJsIqAi62H9D-0-177350d1fa305304d8a68fdf6b41b3a0)
(二)DFT的共轭对称性
与时域对称性质的表达方式类似。在频域,用Xep(k)和Xop(k)分别表示离散傅里叶变换的共轭对称性和共轭反对称性。它们满足下面的公式:
![](https://epubservercos.yuewen.com/688291/21570843601309506/epubprivate/OEBPS/Images/66_02.jpg?sign=1739337248-5F4qEfG3eMIQresoihZy0De9mwWddF5x-0-802e40837bba1167ab93b4dc075ec348)
如果序列x(n)的DFT用X(k)表示,那么X(k)可表示成
X(k)=Xep(k)+Xop(k)
与时域的推导方法类似,可得到下面两个公式:
![](https://epubservercos.yuewen.com/688291/21570843601309506/epubprivate/OEBPS/Images/66_03.jpg?sign=1739337248-2VrLSuaQENOv4wUhOjbiyfPUyPYXEU5u-0-b6ca22a3edbb3ca42394b40290eaa04f)
(1)将x(n)写成实部和虚部相加的形式,即
![](https://epubservercos.yuewen.com/688291/21570843601309506/epubprivate/OEBPS/Images/66_04.jpg?sign=1739337248-S5e6wCbrKNo1NGPXoM3sAovaaVdxPTWh-0-962872b11d3ec539dffd190a2d66d754)
式中
![](https://epubservercos.yuewen.com/688291/21570843601309506/epubprivate/OEBPS/Images/66_05.jpg?sign=1739337248-UlbICWeBdZA49VD6EmrEt0jUnRoTjDio-0-27e682f28f6ca8bda9a1ed5777d9751f)
对以上两式左右两边进行DFT,有
![](https://epubservercos.yuewen.com/688291/21570843601309506/epubprivate/OEBPS/Images/66_06.jpg?sign=1739337248-jJrdFwQ5yPrBqHBhunW2DJrxbyfuZNUP-0-a2488a9811f9a80c0d3f77e4ba54294b)
以上两式说明,序列x(n)的实部xr(n)的傅里叶变换具有共轭对称的性质,虚部xi(n)和j的乘积的傅里叶变换具有共轭反对称的性质。
(2)将x(n)写成共轭对称分量和共轭反对称分量相加的形式,即
x(n)=xep(n)+xop(n)
对xep(n)和xop(n)分别进行DFT,可得
![](https://epubservercos.yuewen.com/688291/21570843601309506/epubprivate/OEBPS/Images/66_07.jpg?sign=1739337248-6ZMQ8kfHdx5vQRqahggZglZa3ewWRta4-0-65e011861a4f2348f5b59fbec1d5d932)
以上两式说明,如果一个有限长序列写成共轭对称分量和反对称分量相加的形式,其共轭对称分量xep(n)的离散傅里叶变换是原序列x(n)的离散傅里叶变换的实部,其共轭反对称分量xop(n)的离散傅里叶变换是原序列x(n)的离散傅里叶变换的虚部乘以j。
如果x(n)是实序列,那么x(n)=x*(n),即序列只有实部。因此,可得
![](https://epubservercos.yuewen.com/688291/21570843601309506/epubprivate/OEBPS/Images/67_01.jpg?sign=1739337248-cLn0ixkoW1cwedVXAe2gYO1FmvOVNx6p-0-3d9bd187033db82e5ff97cc50103ac10)
下面用一个MATLAB的例子说明序列的对称性质。
例3-11 设x(n)=10(0.5)n,0≤n≤7,绘出x(n)的共轭对称分量xep(n)和共轭反对称分量xop(n)的图形,验证式(3-63)和式(3-64)。
解:MATLAB参考脚本如下:
![](https://epubservercos.yuewen.com/688291/21570843601309506/epubprivate/OEBPS/Images/67_02.jpg?sign=1739337248-ikt0jI4vDxYzQRSoHigvyY1YPJAzmSaq-0-ebbce4bf9a50826ca9211cd358d542d4)
![](https://epubservercos.yuewen.com/688291/21570843601309506/epubprivate/OEBPS/Images/68_01.jpg?sign=1739337248-F41TTDJuwa6gEcfqN5rzxD6slsTexEca-0-ffc1e477428e4b1fdccc80c91dffe82a)
运行结果如图3-6所示。
![](https://epubservercos.yuewen.com/688291/21570843601309506/epubprivate/OEBPS/Images/68_02.jpg?sign=1739337248-SesYDO6xzLBkLMsOb6lsxVVbytdaqhTn-0-5332e26c464025b970e0268f0a3298df)
图3-6 例3-11运行结果
由图3-6b可知,共轭对称分量的DFT等于原序列DFT的实部,共轭反对称分量的DFT等于原序列DFT的虚部乘以j。