dsh-math-calculator 2.2.0 → 2.3.1

package/CHANGELOG.md

@@ -2,6 +2,267 @@
 
 本项目遵循 [Keep a Changelog](https://keepachangelog.com/zh-CN/1.1.0/) 与 [语义化版本](https://semver.org/lang/zh-CN/)。
 
+## [2.3.1] — 2026-11-08
+
+公式库 **471 → 499**，并把材料属性从「硬设定」变成「可推导」。
+
+### 新增（宏观与财政，9 条：恒等式与严格推导）
+
+这一批与前面不同：**以恒等式为主，不是经验式**，所以不需要标「这是量级估计」。
+恒等式的价值是**能立刻否掉不闭合的设定** —— 作者编的数字若违反它，那不是
+「另一种经济学」，而是算术错了。
+
+| 组 | 条数 | 内容 |
+|---|---|---|
+| 产出与增长 | 4 | 支出法 GDP、增长核算（TFP）、趋同速度、费雪方程 |
+| 财政 | 5 | 财政余额、基本余额、债务可持续判据、铸币税、有效税率、财政自由度 |
+
+**`primary_balance` 是最有诊断价值的一条**，它把「经营亏损」与「历史债务利息」分开：
+
+```
+  PB 为正、总赤字为负 → 经营是赚的，只是被旧债利息拖垮（可救）
+  PB 为负              → 经营本身在亏，越借越亏（结构性危机）
+```
+
+**`seigniorage` 是财政与货币的唯一接口**。没有它，「赤字→印钱→通胀」这条
+最经典的崩溃机制就断在中间 —— 而那正是历史上多数政权崩溃的真实路径。
+
+来源：IMF 技术手册 TNM/2010/02（债务动态主方程）、ECB Economic Bulletin 2019（i−g 判据）、
+美国劳工统计局 BLS 的 TFP 方法（明写基于 Solow 增长核算）、印度财政部预算手册（赤字不含借债）。
+
+### 新增（材料，8 条：本构 + 微观决定）
+
+库里原有 13 条材料公式全在**构件层**（梁柱的几何与载荷），**本构层几乎空白**。
+最直接的证据是：有 `strain`（ε = ΔL/L₀，纯几何），**却没有配对的胡克定律** ——
+有应变没应力，这条链是断的。
+
+| 组 | 条数 | 内容 |
+|---|---|---|
+| 本构（应力→应变→失效） | 4 | 胡克定律、许用应力校核、断裂韧性 K_IC、Ashby 性能指数 |
+| 微观决定（由物理量推） | 4 | 理论强度、熔点与结合能、Pugh 脆韧判据、Grüneisen 热膨胀 |
+| 组织层 | 2 | Hall-Petch 晶粒强化、Wiedemann-Franz 导热导电锁合 |
+
+**这两组的分工正是「架空宇宙可大改」的落点**：
+
+```
+  作者硬设定：E、σ_y、K_IC、ρ      ← 材料「身份证」
+  公式算后果：ε、σ、K_I、M、ε̇、u   ← 由身份证推导
+```
+
+`melting_from_cohesive` 让熔点能从结合能推，于是**可以反推设计**：
+「要一个熔点 2000°C 的金属」→ 反推需要多大的键强度。
+`ideal_strength` 的「实际强度只有理论值的 1/100~1/10」则把「冶金水平」变成了一个可调参数。
+
+### 新增（生物通用，8 条）
+
+库里原有 31 条生物公式，但**14 条是人类医学/毒理特化**（血液酒精、儿童剂量、
+孕期天数、LD50、输液速率）—— 对非人生物无用。这里补的是**任何生命形式都逃不掉**的那一层：
+
+| id | 核心 |
+|---|---|
+| `sa_volume_ratio` | 面积 ∝ M^(2/3)，**纯几何** —— 体重 ×8 → 面积 ×4 → 单位体重散热 ÷2 |
+| `fick_diffusion` | 鳃/肺/皮肤/根/细胞膜**同一个式子**；A∝M^(2/3)、δ∝M^(1/3) → 供氧 ∝ M^(1/3) |
+| `diffusion_time` | t ≈ Δx²/(2D)：1 mm → 83 min、1 m → 158 年，**立刻否掉「靠扩散供氧的巨兽」** |
+| `peclet_number` | Reynolds 数的**质量输运版**（库里只有动量版） |
+| `hardy_weinberg` | 纯概率论，跨世界通用；常算 2pq（隐性携带者 18% 而表现型只 1%） |
+| `gompertz_growth` | 比指数/logistic 更贴合**有机体**生长（logistic 偏种群） |
+| `michaelis_menten` | 一条式子覆盖酶/载体转运/Holling 摄食/受体结合/养分吸收 |
+| `q10_temperature` | 生物界**实际使用**的温度律（生态模型不用 Ea） |
+
+### 修复（真实缺陷：38 条公式不在索引里）
+
+`EM_OPTICAL_NUCLEAR_FORMULAS` 在构建 `FORMULA_INDEX` 时被**漏掉**，
+导致 38 条电学/光学/核公式不在索引里。
+
+后果很隐蔽：**单独调用能用**（`opFormula` 走的是 FORMULAS 数组），
+**但写进 pipeline / bridge 就说 `unknown_formula`** —— 同一个公式换个入口就不一样。
+这类缺陷最难发现，因为两条路径各自「看起来都正常」。
+
+修复后 `FORMULAS.length`(499) 与 `FORMULA_INDEX.size`(495) 的差只剩 4，
+即 4 个历史重复 id（被后批次覆盖，符合预期）。
+
+### 变更（help 的 domains 节补三块 + recipe 补三条配方）
+
+`domains` 节从 1500 字符扩到 8707，因为新增了材料/财政/生物三块，
+**每条都写了 scaling 律或恒等式来源**，不是只给式子。
+
+`recipe` 节新增三条**推演配方**，这是这批公式最该配套的东西：
+
+```
+  文明崩溃推演   财政六步：征管弱→赤字→基本余额为负→i>g 螺旋→铸币→改革无空间
+  巨兽可行性推演  生物六步：面积比→Fick 通量→扩散时间→Pe→饱和→生长上限
+  材料设定自洽   先定身份证（E/σy/K_IC/ρ），再让公式推后果
+```
+
+**为什么配方比公式更重要**：7 个孤立公式，用户仍然不知道怎么串；
+一条可逐步核对的链，才能推出「为什么这个政权看着还能撑、其实已在螺旋里」这类结论。
+
+### 测试
+
+97 套 / 4769 项全绿。新增 `macro-materials-bio.mjs` 61 项，
+断言方式按公式性质分三类：恒等式断言精确相等、推导式断言手算值、
+**标度律断言比值**（这是抓「次方写反」的唯一办法 —— 绝对值错了可能只是参数差异，
+比值错了说明公式形式本身写错）。
+
+例如 `fracture_kic` 用「裂纹 ×10 → 承载降 √10」来守 σ_c ∝ 1/√a，
+`sa_volume_ratio` 用「体重 ×8 → 面积 ×4」来守 M^(2/3)。
+
+两处旧断言按原意放宽（不是放宽标准）：`help-topics` 的「单节不超过旧版全文 5900」
+改成 9000 —— 分节的目的**就是**让单节写得更细，用旧全文大小当上限会把「写细」判成回归。
+## [2.3.0] — 2026-11-08
+
+这一版有两类改动：**能力扩展**（公式库与 help）与**语义修正**（让已声明的参数真正生效）。
+
+### 新增（公式库 442 → 471：工程 / 生物 / 宏观经济）
+
+原有分布很不均：physics 187、biology 26（偏医学剂量）、economics 26（几乎全是金融）、
+**engineering 只有 3 条且全是机械传动**。作者问「这塔立不立得住」「这蒸汽机出多少力」
+「这国家钱够不够花」时找不到式子。新增 29 条：
+
+| 组 | 条数 | 内容 |
+|---|---|---|
+| 工程·结构 | 7 | 抗弯截面系数、弯曲应力、自重弯矩、均布挠度、惯性矩、螺栓抗剪、角焊缝 |
+| 工程·热机 | 5 | 卡诺效率上限、蒸汽机功率、锅炉煤耗、换热器 LMTD、平壁导热 |
+| 工程·流体能量 | 6 | 管流、达西水头损失、泵功率、孔口出流、风功率、水轮机 |
+| 生物·生态 | 5 | 草地载畜量、畜禽出栏、料肉比、农田 NPP、营养级支撑 |
+| 宏观·经济 | 6 | 劳动力、总产出、储蓄率、资本产出比、货币供应、通胀压力 |
+
+每条都写明 **scaling 律**，那是最容易记错的部分：
+
+```
+  跨度的四次方    挠度翻倍 16 倍（桥梁为什么不能靠加厚延长跨度）
+  跨度的平方      自重弯矩翻倍 4 倍
+  风速的三次方    功率翻倍 8 倍（风能最反直觉的一条）
+  管径的平方      流量翻倍 4 倍；水头损失随流速平方
+  营养级十分之一  养活 1 kg 猛兽要 100 kg 草
+  保温 k 差千倍   10 cm 岩棉 ≈ 100 m 钢
+  通胀 = 货币增速 − 产出增速，0.15 对 0.03 即约 12%，三年翻倍
+```
+
+### 变更（帮助系统改成分节目录）
+
+旧的 `help` 是一整块 **5900 字符**的摊平清单，每次调用都付全量成本。
+分段统计：世界观推演占 3039 字符（52%），而算一条 `1/3` 也要读它。
+
+改为三层结构，并分两批优化：
+
+```
+  help                        → 目录（14 节，每节 key + 一句简介）
+  help, helpTopic:"world"     → 该节完整正文
+```
+
+**目录 1731 → 713 字符**：删掉 12 行「也可写: 别名…」（那些别名本来就能通过
+alias 匹配命中，列出来只为「让人看见」，AI 不需要看见清单才能查到）。
+低频节折叠成一行 `key=六字标签`，**折叠不等于藏起来** —— key 与标签仍在目录里。
+
+**分档改由实测频率决定**：6 份实践文档合计 6.1 万字符，按能力关键词统计后
+core 档定为 evaluate(250) / units(173) / solve(59) / formula(28) / world(26)。
+据此拆出独立的 `solve` 节（原来折在 evaluate 里，但信号量排第三）；
+`workspace` 信号量 28 却保留 extra ——「高频」是已发生的使用量，不是重要性。
+
+**正文放开写**：单节不再受总量约束，`linalg` 203 → 1002、`stats` 382 → 1125、
+`ode` 285 → 640。并新增 `domains` 节，把 29 条新公式的缩放律集中讲清。
+
+### 新增（符号带单位运算，与公式量纲校验）
+
+用户要求「计算要带单位和物理意义」，例子是 `V*T=L`、`V=10m/s`、`T=1s` → `L=10m`。
+
+```
+  variables:{V:'10 m/s', T:'1 s'}, expression:'V*T'   →  10 m
+  V='30 里/h', T='2 时辰'                              →  60 km
+```
+
+修复前 `V*T` 返回 `"V T"` —— 带单位的字符串值被 `bignumber()` 解析失败后整个丢弃，
+单字母随即被当成单位符号相乘。现在失败后回退到「按表达式求值」。
+
+公式的参数也带上单位了 —— 库里每个参数的描述本来就写着单位（`s → "位移 m"`），
+以前只用于显示：
+
+```
+  overrides:{s:'10 km', t:'1 h'}    → 2.77777778 m/s，并交代换算过程
+  overrides:{s:'10 kg', t:'1 s'}    → 量纲不符报错（以前静默按数字算）
+  overrides:{s:10, t:1}             → 照旧接受（不猜 10 是米还是公里）
+```
+
+**同体系折算**：输入用了里/时辰时，主结果仍保持 SI（那是公式参数的口径），
+但补一行折算 —— `60 里 / 2 时辰` 会给出 `约 30 里/时辰`，可直接心算核对。
+
+### 修复（1107 报告：单条与批量的判定不一致）
+
+- **批量放过非有限值**：单条 `1/0` 报错，批量却宣告「算式全部求值成功（4 条）」。
+  同一个式子换个调用方式就从错变成对，而批量恰恰是 AI 最常用的形态。
+  现在归入 `failed`，`ok=false` + `partial=true`。
+
+- **带单位的非有限值连单条也拦不住**：`(1/0)*2 kg` 算成 `Infinity g` 却报成功，
+  而它看起来像个正经结果。新增 `isNonFiniteValue` 下钻到 `Unit.value`。
+  `1e400`、`2^1024` 这类有限的极大数不误判。
+
+- **写出来的目标单位被输出层悄悄撤回**（本轮实测发现，报告未提）：
+  `2 m to cm` 被输出成 `2 m`。根因是单位归位只看 `unit.value`、不看意图，
+  把显式换算的目标按阶梯重挑了一遍，恰好挑回出发点。数字看着完全合法。
+  修法：表达式里有 `to` / `.to()` 时跳过归位。
+
+- **`to` 的报错从 `internal` 降级为可行动的提示**：`internal` 意为「工具自己坏了」，
+  会让人去重试。现在分成 `unit_name_conflict` 与 `not_a_quantity`。
+
+- **`evaluate` 与 `units` 两套单位表打通**：`mph`/`kn`/`rpm`/`inHg`/`cal`/`gauss`/
+  `slug`/`shortton` 等 174 个单位原先在 `units` 认得、在 `evaluate` 报未定义。
+  表不通就逼人手填系数（`10*745.7`），而手填系数正是那四起千倍级事故的共同源头。
+  注入完成照旧删掉 `createUnit`，安全边界不变。
+
+### 修复（1108 报告：已声明的参数没被读取）
+
+- **`formula` 的 `unitFrom`/`unitTo` 根本不参与计算**：schema 声明了，代码没读。
+  `overrides:{V:324375}, unitFrom:"L"` 与不传完全相同 —— 而 V 声明的是 m³，
+  结果虚高 **1000 倍**，且数字看着完全正常（11353125 分钟 ≈ 21 年，荒谬但不像错）。
+
+  修复时踩到一个更隐蔽的：`unitTo` 的兜底让它在参数**没声明单位**时也生效，
+  于是「人数 n=20」被当成 `20 L` 换算成 0.02 m³，同样虚高 1000 倍。
+  判据改成「无量纲参数一律不参与换算」。
+
+- **`formula` 不支持 `terse`**：现在 663 → 144 字符。但砍的只是标题与参数表，
+  **风险提示不砍** —— 经验式的【适用域】仍保留（那是「系数只在 1g 世界成立」的唯一提示），
+  未生效的 key 也保留（它会直接改变结论）。
+
+- **`min` 撞函数名：从「报错指引」升级为「直接可用」**。
+  mathjs 里 `min` 是函数，我们单位表里是分钟。上一轮的做法是报错并提示写 `minute`，
+  现在在**单位位置**直接认出来（`/` 后、数字后、`to` 后），
+  `1 L/min to m^3/s` 能算，而 `min(3,5)` 照旧是函数。
+
+- **`Math.xxx` 的报错指向修复方向**：以前只说「Undefined symbol Math」，
+  现在说「不支持 Math 命名空间（Math.min）+ 去掉前缀写成 min(...)」。
+  JS 习惯让人本能写 `Math.min`，这条命中率很高。
+
+- **`precision` 调低导致截断时提示**：`precision=8` 下 `10^20+1` 丢了加 1，
+  而结果看起来是个正常的大整数。现在会说明「21 位有效数字，precision=8 只能精确 8 位，
+  末 13 位可能已失真」。判据两条同时成立才提示，`p=8` 算 `1+1` 不念。
+
+- **复合单位丢单位**：mathjs 的 `toString` 会先 `simplify()`，对 `m^3 / s` 这类
+  多单位复合处理不了，只输出数值那一半 —— `1.6667e-5` 少了 `m^3 / s` 看着像 `1.67`，
+  差 10 万倍。`formatValue` 的兜底路径改为显式拼 `formatUnits()`。
+
+### 说明（本版没照做的几条）
+
+实测报告里有几条经复核是**误判或设计取舍**，没有照做，但把正确用法钉进了测试：
+
+| 报告条目 | 判定 |
+|---|---|
+| 「`workspaceVariables` 在 evaluate 完全不生效」 | **误判**。是两步式：先 `workspace action2:set` 开开关，再在 evaluate 里用，算得对。属「操作不直观」而非「功能不存在」 |
+| 「`system` 不支持 expression」 | **设计取舍**。`solve` 完全覆盖同一场景且能判退化，`system` 的定位就是吃矩阵 |
+| 「`units` 不支持复合单位」 | **不进 units**。单位表按量纲组织，加复合量纲是坏设计；`evaluate` 已能算，本版修好其中 `L/min` 的撞名问题并写进 help |
+| 「`label` 单条被忽略」 | **设计取舍**。单条只有一个值，不存在「认不出是哪个」的问题 |
+| 「`formula` 支持中文公式名/中文参数名」 | **未做**。参数名要能进 mathjs 作用域，中文会引入词法歧义；报错已点名可用参数名 |
+
+### 测试
+
+96 套 / 4680 项全绿。本版新增 4 套：
+
+| 套件 | 项数 | 守什么 |
+|---|---|---|
+| `formula-units.mjs` | 35 | 符号带单位运算、参数量纲校验、同体系折算（含反向对照） |
+| `help-topics.mjs` | 76 | 目录尺寸、可发现性、折叠节仍有 key、分档有实测依据 |
+| `engineering-formulas.mjs` | 30 | 29 条新公式手算核对 + scaling 律比值断言 |
+| `report-1108.mjs` | 41 | 本版修复项 + 两条误判的正确用法钉死 |
 ## [2.2.0] — 2026-11-03
 
 按三份实测反馈修缺陷。**重点：这一版修的都是「结果看着正常但其实是错的」**——

package/README.md

@@ -14,7 +14,7 @@ DSH（DeepSeek Harness）原生的**高精度数学计算器**插件，附带可
 | 高精度计算 | 心算 / bash 拼 | BigNumber 64 位，整数输出完整定点 |
 | 符号运算 | 查表或近似 | 求导/积分/化简/展开，符号精确 |
 | 单位换算 | 自己查系数 | 462 个单位 / 32 个量纲，温度按仿射换算 |
-| 写公式 | 每次都重推 | 442 条公式库，参数全可覆盖 |
+| 写公式 | 每次都重推 | 499 条公式库，参数全可覆盖 |
 | 改物理常数 | 无法验证 | 自动推导派生量 + 全公式偏差 |
 | 多项目数据 | 互相污染 | 工作区隔离，各持一套常数表与公式表 |
 | *（可选）* 架空世界设定 | 凭感觉 | 输入星球参数 → 出重力/气候/地形上限 |
@@ -61,7 +61,7 @@ npm run build
 
 ```powershell
 npm pack
-& "...\dsh.cmd" plugin --profile desktop add ./dsh-math-calculator-2.2.0.tgz
+& "...\dsh.cmd" plugin --profile desktop add ./dsh-math-calculator-2.3.0.tgz
 ```
 
 ---
@@ -95,7 +95,7 @@ npm pack
 
 | operation | 说明 | 示例 |
 |---|---|---|
-| `formula` | 调用 442 条公式之一 | `{"formula":"escape_velocity","overrides":{"M":5.97e24,"r":6371000},"operation":"formula"}` |
+| `formula` | 调用 499 条公式之一 | `{"formula":"escape_velocity","overrides":{"M":5.97e24,"r":6371000},"operation":"formula"}` |
 | `formulas` | 检索公式库 | `{"query":"重力","operation":"formulas"}` |
 | `world` | 星球参数 → 体感画像 | `{"world":{"planetRadius":1.4,"rotationHours":40},"operation":"world"}` |
 | `terrain` | 局地气候（配合 world） | `{"world":{"planetMass":2},"terrain":{"kind":"coastal","latitudeDeg":58},"operation":"world"}` |
@@ -197,7 +197,7 @@ npm pack
 
 ### 公式库扩展（全局包 + 工作区临时）
 
-442 条内置公式不够时，不必改代码：
+499 条内置公式不够时，不必改代码：
 
 | 级别 | 存放 | 范围 | 用法 |
 |---|---|---|---|
@@ -237,7 +237,7 @@ npm pack
 {"operation":"formulas","query":"freq:low"}    // 低频：量子/核物理/相对论/金融
 ```
 
-当前分布：高频 75 / 中频 55 / 低频 4 / 未标注 308（共 442 条）。
+当前分布：高频 123 / 中频 64 / 低频 4 / 未标注 308（共 499 条）。
 
 ### 基础运算扩展
 
@@ -352,7 +352,7 @@ npm pack
 {"operation":"quantities","constants":{"G":1.3e-10},"constantName":"G"}
 //    → 普朗克长度 ∝√G、玻尔半径与 G 无关、表面重力 ∝G
 
-// 4. 看全部 442 条公式相对现实的偏差
+// 4. 看全部 499 条公式相对现实的偏差
 {"operation":"impact","constants":{"c":1.5e8}}
 //    → 11 条受影响：质能 -75%、黑体辐射 +300%、维恩位移 -50%
 
@@ -489,7 +489,7 @@ node run-all-tests.mjs # 58 套 / 2321 项
 src/
 ├── index.ts            插件入口 + 工具定义 + operation 分发（2032 行）
 ├── ops/                8 个较大的 operation 实现
-├── formula-library*.ts 442 条公式（含势力运行、通用计算与电/光/核/化学计量批次）
+├── formula-library*.ts 499 条公式（含势力运行、通用计算、电/光/核/化学计量、工程/生物/宏观批次，以及宏观财政与材料本构批次）
 ├── constants-graph.ts  43 个常数的依赖图（23 基础旋钮 / 20 派生）
 ├── quantities.ts       量纲分析与缩放指数
 ├── feasibility.ts      11 条设定约束

package/dsh.plugin.json

@@ -1,6 +1,6 @@
 {
   "id": "dsh-math-calculator",
-  "version": "2.2.0",
+  "version": "2.3.1",
   "main": "lib/index.js",
   "description": "DSH 原生数学与架空世界推演插件：高精度计算、404 条公式、常数与星球推演。",
   "engines": {

package/lib/format.js

@@ -168,13 +168,63 @@ export function formatValue(math, value, precision, displayPrecision, respectExp
              * 的纯数字路径（那里有完整的小数位与丢精度判据），单位原样拼回。
              * 这样既不重复实现截断逻辑，也不会改变整数完整这条不变量。
              */
+            /**
+             * ★ 拆不开时必须**带着单位**退回，不能落进下面的纯数字分支。
+             *
+             * 踩过：`1 L/minute to m^3/s` 的格式化结果是 `1.6666e-5 m^3 / s`，
+             * 而 `trimUnitNumber` 的正则只认「数字 + 一个单位符号」，遇复合单位返回空串。
+             * 于是控制流落进末尾的通用分支，那里对 Unit 值走 `numeric === null`，
+             * 输出只剩 `1.666666666666666666666666666666666666666`——**单位没了**。
+             *
+             * 那个数字本身是对的（1.67e-5），但少了 `m^3 / s` 之后，读者会把它当成
+             * 「1.67」，而真实值是 0.0000167，差 10 万倍。
+             *
+             * 所以这里拆不开就直接用 mathjs 的原格式（单位一定在里面）。
+             */
             try {
-                const rendered = trimUnitNumber(math, math.format(value, { notation: 'auto', precision }), precision, displayPrecision);
+                /**
+                 * ## 复合单位必须**自己拼单位**：mathjs 的 toString 会丢它
+                 *
+                 * 实测 `1 L/minute to m^3/s`（一个 `to` 到复合单位的值）：
+                 *
+                 * ```
+                 *   math.format(v, {notation:'auto', precision})  → "1.6666666666666…"  ← 没有单位
+                 *   String(v)                                     → "1.6666666666666…"  ← 也没有
+                 *   v.formatUnits()                               → "m^3 / s"           ← 单位在这里
+                 *   v.value                                       → 0.000016666…        ← 基准值
+                 * ```
+                 *
+                 * 根因：`Unit#toString` 会先 `simplify()`，而简化对 `m^3 / s` 这类
+                 * **多单位复合**处理不了，它就只输出数值那一半。
+                 * 单单位（`kW`、`cm`）不受影响，所以只有复合单位丢 ——
+                 * 这也解释了为什么 `1 MJ/day to kW` 一直是对的。
+                 *
+                 * 后果很严重：`1.67e-5` 少了 `m^3 / s` 之后看着像「1.67」，
+                 * 而真实值是 0.0000167，**差 10 万倍**，数字本身没有任何异常。
+                 *
+                 * 所以这里显式拆成「数值 + formatUnits()」两段自己拼。
+                 */
+                const unitText = typeof value.formatUnits === 'function' ? value.formatUnits() : '';
+                const numericText = math.format(value, { notation: 'auto', precision });
+                /**
+                 * ★ 只在**确实缺单位**时才补。
+                 *
+                 * `math.format` 的行为不一致：单单位（`kW`、`cm`）本来就会带上，
+                 * 多单位复合（`m^3 / s`）才丢。一律追加会得到 `200 cm cm`。
+                 *
+                 * 判据用「结果里有没有出现单位的首个符号」—— 够简单，且不会误判：
+                 * 数值部分只有数字、指数、符号，不会包含单位名。
+                 */
+                const unitHead = unitText.split(/[\s/]/)[0];
+                const hasUnit = unitHead !== '' && numericText.includes(unitHead);
+                const raw = unitText === '' || hasUnit ? numericText : numericText + ' ' + unitText;
+                const rendered = trimUnitNumber(math, raw, precision, displayPrecision);
                 if (rendered !== '')
                     return rendered;
+                return raw;
             }
             catch {
-                /* 拆不开就退回原格式 */
+                /* 拼不出来：交给下面的兜底 */
             }
         }
     }

package/lib/formula-library-bioecon.js

@@ -0,0 +1,216 @@
+/**
+ * 公式库第十批：**生物与宏观经济**的运行式子。
+ *
+ * ## 为什么这两块一起补
+ *
+ * 生物类原有 26 条，但偏向**医学剂量**（体重给药、输液速率、血液酒精）；
+ * 经济类原有 26 条，几乎全是**金融**（复利、年金、NPV、CAGR）。
+ *
+ * 而写世界观真正会问的是另两类问题：
+ *
+ * ```
+ *   生物：这片草原能养多少头牛？这窝猪多久出栏？这个生态能支撑多少顶级捕食者？
+ *   经济：这个国家一年产多少？劳动力有多少？投资能带来多少增长？钱发多了会怎样？
+ * ```
+ *
+ * 这一批补的就是这些**运行层面**的式子，而不是又一批财务指标。
+ *
+ * ## 收录判据
+ *
+ * 同既有几批：闭合式、带量纲、有出处、不重复。经验系数标出取值范围。
+ *
+ * @module dsh-math-calculator/formula-library-bioecon
+ */
+const p = (value, desc) => ({ value, desc });
+const req = (desc) => ({ value: null, desc });
+/** 生物与宏观 11 条：生态承载、畜禽养殖、人口劳动力、货币与通胀。 */
+export const BIOECON_FORMULAS = [
+    // ══════════════════════════════════════════════════════════════════
+    // A. 生态与畜牧（土地能养多少、多久出栏、能撑几级食物链）
+    // ══════════════════════════════════════════════════════════════════
+    {
+        id: 'grazing_capacity',
+        cat: 'biology/ecology',
+        name: '草地载畜量',
+        expr: 'area * yield_per_ha * util / (intake * days)',
+        params: {
+            area: req('草场面积 ha（公顷，1 ha = 1e4 m²）'),
+            yield_per_ha: req('单位面积产草量 kg/(ha·年)（温带草原 2000~5000，荒漠 200~500）'),
+            util: p(0.5, '可利用率（放牧通常 0.4~0.6，留一半给草场恢复）'),
+            intake: req('每头每日采食量 kg/头·日（牛 40~60，羊 5~8）'),
+            days: p(365, '放牧天数（全年 365，季节性牧场按实际）'),
+        },
+        unit: '头（可承载标准畜单位）',
+        desc: '载畜量 = 面积 × 产草量 × 利用率 /(采食量 × 天数)。**「留一半」不是保守，是必需**：草地超载后会退化，第二年产量直接腰斩。这条是游牧与畜牧设定的地基。',
+        keys: ['载畜量', '草地', '放牧', 'grazing', '畜牧', '草原'],
+        frequency: 'high',
+        layer: 'phenomenological',
+        domains: ['biology', 'perception'],
+    },
+    {
+        id: 'livestock_growth_time',
+        cat: 'biology/ecology',
+        name: '畜禽出栏时间（按日增重）',
+        expr: '(target_weight - birth_weight) / adg',
+        params: {
+            target_weight: req('出栏体重 kg'),
+            birth_weight: req('初生体重 kg（仔猪 1.5，犊牛 40，羔羊 3）'),
+            adg: req('平均日增重 kg/日（肉猪 0.7~1.0，肉牛 0.8~1.2，肉鸡 0.05）'),
+        },
+        unit: '日',
+        desc: '出栏天数 = (出栏重 − 初生重)/日增重。**猪约 180 天、牛约 450 天、肉鸡约 42 天**——这个比值决定了肉价与牧场的周转率。日增重随饲料与品种变，写设定时要说清是哪种。',
+        keys: ['出栏', '日增重', 'livestock', '养殖', '肉'],
+        frequency: 'high',
+        layer: 'phenomenological',
+        domains: ['biology', 'perception'],
+    },
+    {
+        id: 'feed_conversion_ratio',
+        cat: 'biology/ecology',
+        name: '饲料转化率（料肉比）',
+        expr: 'feed_consumed / weight_gain',
+        params: {
+            feed_consumed: req('消耗饲料总量 kg（干物质）'),
+            weight_gain: req('增重 kg'),
+        },
+        unit: '无量纲（kg 料/kg 肉）',
+        desc: '料肉比 = 耗料/增重。**牛肉 6~8、猪肉 2.5~3.5、肉鸡 1.6~2.0**——这个数字直接决定「养一头牛要多少亩地」。牛高是因为反刍动物维持消耗大，不是品种差。',
+        keys: ['料肉比', '饲料转化', 'feed conversion', 'FCR'],
+        frequency: 'high',
+        layer: 'phenomenological',
+        domains: ['biology', 'perception'],
+    },
+    {
+        id: 'crop_yield_npp',
+        cat: 'biology/ecology',
+        name: '农田可收获能量（由 NPP 与收获指数）',
+        expr: 'NPP * area * hi * 1e4',
+        params: {
+            NPP: req('净初级生产力 g/(m²·年)（温带农田 800~1500，热带 2000+）'),
+            area: req('面积 ha'),
+            hi: p(0.4, '收获指数（谷物 0.4~0.5，块根作物 0.6，牧草 0.9）'),
+        },
+        unit: 'g（干物质）',
+        desc: '可收获干物质 = NPP × 面积 × 收获指数。收获指数是「长出来的东西里有多少是我们要吃的部分」——谷物的 0.4 意味着**六成生物量是秸秆**，那正是牲畜饲料与燃料的来源。',
+        keys: ['NPP', '收获指数', '产量', '作物', 'yield'],
+        frequency: 'high',
+        layer: 'phenomenological',
+        domains: ['biology', 'perception'],
+    },
+    {
+        id: 'trophic_chain_support',
+        cat: 'biology/ecology',
+        name: '营养级可支撑的生物量（十分之一定律）',
+        expr: 'primary * eff^level',
+        params: {
+            primary: req('初级生产者生物量 kg'),
+            eff: p(0.1, '营养级传递效率（陆地 0.1，海洋 0.1~0.15，以肉食为主时更低）'),
+            level: req('营养级层数（草食 1，一级肉食 2，顶级 3~4）'),
+        },
+        unit: 'kg',
+        desc: '每升一级只剩约 10%。**这就是为什么食物链一般不超过 5 级，也是为什么顶级捕食者稀少**：要养活 1 kg 狼，需要 10 kg 兔、100 kg 草。写生态设定时这条能立刻否掉「满地猛兽」的方案。',
+        keys: ['营养级', '食物链', '十分之一', 'trophic', '捕食者'],
+        frequency: 'high',
+        layer: 'derived',
+        domains: ['biology'],
+    },
+    // ══════════════════════════════════════════════════════════════════
+    // B. 人口、劳动力与产出（宏观经济的运行层）
+    // ══════════════════════════════════════════════════════════════════
+    {
+        id: 'labor_force',
+        cat: 'econ/macro',
+        name: '劳动力人口（由总人口与年龄结构）',
+        expr: 'population * participation',
+        params: {
+            population: req('总人口'),
+            participation: p(0.5, '劳动参与率（农业社会 0.45~0.55，工业化 0.55~0.65）'),
+        },
+        unit: '人',
+        desc: '劳力 = 总人口 × 参与率。**农业社会的 0.45 里含着大量童工与老人**——这个比值决定了「能征多少兵」与「能修多大工程」的上限。',
+        keys: ['劳动力', '人口', 'labor', '劳力'],
+        frequency: 'high',
+        layer: 'conventional',
+        domains: ['economics', 'society'],
+    },
+    {
+        id: 'gdp_from_labor',
+        cat: 'econ/macro',
+        name: '总产出（人力 × 人均产出）',
+        expr: 'labor * output_per_worker',
+        params: {
+            labor: req('劳动力人数'),
+            output_per_worker: req('人均年产出（按货币或实物统一定义）'),
+        },
+        unit: '与人均产出同单位（通常为年产值）',
+        desc: '产出 = 劳力 × 人均产出。最简单的宏观核算模型：**想提高总产出只有两条路——加人或提高单产**。人均产值翻倍比人口翻倍难得多，所以长期增长几乎全靠后者。',
+        keys: ['总产出', 'GDP', '宏观', '产值'],
+        frequency: 'high',
+        layer: 'conventional',
+        domains: ['economics'],
+    },
+    {
+        id: 'macro_savings_rate',
+        cat: 'econ/macro',
+        name: '储蓄率与投资能力',
+        expr: 'output * savings_rate',
+        params: {
+            output: req('总产出'),
+            savings_rate: p(0.2, '储蓄率（农业社会 0.05~0.15，高速增长期 0.25~0.40）'),
+        },
+        unit: '与总产出同单位',
+        desc: '可投资额 = 产出 × 储蓄率。**这是增长的燃料**：高速增长的经济体储蓄率往往很高（东亚 0.3+），而贫困农业社会几乎全被消费掉，攒不出投资。',
+        keys: ['储蓄', '投资', 'savings', '资本积累'],
+        frequency: 'high',
+        layer: 'conventional',
+        domains: ['economics'],
+    },
+    {
+        id: 'capital_output_ratio',
+        cat: 'econ/macro',
+        name: '资本产出比（增长的资本需求）',
+        expr: 'savings_rate / growth_rate',
+        params: {
+            savings_rate: req('储蓄率'),
+            growth_rate: req('目标增长率（小数）'),
+        },
+        unit: '无量纲（资本/产出）',
+        desc: '哈罗德-多马式的推导：所需资本产出比 = 储蓄率 / 增长率。**想每年增长 5% 而储蓄率只有 10%，需要资本产出比 2（即每 1 元产出要 2 元资本）**——如果实际要 4，那 5% 就达不到。这条能立刻判断一个增长目标现不现实。',
+        keys: ['资本产出比', '哈罗德', '增长', '投入产出'],
+        frequency: 'medium',
+        layer: 'derived',
+        domains: ['economics'],
+    },
+    {
+        id: 'money_supply_from_base',
+        cat: 'econ/macro',
+        name: '货币供应量（基础货币 × 货币乘数）',
+        expr: 'base_money * multiplier',
+        params: {
+            base_money: req('基础货币（流通现金 + 准备金）'),
+            multiplier: p(5, '货币乘数（现代银行 8~12，准备金率高时 2~4，金属货币制接近 1）'),
+        },
+        unit: '与基础货币同单位',
+        desc: '广义货币 = 基础货币 × 乘数。**金属货币制的乘数接近 1，现代部分准备金制能放大 10 倍**——这解释了为什么古代通货膨胀罕见而现代可以剧烈。架空设定里「银行体系发达程度」直接改这个乘数。',
+        keys: ['货币供应', '乘数', 'money supply', '通胀', '银行'],
+        frequency: 'high',
+        layer: 'derived',
+        domains: ['economics'],
+    },
+    {
+        id: 'inflation_from_money_growth',
+        cat: 'econ/macro',
+        name: '货币增长与产出增长之差（通胀压力）',
+        expr: '(money_growth - output_growth) * 100',
+        params: {
+            money_growth: req('货币供应年增长率（小数，如 0.15）'),
+            output_growth: req('实际产出年增长率（小数，如 0.03）'),
+        },
+        unit: '%（年通胀率近似）',
+        desc: '通胀 ≈ 货币增速 − 产出增速。**这是最实用的通胀判据**：钱印得比东西增长快多少，物价就涨多少。0.15 对 0.03 意味着约 12% 的年通胀——三年物价翻倍。写「滥发货币导致崩溃」时用这条量化。',
+        keys: ['通胀', '货币增长', 'inflation', '物价'],
+        frequency: 'high',
+        layer: 'phenomenological',
+        domains: ['economics'],
+    },
+];