3. 多区加热系统状态空间建模

多区加热系统的控制难点，根源在于“每个加热区并不是孤立的”。当某一区升温时，热量会通过固体导热、对流与辐射向相邻区扩散，导致相邻区温度被动上升，形成区域间热耦合。若仍按单回路 PID 独立整定，就会出现“调一区、动多区”的反复振荡。状态空间建模的价值，正是把这种多输入多输出（MIMO）耦合关系写成统一的矩阵形式，为后续解耦控制、模型预测控制提供可计算的基础。

一、从能量平衡到状态方程

设系统划分为 n 个加热区，第 i 区的温度为状态变量 T_i。对每个区做能量平衡：单位时间内区内的热量变化率，等于加热器输入功率减去向环境散失的热量，再减去流向相邻区的热量。以第 i 区为例：

C_i · dT_i/dt = P_i(t) − (T_i − T_amb)/R_i,amb − Σ_j (T_i − T_j)/R_ij

其中 C_i 为该区等效热容（J/K），P_i 为加热功率（W），R_i,amb 为对环境的热阻（K/W），R_ij 为第 i 区与第 j 区之间的耦合热阻（K/W）。把耦合项展开并整理，可写成标准状态方程：

Ṫ = A·T + B·u + E·d
y = C·T

这里 T∈Rⁿ 是温度状态向量，u∈R^m 是加热功率输入向量，d 为环境温度等可测扰动。A 矩阵的对角元为 −(1/(C_iR_i,amb) + Σ_j1/(C_iR_ij))，非对角元为 +1/(C_iR_ij)，正是耦合的数学体现。

二、区域间热耦合矩阵的物理含义

A 矩阵的非对角项决定了耦合强度。若 R_ij 很小（区与区之间导热良好），非对角元绝对值大，耦合强；若两区之间加了隔热层，R_ij 很大，非对角元趋近于零，系统近似解耦。工程上常用“耦合增益比”来评估：把第 j 区输入到第 i 区输出的稳态增益，与第 i 区自身增益相比，比值越大说明越不能忽略耦合。

需要特别说明符号约定：本文统一取“温度升高为正、热量流出该区为负”，因此 A 的对角元恒为负（自稳定），非对角元为正（邻区热量流入使本区升温）。若采用相反的散热符号约定，矩阵符号会整体变化，但物理结论不变。全文与图示均遵循这一约定。

三、输入输出变量的选择

输入变量通常取各区的加热功率或占空比。若执行器为 PWM 驱动的加热棒，输入 u_i 可取占空比 D_i∈[0,1]，此时 B 矩阵元素为 P_i,max/(C_i)，把最大功率折算进增益。若驱动器带饱和，需在模型外单独处理限幅。

输出变量按控制目标选择。若目标是各区温度均匀，则 y 直接取全部 T_i，C 为单位阵；若只关心关键区或平均温度，可令 C 为选择矩阵或行向量。例如三区系统取“区2温度 + 三区平均温度”作为输出：

C = [0 1 0
     1/3 1/3 1/3]

这种选择直接影响可观测性与后续解耦设计：输出若不含某区信息，该区状态就难以被反馈校正。

四、模型线性化与工作点

辐射散热项含 T⁴，是典型非线性来源。工程上在目标工作点 T₀ 附近做一阶泰勒展开：把辐射热流 q_rad=εσA(T⁴−T_amb⁴) 对 T 求导，得到等效辐射热导 h_rad=4εσA T₀³，并入 R_i,amb 的倒数中。这样非线性模型就转化为工作点附近的线性状态空间模型，A、B 均为常数矩阵。

线性化的有效性取决于工作点范围。若系统需在 80℃ 到 200℃ 宽范围运行，单一工作点误差会偏大，可按温度区间分段线性化，或采用增益调度。判断依据是：在目标区间内比较线性模型与非线性模型的阶跃响应，若稳态误差与超调差异在允许范围内（例如温度偏差小于 ±1.5℃），即可接受。

五、一个可复现的三区建模示例

假设三区加热板，参数取示例值：C₁=C₂=C₃=120 J/K，R_i,amb=2.5 K/W，R₁₂=R₂₃=4 K/W，R₁₃ 近似无穷（不相邻）。环境温度 25℃。则 A 矩阵对角元约为 −(1/(120×2.5)+1/(120×4)+1/(120×4))≈−0.0075 s⁻¹，相邻非对角元为 1/(120×4)≈0.00208 s⁻¹。B 取 P_max=60 W 时，元素为 60/120=0.5 K/(s·占空比)。

用这组参数做仿真，给区1 占空比 0.5、其余为 0，可观察到区1 上升的同时区2 被带动上升，区3 上升更慢——这正是耦合矩阵非对角项在起作用。若把 R₁₂ 改为 20 K/W（加隔热），区2 的被动温升会明显减小，验证了耦合强度对 R_ij 的敏感性。

完成状态空间建模后，A、B、C 三个矩阵就成为解耦器与控制器设计的直接输入。下一步的静态解耦，本质上就是求 B 的逆或伪逆来抵消耦合；而动态解耦则依赖 A 矩阵的结构信息。因此，本节建立的模型质量，直接决定后续温度均匀性能达到什么水平。