更新记录

1.0.0（2026-11-06）

首个版本：威布尔最大似然拟合，得分方程二分法解形状参数 β，支持右删失样本，出 β、η、B1/B10/B50、MTTF，并用中位秩回归给拟合质量 R²；按 β 区间判读失效模式（β<0.8 早期失效、0.8–1.2 随机失效、>1.2 磨损型），提示该做筛选还是改设计；寿命分位按 −ln(1−F) 反解，可查任意累积失效比例对应的寿命；阿伦尼乌斯加速因子按 Ea/k(1/T_use−1/T_stress) 计算，支持摄氏与开尔文双入口，反推活化能并给出常见失效机理的 Ea 参考区间；TM-21 光衰外推以 Φ(t)=B·e^(−αt) 拟合，出 L70/L90 与 α、B，外推时长按实测最后一点六倍封顶，超限时明确标注已封顶；多组数据并排对比与文本解析（逗号/制表符/分号/竖线，表头注释自动跳过），三份 CSV 带 BOM 导出。UTS→Kotlin 编译通过（index.kt 1441 行、classes.dex 63724 字节），116 项回归测试全部通过；真机未执行（测试机锁屏，adb 报 deviceLocked=1），市场截图为 H5 渲染并标注示例数据。

平台兼容性

uni-app

可靠性寿命推算

一批货到底能用多久？B10 是多少？85℃ 跑 1000 小时能折成常温几年？灯珠光衰到 70% 要多久？

三个问题分别对应威布尔拟合、阿伦尼乌斯加速、TM-21 光衰外推。都是纯算法，数据你给，结论它算。

纯本地计算，不联网、不上传、不需要任何权限。仅支持 Android。

购买前必读

• 本插件算寿命，不做数据采集：不接老化架、不接光谱仪、不接 PLC。试验数据由你的设备或台账提供。
• 威布尔用最大似然（MLE），不是最小二乘；支持右删失（到截止时间还没坏的样本）。 与 Excel 里用「中位秩 + 线性回归」算出来的 β 会有差异，这是方法差异，不是算错。 插件同时给出回归法的 β 与 R² 作参照，两个口径差得远时先查混批。
• 时间单位由你决定（小时、循环次数、公里都行），插件不做换算；三个功能的单位必须一致。
• 加速因子与 TM-21 都是模型外推，结论只对模型成立。Ea 取错，寿命就错一个数量级 —— 报告里必须写清 Ea 来源（实测还是经验值）与试验条件。
• TM-21 的 6 倍限制是硬规则：外推超过最后实测点 6 倍时，插件只报到 6 倍并标注「受限」， 不会给你一个好看但站不住的整数。

解决什么问题

1. 威布尔：这批货的寿命分布

平均寿命是最容易被误用的指标 —— 一批货里 80% 在 2000 小时坏了、20% 撑到 20000 小时， 平均值可能还有 5600 小时，但客户在第 2000 小时就开始退货。

威布尔给的是分布：形状参数 β 说明失效是「早期夭折」「随机」还是「磨损」， 尺度参数 η 是 63.2% 失效对应的寿命，B10 是 10% 失效对应的寿命 —— 后者才是客户体验到的那个数。

2. 阿伦尼乌斯：加速试验折几年

85℃ 跑 1000 小时，等于常温跑多久？取决于激活能 Ea：

AF = exp( Ea/k × (1/T_use − 1/T_stress) )，k = 8.617333262e-5 eV/K

Ea = 0.45 eV、25℃ 使用、85℃ 试验时 AF ≈ 18.8，1000 小时只折 2.1 年。 常听到的「温度每升 10℃ 寿命减半」只在 Ea ≈ 0.7 eV 时成立，照搬会高估加速效果。

3. TM-21：光衰外推到 L70

对 LM-80 的光通维持率序列做指数拟合 Φ(t) = B·e^(−αt)，解出 L70 = ln(B/0.7)/α。 TM-21 规定外推不得超过最后实测点的 6 倍，插件按这条规则封顶。

与同类的边界

顺序建议：先 szhzy-msa-grr 证明量具可信 → szhzy-spc-capability 看过程能力 → 本插件回答「能撑多久」。三者数据不重叠，算法也不重叠。

快速上手

威布尔拟合

import { weibullFit } from '@/uni_modules/szhzy-life-predict'

const out = weibullFit({
  name: '驱动电源电解电容',
  times: [520, 680, 900, 1100, 1500, 6000, 6000],   // 小时
  censored: [0, 0, 0, 0, 0, 1, 1]                    // 1 = 到截止时间还没坏
})

if (out.ok) {
  console.log(out.beta, out.eta)   // 形状参数、尺度参数
  console.log(out.b10)             // 10% 失效寿命
  console.log(out.modeText)        // 早期失效 / 随机失效 / 磨损失效
  console.log(out.r2)              // 中位秩回归拟合优度
}

阿伦尼乌斯

import { arrheniusCalc, activationEnergy } from '@/uni_modules/szhzy-life-predict'

const a = arrheniusCalc({ ea: 0.45, tUseC: 25, tStressC: 85, stressHours: 1000 })
// a.af → 18.8，a.equivalentHours → 18807，a.equivalentYears → 2.15

// 有两个温度下的加速因子时，反算激活能
const ea = activationEnergy({ af: 18.8, t1C: 25, t2C: 85 })

TM-21 光衰

import { tm21Fit } from '@/uni_modules/szhzy-life-predict'

const m = tm21Fit({
  name: '灯珠 2835',
  times: [0, 500, 1000, 2000, 3000, 4000, 5000, 6000],
  lumen: [1.0, 0.995, 0.988, 0.976, 0.965, 0.953, 0.942, 0.931]  // 0.98 或 98 都认
})
// m.l70 → 外推值，m.l70Report → 受 6 倍上限约束后的可报值，m.capped70 → 是否受限

多组对比

import { lifeBatch } from '@/uni_modules/szhzy-life-predict'

const b = lifeBatch({
  items: [
    { name: '电解电容', times: [...], censored: [...] },
    { name: '焊点', times: [...], censored: [...] }
  ]
})
// b.worstName / b.worstB10 —— 直接点名最短的那组

错误码与提示

β 触到搜索边界（0.05 或 50）时不会假装收敛，会在 warn 里写明「数据不足或分布退化」。

验证状态

• UTS → Kotlin 编译通过：Project demo-life-predict compiled successfully，产物 index.kt 1441 行、 classes.dex 63724 字节；
• 纯逻辑回归测试 116 项全部通过：Γ 函数、B10/B50 解析式、已知参数反拟合、右删失、 加速因子与手算逐位对照、Ea 反算、TM-21 参数精确还原与 6 倍封顶、解析与 CSV 结构；
• 真机未验证：本轮测试机（HUAWEI CET-AL60 / Android 12）锁屏（deviceLocked=1）， 无法安装与操作界面，页面交互与真机文件读写均未实测；
• 市场截图为 H5 构建产物在浏览器中的真实渲染，页面标注「示例数据」；
• iOS / 鸿蒙 / 小程序 / Web 均未实现，平台声明一律为不支持。