flowchart TD

 A\["模型已建立\\n状态 x=\[z,v\_p,T\_b,S,w]\\n输入 u=\[L,C,N,Ψ,T\_a]\\n终止: V\_term<=V\_cut 或 z<=0"] --> B\["0. 准备：场景/数据\\n给定 u(t)、初值 x(0)、参数 Θ"]



 %% ---- 主仿真（确定性/单轨迹） ----

 B --> C\["1. 主仿真循环：t=0→…"]

 C --> C1\["1.1 计算总功耗 P\_tot(u,x)\\n(背景+屏幕+CPU+网络/尾效应)"]

 C1 --> C2\["1.2 计算构成关系\\nV\_oc(z\_eff), R0(T\_b,S), Q\_eff(T\_b,S)"]

 C2 --> C3{"1.3 CPL 代数闭合求 I\\nΔ>=0 ?"}

 

 C3 -- "是" --> C4\["I\_CPL 分支"]

 C3 -- "否" --> C4b\["记录 t\_Δ（坍塌风险）\\n进入降级/限流策略"]

 

 C4 --> C5\["1.4 限流/降级：\\nI=min(I\_CPL, I\_max(T\_b))\\n(得到实际 I 与 P\_del)"]

 C4b --> C5

 

 C5 --> C6\["1.5 ODE RHS：\\nż, ṽ\_p, Ṫ\_b, Ṡ, ŵ"]

 C6 --> C7\["1.6 RK4 推进（每个 stage 都嵌套求 I）\\n中间时刻 u(t) 插值/生成"]

 C7 --> C8\["1.7 投影/物理约束：\\nz,S,w 截断到可行区间\\nQ\_eff>=0 等保护"]

 C8 --> C9{"1.8 事件检测：\\nV\_term<=V\_cut 或 z<=0 ?"}

 

 C9 -- "否" --> C1

 C9 -- "是" --> D\["2. 输出：轨迹 + TTE\\n(可线性插值得到 t\*)"]



 %% ---- 鲁棒性检验 ----

 D --> R{"3. 鲁棒性检验是否通过？"}



 R --> R1\["3.1 步长/收敛鲁棒性：\\nΔt vs Δt/2 的 step-halving\\n||zΔt - zΔt/2||∞ < 1e-4\\n且 TTE 相对误差 < 1%"]

 R --> R2\["3.2 事件定位鲁棒性：\\n电压/ SOC 过阈值的插值稳定\\n(比较不同插值/求根策略)"]

 R --> R3\["3.3 CPL/低电压鲁棒性：\\nΔ<0 记录 + I 饱和\\n避免非物理大电流/数值爆炸"]

 R --> R4\["3.4 约束/奇异性鲁棒性：\\nz\_eff 保护(避免 1/z 奇异)\\n状态投影抑制漂移"]



 R1 \& R2 \& R3 \& R4 --> ROK{"全部通过？"}



 ROK -- "否" --> Fix\["回到修正：\\n减小 Δt / 调整容差\\n检查限流策略/保护阈值\\n必要时回到参数标定 Θ"]

 Fix --> C



 ROK -- "是" --> E\["4. 形成确定性结论\\n场景对比、指标归因\\n(谁最耗电/何时风险高)"]



 %% ---- 敏感性 + UQ ----

 E --> S\["5. 全局敏感性（Sobol/Saltelli）\\nY=TTE 或 E\[TTE]\\n得到主效应/总效应排序"]

 S --> U\["6. 不确定性量化 UQ：\\n生成随机输入路径 u(t,ω)\\n(OU 或 切换OU)\\n对每条路径重复主仿真→TTE样本"]

 U --> U1\["6.1 统计推断：\\n经验CDF/均值方差/CI/分位数\\n生存曲线 P(TTE>t)"]

 U1 --> G\["7. 输出与建议：\\n(1) 预测 TTE(点估计+区间)\\n(2) 风险指标 t\_Δ\\n(3) 敏感因子→策略建议"]