# 公式表
# 傅里叶变换性质
# 基本傅里叶变换对
连续时间非周期信号傅里叶变换的推导基本思想
- 非周期函数可以通过周期性延拓得到周期函数
- 周期函数可以令基波周期趋于无穷大得到连续的非周期函数
对于周期性矩形脉冲信号:
akTak=T2T1Sa(kπT2T1)=2T1Sa(ωT1)
其中:ω=kω0=T2kπ。
两点之间的距离:
Δω=(k+1)ω0−kω0=ω0=T2π
当 T→∞ 时有 ω0=T2π→0 ,即两点之间距离趋于 0 ,函数包络线趋于连续函数。
# 非周期信号
# FT 公式推导
连续时间非周期信号傅里叶变换的推导
周期函数 x~(t) 以 T 为周期,在 [−2T,2T] 内与 x(t) 相等,它的傅里叶级数表示为:
x~(t)=k=−∞∑+∞akejkω0t
当 T→∞ 时,在 (−∞,∞) 内 x~(t)=x(t) ;
ak=T1∫−2T2Tx~(t)e−jkω0tdt=T1∫−2T2Tx(t)e−jkω0tdt=T1∫−∞∞x(t)e−jkω0tdt
定义频谱:(傅里叶正变换)
X(jω)=∫−∞∞x(t)e−jkω0tdt
从而 ak 可以表示为:
ak=T1X(jkω0)
因此,在 [−2T,2T] 内有:
x(t)=x~(t)=k=−∞∑+∞akejkω0t=T1k=−∞∑+∞X(jkω0)ejkω0t=2π1k=−∞∑+∞X(jkω0)ejkω0tω0=2π1k=−∞∑+∞X(jω)ejωtΔω
当 T→∞ 时,Δω=ω0→0,即 Δω→dω 则得到(傅里叶反变换):
x(t)=x~(t)=2π1∫−∞+∞X(jω)ejωtdω
从而得到傅里叶变换对:
这两个公式是 “一对” 的,要么同时可积,要么都不可积
傅里叶正变换,分析公式,(ejωt 的线性组合)
X(jω)=∫−∞∞x(t)e−jωtdt
傅里叶反变换,综合公式,(获得 x(t) 的频谱)
x(t)=2π1∫−∞+∞X(jω)ejωtdω
可以表示为:
x(t)⟷CTFTX(jω)F{x(t)}=X(jω)
周期信号傅里叶级数的系数为其一个周期内信号傅里叶变换的等间隔采样:
{x(t)↔X(jω)x~(t)↔ak⇒ak=T1X(jkω0)
# 收敛性
既然傅里叶变换是从周期信号的傅里叶级数引出的,傅里叶变换的收敛问题就应该和傅里叶级数的收敛相一致。所以也有对应的两组条件;同样地两组条件并不完全等价。它们都是傅里叶变换存在的充分条件。
和周期信号的情况一样,当 x(t) 的傅里叶变换 X(jω) 存在时,其反变换在间断点附近会产生 Gibbs 现象
# 平方可积条件
# Dirichlet 条件
# 常见信号的 FT
# 低通信号 (单侧)
x(t)=e−atu(t), a>0
有正变换:
X(jω)∣X(jω)∣dωd∣X(jω)∣=∫0+∞e−ate−jωtdt=a+jω1=a2+ω21=−arctanaω
# 低通信号 (双侧)
x(t)=e−a∣t∣, a>0
有正变换:
X(jω)∣X(jω)∣=∫−∞0eate−jωtdt+∫0+∞e−ate−jωtdt=a−jω1+a+jω1=a2+ω22a=X(jω)
# 冲激信号
x(t)=δ(t)
有正变换:
X(jω)=∫−∞+∞δ(t)e−jωtdt=1
这表明 δ(t) 中包括了所有的频率成分,且所有频率分量的幅度、相位都相同。因此,系统的单位冲激响应 h(t) 才能完全描述一个 LTI 系统的特性,δ(t) 才在信号与系统分析中具有如此重要的意义。
# 常数信号
* 这个名字我乱起的
x(t)=1
因为:
2π1∫−∞+∞2π δ(ω) ejωtdω=1
所以有反变换:
X(jω)=2πδ(ω)
冲激信号和双侧阶跃信号 (我乱起名) 的对偶关系:
δ(t)⟷F1⟷F2πδ(ω)
# 矩形信号
x(t)={1, ∣t∣<T10, ∣t∣>T1
有正变换:
X(jω)ak=TX(jω)=∫−T1T1e−jωtdt=ω2sinωT1=ωT12T1sinωT1=2T1Sa(ωT1)=2T1sinc(πωT1)=T2T1Sa(kω0T1)=T2T1Sa(T2T1)
不同脉冲宽度对频谱的影响:
可见,信号宽度在时域和频域之间有一种相反的关系。
# 采样信号
X(jω)={1, ∣ω∣<W0, ∣ω∣>W
有反变换:
x(t)=2π1∫−W+We−jωtdω=πtsinWt=πWSa(Wt)=πWSa(πWt)
可见 x(t) 是采样信号
矩形信号和采样信号的对偶关系:(上矩形信号、下采样信号)
# 周期信号
Tak 无论在 T 有限还是无限的时候都可以是有限的,所以我们可以将 T 有限时的周期信号和 T 无限时的非周期信号用傅里叶变换统一起来
但是 ak 却不是,所以我们无法将周期信号和非周期信号都用傅里叶级数表示。
# FT 公式推导
考虑 X(jω)=2πδ(ω−ω0) 所对应的信号:
x(t)=2π1∫−∞+∞X(jω) ejωtdω=∫−∞+∞δ(ω−ω0) ejωtdω=ejω0t
于是对于表示为傅里叶级数的周期信号 x(t):
x(t)=k=−∞∑+∞ak ejkω0t
就有频谱即正变换:
X(jω)=2πk=−∞∑+∞ak δ(ω−ω0)
这表明周期信号的傅里叶正变换由一系列位于谐波频率处的冲激组成,其中 kω0 处的冲激强度为对应傅里叶级数系数 ak 的 2π 倍。
# 常见信号的 FT
# 正弦信号
x(t)=sin(ω0t)=2j1(ejω0t−e−jω0t)
有正变换:
X(jω)=jπ[δ(ω−ω0)−δ(ω+ω0)]
# 余弦信号
x(t)=cos(ω0t)=21(ejω0t+e−jω0t)
有正变换:
X(jω)=π[δ(ω−ω0)+δ(ω+ω0)]
# 均匀冲激串
x(t)=n=−∞∑+∞δ(t−nT)
其中 T=ω02π;
其傅里叶系数:
ak=T1∫−2T2Tδ(t) e−jT2πktdt=T1∫−2T2Tδ(t)dt=T1
可得傅里叶正变换:
X(jω)=T2πk=−∞∑+∞δ(ω−T2πk)
# 周期性矩形脉冲
其傅里叶系数:
ak=T2T1Sa(πkT2T1)
有傅里叶正变换:
X(jω)=T4πT1k=−∞∑+∞Sa(πkT2T1) δ(ω−T2πk)=2T1ω0k=−∞∑+∞Sa(ω0T1⋅k) δ(ω−ω0⋅k)
其中 ω0=T2π 为谱线间隔。
# 特殊信号 FT
# 符号函数
可以将符号函数定义为:
Sgn(t)=a→0lim[e−atu(t)−eatu(−t)]
其傅里叶变换为:
F{Sgn(t)}=a→0lim[a+jω1−a−jω1]=a→0lima2+(jω)2−j2ω=jω2
# 单位阶跃函数
将 u(t) 分解为偶部和奇部有:
u(t)=ue(t)+uo(t)ue(t)=21uo(t)=21Sgn(t)
而 1⟷CTFT2πδ(ω) ,Sgn(t)⟷CTFTjω2 ,因此有:
u(t)⟷Fjω1+πδ(ω)
# 性质
# 线性
x(t)⟷FX(jω) , y(t)⟷FY(jω)⇓ax(t)+by(t)⟷FaX(jω)+bY(jω)
# 时移特性
信号的时移会导致其频谱产生一个线性相移
x(t)⟷FX(jω)⇓x(t−t0)⟷FX(jω)e−jωt0
证明
F{x(t−t0)}=∫−∞+∞x(t−t0) e−jωtdt=∫−∞+∞x(t′) e−jω(t′+t0)dt′=e−jωt0∫−∞+∞x(t′) e−jωt′dt′=e−jωt0X(jω)
# 共轭对称性
x(t)⟷FX(jω)⇓x∗(t)⟷FX∗(−jω)
特别地,当 x(t) 是实信号 (实函数),则有频谱: X(jω)=X∗(−jω) 为共轭偶函数(实幅偶,虚相奇)
证明
由 x(t) 的傅里叶正变换两侧同时取共轭有:
X∗(jω)=∫−∞∞[x(t)e−jkω0t]∗dt=∫−∞∞x(t)∗ejkω0tdt
所以:
X∗(−jω)=∫−∞∞x(t)∗e−jkω0tdt
即得证:
x∗(t)⟷FX∗(−jω)
当 x(t) 是实信号时有 x(t)=x∗(t) ,于是有:
X(jω)=X∗(−jω)=Re[X(−jω)]+j Im[X(−jω)]
共轭偶函数:
- 实部是偶函数:Re[X(jω)]=Re[X(−jω)];
- 虚部是奇函数:Im[X(jω)]=−Im[X(−jω)];
- 幅度是偶函数:∣X(jω)∣=∣X(−jω)∣;
- 相位是奇函数:ArgX(jω)=−ArgX(−jω);
# 奇偶对称性
- 偶信号的傅里叶 (正) 变换还是偶函数
- 实偶信号的傅里叶 (正) 变换还是实偶函数
- 奇信号的傅里叶 (正) 变换还是奇函数
- 实奇信号的傅里叶 (正) 变换是虚奇函数
用表格总结奇偶性、是否为实函数的不同情况下的傅里叶正变换奇偶情况:
x(t) | 偶函数:x(t)=x(−t) | 奇函数:x(t)=−x(−t) |
---|
一般的 | X(jω)=X(−jω) | X(jω)=−X(−jω) |
x(t) 是实函数:x(t)=x∗(t) | X(jω)=X(−jω)=X∗(jω)=X∗(−jω) | X(jω)=−X(−jω)=X∗(−jω) |
# 奇偶分解
实信号 x(t) 若有奇偶分解 x(t)=xe(t)+xo(t) ,则有:
X(jω)=Xe(jω)+Xo(jω)xe(t)⟷FXe(−jω)xo(t)⟷FXo(−jω)
# 时域微分特性
可将微分运算转变为代数运算
x(t)⟷FX(jω)⇓dtdx(t)⟷FjωX(jω)
可以推广到高次微分,在后面的微分方程分析 LTI 系统中会用到:
dtkdkx(t)⟷F(jω)kX(jω)
证明:傅里叶反变换两侧同时对 t 微分即可
x(t)dtdx(t)=2π1∫−∞+∞X(jω)ejωtdω=2π1∫−∞+∞X(jω)dtdejωtdω=2π1∫−∞+∞jωX(jω) ejωtdω
# 时域积分特性
可将积分运算转变为代数运算
x(t)⟷FX(jω)⇓∫−∞tx(τ)dτ⟷Fjω1X(jω)+πX(0)δ(ω)
证明:傅里叶正变换两侧同时对 τ 积分即可
F{∫−∞tx(τ)dτ}=∫−∞+∞[∫−∞tx(τ)dτ]e−jωtdt=∫−∞+∞[∫−∞+∞x(τ)u(t−τ)dτ]e−jωtdt=∫−∞+∞x(τ)[∫−∞+∞e−jωtu(t−τ)dt]dτ=∫−∞+∞x(τ)[∫−∞+∞e−jω(t−τ)−jωτu(t−τ)d(t−τ)]dτ=∫−∞+∞x(τ)e−jωτ[∫−∞+∞e−jωsu(s)ds]dτ=[jω1+πδ(ω)]∫−∞+∞x(τ)e−jωτdτ=[jω1+πδ(ω)]X(jω)=jω1X(jω)+πX(0)δ(ω)
其中 u(s)⟷CTFTjω1+πδ(ω) 是单位阶跃函数的傅里叶变换
# 时域和频域的尺度变换
信号如果在时域扩展 a 倍,则其带宽相应压缩 a 倍,反之亦然
x(t)⟷FX(jω)⇓x(at)⟷F∣a∣1X(ajω)
证明:
F{x(at)}=∫−∞+∞x(at)e−jωtdt=∣a∣1∫−∞+∞x(τ)e−jaωτdτ=∣a∣1X(ajω)
# 对偶性
由对偶性可以方便地将时域的某些特性对偶到频域
x(t)⟷CTFTX(jω)⇓X(jt)⟷CTFT2πx(−ω)
证明:
x(t)2πx(t)2πx(ω)2πx(−ω)=2π1∫−∞+∞X(jω)ejωtdω=∫−∞+∞X(jω)ejωtdω=∫−∞+∞X(jt)ejωtdt=∫−∞+∞X(jt)e−jωtdt
# 移频特性
将时域的时移性对偶可以得到频域的移频特性
x(t)⟷FX(jω)⇓x(t)ejω0t⟷FX[j(ω−ω0)]
证明:
x(t)⟷FX(jω)⇓对偶性X(jt)⟷F2πx(−ω)⇓时移特性X(j(t−ω0))⟷F2πx(−ω)e−jω0t⇓对偶性2πx(−t)e−jω0t⟷F2πX[j(−ω−ω0)]⇓x(t)ejω0t⟷FX[j(ω−ω0)]
# 频域微分特性
x(t)⟷FX(jω)⇓−jtx(t)⟷FdωdX(jω)
# 频域积分特性
x(t)⟷FX(jω)⇓−jtx(t)+πx(0)δ(t)⟷F∫−∞ωX(jτ)dτ
# Parseval 定理
x(t)⟷FX(jω)⇓∫−∞+∞∣x(t)∣2dt=2π1∫−∞+∞∣X(jω)∣2dω
这表明:信号的能量既可以在时域求得,也可以在频域求得。由于 ∣X(jω)2∣ 表示了信号能量在频域的分布,因而称其为能量谱密度函数。
# ⭐️卷积特性
x(t)⟷FX(jω) , h(t)⟷FH(jω)⇓x(t)∗h(t)⟷FX(jω)H(jω)
本质上,卷积特性的成立正是因为复指数信号 ejωt 是一切 LTI 系统的特征函数,H(jω) 则是对应的特征值
证明
令 y(t)=x(t)∗h(t)=∫−∞+∞x(τ) h(t−τ) dτ ,求其傅里叶正变换为:
F{y(t)}=Y(jω)=∫−∞+∞e−jωt∫−∞+∞x(τ) h(t−τ) dτ dt=∫−∞+∞e−jω(t−τ)−jωτ∫−∞+∞x(τ) h(t−τ) dτ d(t−τ)=∫−∞+∞e−jωs−jωτ∫−∞+∞x(τ) h(s) dτ ds=∫−∞+∞e−jωse−jωτ∫−∞+∞x(τ) h(s) dτ ds=∫−∞+∞x(τ) e−jωτ dτ∫−∞+∞h(s) e−jωs ds=X(jω)H(jω)
另外,可将 x(t) 视为复指数分量 ejωt 的线性组合,每个 ejωt 通过 LTI 系统时都要受到系统对应的特征值 H(jω) 的加权
LTI 系统的单位冲激响应 h(t) 与频率响应 H(jω) 一一对应,因而 LTI 系统可以由其频率响应完全表征。
⚠️ 并非任何系统的频率响应都存在。
用频率响应表征系统时,一般都限于对稳定系统。
稳定性保证了 h(t) 的绝对可积,结合 Dirichlet 的另外两个条件(实际系统一般满足),则频率响应必然存在
# ⭐️相乘特性
s(t)⟷FS(jω) , p(t)⟷FP(jω)⇓s(t)p(t)⟷F2π1S(jω)∗P(jω)
证明:利用对偶性可以从卷积性质得出相乘性质
由对偶性有:
s(t)⟷FS(jω) , p(t)⟷FP(jω)⇓S(jt)⟷F2πS(−ω) , P(jt)⟷F2πP(−ω)
由卷积特性有:
S(jt)∗P(jt)⟷F4π2s(−ω)p(−ω)
由对偶性有:
4π2s(−t)p(−t)s(−t)p(−t)⟷F2πS(−jω)∗P(−jω)⟷F2π1S(−jω)∗P(−jω)
# 连续时间 LTI 系统的分析
利用傅里叶变换及其特性对 LTI 系统分析
# 频域分析法
根据卷积特性,可以对 LTI 系统进行频域分析,其过程为:
- 求出(或已知、根据系统描述直接得出)系统单位冲激响应 h(t);
- 求输入 x(t) 的傅里叶正变换:x(t)⟷CTFTX(jω);
- 求单位冲激响应 h(t) 的傅里叶正变换:h(t)⟷CTFTH(jω);
- 若需要求 Y(jω) ,由傅里叶变换卷积特性得:Y(jω)=X(jω)H(jω);
- 若需要求 y(t):
- x(t) 是非周期的(或非周期部分):
- 往往由傅里叶反变换求得:y(t)=F−1{Y(jω)}=2π1∫−∞+∞Y(jω)ejωtdω;
- 对有理函数求傅里叶反变换通常采用部分分式展开和利用常用变换对进行
- x(t) 是周期的(或周期部分):
- ∣Y(jω)∣=∣X(jω)∣×∣H(jω)∣;
- argY(jω)=argX(jω)+argH(jω);
# 线性常微分方程
* 线性常系数微分方程
对于由线性常系数微分方程描述的 LTI 系统:
k=0∑Nakdtkdky(t)=k=0∑Mbkdtkdkx(t)
对方程两边进行傅里叶变换,有:
k=0∑Nak(jω)kY(jω)Y(jω)k=0∑Nak(jω)kX(jω)Y(jω)=k=0∑Mbk(jω)kX(jω)=X(jω)k=0∑Mbk(jω)k=∑k=0Nak(jω)k∑k=0Mbk(jω)k
由 CTFT 卷积特性有:
Y(jω)=X(jω)H(jω)
可得:
H(jω)=X(jω)Y(jω)=∑k=0Nak(jω)k∑k=0Mbk(jω)k
可见由线性常系数微分方程描述的 LTI 系统的频率响应是一个有理函数。
若要求 h(t) ,也就是系统的单位冲激响应,则往往可以对有理函数(有理分式)形式的 H(jω) 变形化得到常用变换对的表示形式,再利用常用变换对反变换得到 h(t) 。
反之,若已知一个 LTI 系统的频率响应 H(jω)=X(jω)Y(jω) ,则可以列出两侧分别是关于 Y(jω) 和X(jω) 的方程,两侧同时进行傅里叶反变换,(注意将 jω 视为一个整体)则可得到描述该 LTI 系统的线性常系数微分方程