数学基础
傅里叶分析:时域与频域的桥梁
-
核心思想:把函数表示为多个正弦或余弦函数(复指数函数)的叠加,即傅里叶分解。
- 联想:幂级数分解(泰勒展开)。
- 核心价值:实现了时域视角向频域视角的转换。我们不只看信号随时间如何变化(时域),还看它包含什么样的频率成分(频域)。
- 线性代数视角:任何能量有限的信号(平方可积函数),都可视为无穷维函数空间中的一个向量。它等于一组正交基(不同频率的复指数函数)的线性组合。
傅里叶级数
将一个周期信号,分解成一系列频率为基频整数倍的谐波之和。此时频谱离散。
- ,其中 为基频, 为傅里叶系数,其取值亦为复数,能同时表示谐波的幅度和相位。
-
,为傅里叶系数的求法(从时域函数中提取频域信息)。
- 从复指数函数的角度理解,只有指数互为相反数的两个复指数函数相乘,才能使周期内积分不为0;从三角函数角度理解,只有相同的三角函数相乘,才能使周期内积分不为0。
- 本质是利用正交性进行匹配滤波。让信号与基向量做内积,由基的正交性,只会留下相应方向上的坐标。
傅里叶变换
可以视为傅里叶级数中周期趋于无穷。针对非周期信号,变换为连续的频谱密度函数。
- .
- 逆变换 .
实际应用
- 离散傅里叶变换(DFT):针对数字信号。时域和频域都是离散且有限长的。
- 快速傅里叶变换(FFT):DFT的高效算法,利用周期性和对称性,极大降低了运算的复杂度。
-
实时频谱分析仪:对音频分帧、加窗、FFT、取模后显示。即计算信号在各频率基上的坐标。
- 加窗:在时域乘以窗函数,让帧的两端平滑过渡到零,避免频谱泄露。
-
不确定性原理:时域宽度与频域宽度的乘积具有下界,二者不可兼得。
- 直观理解:短暂的音符难以分辨音高,更长的时间才能带来更精准的频率。
拉普拉斯变换:复频域
理解指数函数:
-
本质特征:函数与导数成线性关系
.𝑓 ′ ( 𝑡 ) = 𝑠 𝑓 ( 𝑡 ) - 比例系数
为复数时,即得到复指数函数。𝑠 - 速度与位矢成线性关系时(固定倍数、夹角),质点的轨迹为复指数函数。
- 比例系数
-
复频率
𝑠 = 𝜎 + 𝑖 𝜔 - 实部
为衰减常数,正负控制增长或衰减,大小控制速度,仅含实部时(速度与位矢共线时)为直线运动。𝜎 - 虚部
为角频率,正负控制旋转方向,大小控制旋转(分量上表现为振荡)速度,仅含虚部时(速度与位矢垂直时)为匀速圆周运动。𝜔
- 实部
- 复振幅
,其中 ˜𝑈 = 𝐴 𝑒 𝑖 𝜑 0 为振幅,𝐴 为初始相位。𝜑 0
常系数齐次线性微分方程的解法:指数函数的运用
-
线性微分方程:未知函数
及其各阶导数𝑦 都是一次的微分方程,满足叠加原理。𝐷 𝑛 𝑦 - 微分算子
.(算子的乘法/幂次定义为算子的复合)𝐷 = 𝑑 𝑑 𝑡 , 𝐷 𝑛 = 𝑑 𝑛 𝑑 𝑡 𝑛 - 线性微分算子
.𝐿 = ∑ 𝑖 = 0 𝑛 𝑎 𝑖 ( 𝑡 ) 𝐷 𝑖 - 线性微分方程可记作
,有𝐿 [ 𝑦 ] = 𝑓 ( 𝑡 ) .𝐿 ( 𝑐 1 𝑦 1 + 𝑐 2 𝑦 2 ) = 𝑐 1 𝐿 ( 𝑦 1 ) + 𝑐 2 𝐿 ( 𝑦 2 )
- 微分算子
-
齐次线性微分方程:
,其解空间是一个线性空间。𝐿 [ 𝑦 ] = 0 - 非齐次线性微分方程的通解 = 对应齐次线性微分方程的通解 + 非齐次线性微分方程的一个特解。
-
常系数线性微分方程:所有系数都是常数,此时
,(𝐿 = 𝑃 ( 𝐷 ) = ∑ 𝑖 = 0 𝑛 𝑎 𝑖 𝐷 𝑖 与常数系数可交换)𝐷 是𝐿 的多项式,可以因式分解,将高阶方程降阶为一系列一阶方程求解。𝐷 - 指数函数
是𝑒 𝜆 𝑡 的特征函数,满足𝑃 ( 𝐷 ) ,在算子作用下只是乘上一个常数。𝑃 ( 𝐷 ) 𝑒 𝜆 𝑡 = 𝑃 ( 𝜆 ) 𝑒 𝜆 𝑡
- 指数函数
-
常系数齐次线性微分方程:
.𝑃 ( 𝐷 ) 𝑦 = 0 - 设
,微分方程转化为其特征方程𝑦 = 𝑒 𝜆 𝑡 ,求出特征根𝑃 ( 𝜆 ) = 0 ,𝜆 𝑖 即为解空间的一组基。特别地,若有𝑒 𝜆 𝑖 𝑡 重根𝑘 ,则𝜆 都为解空间的基。𝑒 𝜆 𝑡 , 𝑡 𝑒 𝜆 𝑡 , … , 𝑡 𝑘 − 1 𝑒 𝜆 𝑡 - 速度与位矢成线性关系的物理系统,对应一阶常系数齐次线性微分方程,其解为一个复指数函数。
- 阻尼振子对应二阶常系数齐次线性微分方程,其一般解可由两个复指数函数线性表出。特征方程在欠阻尼时有两个共轭复根,临界阻尼时有重根,过阻尼时有两个负实根。
- 设
常系数线性微分方程的解法:拉普拉斯变换
拉普拉斯变换:将时域函数映射到复频域,定义为
时,𝑓 ( 𝑡 ) = 𝑒 𝑎 𝑡 ,在极点𝐹 ( 𝑠 ) = ∫ 0 ∞ 𝑒 ( 𝑎 − 𝑠 ) 𝑡 𝑑 𝑡 = 1 𝑠 − 𝑎 处积分发散。(解析延拓后,𝑠 = 𝑎 )𝑠 ∈ ℂ 的每个𝑓 ( 𝑡 ) 分量都表现为𝑒 𝑎 𝑖 𝑡 的一个极点𝐹 ( 𝑠 ) .𝑎 𝑖 ,拉普拉斯变换将微分转化为乘法。ℒ︀ { 𝑓 ′ ( 𝑡 ) } = 𝑠 𝐹 ( 𝑠 ) − 𝑓 ( 0 ) - 对常系数线性微分方程
作拉普拉斯变换,在复频域转化为代数方程𝑃 ( 𝐷 ) 𝑦 = 𝑓 ( 𝑡 ) ,解出𝑃 ( 𝑠 ) 𝑌 ( 𝑠 ) − 𝑄 ( 𝑠 ) = 𝐹 ( 𝑠 ) 后再通过逆变换得到𝑌 ( 𝑠 ) 的时域表达式𝑦 .𝑦 ( 𝑡 ) = ℒ︀ − 1 { 𝑌 ( 𝑠 ) }
几种典型函数的拉普拉斯变换
- 单位脉冲函数
,拉普拉斯变换为𝛿 ( 𝑡 ) .1 - 单位阶跃函数
,拉普拉斯变换为1 ( 𝑡 ) = 1 , ( 𝑡 ≥ 0 ) .1 𝑠 - 单位斜坡函数
,拉普拉斯变换为𝑓 ( 𝑡 ) = 𝑡 , ( 𝑡 ≥ 0 ) .1 𝑠 2 - 单位加速度函数
,拉普拉斯变换为𝑓 ( 𝑡 ) = 𝑡 2 2 , ( 𝑡 ≥ 0 ) .1 𝑠 3 - 指数函数
,拉普拉斯变换为𝑓 ( 𝑡 ) = 𝑒 𝑎 𝑡 , ( 𝑡 ≥ 0 ) .1 𝑠 − 𝑎 - 正弦函数
,拉普拉斯变换为𝑓 ( 𝑡 ) = sin ( 𝜔 𝑡 ) , ( 𝑡 ≥ 0 ) .𝜔 𝑠 2 + 𝜔 2 - 余弦函数
,拉普拉斯变换为𝑓 ( 𝑡 ) = cos ( 𝜔 𝑡 ) , ( 𝑡 ≥ 0 ) .𝑠 𝑠 2 + 𝜔 2
控制系统的数学模型
基本概念
- 线性系统:满足叠加原理的系统。
-
线性时不变系统(线性定常系统):系统的特性不随时间改变,可以用常系数线性微分方程描述。
,其中𝑃 ( 𝐷 ) 𝑦 ( 𝑡 ) = 𝑄 ( 𝐷 ) 𝑥 ( 𝑡 ) 是输入(驱动函数),𝑥 ( 𝑡 ) 是输出(响应函数)。𝑦 ( 𝑡 )
-
传递函数:零初始条件下,系统输出与输入的拉普拉斯变换之比。
-
𝐺 ( 𝑠 ) = 𝑌 ( 𝑠 ) 𝑋 ( 𝑠 ) = 𝑄 ( 𝑠 ) 𝑃 ( 𝑠 ) = 𝑏 0 𝑠 𝑚 + 𝑏 1 𝑠 𝑚 − 1 + … + 𝑏 𝑚 𝑎 0 𝑠 𝑛 + 𝑎 1 𝑠 𝑛 − 1 + … + 𝑎 𝑛 -
传递函数通常为有理分式,即
都为实系数多项式。𝑃 、 𝑄 -
分母多项式
的次数𝑃 ( 𝑠 ) 称为系统的阶数,称系统为𝑛 阶系统。𝑛 -
特征方程
的根称为系统的极点,方程𝑃 ( 𝑠 ) = 0 的根称为系统的零点。𝑄 ( 𝑠 ) = 0
-
-
脉冲响应函数:系统对单位脉冲输入的响应输出。
- 单位脉冲输入的拉普拉斯变换为1,因此系统的传递函数等于系统脉冲响应函数的拉普拉斯变换。
- 通过用脉冲信号激励系统,测量系统的输出响应,就可以获得有关系统动态特性的全部信息。
自动控制系统
- 方框图的组成:方框单元、比较点、分支点。
-
:参考输入量 (Reference) .𝑅 ( 𝑠 ) -
:作用误差信号 (Error),𝐸 ( 𝑠 ) .𝐸 ( 𝑠 ) = 𝑅 ( 𝑠 ) − 𝐵 ( 𝑠 ) -
:控制系统/被控对象输出量 (Control output) .𝐶 ( 𝑠 ) -
:反馈信号 (Back) .𝐵 ( 𝑠 ) -
:前向传递函数,𝐺 ( 𝑠 ) .𝐺 ( 𝑠 ) = 𝐶 ( 𝑠 ) 𝐸 ( 𝑠 ) -
:反馈传递函数,𝐻 ( 𝑠 ) .𝐻 ( 𝑠 ) = 𝐵 ( 𝑠 ) 𝐶 ( 𝑠 ) -
开环传递函数 =
,描述反馈与误差(偏移量)之间的关系,刻画系统的固有特性。𝐵 ( 𝑠 ) 𝐸 ( 𝑠 ) = 𝐺 ( 𝑠 ) 𝐻 ( 𝑠 ) -
闭环传递函数 =
,直接描述系统输出与输入之间的关系,决定了系统的动态特性。𝐶 ( 𝑠 ) 𝑅 ( 𝑠 ) = 𝐺 ( 𝑠 ) 1 + 𝐺 ( 𝑠 ) 𝐻 ( 𝑠 ) -
可以导出
.𝐸 ( 𝑠 ) = 𝐶 ( 𝑠 ) 𝐺 ( 𝑠 ) = 𝑅 ( 𝑠 ) 1 + 𝐺 ( 𝑠 ) 𝐻 ( 𝑠 ) - 闭环特征方程直接源于开环传递函数。通过分析开环传递函数,也能研究闭环系统的动态特性。
-
可以导出
- 控制器:将被控对象输出量
与参考输入量𝑐 ( 𝑡 ) 进行比较,确定误差𝑟 ( 𝑡 ) ,并产生控制信号𝑒 ( 𝑡 ) .𝑢 ( 𝑡 ) - 执行器:动力装置,根据控制器的输出信号,产生被控对象的输入量。
- 被控对象:系统的主体,接受执行器的输入量,并产生输出量
.𝑐 ( 𝑡 ) - 传感器:测量被控对象的输出量,将其转化为特定信号反馈给控制器,用来与参考输入量进行比较。
PID控制器(比例-积分-微分)
控制器方程
传递函数
-
:控制器输出信号,𝑈 ( 𝑠 ) :作用误差信号,𝐸 ( 𝑠 ) :比例增益,𝐾 𝑝 :积分增益,𝐾 𝑖 :微分增益𝐾 𝑑 -
:积分时间常数,𝑇 𝑖 :微分时间常数,二者量纲均为时间。𝑇 𝑑 越大积分作用越弱,𝑇 𝑖 越大微分作用越强𝑇 𝑑
控制系统的分析
瞬态响应和稳态响应分析
控制系统的时间响应由瞬态响应和稳态响应两部分组成:
- 瞬态响应:系统从初始状态过渡到最终稳定状态的响应过程,会随时间衰减消失,由系统的初始状态决定,对应微分方程的通解。
- 稳态响应:系统经过足够长时间后的稳定输出状态,由系统的输入信号决定,对应微分方程的特解。
控制系统的动态特性
-
绝对稳定性:稳定指系统在任何有界输入下都能产生有界输出,不稳定则指系统在扰动下输出发散。
- 系统稳定的充要条件:系统的所有闭环极点均位于复平面的左半部。
- 系统在虚轴上的极点为单根,且其余极点均位于复平面的左半部时,系统为临界稳定,输出作等幅振荡,既不衰减,也不发散;若虚轴上的极点为重根,输出将发散。临界稳定无法物理实现,它不属于稳定。
- 相对稳定性:在系统已经绝对稳定的前提下,进一步衡量其稳定程度和抗干扰能力。
-
稳态误差:对于稳定系统,稳态输出与参考输出之间的差值,可以用来衡量系统精度。
- 对于单位反馈控制系统,反馈传递函数
,系统的期望输出等于输入信号𝐻 ( 𝑠 ) = 1 .𝑟 ( 𝑡 ) - 则误差
,稳态误差𝑒 ( 𝑡 ) = 𝑟 ( 𝑡 ) − 𝑐 ( 𝑡 ) .𝑒 𝑠 𝑠 = lim 𝑡 → ∞ 𝑒 ( 𝑡 ) = 𝑟 ( ∞ ) − 𝑐 ( ∞ )
- 对于单位反馈控制系统,反馈传递函数
高阶系统的瞬态响应
对系统的闭环传递函数的分子分母多项式进行因式分解:
若闭环极点互不相同,对
其中
- 如果闭环零点靠近某个闭环极点,则该极点的留数较小,相近的极点和零点相互抵消。
- 距离原点很远的极点,留数会很小,对应的瞬态响应项也会很小,且持续时间也很短。
- 具有很小留数的项可以忽略,则高阶系统可以用低阶系统近似表示。
- 将一对共轭复数极点写成一个
的二阶项,改写展开式,可知稳定高阶系统的响应曲线是一系列指数曲线与阻尼振荡曲线的叠加。𝑠
劳斯判据
劳斯判据给出系统稳定的充要条件,而不必实际求解特征方程的根。
设特征方程为
-
是系统稳定(即所有特征根都位于复平面左半部)的必要条件。𝑎 𝑖 > 0 -
的条件下,构造𝑎 𝑖 > 0 行的劳斯阵列:𝑛 + 1 其中前两行由特征方程的系数组成,且阵列满足递推关系:
𝑏 1 = − 1 𝑎 1 | 𝑎 0 𝑎 2 𝑎 1 𝑎 3 | , 𝑏 2 = − 1 𝑎 1 | 𝑎 0 𝑎 4 𝑎 1 𝑎 5 | , 𝑏 3 = − 1 𝑎 1 | 𝑎 0 𝑎 6 𝑎 1 𝑎 7 | , … 𝑐 1 = − 1 𝑏 1 | 𝑎 1 𝑏 1 𝑎 3 𝑏 2 | , 𝑐 2 = − 1 𝑏 1 | 𝑎 1 𝑏 1 𝑎 5 𝑏 3 | , 𝑐 3 = − 1 𝑏 1 | 𝑎 1 𝑏 1 𝑎 7 𝑏 4 | , …
则特征方程具有正实部的根的个数等于第一列元素符号变化的次数。
系统稳定的充要条件是:劳斯阵列第一列元素全为正数。 -
将
代入特征方程,并对以𝑠 = ˆ𝑠 − 𝜎 为变量的特征方程应用劳斯判据,可以得到位于 ˆ𝑠 右侧的根的个数,作为评估系统相对稳定性的一个指标。𝑠 = − 𝜎
积分和微分控制作用对系统性能的影响
如果被控对象的传递函数中不存在积分器
微分控制作用可以预测未来的误差,提前进行修正,从而提高系统的相对稳定性。它虽然不直接影响稳态误差,但增加了系统的阻尼,从而允许比较大的比例增益
单位反馈控制系统的稳态误差
终值定理给出了求系统稳态误差的简便方法:
设单位反馈控制系统具有开环传递函数:
开环传递函数位于原点处的极点个数
-
系统的型别
直接决定了系统跟踪典型输入信号(如阶跃、斜坡、加速度信号)的稳态精度。𝑁 - 要无差跟踪阶跃信号,系统至少是 1 型。定义静态位置误差常数
,则系统对单位阶跃输入的稳态误差𝐾 𝑝 = lim 𝑠 → 0 𝐺 ( 𝑠 ) = 𝐺 ( 0 ) 。对于 0 型系统,𝑒 𝑠 𝑠 = lim 𝑠 → 0 𝑠 1 + 𝐺 ( 𝑠 ) 1 𝑠 = 1 1 + 𝐾 𝑝 有限,稳态误差不为 0 ;对于 1 型及以上的系统,𝐾 𝑝 = 𝐾 ,稳态误差为 0 。𝐾 𝑝 = ∞ - 要无差跟踪斜坡信号,系统至少是 2 型。定义静态速度误差常数
,则系统对单位斜坡输入的稳态误差𝐾 𝑣 = lim 𝑠 → 0 𝑠 𝐺 ( 𝑠 ) 。对于 0 型系统,𝑒 𝑠 𝑠 = lim 𝑠 → 0 𝑠 1 + 𝐺 ( 𝑠 ) 1 𝑠 2 = 1 𝐾 𝑣 ,稳态误差为无穷大;对于 1 型系统,𝐾 𝑣 = 0 有限,稳态误差不为 0 ;对于 2 型及以上的系统,𝐾 𝑣 = 𝐾 ,稳态误差为 0 。𝐾 𝑣 = ∞ - 同理,可定义静态加速度误差常数等。静态误差常数是控制系统的品质指标,数值越大表示系统的稳态精度越高。
- 要无差跟踪阶跃信号,系统至少是 1 型。定义静态位置误差常数
- 原点处的极点本质上是一个积分器
,提高系统的型别,相当于使系统多一层累积,从而能够消除更高阶输入信号的稳态误差。1 𝑠 - 增加系统的型别能够提高系统的稳态精度,但也会降低系统的相对稳定性,导致更慢的瞬态响应和更大的超调。
根轨迹分析
根轨迹:系统某一参数(通常为开环增益
根轨迹的绘制规则
-
根轨迹起始于开环极点,终止于开环零点(或无穷远)。
- 设开环传递函数
,则特征方程𝐺 ( 𝑠 ) 𝐻 ( 𝑠 ) = 𝐾 𝑁 ( 𝑠 ) 𝐷 ( 𝑠 ) ,即1 + 𝐺 ( 𝑠 ) 𝐻 ( 𝑠 ) = 0 。𝐷 ( 𝑠 ) + 𝐾 𝑁 ( 𝑠 ) = 0
当 时,特征根为开环极点;当𝐾 = 0 时,特征根趋向于开环零点或无穷远。𝐾 → ∞
- 设开环传递函数
- 根轨迹关于实轴对称。这是实系数多项式方程的根的特性。
- 根轨迹的渐近线、分离点、与虚轴交点等可以分析确定。
频率响应分析
系统对正弦输入信号的稳态响应称为频率响应。频率响应法通过分析系统对不同频率正弦输入的响应,来评估系统的稳定性和性能。该方法可以利用对物理系统的测量数据,而不必推导出系统的数学模型。
可以证明,给稳定的线性时不变系统输入正弦信号
的模长表示系统对正弦输入信号幅值的增益,称为幅频特性;𝐻 ( 𝜔 ) 的幅角表示系统对正弦输入信号引入的相移,称为相频特性。频率响应函数是频率特性的复数表示。𝐻 ( 𝜔 ) - 系统的频率响应由系统的传递函数
唯一确定,而与系统的初始状态无关。𝐺 ( 𝑠 ) - 相比频率响应分析,时域瞬态响应分析更直观,但对高阶系统分析复杂,二者可以结合使用。