系统稳定性判定分析(四)---- 传递函数前情提要

发布于:2024-05-19 ⋅ 阅读:(150) ⋅ 点赞:(0)

传递函数前情提要

傅里叶变换与拉普拉斯变换

傅里叶变换的目的是可将时域(即时间域)上的信号转变为频域(即频率域)上的信号,随着域的不同,对同一个事物的了解角度也就随之改变,因此在时域中某些不好处理的地方,在频域就可以较为简单的处理。

拉普拉斯变换可以说是现代工程学使用最广泛的数学工具,它可以通过数学变换将微积分方程转化成代数方程。(重要人物:赫维赛德)

在这里插入图片描述

傅里叶定律 只要一个函数满足狄利赫里条件,则其可以分解为复指数函数之和。

狄利赫里条件
(1) 函数在任意有限区间内连续,或只有有限个第一类间断点;
(2) 在一个周期内,函数有有限个极大值或极小值;
(3) 函数 f ( t ) f(t) f(t)绝对可积,即有 ∫ 0 + ∞ ∣ f ( t ) ∣ d t < ∞ \int_0^{+\infty}|f(t)|dt<\infty 0+f(t)dt<.

当函数 f ( t ) f(t) f(t)满足狄利赫里条件时,其傅里叶变换及其逆变换为:
F ( ω ) = ∫ − ∞ + ∞ f ( t ) e − i ω t d t , (1) F(\omega)=\int_{-\infty}^{+\infty}f(t)e^{-i\omega t}dt,\tag{1} F(ω)=+f(t)etdt,(1) f ( t ) = 1 2 π ∫ − ∞ + ∞ F ( ω ) e i ω t d ω , (2) f(t)=\frac{1}{2\pi}\int_{-\infty}^{+\infty}F(\omega)e^{i\omega t}d\omega,\tag{2} f(t)=2π1+F(ω)etdω,(2)
添加衰减函数 e − σ t e^{-\sigma t} eσt,并将时间域的范围缩减为正半轴,则傅里叶变换变为:
F ( ω ) = ∫ 0 + ∞ f ( t ) e − σ t e − i ω t d t = ∫ 0 + ∞ f ( t ) e − ( σ + i ω ) t d t F(\omega)=\int_0^{+\infty}f(t)e^{-\sigma t}e^{-i\omega t}dt=\int_0^{+\infty}f(t)e^{-(\sigma+i\omega) t}dt F(ω)=0+f(t)eσtetdt=0+f(t)e(σ+)tdt s = σ + i ω s=\sigma+i\omega s=σ+,则得拉普拉斯变换
F ( s ) = ∫ 0 + ∞ f ( t ) e − s t d t , (3) F(s)=\int_0^{+\infty}f(t)e^{-st}dt,\tag{3} F(s)=0+f(t)estdt,(3)

线性时不变系统的优点

线性时不变系统是一种特殊的线性系统,其继承了线性系统的叠加性,同时系统的参数不随时间变化,亦即信号作用的前后只影响响应输出的先后而不影响形状。

控制,要求我们精确知道每个时刻系统的状态是什么,而系统当前的状态又由系统的初始状态以及之前所有时刻的输入决定。对于线性系统其满足叠加定理,因此要确定当前时刻系统的状态,我们可以将之前所有时刻的输入带来的当前时刻响应叠加,再加上系统初始状态即可。

卷积

已知系统的初始状态以及之前所有时刻的输入信号,可通过卷积计算线性时不变系统的任意时刻的输出响应。两函数的卷积定义如下:
f 1 ( t ) ∗ f 2 ( t ) = ∫ − ∞ + ∞ f 1 ( τ ) f 2 ( t − τ ) d τ , (4) f_1(t)*f_2(t)=\int_{-\infty}^{+\infty}f_1(\tau)f_2(t-\tau)d\tau,\tag{4} f1(t)f2(t)=+f1(τ)f2(tτ)dτ,(4)其中 f 1 ( t ) ∗ f 2 ( t ) f_1(t)*f_2(t) f1(t)f2(t)代表函数 f 1 ( t ) f_1(t) f1(t)与函数 f 2 ( t ) f_2(t) f2(t)的卷积。关于卷积的更多介绍及解释可参见如何通俗易懂地解释卷积下的回答及互动内容。

关于卷积的自我小思想: 假定 f ( t ) = f 1 ( t ) ∗ f 2 ( t ) f(t)=f_1(t)*f_2(t) f(t)=f1(t)f2(t)表示 t t t时刻一植株的生长高度, f 1 ( τ ) f_1(\tau) f1(τ)表示在 τ \tau τ时刻的生长速度, f 2 ( t − τ ) f_2(t-\tau) f2(tτ)表示 τ \tau τ时刻植株生长速度对于 t t t时刻植株生长高度的影响比重。植株持续性生长,那么植株 t t t时刻的生长高度由植株 t t t时刻之前的所有时刻的生长速度决定,每一时刻的生长速度对植株 t t t时刻的生长高度影响程度不同,其影响程度对应卷积距离,具体表现为: 离得近的信号,卷积距离短;离得远的信号,卷积距离长。

卷积定理

卷积定理是傅里叶变换满足的一个重要性质。卷积定理指出,函数卷积的傅里叶变换是函数傅里叶变换的乘积。具体分为时域卷积定理和频域卷积定理。这一定理对拉普拉斯变换双边拉普拉斯变换Z变换Mellin变换Hartley变换等各种傅里叶变换的变体同样成立。

定义 F \mathcal{F} F表示傅里叶变换,若函数 f 1 ( t ) f_1(t) f1(t)的傅里叶变换为 F 1 ( ω ) \mathcal{F}_1(\omega) F1(ω), 函数 f 2 ( t ) f_2(t) f2(t)的傅里叶变换为 F 2 ( ω ) \mathcal{F}_2(\omega) F2(ω),则

时域卷积定理 两信号在时域的卷积积分对应于在频域中该两信号的傅里叶变换的乘积,即有 F [ f 1 ( t ) ∗ f 2 ( t ) ] = F 1 ( ω ) × F 2 ( ω ) . \mathcal{F}[f_1(t)*f_2(t)]=\mathcal{F}_1(\omega)\times \mathcal{F}_2(\omega). F[f1(t)f2(t)]=F1(ω)×F2(ω).

频域卷积定理 两信号在时域的乘积对应于该两信号傅里叶变换的卷积的傅里叶逆变换除以 2 π 2\pi 2π,即有 f 1 ( t ) × f 2 ( t ) = 1 2 π I F [ F 1 ( ω ) ∗ F 2 ( ω ) ] . f_1(t)\times f_2(t)=\frac{1}{2\pi}\mathcal{IF}[\mathcal{F}_1(\omega)*\mathcal{F}_2(\omega)]. f1(t)×f2(t)=2π1IF[F1(ω)F2(ω)].

证明

时域卷积定理证明 已知 F \mathcal{F} F为傅里叶变换,则 F [ f 1 ( t ) ∗ f 2 ( t ) ] = ∫ − ∞ + ∞ [ ∫ − ∞ + ∞ f 1 ( τ ) f 2 ( t − τ ) d τ ] e − i w t d t = ∫ − ∞ + ∞ f 1 ( τ ) [ ∫ − ∞ + ∞ f 2 ( t − τ ) e − i w t d t ] d τ = ∫ − ∞ + ∞ f 1 ( τ ) e − i w τ [ ∫ − ∞ + ∞ f 2 ( t − τ ) e − i w ( t − τ ) d ( t − τ ) ] d τ = ∫ − ∞ + ∞ f 1 ( τ ) e − i w τ F 2 ( ω ) d τ = F 2 ( ω ) ∫ − ∞ + ∞ f 1 ( τ ) e − i w τ d τ = F 2 ( ω ) × F 1 ( ω ) . \begin{align*}\mathcal{F}[f_1(t)*f_2(t)]&=\int_{-\infty}^{+\infty}[\int_{-\infty}^{+\infty}f_1(\tau)f_2(t-\tau)d\tau]e^{-iwt}dt\\&=\int_{-\infty}^{+\infty}f_1(\tau)[\int_{-\infty}^{+\infty}f_2(t-\tau)e^{-iwt}dt]d\tau\\&=\int_{-\infty}^{+\infty}f_1(\tau)e^{-iw\tau}[\int_{-\infty}^{+\infty}f_2(t-\tau)e^{-iw(t-\tau)}d(t-\tau)]d\tau\\&=\int_{-\infty}^{+\infty}f_1(\tau)e^{-iw\tau}\mathcal{F}_2(\omega)d\tau\\&=\mathcal{F}_2(\omega)\int_{-\infty}^{+\infty}f_1(\tau)e^{-iw\tau}d\tau\\&=\mathcal{F}_2(\omega)\times\mathcal{F}_1(\omega).\end{align*} F[f1(t)f2(t)]=+[+f1(τ)f2(tτ)dτ]eiwtdt=+f1(τ)[+f2(tτ)eiwtdt]dτ=+f1(τ)eiwτ[+f2(tτ)eiw(tτ)d(tτ)]dτ=+f1(τ)eiwτF2(ω)dτ=F2(ω)+f1(τ)eiwτdτ=F2(ω)×F1(ω).证毕。

频域卷积定理证明 函数 f 1 ( t ) f_1(t) f1(t)的傅里叶变换为 F 1 ( ω ) \mathcal{F}_1(\omega) F1(ω), 函数 f 2 ( t ) f_2(t) f2(t)的傅里叶变换为 F 2 ( ω ) \mathcal{F}_2(\omega) F2(ω),根据公式(2)有
I F [ F 1 ( ω ) ∗ F 2 ( ω ) ] = 1 2 π ∫ − ∞ + ∞ [ ∫ − ∞ + ∞ F 1 ( μ ) F 2 ( ω − μ ) d μ ] e i ω t d ω = 1 2 π ∫ − ∞ + ∞ ∫ − ∞ + ∞ F 1 ( μ ) e i μ t F 2 ( ω − μ ) e i ( ω − μ ) t d ω d μ = 1 2 π ∫ − ∞ + ∞ F 1 ( μ ) e i μ t d μ ∫ − ∞ + ∞ F 2 ( ω − μ ) e i ( ω − μ ) t d ( ω − μ ) = f 1 ( t ) × 2 π × f 2 ( t ) = 2 π f 1 ( t ) × f 2 ( t ) . \begin{align*}\mathcal{IF}[\mathcal{F}_1(\omega)*\mathcal{F}_2(\omega)]&=\frac{1}{2\pi}\int_{-\infty}^{+\infty}[\int_{-\infty}^{+\infty}\mathcal{F}_1(\mu)\mathcal{F}_2(\omega-\mu)d\mu]e^{i\omega t}d\omega\\&=\frac{1}{2\pi}\int_{-\infty}^{+\infty}\int_{-\infty}^{+\infty}\mathcal{F}_1(\mu)e^{i\mu t}\mathcal{F}_2(\omega-\mu)e^{i(\omega-\mu) t}d\omega d\mu\\&=\frac{1}{2\pi}\int_{-\infty}^{+\infty}\mathcal{F}_1(\mu)e^{i\mu t}d\mu\int_{-\infty}^{+\infty}\mathcal{F}_2(\omega-\mu)e^{i(\omega-\mu) t}d(\omega-\mu)\\&=f_1(t)\times2\pi\times f_2(t)\\&= 2\pi f_1(t)\times f_2(t).\end{align*} IF[F1(ω)F2(ω)]=2π1+[+F1(μ)F2(ωμ)dμ]etdω=2π1++F1(μ)eiμtF2(ωμ)ei(ωμ)tdωdμ=2π1+F1(μ)eiμtdμ+F2(ωμ)ei(ωμ)td(ωμ)=f1(t)×2π×f2(t)=2πf1(t)×f2(t).证毕。

学得两句英语哲理

  1. A tough life needs no explanation.
  2. To one who has faith, no explanation is necessary. To one without faith, no explanation is possible.

参考文献

[1] 如何通俗易懂地解释卷积?
[2] J Pan:从另一个角度看拉普拉斯变换
[3] 傅里叶变换概念及公式推导


网站公告

今日签到

点亮在社区的每一天
去签到