3. 船舶操纵运动响应模型

在船舶航向控制系统的设计与仿真中，我们并不需要一开始就建立包含六自由度、附加质量、粘性水动力导数的完整操纵运动方程。对于绝大多数航向保持与航向改变问题，工程上更常用的是以舵角为输入、艏向角为输出的低阶传递函数模型——这就是Nomoto模型。本节围绕一阶与二阶Nomoto模型展开，重点讲清它们的推导来源、KT参数的物理意义，以及在实际仿真中如何选型与取值。

从操纵运动方程到Nomoto模型

先约定坐标系与符号：采用惯性坐标系下的艏向角 ψ（以正北或初始航向为基准，右转为正），舵角 δ 以右舵为正。船舶平面操纵的线性化方程为：

m(u̇ - vr) = X
m(v̇ + ur) = Y
I_z ṙ = N

其中 v 为横漂速度，r 为艏摇角速度，Y、N 为横向水动力与艏摇力矩。消去横漂速度 v 后，可得到 r 与 δ 之间的二阶微分方程。Nomoto 指出，对于航向稳定性较好的船舶，该方程可近似降阶为一阶形式：

T·ṙ + r = K·δ        （一阶Nomoto，r 为艏摇角速度）
T·ψ̈ + ψ̇ = K·δ        （等价地，以艏向角 ψ 为输出）

对应的传递函数为 G(s) = K / (T·s + 1)。当船舶航向稳定性较弱、或需要更精确描述高频段相位特性时，则保留二阶项：

G(s) = K / [(T₁·s + 1)(T₂·s + 1)]

这就是二阶Nomoto模型。它相当于在一阶模型基础上增加了一个时间常数 T₂，用来刻画舵效响应的初始延迟与高频衰减。

KT参数的物理意义

K 称为回转性指数（增益），单位通常为 1/s。它表示单位舵角下船舶最终能达到的稳态艏摇角速度：当 δ 为常值时，稳态 r_ss = K·δ。K 越大，船舶对舵的响应越"灵敏"，回转圈越小。T 称为追随性指数（时间常数），单位为秒，反映船舶达到稳态回转的滞后程度：T 越大，响应越迟缓，操舵后需要更长时间才能建立稳定回转。二阶模型中的 T₁、T₂ 分别对应主导时间常数与次要时间常数，一般 T₁ > T₂，且 T₁+T₂ ≈ T（一阶等效值）。

需要强调，K 与 T 都不是常数，而是随航速变化的。工程上常用无量纲形式 K' = K·(L/U)、T' = T·(U/L)（L 为船长，U 为航速），在傅汝德数相近时近似为常数。仿真时若航速改变，必须重新换算 K、T，否则模型失真。

参数取值与仿真实现

在MATLAB中，一阶模型可直接用 tf(K,[T 1]) 建立，二阶模型用 tf(K,conv([T1 1],[T2 1]))。以某中型货船为例（以下为教学示例假设值，非实测）：船长 L=150 m，航速 U=7.7 m/s（约15 kn），取无量纲 K'=0.15、T'=2.5，则 K = K'·U/L ≈ 0.0077 1/s，T = T'·L/U ≈ 48.7 s。若采用二阶模型，可设 T₁=40 s、T₂=8.7 s，使 T₁+T₂≈48.7 s，保持与一阶等效的一致性。

L = 150; U = 7.7;
Kp = 0.15; Tp = 2.5;
K = Kp*U/L;      % ≈ 0.0077 1/s
T = Tp*L/U;      % ≈ 48.7 s
G1 = tf(K, [T 1]);              % 一阶
G2 = tf(K, conv([40 1],[8.7 1])); % 二阶
step(G1, G2, 200); grid on;
legend('一阶','二阶');

运行阶跃响应后应观察两点：一是稳态值是否等于 K·δ（例如 δ=10°=0.1745 rad，则 r_ss≈0.00134 rad/s）；二是二阶曲线在初始段是否比一阶更平缓——这正是 T₂ 带来的初始延迟效应。若两者稳态值一致而动态过程不同，说明参数换算正确。

适用场景与选型建议

一阶模型结构简单、参数少，非常适合航向保持、航迹保持等低速常规控制器的初步设计与PID参数整定，也是教学与快速验证的首选。但当船舶尺度大、航向稳定性弱（如大型油轮、集装箱船），或控制精度要求高、需要分析高频舵机动态时，一阶模型会低估相位滞后，此时应改用二阶模型。实际工程中，常先用一阶模型完成控制器框架设计，再用二阶模型校核鲁棒性，形成"先简后繁"的仿真流程。