4. PID控制原理回顾

自动舵的核心任务，是让船舶的实际艏向角 ψ 尽可能平稳、准确地跟随设定艏向 ψ_r。操舵产生的舵角 δ 会改变船舶的转艏力矩，而 PID 控制器正是根据“偏差 e = ψ_r − ψ”来计算舵角指令的经典算法。本节把 PID 的三项作用、时域响应特征、稳定性判断以及离散化实现串起来，为后续在仿真中整定参数打基础。

一、三项作用各自解决什么问题

PID 的输出可写成连续形式：

δ(t) = Kp·e(t) + Ki·∫e(τ)dτ + Kd·de(t)/dt

比例项 Kp·e 提供“偏差越大、打舵越狠”的即时响应，是控制作用的主体。Kp 偏小，船舶转向迟缓、回舵无力；Kp 偏大，舵角动作猛烈，容易引起艏向过冲和来回摆动。对自动舵而言，比例项还常与“舵角限幅”配合，避免指令超出实际舵机行程（一般商船舵角限制在 ±35° 左右）。

积分项 Ki·∫e 负责消除稳态误差。船舶在风、流、不对称负载作用下会产生持续偏航力矩，仅靠比例项会留下一个固定静差，船始终“偏一点”。积分项累积历史偏差，逐步把舵角推向能抵消扰动的方向。但积分过强会带来积分饱和与滞后，表现为艏向缓慢爬升后突然超调，因此工程上常加积分限幅或抗饱和处理。

微分项 Kd·de/dt 感知偏差变化趋势，起“提前刹车”作用。船舶惯性大、响应滞后，微分项能在艏向接近设定值前提前回舵，抑制过冲、增加阻尼。代价是对测量噪声敏感——罗经艏向信号若有抖动，微分项会放大成舵机的高频动作，所以实际自动舵常对微分项做滤波，或采用“微分先行”结构。

二、时域响应特性与稳定性概念

给自动舵一个阶跃设定（例如从 000° 改向到 030°），典型的闭环响应包含上升时间、超调量、调节时间和稳态误差几个指标。比例增益提高会缩短上升时间但增大超调；积分增益提高能压掉静差却延长调节时间；微分增益适当增加可减小超调、改善平稳性。三者相互耦合，这正是需要整定的原因。

稳定性是比“响应快慢”更根本的要求。若闭环极点位于复平面右半部，系统发散；临界稳定则表现为等幅振荡。对船舶这种大惯性、长滞后对象，过高的 Kp 与 Ki 组合最容易把系统推向持续摆艏。工程判断上，除了看阶跃响应是否收敛，还可观察艏向振荡的衰减比：相邻两个波峰的幅值比明显小于 1 且逐步衰减，说明阻尼充足；若接近 1，则已逼近稳定边界。

三、从连续到离散：自动舵里的 PID 实现

实际自动舵以固定采样周期 T 运行，控制器必须写成离散递推形式。设第 k 个采样时刻的偏差为 e(k)，常用位置式与增量式两种写法。位置式直接计算舵角指令：

e(k)   = ψr(k) − ψ(k)
P(k)   = Kp · e(k)
I(k)   = I(k−1) + Ki · T · e(k)      // 积分累加，需限幅
D(k)   = Kd · (e(k) − e(k−1)) / T
δ(k)   = P(k) + I(k) + D(k)          // 再经舵角限幅 ±35°

增量式只计算舵角的变化量，天然具备抗积分饱和的优点，也更适合舵机按“增量打舵”的机构：

Δδ(k) = Kp·[e(k)−e(k−1)] + Ki·T·e(k)
        + (Kd/T)·[e(k)−2e(k−1)+e(k−2)]
δ(k)  = δ(k−1) + Δδ(k)

离散化时有几个工程要点必须注意。第一，采样周期 T 要与船舶操纵性匹配：船舶艏向响应时间常数通常在数十秒量级，自动舵采样周期取 0.5～2 s 已足够，过小会放大微分噪声，过大则微分项失真。第二，微分项对量化噪声敏感，可改用“对测量值微分”并加一阶低通滤波，例如 D(k) = α·D(k−1) + (1−α)·Kd·(ψ(k−1)−ψ(k))/T。第三，积分项必须限幅，否则长时间大偏差后舵角指令会“卡死”在极限位置，回舵迟缓。

四、一个可复现的整定观察示例

假设某船简化 Nomoto 模型为 T·ψ̇ + ψ = K·δ，取 T = 20 s、K = 0.8，采样周期 T_s = 1 s，设定从 000° 阶跃到 030°。可先只开比例项，取 Kp = 2.0，观察响应：艏向上升较慢，最终稳定在约 026°～027°，存在明显静差。再加入 Ki = 0.05，静差逐步被消除，但若 Ki 加到 0.3，会出现明显超调并伴随缓慢摆动。此时加入 Kd = 8～15，超调显著减小，回舵更提前。上述数值属于教学示例假设，实际整定需结合具体船型、装载与海况，通过仿真或实船试验确认。