1. 从微分方程到传递函数为什么这是控制理论的“通用语言”搞自动控制无论是调一个电机转速还是稳定一架无人机的姿态我们最终面对的都是系统的数学模型。最开始大家习惯用微分方程来描述系统动态。比如一个简单的RC低通滤波电路其输入电压u(t)和输出电压y(t)的关系用微分方程写出来就是RC * dy(t)/dt y(t) u(t)。这个方程很直观物理意义明确但有个大问题求解和运算太麻烦。特别是对于复杂系统多个微分方程联立时域求解简直是噩梦更别提在此基础上分析稳定性、设计控制器了。传递函数的出现就是为了解决这个痛点。它的核心思想是将复杂的时域微积分运算转化为相对简单的复频域代数运算。具体操作是对微分方程两边进行拉普拉斯变换并假设系统的初始条件为零这是传递函数定义的前提。对于上面的RC电路方程两边取拉普拉斯变换利用导数性质L{dy/dt} sY(s) - y(0)并令初始条件y(0)0就得到RC * sY(s) Y(s) U(s)。接下来我们就能轻松地写出输出Y(s)与输入U(s)的比值这就是传递函数G(s)G(s) Y(s) / U(s) 1 / (RC * s 1)你看原来是一个微分方程现在变成了一个关于复变量s的有理分式。这个G(s)就是系统的传递函数。它完全表征了系统本身的动态特性与输入信号的具体形式无关。传递函数成了描述线性定常系统的“标准简历”我们通过分析这份“简历”比如看它的极点、零点就能预判系统对各种输入比如突然的指令、持续的干扰会作出何种反应。在实际工程中尤其是用到像MATLAB这样的工具时传递函数更是不可或缺。你想用bode函数画频率特性图用step函数看阶跃响应或者用rlocus画根轨迹输入参数基本都是传递函数模型。所以掌握传递函数就等于拿到了使用现代控制设计与分析工具的钥匙。2. 传递函数的“武功秘籍”形式、性质与内在限制理解了传递函数为什么重要我们还得摸清它的“脾气秉性”知道它能干什么不能干什么。这就像拿到一个功能强大的软件你得先读读说明书避免用错地方。2.1 标准形式与核心要素一个单输入单输出SISO系统的传递函数通常可以写成如下标准形式G(s) (b_m * s^m b_{m-1} * s^{m-1} ... b_1 * s b_0) / (a_n * s^n a_{n-1} * s^{n-1} ... a_1 * s a_0)其中n m这保证了系统的物理可实现性即系统不能预测未来。为了更清晰地分析我们常将其转化为两种更有用的形式零极点增益形式G(s) K * (s - z_1)(s - z_2)...(s - z_m) / [(s - p_1)(s - p_2)...(s - p_n)]K系统增益影响输出的整体幅度。z_i零点。使传递函数分子为零的s值。零点会影响系统响应的“形状”例如它可能削弱或抵消某些频率成分的影响。在复平面上用“○”表示。p_j极点。使传递函数分母为零的s值。极点决定了系统自由运动的模态如指数增长、衰减、振荡和稳定性。如果有一个极点实部大于零位于s平面右半部分系统就不稳定。在复平面上用“×”表示。这种形式是进行根轨迹分析、判断系统稳定性的直接依据。时间常数形式G(s) K * (τ_{z1}s 1)(τ_{z2}s 1)... / [(τ_{p1}s 1)(τ_{p2}s 1)...]τ时间常数具有时间量纲秒。它直观地反映了系统环节响应的快慢。例如1/(τs1)代表惯性环节其阶跃响应上升到稳态值的63.2%所需的时间就是τ。这种形式在分析系统频率特性伯德图时特别方便因为每个(τs1)项在伯德图上都有标准的渐近线画法。2.2 必须牢记的性质与前提传递函数不是万能的它的成立有几个关键前提忽略这些是初学者常踩的坑线性定常系统这是传递函数的“适用边界”。系统必须是线性的满足叠加性和齐次性且参数不随时间变化。如果你的系统有饱和、死区、齿轮间隙等非线性特性传递函数模型只在某个小范围内近似有效大范围分析必须用其他方法如描述函数法、相平面法。零初始条件定义传递函数时我们假设系统在t0时刻之前是“静止”的所有储能元件电容电压、电感电流、弹簧形变、质量块速度的初始状态为零。如果初始条件不为零拉氏变换方程中会多出一项初始条件项此时Y(s)/U(s)就不再是单纯的传递函数了。在实际应用中我们通常关注系统在平衡点附近受到扰动后的行为这个前提多数情况下是合理的。单输入单输出上述形式是针对SISO系统的。对于多输入多输出MIMO系统我们需要使用传递函数矩阵矩阵中的每个元素G_ij(s)表示第j个输入对第i个输出的传递函数。注意很多人在推导电路传递函数时容易忽略初始条件。例如推导一个二阶RLC电路的传递函数如果电容上有初始电压电感上有初始电流那么直接对微分方程进行拉氏变换得到的式子会比零初始条件下复杂得多。在绝大多数控制系统分析和设计场景中我们默认分析的是“零状态响应”因此使用零初始条件的传递函数是标准做法。2.3 从传递函数看系统性能拿到一个传递函数我们能一眼看出什么系统阶数分母s的最高次幂n就是系统的阶数。一阶系统响应无超调二阶系统可能振荡高阶系统行为更复杂但通常由主导极点决定。系统类型看分母中s的幂次。s^0即常数项对应0型系统s^1对应I型系统含有一个积分环节s^2对应II型系统。系统类型决定了它跟踪不同输入信号阶跃、斜坡、抛物线的稳态误差。稳定性如前所述只要有一个极点位于s平面右半部分实部0系统就不稳定。所有极点都在左半平面则系统稳定。极点在虚轴上实部0则临界稳定等幅振荡工程上通常也视为不稳定。动态响应速度极点离虚轴越远负实部绝对值越大对应的模态衰减越快对系统动态影响越小。最靠近虚轴的极点主导极点主要决定了系统的响应速度。3. 庖丁解牛七大典型环节的传递函数与物理意义复杂的控制系统往往可以分解为一些基本“积木块”的组合。这些“积木块”就是典型环节。熟练掌握它们就像熟练掌握了各种基础电路元件电阻、电容、电感一样任何复杂系统的传递函数你都能一眼看出它是由哪些环节以何种方式构成的。这是进行系统分析、设计和调试的基本功。3.1 比例环节最直接的放大传递函数G(s) K时域方程y(t) K * u(t)特点输出瞬间、成比例地复现输入。无相位滞后无失真。实例理想放大器、杠杆、减速器忽略惯性、电阻分压电路。在控制系统中控制器本身常常包含一个比例系数Kp这就是一个比例环节。伯德图幅频特性是一条高度为20lgKdB的水平线相频特性是0度线。3.2 积分环节累积的效应传递函数G(s) K / s时域方程y(t) K * ∫ u(τ) dτ(从0到t积分)特点输出是输入的积分累积。它能消除阶跃输入的稳态误差这是提高系统稳态精度的关键。但会带来-90度的相位滞后可能影响稳定性。实例水箱的水位h是进水流量q的积分h (1/A)∫ q dtA为底面积电容的电压u_C是电流i的积分u_C (1/C)∫ i dt运动物体的位移是速度的积分。伯德图幅频特性是一条斜率为-20dB/十倍频程的直线相频特性是恒为-90度。实操心得在PID控制器中I积分项就是用来消除静差的。但积分作用太强Ki太大容易导致系统响应初期超调大甚至引起振荡。实际调试时通常先调P再慢慢加入I并观察系统响应。3.3 微分环节预测的趋势传递函数G(s) K * s时域方程y(t) K * du(t)/dt特点输出与输入的变化率成正比。它能“预见”输入的变化趋势从而提供超前校正作用提高系统响应速度、改善稳定性。但对噪声极其敏感。实例理想的微分环节在物理世界中难以独立实现常与其他环节结合。测速发电机输出电势与电机转速/角速度成正比而转速是角位移的微分近似为一个微分环节。伯德图幅频特性是一条斜率为20dB/十倍频程的直线相频特性是恒为90度。实操心得纯微分环节会放大高频噪声在实际控制器如PID中使用的是实际微分环节G(s) Kd * s / (Td * s 1)。它只在一定的频率范围内具有微分效应高频段增益被限制避免了噪声放大问题。参数Td微分时间常数的选择至关重要。3.4 惯性环节缓慢的跟随传递函数G(s) K / (T s 1)时域方程T * dy(t)/dt y(t) K * u(t)特点当输入突变时输出不能立刻跟上而是按指数规律逐渐变化。时间常数T直接决定了响应速度T越大响应越慢。实例RC低通滤波电路TRC、一阶热系统如温度计测温、忽略电枢电感后的直流电机转速对电枢电压的响应。阶跃响应从0上升到稳态值的63.2%所需时间就是T经过3T~4T时间可认为达到稳态。伯德图幅频特性在转折频率ω1/T处开始以-20dB/十倍频程下降相频特性从0度趋向-90度在ω1/T时为-45度。3.5 振荡环节有阻尼的波动传递函数G(s) ω_n^2 / (s^2 2ζ ω_n s ω_n^2)特点这是二阶系统的标准形式。ω_n为无阻尼自然振荡频率ζ为阻尼比。阻尼比ζ决定了系统的振荡特性ζ 0无阻尼等幅振荡。0 ζ 1欠阻尼衰减振荡。这是最常见的工况。ζ 1临界阻尼无超调的最快响应。ζ 1过阻尼响应缓慢无振荡。实例质量-弹簧-阻尼系统、RLC串联电路、直流电机的完整模型考虑电枢电感。性能指标对于欠阻尼系统上升时间、峰值时间、超调量、调节时间等都与ω_n和ζ有直接的计算关系。例如超调量σ% e^{-πζ/√(1-ζ^2)} * 100%调节时间t_s ≈ 4/(ζ ω_n)按2%误差带。实操心得在运动控制中我们常常希望系统是轻微欠阻尼的如ζ0.7~0.8这样既有较快的响应速度超调量也在可接受范围内约1.5%~5%。通过调整控制器参数本质上就是在调整闭环系统的等效ω_n和ζ。3.6 一阶微分与二阶微分环节一阶微分G(s) τ s 1。它不是一个独立的物理环节通常与惯性环节成对出现构成超前校正网络用于提供相位超前改善系统动态性能。二阶微分G(s) τ^2 s^2 2ζ τ s 1。同样常与振荡环节组合用于更复杂的校正。3.7 延迟环节信息的迟到传递函数G(s) e^{-τ s}时域方程y(t) u(t - τ)特点输出是输入的纯时间延迟不改变信号形状只滞后时间τ。它对系统稳定性非常不利因为引入了额外的相位滞后。实例物料传输皮带、长管道流体输送、网络控制系统中的通信延迟、数字控制中的计算延迟。近似处理由于e^{-τ s}不是有理分式在经典控制理论中常采用帕德Pade近似将其转化为有理传递函数进行分析。例如一阶帕德近似为e^{-τ s} ≈ (1 - τ s / 2) / (1 τ s / 2)。伯德图幅频特性恒为0dB相频特性φ(ω) -τ ω弧度相位滞后随频率线性增加无边界。4. 实战推演从物理系统到传递函数的建模案例理论说得再多不如亲手算一遍。我们结合热词里的几个具体问题来看看如何从物理定律出发推导出传递函数。这个过程是控制工程师的基本功。4.1 案例一二阶低通滤波器传递函数推导这是一个非常经典的电路问题。假设我们有一个由电阻R、电感L、电容C串联组成的二阶低通滤波器输入是电压U_in(s)输出是电容两端的电压U_out(s)。列写元件方程电阻U_R(s) R * I(s)电感U_L(s) L * s * I(s)零初始条件下电容U_out(s) I(s) / (C * s)根据基尔霍夫电压定律KVLU_in(s) U_R(s) U_L(s) U_out(s)代入并整理U_in(s) R * I(s) L * s * I(s) I(s) / (C * s) I(s) * [R L s 1/(C s)]用输出表示电流从电容方程知I(s) C * s * U_out(s)代入上式U_in(s) C * s * U_out(s) * [R L s 1/(C s)] U_out(s) * [R C s L C s^2 1]得到传递函数G(s) U_out(s) / U_in(s) 1 / (L C s^2 R C s 1)化为标准振荡环节形式 令ω_n 1 / √(LC)ζ (R/2) * √(C/L)则G(s) ω_n^2 / (s^2 2ζ ω_n s ω_n^2)看这就是一个标准的振荡环节。ζ由R决定R越大阻尼越大滤波器的响应越平缓越不容易在截止频率附近产生谐振峰。4.2 案例二Boost变换器占空比到电感电流的传递函数推导CCM模式这是电力电子中的一个核心建模问题。推导需要用到状态空间平均法这里简述关键思路和结果让大家理解其复杂性。在连续导通模式CCM下Boost变换器在一个开关周期内有两个工作状态开关管导通和开关管关断。我们需要对这两个状态下的电路分别列写微分方程然后进行平均化处理。定义状态变量通常选择电感电流i_L和电容电压v_C即输出电压作为状态变量。分段线性建模当开关导通时持续时间d*Td为占空比输入电源V_in给电感充电电容向负载放电。当开关关断时持续时间(1-d)*T电感储能与输入电源一起向电容和负载释放能量。列写两个阶段的状态方程形式为dx/dt A1*x B1*u导通和dx/dt A2*x B2*u关断其中x[i_L; v_C]^TuV_in。状态空间平均得到平均状态方程dx/dt [d*A1 (1-d)*A2] * x [d*B1 (1-d)*B2] * u。注意此时方程中含有状态变量和占空比d的乘积项是非线性的。小信号线性化假设系统在某个稳态工作点(I_L, V_C, D)附近有微小扰动即令i_L I_L î_Lv_C V_C v̂_Cd D d̂。将这些代入平均方程忽略二阶及以上的扰动乘积项小信号假设可以将方程分解为直流稳态方程和交流小信号方程。求解传递函数对交流小信号方程进行拉普拉斯变换并消去中间变量最终可以得到占空比扰动d̂(s)到电感电流扰动î_L(s)的传递函数。这个传递函数通常形式比较复杂包含了右半平面零点RHPZ这是Boost等升压拓扑的固有特性会导致控制上的特殊挑战例如不能用单纯的积分环节否则可能不稳定。其最终形式大致为G_id(s) î_L(s) / d̂(s) [某个包含s的表达式] / [s^2 LC s(L/R) (1-D)^2]的具体变体。详细的推导涉及大量代数运算但核心思想就是上述的“分段线性 - 平均 - 线性化”。这个过程深刻体现了如何将一个非线性、时变的开关系统在特定假设下CCM小信号简化为一个可用于经典控制理论分析的线性定常传递函数模型。4.3 案例三利用MATLAB计算复杂电路传递函数对于复杂电路手动推导传递函数非常繁琐且容易出错。此时可以借助MATLAB的符号数学工具箱Symbolic Math Toolbox来辅助完成。假设我们有一个稍复杂的无源网络包含多个R、L、C。我们可以列写s域阻抗方程将电阻R、电感sL、电容1/(sC)代入电路。使用基尔霍夫定律或节点电压法、网孔电流法列方程。在MATLAB中定义符号变量syms s R1 R2 L C Vin Vout % 定义符号变量列写方程并求解% 假设使用节点电压法设某个节点电压为Vx eq1 (Vin - Vx)/R1 Vx/(1/(s*C)) (Vx - Vout)/s*L); eq2 (Vx - Vout)/(s*L) Vout/R2; % 求解方程组得到Vout关于Vin的表达式 sol solve([eq1, eq2], [Vx, Vout]); G simplify(sol.Vout / Vin) % 化简得到传递函数G(s)结果MATLAB会输出化简后的传递函数表达式G(s)。你还可以使用pretty(G)命令让显示更直观。这种方法特别适合验证手算结果或者处理那些结构复杂、难以一眼看出的电路。得到符号传递函数后你可以代入具体的元件参数值用tf函数生成数值传递函数模型进而进行时域仿真或频域分析。5. 典型环节的组合与系统辨识如何看懂一个复杂传递函数在实际工程中我们拿到一个复杂系统的传递函数或者从频域测试数据中拟合出一个传递函数模型如何快速理解它这就需要我们具备“拆解”的能力。5.1 串联、并联与反馈连接复杂系统通常由典型环节通过三种基本方式连接而成串联总传递函数等于各环节传递函数之积。G_total(s) G1(s) * G2(s) * ...物理意义信号依次通过各个环节。伯德图总幅频特性 各环节幅频特性相加dB值总相频特性 各环节相频特性相加。并联总传递函数等于各环节传递函数之和。G_total(s) G1(s) G2(s) ...物理意义多个环节对同一个输入信号进行处理其输出叠加在一起。注意在复频域直接相加在伯德图上不能简单叠加需要计算。反馈连接这是控制系统中最核心的结构。基本负反馈结构的总传递函数闭环传递函数为Φ(s) G(s) / [1 G(s)H(s)]其中G(s)是前向通道传递函数H(s)是反馈通道传递函数。如果H(s)1则为单位负反馈。5.2 系统辨识从伯德图反推传递函数很多时候我们通过实验测得了系统的频率响应数据伯德图需要从中估计出传递函数。这个过程叫系统辨识。观察幅频特性渐近线从低频段开始看渐近线的斜率。斜率为0 dB/dec - 比例或振荡环节的低频段。斜率变为-20 dB/dec - 遇到了一个惯性环节或积分环节的转折频率。斜率变为-40 dB/dec - 可能遇到了两个惯性环节转折频率很接近或者一个振荡环节。斜率变为20 dB/dec - 遇到了一阶微分环节。斜率在某个频点急剧下降 - 可能遇到了高阶环节或谐振峰。确定环节类型和参数积分环节低频段斜率为-20 dB/dec的直线延长线与0dB线的交点频率即为KωK。惯性环节在转折频率ω1/T处斜率增加-20 dB/dec。从幅频曲线可以找到该转折点。振荡环节在自然频率ω_n附近可能出现谐振峰。通过峰值高度可以估算阻尼比ζM_r 1/(2ζ√(1-ζ^2))对于ζ0.707。一阶微分环节在转折频率ω1/τ处斜率增加20 dB/dec。结合相频特性验证惯性环节在转折频率处相位滞后约45度高频段滞后90度。积分环节始终滞后90度。一阶微分环节在转折频率处相位超前约45度高频段超前90度。振荡环节在ω_n附近相位变化剧烈从0度快速下降到-180度。例如你测得的伯德图显示低频段斜率为0在ω1处斜率变为-20 dB/dec在ω2处斜率变为-40 dB/dec在ω3处斜率又变回-20 dB/dec。对应的相频特性在低频接近0度随频率增加而滞后在ω2附近滞后达到最大可能超过-90度。那么一个合理的传递函数猜测可能是G(s) K / [(T1 s 1)( (s/ω_n)^2 2ζ(s/ω_n) 1 )]其中1/T1 ≈ ω1ω_n ≈ ω2而ω3可能是振荡环节的第二个转折点或另一个惯性环节。这个过程需要反复尝试和验证可以借助MATLAB的tfest或系统辨识工具箱来辅助完成曲线拟合。6. 传递函数在控制系统设计与分析中的核心应用掌握了传递函数的表示和典型环节我们最终要把它用起来。传递函数是连接系统建模、分析和设计的桥梁。6.1 稳定性分析劳斯判据与奈奎斯特判据直接法劳斯判据给定闭环系统的特征方程即传递函数分母多项式等于零列出劳斯表。第一列符号改变的次数就是系统在右半平面的极点个数。这是基于传递函数系数进行代数判断的方法非常直接。频域法奈奎斯特判据根据开环传递函数G(s)H(s)的奈奎斯特图环绕(-1, j0)点的圈数等于开环右极点数与闭环右极点数的差值。这不仅能判断稳定性还能定性地看出稳定裕度相位裕度、幅值裕度。6.2 稳态误差计算利用终值定理和系统类型可以直接从传递函数计算系统对典型输入阶跃、斜坡、抛物线的稳态误差。e_ss lim_{s-0} s * R(s) / [1 G_open(s)]其中G_open(s)是开环传递函数。0型系统对阶跃输入有静差对斜坡输入误差无穷大。I型系统对阶跃输入无静差对斜坡输入有恒值静差。II型系统对阶跃和斜坡输入均无静差对抛物线输入有恒值静差。6.3 根轨迹设计法根轨迹是开环增益K从0变化到∞时闭环系统极点在s平面上的运动轨迹。绘制根轨迹的规则完全基于开环传递函数的零极点分布。通过调整增益K或在开环中增加零极点即设计校正环节我们可以让闭环极点移动到期望的位置从而获得理想的动态性能如指定的阻尼比ζ和自然频率ω_n。6.4 频域校正设计法伯德图法这是工程上最常用的方法之一。我们画出未校正系统的开环伯德图根据性能指标稳态误差、相位裕度、截止频率要求设计串联校正装置如超前校正、滞后校正、滞后-超前校正。超前校正本质是一阶微分环节为主提供相位超前提高系统稳定性和响应速度。通过在伯德图上抬高中频段增加相位裕度来实现。滞后校正本质是惯性环节为主提高低频增益减小稳态误差但同时会降低中频段可能使响应变慢。通常利用其高频衰减特性压低高频噪声。设计过程就是基于开环传递函数的频率特性反复试探和调整校正环节的参数直到满足所有指标。6.5 在MATLAB/Simulink中的实现现代控制设计几乎离不开这些工具。传递函数是这些工具的“母语”。定义模型sys tf(num, den)或sys zpk(z, p, k)。时域分析step(sys),impulse(sys),lsim(sys, u, t)。频域分析bode(sys),nyquist(sys),margin(sys)。根轨迹rlocus(sys)。连接系统series,parallel,feedback函数。控制器设计pidtune,sisotool等工具箱其内部运算的核心对象就是传递函数。从我个人的经验来看很多初学者在学自动控制原理时容易陷入公式推导的细节而忽略了传递函数作为一个“工具”的直观性和实用性。我的建议是在学习每个典型环节时一定要在MATLAB里把它画出来看看它的阶跃响应长什么样看看它的伯德图有什么特征。把抽象的数学表达式和直观的图形响应对应起来这种形象的理解比死记硬背公式要牢固得多。当你拿到一个实际系统无论是电路、机械还是热力学系统第一反应就应该是“它的主导环节是什么是惯性为主还是振荡为主时间常数大概多大” 有了这种直觉无论是分析问题还是设计控制器方向都会清晰很多。 
长河君 · 2026/8/2 4:20:05 · 阅读 1.2k 从微分方程到传递函数：控制理论核心建模与典型环节解析
1. 从微分方程到传递函数为什么这是控制理论的“通用语言”搞自动控制无论是调一个电机转速还是稳定一架无人机的姿态我们最终面对的都是系统的数学模型。最开始大家习惯用微分方程来描述系统动态。比如一个简单的RC低通滤波电路其输入电压u(t)和输出电压y(t)的关系用微分方程写出来就是RC * dy(t)/dt y(t) u(t)。这个方程很直观物理意义明确但有个大问题求解和运算太麻烦。特别是对于复杂系统多个微分方程联立时域求解简直是噩梦更别提在此基础上分析稳定性、设计控制器了。传递函数的出现就是为了解决这个痛点。它的核心思想是将复杂的时域微积分运算转化为相对简单的复频域代数运算。具体操作是对微分方程两边进行拉普拉斯变换并假设系统的初始条件为零这是传递函数定义的前提。对于上面的RC电路方程两边取拉普拉斯变换利用导数性质L{dy/dt} sY(s) - y(0)并令初始条件y(0)0就得到RC * sY(s) Y(s) U(s)。接下来我们就能轻松地写出输出Y(s)与输入U(s)的比值这就是传递函数G(s)G(s) Y(s) / U(s) 1 / (RC * s 1)你看原来是一个微分方程现在变成了一个关于复变量s的有理分式。这个G(s)就是系统的传递函数。它完全表征了系统本身的动态特性与输入信号的具体形式无关。传递函数成了描述线性定常系统的“标准简历”我们通过分析这份“简历”比如看它的极点、零点就能预判系统对各种输入比如突然的指令、持续的干扰会作出何种反应。在实际工程中尤其是用到像MATLAB这样的工具时传递函数更是不可或缺。你想用bode函数画频率特性图用step函数看阶跃响应或者用rlocus画根轨迹输入参数基本都是传递函数模型。所以掌握传递函数就等于拿到了使用现代控制设计与分析工具的钥匙。2. 传递函数的“武功秘籍”形式、性质与内在限制理解了传递函数为什么重要我们还得摸清它的“脾气秉性”知道它能干什么不能干什么。这就像拿到一个功能强大的软件你得先读读说明书避免用错地方。2.1 标准形式与核心要素一个单输入单输出SISO系统的传递函数通常可以写成如下标准形式G(s) (b_m * s^m b_{m-1} * s^{m-1} ... b_1 * s b_0) / (a_n * s^n a_{n-1} * s^{n-1} ... a_1 * s a_0)其中n m这保证了系统的物理可实现性即系统不能预测未来。为了更清晰地分析我们常将其转化为两种更有用的形式零极点增益形式G(s) K * (s - z_1)(s - z_2)...(s - z_m) / [(s - p_1)(s - p_2)...(s - p_n)]K系统增益影响输出的整体幅度。z_i零点。使传递函数分子为零的s值。零点会影响系统响应的“形状”例如它可能削弱或抵消某些频率成分的影响。在复平面上用“○”表示。p_j极点。使传递函数分母为零的s值。极点决定了系统自由运动的模态如指数增长、衰减、振荡和稳定性。如果有一个极点实部大于零位于s平面右半部分系统就不稳定。在复平面上用“×”表示。这种形式是进行根轨迹分析、判断系统稳定性的直接依据。时间常数形式G(s) K * (τ_{z1}s 1)(τ_{z2}s 1)... / [(τ_{p1}s 1)(τ_{p2}s 1)...]τ时间常数具有时间量纲秒。它直观地反映了系统环节响应的快慢。例如1/(τs1)代表惯性环节其阶跃响应上升到稳态值的63.2%所需的时间就是τ。这种形式在分析系统频率特性伯德图时特别方便因为每个(τs1)项在伯德图上都有标准的渐近线画法。2.2 必须牢记的性质与前提传递函数不是万能的它的成立有几个关键前提忽略这些是初学者常踩的坑线性定常系统这是传递函数的“适用边界”。系统必须是线性的满足叠加性和齐次性且参数不随时间变化。如果你的系统有饱和、死区、齿轮间隙等非线性特性传递函数模型只在某个小范围内近似有效大范围分析必须用其他方法如描述函数法、相平面法。零初始条件定义传递函数时我们假设系统在t0时刻之前是“静止”的所有储能元件电容电压、电感电流、弹簧形变、质量块速度的初始状态为零。如果初始条件不为零拉氏变换方程中会多出一项初始条件项此时Y(s)/U(s)就不再是单纯的传递函数了。在实际应用中我们通常关注系统在平衡点附近受到扰动后的行为这个前提多数情况下是合理的。单输入单输出上述形式是针对SISO系统的。对于多输入多输出MIMO系统我们需要使用传递函数矩阵矩阵中的每个元素G_ij(s)表示第j个输入对第i个输出的传递函数。注意很多人在推导电路传递函数时容易忽略初始条件。例如推导一个二阶RLC电路的传递函数如果电容上有初始电压电感上有初始电流那么直接对微分方程进行拉氏变换得到的式子会比零初始条件下复杂得多。在绝大多数控制系统分析和设计场景中我们默认分析的是“零状态响应”因此使用零初始条件的传递函数是标准做法。2.3 从传递函数看系统性能拿到一个传递函数我们能一眼看出什么系统阶数分母s的最高次幂n就是系统的阶数。一阶系统响应无超调二阶系统可能振荡高阶系统行为更复杂但通常由主导极点决定。系统类型看分母中s的幂次。s^0即常数项对应0型系统s^1对应I型系统含有一个积分环节s^2对应II型系统。系统类型决定了它跟踪不同输入信号阶跃、斜坡、抛物线的稳态误差。稳定性如前所述只要有一个极点位于s平面右半部分实部0系统就不稳定。所有极点都在左半平面则系统稳定。极点在虚轴上实部0则临界稳定等幅振荡工程上通常也视为不稳定。动态响应速度极点离虚轴越远负实部绝对值越大对应的模态衰减越快对系统动态影响越小。最靠近虚轴的极点主导极点主要决定了系统的响应速度。3. 庖丁解牛七大典型环节的传递函数与物理意义复杂的控制系统往往可以分解为一些基本“积木块”的组合。这些“积木块”就是典型环节。熟练掌握它们就像熟练掌握了各种基础电路元件电阻、电容、电感一样任何复杂系统的传递函数你都能一眼看出它是由哪些环节以何种方式构成的。这是进行系统分析、设计和调试的基本功。3.1 比例环节最直接的放大传递函数G(s) K时域方程y(t) K * u(t)特点输出瞬间、成比例地复现输入。无相位滞后无失真。实例理想放大器、杠杆、减速器忽略惯性、电阻分压电路。在控制系统中控制器本身常常包含一个比例系数Kp这就是一个比例环节。伯德图幅频特性是一条高度为20lgKdB的水平线相频特性是0度线。3.2 积分环节累积的效应传递函数G(s) K / s时域方程y(t) K * ∫ u(τ) dτ(从0到t积分)特点输出是输入的积分累积。它能消除阶跃输入的稳态误差这是提高系统稳态精度的关键。但会带来-90度的相位滞后可能影响稳定性。实例水箱的水位h是进水流量q的积分h (1/A)∫ q dtA为底面积电容的电压u_C是电流i的积分u_C (1/C)∫ i dt运动物体的位移是速度的积分。伯德图幅频特性是一条斜率为-20dB/十倍频程的直线相频特性是恒为-90度。实操心得在PID控制器中I积分项就是用来消除静差的。但积分作用太强Ki太大容易导致系统响应初期超调大甚至引起振荡。实际调试时通常先调P再慢慢加入I并观察系统响应。3.3 微分环节预测的趋势传递函数G(s) K * s时域方程y(t) K * du(t)/dt特点输出与输入的变化率成正比。它能“预见”输入的变化趋势从而提供超前校正作用提高系统响应速度、改善稳定性。但对噪声极其敏感。实例理想的微分环节在物理世界中难以独立实现常与其他环节结合。测速发电机输出电势与电机转速/角速度成正比而转速是角位移的微分近似为一个微分环节。伯德图幅频特性是一条斜率为20dB/十倍频程的直线相频特性是恒为90度。实操心得纯微分环节会放大高频噪声在实际控制器如PID中使用的是实际微分环节G(s) Kd * s / (Td * s 1)。它只在一定的频率范围内具有微分效应高频段增益被限制避免了噪声放大问题。参数Td微分时间常数的选择至关重要。3.4 惯性环节缓慢的跟随传递函数G(s) K / (T s 1)时域方程T * dy(t)/dt y(t) K * u(t)特点当输入突变时输出不能立刻跟上而是按指数规律逐渐变化。时间常数T直接决定了响应速度T越大响应越慢。实例RC低通滤波电路TRC、一阶热系统如温度计测温、忽略电枢电感后的直流电机转速对电枢电压的响应。阶跃响应从0上升到稳态值的63.2%所需时间就是T经过3T~4T时间可认为达到稳态。伯德图幅频特性在转折频率ω1/T处开始以-20dB/十倍频程下降相频特性从0度趋向-90度在ω1/T时为-45度。3.5 振荡环节有阻尼的波动传递函数G(s) ω_n^2 / (s^2 2ζ ω_n s ω_n^2)特点这是二阶系统的标准形式。ω_n为无阻尼自然振荡频率ζ为阻尼比。阻尼比ζ决定了系统的振荡特性ζ 0无阻尼等幅振荡。0 ζ 1欠阻尼衰减振荡。这是最常见的工况。ζ 1临界阻尼无超调的最快响应。ζ 1过阻尼响应缓慢无振荡。实例质量-弹簧-阻尼系统、RLC串联电路、直流电机的完整模型考虑电枢电感。性能指标对于欠阻尼系统上升时间、峰值时间、超调量、调节时间等都与ω_n和ζ有直接的计算关系。例如超调量σ% e^{-πζ/√(1-ζ^2)} * 100%调节时间t_s ≈ 4/(ζ ω_n)按2%误差带。实操心得在运动控制中我们常常希望系统是轻微欠阻尼的如ζ0.7~0.8这样既有较快的响应速度超调量也在可接受范围内约1.5%~5%。通过调整控制器参数本质上就是在调整闭环系统的等效ω_n和ζ。3.6 一阶微分与二阶微分环节一阶微分G(s) τ s 1。它不是一个独立的物理环节通常与惯性环节成对出现构成超前校正网络用于提供相位超前改善系统动态性能。二阶微分G(s) τ^2 s^2 2ζ τ s 1。同样常与振荡环节组合用于更复杂的校正。3.7 延迟环节信息的迟到传递函数G(s) e^{-τ s}时域方程y(t) u(t - τ)特点输出是输入的纯时间延迟不改变信号形状只滞后时间τ。它对系统稳定性非常不利因为引入了额外的相位滞后。实例物料传输皮带、长管道流体输送、网络控制系统中的通信延迟、数字控制中的计算延迟。近似处理由于e^{-τ s}不是有理分式在经典控制理论中常采用帕德Pade近似将其转化为有理传递函数进行分析。例如一阶帕德近似为e^{-τ s} ≈ (1 - τ s / 2) / (1 τ s / 2)。伯德图幅频特性恒为0dB相频特性φ(ω) -τ ω弧度相位滞后随频率线性增加无边界。4. 实战推演从物理系统到传递函数的建模案例理论说得再多不如亲手算一遍。我们结合热词里的几个具体问题来看看如何从物理定律出发推导出传递函数。这个过程是控制工程师的基本功。4.1 案例一二阶低通滤波器传递函数推导这是一个非常经典的电路问题。假设我们有一个由电阻R、电感L、电容C串联组成的二阶低通滤波器输入是电压U_in(s)输出是电容两端的电压U_out(s)。列写元件方程电阻U_R(s) R * I(s)电感U_L(s) L * s * I(s)零初始条件下电容U_out(s) I(s) / (C * s)根据基尔霍夫电压定律KVLU_in(s) U_R(s) U_L(s) U_out(s)代入并整理U_in(s) R * I(s) L * s * I(s) I(s) / (C * s) I(s) * [R L s 1/(C s)]用输出表示电流从电容方程知I(s) C * s * U_out(s)代入上式U_in(s) C * s * U_out(s) * [R L s 1/(C s)] U_out(s) * [R C s L C s^2 1]得到传递函数G(s) U_out(s) / U_in(s) 1 / (L C s^2 R C s 1)化为标准振荡环节形式 令ω_n 1 / √(LC)ζ (R/2) * √(C/L)则G(s) ω_n^2 / (s^2 2ζ ω_n s ω_n^2)看这就是一个标准的振荡环节。ζ由R决定R越大阻尼越大滤波器的响应越平缓越不容易在截止频率附近产生谐振峰。4.2 案例二Boost变换器占空比到电感电流的传递函数推导CCM模式这是电力电子中的一个核心建模问题。推导需要用到状态空间平均法这里简述关键思路和结果让大家理解其复杂性。在连续导通模式CCM下Boost变换器在一个开关周期内有两个工作状态开关管导通和开关管关断。我们需要对这两个状态下的电路分别列写微分方程然后进行平均化处理。定义状态变量通常选择电感电流i_L和电容电压v_C即输出电压作为状态变量。分段线性建模当开关导通时持续时间d*Td为占空比输入电源V_in给电感充电电容向负载放电。当开关关断时持续时间(1-d)*T电感储能与输入电源一起向电容和负载释放能量。列写两个阶段的状态方程形式为dx/dt A1*x B1*u导通和dx/dt A2*x B2*u关断其中x[i_L; v_C]^TuV_in。状态空间平均得到平均状态方程dx/dt [d*A1 (1-d)*A2] * x [d*B1 (1-d)*B2] * u。注意此时方程中含有状态变量和占空比d的乘积项是非线性的。小信号线性化假设系统在某个稳态工作点(I_L, V_C, D)附近有微小扰动即令i_L I_L î_Lv_C V_C v̂_Cd D d̂。将这些代入平均方程忽略二阶及以上的扰动乘积项小信号假设可以将方程分解为直流稳态方程和交流小信号方程。求解传递函数对交流小信号方程进行拉普拉斯变换并消去中间变量最终可以得到占空比扰动d̂(s)到电感电流扰动î_L(s)的传递函数。这个传递函数通常形式比较复杂包含了右半平面零点RHPZ这是Boost等升压拓扑的固有特性会导致控制上的特殊挑战例如不能用单纯的积分环节否则可能不稳定。其最终形式大致为G_id(s) î_L(s) / d̂(s) [某个包含s的表达式] / [s^2 LC s(L/R) (1-D)^2]的具体变体。详细的推导涉及大量代数运算但核心思想就是上述的“分段线性 - 平均 - 线性化”。这个过程深刻体现了如何将一个非线性、时变的开关系统在特定假设下CCM小信号简化为一个可用于经典控制理论分析的线性定常传递函数模型。4.3 案例三利用MATLAB计算复杂电路传递函数对于复杂电路手动推导传递函数非常繁琐且容易出错。此时可以借助MATLAB的符号数学工具箱Symbolic Math Toolbox来辅助完成。假设我们有一个稍复杂的无源网络包含多个R、L、C。我们可以列写s域阻抗方程将电阻R、电感sL、电容1/(sC)代入电路。使用基尔霍夫定律或节点电压法、网孔电流法列方程。在MATLAB中定义符号变量syms s R1 R2 L C Vin Vout % 定义符号变量列写方程并求解% 假设使用节点电压法设某个节点电压为Vx eq1 (Vin - Vx)/R1 Vx/(1/(s*C)) (Vx - Vout)/s*L); eq2 (Vx - Vout)/(s*L) Vout/R2; % 求解方程组得到Vout关于Vin的表达式 sol solve([eq1, eq2], [Vx, Vout]); G simplify(sol.Vout / Vin) % 化简得到传递函数G(s)结果MATLAB会输出化简后的传递函数表达式G(s)。你还可以使用pretty(G)命令让显示更直观。这种方法特别适合验证手算结果或者处理那些结构复杂、难以一眼看出的电路。得到符号传递函数后你可以代入具体的元件参数值用tf函数生成数值传递函数模型进而进行时域仿真或频域分析。5. 典型环节的组合与系统辨识如何看懂一个复杂传递函数在实际工程中我们拿到一个复杂系统的传递函数或者从频域测试数据中拟合出一个传递函数模型如何快速理解它这就需要我们具备“拆解”的能力。5.1 串联、并联与反馈连接复杂系统通常由典型环节通过三种基本方式连接而成串联总传递函数等于各环节传递函数之积。G_total(s) G1(s) * G2(s) * ...物理意义信号依次通过各个环节。伯德图总幅频特性 各环节幅频特性相加dB值总相频特性 各环节相频特性相加。并联总传递函数等于各环节传递函数之和。G_total(s) G1(s) G2(s) ...物理意义多个环节对同一个输入信号进行处理其输出叠加在一起。注意在复频域直接相加在伯德图上不能简单叠加需要计算。反馈连接这是控制系统中最核心的结构。基本负反馈结构的总传递函数闭环传递函数为Φ(s) G(s) / [1 G(s)H(s)]其中G(s)是前向通道传递函数H(s)是反馈通道传递函数。如果H(s)1则为单位负反馈。5.2 系统辨识从伯德图反推传递函数很多时候我们通过实验测得了系统的频率响应数据伯德图需要从中估计出传递函数。这个过程叫系统辨识。观察幅频特性渐近线从低频段开始看渐近线的斜率。斜率为0 dB/dec - 比例或振荡环节的低频段。斜率变为-20 dB/dec - 遇到了一个惯性环节或积分环节的转折频率。斜率变为-40 dB/dec - 可能遇到了两个惯性环节转折频率很接近或者一个振荡环节。斜率变为20 dB/dec - 遇到了一阶微分环节。斜率在某个频点急剧下降 - 可能遇到了高阶环节或谐振峰。确定环节类型和参数积分环节低频段斜率为-20 dB/dec的直线延长线与0dB线的交点频率即为KωK。惯性环节在转折频率ω1/T处斜率增加-20 dB/dec。从幅频曲线可以找到该转折点。振荡环节在自然频率ω_n附近可能出现谐振峰。通过峰值高度可以估算阻尼比ζM_r 1/(2ζ√(1-ζ^2))对于ζ0.707。一阶微分环节在转折频率ω1/τ处斜率增加20 dB/dec。结合相频特性验证惯性环节在转折频率处相位滞后约45度高频段滞后90度。积分环节始终滞后90度。一阶微分环节在转折频率处相位超前约45度高频段超前90度。振荡环节在ω_n附近相位变化剧烈从0度快速下降到-180度。例如你测得的伯德图显示低频段斜率为0在ω1处斜率变为-20 dB/dec在ω2处斜率变为-40 dB/dec在ω3处斜率又变回-20 dB/dec。对应的相频特性在低频接近0度随频率增加而滞后在ω2附近滞后达到最大可能超过-90度。那么一个合理的传递函数猜测可能是G(s) K / [(T1 s 1)( (s/ω_n)^2 2ζ(s/ω_n) 1 )]其中1/T1 ≈ ω1ω_n ≈ ω2而ω3可能是振荡环节的第二个转折点或另一个惯性环节。这个过程需要反复尝试和验证可以借助MATLAB的tfest或系统辨识工具箱来辅助完成曲线拟合。6. 传递函数在控制系统设计与分析中的核心应用掌握了传递函数的表示和典型环节我们最终要把它用起来。传递函数是连接系统建模、分析和设计的桥梁。6.1 稳定性分析劳斯判据与奈奎斯特判据直接法劳斯判据给定闭环系统的特征方程即传递函数分母多项式等于零列出劳斯表。第一列符号改变的次数就是系统在右半平面的极点个数。这是基于传递函数系数进行代数判断的方法非常直接。频域法奈奎斯特判据根据开环传递函数G(s)H(s)的奈奎斯特图环绕(-1, j0)点的圈数等于开环右极点数与闭环右极点数的差值。这不仅能判断稳定性还能定性地看出稳定裕度相位裕度、幅值裕度。6.2 稳态误差计算利用终值定理和系统类型可以直接从传递函数计算系统对典型输入阶跃、斜坡、抛物线的稳态误差。e_ss lim_{s-0} s * R(s) / [1 G_open(s)]其中G_open(s)是开环传递函数。0型系统对阶跃输入有静差对斜坡输入误差无穷大。I型系统对阶跃输入无静差对斜坡输入有恒值静差。II型系统对阶跃和斜坡输入均无静差对抛物线输入有恒值静差。6.3 根轨迹设计法根轨迹是开环增益K从0变化到∞时闭环系统极点在s平面上的运动轨迹。绘制根轨迹的规则完全基于开环传递函数的零极点分布。通过调整增益K或在开环中增加零极点即设计校正环节我们可以让闭环极点移动到期望的位置从而获得理想的动态性能如指定的阻尼比ζ和自然频率ω_n。6.4 频域校正设计法伯德图法这是工程上最常用的方法之一。我们画出未校正系统的开环伯德图根据性能指标稳态误差、相位裕度、截止频率要求设计串联校正装置如超前校正、滞后校正、滞后-超前校正。超前校正本质是一阶微分环节为主提供相位超前提高系统稳定性和响应速度。通过在伯德图上抬高中频段增加相位裕度来实现。滞后校正本质是惯性环节为主提高低频增益减小稳态误差但同时会降低中频段可能使响应变慢。通常利用其高频衰减特性压低高频噪声。设计过程就是基于开环传递函数的频率特性反复试探和调整校正环节的参数直到满足所有指标。6.5 在MATLAB/Simulink中的实现现代控制设计几乎离不开这些工具。传递函数是这些工具的“母语”。定义模型sys tf(num, den)或sys zpk(z, p, k)。时域分析step(sys),impulse(sys),lsim(sys, u, t)。频域分析bode(sys),nyquist(sys),margin(sys)。根轨迹rlocus(sys)。连接系统series,parallel,feedback函数。控制器设计pidtune,sisotool等工具箱其内部运算的核心对象就是传递函数。从我个人的经验来看很多初学者在学自动控制原理时容易陷入公式推导的细节而忽略了传递函数作为一个“工具”的直观性和实用性。我的建议是在学习每个典型环节时一定要在MATLAB里把它画出来看看它的阶跃响应长什么样看看它的伯德图有什么特征。把抽象的数学表达式和直观的图形响应对应起来这种形象的理解比死记硬背公式要牢固得多。当你拿到一个实际系统无论是电路、机械还是热力学系统第一反应就应该是“它的主导环节是什么是惯性为主还是振荡为主时间常数大概多大” 有了这种直觉无论是分析问题还是设计控制器方向都会清晰很多。