在之前的 B l o g Blog Blog 里面,我们都已经进行过了推导。我们这里再一次简洁地推一下下:
对于下面这样的周期矩形信号,我们可以有下面的结论:
a k = 1 T ∫ − T 1 + T 1 x ( t ) ~ e − j k ω 0 t d t (1) a_k = \frac{1}{T}\int_{-T_1}^{+T_1}\widetilde{x(t)}e^{-jkω_0t}dt\tag{1} ak=T1∫−T1+T1x(t) e−jkω0tdt(1)
下面我们考虑这样一个非周期信号:
我们发现,其实在区间 − T 1 -T_1 −T1 ~ T 1 T_1 T1 内,有: x ( t ) = x ( t ) ~ x(t) = \widetilde{x(t)} x(t)=x(t) 。那么我们把这个式子带入(1)得到: a k = 1 T ∫ − T 1 + T 1 x ( t ) e − j k ω 0 t d t T a k = ∫ − T 1 + T 1 x ( t ) e − j k ω 0 t d t \begin{aligned} &a_k = \frac{1}{T}\int_{-T_1}^{+T_1}x(t)e^{-jkω_0t}dt\\ &Ta_k = \int_{-T_1}^{+T_1}x(t)e^{-jkω_0t}dt\\ \end{aligned} ak=T1∫−T1+T1x(t)e−jkω0tdtTak=∫−T1+T1x(t)e−jkω0tdt
我们发现,现在的式子就变成针对 x ( t ) x(t) x(t) 的了,那么对于 x ( t ) x(t) x(t) ,他不是只在 区间 − T 1 -T_1 −T1 ~ T 1 T_1 T1 内才有值嘛,那么其实在区间 − T 1 -T_1 −T1 ~ T 1 T_1 T1 内积分就可以等价于在 − ∞ -∞ −∞ ~ + ∞ +∞ +∞ 内积分了: T a k = ∫ − ∞ + ∞ x ( t ) e − j k ω 0 t d t Ta_k = \int_{-∞}^{+∞}x(t)e^{-jkω_0t}dt Tak=∫−∞+∞x(t)e−jkω0tdt
当 T T T 趋近于无穷时,谱线的间隔趋于0,那么 k ω 0 kω_0 kω0(本来代表的就是一根一根的谱线的横坐标,因为现在谱线间隔没有了,所以这个 k ω 0 kω_0 kω0 就可以用一个连续的 ω ω ω 来表示了: X ( j ω ) = ∫ − ∞ + ∞ x ( t ) e − j ω t d t X(jω) = \int_{-∞}^{+∞}x(t)e^{-jωt}dt X(jω)=∫−∞+∞x(t)e−jωtdt
傅里叶逆变换的推导思路也很清晰:
还记得我们在上一篇 B l o g Blog Blog 里面说过的吗?连续时间周期矩形信号的频谱可以看成是对连续时间非周期矩形信号频谱的采样。 所以我们有下面的式子: a k = 1 T X ( j k ω 0 ) a_k = \frac{1}{T}X(jkω_0) ak=T1X(jkω0)
下面我们依然考虑一个连续时间周期矩形信号的傅里叶级数表示: x ( t ) ~ = ∑ k = − ∞ + ∞ a k e j k ω 0 t = ∑ k = − ∞ + ∞ 1 T X ( j k ω 0 ) e j k ω 0 t = 1 2 π ∑ k = − ∞ + ∞ X ( j k ω 0 ) e j k ω 0 t ω 0 \begin{aligned} \widetilde{x(t)} &= \sum_{k=-∞}^{+∞}a_ke^{jkω_0t}\\ &=\sum_{k=-∞}^{+∞}\frac{1}{T}X(jkω_0)e^{jkω_0t}\\ &=\frac{1}{2π}\sum_{k=-∞}^{+∞}X(jkω_0)e^{jkω_0t}ω_0 \end{aligned} x(t) =k=−∞∑+∞akejkω0t=k=−∞∑+∞T1X(jkω0)ejkω0t=2π1k=−∞∑+∞X(jkω0)ejkω0tω0
当 T T T 趋于无穷时,我们用求和去极限的方式表示(就引入了积分):
x ( t ) = lim T → ∞ x ( t ) ~ = lim T → ∞ 1 2 π ∑ k = − ∞ + ∞ X ( j k ω 0 ) e j k ω 0 t ω 0 = 1 2 π lim T → ∞ ∑ k = − ∞ + ∞ X ( j k ω 0 ) e j k ω 0 t ω 0 = 1 2 π ∫ − ∞ + ∞ X ( j ω ) e j ω t d ω \begin{aligned} x(t) = \lim_{T\to ∞} \widetilde{x(t)} &=\lim_{T\to ∞}\frac{1}{2π}\sum_{k=-∞}^{+∞}X(jkω_0)e^{jkω_0t}ω_0\\ &=\frac{1}{2π}\lim_{T\to ∞}\sum_{k=-∞}^{+∞}X(jkω_0)e^{jkω_0t}ω_0\\ &=\frac{1}{2π}\int_{-∞}^{+∞}X(jω)e^{jωt}dω\\ \end{aligned} x(t)=T→∞limx(t) =T→∞lim2π1k=−∞∑+∞X(jkω0)ejkω0tω0=2π1T→∞limk=−∞∑+∞X(jkω0)ejkω0tω0=2π1∫−∞+∞X(jω)ejωtdω
因此,我们就推出了傅里叶逆变换公式。
下面,我们把傅里叶正变换和逆变换写到一起方便记忆: X ( j ω ) = ∫ − ∞ + ∞ x ( t ) e − j ω t d t X(jω) = \int_{-∞}^{+∞}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π}\int_{-∞}^{+∞}X(jω)e^{jωt}dω x(t)=2π1∫−∞+∞X(jω)ejωtdω
有两个条件:满足其中之一即可:
【1】信号的能量有限: ∫ − ∞ + ∞ ∣ x ( t ) ∣ 2 d t < ∞ \int_{-∞}^{+∞}|x(t)|^2dt < ∞ ∫−∞+∞∣x(t)∣2dt<∞
【2】信号绝对可积: ∫ − ∞ + ∞ ∣ x ( t ) ∣ d t < ∞ \int_{-∞}^{+∞}|x(t)|dt < ∞ ∫−∞+∞∣x(t)∣dt<∞
连续时间周期信号的频谱是其对应的非周期信号频谱的样本。(注:这个“对应的” 应该怎么解释呢?因为任何一个周期信号都可以是其一个周期的信号做周期延拓得来的)
非周期信号的频谱就是其对应的周期信号频谱的包络
在记录之前,我们打算说明一下广义采样函数的定义。我们之前的 s i n c sinc sinc 函数是狭义的采样函数,那么为了和老师授课所用的符号一致,我们以后在某些场合将采用广义的采样函数 S a Sa Sa 函数,其函数表达式为: S a ( x ) = s i n ( x ) x Sa(x) = \frac{sin(x)}{x} Sa(x)=xsin(x)
当 x = k π ( k = ± 1 , ± 2 , ⋯ ) x = kπ(k = ±1, ±2, \cdots) x=kπ(k=±1,±2,⋯)时, S a ( x ) = 0 Sa(x) = 0 Sa(x)=0
对于狭义的 s i n c ( x ) sinc(x) sinc(x) 函数,当 x = ± k x = ±k x=±k( k k k 为正整数时), s i n c ( x ) = 0 sinc(x) = 0 sinc(x)=0
对于非周期矩形信号的频谱,大家可以按照我以前那个 B l o g Blog Blog 里面的方法记忆(更容易记忆一些):
若这个非周期矩形信号的脉冲宽度是 τ τ τ,那么其频谱函数就是: τ s i n c ( τ f ) τsinc(τf) τsinc(τf)
我们看看下面的例子:
你看,脉冲宽度是 2 T 1 2T_1 2T1,那么频谱直接出来啦,就是: 2 T 1 s i n c ( 2 T 1 f ) 2T_1sinc(2T_1f) 2T1sinc(2T1f),你想把它写成 ω ω ω 的表达式也很方便,就是: 2 T 1 s i n c ( 2 T 1 ω 2 π ) = 2 T 1 s i n c ( T 1 ω π ) 2T_1sinc(2T_1\frac{ω}{2π}) = 2T_1sinc(T_1\frac{ω}{π}) 2T1sinc(2T12πω)=2T1sinc(T1πω)
根据狭义采样函数的定义,在: s i n c = 0 ( i f T 1 ω π = ± 1 , ± 2 , ⋯ ) sinc = 0\space(if\space T_1\frac{ω}{π} =±1, ±2, \cdots) sinc=0 (if T1πω=±1,±2,⋯)
我们可以得到频谱的第一个零点的坐标就是: ω = π T 1 ω = \frac{π}{T_1} ω=T1π
如果给出的是一个矩形的频谱,要你求原时域函数呢?如下图:
我们先直接用公式求解: x ( t ) = 1 2 π ∫ − ∞ + ∞ X ( j ω ) e j ω t d ω = 1 2 π ∫ − ω + ω e j ω t d ω = s i n ( ω t ) π t = ω π s i n ( ω t ) ω t = ω π S a ( ω t ) \begin{aligned} x(t) &= \frac{1}{2π}\int_{-∞}^{+∞}X(jω)e^{jωt}dω\\ &=\frac{1}{2π}\int_{-ω}^{+ω}e^{jωt}dω\\ &=\frac{sin(ωt)}{πt}\\ &=\frac{ω}{π}\frac{sin(ωt)}{ωt} = \frac{ω}{π}Sa(ωt) \end{aligned} x(t)=2π1∫−∞+∞X(jω)ejωtdω=2π1∫−ω+ωejωtdω=πtsin(ωt)=πωωtsin(ωt)=πωSa(ωt)
对于 Sa(x) 函数而言,我们之前也提到过了,当 x = k π ( k = ± 1 , ± 2 , ⋯ ) x = kπ(k = ±1, ±2, \cdots) x=kπ(k=±1,±2,⋯)时, S a ( x ) = 0 Sa(x) = 0 Sa(x)=0
所以, x ( t ) x(t) x(t) 在正半时间轴的第一个零点就是: x = π ω x = \frac{π}{ω} x=ωπ
这个结论我们直接记忆:单位冲激信号的傅里叶变换是幅度为 1 的频谱:
大家看分析一下:单位冲激信号的频谱居然包括了所有频率范围,并且各个频率成分的幅度都是 1 !相位也相同! 这也是为什么 δ δ δ 信号在这门课如此重要的原因!
如果某一个信号的频谱是单位冲激,那么做一下傅里叶逆变换: x ( t ) = 1 2 π ∫ − ∞ + ∞ δ ( j ω ) e j ω t d ω = 1 2 π ∫ − ∞ + ∞ δ ( j ω ) e 0 d ω = 1 2 π \begin{aligned} x(t) &= \frac{1}{2π}\int_{-∞}^{+∞}δ(jω)e^{jωt}dω\\ &=\frac{1}{2π}\int_{-∞}^{+∞}δ(jω)e^0dω\\ &=\frac{1}{2π} \end{aligned} x(t)=2π1∫−∞+∞δ(jω)ejωtdω=2π1∫−∞+∞δ(jω)e0dω=2π1
【注:上面推导的第二行用了单位冲击信号的采样性】
由此可见,幅度为 1 2 π \frac{1}{2π} 2π1 的直流信号的频谱是单位冲激
x ( t ) = e − a t u ( t ) ( a > 0 ) x(t) = e^{-at}u(t)\space(a > 0) x(t)=e−atu(t) (a>0)
其频谱就是: X ( j ω ) = 1 a + j ω X(jω) = \frac{1}{a + jω} X(jω)=a+jω1