K3-10 · 控制、软件与标定 / 控制策略与算法
待认领listed
→
众筹中recruiting
→
已开课open · Course 售卖
一课题可并存多讲师版本
课程大纲
学完能做什么
- 能按「模型信息够不够、被估量带宽、域控算力预算」三条判据，在守恒式直算／部件特性曲线反解／一阶 RC 网络／龙贝格观测器／互补滤波／KF-EKF-UKF 六类结构中选定一种，并说清落选结构各自的失效点
- 能算可观性矩阵的秩与条件数，判定某个不可测量到底估不估得出来；不可观或数值病态时给出「加测点／降阶」而不是「换更复杂滤波器」的处置
- 能从静置段数据量出 R、按可接受跟踪滞后反推 Q，整定龙贝格增益 L，并算出常值模型失配下的稳态估计偏差 e_ss
- 能把残差带做成工况化的可标定量（ε 工况分段表 + 权重调度 + 有效工况窗口三张表），并设计发散防护与重初始化触发清单
- 能把设计期的制冷量／换热量公式改造成域控每周期可算的在线量，交付「值 + 置信度 + 有效窗口」三件套，并给出误差传播与拒绝出值的判据
- 能设计置信度三档分级与退出状态机，与 fault flag 三级、降级等级体系对接，并判定某测点的软测量属「可替代／只可交叉校验／不可替代」硬件冗余哪一档
- 能在壁温类／量测重构类／源项反解类三种估计结构之间做迁移，对一个新的不可测量走完「定被估量→选结构→整定→写退出条件」四步
内容大纲
- 估计器结构谱系与车规选型
- 六类结构按「需要什么模型信息／吃多少算力／抗什么误差」排成一张选型表：守恒式直算（只需一条平衡式、无状态，误差全从输入进来）→ 部件特性曲线反解（需要 map，误差随 map 外推放大）→ 一阶 RC 网络（一个状态，抗量测噪声但慢）→ 龙贝格观测器（需要 A/C，增益按极点配置，不需要噪声统计）→ 互补滤波（两路频段互补，不需要协方差，算力最省）→ KF/EKF/UKF（需要 Q/R 噪声统计，算力最贵）。选型判据不是「谁最先进」，是这三栏对不对得上
- 选结构之前先判可观性：rank(O)=n 是「理论上能不能估」的是非题，条件数 κ(O) 才是「工程上估不估得住」的程度量；秩满而 κ 偏大时量测噪声会被放大到估计量失去意义，正确动作是加测点或把不可观那部分降阶掉。B4-04 只把「哪些温度必须实测、哪些可以观测器估」留到选型口径，定量判据在本讲
- 守恒式直算与观测器的分界线：被估量若能由同一时刻的量测代数解出（如 Q̇=ṁ·c_p·ΔT），就不要引入状态——多一个状态就多一套整定与发散风险；只有当被估量有显著热惯性、或量测比对象动态更慢/更噪时才需要状态与增益
- EKF 与 UKF 的取舍：EKF 需要工作点雅可比，强非线性或大初始误差下线性化误差直接进协方差且不体现在 Q 里；UKF 免雅可比但采样点随维数增长（单步传播量级约 2n+1 次）。热管理里状态维数低、非线性主要来自 map 与饱和曲线，多数场合 EKF 够用，UKF 留给饱和线附近强弯曲的情形
- 互补滤波是本领域被低估的一档：凡是「有一路慢但无漂的基准 + 一路快但会漂的积分」的场合（冷却液入口温度配损耗积分、长静置读数配热流积分），一个交叉时间常数 τ 就是全部整定量，不需要协方差与矩阵运算——车规实现里这是最先该试的结构
- 车规约束会把选型压回定点/查表实现：无动态内存分配、矩阵求逆维数受限（实做常把 KF 拆成逐量测的标量更新以回避矩阵反演）、单步耗时须满足 t_solve ≤ Δt。实时性判据本体归 J7-02/J7-04，本讲只讲它如何反过来限制结构选择
- 本讲接住 K1-03 的书面缺口：K1-03 只列出这串结构名并声明「不展开算法推导」，冗余语境下拿残差当阈值用的取舍仍归 K1-03，本课不重复讲冗余表决与独立性
- 整定量的确定：R/Q、增益 L、稳态偏差 e_ss 与工况化残差带
- Q/R 不是拟合出来的两个数，是「信模型还是信量测」的显式声明：R 可测——用长静置段数据估方差，把量化台阶一并计入；Q 不可测——它是留给未建模项的旷量，只能按可接受的跟踪滞后反推。先定 R 再调 Q/R 比值，调绝对值没有意义
- 龙贝格增益 L 的整定判据：观测器极点比被观测动态快 3~5 倍即可，再快只是把量测噪声按 L 的倍数灌进估计——同一句话在 KF 里表现为 Q/R 抬高导致 K 变大，两条路整定的是同一个 trade-off
- 常值模型失配下的稳态偏差 e_ss = −(A−LC)⁻¹·E·d 把整定的全部内容压成一对矛盾：增大增益缩小 e_ss 但放大噪声。结论是别指望靠加大增益吃掉系统性失配——那要靠在线参数辨识，或把这段偏差如实写进误差带交给下游
- 残差带必须工况化：用创新协方差 √(H·P⁻·Hᵗ + R) 乘覆盖系数得 ε(工况)，而不是全工况拍一个常数。这就是把 K1-03 那句定性的「置信度随工况动态调整」落成三张可标定表——ε 的工况分段、权重调度、有效工况窗口
- 有效工况窗口要显式写出来并作为退出条件的第一道闸：小 ΔT、瞬态、模式切换过渡段、输入信号缺失——窗外一律不出值，不许硬算（与 G12-01 拒绝小 ΔT 热平衡反算同源）
- 运行期自保与在线参数辨识：发散防护、重初始化与冻结回退
- 发散防护三件套：协方差上下钳位（防 P 收敛到 0 后滤波器再也不吃量测）、对称化与正定性维护（P←(P+Pᵗ)/2）、创新门限（|r| 超 k·σ 的量测先拒收再计数）。EKF 发散最典型的实机表现是长时间稳定运行后某次工况突变一去不回，而不是一上车就飘
- 重初始化的触发清单：模式切换、执行器控制权交接、长时间无有效量测后恢复、置信度跌到不可用——一律以慢变但无漂的边界量测重置状态与协方差。E2-05 用冷却液入口温度重对齐 Cauer 网络就是这条的一个实例，本课给通用触发清单，结温那一路不重讲
- 在线参数辨识与状态估计是两件事：辨识改的是模型参数（UA 衰减、泵特性下移、map 偏置），估计重构的是不可测状态。辨识须带遗忘因子，并规定冻结条件（激励不足、可辨识性病态、正在降级）与回退到出厂值的路径；同一残差不能既归参数又归状态，实做上让辨识跑慢时间尺度、状态估计跑快时间尺度
- 本讲只到「辨得出来、且辨错了能回退」：K7-03 从能效自学习侧前指本课、K3-03 的阈值老化在线修正也前指本课；衰减趋势判定与服务决策归 K5-05。传感器级在线校零（被测量没变、单只传感器零点漂了）与系统级参数辨识（被测系统真的变了）的判别依据是残差在多通道上的一致性，传感器侧口径归 K1-01/K1-03，本课只接系统级那一半
- 在线热负荷、制冷量与换热量估计
- 设计期公式与在线量的差别不在公式本身，在参数怎么来：η_v、c_p、h(p,T) 在设计期是给定值，在域控里要么是 map 外推、要么随乙二醇浓度与液温漂移、要么本身就是另一个估计量——所以每一个在线量都必须挂一条「参数来源 + 误差账」，不能只交一个数
- 制冷量在线化 Q̇_e = ṁ_ref·(h₁−h₄)、ṁ_ref = ρ_s·V_disp·n·η_v（公式本体承 A1-05 与 C1-01/G1-01）：三个误差源按量级排——η_v map 在低转速/大压比区的外推最差、吸气密度 ρ_s 依赖 (p_s,T_s) 两路的时刻一致性、吸气带液时 h₁ 的单相反算前提直接失效。两相判别与 P-T 联测口径归 G12-02，本课只管把它做成带置信度的在线量
- 水侧换热量 Q̇_w = ṁ_w·c_p·ΔT_w 与 K5-06 方向相反——K5-06 已知 Q̇ 求 ṁ，本课已知 ṁ 求 Q̇。两者不可循环互喂：用热平衡估出的 ṁ 再回代算 Q̇ 是恒等式，看着自洽实际零信息，且残差恒为零、诊断抓不到任何异常。这是最容易犯的实现级错误
- 空气侧用焓差法 Q̇_air = ṁ_air·(h_in − h_out)（承 A2-04 的湿空气口径）：湿工况下必须用焓差而不是显热温差，否则蒸发器的潜热被整体漏掉；ṁ_air 在整车上通常也是由鼓风机工作点反推的估计量，它的误差要串进这条链而不是当已知量
- 误差传播按相对量做：δQ/Q ≈ √((δṁ/ṁ)² + (δΔT/ΔT)²)，因此 ΔT 小的工况这个量先失效。判据给量级口径——ΔT 须大于温度传感器误差带若干倍才出值——而不是给一个假装精确的阈值，各车型传感器误差带不同
- 三处悬空输入的交付口径统一为三件套「值 + 置信度档 + 有效窗口标志」：K3-05 前馈通道要的负荷突变量、K3-06 状态空间里的不可测扰动 d、K3-07 的各路 Q_demand,i 与 Q_avail。下游按置信度选分支（置信不足→退保守分支），而不是无条件消费一个裸值
- 估计误差必须折进约束裕度，不能默默忽略：MPC 按估计值求解、实车按真值运行，硬约束会在边缘工况被悄悄吃掉。误差带→裕度的折算关系本体在 B4-04，本课只交误差带本身及它随工况的变化
- 能力量还有第二路来源——执行器回传的能力上限，它比模型估算更直接但受各供应商状态字语义不统一之累。两路误差性质不同，不可简单取平均；两路怎么进仲裁归 K5-02/K3-01，本课只交模型估算这一路的误差性质
- 置信度分级与退出状态机
- 置信度不是一个 0~1 的连续数丢给下游，而是三档离散状态 {可信／降级可用／不可用} 加一个有效窗口标志，与 K1-01 的 fault flag 三级同一套语法。K1-03 已要求置信度分级对齐这三档并显式输出，本讲交的是状态机本体与各档判据
- 每档的进入判据分四层：残差超带的持续时长（|r| > ε 持续 > t_hold）、是否落在有效工况窗口之外、输入自身的置信度（估计量以别的估计量为输入时置信度要串乘，不能取最好那一路）、协方差或条件数是否越界
- 退出（降级为不可用）必须绑定下游动作，否则等于没退出：转保守常量假设／切硬件冗余通道／再降一级，三选一的判据是「这个量在闭环里承担什么角色」。参与安全限值的量禁止「保持上次值」——K1-01 已写死此条，本课引用不重讲
- 反向命题正是全库空白：现有课只有正向的「切到软测量估计」，本讲交「软测量何时不可信」。K5-02 已把「估计失效、置信度跌破门限」列为独立于硬件失效的一类降级触发源（看不见了 ≠ 坏了，两条降级路径不同），本课交它的判据，降级等级体系与资源仲裁归 K5-02
- 何时必须让位给硬件冗余：解析冗余的模型误差随工况漂移，且当模型与被污染的量测共享同一输入时残差恒为零（共模盲区，K1-03 已论证）。本讲给「可替代／只可交叉校验／不可替代」的分档判据；冗余架构、表决结构与独立性论证本体归 K1-03
- 置信度档位切换要带迟滞与最小驻留时间，否则残差在带边抖动会把下游降级状态机一起带抖；不许靠抬高 ε 应付——抬 ε 同时抬掉真失效的检出能力（与 K1-03 对 ε_th 的告警是同一句话）
- 回切（不可用→可信）不是退出条件的镜像：须满足连续 N 帧有效 + 状态已重初始化 + 首拍不得让执行器跳。无扰接管的手法本体归 K3-05，降级态回交自动态那一拍的责任归 K5-02
- 置信度状态要与 fault flag 分开报：估计不可信不一定是硬件故障，不该置 DTC、不该点故障灯，但必须有一组独立于 DTC 的可读状态量（触发源、生效档位、进入/退出时刻）。这与 K5-02 对「非故障策略触发」的规定同源，字段承载归 K5-01/K5-02
- 融合置信度被用作检测判决时另有一层口径——K5-04 用多特征融合置信度做热失控早期预警的误报/漏报权衡，那是判决层的事；本课只交估计层的置信度分档与退出，不重讲预警判决的权衡
- 边界与三类结构的迁移实操
- 壁温类估计（范例：玻璃内表面温度）：对象特征是低热容、大面积围护壁面，输入含日照短波吸收与风窗送风量/送风温度，估计器形态是单节点集总或两节点 RC。物理机理、最冷点偏置（风窗下沿、A 柱根部比传感点更冷）与标定口径全部归 C5-03，本课只讲这类结构的整定与残差带
- 量测重构类估计（范例：虚拟内温）：对象特征是量测点本身被污染——仪表台辐射与滞留空气使内温读数系统性偏高，须由多点温度 + 送风温度 + 风量 + 车速加权重构，权重随风量与车速调度（低风量时传感点代表性最差、权重最低）。输入清单与权重调度口径归 C6-01，传感器误差机理与抽吸取样归 G12-04，本课只讲重构结构与权重的置信度含义
- 源项反解类估计（范例：产湿率）：对象特征是被估量根本不在状态里、而是守恒式里的源项，由座舱湿平衡反解。可辨识性边界必须写死——新风量未知或内外循环比例正在变化时该反解退化，只能给「湿源在增/在减」的趋势、不能给绝对值，此时防雾侧一律取保守值。对象侧归 C4-04
- 与既有两个应用样板对照即凑齐五种形态：冷却液流量（多路互校、逐路降权而非一超阈就踢掉，构造路径见 G12-01 讲义 5.1，估计器与流量/内漏诊断本体归 K5-06）、结温（Cauer 网络在线观测器全流程，含实时性与模型自校正，样板归 E2-05）。本课不重讲任何一个，只用它们做结构对照
- 分工总账一句话说清：对象侧物理机理、冷点偏置与标定口径归 C5-03/C6-01/C4-04；模型结构与参数辨识（座舱多节点划分、UA 与太阳有效开口的反辨识）归 J1-08；孪生侧融合架构、e_drift 漂移监测与 t_solve 实时约束归 J7-04；降阶与状态空间形式归 J7-02；本课只交估计器形态差异、整定、残差带与退出条件
- A4-08 原先引的「K3-05 状态估计」是错引——K3-05 实为《PID/前馈/增益调度控制实战》，与状态估计无关；t_eq 软传感器的整定、残差带与退出条件改指本课，本课须接住。对象侧的舒适量定义与分区送风策略仍归 A4-08
- 迁移练习固定四步走，任何一个新的不可测量都用同一套流程：定被估量与可观性 → 按第 1 讲选型表选结构 → 整定 Q/R 或 L 并给出残差带的工况分段 → 写退出条件与对应的下游降级动作
- 本课自身的深度边界也要讲明：滤波理论的最优性与收敛性数学证明不在范围内（只给工程整定判据与失效边界）；数据驱动/黑箱的能效自学习归 K7-03；部件级传感器误差机理归 G12-04/K1-01；估计失准属「功能不足而非失效」，其安全论证方法学不在本课
关键公式
O = [C; CA; CA²; …; CA^(n−1)]，可估条件 rank(O) = n；数值可估程度用条件数 κ(O) = σ_max/σ_min 衡量
选结构之前的第一道闸——判某个不可测量到底估不估得出来；不可观或数值病态时的正确动作是加测点或把不可观部分降阶掉，不是换更复杂的滤波器
dx̂/dt = A·x̂ + B·u + L·(y − C·x̂)，增益 L 按「观测器极点比被观测动态快 3~5 倍」配置
模型信息够、噪声统计说不清时的默认结构；热管理里壁温、绕组热点这类单点低维被估量最常用
e = x − x̂ 满足 de/dt = (A − L·C)·e + E·d ⇒ 常值失配下 e_ss = −(A − L·C)⁻¹·E·d
回答「模型少建了一项，估出来会偏多少」，是下游误差带与约束裕度的直接来源（误差带→裕度的折算归 B4-04）
x̂ = LPF_ωc(x_slow) + HPF_ωc(x_fast)，一阶实现 x̂[k] = α·(x̂[k−1] + Δx_fast[k]) + (1−α)·x_slow[k]，α = τ/(τ + Δt)，τ = 1/ω_c
在「一路慢基准 + 一路快积分」的场合替代 KF；τ 是唯一整定量，取值落在慢基准的漂移时间尺度与快通道的积分漂移时间尺度之间
P⁻ = A·P·Aᵗ + Q；K = P⁻·Hᵗ·(H·P⁻·Hᵗ + R)⁻¹；P = (I − K·H)·P⁻；整定量是 Q/R 的比值
状态更新式的框架侧口径已在 J7-04，本课交的是这三式里的整定量——Q/R 怎么定、P 怎么防发散
r = y − H·x̂⁻；ε(工况) = k·√(H·P⁻·Hᵗ + R)；判据 |r| > ε 持续 > t_hold ⇒ 置信度降档
把 K1-03 那句定性的「置信度随工况动态调整」落成可标定量：ε 的工况分段表 + 权重调度 + 有效工况窗口三张表
θ̂[k] = θ̂[k−1] + K_θ[k]·(y[k] − φᵗ[k]·θ̂[k−1])，带遗忘因子 λ ∈ (0,1]，有效记忆长度 ≈ 1/(1−λ)
处理「模型参数真的变了」那一类失配（UA 衰减、泵特性下移、map 偏置），即稳态偏差式吃不掉的部分；同时规定冻结条件与回退到出厂值的路径
Q̇_e = ṁ_ref·(h₁ − h₄)，ṁ_ref = ρ_s·V_disp·n·η_v（公式本体承 A1-05 与 C1-01/G1-01，本课只做在线化）
把设计期的循环性能公式做成域控每周期可算的在线能力量，补上 K3-07 的 Q_avail 与各路 Q_demand,i 的悬空输入
Q̇_w = ṁ_w·c_p·ΔT_w；Q̇_air = ṁ_air·(h_in − h_out)（湿工况必用焓差，承 A2-04）；δQ/Q ≈ √((δṁ/ṁ)² + (δΔT/ΔT)²)
各路热需求与实际交付量的在线口径，兼作两侧互校（冷媒侧 Q̇_e 与水侧 Q̇_w 交叉校验）
关键概念
估计器结构谱系（守恒式直算／特性曲线反解／一阶 RC／龙贝格／互补滤波／KF-EKF-UKF）可观性矩阵的秩与条件数 κ(O)观测器极点配置与增益 L互补滤波交叉频率 ω_cQ/R 整定与创新协方差残差带工况化 ε(工况)覆盖系数与消抖时长 t_hold发散防护（协方差钳位、对称化、创新门限）观测器重初始化触发常值失配下的稳态估计偏差 e_ss在线参数辨识（遗忘因子、冻结条件、回退路径）有效工况窗口置信度三档分级（可信／降级可用／不可用）退出状态机与回切条件置信度串乘（估计量以估计量为输入）误差传播（相对量平方和）软测量对硬件冗余的可替代性分档壁温类／量测重构类／源项反解类三种估计结构值 + 置信度 + 有效窗口三件套
推荐工具与标准
MATLAB/Simulink（观测器与 KF/EKF 原型；Control System Toolbox 算可观性矩阵秩、条件数与极点配置） Simulink Coder 与定点定标工具链（估计器的车规实现：无动态分配、逐量测标量更新） dSPACE / Speedgoat 等实时台架（t_solve ≤ Δt 的实时性核算与注入试验） CANape / INCA（残差带工况分段表、置信度门限的实车标定与数采） Python / MATLAB（用静置段数据估 R、创新序列白化检验、协方差与条件数的离线诊断）
ISO 26262（估计值进入安全相关闭环时的失效响应与降级路径；具体部与版本按现行版本确认） ISO 21448 / SOTIF（估计失准属「功能不足而非失效」，其安全论证方法学归 K2-05；版本按现行版本确认） 注：本课不引用任何规定观测器算法或整定量的标准——此类量值一律来自平台标定，无标准可依
工程案例
某纯电平台用热平衡法在线估算电池回路换热量，喂给功率仲裁当 Q_avail 用。台架与常温工况都对得上，冬季小负荷段整车却间歇出现能力误判：电池冷却过早退出、随后温度追上来又急切回，用户端表现为快充功率反复台阶。查下来两条根因叠在一起——一是小负荷段 ΔT 落到与温度传感器误差带同量级，δΔT/ΔT 把相对误差放大到不可用，而实现里根本没有有效工况窗口，照样每周期出值；二是这个 Q̇ 又被回代去反算流量、流量再回来算 Q̇，构成循环互喂，残差恒为零，所有残差类诊断全部放行，看起来「自洽得很好」。整改三条：加有效工况窗口（ΔT 低于误差带若干倍即置不可用，窗外不出值）；把置信度做成三档显式输出给仲裁，不可用时仲裁退保守分支而不是继续消费裸值；切断 Q̇ 与 ṁ 的循环，流量改由泵特性与电功率两路独立提供、逐路降权互校（那一路的估计器本体归 K5-06）。
动手做
交付物 · 在三类范例中任选一个被估量（壁温类=玻璃内表面温度／量测重构类=虚拟内温／源项反解类=产湿率），走完四步并交付四件材料：① 写出状态方程与量测方程，算可观性矩阵的秩与条件数，给出「可估／数值病态／不可估」的判定及对应处置（加测点还是降阶）；② 在同一被估量上并行实现两种结构（如一阶 RC + 互补滤波，或龙贝格 + EKF），用同一段工况数据对比跟踪滞后与噪声放大倍数，说明为什么最终选它；③ 用静置段数据量 R、按可接受滞后反推 Q（或配置 L 的极点），并给出残差带 ε 的工况分段表（至少覆盖大负荷、小负荷、瞬态、输入缺失四段）；④ 写置信度三档的进入/退出/回切判据与每档对应的下游动作，然后做一次注入试验（把某路输入置为丢失或加恒定漂移），验证退出真的发生、且回切时执行器没有阶跃。交付物：可观性判定表 + 两结构对比曲线 + 残差带工况分段表 + 置信度状态机与注入试验记录。对象侧的物理机理与标定口径不必自己推，直接引 C5-03／C6-01／C4-04。
常见误区
- 以为选更高级的滤波器就能把估不出来的量估出来——其实不可观是结构问题，κ(O) 已经很大时 EKF/UKF 只是把噪声放大得更好看；正确动作是加测点或降阶
- 以为 Q 和 R 是两个待拟合参数——其实 R 可以从静置段数据量出来，Q 是留给未建模项的旷量、只能按可接受滞后反推；把两个都当自由参数去拟合，会得到一个只在标定工况成立的滤波器
- 以为残差带定一个常数就够——其实创新协方差随工况变，常数带在低信噪工况必误报、高信噪工况必漏报，两头都错
- 以为残差抖动可以靠抬高阈值压住——其实抬 ε 同时抬掉真失效的检出能力（与 K1-03 对 ε_th 的告警同源）；该做的是迟滞、最小驻留时间与有效工况窗口
- 以为置信度低就让下游「少信一点、按权重打个折」——其实下游需要的是离散档位加明确动作；一个打了折的估计值进了安全限值闭环，等于把保护也打了折
- 以为有了软测量就可以省掉硬件冗余——其实两者防的失效不同，且解析冗余若与被污染的量测共享输入，残差恒为零（K1-03 的共模盲区），安全相关测点不可替代
- 以为热平衡式既能算流量又能算换热量——其实一条式子里两个未知量只能定一个；用热平衡估出的 ṁ 回代算 Q̇ 是恒等式，零信息且残差恒为零
- 以为估计误差可以留给下游「自己留点裕度」——其实不把误差带写出来，MPC 就会按估计值当真值求解，硬约束在边缘工况被悄悄吃掉（折算关系见 B4-04）
- 以为模式切换后观测器会自己收敛回来——其实积分历史带着切换前的状态，切换后若干秒是系统性偏移；须以慢变无漂的边界量测重初始化
- 以为参数在线辨识就是让它一直跑着自适应——其实激励不足时必须冻结，否则参数被噪声牵走，且错值会跨上电周期一路信下去；冻结条件与回退路径是辨识功能的必要组成
- 以为传感器在线校零和系统参数辨识差不多，一起做更省事——其实前者假设被测量没变、后者假设系统真变了，混做会把真实劣化当传感器漂移校掉；判别依据是残差在多通道上的一致性
相关课题