← Archive
lm-005217 · 2026-11

GMDIS-11|規則勢態勢換位與物理反事實:狀態—規則—邊界—觀察的四層干預與可辨識世界線

下載 MD 檔 ⬇
title
GMDIS-11|規則勢態勢換位與物理反事實:狀態—規則—邊界—觀察的四層干預與可辨識世界線
english_title
Rule-Shi Position Shifts and Physical Counterfactuals: Four-Level Interventions over State, Rule, Boundary, and Observation, and Identifiable Future Trajectories
series
Global Multiscale Dynamics and Informational Sufficiency Theory (GMDIS)
paper
11 / 14
version
v0.1
date
2026-10-17
timezone
Asia/Taipei
author
Neo.K
research_collaboration
Aletheia / GPT-6
institution
EveMissLab/一言諾科技有限公司
language
zh-TW
status
Detailed research manuscript; model-conditional analytical propositions, synthetic reproducibility checks; real-world experimental validation pending
canonical_format
UTF-8 Markdown; mathematics delimited only by $...$ and $$...$$
predecessors
["GMDIS-01", "GMDIS-02", "GMDIS-03", "GMDIS-04", "GMDIS-05", "GMDIS-06", "GMDIS-07", "GMDIS-08", "GMDIS-09", "GMDIS-10"]
next_paper
GMDIS-12|節拍器矩陣的多尺度基準世界
license_note
研究原稿;公開版本及授權條件由權利人另行決定

GMDIS-11|規則勢態勢換位與物理反事實

狀態—規則—邊界—觀察的四層干預、可比較未來與可辨識實驗

Rule-Shi Position Shifts and Physical Counterfactuals: Four-Level Interventions over State, Rule, Boundary, and Observation, and Identifiable Futures

系列: 全域多尺度動力學與資訊充分性理論(GMDIS),第 11/14 篇
作者: Neo.K;AI 研究協作: Aletheia/GPT-6
機構: EveMissLab/一言諾科技有限公司
版本: v0.1;2026 年 10 月 17 日(Asia/Taipei)
主要學科: 非線性動力學、數學物理、因果干預、偏微分方程、結構辨識、反事實實驗、全域規則勢、資訊自由度。
學術聲明: 本文直接證明的命題僅在所列有限維模型或明確的線性 PDE 假設下成立;合成數值程序不是實物測量;GRST/SIO 的語義接縫是提出中的研究方法,不等於已證成新的基本物理定律。本篇不宣稱能人為更換自然宇宙的基本法則,更不宣稱「規則勢」本身是能量、牛頓力或額外物理場。

日期: 2026-10-17


摘要

GMDIS-01 至 GMDIS-10 已從物理振盪系統的相位與吸引子投影不足問題出發,逐步建立全域多尺度世界的狀態分層、GRST—SIO 接縫、廣義張力的物理語義校準、多場泛函分析、固體—聲學—熱學耦合、持續驅動與穩定性、隱藏自由度的記憶核、24/72 範式的有型別模型分類,以及多不變量觀測充分性稽核(MIOSA)。但至第十篇為止,仍存在一個方法論缺口:即使已經觀察到兩個世界的長期運動不同,差異究竟來自初始狀態、物理生成規則、邊界條件,還是觀測方式本身? 若四者不分,則「規則勢換位」容易退化成事後描述,也無法區分預測增益、干預因果增益與模型選取偏差。

本篇依循《GRST-07:因果勢與雙層反事實》的四種干預 doX,doR,doB,doOdo_X,do_R,do_B,do_O,重新建立物理世界的可執行語義:以同一可比較世界狀態或明示的狀態橋接,凍結應保留的物理單位、時間窗、外部驅動及觀測任務,對規則、邊界、狀態、感測依次做控制干預。對規則勢換位,本篇不將它定義為神祕的「勢量增加」,而給出一組具型別的可行未來剖面:局部線性穩定、有限時間擾動放大、可達性與作用傳播、能量交換與耗散、分岔位置、干預敏感度、跨尺度觀測充分性。只有在可比性契約成立、且至少一項預註冊的未來相關量呈現可重現差異時,才把該結果記為在指定任務與語境下的「態勢換位」。

本篇建立四組可獨立核對的解析案例。第一,對受控的超臨界叉式分岔正規形,同一平衡點在參數由負變正時,未受擾的觀察可完全相同,但局部穩定性與擾動後未來發生改變。第二,對同一根熱桿、同一初始溫度分布與同一擴散率,只改絕熱 Neumann 或固定端溫 Dirichlet 邊界,即可證明長期總熱含量與空間形狀顯著不同。第三,兩個相同機械振盪器在同一初始位移、速度與空間位置下,僅變更它們間的彈性耦合,衝量造成的遠端響應由嚴格零變為非零。第四,利用確定性線性流,直接推導規則與狀態的非加性交互項,說明只有四因子或多因子的設計才能合理分離共同作用。

本文同時給出「同樣被動資料,干預後可能不同」的參數不可辨識反例、不同模型無法直接比較狀態時的橋接條件、物理守恆與能量帳本要求,以及兩種避免誤判的負例:單純更換感測器不必改變物理規則勢;在實際研究任務與所有合法干預下等價的規則,不應硬被判為不同的可觀測態勢。 以此形成 GMDIS-12/13 的物理干預矩陣、合成重現測試和對現有因果推論、系統辨識、混合動力學方法的公平競爭條件。

關鍵詞: physical counterfactual、rule-level intervention、boundary intervention、observational intervention、GRST、SIO、Shi position shift、intervention identifiability、future-structure profile、bifurcation、heat boundary conditions、oscillator networks、state–rule interaction、MIOSA、Rashomon rule family、physical conservation audit。


0. 研究責任與命題等級:不把新名詞當結果

0.1 前十篇的直接繼承

本篇接受以下已在原稿中形成的研究契約:

  1. GMDIS-01/02: 只呈現相位差、吸引子形狀或同步程度,不必然足以識別完整物理機制;但在適當可觀測與重建條件下,低維或延遲表示可以足夠。
  2. GMDIS-03: SIO 的硬約束自由度陰影與物理交互作用的動態影響陰影不能混同;即使局部合法自由度數目不變,未來穩定性及耦合仍可不同。
  3. GMDIS-04: 廣義張力是概念性關係語法,不能不經單位校準直接當作牛頓力;物理的功、功率、能量、應力、熵和規則勢均需分開。
  4. GMDIS-05/06: 場的方程與界面必須滿足其模型級別的適定性、傳播性和能量帳本;反事實不是免除守恆與邊界相容要求的許可。
  5. GMDIS-07: 穩定運動可以是持續驅動與耗散平衡;同步不等於零能流,模態穩定不保證唯一吸引域。
  6. GMDIS-08: 被粗粒化的自由度可藉由記憶項及初始隱藏態繼續影響可見系統;需區分新物理介入與模型換表示。
  7. GMDIS-09: 24/72 為帶語境的計算/模型表示分類,不是自然界 72 條作用律或 72 種力量。
  8. GMDIS-10: MIOSA 為指定任務、成本、誤差與合法干預集下的觀測充分性稽核; XX (拒單測)是條件性判定,不是「多變數便不可測」。

因此,本篇不是另一篇同步現象回顧,而是將「世界可行未來的結構性非均勻」轉為可以被干預設計檢查的跨規則比較問題。

0.2 D、M、C、H、N 的證據記號

  • D-11-xx|定義: 在本研究中明確約定的數學或型別語義,不需以實驗冒充其已存在。
  • M-11-xx|模型內命題: 在全部前提成立時可由既有數學推導的結果;不能自動升格為普遍物理新定律。
  • C-11-xx|構造性反例: 用已知微分方程、偏微分方程或受控機械模型展示某個不充分敘述的失敗。
  • H-11-xx|待驗證研究假說: 需要實物、不同參數族或跨資料集外部測試。
  • N-11-xx|合成數值一致性檢查: 附錄可執行程式驗證解析預測;不是新實驗材料。

任何數值、證明或文獻若不滿足上述對應,不得混入更高證據層級。

0.3 概念中最危險的六種混淆

容易混淆的對象 必須維持的區分
規則 RR 與狀態 xx 更改生成法則不是單純替換一個觀測值;但若把參數納入擴展狀態,也可以在統一世界模型中表示
規則 RR 與邊界 BB 物理上都可能由方程與邊界條件共同定義;本篇採實驗設計的型別分工,不宣稱存在普遍唯一劃線
邊界 BB 與幾何排列 GG 幾何相同不表示接觸、固定方式、材料剛度與介質邊界相同
被動觀測與實際干預 看見兩個系統不同,不證明是哪個物理因素造成它們不同
觀察方式 OO 與世界物理 換相機或讀數通常只改觀察;但加裝有反作用的感測器會改物理環境
物理態勢換位與計算分類換位 C ⁣→ ⁣DC\!\to\!D 的標籤改變可只來自離散化,不代表世界本身已失穩

1. 從 GRST-07 回收四層干預,但改用可實驗的物理語義

1.1 GRST-07 的原始研究問題

《GRST-07》將世界內的狀態干預與世界生成規則干預分開,並額外區分邊界與觀察。四類寫為:

I={doX,doR,doB,doO}.\mathcal I=\{do_X,do_R,do_B,do_O\}.

其核心的哲學洞察並不是「人類能讓宇宙法則改寫」,而是:若只改狀態和若改規則可以得到不同的未來結構,那麼只使用一個普通狀態差分不足以描述全部反事實差異。以節拍器而言,移動擺錘初態、修改彈簧常數、固定或釋放平台、換攝影讀出,分別是不同研究干預層。

1.2 D-11-01:四層的操作型別

狀態干預 doXdo_X: 在同一生成模型下改變初始位移、速度、擒縱內部狀態或已知介質擾動;需要報告狀態準備所耗能量與合法性。

規則干預 doRdo_R: 改變有效模型的演化參數或機制,如彈簧常數、阻尼、接觸力律、擒縱驅動策略、作用網絡的耦合律、確定與隨機模型假設。這是實驗模型級的改變,不表示宇宙基本力改寫。

邊界干預 doBdo_B: 修改周圍材料、固定方式、聲腔反射條件、溫度梯度、容器幾何或與外部交換的功率埠。這些也能被包含在更大形式的 RR 內;本篇將其另列,是為了可控制、可歸因的實驗分解。

觀察干預 doOdo_O: 改變觀察通道、相機角度、感測頻率、測量誤差、分析尺度與壓縮投影。理想的被動觀察干預不改物理軌道;若感測負載有物理作用,必須把它額外記為 doXdo_X 、 doRdo_R 或 doBdo_B,不許悄悄隱藏。

1.3 D-11-02:實驗幹預不是哲學上的「改真理」

要比較兩條物理軌跡,需將外界物理定律(例如牛頓定律、流體動量守恆)維持為模型的共同假設,只在可操控機構、材料、幾何、邊界與輸入上干預。對某種真實但不可操控的「自然常數改變」只能提出模型內反事實,不可聲稱完成實驗介入。例如「若重力加速度不同,鐘擺週期如何變化」是明確的模型比較,但不等於在實驗室直接改變基本引力常數。

1.4 干預必須記錄資源帳本

狀態被瞬時重置時,動量或機械能可能改變;更換彈簧或更換邊界亦可能造成外部施力、重新裝配與初態準備。合法反事實應明列:

ΔE=∫0TPin(t) dt−∫0TPout(t) dt\Delta E=\int_0^T P_{\mathrm{in}}(t)\,dt-\int_0^T P_{\mathrm{out}}(t)\,dt

此式是所選系統控制體的能量帳本;若還有內能、熱、衝擊儲能或其他介質,必須增加其對應項。將更換設備的物理代價忽略,會讓所謂「只改規則」其實暗中同時改了狀態與能量。


2. 世界模型、共同狀態橋接與可比較的未來

2.1 D-11-03:一般世界模型

令待研究的有效物理世界為:

W=(X,R,B,U,O,Γ).W=(\mathcal X,R,B,\mathcal U,\mathcal O,\Gamma).

其中:

  • X\mathcal X:明定單位、可行域及約束的完整狀態空間;在多場情形可為 Hilbert 或 Banach 積空間,在混合機械系統可含離散模式。
  • RR:狀態演化的有效生成規則或算子,包括擒縱、材料與非線性機制;需標明是否自治、時變、隨機或含事件重置。
  • BB:容器及外界邊界條件;不可只給圖形尺寸而省略牽引、熱流及功率埠。
  • U\mathcal U:合法外部輸入與干預的集合;不可為了找到差異而使用實驗上做不到的無窮脈衝。
  • O\mathcal O:合法量測/讀出族,含校準、時間與解析度。
  • Γ\Gamma:研究語境,記錄觀察尺度、時間窗、誤差限、資源限制、初態分布和決策/預測任務。

在確定性且適定的情況下,世界線可寫為:

ΦR,Bu(t,x0).\Phi_{R,B}^{u}(t,x_0).

這裡 uu 表示外部輸入或試驗協定。若是隨機系統,應改為一個有條件的軌跡機率律或可測轉移核;若包含碰撞、滯回及離散模式,應附守衛、重置與事件時間。本文不強迫不同動力學都由單一連續微分方程表示。

2.2 D-11-04:同一完整狀態的最低條件

比較 R0R_0 和 R1R_1 時,不可隨意宣布 x0x_0 「相同」。最低條件是:

x0∈XR0,B0∩XR1,B1.x_0\in\mathcal X_{R_0,B_0}\cap\mathcal X_{R_1,B_1}.

而且每一分量必須共享相同的單位與參考架構。同為「擺角 15∘15^\circ 」並不代表擒縱儲能或空氣壓力場相同。同為「同一張照片」亦不代表具有相同的速度及隱藏邊界態。

如果兩個模型根本不具相同狀態空間,例如一個是有限維振子、一個是聲場與固體耦合 PDE,應使用兩個合法嵌入:

j0:Z→X0,j1:Z→X1.j_0:\mathcal Z\to\mathcal X_0, \qquad j_1:\mathcal Z\to\mathcal X_1.

比較的共同條件是 z∈Zz\in\mathcal Z。不在橋接 j0,j1j_0,j_1 中說明的隱藏態,不能假裝已控制一致。

2.3 D-11-05:共同輸出型別與反事實核

讓各模型結果都落在可比的任務輸出空間 Y\mathcal Y:

Yi(t;z,u)=hi ⁣(ΦRi,Biu(t,ji(z))),i=0,1.Y_i(t;z,u)=h_i\!\left(\Phi_{R_i,B_i}^{u}(t,j_i(z))\right),\qquad i=0,1.

hih_i 必須單位相同,且物理意義可以比較。允許不同內部座標,但不能把一邊的攝影像素當作另一邊的焦耳或牛頓。若任務比較未來事件、穩定性或風險分布, Y\mathcal Y 也可為向量或分布空間。跨世界比較屬於有型別、受限語境的比較;不是聲稱世界本體彼此同一。

2.4 三種不相等的「相同」

  • 幾何相同: 外觀與座標一致;可以仍有不同應力、接觸剛度或介質作用。
  • 完整狀態相同: 在指定模型類中,全部狀態分量及隱藏初態一致。
  • 被動觀察相同: 同一個觀測算子無法辨認差異;未必表示內部狀態或規則相同。

尤其把「同樣相位」直接當「同樣完整狀態」是 GMDIS-02 已指出的典型錯誤。

2.5 四種可比較的反事實成對條件

實驗條件 必須固定 只改變 首先評估
doXdo_X R,B,O,uR,B,O,u 合法初態或外部擾動 方向敏感度、穩定域、回復時間
doRdo_R x0,B,O,ux_0,B,O,u 有效生成參數/機制 模態、分岔、動力影響陰影、未來分布
doBdo_B x0,R,O,ux_0,R,O,u 物理邊界/介面 反射、熱交換、可達性、端口功率
doOdo_O x0,R,B,ux_0,R,B,u 被動感測與表示 觀測充分性、噪聲與識別能力

當「只改變」無法物理實現時,必須同時保留干預帶來的狀態準備改變;這在物理研究中是不可消除的設計問題,不是記帳上的失敗。


3. GRST 勢、局部態勢與「換位」的研究級定義

3.1 GRST-12 的母鏈保持不動

本篇沿用原始順序:

T⟶R⟶FR⟶ΘR⟶ΣR⟶σR(x)⟶ΓR,x⟶O.T\longrightarrow R\longrightarrow\mathcal F_R\longrightarrow\Theta_R\longrightarrow\Sigma_R\longrightarrow\sigma_R(x)\longrightarrow\Gamma_{R,x}\longrightarrow O.
  • 可行自由集合 FR\mathcal F_R 是仍然合法的方向或狀態配置。
  • 勢賦形 ΘR\Theta_R 描述這些合法方向在成本、放大、穩定、延遲與相互耦合上的非均勻。
  • 全域規則勢 ΣR\Sigma_R 是研究規則誘導的未來結構,而不是新增一個可以與機械功率相加的勢能。
  • 局部態勢 σR(x)\sigma_R(x) 是指定尺度/時間窗下的未來結構位置,而非世界空間裡某個具有 SI 單位的物理座標。

3.2 D-11-06:物理可操作的未來態勢剖面

為避免在未驗證之前宣稱存在一個普遍、完整的 ΣR\Sigma_R 數值,本篇只定義選定任務族 Q\mathcal Q 與時間窗 [0,T][0,T] 的剖面:

SΓ(W,x)=(AT,ΛT,BT,ET,CT,HT).S_{\Gamma}(W,x)=\bigl(\mathcal A_T,\Lambda_T,\mathcal B_T,\mathcal E_T,\mathcal C_T,\mathcal H_T\bigr).

六個分量分別為:

  1. AT\mathcal A_T |可達性: 在指定合法控制、能量預算、時間窗和邊界下,可到達哪些目標集合;單純可達集合大小不足以代表作用成本。
  2. ΛT\Lambda_T |擾動放大/恢復性: 局部線性化最大增長率、有限時擴大因子、Floquet 乘數、對穩定集合的回復時間等,須標註所用範數。
  3. BT\mathcal B_T |吸引域/分岔: 指定模態的吸引盆、存在條件及臨界參數;有限時間資料不一定足以確認真正的漸近吸引域。
  4. ET\mathcal E_T |物理帳本: 機械儲能、聲能、熱量、外部功率、衝擊耗散的分離記錄;本篇不將這些組合成自稱「總勢能」的量。
  5. CT\mathcal C_T |耦合與因果影響: 在指定合法干預下,某單位或介質的擾動如何影響其他單位及其延遲、方向與恢復性。
  6. HT\mathcal H_T |表示/觀察充分性: 依照 MIOSA 記錄原觀察與增量觀察是否足以區分候選模型,這是我們知道多少,不是物理世界自身的力量。

這些型別可以同時存在,不能在未定義物理單位、正規化、目標和誤差容忍之前直接加權累計為任意純量。

3.3 D-11-07:有限任務的態勢換位證書

給定兩個完成合法對齊的世界方案 W0,W1W_0,W_1 和同一共同初態 zz。先預註冊一族可比較的未來任務 Q\mathcal Q 、指標 JqJ_q 、誤差容限 εq\varepsilon_q 與測試時間窗。若存在 q∈Qq\in\mathcal Q 使:

∣Jq(W1,z)−Jq(W0,z)∣>εq,\left|J_q(W_1,z)-J_q(W_0,z)\right|>\varepsilon_q,

且差異超過統計與數值不確定性(對解析模型可由恆等式直接核對),則給予:

Shift⁡Γ(W0,W1;z)=1.\operatorname{Shift}_{\Gamma}(W_0,W_1;z)=1.

這是一種研究語境下的差異證書,不是普遍斷言「兩個世界的本體必然不相等」。若差異不超過門檻,標記為未發現或任務等價,而非武斷宣稱全世界完全同勢。

3.4 D-11-08:態勢距離不能比對無關數值

若要計算距離,須先使同名指標具有共同定義、尺度與 SI 單位。例如:

DΓ2=∑q∈Qwq ⁣(Jq(W1,z)−Jq(W0,z)sq)2,D_{\Gamma}^{2}=\sum_{q\in\mathcal Q}w_q\!\left(\frac{J_q(W_1,z)-J_q(W_0,z)}{s_q}\right)^2,

其中 wq≥0w_q\geq0 、 sq>0s_q>0 均需預先固定。這是任務感知特徵距離,不會因形式長得像歐氏範數便自動成為所有 GRST 勢的普遍度量。若跨世界狀態空間沒有合法的輸出橋接,距離應標為未定義。

3.5 換位不必伴隨已發生的軌跡變化

最重要的情形是:同一平衡點在兩套規則下仍保持原位,甚至在被動觀測下永遠為常數,但擾動後穩定性反轉。這時「實際觀察的歷史」可以完全相同,卻仍有不同的未來可能性與脆弱性。GRST 要求比較的是後者,但要產生實證價值,必須找到能夠區分這兩種世界的合法干預。


4. 第一個核心命題:同一平衡點、不同規則、不同局部未來

4.1 C-11-01:受控的超臨界叉式分岔模型

取一維無量綱有效動力系統:

x˙=μx−x3,x(0)=x0.\dot x=\mu x-x^3, \qquad x(0)=x_0.

其中 μ\mu 是可調的有效參數;方程是正常形 toy model,不是節拍器的完整運動方程。此模型在 μ<0\mu<0 與 μ>0\mu>0 具有不同的平衡點穩定結構。固定幾何、座標、狀態空間與未干預初態:

x0=0.x_0=0.

對所有 μ\mu,皆有 x(t)=0x(t)=0,因此只觀測自然放置的零軌道,無法看出規則差異。

4.2 M-11-01:零軌道相同但穩定性換位

命題: 上述模型在 μ<0\mu<0 時原點局部漸近穩定,在 μ>0\mu>0 時原點不穩定。對 μ>0\mu>0 還有兩個局部漸近穩定平衡點:

x±=±μ.x_{\pm}=\pm\sqrt{\mu}.

證明: 向量場 fμ(x)=μx−x3f_\mu(x)=\mu x-x^3。原點之線性化導數 fμ′(0)=μf_\mu'(0)=\mu:負值表示局部指數收縮,正值表示局部指數擴張。非零平衡點滿足 x2=μx^2=\mu (只在正參數下存在),其導數 fμ′(x±)=μ−3μ=−2μ<0f_\mu'(x_{\pm})=\mu-3\mu=-2\mu<0,故各自局部漸近穩定。當 μ=0\mu=0,線性化導數為零,不能用雙曲平衡點論證;此時 x˙=−x3\dot x=-x^3 仍為漸近穩定但非指數穩定。證畢。

這是一個精確的 GRST 物理語義候選:當前位置沒有改變,合法自由維度仍為一維,甚至未受擾動的實際歷史可以完全一致,但未來擾動的放大、恢復與分岔結構已改變。

4.3 M-11-02:正初態的精確解與反事實差

對 x0>0x_0>0 、 μ≠0\mu\ne0,令 v=x2v=x^2,則:

v˙=2μv−2v2.\dot v=2\mu v-2v^2.

解得:

xμ(t;x0)=x0eμt1+x02μ(e2μt−1).x_\mu(t;x_0)=\frac{x_0e^{\mu t}}{\sqrt{1+\frac{x_0^2}{\mu}(e^{2\mu t}-1)}}.

此式適用於 μ<0\mu<0 與 μ>0\mu>0,分母在 t≥0t\ge0 的本例範圍保持正值。當 μ=0\mu=0,解為:

x0(t;x0)=x01+2x02t.x_0(t;x_0)=\frac{x_0}{\sqrt{1+2x_0^2t}}.

因此以相同的小初始偏移 x0=δx_0=\delta 做 doXdo_X,再比較 μ−<0\mu_-<0 與 μ+>0\mu_+>0 的 doRdo_R,兩個世界立即產生不同的擾動回應。早期當 δ\delta 很小且時間不太長,增長率近似由 eμte^{\mu t} 決定;當 μ>0\mu>0 時非線性 −x3-x^3 阻止無限增長,後來靠近 +μ+\sqrt{\mu},而不是無限制爆炸。

4.4 不可把本例直接當作新的「力的萬有定律」

此正常形已是成熟的非線性動力學模型。GMDIS 的候選新意只在研究程序:使用共同狀態、規則、邊界和觀察的型別化干預,判定哪些現象屬於物理態勢差異、哪些只是讀出差異。若這套程序沒有增加設計或識別價值,應由既有分岔與系統辨識方法吸收。

4.5 H-11-01:實物中能否識別對應的局部換位?

在真實節拍器世界,改變底座阻尼或耦合強度可能改變不同同步模式的橫向恢復性,某些模態可能失穩。不能事先假設一定具有叉式分岔;本例只提示一類可測的干預差:固定觀察到的同步狀態,施加可重複的微小偏移,看它是恢復、分群、轉模態還是停止。


5. 第二個核心命題:同一溫度世界,只換邊界,改變整個可行未來

5.1 C-11-02:一根熱桿的兩種世界

令材料佔據同一區間 0<x<L0<x<L,導熱擴散率 κ>0\kappa>0,取線性熱方程:

∂tu=κ∂xxu,0<x<L.\partial_t u=\kappa\partial_{xx}u, \qquad 0<x<L.

固定同一初始溫度相對基準場:

u(x,0)=sin⁡2(πx/L).u(x,0)=\sin^2(\pi x/L).

考慮兩個實驗世界:

  • BNB_N:兩端絕熱, ux(0,t)=ux(L,t)=0u_x(0,t)=u_x(L,t)=0 (Neumann)。
  • BDB_D:兩端固定為基準溫度, u(0,t)=u(L,t)=0u(0,t)=u(L,t)=0 (Dirichlet)。

初始函數在 x=0,Lx=0,L 既為零,其導數亦為零,因此對兩種標準邊界均具有初始相容性。材料、局部熱傳規則、幾何和初始分布都一樣,只有邊界交換條件不同。若把邊界規則全數納入生成元,它們自然具有不同的生成元定義域;但為了物理實驗可歸因,本篇刻意把差異記為 doBdo_B。

5.2 M-11-03:絕熱世界的平均溫度保持

定義總相對熱含量的無量綱代表量:

M(t)=1L∫0Lu(x,t) dx.M(t)=\frac{1}{L}\int_0^L u(x,t)\,dx.

對足夠平滑解,在上述方程下:

M˙(t)=κL[ux(L,t)−ux(0,t)].\dot M(t)=\frac{\kappa}{L}\left[u_x(L,t)-u_x(0,t)\right].

在 Neumann 邊界條件,右側嚴格為零,因而:

MN(t)=MN(0)=12.M_N(t)=M_N(0)=\frac12.

而且可以直接寫出完整解析解:

uN(x,t)=12−12e−4π2κt/L2cos⁡(2πx/L).u_N(x,t)=\frac12-\frac12e^{-4\pi^2\kappa t/L^2}\cos(2\pi x/L).

因此時間趨向無窮大時,整根桿的溫度趨向常數 1/21/2。這裡的「保持」不是 GRST 勢的保存,而是由絕熱熱通量邊界導致的總熱量平衡。

5.3 M-11-04:固定端溫世界的總熱量趨零

在 Dirichlet 世界,對 u(0,t)=u(L,t)=0u(0,t)=u(L,t)=0 的熱方程,使用正弦本徵函數分解:

uD(x,t)=∑n=1∞bnsin⁡(nπx/L)e−κn2π2t/L2,u_D(x,t)=\sum_{n=1}^{\infty}b_n\sin(n\pi x/L)e^{-\kappa n^2\pi^2t/L^2},

其中:

bn=2L∫0Lsin⁡2(πx/L)sin⁡(nπx/L) dx.b_n=\frac2L\int_0^L\sin^2(\pi x/L)\sin(n\pi x/L)\,dx.

計算可得:

bn={−8πn(n2−4),n 為奇數,0,n 為偶數.b_n=\begin{cases} -\dfrac{8}{\pi n(n^2-4)}, & n\ \text{為奇數},\\ 0, & n\ \text{為偶數}. \end{cases}

對每個 t>0t>0,上述展開快速收斂;最大本徵模態的衰減至少受第一個正特徵值 κπ2/L2\kappa\pi^2/L^2 控制,所以:

lim⁡t→∞∥uD(⋅,t)∥L2(0,L)=0,lim⁡t→∞MD(t)=0.\lim_{t\to\infty}\|u_D(\cdot,t)\|_{L^2(0,L)}=0, \qquad \lim_{t\to\infty}M_D(t)=0.

同時:

MD(t)=∑n≥1nodd2bnnπe−κn2π2t/L2.M_D(t)=\sum_{\substack{n\ge1\\n\,\mathrm{odd}}}\frac{2b_n}{n\pi}e^{-\kappa n^2\pi^2t/L^2}.

所以在長時間下,兩個世界的平均溫度終值差是 1/21/2。若溫度尺度為 ΔT\Delta T 、長度為 LL 、密度和比熱已知,乘上熱容對應的權重才能還原總熱能;不可把無量綱 MM 直接稱為焦耳。

5.4 為什麼這不是「換個感測器」?

兩種邊界不是感測投影,而是實際改變與外界交換熱量的物理條件。相同的熱方程在兩種邊界下生成不同的半群;在無量綱表示上:

TN(t)≠TD(t).T_N(t)\ne T_D(t).

即便把所有時刻的中間位置溫度都測得非常精確,也不能用觀察技巧讓兩種邊界下的未來變成相同。這是最乾淨的物理邊界造成未來結構換位的解析例子。

5.5 SIO 的兩種「陰影」在本例怎麼區分?

  1. 硬約束陰影: 固定端溫要求邊界值不能自由變動,兩端溫度自由受到硬限制;此時候選解空間與相容初態空間的性質改變。
  2. 動態影響陰影: 在桿中間沒有直接碰到外界散熱端的點,其未來溫度仍會因傳熱而受邊界改變影響;這種傳播並非可用「幾何接觸範圍」直接替代。

熱方程在數學上具有任意正時間的非零核支撐,這是擴散 PDE 的理想化性質,不應直接宣稱熱擾動在真實物理中以無限速度傳播;需在所要求的時間/空間尺度下改用具有適當微觀或有限傳播描述的物理模型。

5.6 物理反事實不保證「合法空間數目相同」

本例恰恰顯示某些 doBdo_B 對 admissibility domain 的限制會改變。不能為了表現 GRST 自由—勢分離而硬要求這個案例也具有相同自由纖維。第 4 節固定自由維度而改規則;本節則是邊界與自由結構可以同時改變。兩者處理的是不同的對照問題。


6. 第三個核心命題:固定幾何與初始能量,只改耦合,作用傳遞由零變非零

6.1 C-11-03:兩振子共同物理世界

令 u1,u2u_1,u_2 是兩個質量已歸一的振子位移; ω>0\omega>0 是其相同的局部自然頻率; kc≥0k_c\ge0 是連接兩者的有效耦合剛度:

u¨1+ω2u1+kc(u1−u2)=0,\ddot u_1+\omega^2u_1+k_c(u_1-u_2)=0, u¨2+ω2u2+kc(u2−u1)=0.\ddot u_2+\omega^2u_2+k_c(u_2-u_1)=0.

初始條件完全相同:

u1(0)=u2(0)=0,u˙1(0)=v,u˙2(0)=0.u_1(0)=u_2(0)=0, \qquad \dot u_1(0)=v, \qquad \dot u_2(0)=0.

這表示對第一個振子施加相同的速度脈衝,之後不再提供外部輸入。兩個世界具有同一排列、相同初始位移及速度,以及相同的初始機械能 E(0)=v2/2E(0)=v^2/2。在 kc=0k_c=0 與 kc>0k_c>0 兩種情形,系統的合法狀態空間維度皆相同。差別完全由耦合律引起。

6.2 M-11-05:正常模態與第二振子的遠端響應

定義同向與相對坐標:

s=u1+u22,d=u1−u22.s=\frac{u_1+u_2}{2}, \qquad d=\frac{u_1-u_2}{2}.

則兩種模態彼此解耦:

s¨+ω2s=0,d¨+(ω2+2kc)d=0.\ddot s+\omega^2s=0, \qquad \ddot d+(\omega^2+2k_c)d=0.

定義 Ω=ω2+2kc\Omega=\sqrt{\omega^2+2k_c}。由初態 s˙(0)=d˙(0)=v/2\dot s(0)=\dot d(0)=v/2,得到:

u2(t)=v2[sin⁡(ωt)ω−sin⁡(Ωt)Ω].u_2(t)=\frac v2\left[\frac{\sin(\omega t)}{\omega}-\frac{\sin(\Omega t)}{\Omega}\right].
  • 當 kc=0k_c=0, Ω=ω\Omega=\omega,所以 u2(t)=0u_2(t)=0 對所有時間成立。
  • 當 kc>0k_c>0,上述函數一般不恆為零;例如小時間展開:
u2(t)=kcv6t3+O(t5).u_2(t)=\frac{k_c v}{6}t^3+O(t^5).

因此同一初始幾何、初始自由、初始能量與輸入脈衝,可以因連接規則不同而產生不同的遠端因果作用。證畢。

6.3 這不是宣稱「遠端瞬時作用」

兩個集中振子由無質量彈簧連接,是已粗粒化的理想模型;解析式中遠端對 t>0t>0 的非零響應,不應被解讀為真實長彈簧或連續材料的影響可以無限快傳播。若研究距離及延遲,需要改用具波速的彈性桿、薄板或介質 PDE。物理作用的傳播邊界須保留在與解析度一致的描述裡。

6.4 初始總能量一樣,不表示各規則下的能量函數一樣

此模型的保守能量為:

Ekc(t)=12(u˙12+u˙22)+ω22(u12+u22)+kc2(u1−u2)2.E_{k_c}(t)=\frac12(\dot u_1^2+\dot u_2^2) +\frac{\omega^2}{2}(u_1^2+u_2^2) +\frac{k_c}{2}(u_1-u_2)^2.

在初始位移全為零時,兩種耦合強度下皆有 Ekc(0)=v2/2E_{k_c}(0)=v^2/2。然而 kck_c 不同時,能量函數對一般非零位移狀態的取值可能不同,故不能把不同規則的能量值差直接視為能源憑空流入或流出。兩個世界分別滿足自己的保守能量恆等式;真實切換 kck_c 的瞬間仍須計入裝配或外力所作的功。

6.5 對節拍器矩陣的延伸

同一個平面格點排列下,改變底座是否允許轉動、彈性形變、板材厚度、支承接觸及介質通道,可能改變相應耦合算子及模態。可被檢驗的不是「1000 個格點都要同步」,而是固定同一初態/輸入的反事實響應核如何變化。在只允許剛體平移、質量與方向皆等同的模型中,位置改動可能根本無影響;那應列為必要負例,而不是硬把幾何效應塞回模型。


7. 第四個核心命題:規則與狀態不是可任意線性相加的兩項

7.1 D-11-09:二因子反事實設計

先選一個任務輸出 Y(R,x)Y(R,x),不能先把未知的物理效果定義成「規則差 + 狀態差」。真正的兩因子交互作用可用差中之差:

IXR=Y(R1,x1)−Y(R1,x0)−Y(R0,x1)+Y(R0,x0).I_{XR}=Y(R_1,x_1)-Y(R_1,x_0)-Y(R_0,x_1)+Y(R_0,x_0).

此數值的物理單位與輸出 YY 相同。它依賴參數刻度、干預範圍和研究任務;不同參數化下,「交互主導」不必保持同一數值,因此不能稱為宇宙普遍的交互常數。

7.2 M-11-06:即使最簡線性流也有非零狀態—規則交互

取:

x˙=ax,x(0)=x0.\dot x=ax, \qquad x(0)=x_0.

指定時間 T>0T>0 的輸出:

Y(a,x0)=x0eaT.Y(a,x_0)=x_0e^{aT}.

取 a0≠a1a_0\ne a_1 、 x0≠x1x_0\ne x_1,則:

IaX=(ea1T−ea0T)(x1−x0).I_{aX}=\bigl(e^{a_1T}-e^{a_0T}\bigr)(x_1-x_0).

除特殊退化條件外,此式一般非零。這是動力學流程的時間演化造成乘性關係的結果,不需發明新的物理力才能存在。若只報告「改初態效果」與「改參數效果」而不考察交互項,便可能誤判何者支配某次未來變化。

7.3 有邊界與觀測時的完整設計

若多做 R×X×B×OR\times X\times B\times O 四因子設計,所有主效應和交互項可能需要按照預先定義的效果編碼及階層模型分離;其重點不是每個實驗都必須做完所有組合,而是未分離的作用不可被事後任意歸因。對物理軌跡,觀察 OO 通常是輸出測量因子;若被動換讀出,物理演化的主效應應仍為零。對事件指標及資料可辨識度, OO 可以產生顯著觀測效果。

7.4 交互項不是對「同時改很多東西」的赦免

二因子差中之差雖能識別某種模型內交互,但若更換耦合剛度同時無法避免改動平台重量、摩擦或固定條件,這些都是未控制的伴隨因子,不能因已報 IXRI_{XR} 就當不存在。實物設計需要校準、共變量、隨機化與未量測混淆分析。


8. 被動觀測的相同歷史,與干預可辨識性

8.1 C-11-04:不可辨識的「零輸入、零初態」

令不同參數的有效系統為:

x˙=ax+bu,y=x.\dot x=ax+bu, \qquad y=x.

若長期都使用 u(t)=0u(t)=0 、 x(0)=0x(0)=0,則無論 a,ba,b 是什麼:

y(t)=0.y(t)=0.

因此單靠這些被動資料無法識別 aa 或 bb。這不是 GRST 的獨特新發現,而是系統辨識中已知的激勵不足問題。即使資料取樣無限頻繁,對這個特定實驗世界也不會平白產生缺失的因果識別能力。

8.2 M-11-07:合法輸入讓部分規則產生不同響應

若在同樣零初態下施加已知常數輸入 u(t)=1u(t)=1,且 a≠0a\ne0,則:

y(t)=ba(eat−1).y(t)=\frac ba(e^{at}-1).

a=0a=0 時連續極限為 y(t)=bty(t)=bt。因此不同 (a,b)(a,b) 可以在此干預下被區分。當然,若只量有限個噪聲樣本,參數仍可能實務上不可辨識;單一固定輸入未必足夠識別所有多參數模型。應使用可辨識分析、先驗約束、不同頻率或時間形狀的輸入和獨立驗證。

8.3 D-11-10:任務相對的規則等價類

給定合法初態集 ZΓ\mathcal Z_\Gamma 、輸入/邊界干預族 UΓ\mathcal U_\Gamma 、觀察族 OΓ\mathcal O_\Gamma 與時間窗,定義:

R0∼ΓR1R_0\sim_\Gamma R_1

當且僅當對每個合法組合,兩者經共同型別橋接後的任務輸出相同(或在預先指定的容差內任務等價;容差版本一般只是近似關係,未必具嚴格遞移性)。這將 GRST-09 的 Rashomon 規則族與 GMDIS-10 的任務充分性銜接起來。

結論: 在所有可用的合法干預與觀測下仍然不可區分的兩種規則,研究者不能以「我有一個新的 Shi 名詞」直接宣稱其物理機制已被識別。應報告等價類及未識別的自由,並嘗試設計更具鑑別力的實驗。

8.4 對兩個模型的鑑別力最佳化不是免費的

可提出最佳化候選:

u∗∈arg⁡max⁡u∈UΓDΓ(PR0u,PR1u),u^*\in\arg\max_{u\in\mathcal U_\Gamma}D_\Gamma\bigl(P_{R_0}^{u},P_{R_1}^{u}\bigr),

其中 PRiuP_{R_i}^{u} 是指定任務輸出的機率律或 Dirac 測度, DΓD_\Gamma 是合法的分布差異;需加入最大輸入能量、設備安全、取樣成本及不可觀測隱藏狀態等約束。若直接以模型差異最大化而不限制實驗成本和危險,得到的「最佳實驗」可能是破壞設備或根本不可實現的干預。

8.5 觀察等價不等於物理本體等價

兩個規則可以對某個感測器等價、對另一個感測器不等價;也可以對所有現有可用的干預觀察等價,卻在更大、目前不可實施的干預集合中不同。這是可辨識的邊界,不是說世界沒有客觀規則,更不是說凡無法區分者就必須被任意統一。


9. 觀察干預:唯一不應被直接誤認為物理勢換位的操作

9.1 C-11-05:只換觀察器,不改動力世界

固定同一個振子軌跡 x(t)=(q(t),v(t))x(t)=(q(t),v(t)),可以選擇:

O1(x)=q,O2(x)=(q,v).O_1(x)=q, \qquad O_2(x)=(q,v).

兩個感測器給出不同的資訊充分性與重建表現,但若都為理想被動觀測,底座、阻尼、能量交換、吸引域與穩定性本身沒有改變。物理世界中的 ΣR\Sigma_R 不應僅因更換攝影機座標或讀出維度便被宣布發生物理態勢換位。

9.2 D-11-11:物理換位證書與認識換位證書分離

定義兩種不同記錄:

  • Physical-Shift: 固定共同物理可比量及合法實驗,物理未來的任務指標差超出門檻。
  • Observation-Shift: 在同一物理世界下,由於觀測器/表示/成本改變,模型可辨識性或任務充分性證書有所改變。

兩者有時同時發生(因感測器實際施力或溫度影響介質),但因果理由要分開記錄。記錄 doOdo_O 的物理回饋時應新增對應作用埠,不得把觀測施力視為純抽象資訊增量。

9.3 不應把認知優勢當成新的物理能量

觀測通道更細、原本被忽略的振子模式被發現、或 MIOSA 從 X 轉為充分,可以改變研究者的判斷,但本身不是物理場出現新能量。這種分離亦維持 SIO 原稿的 Representation Non-Ownership 原則:表示可以揭露或隱藏關係,但不因此擁有所揭露的物理真值。


10. 規則勢的可檢驗跨尺度因果響應:未來不只是一次輸出

10.1 D-11-12:有限時間介入響應核

對一個具有適當可微流映射的確定動力系統,令 uku_k 為作用在子域 kk 的合法、帶物理單位的外部擾動, yj(t)y_j(t) 為子域 jj 的物理量測。可定義線性介入響應:

Gjk(t,s)=δyj(t)δuk(s)∣u=u0,t≥s.G_{jk}(t,s)=\left.\frac{\delta y_j(t)}{\delta u_k(s)}\right|_{u=u_0},\qquad t\ge s.

若數學條件不支援 Fréchet 導數或存在碰撞/跳躍,則改用有限幅度實驗:

Δjk(t;η)=yj(t;u0+ηek)−yj(t;u0),\Delta_{jk}(t;\eta)=y_j(t;u_0+\eta e_k)-y_j(t;u_0),

並將干預幅度 η\eta 、其單位、作用時長、時刻與介質條件記錄完整。小擾動響應不是所有大幅度突變、強噪聲或多穩態跨盆地作用的充分表示,必須另外報告非線性及事件依賴。

10.2 D-11-13:物理動態影響陰影的任務化定義

採用 GMDIS-03 的「動態影響陰影」概念,但不宣稱其是空間中固定的額外物理場。給定閾值 δ>0\delta>0 、觀測時間窗及合法擾動集,定義:

Sdyn(k;T,δ)={j:sup⁡0≤t≤T∥Δjk(t;η)∥sj>δ},\mathcal S_{\rm dyn}(k;T,\delta) =\left\{j:\sup_{0\le t\le T}\frac{\|\Delta_{jk}(t;\eta)\|}{s_j}>\delta\right\},

sj>0s_j>0 為預先固定的輸出尺度,讓不同單位可公平比較。變更 T,δT,\delta 、觀測器或干預幅度可能改變陰影範圍;這不是「陰影不真實」,而是證書依指定實驗任務定義。如果延遲或波速重要,應另外記錄首達時間,不可將跨距離的小延遲混入一個無時間標記的非局部網絡。

10.3 M-11-08:有限時間影響等價不推出長期穩定性等價

在任意有限 TT 內,兩個模型的觀測輸出可以任意接近,而它們的漸近穩定結構仍可能不同。例如第 4 節叉式分岔模型取 μ0=−ε\mu_0=-\varepsilon 、 μ1=+ε\mu_1=+\varepsilon,初始擾動極小, TT 固定,當 ε→0\varepsilon\to0 時兩條有限時軌道差趨近零;但只要 ε≠0\varepsilon\ne0,零平衡點的線性穩定號數仍不同。因而:

有限時間近似等價⇏漸近穩定類相同.\text{有限時間近似等價}\not\Rightarrow\text{漸近穩定類相同}.

這不是對所有數值方法的否定,而是說實驗時間窗與必要的長期結論之間,應有明確的誤差與辨識邊界。

10.4 D-11-14:反事實世界線族

固定共同狀態 zz 與合法干預族 I\mathfrak I,研究世界族:

GΓ(z;W)={YWι(⋅;z):ι∈I}.\mathscr G_{\Gamma}(z;W)=\bigl\{Y_{W}^{\iota}(\cdot;z):\iota\in\mathfrak I\bigr\}.

這不是宣稱所有世界線均已能觀測,也不自動賦予決定性世界線「客觀機率分布」。它是在預註冊試驗集合下可展開的未來軌道束;其中每條具有自身初始條件、作用路徑、能量限制與感測約束。對隨機模型,應改用帶條件的路徑分布族。

10.5 GRST 勢剖面與作用核的關係

GjkG_{jk} 、有限時間 Lyapunov 指數、穩定盆地及能量流能夠構成 GRST 候選未來剖面的物理特徵,但這不證明它們與抽象 ΣR\Sigma_R 存在唯一同構。不同規則可能給出同一組有限觀察特徵;同一規則在不同時間窗可被標成不同操作性態勢。研究的原則應是:保留特徵、證書、出處與失敗邊界,而非宣布已測量到「Shi 粒子」或「Shi 的物理量綱」。


11. 規則與邊界何時應被視為真正不同的干預型別?

11.1 物理數學中的邊界通常本來就是算子的一部分

一個 PDE 的生成元不只取決於形式微分算子,也取決於定義域及邊界條件。例如熱方程雖都使用 κ∂xx\kappa\partial_{xx},在 Dirichlet 與 Neumann 情形其算子定義域、譜、長期狀態不同。因此嚴格地說,若將「完整規則」定為生成元及其定義域,則改邊界也是改 RR。GMDIS 的 R/BR/B 分離是可審計的研究設計因子,而不是所有物理描述下不可再分的本體二元。

11.2 D-11-15:統一標記與型別保留

允許完整實施機制以結構對表示:

M=(R,B),\mathsf M=(R,B),

其中 RR 為材料、內部交互作用與局部轉移方程, BB 為容器幾何與外部交換條件。某些模型化程度下 BB 可以嵌入 RR,但反事實記錄仍保留「研究者有意改變了哪個實驗因子」的欄位。這讓兩種系統描述可以互譯,而不致因符號習慣不同產生假性發現。

11.3 規則干預也可被擴展狀態吸收

若一個參數 θ\theta 在單次實驗中保持常數,原系統為:

x˙=f(x;θ).\dot x=f(x;\theta).

可建立增廣狀態 z=(x,θ)z=(x,\theta),使:

z˙=(f(x;θ),0).\dot z=\bigl(f(x;\theta),0\bigr).

在此增廣模型中,「改參數」能被表述為「準備不同初始 θ\theta 」。因此 doRdo_R 與 doXdo_X 並非任何數學表達下都絕對不能互換的本體範疇。有用的區分是:研究者究竟只替換某次狀態,還是更換在整類狀態、輸入及未來中都持續生效的生成機制;該更換能否由可實施的物理控制實現。若增廣模型已經提供所有所需比較而更簡潔,GMDIS 沒有理由拒絕它。

11.4 這是本篇的重要讓位條件

GRST-07 的雙層反事實優勢若完全可由一個標準的增廣狀態/結構因果模型無損表達,新增符號應縮減為研究註記規約;反之若顯式型別分離能系統性減少錯誤歸因或改善設計效率,才具有可測工程增益。新系列不把哲學語詞本身當作免於比較的證據。


12. 系統辨識、Rashomon 規則族與未來反事實的可識別性

12.1 D-11-16:結構性與實務性辨識分開

  • 結構性可辨識: 在理想無噪、充分長或理論允許的觀測下,參數/規則是否唯一。
  • 實務性可辨識: 在現有取樣率、噪聲、干預幅度與資源限制下,參數是否能被精確估計。
  • 任務性可辨識: 即使參數不唯一,不同候選規則是否對目標物理決策產生不同結果。

三者不可互相取代。Raue 等對部分觀測動力模型的 profile likelihood 工作(文獻 E-11-04)顯示,模型看似貼合資料,參數仍可能無法被穩健識別。

12.2 D-11-17:反事實區分試驗

若有候選 R1,R2R_1,R_2,應比較它們在同樣的合法輸入下所產生的可見輸出:

δI(R1,R2)=sup⁡ι∈I dY(YR1ι,YR2ι).\delta_{\mathcal I}(R_1,R_2)=\sup_{\iota\in\mathcal I}\,d_{\mathcal Y}\bigl(Y_{R_1}^{\iota},Y_{R_2}^{\iota}\bigr).

若 δI=0\delta_{\mathcal I}=0,則在所選干預/觀測族下無法區分。若大於零但小於儀器不確定度或模型結構誤差,也不應宣稱實驗已識別。研究應報告可區分度的誤差下界與置信範圍,而不是只給單一 p 值。

12.3 MIOSA 必須擴展到跨規則干預

第十篇 MIOSA 主要審計一組感測是否足以支援任務;第十一篇增加:同一感測方案是否能區分模型世界族在合法干預下的未來。若沒有干預,觀測方案可能對被動任務充分,對機制判定不充分。若在新方案中單一便宜感測器加適當輸入便足夠,也應撤回「需要高昂多感測器」的事先偏見。

12.4 從殘差到機制,禁止跨越的推理

若模型 A 和 B 的預測殘差不同,只能說某個模型在指定資料與評分上更好。不能直接把較小殘差當成「模型 A 所有物理作用都是真的」。兩套方程可能在短資料窗內提供等價預測;對模型選擇還需結構辨識、守恆審查、參數合理性、適用域與獨立實驗。SINDy 類稀疏辨識工具是強基線,GMDIS 必須與它正面比較,而不能只以自創名詞勝出。

12.5 不同感測足以驗證動力,未必足以驗證因果

僅憑時間相干、同步率或高互資訊,不足以排除共同外部驅動或隱藏介質。研究必須能夠施加或利用可信的外源性擾動、建立因果時序、比較實際響應方向與能量路徑,才可以把某種傳播判為物理因果,而不是相關投影。


13. 24/72、SIO、GRST 在本篇的分工,不得混線

13.1 原始 24/72 分類坐標不變

24 範式的底空間軸 B2B_2 (連續/離散)、更新軸 U4U_4 (序列、選擇/跳躍、並行、識別/存取)、觀察軸 O3O_3 (連續、離散、條件性拒單測 X)全部保持原始語義。72 格新增最低充分轉移律軸 L3L_3:F(確定函數)、K(古典機率核)、Q(必要的量子通道)。24/72 用於標記模型事件及求解/觀察方式,不能用其代碼替代物理定律。

13.2 反事實事件的多層標籤

同一研究可能包含:

層 事件 分類/物理意義
W/世界 材料、擺錘、底座、聲腔 物理世界不是一格 72 的「硬體物種」
M/模型 兩振子連續確定微分方程 F 類轉移;連續底空間表示
N/數值 時間離散、網格、向量化 所用求解算法可標序列/並行等,但不改世界的時間本體
I/干預 改 kck_c 、固定底座、局部脈衝 明確的 doRdo_R 、 doBdo_B 、 doXdo_X
O/觀察 位移、相位、音壓、恢復率 依任務決定 C/D/X 的充分性證書
A/稽核 能量核對、反事實差、MIOSA 保留模型條件、數值與不確定度證據

量子 Q 僅在物理量子態與相干通道在該任務不可約時使用;不能因為「節拍器有相位」或「物理波有相干」就將經典振子分類為量子 Q。

13.3 SIO:先看自由,再看可行未來如何不同

在第 4 節同一叉式模型, μ\mu 變化不改一維狀態空間,但改變原點的穩定性及非零平衡點。這展示「自由度相同」不足以決定未來動力構形。第 5 節在不同邊界下,合法的演化函數空間及長期熱含量都可能改變;故先判定是否屬於硬約束變化,再討論規則勢剖面。第 6 節耦合強度改變了動態影響陰影,但並無理由宣稱局部物理粒子數或物理位移自由度因此被消除。

13.4 GRST:把規則、自由、賦形與未來保持為不同語義層

在某個測試世界,若 SΓ(W0,x)≠SΓ(W1,x)S_{\Gamma}(W_0,x)\ne S_{\Gamma}(W_1,x),只能給出指定研究條件下的態勢差證書。不能由此倒推唯一真實規則;有許多參數與邊界組合可能得出相同剖面。GRST-09 的規則族/不可辨識觀點在此仍有效。

13.5 與早期萬有張力的距離

萬有張力提供「交互作用持續交換與競爭」的啟發。但物理比較必須分離:接觸力的牛頓值、能量/功率的 SI 單位、應力與壓力、阻尼、波傳與延遲,以及本篇任務相對的未來勢剖面。跨介質「張力總和」若沒有單位、狀態方程與守恆設計,就不能視為可驗證數學物理量。


14. 針對多介質真實系統的因果設計

14.1 最小目標並非一開始涵蓋整個宇宙

物理實驗應從兩振子及已知邊界的熱模型開始,逐步擴展到少數機械節拍器、共享可移動底座、可轉動板、彈性平台、聲學介質、空氣耗散與擒縱非線性。每增加一層,必須保留既有基線的守恆、數值收斂與資料對齊。一次把固體、流體、熱、量子、觀測者全部塞進一個尚未校準的「全域張力總算子」,無法成為可靠的實驗研究。

14.2 D-11-18:節拍器矩陣的四類可執行干預

干預 操作示例 主要保持量 預註冊讀出 風險
狀態 doXdo_X 對單擺施加校準的微小角度偏移 R,B,OR,B,O 、能源裝置設定 恢復時間、相位差、模態變化、音壓響應 擾動帶入不同動能
規則 doRdo_R 更換等效彈性/阻尼參數或擒縱設定 幾何與初態、外部驅動政策 轉移核、局部穩定譜、吸引域、能耗 改裝可能同時改質量及接觸
邊界 doBdo_B 固定/釋放底座,改支點或反射聲腔 同一節拍器與驅動參數 底座速度、反射、音場、功率 邊界重設連帶改初態
觀察 doOdo_O 替換追蹤相機、音壓計或採樣率 理想情形下全部物理參數 MIOSA、可辨識度、噪聲 感測負載與時間同步誤差

在真實實驗中,同一初態應至少同時記錄所有擺角、角速度、底座位姿/速度及具有顯著影響的隱藏能量狀態;若不能控制空氣聲場的初態,應隨機化、暖機及追蹤背景噪聲或將其列為不確定度,而非默認零。

14.3 幾何的零效應是必需的負例

對只允許共用底座一維剛體平移、各節拍器質量與運動方向完全相同的理想模型,單純把節拍器在平面上重新排列而不改任何力臂與約束,通常不會改變耦合方程。該幾何零效應不得被修改成「理論也能解釋,因為 GRST 剪枝」。它是模型正確性檢查。只有允許旋轉力矩、板形變、非均勻質量、實際聲場與邊界時,排列位置才具有可研究的幾何因果途徑。

14.4 多模態的真實同步不是單一吸引子標籤

至少分開記錄:完全同向、反相、群聚/簇、行波、漂移、無穩定鎖定和振盪死亡,並區分「彼此同步」和「對一個外部時鐘鎖定」。相同同步序參量可能隱藏不同能流、外部驅動及橫向擾動恢復。GMDIS-07 的持續驅動模型與 GMDIS-08 的記憶核在此為強基線。

14.5 只有條件確定、資料夠多才可給物理態勢差證書

物理實驗報告不能只展示某次影片的兩種結局。至少需要:多個獨立試次、跨初態與干預幅度重複、隨機或平衡操作次序、記錄裝置老化、機械摩擦改變和能源補給狀態。若裝置間差異會隨時間變化,必須考慮漂移模型和配對設計,而不是直接套單一同分布假設。


15. 實驗計畫:反事實矩陣、預註冊指標與最小比較基線

15.1 D-11-19:四因素矩陣的最小形式

取兩個規則設定 R0,R1R_0,R_1 、兩個初態設定 X0,X1X_0,X_1 、兩種邊界 B0,B1B_0,B_1 、兩種被動觀測 O0,O1O_0,O_1,形成候選 24=162^4=16 個格子。但不宣稱所有 16 種組合都物理合法:若更換規則導致某個初態不可實現,該組合應標「inapplicable」並建立橋接或重設研究範圍。

真正實驗可採分段設計:先做 R×XR\times X 的 4 格,確認基本差分;再加入 BB,定位邊界與規則交互;最後以相同物理試驗同步收集多種被動觀測。這可以避免把換感測當作「再做一組新的物理世界」。

15.2 反事實矩陣的示例

對第 7 節簡單流 x˙=ax\dot x=ax 、指定 T=2T=2,取:

a0=−0.6,a1=0.4,x0=0,x1=0.1.a_0=-0.6,\quad a_1=0.4,\quad x_0=0,\quad x_1=0.1.

四格數值由精確式 Y=xeaTY=x e^{aT} 生成:

狀態 x0=0x_0=0 狀態 x1=0.1x_1=0.1
規則 a0=−0.6a_0=-0.6 00 0.1e−1.20.1e^{-1.2}
規則 a1=0.4a_1=0.4 00 0.1e0.80.1e^{0.8}

因而 IaX=0.1(e0.8−e−1.2)I_{aX}=0.1(e^{0.8}-e^{-1.2}),明確非零。這不是四次實物測量,而是供程式核對的模型輸出。實驗中須給每格重複測試與不確定度。

15.3 共同最低輸出剖面

每組實驗建議至少記錄以下欄位:

型別 量 典型單位 檢驗理由
物理狀態 擺角、角速度、平台位移/速度 rad、rad/s、m、m/s 避免把同相位當同全態
力學/聲學 加速度、接觸力、聲壓、界面功率 m/s²、N、Pa、W 跨介質能量與作用通路
穩定性 橫向恢復率、模態轉換、盆地入界 1/s、事件、無量綱 直接比較未來結構
介面/能量 外部功率、散逸熱、儲能、殘差 W、J、J、J 阻止重複計能及虛假能量
觀察 採樣率、感測漂移、噪聲、畫素追蹤誤差 Hz、特定 SI 單位 MIOSA 及儀器限制
模型 規則參數、邊界版本、求解器、步長 對應 SI、版本 ID 讓反事實有可重現實體
推斷 效應量、區間估計、等價類、失效標籤 任務指定 區分資料支持與模型詮釋

15.4 必須設的強基線

對 GMDIS-12/13,應至少同時比較:

  1. B0: 傳統有限維經典振盪器模型,直接以機械參數進行系統辨識。
  2. B1: 帶振幅的相位模型或可用的相位—等穩幅近似,而非只用最弱的 Kuramoto 形式。
  3. B2: 真實可行的底座移動、轉動及形變模型,加上必要的聲場/熱場修正。
  4. B3: 完整輸入輸出辨識方法,例如狀態空間辨識、SINDy 或適當的混合系統辨識。
  5. B4: GMDIS 的 24/72、SIO、GRST 反事實剖面與型別稽核。

核心公平規則: 若 B4 只是 B2/B3 的所有物理變數重新命名,而沒有額外預測、比較、建模錯誤辨識或操作成本增益,則 B4 在該任務上的獨立工程價值不成立。

15.5 原則性的訓練/測試切分

不要把同一連續模擬軌跡相鄰的 1000 幀當作 1000 組獨立試次。建議以獨立初態、實驗日、材料、外部溫度與平台構造進行 hold-out。最有鑑別力的測試是在訓練階段沒看過的參數/邊界干預下預測輸出,例如訓練在 kck_c 低值,測試在中等與高值;或訓練在固定底座,測試在釋放/轉動底座。必須明示這類外推是否超過模型有效範圍。

15.6 預註冊範例

至少在資料取得前寫下:主要終點、替代終點、接受效應的最低尺度、儀器精度、樣本量規劃、隨機化單位、共變量調整、實驗撤回理由、外推禁區、固定軟體版本與不利結果的保存方式。若實驗中期改了定義,必須另存為探索性分析,不應回寫成預先成功的理論預測。


16. 候選新貢獻與明確判敗:不能讓 GRST 成為萬能解釋詞

16.1 H-11-01:規則勢換位可由可重複物理干預辨認

假說: 在部分節拍器多介質實驗,存在保持物理幾何和初態近似相同的設定變更,使指定未來剖面中至少一項出現可重複顯著差異;且某些差異在單純相位指標下並不易辨認。判敗: 若成熟的機械模型已完整解釋,而且新指標沒有超出其合理衍生特徵的效果,便不得宣稱增量物理新律。

16.2 H-11-02:跨規則干預有助於釐清假性同步因果

假說: 對共同外部驅動、底座機械耦合、聲音耦合的不同情形,四層干預協議比單純比較相位相關能更穩健地識別主作用通路。判敗: 若標準控制理論與現有因果實驗設計已在同等成本下給出相同或更好表現,新框架應退位為跨領域記錄慣例。

16.3 H-11-03:SIO 陰影指標對場效應具有跨設定可遷移性

假說: 在校準過的時間窗與介質條件下,物理動態影響陰影可改善對干預後因果可達範圍或能量傳播延遲的預測。判敗: 若影響陰影高度依賴閾值、測量噪聲和任意正規化,而不能穩定跨樣本外重現,則它只保留為視覺化/稽核輔助。

16.4 H-11-04:GRST 多不變量剖面改善規則族辨識

假說: 在同一成本和同一物理模型類下,包含能量通路、模態回復與干預敏感度的剖面比只有相位或單一序參量更能分開對任務重要的候選規則族。判敗: 如更少、更標準的可觀測性 Gramian、轉移函數或系統參數已達同等能力,應降低 GRST 名義的新穎性宣稱。

16.5 H-11-05:高精度物理模型可能證明「新框架」不必要

這必須作為真正允許的結果:完整機械與多場模型可以在所有指定實驗中成功預測同步、反相、分群及耗散,而將其重新標為 GRST/SIO/72 格並無改善。屆時研究價值可能是清晰化的跨領域模型分類,而非獨立科學定律。可被更簡潔既有方法吸收是此系列保持科學誠實的正式路徑,不是該掩蓋的失敗。

16.6 M-11 結果與 H-11 假說不得互相挪用

第 4–8 節解析證明的是指定方程中的結構差異,不能作為 H-11-01 至 H-11-04 對實際節拍器成立的直接實物證據。反過來,真實節拍器若出現某些非線性模式,也不會因此自動證明全部 GRST 世界本體論命題。


17. 本篇已完成/尚未完成的研究邊界

17.1 已完成的可審核工作

  • 清晰地分離四種干預的研究型別及其非唯一的本體界線。
  • 在一個可解的一維非線性模型中給出「同軌道、異穩定性」的解析證明。
  • 在相同初態熱擴散 PDE 中給出 Neumann 與 Dirichlet 未來差異的解析證明。
  • 在兩振子模型中推導相同初始能量、相同幾何、不同耦合產生不同遠端影響的確定式。
  • 給出規則—狀態交互項的精確差中之差公式。
  • 提出規則族相對等價、觀測換位、最小充分干預與多不變量物理稽核。
  • 提供後附可重現的純 NumPy 模型核對程式及文獻/來源審計。

17.2 尚未完成、不應冒稱已完成的工作

  • 100、1000 台實體節拍器的物理反事實測量。
  • 全固體—流體—熱—擒縱 PDE/ODE/事件閉合的全域存在唯一性定理。
  • 跨全部物理尺度通用的「唯一 Shi 距離」或「新資訊守恆律」。
  • GRST 對經典系統辨識方法的實測優越性。
  • 將所有有效模型參數干預唯一歸因於世界基本法則的能力。
  • 任何由上述 toy model 推出的宇宙普遍同步、臨界或必然收斂定律。

17.3 物理測量上的判定應誠實保留三種中間狀態

「支持特定效應」「拒絕特定效應」「資料不足而未判定」必須相互區分。測不到效果不等於證明完全無效,有限範圍觀察有差異也不等於證明所有條件下均有差異。對不可比較情況,還需要第四種「未合法對齊」狀態,不應放進顯著性計算。


18. 對後續三篇的正式交付

18.1 GMDIS-12:節拍器矩陣的多尺度基準世界

GMDIS-12 必須接收本篇的「四層干預資料契約」,優先用兩振子機械耦合、固定/可動/可轉底座、聲場、擒縱驅動、熱耗散等已校準子域組成逐層擴充的基準世界。每個實驗應提供兩份資料:反事實設計與實際可重現輸出。如模型細節未定而無法計算某反事實,應標未支持,不得虛構實驗結果。

18.2 GMDIS-13:基線競爭、消融與可反駁性

GMDIS-13 必須檢驗:不同研究組別的初態可比性、真實能量帳本、同模型參數外推、觀測不足對因果的影響、強基線的公平資源分配,以及 GMDIS 特徵剖面的實際增益。應允許「所有新符號都被成熟系統辨識方法吸收」的負面結果。

18.3 GMDIS-14:全域多尺度動力學統合

GMDIS-14 不能因前 13 篇名詞已形成完整母鏈,就提前宣稱具有新的統一物理定律。只有被解析與數值/真實實驗雙重支持的機制可以進入已證結果;其餘命題必須保留「數學候選」「研究性詮釋」「未驗證實物」或「失敗邊界」的明確地位。系列收束應強調可生成、可干預、可稽核、可反駁的多尺度世界模型。


19. 結論:未來不同,究竟是誰改變了世界?

若只是把兩張節拍器同步影片放在一起,看到一張同向擺、一張反相擺,仍無法判定:兩者差異究竟源於起始相位、某擺的動能、底座耦合、空氣聲場、外部電源、邊界固定方式,還是攝影機採樣偏差。GRST 的「態勢」提示我們不要只看當前外觀,SIO 的自由度觀提示我們不要把幾何佔據當作全部作用;但要成為可研究科學,就必須再往前走一步——將候選生成機制及邊界以合法、可重複的方式實際干預,並從共同物理世界的未來對照中排除混淆。

本篇最重要的結論不是「物理世界必然走向強勢穩定同步」,而是:

相同位置、相同初態甚至相同被動歷史,都不足以無條件保證相同未來。真正可比較的未來位置,必須相對於生成機制、邊界、合法干預、觀察與任務而定義;若不同世界的未來差異能在控制條件下穩定被辨識,才形成一份可審計的物理態勢換位證書。

與此相反,若更換的只是理想被動觀測器、或兩套候選規則在全部合法實驗下等價,則不得硬宣稱物理世界已換位。這個正例與負例的共同存在,才使 GRST 與 SIO 的新應用有接受反駁、降維與讓位的科學條件。


附錄 A|解析證明、可比較性與物理單位的覆核

A.1 叉式分岔之精確解

對 x˙=μx−x3\dot x=\mu x-x^3,若 x≠0x\ne0,置 v=x2v=x^2,得到 v˙=2μv−2v2\dot v=2\mu v-2v^2。對 μ≠0\mu\ne0 的 logistic 類方程,其顯式解為:

v(t)=μv0e2μtμ+v0(e2μt−1).v(t)=\frac{\mu v_0e^{2\mu t}}{\mu+v_0(e^{2\mu t}-1)}.

因 v0=x02v_0=x_0^2,若 x0>0x_0>0,取正平方根即正文公式;若 x0<0x_0<0,取負平方根。對 μ<0\mu<0 且 t≥0t\ge0,分母分子同號, v(t)v(t) 保持非負;不應把分母的未整理符號誤讀為負數平方根。若 μ=0\mu=0,則 v˙=−2v2\dot v=-2v^2,故 v(t)=v0/(1+2v0t)v(t)=v_0/(1+2v_0t)。 x=0x=0 為每個 μ\mu 的恆定解,但它的擾動敏感度隨 μ\mu 變號,這就是正文「被動同歷史,局部未來換位」的最小證明。

A.2 熱方程 Fourier 係數的精確核對

令 L=1L=1 作無量綱基準。初態 u0(x)=sin⁡2(πx)=(1−cos⁡2πx)/2u_0(x)=\sin^2(\pi x)=(1-\cos 2\pi x)/2。其 Dirichlet 正弦係數:

bn=2∫01sin⁡2(πx)sin⁡(nπx) dx.b_n=2\int_0^1\sin^2(\pi x)\sin(n\pi x)\,dx.

因 u0(1−x)=u0(x)u_0(1-x)=u_0(x) 而 sin⁡(nπ(1−x))=(−1)n+1sin⁡(nπx)\sin(n\pi(1-x))=(-1)^{n+1}\sin(n\pi x),偶數 nn 的積分為零;奇數時由積化和差:

∫01cos⁡(2πx)sin⁡(nπx) dx=2nπ(n2−4).\int_0^1\cos(2\pi x)\sin(n\pi x)\,dx=\frac{2n}{\pi(n^2-4)}.

因此:

bn=2nπ−2nπ(n2−4)=−8πn(n2−4).b_n=\frac2{n\pi}-\frac{2n}{\pi(n^2-4)}=-\frac8{\pi n(n^2-4)}.

n=1n=1 給 b1=8/(3π)>0b_1=8/(3\pi)>0 ; n=3n=3 給 b3=−8/(15π)<0b_3=-8/(15\pi)<0。有限截斷在 t=0t=0 的總量重建具有非零尾項誤差,因此附錄程式對初時積分使用比浮點捨入更寬的、明示的 Fourier 截斷容忍度;對 t>0t>0 則因高頻指數衰減而迅速變精確。

邊界與熱帳本: 若 ρ,c\rho,c 不變且截面積為 AA,則物理相對熱含量 Q(t)=ρcA∫0LΔT u(x,t) dxQ(t)=\rho c A\int_0^L \Delta T\,u(x,t)\,dx,其變化由兩端的導熱通量決定。Dirichlet 邊界是對外部熱庫的理想有效條件;熱量從桿離開並不代表總能量在封閉大世界中消失。

A.3 兩振子衝量的局部時間展開

有 Ω2=ω2+2kc\Omega^2=\omega^2+2k_c,且:

sin⁡(at)a=t−a2t36+a4t5120+O(t7).\frac{\sin(at)}a=t-\frac{a^2t^3}{6}+\frac{a^4t^5}{120}+O(t^7).

故:

sin⁡(ωt)ω−sin⁡(Ωt)Ω=Ω2−ω26t3+O(t5)=kc3t3+O(t5).\frac{\sin(\omega t)}\omega-\frac{\sin(\Omega t)}\Omega =\frac{\Omega^2-\omega^2}{6}t^3+O(t^5) =\frac{k_c}{3}t^3+O(t^5).

乘上 v/2v/2 得 u2(t)=(kcv/6)t3+O(t5)u_2(t)=(k_cv/6)t^3+O(t^5)。在此簡化模型中 kck_c 的單位為 s−2\mathrm s^{-2} (已除以單位質量);若改回有量綱彈簧剛度 KK 與質量 mm,應使用 kc=K/mk_c=K/m,不可直接將兩種不同單位的參數互換。

A.4 功率與動量:不容許的「全球力守恆」捷徑

在牛頓力學中,力是動量變化率;同時存在多個接觸力、外力、慣性參考效應及約束反力時,應逐控制體建立平衡式,不可直接宣稱「總力是一個普遍守恆量」。對兩振子,互連彈簧作用大小相等、方向相反,當作內力在總動量方程中相消;但各自的能量仍可能彼此交換。對固定外部彈簧,還存在外界支承反力,不能由一對彈簧耦合的相消推出整個開放系統動量必定恆定。與 GRST 勢特徵的任何比較都必須維持此分離。

A.5 反事實交互項的單位與依賴

對 Y(a,x)=xeaTY(a,x)=xe^{aT}, aa 單位 1/s1/\mathrm s, TT 單位秒, aTaT 無量綱;差中之差 IaXI_{aX} 與 xx 同單位。若選擇 ln⁡x\ln x 作輸出、或更換初始狀態的參數化,交互項的數值可能不同。故「非零交互」只是一個相對輸出與實驗因子的具體結果,不是不存在座標依賴的宇宙力項。

A.6 何種情況「同一狀態」不合法?

  • 幾何尺寸變更後,原始位置位於新材料之外;此時初態無法直接傳送。
  • 無滑動接觸改為允許滑動,速度切向約束被改變,原初態的相容性必須重查。
  • Dirichlet 邊界改為 Neumann,若初態違反任一邊界的強相容條件,不能直接比較同一類古典解;可降到適當弱解框架,並如實標註。
  • 熱擾動與機械擾動未匹配時,初始壓力場或隱藏聲腔態可能已經不同。
  • 量子世界與古典世界不能只因都用「相位」便假定具同一狀態;應明定可觀察代數與合法介面映射。
  • 一個世界從 100 個單位改為 1000 個單位,不具有自動的一對一全態等號;需縮放、嵌入、統計對齊或共同可觀察量。

如果狀態橋接不可構造,則正確記錄是 NOT_COMPARABLE,而不是理論只要足夠宏大就可以任意跨越。


附錄 B|可獨立執行的八組 NumPy 合成模型驗證

B.1 測試責任

本附錄為 N-11-01 至 N-11-08 的完整可執行程式,需 Python 3.10+ 與 NumPy,不依賴外部網路或自研計算框架。其工作只限於:確認正文中的解析解、係數、反事實差、能量帳本與同觀測不同規則的模型內行為一致。程式沒有測量自然物理,不接受用於宣稱真實 100/1000 節拍器實驗成功。

執行方式:將下方程式碼區塊原樣存為 gmdis11_check.py,執行 python gmdis11_check.py。本篇的正式研究來源仍是本 Markdown;程式可作為日後產生測試報告的附件來源。

B.2 完整測試程式

# GMDIS-11 reproducible synthetic checks; Python 3.10+, NumPy only.
import math
import numpy as np

RTOL, ATOL = 1e-11, 1e-12

def assert_close(x, y, msg, rtol=RTOL, atol=ATOL):
    if not np.allclose(x, y, rtol=rtol, atol=atol):
        raise AssertionError(f'{msg}: got {x}, expected {y}')

def pitchfork(mu, delta, t):
    if abs(mu) < 1e-14:
        return delta / math.sqrt(1 + 2 * delta**2 * t)
    den = 1 + delta**2 * math.expm1(2 * mu * t) / mu
    return delta * math.exp(mu * t) / math.sqrt(den)

def heat_dirichlet_mass(kappa, t, n_odd=300):
    n = np.arange(1, 2*n_odd, 2, dtype=float)
    b = -8/(np.pi*n*(n*n-4))
    return float(np.sum(2*b/(np.pi*n)*np.exp(-kappa*(np.pi*n)**2*t)))

def heat_neumann_mass(kappa,t):
    return 0.5

def heat_neumann(x, kappa, t):
    return .5 - .5*np.exp(-4*np.pi**2*kappa*t)*np.cos(2*np.pi*x)

def impulse_second(kc, omega, v, t):
    Om = math.sqrt(omega**2+2*kc)
    return .5*v*(math.sin(omega*t)/omega - math.sin(Om*t)/Om)

def exponential(a,x,t):
    return x*math.exp(a*t)

def forced_zero(a,b,t):
    return b*t if abs(a)<1e-14 else b*math.expm1(a*t)/a

def tests():
    # N-11-01: same unperturbed equilibrium, distinct stability and perturbed futures
    mu_neg, mu_pos, delta, T = -.4, .4, .01, 5.
    y_neg, y_pos = pitchfork(mu_neg, delta, T), pitchfork(mu_pos, delta, T)
    assert 0 < y_neg < delta < y_pos < math.sqrt(mu_pos)
    assert_close(pitchfork(mu_neg,0,T),pitchfork(mu_pos,0,T),'same zero trajectory')
    assert_close(pitchfork(.4,delta,0),delta,'pitchfork initial')
    # Exact solution derivative equals right-hand side for nonzero parameter
    for mu in [-.4,.4]:
        t=.7; hh=1e-5
        diff=(pitchfork(mu,delta,t+hh)-pitchfork(mu,delta,t-hh))/(2*hh)
        yy=pitchfork(mu,delta,t)
        assert_close(diff, mu*yy-yy**3, 'pitchfork ODE', rtol=1e-8, atol=1e-9)
    print(f'N11-01 PASS pitchfork: negative_mu={y_neg:.10f}, positive_mu={y_pos:.10f}')
    # N-11-02: same heat initial field, different boundary -> different mass
    kappa, t = .1, 1.0
    assert_close(heat_dirichlet_mass(kappa,0,300),.5,'initial Dirichlet mass (truncated series)',atol=2e-9)
    MN=heat_neumann_mass(kappa,t); MD=heat_dirichlet_mass(kappa,t)
    assert 0 < MD < MN
    x=np.linspace(0,1,2001)
    nx=heat_neumann(x,kappa,t)
    numerical_mean=float(np.sum(.5*(nx[1:]+nx[:-1])*np.diff(x)))
    assert_close(numerical_mean,MN,'Neumann mass quadrature',atol=1e-9)
    assert abs(heat_dirichlet_mass(kappa,5))<abs(MD)
    print(f'N11-02 PASS heat: Neumann_mass={MN:.10f}, Dirichlet_mass={MD:.10f}, contrast={MN-MD:.10f}')
    # N-11-03: fixed geometry, impulse and energy, coupling change
    omega, v, k1,t=.9,1.,.6,.05
    assert_close(impulse_second(0,omega,v,t),0,'decoupled impulse',atol=1e-16)
    resp=impulse_second(k1,omega,v,t)
    leading=k1*v*t**3/6
    assert abs(resp-leading)/abs(leading) < .005
    print(f'N11-03 PASS impulse: coupled={resp:.11e}, leading={leading:.11e}, relative_error={abs(resp-leading)/abs(leading):.7f}')
    # N-11-04: finite factorial contrast
    a0,a1,x0,x1,T=-.6,.4,0.,.1,2.
    dd=exponential(a1,x1,T)-exponential(a1,x0,T)-exponential(a0,x1,T)+exponential(a0,x0,T)
    expected=(math.exp(a1*T)-math.exp(a0*T))*(x1-x0)
    assert_close(dd,expected,'double difference')
    print(f'N11-04 PASS factorial: double_difference={dd:.10f}')
    # N-11-05: passive non-identifiable, active distinguishability
    for a,b in [(-1.,.7),(.8,2.),(0,5.)]:
        assert_close(exponential(a,0,2),0,'passive zero trajectory')
    yA,yB=forced_zero(-.5,1,1),forced_zero(.5,1,1)
    assert abs(yA-yB)>.1
    assert_close(forced_zero(0,1,1),1,'a=0 continuous limit')
    print(f'N11-05 PASS identification: forced_minus={yA:.10f}, forced_plus={yB:.10f}')
    # N-11-06: observer change affects data, not physical energy or trajectory
    t=.67; q=math.cos(t); velocity=-math.sin(t)
    obs_1=q;obs_2=np.array([q,velocity]);E=.5*(q*q+velocity*velocity)
    assert_close(obs_1,obs_2[0],'passive observation same state')
    assert_close(E,.5,'harmonic energy')
    print(f'N11-06 PASS passive observer: scalar={obs_1:.10f}, two_channel_energy={E:.10f}')
    # N-11-07: rule-dependent energy is constant within each world
    for k in [1.,4.]:
        for t in np.linspace(0,5,101):
            q=math.cos(math.sqrt(k)*t)
            vel=-math.sqrt(k)*math.sin(math.sqrt(k)*t)
            assert_close(.5*(vel**2 + k*q*q),.5*k,'constant rule-specific E')
    print('N11-07 PASS rule-specific conservation: E(k=1)=0.5; E(k=4)=2.0')
    # N-11-08: intervention sign and time scale sensitivity
    for t in [0.01,0.1,0.5,1.0]:
        a=impulse_second(.6,.9,1.,t)
        assert a>0, 'coupled response should be positive on chosen finite horizon'
    print('N11-08 PASS finite-time directional response for chosen parameter window')
    print('N-11 CHECKS: 8/8 PASS; all synthetic model checks only, not physical experiments.')
if __name__=='__main__':
    tests()

B.3 合成輸出的數值回報

在 Python 3.11+/NumPy(實際執行工具環境)下,合成檢查輸出如下;均僅對已指定模型與程式設定有效:

N11-01 PASS pitchfork: negative_mu=0.0013531868, positive_mu=0.0734004315
N11-02 PASS heat: Neumann_mass=0.5000000000, Dirichlet_mass=0.2013987306, contrast=0.2986012694
N11-03 PASS impulse: coupled=1.24955943382e-05, leading=1.25000000000e-05, relative_error=0.0003525
N11-04 PASS factorial: double_difference=0.1924346717
N11-05 PASS identification: forced_minus=0.7869386806, forced_plus=1.2974425414
N11-06 PASS passive observer: scalar=0.7838216659, two_channel_energy=0.5000000000
N11-07 PASS rule-specific conservation: E(k=1)=0.5; E(k=4)=2.0
N11-08 PASS finite-time directional response for chosen parameter window
N-11 CHECKS: 8/8 PASS; all synthetic model checks only, not physical experiments.

這裡的熱模型在 Dirichlet 邊界採用有限正弦級數(300 個奇數模態)計算 MDM_D ;在 t=0t=0 使用 2×10−92\times10^{-9} 截斷容差。這不是「數值解精度已達自然宇宙」,更不是連續 PDE 的嚴格浮點誤差界;若需計算機輔助正式證明,應進一步採用區間算術或譜截斷餘項的解析上界。

B.4 必需的迴歸負例

  1. 若 kc=0k_c=0 而遠端位移不為零,應標記模型或程式錯誤。
  2. 若 μ=0\mu=0 的叉式分岔公式被直接以分母 μ\mu 算到除零,應標記失敗。
  3. 若 Neumann 兩端絕熱且無內熱源,總平均熱量不守恆,應拒絕數值結果。
  4. 若 Dirichlet 固定端溫且初始場非零,長時間仍趨向 1/21/2,應拒絕模型或邊界實作。
  5. 若不同規則世界的保守能量函數都被硬要求給出同一數值,而未控制初始位移的彈性勢能,應指出比較契約錯誤。
  6. 若單純換被動觀察器造成物理能量變化,必須回查感測反作用或程式耦合錯置。
  7. 若相位資料在目標任務已充分,就不應硬標 X 或增加無必要感測器。
  8. 若新的 GRST 指標只是已知力學量的改名,應降級為語義標註,不應宣稱產生新的物理量。

附錄 C|可跨研究交接的反事實資料契約

C.1 實驗資料與模型資料的命名

本篇建議 GMDIS-12 將每項研究至少拆成五個檔案層:

  • world_spec:物理介質、設備、材料、時間座標、邊界、單位、合法狀態域、版本、測量日期。
  • intervention_spec:干預型別、改變內容、固定內容、可實施性、操作與物理狀態準備成本。
  • observation_spec:傳感器、校準、頻率、誤差、場域、尺度、壓縮與資料品質。
  • outcomes:實測或合成時間序列、能量帳本、模式與有限時間統計、反事實比較。
  • audit:使用的模型、數值算法、假設、失敗條件、基線、單元測試及評估結果。

不得把 simulation 標成 physical_measurement;不得把 explanatory_hypothesis 標成 verified_law。

C.2 JSON 範例:第 6 節兩振子的規則級反事實

{
  "schema": "GMDIS-CF-11-v0.1",
  "experiment_id": "synthetic_pair_coupling_001",
  "evidence_level": "synthetic_analytical_and_numerical",
  "world": {
    "state": "u1,u2,du1_dt,du2_dt",
    "units": {"u1": "m", "u2": "m", "du1_dt": "m/s", "du2_dt": "m/s"},
    "mass": "1 kg per oscillator (model normalization)",
    "angular_frequency": {"value": 0.9, "unit": "1/s"},
    "geometry": "fixed; two oscillator locations unchanged",
    "boundary": "unchanged; no acoustic or thermal channel in toy model"
  },
  "initial_conditions": {
    "u1": 0, "u2": 0, "du1_dt": 1, "du2_dt": 0,
    "initial_energy": {"value": 0.5, "unit": "J"}
  },
  "intervention": {
    "kind": "do_R",
    "changed_parameter": "effective coupling K/m",
    "old_value": 0,
    "new_value": 0.6,
    "unit": "1/s^2",
    "unchanged": ["initial state", "geometry", "local omega", "observables"]
  },
  "observation": {
    "signal": "u2(t)", "unit": "m", "time": {"value": 0.05, "unit": "s"},
    "kind": "synthetic exact expression"
  },
  "expected": {
    "no_coupling": 0,
    "coupled_u2": 0.0000124955943382,
    "first_nonzero_order": "t^3"
  },
  "comparator": {"normalization": "not required; output is u2 in meters", "tolerance": 1e-10},
  "labels": ["physical_shift_model_conditional", "not_a_real_experiment"],
  "outcome": "analytic_and_synthetic_consistency_only"
}

上述 1 kg1\,\mathrm{kg} 與 1 m/s1\,\mathrm{m/s} 是數值實驗的單位化示例,並不代表任意真實機械節拍器真的以 1 kg 集中質點近似。若改以轉動擺錘,請改為轉動慣量、角速度與力矩,不能把位移與擺角混用。

C.3 模型世界比較必須強制列出缺漏

當某次實驗無法控制完整初始聲場或所有接觸微結構時,請在 uncontrolled_variables、assumptions 和 uncertainty_sources 明示。若其中任一項足以解釋觀察差異,就不能報告「純 doRdo_R 已識別」。研究應改稱帶伴隨擾動的介入比較,並保留後續補證計畫。

C.4 態勢換位的建議證書

{
  "certificate_version": "GMDIS-SHIFT-11-v0.1",
  "certificate_kind": "physical_shift_or_observational_shift",
  "world_pair": ["W0", "W1"],
  "intervention_layer": "do_R | do_X | do_B | do_O",
  "alignment": {"status": "comparable", "state_bridge": "identity_on_common_state", "unit_consistency": true},
  "task_horizon": {"start": 0, "end": 5, "unit": "s"},
  "pre_registered_metrics": ["transverse_recovery_rate", "impulse_transfer", "energy_residual"],
  "thresholds": {"effect_min": "task_defined", "measurement_error": "instrument_specific"},
  "evidence_level": "model_proof | numerical_check | real_measurement | unverified",
  "decision": "shift_supported | no_detectable_shift | unresolved | not_comparable",
  "counterexamples_preserved": true,
  "reproduction_reference": "model_hash_or_experiment_record_uri"
}

這只是契約示例,不代表 JSON 欄位中的字串已被實驗或外部標準正式採用。GMDIS-12 應為實際輸入資料定義 JSON Schema、校準版本與統一 SI 單位,並在系統實作後再進行 schema validation。


附錄 D|來源與參考文獻(原始理論與外部科學分離)

D.1 EveMissLab/GMDIS 內部來源

這些文件用於延續定義及工作歷史,不能當作本文物理模型受到獨立實物驗證的外部證據。

[I-11-01] Neo.K,2026-10-05,《GRST-12:全域規則勢統合理論——從底空間到新科學方法論》,v0.1。原始來源:無界智庫公開版。支持:規則、可行自由、勢賦形、未來世界線、觀察的分層定義;不支持「新物理力已證實」。

[I-11-02] Neo.K,2026-10-04,《GRST-07:因果勢與雙層反事實——改狀態與改宇宙規則》,v0.1。原始來源:無界智庫公開版。支持: doXdo_X 、 doRdo_R 、 doBdo_B 、 doOdo_O,規則—狀態交互與相對辨識的候選語法。

[I-11-03] Neo.K,2026-10-05,《GRST-11:最小全域規則勢實驗》,v0.1。原始來源:無界智庫公開版。支持:控制同等初態、規則與資源的對照原則,並非連續物理世界的完整驗證。

[I-11-04] Neo.K,2026-10-02,《空間資訊佔據、構成纖維與容器自由度:一個主客邊—ACTM 統合框架》,v0.1。原始來源:無界智庫公開版。支持:容器、構成位、自由纖維、非局部自由陰影、表示與世界分離;其「纖維」不自動是真實物理粒子或量子場。

[I-11-05] Neo.K,2026-07-19,《計算的二十四重範式》,v4.0。原始來源:無界智庫公開版。支持:B/U/O 三軸、X 的任務相對拒單測。

[I-11-06] Neo.K × Aletheia,2026-09-17,《從二十四重形態到七十二格動力學:最低充分轉移律、混合格與計算路徑》,v0.1。原始來源:無界智庫公開版。支持:F/K/Q 轉移律及混合事件分解。

[I-11-07] GMDIS-01 至 GMDIS-10 既有研究原稿(2026-10-16 至 2026-10-17),本篇承接其數學建模、理論接縫、語義分層、計算測試與學術邊界。直接上篇檔名:GMDIS_10_Dynamical_Multi_Invariant_Profiles_and_Non_Single_Measurement_Observation_v0.1_2026-10-17.md。注意:系列原稿的人工審核與外部同行評審狀態不可因篇數增加而擅自上升。

D.2 外部成熟科學來源(已核對書目與公開摘要/書目頁)

[E-11-01] Peters, J., Janzing, D., & Schölkopf, B. (2017). Elements of Causal Inference: Foundations and Learning Algorithms. MIT Press. Publisher。參照:觀察資料、介入分布與因果可識別區分。不推導 GRST 真值。

[E-11-02] Henzinger, T. A. (1996). The Theory of Hybrid Automata. Proceedings of LICS 1996, 278–292. DOI。參照:離散事件與連續動力系統的可驗證形式框架。

[E-11-03] Brunton, S. L., Proctor, J. L., & Kutz, J. N. (2016). Discovering governing equations from data by sparse identification of nonlinear dynamical systems. PNAS, 113(15), 3932–3937. DOI。參照:資料驅動模型辨識之強基線。

[E-11-04] Raue, A., Kreutz, C., Maiwald, T., Bachmann, J., Schilling, M., Klingmüller, U., & Timmer, J. (2009). Structural and practical identifiability analysis of partially observed dynamical models by exploiting the profile likelihood. Bioinformatics, 25(15), 1923–1929. DOI。參照:結構與實務可辨識性。

[E-11-05] Sobol', I. M. (2001). Global sensitivity indices for nonlinear mathematical models and their Monte Carlo estimates. Mathematics and Computers in Simulation, 55(1–3), 271–280. DOI。參照:模型參數的全域敏感度方法;須注意 Sobol 指數為依參數分布與目標輸出定義的統計對象。

[E-11-06] Bongard, J., & Lipson, H. (2007). Automated reverse engineering of nonlinear dynamical systems. PNAS, 104(24), 9943–9948. DOI。參照:由觀察重建動力方程的歷史研究與辨識限制。

[E-11-07] Scheffer, M., Bascompte, J., Brock, W. A., et al. (2009). Early-warning signals for critical transitions. Nature, 461, 53–59. DOI。參照:分岔與臨界轉變的先兆;並非任何轉變都保證能由單一先兆預測。

[E-11-08] Holling, C. S. (1973). Resilience and Stability of Ecological Systems. Annual Review of Ecology and Systematics, 4, 1–23. DOI。參照:穩定與恢復性不同的比較視角;不是振盪器精確方程。

[E-11-09] van der Schaft, A., & Jeltsema, D. (2014). Port-Hamiltonian Systems Theory: An Introductory Overview. Foundations and Trends in Systems and Control, 1(2–3), 173–378. DOI。參照:多物理場連接與能量功率一致性。

[E-11-10] Pantaleone, J. (2002). Synchronization of metronomes. American Journal of Physics, 70(10), 992–1000. DOI。參照:共同移動平台造成節拍器同步的經典實驗與模型基線。

[E-11-11] de Silva, B. M., Champion, K., Quade, M., Loiseau, J.-C., Kutz, J. N., & Brunton, S. L. (2020). PySINDy: A Python package for the sparse identification of nonlinear dynamical systems from data. JOSS, 5(49), 2104. DOI。參照:公開可用的系統辨識工具。

[E-11-12] Kaptanoglu, A. A., de Silva, B. M., Fasel, U., et al. (2022). PySINDy: A comprehensive Python package for robust sparse system identification. JOSS, 7(69), 3994. DOI。參照:更穩健的稀疏辨識方法。

[E-11-13] Chorin, A. J., Hald, O. H., & Kupferman, R. (2002). Optimal prediction with memory. Physica D, 166(3–4), 239–257. DOI。參照:未解析自由度造成的有效記憶項,作為 GMDIS-08 的接續基線。

[E-11-14] Zwanzig, R. (1973). Nonlinear generalized Langevin equations. Journal of Statistical Physics, 9, 215–220. DOI。參照:選定熱浴模型下的廣義 Langevin 與記憶形式;不宣稱所有物理記憶都遵循該特例。

[E-11-15] Takens, F. (1981). Detecting strange attractors in turbulence. Lecture Notes in Mathematics, 898, 366–381. DOI。參照:某些觀測歷程可以重建動力系統的重要反方;不代表延遲嵌入能自動識別唯一微觀機制。

D.3 書目使用邊界與可稽核性

本稿外部資料主要以出版社、學術期刊、學校書目與 DOI 登記頁核對標題、作者、年份、頁碼和研究主題;不聲稱逐頁閱讀、重新實驗或獨立重證上述所有研究。所引文獻支持的是方法與科學基線的存在,不是 GMDIS 全部命題已獲外部支持。在投稿前應補上精確引用頁碼、與主張對應的原文證據、授權/著作權及必要的最新文獻審查。


附錄 E|獨立審稿的反方問題與下一階段品質閘門

E.1 對定理及模型的反方

Q1: 只要改正規形的 μ\mu,不就能直接得知穩定性?為什麼需要 GRST?
答: 若僅為識別該 toy model 的局部穩定性,既有分岔理論已足夠。GMDIS 的附加主張只可能存在於跨模型的型別化比較、觀察成本、因果設計與不變量追蹤;必須實測才成立。

Q2: 只要 Dirichlet 與 Neumann 邊界不同,就一定不同,這不是廢話嗎?
答: 此例有意作為模型校準與負例對照:它直接證明不能把同一 PDE 形式當作同一完整規則世界,也迫使反事實記錄保留邊界及作用帳本。它不是世界普遍新定律。

Q3: 不同耦合剛度本來就有不同遠端響應,何來創新?
答: 響應本身屬於成熟線性力學;新研究問題是以同一初態、同一幾何、同一初始能量控制條件比較規則級干預,並證明相位/幾何投影不應冒充完整動態影響。

Q4: 若 doRdo_R 可以靠增廣狀態做到,那還有獨立意義嗎?
答: 在數學語義上不一定獨立。其價值是能否降低研究中的型別混淆、錯誤實驗設計、不可比較結果和無意義的跨規則斷言;若無益,應讓位。

Q5: 所謂規則勢距離是否只是任意加權函數?
答: 本篇只給出任務相關、多不變量、有單位校準的候選距離,不宣稱唯一自然度量;若對權重極度敏感或無樣本外預測增益,該標量距離應撤回。

Q6: 觀察者本身改變系統怎麼辦?
答: 儀器負載必須進入物理模型與能源帳本,否則 doOdo_O 的理想被動假設不成立;應以混合干預報告。

Q7: 熱方程是否隱含非物理的瞬時無限傳播?
答: 擴散近似在極短時/極小尺度有適用範圍;不能把 PDE 核的數學支撐性質當作超光速物理證據。

Q8: 如果多個規則擬合同一資料,誰是對的?
答: 應保留規則等價類或信賴區間,設計具有辨識力的合法干預;不可從有限資料直接宣稱唯一真實底層規則。

E.2 投稿或提交 GMDIS-12 的必要閘門

  • 每組世界的模型、初態、邊界、外部輸入與觀測版本分開記載。
  • 反事實比較的單位、時間、尺度、符號、狀態橋接合法,並通過初始相容性檢查。
  • 物理力、功率、熱、能量與 GRST 勢指標分開記帳。
  • 所有可測結果明列感測精度、噪聲、資料缺漏與實驗重複性。
  • 近似模型、解析恆等式、合成結果、實物實驗各有獨立證據標籤。
  • 原始 24/72 分類軸及 F/K/Q 語義沒有被物理詞彙悄悄取代。
  • 凡自稱額外預測增益,均提供同成本強基線與樣本外檢驗。
  • 保存幾何零效應、純觀測換位、無法辨識與模型本身失效等負例。
  • 所有定理有精確假設;所有合成程式可從本文直接擷取再執行。
  • 研究若涉及無法物理實施的「改宇宙基本規則」,明確標為模型反事實。
  • 原始 Markdown 保留 UTF-8、完整 LaTeX 區塊及程式,不透過不可信 UI 自動改寫公式。

E.3 本篇最小交付閉包

數學閉包: 四組模型內解析命題與導出邊界清楚;至少八項合成數值檢查可重現。
語義閉包: 可比初態、四層干預、世界—模型—求解—感測分離、GRST/SIO 與 24/72 的接口有明確型別。
物理閉包: 不宣稱新物理守恆定律、不把觀測切換誤認真實能量流、所有實物未測條件有註記。
後續閉包: GMDIS-12 接手可執行節拍器矩陣世界;GMDIS-13 接手盲測、消融、基線與判敗;GMDIS-14 僅依已獲支持證據統合。


文件最終聲明

本研究是研究性完整原稿,不是已完成自然界統一理論的宣告。任何不同法則之比較都依賴明確的模型、合法干預、觀察條件與時間尺度。文內數學恆等式可被獨立重算,附錄合成程式可被重新執行,真正的 100/1000 節拍器實物研究及更複雜多場模型驗證仍屬後續工作。若成熟既有方法能以更簡潔的方式完全吸收某項 GRST/SIO 實作,應主動記錄並讓位。