dsh-math-calculator 2.3.0 → 2.3.1

package/CHANGELOG.md

@@ -2,6 +2,113 @@
 
 本项目遵循 [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）与**语义修正**（让已声明的参数真正生效）。

package/README.md

@@ -14,7 +14,7 @@ DSH（DeepSeek Harness）原生的**高精度数学计算器**插件，附带可
 | 高精度计算 | 心算 / bash 拼 | BigNumber 64 位，整数输出完整定点 |
 | 符号运算 | 查表或近似 | 求导/积分/化简/展开，符号精确 |
 | 单位换算 | 自己查系数 | 462 个单位 / 32 个量纲，温度按仿射换算 |
-| 写公式 | 每次都重推 | 471 条公式库，参数全可覆盖 |
+| 写公式 | 每次都重推 | 499 条公式库，参数全可覆盖 |
 | 改物理常数 | 无法验证 | 自动推导派生量 + 全公式偏差 |
 | 多项目数据 | 互相污染 | 工作区隔离，各持一套常数表与公式表 |
 | *（可选）* 架空世界设定 | 凭感觉 | 输入星球参数 → 出重力/气候/地形上限 |
@@ -95,7 +95,7 @@ npm pack
 
 | operation | 说明 | 示例 |
 |---|---|---|
-| `formula` | 调用 471 条公式之一 | `{"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
 
 ### 公式库扩展（全局包 + 工作区临时）
 
-471 条内置公式不够时，不必改代码：
+499 条内置公式不够时，不必改代码：
 
 | 级别 | 存放 | 范围 | 用法 |
 |---|---|---|---|
@@ -237,7 +237,7 @@ npm pack
 {"operation":"formulas","query":"freq:low"}    // 低频：量子/核物理/相对论/金融
 ```
 
-当前分布：高频 99 / 中频 60 / 低频 4 / 未标注 308（共 471 条）。
+当前分布：高频 123 / 中频 64 / 低频 4 / 未标注 308（共 499 条）。
 
 ### 基础运算扩展
 
@@ -352,7 +352,7 @@ npm pack
 {"operation":"quantities","constants":{"G":1.3e-10},"constantName":"G"}
 //    → 普朗克长度 ∝√G、玻尔半径与 G 无关、表面重力 ∝G
 
-// 4. 看全部 471 条公式相对现实的偏差
+// 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 471 条公式（含势力运行、通用计算、电/光/核/化学计量与工程/生物/宏观批次）
+├── 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.3.0",
+  "version": "2.3.1",
   "main": "lib/index.js",
   "description": "DSH 原生数学与架空世界推演插件：高精度计算、404 条公式、常数与星球推演。",
   "engines": {

package/lib/formula-library-biol-general.js

@@ -0,0 +1,175 @@
+/**
+ * 公式库第十二批：**生物学的通用公式（跨类群/跨尺度/跨世界）**。
+ *
+ * ## 为什么这一批值得单独成文件
+ *
+ * 库里原有 31 条生物公式，但**14 条是人类医学/毒理特化**：
+ * 血液酒精浓度、儿童剂量、孕期剩余天数、LD50、输液速率 ——
+ * 这些对「9000 吨工程生物」这类设定毫无用处，且参数标定与人类解剖强绑定。
+ *
+ * 这里补的是**任何生命形式都逃不掉的那一层**，判据五条全过：
+ *
+ * ```
+ *   ① 跨类群    不依赖人类解剖
+ *   ② 跨尺度    微生物到鲸，跨 10 个数量级
+ *   ③ 跨世界    常数显式暴露，形式与地球无关
+ *   ④ 有根源    物理/几何/化学/数学，不是拟合的经验系数
+ *   ⑤ 参数可得  不需要查专业表就能填
+ * ```
+ *
+ * 第 ⑤ 条最易忽略：**通用幂律数学上能算任何东西，但用户不知道指数填多少，它就成了空壳**。
+ * 所以本批凡是幂律类，指数都**内置理论值**并注明出处，不留给用户填。
+ *
+ * ## 指数从哪来（分清「几何必然」与「网络理论推论」）
+ *
+ * 这个区分很重要 —— 前者在任何世界都成立（含非碳基），后者依赖「存在分形输运系统」这个前提：
+ *
+ * ```
+ *   2/3  ← 面积/体积的纯几何（无需任何生物假设）
+ *   1/4  ← 分形输运网络的最优层级（West-Brown-Enquist）
+ *   3/4  ← 分形网络的能量通量（Damuth）
+ * ```
+ *
+ * @module dsh-math-calculator/formula-library-biol-general
+ */
+const p = (value, desc) => ({ value, desc });
+const req = (desc) => ({ value: null, desc });
+/** 生物通用公式 7 条。 */
+export const BIO_GENERAL_FORMULAS = [
+    {
+        id: 'sa_volume_ratio',
+        cat: 'biology/scaling',
+        name: '表面积-体积比（几何必然）',
+        expr: 'surface_coeff * mass^(2/3)',
+        params: {
+            mass: req('体重 kg'),
+            surface_coeff: p(0.11, '表面积系数（人体约 0.11，即 A ≈ 0.11·M^(2/3)）'),
+        },
+        unit: 'm²（表面积；除以 mass 即单位体重的表面积）',
+        desc: ['A ∝ M^(2/3)，故「面积/体积」∝ M^(−1/3)。**指数 2/3 是纯几何** ——', '任何三维物体体积 ∝ L³、表面积 ∝ L² 都成立，不需要任何生物假设，', '故跨物种、跨世界、甚至非碳基生命都适用。', '★ 一个式子回答巨兽设定的两个核心问题：体重 ×8 → 表面积 ×4 →', '单位体重散热能力 ÷2（为什么会过热）；鳃面积同样只 ×4 → 单位体重摄氧 ÷2（为什么会缺氧）。', '★ Du Bois 那种人体经验式（body_surface_area）外推到 9000 吨必错；', '但这条给的是比例关系 —— 而判断大尺度是否可行，靠的就是比例。'].join(''),
+        keys: ['表面积', '体积比', 'SA:V', '平方立方', '巨兽', '表面'],
+        frequency: 'high',
+        layer: 'fundamental',
+        domains: ['biology'],
+    },
+    {
+        id: 'fick_diffusion',
+        cat: 'biology/transport',
+        name: '菲克扩散通量（呼吸/吸收/供氧）',
+        expr: 'diffusion_coeff * area * delta_conc / thickness',
+        params: {
+            diffusion_coeff: req('扩散系数 D m²/s（O₂ 在水中 2e-9、在空气中 2e-5）'),
+            area: req('交换面积 m²'),
+            delta_conc: req('浓度差或分压差（与 D 单位配套）'),
+            thickness: req('屏障厚度 m'),
+        },
+        unit: '物质通量（与 ΔC·D 的单位配套）',
+        desc: ['J = D·A·ΔC/δ。鳃、肺、皮肤吸收、根吸水、细胞膜运输全是同一个式子。', '★ 纯物理，不含物种假设；它也是鳃氧限制理论（GOLT）的根基。', '★ 关键推论：A ∝ M^(2/3)、δ ∝ M^(1/3)，故供氧能力 ∝ M^(1/3) ——', '这就是「大生物必须换气而不能靠扩散」的定量证明。', '★ 与 diffusion_time 配合：前者给通量，后者给「扩散够不够快」。'].join(''),
+        keys: ['Fick', '菲克', '扩散', '供氧', '鳃', '通量'],
+        frequency: 'high',
+        layer: 'fundamental',
+        domains: ['biology'],
+    },
+    {
+        id: 'diffusion_time',
+        cat: 'biology/transport',
+        name: '扩散时间标度（尺寸上限的判据）',
+        expr: 'distance^2 / (2 * diffusion_coeff)',
+        params: {
+            distance: req('特征长度 m'),
+            diffusion_coeff: req('扩散系数 D m²/s（蛋白在水中约 1e-10）'),
+        },
+        unit: 's（扩散跨越该距离所需时间）',
+        desc: ['t ≈ Δx²/(2D)。布朗运动的数学结论，与是不是生物无关。', '★ 尺寸上限的直接判据（蛋白在水里 D≈1e-10）：', '1 µm → 5 ms ｜ 10 µm → 0.5 s（细胞尺度，够用）｜ 1 mm → 83 min（已不够）｜', '1 m → 158 年（不可能）—— 所以大生物必须有循环系统。', '★ 这条能立刻否掉「靠扩散供氧的巨兽」方案，且用的是硬数学而非估计。'].join(''),
+        keys: ['扩散时间', 'diffusion time', '尺寸上限', '布朗运动', '循环系统'],
+        frequency: 'high',
+        layer: 'derived',
+        domains: ['biology'],
+    },
+    {
+        id: 'peclet_number',
+        cat: 'biology/transport',
+        name: '佩克莱数（扩散 vs 对流）',
+        expr: 'length * velocity / diffusion_coeff',
+        params: {
+            length: req('特征长度 m'),
+            velocity: req('流速 m/s'),
+            diffusion_coeff: req('扩散系数 D m²/s'),
+        },
+        unit: '无量纲（Pe ≪ 1 扩散主导；Pe ≫ 1 必须对流）',
+        desc: ['Pe = L·u/D。判断「靠扩散够不够、必须不必须对流」的无量纲判据。', '★ 它是 Reynolds 数的**质量输运版**：库里已有 reynolds（动量输运），', '但缺质量输运的对应物 —— 两者成对才有意义。', '★ 实例：毛细血管里血蛋白 v≈0.3 mm/s，微秒尺度扩散位移仅 0.3 nm → 扩散主导；', '大血管 u≈m/s → 对流主导。这解释了为什么微循环靠扩散、大循环必须泵。'].join(''),
+        keys: ['Peclet', '佩克莱', '无量纲', '对流扩散', '质量输运'],
+        frequency: 'medium',
+        layer: 'derived',
+        domains: ['biology', 'physics'],
+    },
+    {
+        id: 'hardy_weinberg',
+        cat: 'biology/genetics',
+        name: '哈迪-温伯格基因型频率',
+        expr: 'p^2 + 2 * p * q + q^2',
+        params: {
+            p: req('等位基因 A 频率（0—1）'),
+            q: req('等位基因 a 频率（0—1）'),
+        },
+        unit: '无量纲（应恒为 1，用于自检与分项计算）',
+        desc: ['p² + 2pq + q² = 1。**这是纯概率论结论**（配子随机组合），不需要任何生物机制，', '所以只要「存在等位基因 + 随机交配」它就成立 —— 跨世界完全通用。', '★ 常算的是中间项 2pq：隐性基因的携带者比例。例 p=0.9、q=0.1 → 2pq=18%，', '即每 100 人里有 18 个携带者，而表现型只有 q²=1%。', '★ 若结果明显不等于 1，说明 p+q≠1，输入有误 —— 这是一条自带自检的公式。'].join(''),
+        keys: ['Hardy-Weinberg', '哈迪温伯格', '基因频率', '遗传平衡', '携带者'],
+        frequency: 'high',
+        layer: 'fundamental',
+        domains: ['biology'],
+    },
+    {
+        id: 'gompertz_growth',
+        cat: 'biology/growth',
+        name: 'Gompertz 生长曲线',
+        expr: 'asymptote * exp(-shape * exp(-rate * t))',
+        params: {
+            asymptote: req('渐近上限（最终体型/规模）'),
+            shape: req('形状参数 b（决定曲线起点位置）'),
+            rate: req('早期增长率 c（1/时间）'),
+            t: req('时间'),
+        },
+        unit: '与 asymptote 同单位',
+        desc: ['N(t) = a·exp(−b·exp(−c·t))。描述「越大长得越慢」的通用形式。', '★ 比指数与 logistic 更贴合**有机体**的真实生长（logistic 更偏种群）：', '跨动物、植物、细菌、肿瘤都拟合良好。', '★ 三个参数的物理含义：a = 渐近上限，c = 早期增长率，b = 曲线形状与起点。', '★ 跨世界照用：形式通用，系数按对象标定。'].join(''),
+        keys: ['Gompertz', '冈珀茨', '生长曲线', '渐近', '生长'],
+        frequency: 'high',
+        layer: 'phenomenological',
+        domains: ['biology'],
+    },
+    {
+        id: 'michaelis_menten',
+        cat: 'biology/kinetics',
+        name: '米氏方程（饱和动力学）',
+        expr: 'vmax * substrate / (km + substrate)',
+        params: {
+            vmax: req('最大速率 Vmax'),
+            km: req('半饱和常数 Km（速率达 Vmax/2 时的底物浓度）'),
+            substrate: req('底物浓度 S'),
+        },
+        unit: '与 vmax 同单位（速率）',
+        desc: ['v = Vmax·S/(Km+S)。一条式子覆盖五种现象：酶动力学、载体转运、', 'Holling II 型摄食功能反应、受体结合、土壤养分吸收 —— 因为它们都是同一个', '「先线性、再饱和」的双曲线形式。', '★ 结构性质：S ≪ Km 时 v ≈ (Vmax/Km)·S（线性）；S ≫ Km 时 v → Vmax（饱和）。', '这个「处理能力有限必然饱和」是所有系统的共同特征。', '★ 与 q10_temperature 配合：前者管底物，后者管温度上限。'].join(''),
+        keys: ['Michaelis', '米氏', 'Km', 'Vmax', '酶', '饱和', 'Holling'],
+        frequency: 'high',
+        layer: 'fundamental',
+        domains: ['biology'],
+    },
+    {
+        id: 'q10_temperature',
+        cat: 'biology/kinetics',
+        name: 'Q₁₀ 温度系数（生物代谢的温度律）',
+        expr: 'rate_ref * q10^((temperature - temp_ref) / 10)',
+        params: {
+            rate_ref: req('参考温度下的速率'),
+            q10: p(2.5, 'Q₁₀ 系数（跨类群常见 2~3）'),
+            temperature: req('实际温度'),
+            temp_ref: req('参考温度'),
+        },
+        unit: '与 rate_ref 同单位（速率）',
+        desc: ['R(T) = R₀·Q₁₀^((T−T₀)/10)。温度每升高 10 度，速率变为 Q₁₀ 倍。', '★ 这是**生物系统实际使用的温度律**：生态模型几乎都用它，', '因为不需要活化能 Ea（那要查专业表）。', '★ 与库里已有的 arrhenius（化学式，要 Ea）并存而不替代：', '化学场景用 Arrhenius，生物场景用 Q₁₀ —— 同一条规律的两种学科参数化。', '★ Q₁₀ 跨类群很稳（2~3），这是它比 Arrhenius 更好用的原因。'].join(''),
+        keys: ['Q10', '温度系数', '代谢温度', '温度律', '生态模型'],
+        frequency: 'high',
+        layer: 'phenomenological',
+        domains: ['biology'],
+    },
+];