一般信号应该如何分解,傅里叶级数还是傅里叶变换?基本信号作为激励时,LTI系统如何响应?一般信号作为激励时,LTI系统如何响应?频域分析法的原理、工具、步骤、优缺点是什么?这是本文致力阐述的问题。
1. 基本信号 e j ω t e^{j\omega t} ejωt 作用于LTI系统的响应
设LTI系统的冲激响应为
h
(
t
)
h(t)
h(t), 当激励是角频率
ω
\omega
ω 的基本信号
e
j
ω
t
e^{j\omega t}
ejωt 时,其响应(频域分析中,基本信号的定义域是在负无穷到正无穷区间内的,而在负无穷处可以认为系统的状态为0, 因此这里的响应指的是全响应,也等于零状态响应)为:
y
(
t
)
=
h
(
t
)
∗
e
j
ω
t
y(t)=h(t)*e^{j\omega t}
y(t)=h(t)∗ejωt
根据卷积定义,得:
y
(
t
)
=
∫
−
∞
∞
h
(
τ
)
e
j
ω
(
t
−
τ
)
d
τ
=
∫
−
∞
∞
h
(
τ
)
e
−
j
ω
τ
d
τ
⋅
e
j
ω
t
y(t)=\int_{-\infty}^{\infty} h(\tau) \mathrm{e}^{j \omega(t-\tau)} \mathrm{d} \tau =\int_{-\infty}^{\infty} h(\tau) \mathrm{e}^{-j \omega \tau} \mathrm{d} \tau \cdot \mathrm{e}^{j \omega t}
y(t)=∫−∞∞h(τ)ejω(t−τ)dτ=∫−∞∞h(τ)e−jωτdτ⋅ejωt
h
(
t
)
h(t)
h(t) 的傅里叶变换,记为
H
(
j
ω
)
H(j\omega)
H(jω) ,常称为系统的频率响应函数,也称为传输函数:
H
(
j
ω
)
=
F
[
h
(
t
)
]
=
∫
−
∞
∞
h
(
τ
)
e
−
j
ω
τ
d
τ
H(j \omega)=\mathrm{F}[h(t)]=\int_{-\infty}^{\infty} h(\tau) \mathrm{e}^{-j \omega \tau} \mathrm{d} \tau
H(jω)=F[h(t)]=∫−∞∞h(τ)e−jωτdτ
所以,基本信号
e
j
ω
t
e^{j\omega t}
ejωt 作用于LTI系统的响应为:
y
(
t
)
=
H
(
j
ω
)
⋅
e
j
ω
t
y(t)=H(j\omega)\cdot e^{j\omega t}
y(t)=H(jω)⋅ejωt
显然,
H
(
j
ω
)
H(j\omega)
H(jω) 反映了响应
y
(
t
)
y(t)
y(t) 的幅度和相位。
2. 一般信号 f ( t ) f(t) f(t) 作用于LTI系统的响应
基本信号
e
j
ω
t
e^{j\omega t}
ejωt 作用于LTI系统的响应可表示为:
e
j
ω
t
⟹
H
(
j
ω
)
e
j
ω
t
e^{j \omega t} \Longrightarrow H(j \omega) e^{j \omega t}
ejωt⟹H(jω)ejωt
由齐次性:
1
2
π
F
(
j
ω
)
d
ω
e
j
ω
t
⟹
1
2
π
F
(
j
ω
)
d
ω
H
(
j
ω
)
e
j
ω
t
\frac{1}{2 \pi} {F}(j \omega)\: \mathrm{d} \omega \: e^{j \omega t} \Longrightarrow \frac{1}{2 \pi} {F}(j \omega)\: \mathrm{d} \omega \:H(j \omega) e^{j \omega t}
2π1F(jω)dωejωt⟹2π1F(jω)dωH(jω)ejωt
由可叠加性:
1
2
π
∫
−
∞
∞
F
(
j
ω
)
e
j
ω
t
d
ω
⟹
1
2
π
∫
−
∞
∞
H
(
j
ω
)
F
(
j
ω
)
e
j
ω
t
d
ω
\frac{1}{2 \pi} \int_{-\infty}^{\infty} {F}(j \omega) \: e^{j \omega t} \: \mathrm{d} \omega\Longrightarrow \frac{1}{2 \pi} \int_{-\infty}^{\infty} H(j \omega) {F}(j \omega)\:e^{j \omega t}\: \mathrm{d} \omega
2π1∫−∞∞F(jω)ejωtdω⟹2π1∫−∞∞H(jω)F(jω)ejωtdω
由傅里叶逆变换的公式可知,上式左侧即为
f
(
t
)
f(t)
f(t) ;而上式右侧为
H
(
j
ω
)
F
(
j
ω
)
H(j \omega) {F}(j \omega)
H(jω)F(jω) 的傅里叶逆变换,所以上式可以表示为:
f
(
t
)
⟹
F
−
1
[
H
(
j
ω
)
F
(
j
ω
)
]
f(t) \Longrightarrow F^{-1}[H(j \omega) {F}(j \omega)]
f(t)⟹F−1[H(jω)F(jω)]
因此,一般新号
f
(
t
)
f(t)
f(t) 作用于LTI系统的响应为:
y
(
t
)
=
F
−
1
[
H
(
j
ω
)
F
(
j
ω
)
]
y(t) = F^{-1}[H(j \omega) {F}(j \omega)]
y(t)=F−1[H(jω)F(jω)]
对
y
(
t
)
y(t)
y(t) 进行傅里叶变换有:
Y
(
t
)
=
H
(
j
ω
)
F
(
j
ω
)
Y(t) = H(j \omega) {F}(j \omega)
Y(t)=H(jω)F(jω)
所以可以这么描述:对于一般信号 ,其作用于LTI系统的响应的傅里叶变换为该新号的傅里叶变换与系统的传输函数(频率响应函数)的乘积。这不正是 时域卷积对应频域相乘 吗?
3. 傅里叶变换分析法
从上面可知,在求一个系统的响应时,可以分为下面4个步骤来完成:
- 求输入信号 f ( t ) f(t) f(t) 的傅里叶变换 F ( j ω ) F(j \omega) F(jω) ;
- 求系统的传输函数(频率响应函数);
- 求零状态响应 y ( t ) y(t) y(t) 的傅里叶变换 Y ( j ω ) = H ( j ω ) F ( j ω ) Y(j\omega) = H(j \omega) {F}(j \omega) Y(jω)=H(jω)F(jω) ;
- 求 Y ( j ω ) Y(j\omega) Y(jω) 的傅里叶逆变换 y ( t ) = F − 1 [ H ( j ω ) F ( j ω ) ] y(t) = F^{-1}[H(j \omega) {F}(j \omega)] y(t)=F−1[H(jω)F(jω)] .
4. 傅里叶级数分析法
无论是周期信号还是非周期信号我们均可以采用傅里叶变换来进行频域分析,但是周期信号也可以用傅里叶级数来进行频域分析。因此,如果掌握了傅里叶变换分析法,按道理来说就足够了。但是也不能说傅里叶级数分析法鸡肋,因为有时候工程追求的是最优解法,所以多掌握几种方法肯定是不错的。
4.1 指数形式的傅里叶级数
周期信号
f
T
(
t
)
f_T(t)
fT(t) 的指数形式的傅里叶级数可以表示为:
f
T
(
t
)
=
∑
n
=
−
∞
∞
F
n
e
j
n
Ω
t
f_T(t) =\sum_{n=-\infty}^{\infty}F_n\rm e^{\:\it j n\Omega t}
fT(t)=n=−∞∑∞FnejnΩt
系统的零状态响应可以表示为:
y
(
t
)
=
h
(
t
)
∗
f
T
(
t
)
=
∑
n
=
−
∞
∞
F
n
[
h
(
t
)
∗
e
j
n
Ω
t
]
=
∑
n
=
−
∞
∞
F
n
H
(
j
n
Ω
)
e
j
n
Ω
t
y(t)=h(t) * f_{T}(t)=\sum_{n=-\infty}^{\infty} F_{n}\left[h(t) * e^{j n \Omega t}\right]=\sum_{n=-\infty}^{\infty} F_{n} H(j n \Omega) e^{j n \Omega t}
y(t)=h(t)∗fT(t)=n=−∞∑∞Fn[h(t)∗ejnΩt]=n=−∞∑∞FnH(jnΩ)ejnΩt
上式中,
F
n
H
(
j
n
Ω
)
F_{n} H(j n \Omega)
FnH(jnΩ) 即为系统的零状态响应
y
(
t
)
y(t)
y(t) 的傅里叶系数。
因此,傅里叶级数分析法也可以总结为以下4个步骤:
- 求周期输入信号 f T ( t ) f_{T}(t) fT(t) 的傅里叶系数 F n F_n Fn ;
- 求系统的传输函数(频率响应函数)离散化后的取值序列 H ( j n Ω ) = H ( j ω ) ∣ ω = n Ω H(j n \Omega)=H(j \omega)|_{\omega=n\Omega} H(jnΩ)=H(jω)∣ω=nΩ;
- 求零状态响应 y ( t ) y(t) y(t) 的傅里叶系数 F n H ( j n Ω ) F_{n} H(j n \Omega) FnH(jnΩ) ;
- 求傅里叶级数的展开式 y ( t ) y(t) y(t) .
4.2 三角形式的傅里叶级数
周期信号
f
T
(
t
)
f_T(t)
fT(t) 的三角形式的傅里叶级数可以表示为:
f
T
(
t
)
=
A
0
2
+
∑
n
=
1
∞
A
n
c
o
s
(
n
Ω
t
+
φ
n
)
f_{T}(t)=\frac{A_{0}}{2}+\sum_{n=1}^{\infty} A_{n} \rm cos \left(n \Omega t+\varphi_{n}\right)
fT(t)=2A0+n=1∑∞Ancos(nΩt+φn)
系统的传输函数还可以表示为模值和相角相乘的形式:
H
(
j
ω
)
=
∣
H
(
j
ω
)
∣
e
j
θ
(
ω
)
H(j \omega)=|H(j \omega)| e^{j \theta(\omega)}
H(jω)=∣H(jω)∣ejθ(ω)
那么传输函数在傅里叶级数基频和各谐频处的离散化取值为:
H
(
j
n
Ω
)
=
∣
H
(
j
n
Ω
)
∣
e
j
θ
(
n
Ω
)
=
H
(
j
ω
)
∣
ω
=
n
Ω
H(j n \Omega)=|H(j n \Omega)| e^{j \theta(n \Omega)}=\left.H(j \omega)\right|_{\omega=n \Omega}
H(jnΩ)=∣H(jnΩ)∣ejθ(nΩ)=H(jω)∣ω=nΩ
类比指数形式表达的零状态响应,容易得到当周期信号作为激励时,系统的零状态响应为:
y
(
t
)
=
A
0
2
H
(
0
)
⏟
直流增益
+
∑
n
=
1
∞
A
n
∣
H
(
j
n
Ω
)
∣
⏟
幅值相乘
cos
[
n
Ω
t
+
ϕ
n
+
θ
(
n
Ω
)
⏟
相位相加
]
y(t)=\frac{A_{0}}{2}\underbrace{H(0)}_{\text{直流增益 }} +\sum_{n=1}^{\infty} \underbrace{A_{n} \mid H(j n \Omega)\mid}_{\text {幅值相乘 }} \cos \:[n \Omega t+\underbrace{\phi_{n}+\theta(n \Omega)}_{\text {相位相加 }}]
y(t)=2A0直流增益
H(0)+n=1∑∞幅值相乘
An∣H(jnΩ)∣cos[nΩt+相位相加
ϕn+θ(nΩ)]
5. 实践(做题是掌握概念的有效检验方式)
若某LTI系统的幅频特性 ∣ H ( j ω ) ∣ |H(j\omega)| ∣H(jω)∣ 和相频特性 θ ( ω ) \theta(\omega) θ(ω) 如下图所示,若激励为 f ( t ) = 2 + 4 cos ( 5 t ) + 4 cos ( 10 t ) f(t)=2+4 \cos (5 t)+4 \cos (10 t) f(t)=2+4cos(5t)+4cos(10t) , 如何求系统的响应?
下面将分别用傅里叶级数分析法和傅里叶系数分析法来求解系统响应。
5.1 傅里叶变换分析法解答
系统的传输函数是已知的,为:
H
(
j
ω
)
=
∣
H
(
j
ω
)
∣
e
j
θ
(
ω
)
H(j \omega)=|H(j \omega)| e^{j \theta(\omega)}
H(jω)=∣H(jω)∣ejθ(ω)
激励函数
f
(
t
)
f(t)
f(t) 的傅里叶变换为:
F
(
j
ω
)
=
4
π
δ
(
ω
)
+
4
π
[
δ
(
ω
+
5
)
+
δ
(
ω
−
5
)
]
+
4
π
[
δ
(
ω
+
10
)
+
δ
(
ω
−
10
)
]
F(j \omega)=4 \pi \delta(\omega)+4 \pi[\delta(\omega+5)+\delta(\omega-5)]+4 \pi[\delta(\omega+10)+\delta(\omega-10)]
F(jω)=4πδ(ω)+4π[δ(ω+5)+δ(ω−5)]+4π[δ(ω+10)+δ(ω−10)]
则可以求出响应的傅里叶变换为:(把下面式子中的分数都变成小数)
Y
(
j
ω
)
=
F
(
j
ω
)
H
(
j
ω
)
=
4
π
δ
(
ω
)
+
4
π
⋅
δ
(
ω
+
5
)
⋅
H
(
−
5
j
)
+
4
π
⋅
δ
(
ω
−
5
)
⋅
H
(
5
j
)
+
4
π
⋅
δ
(
ω
+
10
)
⋅
H
(
−
10
j
)
+
4
π
⋅
δ
(
ω
−
10
)
⋅
H
(
10
j
)
=
4
π
δ
(
ω
)
+
4
π
⋅
0.5
⋅
e
0.5
π
⋅
δ
(
ω
+
5
)
+
4
π
⋅
0.5
⋅
e
−
0.5
π
⋅
δ
(
ω
−
5
)
+
0
+
0
=
4
π
δ
(
ω
)
+
j
⋅
2
π
⋅
δ
(
ω
+
5
)
−
j
⋅
2
π
⋅
δ
(
ω
−
5
)
\begin{aligned} Y(j \omega) &=F(j \omega) H(j \omega) \\ &=4 \pi \delta(\omega)+4 \pi\cdot \delta(\omega+5)\cdot H(-5j) +4 \pi\cdot \delta(\omega-5)\cdot H(5j) +4 \pi\cdot\delta(\omega+10) \cdot H(-10j) +4 \pi\cdot \delta(\omega-10)\cdot H(10j) \\ &=4 \pi \delta(\omega)+4 \pi\cdot 0.5 \cdot e^{0.5{\pi}}\cdot \delta(\omega+5)+4 \pi\cdot 0.5 \cdot e^{-0.5{\pi}}\cdot \delta(\omega-5)+0+0\\ &=4 \pi \delta(\omega)+j\cdot 2 \pi \cdot \delta(\omega+5)-j\cdot 2 \pi \cdot \delta(\omega-5) \end{aligned}
Y(jω)=F(jω)H(jω)=4πδ(ω)+4π⋅δ(ω+5)⋅H(−5j)+4π⋅δ(ω−5)⋅H(5j)+4π⋅δ(ω+10)⋅H(−10j)+4π⋅δ(ω−10)⋅H(10j)=4πδ(ω)+4π⋅0.5⋅e0.5π⋅δ(ω+5)+4π⋅0.5⋅e−0.5π⋅δ(ω−5)+0+0=4πδ(ω)+j⋅2π⋅δ(ω+5)−j⋅2π⋅δ(ω−5)
因为正弦函数
sin
(
ω
0
t
)
\sin \left(\omega_{0} t\right)
sin(ω0t) 的傅里叶变换为:
sin
(
ω
0
t
)
=
1
2
j
(
e
j
ω
0
t
−
e
−
j
ω
0
t
)
⟷
j
π
[
δ
(
ω
+
ω
0
)
−
δ
(
ω
−
ω
0
)
]
\sin \left(\omega_{0} t\right)=\frac{1}{2 j}\left(\mathrm{e}^{j \omega_{0} t}-\mathrm{e}^{-j \omega_{0} t}\right) \longleftrightarrow j \pi\left[\delta\left(\omega+\omega_{0}\right)-\delta\left(\omega-\omega_{0}\right)\right]
sin(ω0t)=2j1(ejω0t−e−jω0t)⟷jπ[δ(ω+ω0)−δ(ω−ω0)]
所以,对
Y
(
j
ω
)
Y(j\omega)
Y(jω) 进行傅里叶逆变换,可以得到响应为:
y
(
t
)
=
F
−
1
[
Y
(
j
ω
)
]
=
2
+
2
sin
5
t
y(t)=\mathrm{F}^{-1}[Y(j \omega)]=2+2 \sin 5 t
y(t)=F−1[Y(jω)]=2+2sin5t
5.2 傅里叶级数分析法解答
激励函数的基波角频率为
Ω
=
5
r
a
d
/
s
\Omega=5\: \rm rad/s
Ω=5rad/s , 所以激励函数可以写为:
f
(
t
)
=
2
+
4
cos
(
Ω
t
)
+
4
cos
(
2
Ω
t
)
f(t)=2+4 \cos (\Omega t)+4 \cos (2 \Omega t)
f(t)=2+4cos(Ωt)+4cos(2Ωt)
观察易得:
H
(
0
)
=
1
,
H
(
j
Ω
)
=
0.5
e
−
j
0.5
π
,
H
(
j
2
Ω
)
=
0
{H}(0)=1, \quad H(j \Omega)=0.5 \mathrm{e}^{-j 0.5 \pi}, \quad {H}({j} 2 \Omega)=0
H(0)=1,H(jΩ)=0.5e−j0.5π,H(j2Ω)=0
根据本文公式可以直接求得响应为:
y
(
t
)
=
A
0
2
H
(
0
)
+
∑
n
=
1
∞
A
n
∣
H
(
j
n
Ω
)
∣
cos
[
n
Ω
t
+
ϕ
n
+
θ
(
n
Ω
)
]
=
2
×
1
+
4
×
0.5
cos
(
Ω
t
−
0.5
π
)
+
4
×
0
=
2
+
2
sin
5
t
\begin{aligned} y(t) &=\frac{A_{0}}{2} H(0)+\sum_{n=1}^{\infty} A_{n}|H(j n \Omega)| \cos \left[n \Omega t+\phi_{n}+\theta(n \Omega)\right] \\ &=2 \times 1+4 \times 0.5 \cos (\Omega t-0.5 \pi)+4 \times 0 \\ &=2+2 \sin 5 t \end{aligned}
y(t)=2A0H(0)+n=1∑∞An∣H(jnΩ)∣cos[nΩt+ϕn+θ(nΩ)]=2×1+4×0.5cos(Ωt−0.5π)+4×0=2+2sin5t
本文首发于微信公众号 振动信号研究所