系列文章目录
【学习笔记】奥本海姆第二版《信号与系统》第一章:信号与系统
【学习笔记】奥本海姆第二版《信号与系统》第二章:线性时不变系统
【学习笔记】奥本海姆第二版《信号与系统》第三章:周期信号的傅里叶级数表示
文章目录
一、非周期信号的表示:连续时间傅里叶变换
1.1 非周期信号傅里叶变换的导出
在建立非周期信号的傅里叶变换时,可以把非周期信号当成一个周期信号在周期任意大时的极限来看待,并且研究这个周期信号傅里叶级数表示式的极限特性。现在考虑一个信号 x ( t ) x(t) x(t),它具有有限持续期,即对某个 T 1 T_1 T1,当 ∣ t ∣ > T 1 |t|>T_1 ∣t∣>T1时, x ( t ) = 0 x(t)=0 x(t)=0,如下图所示:
从这个非周期信号出发,可以构成一个周期信号 x ˜ ( t ) \~{x}(t) x˜(t),使 x ( t ) x(t) x(t)就是 x ˜ ( t ) \~{x}(t) x˜(t)的一个周期,如图下图所示:
当把 T T T选得比较大时, x ˜ ( t ) \~{x}(t) x˜(t)就在一个更长的时段上与 x ( t ) x(t) x(t)相一致,并且随着 T → ∞ T→∞ T→∞,对任意有限时间 t t t值而言, x ˜ ( t ) \~{x}(t) x˜(t)就等于 x ( t ) x(t) x(t)。再由第三章的连续时间周期函数的傅里叶级数公式可得:
x ˜ ( t ) = ∑ k = − ∞ + ∞ a k e j k ω 0 t \~{x}(t)=\sum_{k=-∞} ^{+∞}a_ke^{jkω_0t} x˜(t)=k=−∞∑+∞akejkω0t
a k = 1 T ∫ − T / 2 T / 2 x ˜ ( t ) e − j k ω 0 t d t a_k=\frac{1}{T}∫_{-T/2}^{T/2}\~{x}(t)e^{-jkω_0t}dt ak=T1∫−T/2T/2x˜(t)e−jkω0tdt
其中 ω 0 = 2 π T ω_0=\frac{2π}{T} ω0=T2π,由于在 ∣ t ∣ < T 2 |t|<\frac{T}{2} ∣t∣<2T时, x ˜ ( t ) = x ( t ) \~{x}(t)=x(t) x˜(t)=x(t),其它情况下 x ( t ) = 0 x(t)=0 x(t)=0,所以上式可以写成:
a k = 1 T ∫ − T / 2 T / 2 x ( t ) e − j k ω 0 t d t = 1 T ∫ − ∞ ∞ x ( t ) e − j k ω 0 t d t a_k=\frac{1}{T}∫_{-T/2}^{T/2}x(t)e^{-jkω_0t}dt=\frac{1}{T}∫_{-∞}^{∞}x(t)e^{-jkω_0t}dt ak=T1∫−T/2T/2x(t)e−jkω0tdt=T1∫−∞∞x(t)e−jkω0tdt
因此,定义 X ( j ω ) X(jω) X(jω)为:
X ( j ω ) = ∫ − ∞ + ∞ x ( t ) e − j ω t d t X(jω)=∫_{-∞}^{+∞}x(t)e^{-jωt}dt X(jω)=∫−∞+∞x(t)e−jωtdt
此时系数 a k a_k ak可以写成
a k = 1 T X ( j k ω 0 ) a_k=\frac{1}{T}X(jkω_0) ak=T1X(jkω0)
由 x ˜ ( t ) \~x(t) x˜(t)的傅里叶级数公式和上式结合可得:
x ˜ ( t ) = ∑ k = − ∞ + ∞ 1 T X ( j k ω 0 ) e j k ω 0 t \~x(t)=\sum_{k=-∞} ^{+∞}\frac{1}{T}X(jkω_0)e^{jkω_0t} x˜(t)=k=−∞∑+∞T1X(jkω0)ejkω0t
又因为 2 π T = ω 0 \frac{2π}{T}=ω_0 T2π=ω0,所以上式又可以表示为:
x ˜ ( t ) = 1 2 π ∑ k = − ∞ + ∞ X ( j k ω 0 ) e j k ω 0 t ω 0 \~x(t)=\frac{1}{2π}\sum_{k=-∞} ^{+∞}X(jkω_0)e^{jkω_0t}ω_0 x˜(t)=2π1k=−∞∑+∞X(jkω0)ejkω0tω0
随着 T → ∞ T→∞ T→∞, x ˜ ( t ) \~x(t) x˜(t)趋近与 x ( t ) x(t) x(t),上式的极限形式就变成了 x ( t ) x(t) x(t)的表达式。最终:
X ( j ω ) = ∫ − ∞ + ∞ x ( t ) e − j ω t d t X(jω)=∫_{-∞}^{+∞}x(t)e^{-jωt}dt X(jω)=∫−∞+∞x(t)e−jωtdt
x ( t ) = 1 2 π ∫ − ∞ + ∞ X ( j ω ) e j ω t d ω x(t)=\frac{1}{2π}∫_{-∞}^{+∞}X(jω)e^{jωt}dω x(t)=2π1∫−∞+∞X(jω)ejωtdω
称为傅里叶变换对,函数 X ( j ω ) X(jω) X(jω)称为 x ( t ) x(t) x(t)的傅里叶变换。一个非周期信号的傅里叶变换 X ( j ω ) X(jω) X(jω)通常称为 x ( t ) x(t) x(t)的频谱。
1.2 傅里叶变换的收敛
需要满足狄里赫利条件:
- 条件一: x ( t ) x(t) x(t)绝对可积,即:
∫ − ∞ + ∞ ∣ x ( t ) ∣ d t < ∞ ∫_{-∞}^{+∞}|x(t)|dt<∞ ∫−∞+∞∣x(t)∣dt<∞
- 条件二:在任意有限区间内, x ( t ) x(t) x(t)只有有限个最大值和最小值
- 条件三:在任意有限区间内, x ( t ) x(t) x(t)有有限个不连续点,并且在这些不连续点上都必须是有限值。因此,本身是连续的或者只有有限个不连续点的绝对可积信号都存在傅里叶变换。
1.3 几个常见的典型傅里叶变换对
- 单边指数信号:
x ( t ) = e − a t u ( t ) , a > 0 x(t)=e^{-at}u(t) ,a>0 x(t)=e−atu(t),a>0
其傅里叶变换 X ( j ω ) X(jω) X(jω):
X ( j ω ) = ∫ 0 ∞ e − a t e − j ω t d t = − 1 a + j ω e − ( a + j ω ) t ∣ 0 ∞ = 1 a + j ω , a > 0 X(jω)=∫_{0}^{∞}e^{-at}e^{-jωt}dt=-\frac{1}{a+jω}e^{-(a+jω)t}|_{0}^{∞}=\frac{1}{a+jω},a>0 X(jω)=∫0∞e−ate−jωtdt=−a+jω1e−(a+jω)t∣0∞=a+jω1,a>0
- 单位冲击函数:
x ( t ) = δ ( t ) x(t)=δ(t) x(t)=δ(t)
其傅里叶变换 X ( j ω ) X(jω) X(jω):
X ( j ω ) = ∫ − ∞ + ∞ δ ( t ) e − j ω t d t = 1 X(jω)=∫_{-∞}^{+∞}δ(t)e^{-jωt}dt=1 X(jω)=∫−∞+∞δ(t)e−jωtdt=1
- 双边指数信号:
x ( t ) = e − a ∣ t ∣ , a > 0 x(t)=e^{-a|t|} ,a>0 x(t)=e−a∣t∣,a>0
其傅里叶变换 X ( j ω ) X(jω) X(jω):
X ( j ω ) = ∫ − ∞ + ∞ e − a ∣ t ∣ e − j ω t d t = ∫ − ∞ 0 e a t e − j ω t d t + ∫ 0 + ∞ e − a t e − j ω t d t X(jω)=∫_{-∞}^{+∞}e^{-a|t|}e^{-jωt}dt=∫_{-∞}^{0}e^{at}e^{-jωt}dt+∫_{0}^{+∞}e^{-at}e^{-jωt}dt X(jω)=∫−∞+∞e−a∣t∣e−jωtdt=∫−∞0eate−jωtdt+∫0+∞e−ate−jωtdt
= ∫ − ∞ 0 e ( a − j ω ) t d t + ∫ 0 + ∞ e − ( a + j ω ) t d t =∫_{-∞}^{0}e^{(a-jω)t}dt+∫_{0}^{+∞}e^{-(a+jω)t}dt =∫−∞0e(a−jω)tdt+∫0+∞e−(a+jω)tdt
= 1 a − j ω + = 1 a + j ω =\frac{1}{a-jω}+=\frac{1}{a+jω} =a−jω1+=a+jω1
= 2 a a 2 + ω 2 =\frac{2a}{a^2+ω^2} =a2+ω22a
- 抽样函数:
x ( t ) = s i n ( ω c t ) π t x(t)=\frac{sin(ω_ct)}{πt} x(t)=πtsin(ωct)
其傅里叶变换 X ( j ω ) X(jω) X(jω):
X ( j ω ) = { 0 ,其它 1 , − ω c < ω < w c X(jω)=\lbrace_{0,其它}^{1,-ω_c<ω<w_c} X(jω)={0,其它1,−ωc<ω<wc
- 冲击偶函数:
x ( t ) = δ ′ ( t ) x(t)=δ^{'}(t) x(t)=δ′(t)
根据冲击偶函数的取样性质:
∫ − ∞ + ∞ x ( t ) δ ′ ( t ) d t = − x ′ ( 0 ) ∫_{-∞}^{+∞}x(t)δ^{'}(t)dt=-x^{'}(0) ∫−∞+∞x(t)δ′(t)dt=−x′(0)
所以其傅里叶变换 X ( j ω ) X(jω) X(jω):
X ( j ω ) = ∫ − ∞ + ∞ δ ′ ( t ) e − j ω t d t X(jω)=∫_{-∞}^{+∞}δ^{'}(t)e^{-jωt}dt X(jω)=∫−∞+∞δ′(t)e−jωtdt
= − d ( e − j ω t ) d t ∣ t = 0 = - \frac{d(e^{-jωt})}{dt}|_{t=0} =−dtd(e−jωt)∣t=0
= j ω e − j ω t ∣ t = 0 =jωe^{-jωt}|_{t=0} =jωe−jωt∣t=0
= j ω =jω =jω
- 矩形窗函数:
x ( t ) = { 0 ,其它 E , − τ / 2 < ω < τ / 2 x(t)=\lbrace_{0,其它}^{E,-τ/2<ω<τ/2} x(t)={0,其它E,−τ/2<ω<τ/2
其傅里叶变换 X ( j ω ) X(jω) X(jω):
X ( j ω ) = ∫ − τ 2 τ 2 E e − j ω t d t X(jω)=∫_{-\frac{τ}{2}}^{\frac{τ}{2}}Ee^{-jωt}dt X(jω)=∫−2τ2τEe−jωtdt
= E 1 j ω e − j ω t ∣ τ 2 − τ 2 = E \frac{1}{jω}e^{-jωt}|_{\frac{τ}{2}}^{-\frac{τ}{2}} =Ejω1e−jωt∣2τ−2τ
= E 1 j ω 2 j s i n ( ω τ 2 ) =E\frac{1}{jω}2jsin(\frac{ωτ}{2}) =Ejω12jsin(2ωτ)
= 2 E ω s i n ( ω τ 2 ) =\frac{2E}{ω}sin(\frac{ωτ}{2}) =ω2Esin(2ωτ)
= E τ S a ( τ ω 2 ) =EτSa(\frac{τω}{2}) =EτSa(2τω)
- 常数:
x ( t ) = 1 x(t)=1 x(t)=1
其傅里叶变换 X ( j ω ) X(jω) X(jω):
X ( j ω ) = 2 π δ ( ω ) X(jω)=2πδ(ω) X(jω)=2πδ(ω)
- 阶跃函数:
x ( t ) = u ( t ) x(t)=u(t) x(t)=u(t)
其傅里叶变换 X ( j ω ) X(jω) X(jω):
X ( j ω ) = 1 j ω + π δ ( ω ) X(jω)=\frac{1}{jω}+πδ(ω) X(jω)=jω1+πδ(ω)
二、连续时间傅里叶变换性质
为了方便,将 X ( j ω ) X(jω) X(jω)用 Γ { x ( t ) } \Gamma\{x(t)\} Γ{x(t)}表示,将 x ( t ) x(t) x(t)用 Γ − 1 { X ( j ω ) } \Gamma^{-1}\{X(jω)\} Γ−1{X(jω)}表示,符号如下:
x ( t ) ⟷ Γ X ( j ω ) x(t)\stackrel{\Gamma}{\longleftrightarrow}X(jω) x(t)⟷ΓX(jω)
2.1 线性性质
若:
x ( t ) ⟷ Γ X ( j ω ) x(t)\stackrel{\Gamma}{\longleftrightarrow}X(jω) x(t)⟷ΓX(jω)
且:
y ( t ) ⟷ Γ Y ( j ω ) y(t)\stackrel{\Gamma}{\longleftrightarrow}Y(jω) y(t)⟷ΓY(jω)
则:
a x ( t ) + b y ( t ) ⟷ Γ a X ( j ω ) + b Y ( j ω ) ax(t)+by(t)\stackrel{\Gamma}{\longleftrightarrow}aX(jω)+bY(jω) ax(t)+by(t)⟷ΓaX(jω)+bY(jω)
2.2 时移性质
若:
x ( t ) ⟷ Γ X ( j ω ) x(t)\stackrel{\Gamma}{\longleftrightarrow}X(jω) x(t)⟷ΓX(jω)
则:
x ( t − t 0 ) ⟷ Γ e − j ω t 0 X ( j ω ) x(t-t_0)\stackrel{\Gamma}{\longleftrightarrow}e^{-jωt_0}X(jω) x(t−t0)⟷Γe−jωt0X(jω)
2.3 共轭对称性
实偶 ⟷ Γ 实偶 实偶\stackrel{\Gamma}{\longleftrightarrow}实偶 实偶⟷Γ实偶
实奇 ⟷ Γ 实奇 实奇\stackrel{\Gamma}{\longleftrightarrow}实奇 实奇⟷Γ实奇
实函数的傅里叶变换实部为偶数,虚部为奇函数。实函数的幅频特性为偶函数,相频特性为奇函数。
2.4 频移性质
若:
x ( t ) ⟷ Γ X ( j ω ) x(t)\stackrel{\Gamma}{\longleftrightarrow}X(jω) x(t)⟷ΓX(jω)
则:
x ( t ) e j ω 0 t ⟷ Γ X ( j ( ω − ω 0 ) ) x(t)e^{jω_0t}\stackrel{\Gamma}{\longleftrightarrow}X(j(ω-ω_0)) x(t)ejω0t⟷ΓX(j(ω−ω0))
2.5 微分性质
若:
x ( t ) ⟷ Γ X ( j ω ) x(t)\stackrel{\Gamma}{\longleftrightarrow}X(jω) x(t)⟷ΓX(jω)
则:
x ′ ( t ) ⟷ Γ j ω X ( j ω ) x^{'}(t)\stackrel{\Gamma}{\longleftrightarrow}jωX(jω) x′(t)⟷ΓjωX(jω)
2.6 积分性质
若:
x ( t ) ⟷ Γ X ( j ω ) x(t)\stackrel{\Gamma}{\longleftrightarrow}X(jω) x(t)⟷ΓX(jω)
则:
∫ − ∞ t x ( u ) d u ⟷ Γ X ( j ω ) j ω + π X ( 0 ) δ ( ω ) ∫_{-∞}^tx(u)du\stackrel{\Gamma}{\longleftrightarrow}\frac{X(jω)}{jω}+πX(0)δ(ω) ∫−∞tx(u)du⟷ΓjωX(jω)+πX(0)δ(ω)
2.7 时间与频率尺度变换
若:
x ( t ) ⟷ Γ X ( j ω ) x(t)\stackrel{\Gamma}{\longleftrightarrow}X(jω) x(t)⟷ΓX(jω)
则:
x ( a t ) ⟷ Γ 1 ∣ a ∣ X ( j ω a ) x(at)\stackrel{\Gamma}{\longleftrightarrow}\frac{1}{|a|}X(\frac{jω}{a}) x(at)⟷Γ∣a∣1X(ajω)
2.8 帕斯瓦尔定理
若 x ( t ) x(t) x(t) 和 X ( j ω ) X(jω) X(jω)是一对傅里叶变换,则:
∫ − ∞ + ∞ ∣ x ( t ) ∣ 2 d t = 1 2 π ∫ − ∞ + ∞ ∣ X ( j ω ) ∣ 2 d ω ∫_{-∞}^{+∞}|x(t)|^2dt=\frac{1}{2π}∫_{-∞}^{+∞}|X(jω)|^2dω ∫−∞+∞∣x(t)∣2dt=2π1∫−∞+∞∣X(jω)∣2dω
左边是信号 x ( t ) x(t) x(t)的总能量。帕斯瓦尔定理指出,这个总能量既可以按每单位时间内的能量 ∣ x ( t ) ∣ 2 |x(t)|^2 ∣x(t)∣2在整个时间内积分计算出来,也可以按每单位频率内的能量 ∣ X ( j ω ) ∣ 2 |X(jω)|^2 ∣X(jω)∣2在整个频率范围内积分而得到。因此, ∣ X ( j ω ) ∣ 2 |X(jω)|^2 ∣X(jω)∣2常称为信号 x ( t ) x(t) x(t)的能谱密度。
2.9 卷积性质
若:
x ( t ) ⟷ Γ X ( j ω ) x(t)\stackrel{\Gamma}{\longleftrightarrow}X(jω) x(t)⟷ΓX(jω)
且:
y ( t ) ⟷ Γ Y ( j ω ) y(t)\stackrel{\Gamma}{\longleftrightarrow}Y(jω) y(t)⟷ΓY(jω)
则:
x ( t ) ∗ y ( t ) ⟷ Γ X ( j ω ) Y ( j ω ) x(t)*y(t)\stackrel{\Gamma}{\longleftrightarrow}X(jω)Y(jω) x(t)∗y(t)⟷ΓX(jω)Y(jω)
结论:时域卷积等于频域相乘
举例:在上一章我们知道了理想低通滤波的频率响应如下:
H ( j ω ) = { 0 , ∣ ω ∣ > ω c 1 , ∣ ω ∣ < ω c H(jω)=\{\ _{0,|ω|>ω_c} ^{1,|ω|<ω_c} H(jω)={ 0,∣ω∣>ωc1,∣ω∣<ωc
所以,根据傅里叶逆变换可得单位冲击响应 h ( t ) h(t) h(t)为:
h ( t ) = s i n ( ω c t ) π t h(t)=\frac{sin(ω_ct)}{πt} h(t)=πtsin(ωct)
如图所示:
这个例子已经能够开始看到在滤波器设计中所出现的一些问题,滤波器设计中涉及到时域和频域两方面的要求。尽管理想低通滤波器确实有非常完美的频率选择性,但是它的单位冲激响应的某些特性却可能是我们不希望的。首先注意到 h ( t ) h(t) h(t)在 t < 0 t<0 t<0时不是零,其结果就是理想低通滤波器不是因果的,因此在要求因果系统的应用中,就无法采用理想低通滤波器。
2.10 调制(相乘)性质
若:
x 1 ( t ) ⟷ Γ X 1 ( j ω ) x_1(t)\stackrel{\Gamma}{\longleftrightarrow}X_1(jω) x1(t)⟷ΓX1(jω)
且:
x 2 ( t ) ⟷ Γ X 2 ( j ω ) x_2(t)\stackrel{\Gamma}{\longleftrightarrow}X_2(jω) x2(t)⟷ΓX2(jω)
则:
x 1 ( t ) x 2 ( t ) ⟷ Γ 1 2 π X 1 ( j ω ) ∗ X 2 ( j ω ) x_1(t)x_2(t)\stackrel{\Gamma}{\longleftrightarrow}\frac{1}{2π}X_1(jω)*X_2(jω) x1(t)x2(t)⟷Γ2π1X1(jω)∗X2(jω)
三、傅里叶变换性质和基本傅里叶变换对列表
四、小结
这一章建立了连续时间信号的傅里叶变换表示,并研究了许多很有用的性质。特别是在把一个非周期信号看成周期变得任意大时一个周期信号的极限之后。另外,周期信号本身也可以用傅里叶变换来表示,这个傅里叶变换由发生在该周期信号各谐波频率上的冲激串所组成,并且每个冲激串的面积正比于各傅里叶级数系数。傅里叶变换具有一系列重要性质,这些性质表达了不同的信号特性是如何反映到它们的变换中去的。