![数字信号处理及MATLAB实现](https://wfqqreader-1252317822.image.myqcloud.com/cover/677/41398677/b_41398677.jpg)
第一节 离散信号的傅里叶变换
一、定义
连续时间非周期信号的频域分析,是对其进行傅里叶变换,同样,离散时间非周期信号(采样信号)的频域分析,也可进行傅里叶变换。
设原信号为x(t),采样信号为x(n)=x(t)|t=nT,则x(n)的傅里叶变换定义为
![](https://epubservercos.yuewen.com/688291/21570843601309506/epubprivate/OEBPS/Images/48_01.jpg?sign=1738905539-zIsdUTzYZWRUIOKokvbsc05MmBHW9FXp-0-c87d96fde4b744bb31a7fe60ad1e046f)
式(3-1)成立的前提条件是X(ejω)绝对可和,即|X(ejω)|<∞。
因为
所以,|X(ejω)|<∞等价于
其傅里叶逆变换为
![](https://epubservercos.yuewen.com/688291/21570843601309506/epubprivate/OEBPS/Images/48_04.jpg?sign=1738905539-2UsNFDptEZWJ2PJiyZ5vhE7YAR5gDDIj-0-d1676e840c6d84beccfdb89f03ca41e4)
式(3-1)和式(3-2)组成了序列x(n)的傅里叶变换对。表3-1列举了一些常用序列的傅里叶变换,这里仅举几例说明。
表3-1 常用离散信号的傅里叶变换
![](https://epubservercos.yuewen.com/688291/21570843601309506/epubprivate/OEBPS/Images/48_05.jpg?sign=1738905539-NzXgtT8iydeaB3qxbGCPVFAHA8s7vHqy-0-833a665c7ca80dd357d3e05df2bb9333)
(续)
![](https://epubservercos.yuewen.com/688291/21570843601309506/epubprivate/OEBPS/Images/49_01.jpg?sign=1738905539-rsoFWy6qdV5dj3hho4ikMvPVEMssVEJS-0-05796e7907c13758b2166a711096faf9)
例3-1 求单位采样序列δ(n)的傅里叶变换。
解:
例3-2 求指数序列anu(n)的傅里叶变换,其中|a|<1。
解:
二、性质
离散时间信号的傅里叶变换性质有很多,这里对它们作简单介绍,并归纳在表3-2中。
(一)周期性
X(ejω)具有隐含的周期性,通过式(3-3)可以证明。
![](https://epubservercos.yuewen.com/688291/21570843601309506/epubprivate/OEBPS/Images/49_04.jpg?sign=1738905539-GdffjWade1ah5VwjlzpIMR2Iq92iXq0U-0-92fe4da9a2cb6d2a2baabdbce64d32ed)
利用公式e-j2πk=1,可证明X(ejω)是关于ω=2πk的周期函数。通常对X(ejω)的研究只取一个周期内的数据,把ω∈[-π,π]称作主值区间。
表3-2 离散时间信号傅里叶变换的性质
![](https://epubservercos.yuewen.com/688291/21570843601309506/epubprivate/OEBPS/Images/49_05.jpg?sign=1738905539-O9CXX39Cn1xbezgHIM8m5FdoEE41ahWK-0-9ccc627dcf779b7103504b97b6cf3ace)
(续)
![](https://epubservercos.yuewen.com/688291/21570843601309506/epubprivate/OEBPS/Images/50_01.jpg?sign=1738905539-xyuIOGxQn5POCxTZzXvokI2wXzIEfUYl-0-ede1f860e78a46e3da2363ab819c319e)
(二)线性性质
若序列x1(n)和x2(n)的傅里叶变换分别为X1(ejω)=FT[x1(n)]和X2(ejω)=FT[x2(n)],则序列x(n)=ax1(n)+bx2(n)(a,b是任意常数)的傅里叶变换X(ejω)等于a X1(ejω)+b X2(ejω),即
![](https://epubservercos.yuewen.com/688291/21570843601309506/epubprivate/OEBPS/Images/50_02.jpg?sign=1738905539-ILDSUueRK7rUlVOhDCUavUvvvx2mTj3k-0-d8a0cda8dcc0414c8197a7f6c455a4ed)
证明略。
(三)时间反转定理
若y(n)=x(-n),则y(n)的傅里叶变换为X(e-jω),即
![](https://epubservercos.yuewen.com/688291/21570843601309506/epubprivate/OEBPS/Images/50_03.jpg?sign=1738905539-bVAiXNlFpoFbaoarYzMLI06mBflelMeD-0-480e869a308338df67c142b41078066a)
证明如下:
![](https://epubservercos.yuewen.com/688291/21570843601309506/epubprivate/OEBPS/Images/50_04.jpg?sign=1738905539-A53m2ZTB0JsAGCEXZEizSlRiPKvC4wjy-0-3b0ac843ea2dc8b50dc31a6b3ac47e70)
(四)时移定理
延时序列y(n)=x(n-n0)的傅里叶变换为,n0为整数,即
![](https://epubservercos.yuewen.com/688291/21570843601309506/epubprivate/OEBPS/Images/50_06.jpg?sign=1738905539-7gUfeJvZDfB9Ggpow3JAFSKJr6ortV4A-0-112554864c0586368a5f63d39d743d70)
证明如下:
![](https://epubservercos.yuewen.com/688291/21570843601309506/epubprivate/OEBPS/Images/50_07.jpg?sign=1738905539-fBkh78COuLAvwUdxmi4OCLfytKY76zQA-0-6877f25a20e8dbd6babd03ee29e74d08)
例3-3 求序列y(n)=anu(n)-anu(n-M)的傅里叶变换。
解:y(n)=anu(n)-anu(n-M)=anu(n)-aM·an-Mu(n-M)
查表3-1可知anu(n)的傅里叶变换为,又由时移定理可得an-Mu(n-M)的傅里叶变换为
。利用线性性质,y(n)的傅里叶变换为
![](https://epubservercos.yuewen.com/688291/21570843601309506/epubprivate/OEBPS/Images/51_01.jpg?sign=1738905539-rnCU1iouZkiiclAiPX3WcP22HvzH3eLo-0-968216314521a8def7020645ec6a3ebf)
(五)频移定理
序列的傅里叶变换为
,即
![](https://epubservercos.yuewen.com/688291/21570843601309506/epubprivate/OEBPS/Images/51_04.jpg?sign=1738905539-kBtj1Km5tzGHmtPT5JXYygzdNqrhWUzn-0-8406914b218a2d26dfe2fbc3e459c95f)
证明方法可参考时移定理。
例3-4 设x(n)=cos(πn/2),y(n)=ejπ/(4n)x(n),用MATLAB程序验证频移定理。
解:MATLAB参考程序如下:
![](https://epubservercos.yuewen.com/688291/21570843601309506/epubprivate/OEBPS/Images/51_05.jpg?sign=1738905539-lyx0MvyT6jQ1MLtqj6pFMLMoi7evqHCj-0-e8e48d51d928208a36c156a4f87ff3af)
运行结果如图3-1所示。
由图3-1中的幅度和相位图可知,y(n)的傅里叶变换相对于x(n)的傅里叶变换向右平移了π/4,由此证明了频移定理。
例3-5 求序列y(n)=(-1)nanu(n)的傅里叶变换,其中a<1。
解:可将序列y(n)变形为y(n)=ejπnx(n)的形式,其中x(n)=anu(n)。由例3-2的结论,再根据频移定理,y(n)的傅里叶变换为
![](https://epubservercos.yuewen.com/688291/21570843601309506/epubprivate/OEBPS/Images/52_01.jpg?sign=1738905539-c69htfneDBSBB10a85PL0mGmsdUAba8s-0-23c4b0b5c6ba1593db3880e71773a9bf)
图3-1 例3-4的运行结果
![](https://epubservercos.yuewen.com/688291/21570843601309506/epubprivate/OEBPS/Images/52_02.jpg?sign=1738905539-p2gxgeolGTvX4NfdF4NXgi5h3xlhrN2z-0-b4d2601b13fd3c7d455c991793c8bf72)
(六)频域微分定理
序列y(n)=nx(n)的傅里叶变换为,即
![](https://epubservercos.yuewen.com/688291/21570843601309506/epubprivate/OEBPS/Images/52_04.jpg?sign=1738905539-i5KbEPlvtfh3VA93nTdYmIkkwQ0tnghe-0-c5eff5c6161983045a052fe5d1b52aee)
证明略。
例3-6 求序列y(n)=nanu(n)的傅里叶变换。
解:根据频域微分定理可得
![](https://epubservercos.yuewen.com/688291/21570843601309506/epubprivate/OEBPS/Images/52_05.jpg?sign=1738905539-EiDVoqOrS3twfvQKPRsaeoMDFop6yGzq-0-ee4d369347d42005d96ce6329b33751a)
(七)卷积定理
设y(n)=x1(n)*x2(n),则Y(ejω)可表示成
![](https://epubservercos.yuewen.com/688291/21570843601309506/epubprivate/OEBPS/Images/52_06.jpg?sign=1738905539-rHzGl11U1tNRMkFDAMHJbd1t9CfIaL2z-0-6ddd2c32f2808f91b02cc2a1399140c0)
证明:由卷积的定义可知:
![](https://epubservercos.yuewen.com/688291/21570843601309506/epubprivate/OEBPS/Images/52_07.jpg?sign=1738905539-0E4SepA52EDWzH10ltyGFjfcVs5ErW8J-0-61a3f5e4558b15dee0c3546be5b6549c)
对上式两边进行傅里叶变换,得
![](https://epubservercos.yuewen.com/688291/21570843601309506/epubprivate/OEBPS/Images/52_08.jpg?sign=1738905539-tyzJPcwNNBKzWmbaOZ915IfPJf5k39u5-0-bb0cd773963c54d7681a81a2ec984fc3)
将k=n-m代入上式,得
![](https://epubservercos.yuewen.com/688291/21570843601309506/epubprivate/OEBPS/Images/53_01.jpg?sign=1738905539-NtwEv79xab1U2BmG3xvUXOAprqMvucs6-0-c816654b293a448990ab0625c799caf8)
(八)调制定理
设y(n)=x1(n)·x2(n),则Y(ejω)可表示成
![](https://epubservercos.yuewen.com/688291/21570843601309506/epubprivate/OEBPS/Images/53_02.jpg?sign=1738905539-4f5FIq5vmH7Zkg2bbrGXI1wsehSAh4iV-0-9dad91627c19b26e3e0ee1ac8d5605f1)
证明略。
(九)帕塞瓦尔定理
![](https://epubservercos.yuewen.com/688291/21570843601309506/epubprivate/OEBPS/Images/53_03.jpg?sign=1738905539-Y6QU1sEafgzcDW9ndtQUQT1bxQG8z1OU-0-1f1b8f3eae18d4528bbf6dc9b83c5ecc)
证明略。
(十)对称性
在学习对称性之前,先介绍共轭对称序列和共轭反对称序列的定义。
设序列xe(n)满足下列表达式:
![](https://epubservercos.yuewen.com/688291/21570843601309506/epubprivate/OEBPS/Images/53_04.jpg?sign=1738905539-WicmrqW1LARVHlutA85yU2ICZTB9G5zC-0-5e578d5d6cf8e9a2e9b25d875c232f41)
x e(n)称作共轭对称序列。如果将其写成实部与虚部相加的形式,即
![](https://epubservercos.yuewen.com/688291/21570843601309506/epubprivate/OEBPS/Images/53_05.jpg?sign=1738905539-fkHGrGXwAsAklW7eObDb6ZsBf8a6QAZj-0-1ebe9dcda59b1b14c6edf7b9f1bf14e6)
将式(3-13)中的n用-n代替,并取共轭,得
![](https://epubservercos.yuewen.com/688291/21570843601309506/epubprivate/OEBPS/Images/53_06.jpg?sign=1738905539-YlI8Bdg1RnDIgqbVRuU6o3hkSQM8eqSM-0-37b48619c745897e109e1a4663ca2f73)
将式(3-13)和式(3-14)代入式(3-12)中,得
![](https://epubservercos.yuewen.com/688291/21570843601309506/epubprivate/OEBPS/Images/53_07.jpg?sign=1738905539-neW4OH2JPeHxlPnPkmRqPT6DyHahNUlo-0-ec18c9960dae0ef2dd42ce5370f85963)
以上两式表明,共轭对称序列的实部为偶函数,虚部为奇函数。类似地,可得出共轭反对称序列[用xo(n)表示]的定义及性质。
![](https://epubservercos.yuewen.com/688291/21570843601309506/epubprivate/OEBPS/Images/53_08.jpg?sign=1738905539-LfmwadIQtzVdTZgYhFyfniOY9NRLGRzG-0-0c226b7f836f6cc4a90b74366470ffad)
满足式(3-17)的序列称为共轭反对称序列。式(3-18)和式(3-19)表明,共轭反对称序列的实部为奇函数,虚部为偶函数,这与共轭对称序列正好相反。
下面研究一般序列与共轭对称序列和共轭反对称序列之间的关系。
1.将序列写成共轭对称部分和共轭反对称部分相加的形式
一个序列通常可用共轭对称序列与共轭反对称序列之和表示,即
![](https://epubservercos.yuewen.com/688291/21570843601309506/epubprivate/OEBPS/Images/53_09.jpg?sign=1738905539-MKd6kP3f9nf1Jzfr1liBA4q5BhWSZjBc-0-8c2ce83a5a0b9eaf7a4966ab0699a96e)
将式(3-20)中的n用-n代替,并取共轭,得到
![](https://epubservercos.yuewen.com/688291/21570843601309506/epubprivate/OEBPS/Images/54_01.jpg?sign=1738905539-smzUBwah5IfcDc4PKsnKXYFqpjOkiXjD-0-e6871e1f60067b1ea876d5f3321838ae)
对照式(3-20)与式(3-21),有
![](https://epubservercos.yuewen.com/688291/21570843601309506/epubprivate/OEBPS/Images/54_02.jpg?sign=1738905539-DYKzOzl9pmgZjO1xgDYvxufHZfWRQytH-0-8284aca749b4683bedf977724788c3c0)
将式(3-22)和式(3-23)分别进行傅里叶变换得
![](https://epubservercos.yuewen.com/688291/21570843601309506/epubprivate/OEBPS/Images/54_03.jpg?sign=1738905539-uLjxDZ1Yaz7MUJFeqM7xEvUcp92Zr2Ld-0-c9a854979b49a8df3f1eace86c50ba62)
式中,X(ejω)为序列x(n)的傅里叶变换;XR(ejω)和XI(ejω)分别为X(ejω)的实部和虚部。
式(3-24)和式(3-25)说明,如果一个序列写成共轭对称和反对称部分相加的形式,则共轭对称部分的傅里叶变换为原来序列傅里叶变换的实部,共轭反对称部分的傅里叶变换为原来序列傅里叶变换的虚部乘以j。
2.将序列写成实部和虚部相加的形式
如果将序列x(n)写成x(n)=xr(n)+jxi(n)的形式,实部和虚部的傅里叶变换分别为
![](https://epubservercos.yuewen.com/688291/21570843601309506/epubprivate/OEBPS/Images/54_04.jpg?sign=1738905539-QQ7AO3b9njS64TnPchY7uPuFopAUMnHo-0-1714f50dd4f0adea3818d9f61a1e4bc9)
可以证明式(3-26)具有共轭对称的性质,式(3-27)具有共轭反对称的性质,参照时域的共轭对称性,定义
![](https://epubservercos.yuewen.com/688291/21570843601309506/epubprivate/OEBPS/Images/54_05.jpg?sign=1738905539-AYqdfcKYUrf6huJ19Z68210PdfrQxHrL-0-4363cac18465c63152c5da2ead23f987)
式(3-28)和式(3-29)说明,如果一个序列写成实部和虚部相加的形式,则其实部的傅里叶变换Xe(ejω)具有共轭对称的性质,虚部与j相乘的傅里叶变换Xo(ejω)具有共轭反对称的性质。
例3-7 设x(n)=sin(πn/2),-5≤n≤10,用MATLAB程序验证该实序列的对称性质。
解:MATLAB参考脚本如下:
![](https://epubservercos.yuewen.com/688291/21570843601309506/epubprivate/OEBPS/Images/54_06.jpg?sign=1738905539-zFoVn7cBq7Wt2W7pSprib6P40Diu1nVG-0-fa4715981c2bac207772a4362b861fcd)
![](https://epubservercos.yuewen.com/688291/21570843601309506/epubprivate/OEBPS/Images/55_01.jpg?sign=1738905539-lp2aHNfOl5d6rACm1Ake9qSiGSNGpqjM-0-c8dee0456c16cd3389407078ee1470fa)
运行结果如图3-2所示。
由图3-2可看出,如果将序列x(n)写成共轭对称部分和反对称部分相加,其共轭对称部分的傅里叶变换(见图3-2c)等于x(n)的傅里叶变换的实部(见图3-2a),用Re(X)表示。其共轭反对称部分的傅里叶变换(见图3-2d)等于x(n)的傅里叶变换的虚部(见图3-2b),用Im(X)表示。
![](https://epubservercos.yuewen.com/688291/21570843601309506/epubprivate/OEBPS/Images/56_01.jpg?sign=1738905539-yyBxG00Kll9yJFUUBdgyvtwcwMUkqPgk-0-f9caa2be1a100f605ba4378983f06051)
图3-2 例3-7的运行结果