title: GMDIS-07｜持續作用與動態穩定結構：驅動、耗散、橫向穩定與集體模態選擇
english_title: Persistent Actions and Dynamically Stable Structures: Driving, Dissipation, Transverse Stability, and Collective Mode Selection
series: GMDIS — Global Multiscale Dynamics and Informational Sufficiency Theory
paper: 07 / 14
version: v0.1
date: 2026-10-17
author: Neo.K
research_collaboration: Aletheia / GPT-6
institution: EveMissLab / 一言諾科技有限公司
locale: zh-TW
status: 完整研究原稿；模型條件下的解析結果與數值驗證；真實節拍器實驗待進行
canonical_format: UTF-8 Markdown; mathematical delimiters $...$ and $$...$$ only
predecessors: ["GMDIS-01", "GMDIS-02", "GMDIS-03", "GMDIS-04", "GMDIS-05", "GMDIS-06"]
next_paper: GMDIS-08｜跨尺度投影、記憶與自由度傳輸
license_note: 研究原稿；授權方式待權利人決定

GMDIS-07｜持續作用與動態穩定結構

Q: 4.5 C-07-02：改變實驗條件，結論會如何失效？

1. 移除阻尼 $c=0$。 相對模態一般持續振盪，不會因共同驅動而自動衰減；共同模態若恰在共振也可能無界地隨時間增長（在理想線性模型中）。 2. 改為異質個體。 若兩者的質量、剛度或阻尼不同，相同外力不必產生相同長期振幅與相位。 3. 改為相反外力。 當兩端輸入方向相反，差異方程獲得非零受迫項，相對差異一般不會趨零。 4. 加入未解析聲場。 若反射聲壓回饋包含延遲或頻率相依阻抗，原本的單一 $\kappa$ 不能無證據地保證仍為純正彈性耦合。 5. 加入擒縱重設。 非光滑速度重設可能改變穩定性，需以混合系統事件映射重算。

Q: 9.3 為何不能按「吸引最強」排出普遍單一勝者？

至少存在三個不可混同的排序標準：

驅動、耗散、橫向穩定與集體模態選擇

Persistent Actions and Dynamically Stable Structures: Driving, Dissipation, Transverse Stability, and Collective Mode Selection

作者： Neo.K
AI 研究協作： Aletheia / GPT-6
機構： EveMissLab（一言諾科技有限公司）
系列： 全域多尺度動力學與資訊充分性理論（GMDIS）07／14
版本與日期： v0.1，2026 年 10 月 17 日（Asia/Taipei）
文件性質： 非線性動力學、受驅動耗散系統、耦合振盪、穩定性理論、跨尺度物理語義、反事實研究規約。
證據地位： 本文給出限定模型下可直接檢查的公式、解析命題與有限精度數值測試；尚未以實物節拍器、聲腔或流固耦合裝置驗證 GMDIS 的任何額外預測增益。不宣稱發現新基本力、普遍收斂律、普遍同步律或新增熱力學守恆定律。

日期： 2026-10-17
版本： v0.1

摘要

本篇研究的核心問題不是為何多個節拍器「最後看起來一樣」，而是：在同時存在外力、機械約束、持續供能、內部回饋、材料與流體損耗、介質儲能、共振、幾何耦合及隨機擾動的完整世界中，哪些動態作用得以持續存在？它們如何改變不同集體運動模式的可達性、恢復速率、吸引域和長期存續結構？這項提問承接 GMDIS-01 的世界生成觀、GMDIS-02 的投影資訊充分性、GMDIS-03 的 SIO–GRST 接縫、GMDIS-04 的力學型別校準、GMDIS-05 的演化算子與能量適定性，以及 GMDIS-06 的跨介質能量帳本。

本文首先進行七項不可混淆的分離：持續作用與主導作用、強力與高有效增益、物理能量平衡與動力穩定性、穩定平衡點與穩定週期軌道、共同外力鎖定與自主振盪器相互同步、吸引子存在與其吸引域大小，以及有限時間擾動放大與無限時間漸近穩定。以物理能量收支 $\dot H=P_{\mathrm{in}}-D-P_{\mathrm{out}}$ 為起點，指出非零的長期穩定運動必須在所選系統邊界下具有相容的能量來源、交換與耗散；長期平均收支相抵並不保證軌道穩定，對穩定軌道的證明仍須檢查橫向動力、Floquet 乘子、Lyapunov 泛函或其他合適判據。

其次，本文建立兩個可分析的互補世界。第一個是共同受迫雙振子：兩個相同阻尼振盪器由共同週期外力驅動，彼此相對位移遵守自治阻尼方程，能精確證明相對運動衰減，同時共同運動仍維持非零週期響應。此例的意義是顯示「差異收束」與「能量持續流動」相容；但其同步也可能由共同外部驅動造成，不能冒充為自主節拍器的自發同步。第二個是自維持 Stuart–Landau 振盪器網路：在指定耦合形式下，可解析證明正耦合時同向集體極限環的局部軌道穩定、零耦合時相位差不被自動消除，以及兩振盪器在負耦合下反相極限環的局部穩定。這直接顯示持續供能不等於必然同步、穩定不等於只有一種相位結構、相同組件也不等於相同未來。

第三，本文以 Floquet／橫向變分、圖拉普拉斯譜、介面能量流與有限時窗受控響應，建立「持續性—耦合—模態選擇」的最低操作剖面，並將其與 GRST 的未來非均勻性、SIO 的動態影響陰影對接。所謂規則勢仍是規則所誘導的可行未來空間之結構性非均勻，不是焦耳或牛頓單位的全新勢能。本文保留 24／72 作為局部模型的計算與觀察分類，不將其稱為 72 種自然力或 72 種穩定態。

最後，本文提出多階模型與實驗矩陣，要求將共同底座、外部週期驅動、擒縱機制、聲學回饋、材料阻尼、幾何旋轉及隨機噪聲分別增模、消融和對照；明列 GRST／SIO 特徵應如何與傳統相位—振幅、完整 Newton／FSI、master stability function、Koopman、Mori–Zwanzig 與機器學習預測基線公平競爭。若新增語義剖面在可重現的干預預測、穩定模態分類或恢復時間估計上沒有增量，則其獨立物理價值應被降低，而非藉由詮釋語言保留必勝地位。

關鍵詞： 持續作用、動態穩定、受驅動耗散、週期軌道、orbital stability、Floquet、橫向恢復、Stuart–Landau、同步、反相、多穩態、共同強迫、聲學回饋、力學能量帳本、GRST、SIO、規則勢、資訊充分性、24／72 分類。

0. 本篇的研究責任、輸入契約與證據分級

0.1 為什麼第七篇必須獨立撰寫？

GMDIS-06 已回答一個物理一致性問題：在指定聲腔、集中振子及熱節點的簡化模型中，機械—聲學的界面功率如何相消，耗散如何被正確記在熱帳本，以及數值時間積分如何不憑空增減能量。但一份正確的能量帳本尚未判定哪種集體運動會吸引鄰近軌道。原因是能量是一個標量或少量儲能量的族，而穩定性可以依賴時變 Jacobian、非正規暫態、記憶、吸引域、邊界事件與整體耦合拓撲。

相反地，一個漂亮的相位鎖定動畫即使顯示 100 或 1,000 台振盪器一致擺動，也還不能回答驅動功率從何而來、聲學與材料的損耗如何處理、穩定性是否對微小擾動有恢復能力、或將底座固定後會發生什麼反事實。

因此 GMDIS-07 必須在能量一致性的前提下，補上長時間穩定性、模態選擇與干預敏感度的可檢查數學語義，並在 GMDIS-08 之前標出無法丟失的隱藏自由度與記憶義務。

0.2 從前六篇接收的契約

上游文件	已提供的層	GMDIS-07 的守界義務
GMDIS-01	世界、規則、觀察、反事實及系列判敗	不把同步指定為宇宙必然終點
GMDIS-02	精確／近似資訊充分性、相位與吸引子強基線	不以低維度本身證明低資訊充分性
GMDIS-03	SIO 的硬約束陰影／動態影響陰影，GRST 勢賦形	不將局部相互作用直接命名為新物理勢能
GMDIS-04	力、應力、功率、作用、儲能與耗散的型別	不將「力量持續」偷換成「做功持續」
GMDIS-05	多場能量空間、演化算子、適定性與非正規暫態	不把有限模型穩定定理外推到任意 PDE
GMDIS-06	跨介質界面、熱帳本、能量相容的離散模型	不重複計算聲學擾動能與完整流體內能

本篇對外延續 GMDIS-01–06 的同一檔案系列。若基準模型採用的是理想化無量綱正規型，必須明示它不是具 SI 單位的節拍器裝置方程。

0.3 四類證據標記

D-07（Definition）：在指定語境下可操作的定義，不自動包含現實自然律的主張。
M-07（Model theorem）：在明示模型與假設下成立的可驗證數學命題，附證明或推導。
N-07（Numerical consistency）：對明確給定的方程、參數與數值方法完成有限精度檢查；不是實物驗證。
H-07（Hypothesis）：待模擬、實驗或異質模型競爭檢驗的預測與方法論主張；F-07 是對應判敗條件。

若某個命題的證明只在局部線性化成立，不能在摘要或結論中升格為所有初態、所有參數、所有物理介質下都成立的全域結論。

0.4 一句話的研究主張

長時間作用所產生的集體秩序，不由力的大小或存在時間單獨決定；它取決於驅動、耦合、耗散、邊界與時空結構共同塑造的穩定演化集合。

此處的「秩序」是一個可分析的動力學結果，不含物理世界具有意圖、目的或必然朝同步演化之意。

1. 七組必須拆開的概念

1.1 持續存在的作用不等於主導作用

一個恆定向下的重力在物體受支持而靜止時，可以由支持力平衡，平均機械做功為零；一個振幅很小的週期激勵，若在高品質因數共振附近反覆輸入能量，卻可能造成大響應。因此作用力持續不代表累積做功為正，作用力的峰值大也不表示其對目標模式的有效增益最大。

1.2 力的平衡不等於能量交換停止

力的向量和可能為零，但界面上仍可能存在不同位置的局部功率流；反之，一個未受淨外力的系統也可以具有內部聲波、相對運動或熱交換。力、動量、能量與熱量必須逐項列帳。

1.3 能量長期平衡不等於動力穩定

假設一條週期軌道每個週期輸入功與總耗散相等，這只是該週期狀態的必要收支要求。它沒有告訴我們附近偏離軌道的擾動會縮小還是被放大；後者是 Floquet 乘子、橫向 Lyapunov 指數或適當能量—變分判據要回答的問題。

1.4 穩定平衡點不等於穩定週期軌道

穩定平衡點可以隨時間保持不變；穩定極限環則可以持續以非零速度運動。對自治週期軌道，沿軌道相位平移方向通常具有中性擾動，不能把這個中性方向誤判成軌道不穩定。相應的穩定概念是模掉共同相位後的軌道穩定，不是所有座標差對固定時間參考都必須趨零。

1.5 共同強迫鎖定不等於自主集體同步

若兩個無耦合、完全相同且受同一正弦外力驅動的阻尼振盪器，穩態可以有完全相同的位移與相位。這說明的是對共同外界時鐘的鎖定，不是兩者透過彼此交換作用而產生的同步。真實節拍器通常具有自己的擒縱供能；共同底座耦合可使不同自治驅動的振盪形成集體關係。實驗必須分開這兩種因果機制。

1.6 吸引子存在不等於吸引域巨大

一個模式可以局部穩定但吸引域很小，僅在精密製備的初態附近可觀察；另一種模式可能在相同參數下吸引較大比例的初態。因此以影片中最終看到哪個模式作為唯一判準，無法得知其他模式、轉換閾值或其吸引域大小。

1.7 長期漸近穩定不等於所有有限時間安全

非正規系統的演化算子可能在其特徵值實部均為負時產生明顯的暫態放大。若研究關注裝置碰撞、聲學峰值或短時間結構破壞，只問 $t\to\infty$ 是否穩定便不足；還應估計特定能量內積下的短時增益。

以下表格歸納研究中最常出現的偷換。

被混合的概念	不可直接等價的理由	最小檢查
持續力／持續功率	力與速度夾角及時間相關性重要	測 $F\cdot v$ 而非僅測 $
能量平衡／穩定性	收支平衡不判定橫向恢復	線性化／Floquet／Lyapunov
個體極限環／群體相位鎖定	各自自維持不等於相互鎖定	耦合移除試驗
同向／全狀態相同	相位相同仍可有振幅與內部模態差異	同步流形距離
相位／完整因果機制	投影可能丟失介質和驅動記憶	相同投影不同未來反例
低波動／高韌性	受測擾動方向可能沒有被激發	系統化注入脈衝
格點相同／耦合相同	底座旋轉、彈性與聲場改變算子	幾何負向基線

2. 完整世界與動力穩定性的基本型別

2.1 D-07-01：帶外部埠的世界演化

固定一個研究語境 $\Gamma$ ，將物理世界的有效狀態寫成：

X(t)=\bigl(q,\dot q,\xi_s,p_a,v_a,T_s,T_a,e,\zeta\bigr)(t).

其中 $q$ 為集中機械位移／擺角， $\xi_s$ 為固體形變場， $p_a,v_a$ 為所選聲學近似中的壓力與速度， $T_s,T_a$ 為溫度場， $e$ 是擒縱或驅動狀態， $\zeta$ 為其他需保留的邊界或內部狀態。並非每個實驗都必須使用所有分量；是否可以刪去某分量需由任務充分性及誤差界決定。

正式寫成受控或混合演化：

\dot X=\mathcal A_R(X)+\mathcal B_R(X)u(t),

並容許在合法事件面 $h_j(X)=0$ 發生：

X^+=\mathcal R_j(X^-).

這裡 $\mathcal A_R$ 是在指定規則、介面與適定性假設下的演化算子， $u$ 是已型別化的外部輸入， $\mathcal R_j$ 是重設算子。無限維 PDE、非光滑接觸或隨機演化不能一律假裝為普通全域 Lipschitz 常微分方程。

2.2 D-07-02：能量—功率帳本

若將系統邊界劃定為機械、聲學與集總熱節點，且外部驅動從邊界進入，則採一般記帳形式：

\dot E_{\mathrm{tot}} =P_{\mathrm{drive}}-Q_{\mathrm{out}}-P_{\mathrm{rad}}+P_{\mathrm{other}}.

各子系統則可以寫為：

\dot E_m=P_{\mathrm{drive}}-D_m-P_{m\to a},

\dot E_a=P_{m\to a}-D_a-P_{a\to\mathrm{rad}},

\dot E_{th}=D_m+D_a-Q_{\mathrm{out}},

其中 $D_m,D_a\ge0$ 是在所採用模型中被解釋為轉入熱節點的功率。相加得：

\dot E_{\mathrm{tot}}=P_{\mathrm{drive}}-Q_{\mathrm{out}}-P_{a\to\mathrm{rad}}.

若某些材料阻尼所轉移的能量其實離開所選容器，而未進入熱節點，則必須改變帳本而不是悄悄把它加進 $E_{th}$ 。若將供能彈簧、電池或發條納入同一封閉世界， $P_{\mathrm{drive}}$ 應改列為內部能量轉移並追蹤其儲能減少，不能一面視為外部供給、一面又在內部重複計算。

2.3 M-07-01：有限週期運動的平均能量必要條件

命題。 若某系統在 $[t,t+T]$ 具有週期 $T>0$ 的完整能量狀態，且上述各功率可積，則：

\frac1T\int_t^{t+T}P_{\mathrm{drive}}\,ds = \frac1T\int_t^{t+T}(Q_{\mathrm{out}}+P_{a\to\mathrm{rad}})\,ds.

證明。 沿一完整週期對總能量方程積分；由 $E_{\mathrm{tot}}(t+T)=E_{\mathrm{tot}}(t)$ ，左側儲能變化為零，故得等式。證畢。

界線。 此命題為指定能量帳本下的必要條件，不是軌道穩定性的充分條件。若外界熱節點也包含在總系統內、溫度持續上升，或驅動所用發條尚未補充能量，則完整狀態不是真正的週期軌道，不能無條件引用。

2.4 D-07-03：動態模式集合與觀察目標

令 $\mathcal M_R\subseteq\mathcal X$ 表示在規則 $R$ 下想研究的不變集合，它可能是平衡點、週期軌道、同步流形、反相流形或更一般吸引集合。對任務距離 $d_\Gamma$ ，考慮：

D_{\mathcal M}(X)=\inf_{Y\in\mathcal M_R}d_\Gamma(X,Y).

若存在鄰域 $\mathcal U$ ，使對 $X_0\in\mathcal U$ 的合法解滿足 $D_{\mathcal M}(X_t)\to0$ ，且其偏離在全時段保持受控，才可以討論局部吸引與穩定。週期軌道的距離應對整條軌道取下確界，不能錯把相位平移造成的差異當作橫向失穩。

此處的 $d_\Gamma$ 可以是無量綱化能量範數、Hilbert 空間距離、可觀測量差異或其他任務相關距離，但若比較模型必須固定型別、尺度及量綱轉換，不能臨時替換最有利的距離。

3. 持續作用的數學形式：力強、作用久，還缺哪些條件？

3.1 D-07-04：有限時窗的因果傳遞算子

對一條合法參考解 $X_\ast(t)$ ，在可微模型且線性化有意義時，令 $\eta$ 表示偏離，並寫：

\dot\eta=A(t)\eta+B_j(t)u_j(t),

其中 $A(t)$ 是沿參考解的變分算子， $u_j$ 是第 $j$ 類受控作用。令 $\Phi(t,s)$ 為對應的狀態轉移算子，則初始偏離為零時：

\eta(T)=\int_0^T\Phi(T,s)B_j(s)u_j(s)\,ds.

對給定觀察算子 $C$ ：

\delta y(T)=C\int_0^T\Phi(T,s)B_j(s)u_j(s)\,ds.

這裡的整個被積核，而不僅是輸入力的瞬時大小，才決定該作用對指定未來觀察量的線性化影響。 作用持續時間、相位與模態方向、系統回授、延遲、邊界與阻尼都藏在 $B_j$ 、 $\Phi$ 與 $C$ 的組合之中。

可提出型別一致的受控增益：

G_{j,\Gamma}(T) =\sup_{\|u_j\|_{U,2}\le1}\|\delta y(T)\|_{Y,\Gamma}.

此增益是指定輸入／輸出規約下的算子範數，不是焦耳單位的能量或全域規則勢的獨立物理實體；如不同輸入有不同 SI 單位，必須先加入明確的尺度轉換或分別報告。

3.2 M-07-02：大而恆定的力可以對穩態無持續功率

令彈簧平衡方程為：

m\ddot q+c\dot q+kq=F_0,

其中 $m,k,c>0$ ， $F_0$ 恆定。唯一平衡點是 $q_\ast=F_0/k$ ， $\dot q_\ast=0$ ，因此：

P_\ast=F_0\dot q_\ast=0.

恆定力確實改變平衡構形並可在過渡過程做功，但到達靜態平衡後沒有持續機械功率輸入。這是對「長期強作用必然累積為主導運動」的最小解析反例。

3.3 C-07-01：小驅動的共振可以產生大響應

令無量綱阻尼振盪器滿足：

\ddot q+2\zeta\omega_0\dot q+\omega_0^2q=\varepsilon\cos(\omega_0t),

其中 $\zeta>0$ 。在精確共振且線性近似成立時，穩態振幅是：

A_{\mathrm{res}}=\frac{\varepsilon}{2\zeta\omega_0^2}.

若 $\omega_0=1$ 、 $\varepsilon=0.01$ 、 $\zeta=0.001$ ，則形式上的 $A_{\mathrm{res}}=5$ 。數值只反映線性模型的高增益；真實裝置可能先觸及撞擊、非線性行程或材料邊界，不得外推至無限振幅。這說明即使驅動很小，只要它長期與能吸收其能量的模式有效配合，就可能有顯著後果。

3.4 不是任何「抵抗」都會耗散

彈性回復力可暫存能量並於後段釋回；理想聲腔也可以透過反射將聲學能回送給活塞。把兩者通稱為「小力被大力吸收」會錯失儲能與回送。真正不可逆轉熱、摩擦耗散或向外輻射，應由有物理量綱與正定性的耗散率或外流功率表示。

3.5 橫向穩定才是長期模態選擇的關鍵之一

對給定集體模式，真正決定鄰近擾動能否衰減的是沿運動歷程的變分動力。某種力可以很大，卻對與同步流形正交的差異方向沒有任何恢復作用；另一個較弱但專門消除相對偏移的耦合，則可能使同步流形穩定。要避免把「能量大」直接推出「吸引子大」或「恢復快」。

4. 可證模型 A：共同週期驅動的雙振子

4.1 D-07-05：模型、參數及物理地位

考慮相同質量 $m>0$ 、自身彈簧剛度 $k>0$ 、黏性阻尼 $c>0$ 與非負互連剛度 $\kappa\ge0$ 的兩個線性振盪器。兩者受到相同的外部週期力 $F\cos(\Omega t)$ ：

m\ddot q_1+c\dot q_1+kq_1+\kappa(q_1-q_2)=F\cos(\Omega t),

m\ddot q_2+c\dot q_2+kq_2+\kappa(q_2-q_1)=F\cos(\Omega t).

這是一個具 SI 單位的簡化受迫機械系統，可用作能量與穩定性命題的透明基線。它不是具有擒縱機構的自維持節拍器模型，也不是用來解釋實物共同底座上所有可能模式。將物理場消元成有效的 $\kappa$ 需要額外的力學推導，而不能因為雙振子可同步就推回聲場／底座的真實耦合係數。

設平均座標與相對座標：

a=\frac{q_1+q_2}{2},\qquad e=q_1-q_2.

將兩式相加、相減可得精確分解：

m\ddot a+c\dot a+ka=F\cos(\Omega t),

m\ddot e+c\dot e+(k+2\kappa)e=0.

其中 $a$ 是共同運動， $e$ 是各振盪器的相對差異。兩者具有不同的演化責任：共同模態持續接受能量，相對模態在指定正阻尼下衰減。

4.2 M-07-03：相對運動指數收斂，非零共同週期運動保留

命題 M-07-03。 假設 $m,c,k>0$ 、 $\kappa\ge0$ 、 $F\neq0$ 、 $\Omega>0$ 。則任意有限初態的相對運動 $e(t),\dot e(t)$ 均以指數速率趨於零；平均運動 $a(t)$ 趨於唯一的受迫週期解 $a_\ast(t)$ ，其振幅為：

A=\frac{|F|}{\sqrt{(k-m\Omega^2)^2+(c\Omega)^2}}.

因此兩者的相對位置與速度差會衰減，同時共有的長期運動仍可保持非零振幅。

證明。 相對運動的特徵多項式為：

m r^2+c r+k+2\kappa=0.

由 $m,c,k+2\kappa>0$ 可知兩個特徵根的實部嚴格為負；若為共軛複根，其實部為 $-c/(2m)$ ，若為實根，兩根均為負。故相對模態指數衰減。共同模態的齊次特徵根亦有負實部，所以所有初態的過渡項衰減，留下頻率 $\Omega$ 的特解。以複指數代入 $a_\ast=\operatorname{Re}(\widehat a e^{i\Omega t})$ ，得：

\widehat a=\frac{F}{k-m\Omega^2+i c\Omega},

振幅即上式。證畢。

重要但容易忽略的負向基線： 若設定 $\kappa=0$ ，命題仍可成立。原因是兩個相同阻尼振盪器受同一週期外力驅動，會各自收斂到同一特解。這不能被當作互連耦合的直接證據。要測自主集體同步，必須另設 $F=0$ 且具有個別自維持供能的模型，或分離共同外力與彼此耦合的獨立作用。

4.3 D-07-06：精確機械能

定義：

H_m=\frac m2(\dot q_1^2+\dot q_2^2) +\frac k2(q_1^2+q_2^2) +\frac\kappa2(q_1-q_2)^2.

此儲能單位為焦耳， $\kappa(q_1-q_2)^2/2$ 是互連彈簧的能量，並非 GRST 的規則勢。

對運動方程各乘相應速度並相加，可得：

\dot H_m =F\cos(\Omega t)(\dot q_1+\dot q_2) -c(\dot q_1^2+\dot q_2^2).

上式沒有顯式 $\kappa$ 功率項，是因為其內部彈簧功率已計入 $H_m$ 的導數，並非耦合不作用。對長期共同週期解， $q_1=q_2=a_\ast$ ，一週期平均帳本為：

\langle P_{\mathrm{in}}\rangle =\langle D\rangle =c\Omega^2A^2>0.

所以穩定運動不等於沒有耗散；它代表受迫振動在指定邊界下持續獲得能量並抵銷平均耗散。

4.4 相對能量與恢復速率不能混入總能量

對差異坐標可定義：

E_e=\frac m2\dot e^2+\frac{k+2\kappa}{2}e^2.

由相對運動方程得到：

\dot E_e=-c\dot e^2\le0.

此非嚴格耗散不等式配合正定剛度與阻尼足以導出相對運動漸近收斂；也可直接由特徵根證得指數收斂。但請注意， $E_e$ 是差異坐標的候選 Lyapunov 函數，並不是整體物理能量 $H_m$ 。它的零點並不表示兩個振盪器停止運動，只表示它們的相對位置與速度差為零。

一個符合 GMDIS 研究目標的報告，應同時呈現 $H_m$ 、 $E_e$ 、相對運動回復率、共同行程與功率，而非只畫「相位差趨近零」。

4.5 C-07-02：改變實驗條件，結論會如何失效？

移除阻尼 $c=0$ 。 相對模態一般持續振盪，不會因共同驅動而自動衰減；共同模態若恰在共振也可能無界地隨時間增長（在理想線性模型中）。
改為異質個體。 若兩者的質量、剛度或阻尼不同，相同外力不必產生相同長期振幅與相位。
改為相反外力。 當兩端輸入方向相反，差異方程獲得非零受迫項，相對差異一般不會趨零。
加入未解析聲場。 若反射聲壓回饋包含延遲或頻率相依阻抗，原本的單一 $\kappa$ 不能無證據地保證仍為純正彈性耦合。
加入擒縱重設。 非光滑速度重設可能改變穩定性，需以混合系統事件映射重算。

此例可證明「某個模式能穩定保持」的最小數學版本，也清楚畫出它不能解釋的物理世界。

5. 可證模型 B：自維持 Stuart–Landau 網路與局部同步

5.1 D-07-07：為什麼要第二個模型？

第一個模型之非零長期運動來自共同外力。要區分「個體自己能持續振盪」與「個體之間能穩定協同」，我們另外使用典型超臨界 Hopf 正規型：

\dot z_j=(\lambda+i\omega-|z_j|^2)z_j +K(\bar z-z_j),

其中：

\bar z=\frac1N\sum_{j=1}^N z_j, \qquad \lambda>0, \quad K\in\mathbb R.

$z_j\in\mathbb C$ 是二維無量綱正規型狀態， $\lambda$ 是自由振幅成長參數， $\omega$ 是內部旋轉頻率， $K$ 是平均場差分耦合係數。 $K>0$ 表示朝群體平均態的吸引型差分耦合， $K<0$ 表示排斥型差分耦合； $K$ 不是牛頓或焦耳單位的基本力，不可直接等同節拍器底座的彈簧、摩擦或聲學耦合係數。

本模型具有非線性振幅飽和：不耦合時，每個 $z_j$ 的振幅趨近 $\sqrt\lambda$ ，相位持續以角速度 $\omega$ 轉動。可將「負阻尼／主動補償直到飽和」作為直覺類比，但在此正規型裡不能把 $|z_j|^2$ 未經轉換直接稱為真實機械能。

5.2 M-07-04：正耦合下全同向極限環的局部軌道穩定

命題。 對 $N\ge2$ 、 $\lambda>0$ 、 $K>0$ 的上述正規型，存在同向集體週期解：

z_j^\ast(t)=\sqrt\lambda\,e^{i(\omega t+\phi_0)}, \qquad j=1,\ldots,N.

其在模掉共同相位平移後對充分小擾動局部漸近穩定。具體地：共同徑向特徵指數為 $-2\lambda$ ；每個與平均場正交的橫向方向，徑向和相角指數分別為：

\gamma_{\perp,r}=-(2\lambda+K), \qquad \gamma_{\perp,\phi}=-K.

共同相位方向的指數為零。

證明。 令 $w_j(t)=z_j(t)e^{-i\omega t}$ ，則旋轉座標下：

\dot w_j=(\lambda-|w_j|^2)w_j +K(\bar w-w_j).

同向解變成平衡 $w_j^\ast=\sqrt\lambda$ 。令 $w_j=\sqrt\lambda+u_j+i v_j$ ，取一階變分，得：

\dot u_j=-2\lambda u_j+K(\bar u-u_j),

\dot v_j=K(\bar v-v_j).

將 $u,v$ 分解成所有分量相同的平均方向與零均值的橫向方向。平均徑向得到 $-2\lambda$ ，平均相角得到 $0$ ；零均值徑向得到 $-(2\lambda+K)$ ，零均值相角得到 $-K$ 。假設 $\lambda,K>0$ ，除共同相位中性方向外，其餘變分指數均為負，由光滑系統的局部正常雙曲週期軌道理論得到局部軌道吸引。證畢。

詮釋。 如果只觀察振幅接近 $\sqrt\lambda$ ，無法知道彼此相位差是否收斂；如果只觀察共同振盪頻率，也無法知道徑向擾動以什麼速率恢復。這兩種資訊由不同的變分方向承載。

5.3 C-07-03：零耦合的正面反例——自維持不必然同步

令 $K=0$ 。每個振盪器仍具有自己的穩定振幅 $\sqrt\lambda$ ，但各自的相位平移是中性方向：

z_j(t)\sim\sqrt\lambda\,e^{i(\omega t+\phi_j)}.

若初態具有不相同的 $\phi_j$ ，差異通常可以永久保留。個體的自維持極限環，並沒有自動生成所有個體共享的單一相位。 此反例與「大量個體長時間都在做一樣的事，最終自然會同步」的無條件敘述直接衝突。

在 GMDIS 語言下，外觀相似、能量流長期存在，只是相應世界的局部描述；要判斷集體模態，仍需檢查耦合結構、橫向恢復與共同環境。

5.4 M-07-05：負耦合可使兩個振盪器的反相週期軌道局部穩定

命題。 對 $N=2$ 、 $\lambda>0$ 、 $K<0$ 的同一模型，存在反相週期解：

z_1^\ast(t)=-z_2^\ast(t) =r\,e^{i(\omega t+\phi_0)},

其中：

r=\sqrt{\lambda-K}.

它在模掉整條反相週期軌道的共同相位後局部軌道穩定。其旋轉座標下的一階特徵指數為：

0,\qquad -2(\lambda-K),\qquad K,\qquad -2\lambda+3K.

由 $\lambda>0$ 、 $K<0$ ，除零的相位指數外均嚴格為負。

證明概要。 在旋轉座標令 $w_1=r$ 、 $w_2=-r$ 。此時 $\bar w=0$ ，每個單位的耦合項為 $-Kw_j$ ；由 $(\lambda-r^2-K)r=0$ 得 $r^2=\lambda-K$ 。設擾動 $\delta w_1,\delta w_2$ ，分解為對稱 $(\delta w_1+\delta w_2)/2$ 與反對稱 $(\delta w_1-\delta w_2)/2$ ，各自再分徑向與切向，得到上述四個特徵指數。其符號保證局部軌道穩定。完整線性化見附錄 B。證畢。

限制。 這是所選兩個複數正規型的解析結果，不是「負的機械摩擦必然反相」或「所有節拍器都只可能同向／反相」的定律。真實機械排斥耦合如何出現，需由底座模型、擒縱與物理接口先行確立。

5.5 這三個模型結果並列的真正意義

相同的個體形式	耦合條件	局部長期運動	正確結論
$\lambda>0$ 、同一 $\omega$	$K>0$	同向週期軌道局部穩定	同向是一個可吸引模式
$\lambda>0$ 、同一 $\omega$	$K=0$	個體振幅穩定、相位差可保留	持續振盪不必同步
$N=2$ 、 $\lambda>0$	$K<0$	反相週期軌道局部穩定	集體穩定性可選擇不同模態

在三個世界裡，個體都可以長期保有振盪能力；但在不同耦合規則下，其集體未來完全不等價。這是 GRST「相同自由能力不保證相同未來結構」非常直接的模型級示例。若有人能用傳統 Hopf 正規型、耦合拓撲和橫向譜完全解釋這些結果，我們也必須承認它們本身已能解釋，不應強行宣稱 GRST 在此提供新的物理預測。

6. 從二元到網路：拓撲如何改變穩定性的時間尺度？

6.1 D-07-08：一般連通網路版本

考慮 $N$ 個全同 Stuart–Landau 振盪器，在無向、非負邊權重圖上以拉普拉斯矩陣 $L$ 耦合：

\dot z_j=(\lambda+i\omega-|z_j|^2)z_j -K\sum_{k=1}^N L_{jk}z_k.

令特徵值滿足：

0=\mu_1\le\mu_2\le\cdots\le\mu_N.

對連通的無向圖， $\mu_2>0$ 。由 $L\mathbf1=0$ ，全同週期軌道依然存在；在其旋轉座標下沿第 $\ell$ 個拉普拉斯本徵模態的徑向與相角指數為：

\gamma_{\ell,r}=-2\lambda-K\mu_\ell,

\gamma_{\ell,\phi}=-K\mu_\ell.

其中 $\ell=1$ 對應共同模式；若 $K>0$ 且圖連通，所有 $\ell\ge2$ 的橫向方向均嚴格穩定。證明就是將第 5 節的橫向變分方程沿 $L$ 的正交本徵基底分解，無須再另立一個普遍同步公理。

6.2 M-07-06：連通度不只改變「能否」同步，也可能改變「多久」

上述模型的最慢橫向相角衰減率為 $K\mu_2$ ；在局部線性化有效的條件下，初始相角偏離約按：

\delta\phi_\perp(t) \approx e^{-K\mu_2t}\delta\phi_\perp(0)

的最慢模態衰減。此處 $K\mu_2$ 具有時間倒數的量綱（在無量綱模型中依所採時間單位決定），是實際可比較的局部恢復率。這個定量結果遠比「有很多節拍器正在互相影響，所以勢很強」更可檢驗。

但這不表示代數連通度單獨決定所有穩定性：有向、時變、延遲、非線性非加性、高階耦合、異質頻率、相位—振幅交互及底座旋轉模態，都可能要求更一般的主穩定函數、Floquet、非正規或非線性分析。

6.3 多尺度幾何與「不要虛構位置效應」

GMDIS-03 已指出：若共同底座只作一維剛體平移，且各振盪器質量、擺向和耦合方式理想同一，單純將振盪器重新排成不同格點，不應憑空改變模型。要觀察排列幾何的物理效應，必須允許底座轉動、柔性彎曲、聲學路徑、接觸分布或其他確實由位置決定的項進入演化算子。

因此研究 $10\times10$ 與 $1\times100$ 的不同時，必須先寫出它們在模型中改變了哪個矩陣、邊界或介面條件。如果只是重新排列圖形，沒有任何方程係數改變，卻得到不同動力結果，那反而可能是程式把繪圖座標錯誤混入物理更新的證據。

6.4 已有的強基線：Master Stability Function

Pecora 與 Carroll（1998）已提出針對耦合系統同步穩定性的主穩定函數方法；更近期非線性耦合 Stuart–Landau 網路的研究亦結合頻譜與 Floquet 分析。因此 GMDIS 不得將「耦合拓撲決定某種橫向穩定」當成自己首次發現的物理定理。新系列所能再進一步研究的是：跨物理介質的可追溯耦合、觀察不足的反事實損失，以及 GRST／SIO 特徵是否提供超越既有強基線的增量任務表現。

7. 吸引子、Floquet、橫向 Lyapunov 與長期存續的精確邊界

7.1 D-07-09：週期軌道之線性變分

對平滑自治系統：

\dot x=f(x),

假設存在週期軌道 $x_\ast(t)$ ，週期 $T>0$ 。沿軌道的變分滿足：

\dot\eta=A(t)\eta, \qquad A(t)=Df(x_\ast(t)), \qquad A(t+T)=A(t).

令基本解矩陣 $\Phi(0)=I$ ，則單週期單值映射的線性部分（monodromy matrix）為 $M=\Phi(T)$ ；其特徵值稱 Floquet 乘子。對自治週期解，由沿流向的時間平移對稱性，通常存在一個乘子 $1$ 。若這個單位乘子為簡單乘子，且其餘乘子的模均小於 $1$ ，則週期軌道在適當局部光滑假設下軌道漸近穩定。

因此不能草率要求所有乘子都嚴格小於 $1$ ，那會把正常的共同相位中性方向誤判成不穩定。

7.2 D-07-10：橫向指數、共同相位與同步距離

對多振盪器的同步流形：

\mathcal M_{\mathrm{sync}} =\{x:x_i=x_j\ \forall i,j\},

若各振盪器帶有需保留的內部擒縱或介質狀態，完整同步條件須對適當相同型別的內部狀態加上限制。局部橫向擾動 $\eta_\perp$ 的指數若為負，代表附近不同振盪器狀態差有恢復傾向；沿共同週期運動的方向則可能中性。

需要特別區分：

Exact full-state synchronization： 選定完整子系統狀態趨於一致。
Phase locking： 選定相角差趨於固定；振幅或隱藏場未必相同。
Frequency locking： 長期平均頻率一致；瞬時相角差可能振盪。
Orbital synchronization： 距離共同軌道收束，容許時間平移自由度。
Cluster / anti-phase synchronization： 分群內一致，或群與群之間具有不同相位關係。

同一實驗可滿足其中幾種，但不能由其中一種直接推出所有其他種。

7.3 Lyapunov 泛函與物理能量要分開

若存在對某不變集合 $\mathcal M$ 的局部 Lyapunov 泛函 $V\ge0$ ，滿足在相應橫向區域內：

a_1 d_{\mathcal M}(x)^2\le V(x) \le a_2d_{\mathcal M}(x)^2,

且：

\dot V(x)\le-\alpha V(x), \quad \alpha>0,

則可推得局部指數恢復估計。但此 $V$ 可以是偏離軌道的數學量，不必等於機械總能量 $H_m$ 。物理系統沿穩定極限環可以一直有非零功率輸入和耗散，而其橫向誤差泛函仍然下降。

對只具有 $\dot V\le0$ 的模型，不能直接宣稱指數穩定；必要時使用 LaSalle 不變集原理、譜條件或其他附加條件檢查是否存在未衰減模態。

7.4 幾個不同的時間尺度

設驅動、阻尼、耦合、聲學傳播、底座彎曲、熱交換與量測時窗的典型時間尺度分別為 $\tau_d,\tau_\gamma,\tau_K,\tau_a,\tau_s,\tau_T,\tau_{\mathrm{obs}}$ 。若沒有把這些尺度明確量測或無量綱化，所謂「長時間」便沒有共同意義。

例如長於 $\tau_a$ 並不表示長於 $\tau_\gamma$ ；相位可能已快速鎖定，但底座慢模態尚未衰減；也可能聲學回饋先改變了瞬時相位，熱場卻在更長時窗逐漸改變材料阻尼。研究必須允許時間尺度的順序與觀察投影本身成為分析變數，而非由幾秒鐘的影片宣告完成永久穩定。

7.5 暫態放大與非正規演化

對線性擾動 $\dot\eta=A\eta$ ，漸近穩定判定不應只看 $\max\operatorname{Re}\operatorname{spec}(A)<0$ ，還應於指定能量內積下檢查：

G(T)=\sup_{\|\eta_0\|=1}\|e^{AT}\eta_0\|.

若某段時間 $G(T)>1$ ，代表存在初始擾動可在該時間放大，即使所有特徵值實部為負。GMDIS-05 已給出此反例；GMDIS-07 將它納入持續作用的判定：某作用最後會消退，不等於它在過程中從未產生結構性風險。

8. 非光滑擒縱、碰撞與事件：Floquet 不是只看平滑微分方程

8.1 真正節拍器的週期運動可能包含事件

物理節拍器中的擒縱裝置、齒輪接觸、棘爪釋放或撞擊可能在特定狀態面觸發瞬時或短時間的速度改變。因此，單純把連續方程的 Jacobian 在一週期積起來，可能忽略事件造成的一階擾動跳變。

典型的混合模型由連續流 $\dot x=f_\pm(x)$ 、事件面 $h(x)=0$ 、以及重設映射 $x^+=R(x^-)$ 組成。若事件是橫向穿越、事件面與重設映射足夠光滑、且 $n^Tf^-(x^-)\neq0$ ，其中 $n=\nabla h(x^-)$ ，則事件的一階擾動變換由 saltation matrix 給出：

S=DR(x^-) +\frac{[f^+(x^+)-DR(x^-)f^-(x^-)]n^T} {n^Tf^-(x^-)}.

此式採事件面與重設規則不顯含時間的版本。若有顯含時間的觸發面、重設或額外事件延遲，必須使用相應擴充式，不可直接沿用。

8.2 事件存在時的週期單值映射

一個週期內若有多段平滑流與多個事件，單週期擾動傳遞可依發生順序寫成：

M=\Phi_{r}S_{r-1}\Phi_{r-1}\cdots S_1\Phi_1.

其乘子須包含事件的貢獻。相同的節拍器擺角相位與擺幅，只要擒縱事件接觸時機、恢復係數或齒形參數不同，可能具有不同的混合 Floquet 乘子和恢復率。

這正是「相位描述可正確描述同步程度，但單獨不足以決定擾動回復」的工程化案例。不過若觀察者已將事件時間、撞擊條件及修正後的有效相位—振幅座標納入降維模型，則其預測可能重新變得充分。研究要與完善後的基線比較，不應只攻擊沒有事件的極簡模型。

8.3 事件理論的失效點

Saltation 的標準線性化在事件橫向條件失效時，例如掠觸（grazing）或短時間內無限多事件（Zeno 行為），可能失去良好有界性。若模型遇到這些問題，不可藉由誤差很小的數值週期結果宣稱局部穩定性已獲證明；必須調整混合系統的假設、採特定非光滑分析或縮小有效參數域。

9. 穩定模式的競爭、吸引域與多穩態

9.1 D-07-11：模式集合，不是單一同步目標

對完整模型 $R$ ，將指定觀察語境中可區分且有實際意義的穩定模式候選記成：

\mathfrak M_R =\{\mathcal M_1,\mathcal M_2,\ldots\}.

每一模式可配備：週期或特徵頻譜、動力流、橫向恢復、吸引域、可達成本、干預敏感度及實驗可識別性。不同模式並不一定互斥地覆蓋所有初態；某些初態可能走向準週期、混沌、停止、碰撞或尚未分類的運動。

9.2 為什麼模式的出現率依賴初態分布？

令 $\mathcal B(\mathcal M_j)$ 是模式 $\mathcal M_j$ 的吸引域， $\nu_0$ 是研究者明確指定的初態分布。則在確定性系統且吸引域定義適用時：

p_j=\nu_0\bigl(\mathcal B(\mathcal M_j)\bigr)

才有可能被解讀成多次隨機製備初態後最終落入模式 $j$ 的理論頻率。這個 $p_j$ 不是模式本身固有的宇宙機率；改變初態製備方式、外部噪聲或有限觀測時間，都可能改變實際估計值。

9.3 為何不能按「吸引最強」排出普遍單一勝者？

至少存在三個不可混同的排序標準：

局部線性擾動恢復率較快；
對指定初態分布具有較大的吸引域；
受到指定類型干預後較不容易離開模式。

它們沒有必然相同的排序。局部恢復快的模式可能吸引域很小；吸引域大但局部收束慢的模式可能在有限時窗更常被判為「未完成收斂」；受到強脈衝時能恢復的模式，未必對緩慢改變參數同樣有韌性。因此若 GMDIS 研究使用「未來勢更高」這類語句，必須立即補上何種可測操作剖面及排序語境，不能使用抽象單一分數代替一切。

9.4 吸引域不能從有限觀察直接唯一恢復

真實實驗只能從有限的初態集合與有限觀察時間估計模式出現率。有限時間內仍在過渡的軌道不應強迫分配給某個長期模式；未知吸引域邊界、混沌鞍及長壽命暫態必須標記為不確定。較好的研究方法是：事前固定初態抽樣設計，記錄「未分類」比例，對模式判定採盲化條件及誤差上下界，並以不同觀測時間窗複核。

9.5 與受迫、耗散及自主模式的關係

對共同週期外力的雙振子，外部時間週期本身是重要規則的一部分；對自治 Stuart–Landau 系統，模式的週期來自內部振盪。將共同外力當作額外「時鐘狀態」可形成增廣自治系統，但增廣後的同步與原物理對象的自主同步仍有不同因果定義，不應透過數學增廣偷換研究問題。

10. GRST × SIO：以可追溯的未來結構取代模糊的「強勢勝出」

10.1 GRST 的原始定義必須原封不動守住

GRST-12 將 Global Rule-Shi 定義為：規則對可行未來空間所誘導出的結構性非均勻（rule-induced structured non-uniformity over admissible future space）。它不是一種額外的牛頓力、自然界已有勢能的替身，也不是將「流動、持續、穩定」三個詞相加便得到的新物理量。

為使本篇研究可檢驗，在固定的模型 $R$ 、容器 $\mathfrak C$ 、觀察語境 $\Gamma$ 、時間窗 $[0,T]$ 、干預類 $\mathcal U$ 及狀態尺度下，提出一個候選操作剖面：

\mathcal S_{R,\Gamma}(x;T,\mathcal U) =\left\langle \mathcal A_T,\mathcal B_T,\gamma_\perp, G_T,\tau_{\mathrm{rec}},\mathcal E_T, \mathcal I_T,\mathcal Q_T \right\rangle.

符號分別表示有限時窗可達集、各集體模式之吸引域估計、橫向恢復率、擾動增益、恢復時間、能量流剖面、干預敏感度與觀察充分性稽核結果。這只是一組可操作的候選特徵，不是 GRST 原稿已證明的唯一或最終勢表示。

各項型別和單位不相同。例如 $\gamma_\perp$ 是時間倒數、 $\tau_{\mathrm{rec}}$ 是時間、 $\mathcal E_T$ 含焦耳和瓦特， $\mathcal B_T$ 可以是集合或相對指定初態分布的機率。沒有共同尺度、目的函數及單位規約時，不能把它們硬加成一個「規則勢分數」。

10.2 SIO 的三種不同「陰影」與本篇位置

GMDIS-03 已把 SIO 的非局部自由度陰影分成兩個不同的概念：

幾何支撐 $\mathrm{Supp}(O)$ ： 物體實際佔據或約束的空間位置。
硬約束自由陰影 $\mathrm{Shadow}_{\mathrm{hard}}(O)$ ： 即使未直接接觸，某個位置的合法配置或合法自由集合仍被改變。
動態影響陰影 $\mathrm{Shadow}_{\mathrm{dyn}}(O;T,\varepsilon)$ ： 局部干預經由合法動力演化，在指定時窗與量測閾值下可被其他位置的輸出察覺。

第三種是本篇的重要對象。它不是某種「超越物理傳播速度」的即時作用。真實聲波、彈性波和熱擴散的適用方程各有自己的因果結構；連續 PDE 的有效遠端響應也須依模型與閾值判定。資訊纖維仍是局部容許狀態的描述結構，不代表已找到新的物理場粒子。

10.3 D-07-12：模式條件化的有限時窗動態陰影

設 $B_j$ 表示對構成位 $j$ 的合法小干預， $C_i$ 是在位置 $i$ 的觀察算子。對沿模式 $\mathcal M$ 的線性化響應 $\Phi_R(t,s)$ ，定義：

K_{i\leftarrow j}^{R,\mathcal M}(t,s) =C_i\Phi_R(t,s)B_j.

在固定輸入測試空間、量測尺度及容許擾動類的前提下，可以令：

\mathcal G_{i\leftarrow j}(T) =\sup_{\|u\|\le1} \left\|\int_0^T K_{i\leftarrow j}(T,s)u(s)\,ds\right\|.

此時，若 $\mathcal G_{i\leftarrow j}(T)>\varepsilon$ ，就可在該語境中將 $i$ 標記為 $j$ 的可觀察動態影響範圍。若 $\mathcal G_{i\leftarrow j}(T)\le\varepsilon$ ，只代表在本時窗、本干預預算與本感測精度下未檢出顯著影響，不代表物理作用絕對不存在。

與單純耦合圖的邊不同，這個指標會納入多步傳播、共振、延遲、放大及阻尼，並清楚區分「拓撲上存在傳播路徑」與「在指定時窗有足夠大觀察響應」。

10.4 D-07-13：模式穩定性的反事實規則比較

固定同一初始實體配置、感測器和任務，考慮兩組合法規則 $R_A,R_B$ ，例如固定底座和可平移底座、自由旋轉和鎖死旋轉，或兩種材料阻尼。若：

\mathcal S_{R_A,\Gamma}(x;T,\mathcal U) \neq \mathcal S_{R_B,\Gamma}(x;T,\mathcal U),

即在指定的至少一個可測分量上有所差異，則稱該比較對研究任務顯示態勢換位的候選操作證據。但不能由上式倒推所有差異均出自 GRST 新效應，因為標準牛頓／流固耦合模型本來就可能完整解釋這些差異。

真正的研究增益須來自：在相同訓練資料、參數量、計算預算及干預分布下，GMDIS 型特徵是否改善未見規則的模式預測、恢復率估計或不確定性校準；若沒有，原有物理量已構成更簡潔的答案。

10.5 三種時間持續性不能混為一談

在 GRST 對物理世界的應用中，「持續」至少可以指：

作用源持續： 外界電池、發條或自維持回饋可不斷提供功率。
動力關係持續： 特定運動模態對小擾動具長期恢復性。
結構剖面持續： 在指定條件下，某些可達、吸引、成本或存續關係在時間窗間可重現。

三者可以同時成立，也可以彼此分離。持續供能的裝置可能永遠不同步；沒有持續供能的被動裝置可以靜止在穩定平衡點；隨參數緩慢變化，原有的模式結構也可能越過分岔而改變。

10.6 從張力語言走向可測物理，不能靠「方向一致」一個條件

早期《OmniTension》強調多維張力與相互作用的語法，具有跨領域問題生成的價值。本篇對其採取受控詮釋：

在物理層，接觸力、牽引、彈性應力、聲壓、機械功率與熱通量各用合法量綱和守恆式表示。
在動力學層，使用因果傳遞、耦合拓撲、參數依賴和穩定性標記不同作用的長期效果。
在規則勢層，以 GRST 剖面比較不同規則對可行未來的非均勻賦形。

這個三層結構不是三份能量，也不是以哲學詞彙代替材料實驗。它的意義是防止跨尺度研究時，先把不可比較的量硬併為單一尺度再宣稱找到物理定律。

11. 計算 24／72 的正確接入：事件分類，而非宇宙物理枚舉

11.1 原始四軸不改動

依既有正式版，24 格的三軸為：

\mathfrak P_{24}=\mathfrak B_2\times\mathfrak U_4\times\mathfrak O_3,

72 格另增轉移律：

\mathfrak P_{72}=\mathfrak P_{24}\times\mathfrak L_3.

其中 B 是連續／離散基底；U 是序列、跳躍／選擇、並行、識別／檢索；O 是連續、離散或指定語境下拒單測的觀察；L 是確定函數、古典隨機核或真正需要量子通道語義的轉移律。

分類時必須先選定「局部模型事件」及其計算／觀察語境。 物理系統的連續時間自主演化不會因為工程師選擇 CPU 序列迴圈就改變自然律；選擇 GPU 平行網格更新，也不表示宇宙在物理上有一個全域同步數位時鐘。

11.2 GMDIS-07 的候選範式路徑

局部研究事件	範式候選	何時成立及限制
連續受迫機械振盪器的數值演化	C-S-C-F	逐步連續值數值更新；F 是確定轉移
連續場的平行數值求解	C-P-C-F	P 描述數值更新組織，不直接描述物理本體
由連續擺角偵測擒縱事件	C-J-D-F	事件觸發／選擇更新並輸出離散事件
同一場中機率噪聲的數值取樣	C-P-C-K	K 是古典機率核，非機械耦合係數
多模態穩定性量測與分類	C-P-X-F	只有在單一容許表示不能保留所需多不變量時標 X
從過往測試庫辨識已知集體模式	D-R-D-F	R 是算法辨識／存取，非物體自己有檢索器

此表是模型與資料處理事件的語義映射，並非已證明所有節拍器必落在哪一格。對真實自適應計算工作流應保留多格路徑或超圖；不能硬將一整個連續場加離散事件、熱隨機及識別演算法的系統標為單一格。

11.3 拒單測 X 的正確使用時機

若研究者只想判斷「是否頻率鎖定」，並且已證明特定單一特徵在指定語境下足以回答該問題，就沒有理由硬標 X。只有當研究任務要求同時保持相位、振幅、恢復率、能量通量、擾動增益、分岔和介質依賴等重要不變量，且給定容許表示類裡的單一讀出無法容差內保留它們時，X 才有操作上的用途。

因此 X 不是「我認為世界很複雜」的標籤，而是對已指定任務、容差、表示類所做的資訊充分性判斷。GMDIS-02 的受控纖維閉合判準及最壞情形預測下界仍是本篇的正式強制依賴。

11.4 分類標記的獨立效益需另行測試

如果 24／72 無法降低跨模型建構錯誤、提高模型路由或減少反事實稽核時的資訊遺失，就不應將「它能成功標註每個局部事件」直接當成物理研究的進展。形式分類閉合性和實際科學預測效用不同；兩者必須分開評分。

12. 反事實與介質增模：如何真正找出模式形成的原因？

12.1 D-07-14：四類干預再次明示

依 GRST-07 的雙層反事實結構，至少區分：

\mathcal I=\{I_X,I_R,I_B,I_O\}.

其中 $I_X$ 改狀態、 $I_R$ 改規則或本構、 $I_B$ 改容器與邊界、 $I_O$ 改觀察。GMDIS-07 不把四者當成可任意替代的同一操作。

以節拍器為例：

狀態干預： 給一台節拍器施加指定角度或速度脈衝、稍微移動底座，保持模型其他規則不變。
規則干預： 改阻尼係數、擒縱供能規律、接觸恢復係數或耦合網路。
邊界干預： 底座鎖死／解鎖、容器開放／封閉、聲學反射面更換、摩擦支承方式改變。
觀察干預： 改感測頻率、加入麥克風、從俯視改側視、加入熱影像或材料形變量測。

若實際實驗只改感測器，就不能將新的觀測模式說成底層物理規則已改變；若換了底座又重新調整初態，則兩種干預必須在分析表中分開標記。

12.2 幾何與底座的三個必需負向基線

基線 G0： 共同底座完全鎖死，個體仍各自持續驅動。若沒有共同外力或其他耦合，理論上不應因為畫面排列就普遍形成彼此同步。

基線 G1： 底座只允許一維剛體平移、各個機械參數完全相同、接觸與質心作用等對應保持不變。單純重新排列顯示格點，在此理想化模型裡不應自動產生幾何排序效應。

基線 G2： 底座允許平面平移與旋轉，或允許柔性彎曲，並將作用點到旋轉中心的力臂與材料模態正確帶入模型。此時不同排列可能確實選擇不同平移／旋轉模態。2026 年 Hoogesteger 等人的平面底座研究已提供可對照的實驗與模型結果；GMDIS 必須將其視為強基線，而不是自己首次提出的觀測。

三個基線將幫助我們把「資訊自由度陰影」轉成對真實作用路徑的模型檢查，而不是隨意把不同排列所造成的所有現象都稱為新幾何效應。

12.3 聲音：從零聲學到回饋再到真實環境

A0： 關閉聲學回饋，保持底座耦合。A1： 加入線性被動聲場，確保界面功率正確相消。A2： 加入聲學損耗、反射邊界、幾何路徑及必要熱黏性近似。A3： 若線性聲學不再充分，改用流體本構與移動邊界求解並重新核對能量；不可將聲學與完整流體能量重複計算。

各級必須記錄新增的自由度、傳播時間、能量通量及數值成本。若 A1 或 A2 相比 A0 對指定穩定性目標沒有可重複增益，就應允許聲場在該任務、條件與精度下被安全刪去；完整不是永遠算得更多，而是保留任務必需的物理機制。

12.4 驅動機制：公共時鐘與個體擒縱分開

D0： 外部共用週期強迫（用於第 4 節的受迫基線）。D1： 各個體具有獨立但參數相同的自維持驅動。D2： 各個體驅動略有異質、具有供能衰退。D3： 加入非光滑擒縱事件、齒輪接觸和供能儲量。

這組拆分可以回答：同步究竟主要來自外部公共節奏、機械耦合，還是兩者交互作用？若未設 D1、D2、D3，便不宜直接把 D0 的結果當成真實節拍器的完成解釋。

12.5 不以模型複雜度作為先驗優越性

每次增加聲場或更細構成位，必須同時比較：

是否改善對未見干預的預測？
是否降低穩定性與恢復時間的系統偏差？
是否改善能量與界面帳本的真實對應？
是否只增加未識別的參數和過擬合能力？
是否能以更便宜的成熟降維模型取得同等精度？

若增加的場無法提高指定任務充分性，就以簡化模型作為工程首選；這不等於宣稱被省略物理量在宇宙本體上不存在。

13. 從模型到實驗：分層基準、量測與數據結構

13.1 D-07-15：建議的八級實驗階梯

層級	系統規格	必測目標	失效／風險
E0	單一受迫或自維持振盪器	能量、極限環、單機橫向恢復	誤把共同強迫當耦合
E1	兩個同參數個體／無耦合	無耦合相位差持續性	數值同步假象
E2	雙振子，共用移動底座	同向／反相模態、機械帳本	材料摩擦未量測
E3	雙振子，加入聲學／熱域	聲學壓力與熱流的增量	聲學雙重計能
E4	10 個異質個體	局部擾動恢復與分群	初態選擇偏差
E5	100 個，行列／方格／環形	幾何—耦合—模式轉換	忽視轉動自由度
E6	1,000 個，分層模擬	收斂時間、計算成本、誤差上界	算法人工同步
E7	至少一種實物裝置配置	校準、反事實重複性、外推	以模擬冒稱實證

E0 至 E6 可以用模擬與標準模型完成；E7 才具有自然科學的實物驗證角色。若尚無 E7，也應公開 E0–E6 的限制與失敗資料，而不能將「模擬通過」寫成已證實某普遍物理定律。

13.2 D-07-16：預先註冊的觀測剖面

最低建議保存：

個體運動： 擺角、角速度、週期、振幅與擒縱事件時間。
共同介質： 底座平移、旋轉、加速度、柔性模態與必要應變。
聲學： 麥克風時序、頻譜、空間位置與反射配置。
能量／功率： 驅動輸入、機械儲能、界面通量、阻尼耗散、熱通量、數值殘差。
穩定性： 同步流形距離、相位差分布、Floquet／局部指數估計、恢復時間與失穩事件。
模式結構： 同向、反相、分群、無法分類、準週期、停止、長壽命暫態。
因果與可識別： 每次脈衝位置、幅度、開始／結束時間、控制組、儀器校正及觀測語境。

不要求每台實物都裝設全部感測器；若受限於預算，可設計多輪可重複的替代測量。但重要輸出應明確區分直接量測、模型估計、數值重建及純假說。

13.3 漸近穩定與有限時窗指標要一起報告

只報告最後一秒是否同步，容易造成多種失真：快但不穩的過渡可能被誤判為完成同步，局部穩定但恢復很慢的模式可能被誤判為不穩定，或出現碰撞後被動重置的機構被誤認為理想相位吸引。因此建議同時報告：

T_{\mathrm{rec}}(\varepsilon),\quad \sup_{t\in[0,T]}\frac{\|\eta_\perp(t)\|}{\|\eta_\perp(0)\|}, \quad \gamma_{\perp},

以及相應的恢復成功率、能量帳本殘差、未分類軌道比例。此處恢復時間門檻 $\varepsilon$ 、範數、擾動幅度與時間窗需在看到結果之前固定。

13.4 24／72 的資料契約

對每一個局部事件 $e_j$ ，記錄：所屬模型、底空間 B、更新 U、觀察 O、轉移律 L、介質與材料型別、狀態尺度、時間步、感測器與數值求解器版本。大型系統採多格組合及有向關係圖表示。這使計算事件的分類可被審查，但是否提升物理預測仍應另外測試。

13.5 預先接受的零結果

以下結果都可能是科學上完全合理的正式結果：聲場對指定尺度的穩定性沒有顯著增益；有限時窗的動態陰影受測量噪聲掩蓋；某幾何重排在剛體平移模型下毫無影響；高階 GRST 特徵被既有流固耦合指標完全吸收；某些參數區間只有停止或混沌沒有同步。研究不得因為實驗結果不符合初始直覺而修改「什麼算成功」的門檻。

14. 強基線與公平競爭：新語義是否真的提高認識品質？

14.1 應至少比較五組成熟方法

Baseline A： 實驗經典相位／相位差模型，含適用條件與校準過的耦合函數，不使用蓄意削弱的理想 Kuramoto 變體作為唯一對照。

Baseline B： 相位—振幅或 phase–isostable 模型，以及適當情況下的延遲座標重建。若模型已能以很少自由度預測擾動恢復，GMDIS 不應把該結果稱為「舊理論不知道」。

Baseline C： 有質量、阻尼、剛度、擒縱與底座自由度的牛頓多體模型；對允許旋轉／彈性基座的裝置使用相容的完整多體模型。

Baseline D： 介質層的完整或可信縮減 FSI／聲學／熱模型，及其 port-Hamiltonian、GENERIC 或可相容的能量—耗散記帳。

Baseline E： Floquet／Master Stability Function、Koopman 或 Mori–Zwanzig／記憶降維等，依任務適用性選用，並將其成本和可用資料量一併記錄。

GMDIS candidate： 在同等可用資料與合理計算預算下，將 SIO 動態影響陰影、GRST 未來剖面及 24／72 路徑型別作為額外特徵或比較文法，研究其獨立增量，不是將完整 PDE 與其他模型不公平地相比。

14.2 資料分割比圖像漂亮更重要

建議把幾何、阻尼、材料、介質與外部邊界組合分成訓練、校準與測試集；至少一組完整幾何類型或邊界規則只在測試時出現。若 GMDIS 只能記住已見配置的穩定性，不應宣稱理解「世界規則誘導的未來結構」。

樣本量與擾動重複次數必須在收集前規劃，對每組輸出不確定性和信賴區間，並記錄感測器失效、模型發散與無法分類。不可事後只挑同步成功的影片或數值種子。

14.3 模型增益指標的建議

依預先固定的任務，可採：模式分類的校準機率或 Brier score、同步建立時間的 MAE、恢復率估計的可信區間覆蓋、外部邊界干預後的動力軌跡預測誤差、跨介質能量殘差、短時擾動增益誤差與執行成本。

對能量、位移與時間等不同型別量，先做無量綱化並保留原始單位結果；若設綜合分數，必須在試驗前固定權重。一組總分較高，不能代替各子任務誤差與失敗模式的公開。

14.4 模型複雜度與資訊增益的 Pareto 比較

更細緻的世界模型通常需要更多網格、材料係數與邊界數據。若其準確度提升非常小、但所需物理參數幾乎無法校準，就可能不值得在指定任務中採用。比較應至少沿「預測誤差—參數辨識成本—計算成本—可解釋性—魯棒性」呈現 Pareto 結構，而非永遠宣布最大模型最好。

15. 研究命題、判敗條件與證據等級

15.1 H-07-01：持續作用的多特徵剖面假說

假說。 在同一受驅動耦合世界、固定觀測預算與相近資料量下，使用持續時間、有效耦合增益、模態方向、能量交換、延遲與阻尼等分量所構造的剖面，可能比單純峰值力、平均輸入功率或原始相位同步指標更好預測未見擾動的模式恢復。

判敗 F-07-01。 若對公平的相位—振幅和完整多體強基線沒有顯著增量，或增量僅由更多參數量帶來，就撤回獨立改進主張。

15.2 H-07-02：多介質必要性的條件性假說

假說。 在特定幾何、聲學反射、底座彈性或熱漂移條件下，省略某個介質會對恢復時間或模式選擇造成超過事前誤差門檻的系統性偏差。

判敗 F-07-02。 若所有受測條件下含介質模型與無該介質模型的預測等價，則只保留其理論上的可能角色，不宣稱它對當前實驗有必要性。

15.3 H-07-03：動態陰影與局部結構的可辨識性假說

假說。 指定時窗的受控響應核能識別單純幾何支撐或直接圖鄰接未能描述的有效影響傳播範圍。

判敗 F-07-03。 若該影響範圍在所有受測場景都被既有 Green 函數、傳遞函數或數值敏感度分析完整取代而無新用途，GMDIS 應明示該指標只是既有工具的標準化語義接口。

15.4 H-07-04：集體模態的可外推規則勢剖面假說

假說。 在規則層與狀態層反事實均事前指定的實驗中，GRST 剖面或其低維近似可協助預測未見幾何與介質變動導致的模式切換，且超過單一相位指標的表現。

判敗 F-07-04。 若成熟的完整多體／Floquet／FSI 模型已達到相同或更好預測，則不能宣稱新的獨立物理效應；此時可以降級為跨模型比較或解釋文法。

15.5 H-07-05：24／72 分類帶來的模型審計效率假說

假說。 使用 24／72 事件標記可以降低多場模型在轉移律、觀察與求解流程上的型別混淆。

判敗 F-07-05。 若專業多物理仿真工具已有同等或更強的模型型別檢查，而本分類無額外降低錯誤率或加速工程協作，則不得宣稱其分類具有新增自然科學發現。

15.6 F-07-06：能源帳本失效

任何聲稱可穩定長時間非零振盪的機械模型，若在既定系統邊界下既沒有可追蹤的供能，也沒有合法的初始儲能與耗散去向，卻同時違反其應守的能量定律，則模型物理判敗，不能靠「規則勢自動注入能量」修復。

15.7 F-07-07：穩定性證明跨域濫用

若從平滑、同參數、局部 Stuart–Landau 模型推出所有真實節拍器、含接觸事件、空氣流動與材料異質性的全域同步定理，屬於定理適用域越界，應撤回外推。不必等實驗失敗，邏輯上已不成立。

15.8 F-07-08：同一結果卻無法區分因果路徑

若兩種機制在現有觀察下對所有測試都不可辨識，不可因研究者偏好 GRST 而宣稱已確認特定因果路徑。應標記 Rashomon 等價機制，或設計新增觀測與干預以求辨識；若不可辨識性在可行資源內無法解除，結果應保持未決。

15.9 證據階梯

等級	可以宣稱的事	不可以宣稱的事
E0 概念定義	已建立詞義與型別	宇宙已證成此本體
E1 解析證明	指定模型與假設下定理成立	所有物理裝置都符合該模型
E2 數值一致性	指定離散模型、演算法與精度通過檢驗	已完成自然實驗
E3 受控實體實驗	指定條件下可重現觀測與因果差分	可推廣到所有尺度
E4 獨立跨域重複	多種裝置、場域與研究者下穩健	不存在更簡潔的競爭理論
E5 理論整合	具有共同、可用、可反駁的形式系統	必然是終極物理學

16. 交付 GMDIS-08 的正式接口

下一篇《跨尺度投影、記憶與自由度傳輸》不應重新從零發明穩定性。它接手的正式對象包括：

一個有完整能量／介面帳本的多場世界狀態 $X$ 。
一個被觀察的低維子狀態 $y=\pi(X)$ 。
在長時間作用下顯現的共同模態、橫向恢復與擾動增益。
指定的物理介質與耦合接口；已知哪種情況可以消去，哪種情況產生記憶或未解析效應。
既有相位、吸引子、Floquet、MSF 及完整多體模型作為強基線。

GMDIS-08 至少應回答：當我們把共同底座柔性模態、聲腔或熱節點從顯式狀態中消去時，哪些變成有效阻尼、哪些變成頻率相依回饋、哪些形成非 Markov 記憶、哪些仍可被相位—振幅模型充分描述？

它並須逐項核對以下命題：

被消去的狀態是被物理耗散，還是只被模型投影隱藏？
消去是否保留被觀察量在未見干預下的閉合預測？
新增記憶核的時間尺度與本篇恢復率、共振與模式切換是否相容？
對相同計算與資料預算，Mori–Zwanzig、Koopman 或傳統縮減是否已提供足夠答案？
24／72 的路徑標記能否在明確語境下降低模型語義錯置？

這些問題形成第八篇的正式上下文，而不是以「更多資訊會更好」這種永遠不可反駁的口號替代科學工作。

17. 結論：穩定的是動力關係，不是萬物停止流動

本篇的兩個解析模型表明：相對差異可以衰減而共同振盪仍然持續；個體可以自維持卻彼此不鎖定；耦合規則可以在相同組件上選擇局部穩定的不同集體模態。 這些都是在限定方程下的標準非線性動力學結果，並非新發現的普遍物理定律。

本篇真正完成的 GMDIS 工作，是把原本直覺中的「持續性強力主導長期秩序」改寫為數個彼此分離、可量測且可失敗的問題：作用源輸入什麼功率？能量跨哪些界面流動？哪些誤差方向受阻尼或恢復？哪些模式具局部吸引性？其吸引域在指定初態分布下有多大？外部幾何、聲場與邊界更換後，模式會不會換位？當低維投影無法辨識這些差異時，可否由成熟的相位—振幅、延遲嵌入、記憶或完整物理模型加以修補？

GRST 可提供「可行未來非均勻結構」的研究語言；SIO 可幫助我們分開容器、自由度、實際佔據與可觀察動態影響；計算 24／72 可以協助標記模型事件的底空間、更新、觀察與轉移律。但這些框架都不替代牛頓方程、介質邊界條件、已證穩定性判據與外部實驗。

最後的科學命題不是「所有持續力終將收斂成同步」，而是「在指定物理規則與觀察語境下，能否完整追溯哪些持續作用、耗散與耦合關係選擇了某種可存續的集體模式，並且以反事實干預檢查這個解釋」？這項研究可以得出同步，也可以得出反相、多穩態、混沌、停止與不收斂；它對所有這些結果都保留可證偽義務。

參考文獻與來源審計

引用紀律。 下列 DOI 與出版資訊已以出版社、PubMed 或可信學術索引核對。此處對部分付費全文僅核對題名、作者、年份、摘要與可公開書目，不聲稱已逐頁核讀所有付費文章。外部研究僅支持其已有的數學、物理及實驗內容；任何文獻都不被視為直接證實 GMDIS、GRST 或 SIO 尚待驗證的新命題。正式投稿前應進一步核對全文、公式定義、版次及頁碼。

A. EveMissLab／GMDIS 內部理論來源

[I-07-01] Neo.K × Aletheia（2026-10-16），《GMDIS-01｜從同步現象到全域多尺度動力學》，v0.1，EveMissLab。系列總綱、觀察語境、實驗責任與理論讓位規則。

[I-07-02] Neo.K × Aletheia（2026-10-16），《GMDIS-02｜相位與吸引子的資訊充分性邊界》，v0.1，EveMissLab。投影纖維一致性、五種充分性與同投影不同未來的反例。

[I-07-03] Neo.K × Aletheia（2026-10-16），《GMDIS-03｜GRST × SIO：底空間、自由度陰影與未來勢賦形》，v0.1，EveMissLab。硬約束／動態陰影分離與同幾何、不同未來勢的機械反例。

[I-07-04] Neo.K × Aletheia（2026-10-16），《GMDIS-04｜萬有張力的物理語義校準》，v0.1，EveMissLab。物理力、功、應力、能量與廣義影響的合法型別。

[I-07-05] Neo.K × Aletheia（2026-10-16），《GMDIS-05｜多場泛函分析底座》，v0.1，EveMissLab。半群、能量適定性、隱藏場及非正規暫態。

[I-07-06] Neo.K × Aletheia（2026-10-16），《GMDIS-06｜多介質耦合動力學》，v0.1，EveMissLab。固體—聲學—熱場、介面功率、離散數值驗證。

[I-07-07] Neo.K（2026-10-05），《全域規則勢論 12：全域規則勢統合理論》，v0.1。原始 GRST 勢定義、自由與勢分離、可反駁邊界。無界智庫公開版。

[I-07-08] Neo.K（2026-10-04），《全域規則勢論 03：自由度不是勢》，v0.1。相同自由空間、不同方向賦形的形式接縫。無界智庫公開版。

[I-07-09] Neo.K（2026-10-05），《全域規則勢論 07：因果勢與雙層反事實》，v0.1。區分狀態干預與規則干預。無界智庫公開版。

[I-07-10] Neo.K × Aletheia（2026-10-02），《空間資訊佔據、構成纖維與容器自由度》，v0.1。幾何支撐、合法自由及非局部自由度陰影。無界智庫公開版。

[I-07-11] Neo.K（2026-07-19），《計算的二十四重範式》正式版 v4.0。三軸分類與拒單測語境。無界智庫公開版。

[I-07-12] Neo.K × Aletheia（2026-09-17），《從二十四重形態到七十二格動力學：最低充分轉移律、混合格與計算路徑》。無界智庫公開版。

[I-07-13] Neo.K（早期內部研究），《OmniTension：基於張力語法的萬象力學統一理論總綱》。本篇借用問題生成語言，不繼承其未證的物理等價公式。

[I-07-14] Neo.K × Aletheia（2026-09-22），《全域耦合計算場論的形式核心》。僅作模型耦合與計算路由概念參照，不等於物理多場定理。

B. 外部非線性動力學、物理與系統研究

[E-07-01] Pantaleone, J. (2002). “Synchronization of metronomes.” American Journal of Physics, 70(10), 992–1000. DOI: 10.1119/1.1501118。共同移動底座上的實際機械同步模型基線。

[E-07-02] Hoogeboom, F. N., Pogromsky, A. Y., & Nijmeijer, H. (2016). “Huygens' inspired multi-pendulum setups: Experiments and stability analysis.” Chaos, 26, 116304. DOI: 10.1063/1.4965031。行列與十字排列的節拍器實驗、同向與反相模式和局部穩定性對照。

[E-07-03] Hoogesteger, M. M., Pogromsky, A. Y., & Nijmeijer, H. (2026). “Huygens' synchrony: Emergent modes of planar symmetrically coupled pendula.” Chaos, 36(9), 093135. DOI: 10.1063/5.0331194。平面平台運動中的平移／旋轉集體模態及幾何轉換的實驗與多體模型強基線；出版資訊另由 PubMed 核對。

[E-07-04] Kada, H., Kondo, K., Sumino, S., & Tokuda, I. T. (2017). “Experimental study on synchronization of sound-coupled electronic metronomes.” Nonlinear Theory and Its Applications, IEICE, 8(1), 15–24. DOI: 10.1587/nolta.8.15。聲音可作為某類電子節拍器的同步耦合媒介；不可直接推廣為所有機械節拍器。

[E-07-05] Pecora, L. M., & Carroll, T. L. (1998). “Master Stability Functions for Synchronized Coupled Systems.” Physical Review Letters, 80, 2109–2112. DOI: 10.1103/PhysRevLett.80.2109。將耦合拓撲與振盪器橫向穩定性分解的成熟方法。

[E-07-06] Nakao, H. (2016). “Phase reduction approach to synchronisation of nonlinear oscillators.” Contemporary Physics, 57(2), 188–214. DOI: 10.1080/00107514.2015.1094987。相位降維與同步分析的適用條件和強基線。

[E-07-07] Shirasaka, S., Kurebayashi, W., & Nakao, H. (2017). “Phase reduction theory for hybrid nonlinear oscillators.” Physical Review E, 95, 012212. DOI: 10.1103/PhysRevE.95.012212。非光滑混合振盪器仍可接受合適相位理論分析的反方參照。

[E-07-08] Segnou Tchio Kuete, W., Muolo, R., Dorchain, M., Nakao, H., & Carletti, T. (2026). “Synchronization of nonlinearly coupled Stuart–Landau oscillators on networks.” Proceedings of the Royal Society A, 482, article 20250907. DOI: 10.1098/rspa.2025.0907。非線性與非對稱網路耦合的橫向穩定、Floquet 及模式分析；本篇的正規型例子為獨立基本模型推導，不是該論文特定非線性結果的複製。

[E-07-09] Willems, J. C. (1972). “Dissipative dynamical systems part I: General theory.” Archive for Rational Mechanics and Analysis, 45, 321–351. DOI: 10.1007/BF00276493。storage function、supply rate、耗散與系統互連的基礎理論。

[E-07-10] 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: 10.1561/2600000002。能量埠及多物理互連的成熟綜述。

[E-07-11] Grmela, M., & Öttinger, H. C. (1997). “Dynamics and thermodynamics of complex fluids. I. Development of a general formalism.” Physical Review E, 56, 6620–6632. DOI: 10.1103/PhysRevE.56.6620。GENERIC 中可逆與不可逆演化的形式分離。

[E-07-12] Chorin, A. J., Hald, O. H., & Kupferman, R. (2000). “Optimal prediction and the Mori–Zwanzig representation of irreversible processes.” PNAS, 97(7), 2968–2973. DOI: 10.1073/pnas.97.7.2968。研究降維之未解析狀態與記憶作用的成熟起點，交接 GMDIS-08。

[E-07-13] Chorin, A. J., Hald, O. H., & Kupferman, R. (2002). “Optimal prediction with memory.” Physica D, 166(3–4), 239–257. DOI: 10.1016/S0167-2789(02)00446-3。記憶與粗粒化預測的實務延伸。

[E-07-14] Trefethen, L. N., & Embree, M. (2005). Spectra and Pseudospectra: The Behavior of Nonnormal Matrices and Operators. Princeton University Press. Publisher / JSTOR。漸近譜穩定與有限時間放大並不等價的權威數學參照。

[E-07-15] Pikovsky, A., Rosenblum, M., & Kurths, J. (2001). Synchronization: A Universal Concept in Nonlinear Sciences. Cambridge University Press。振盪同步、鎖定與集體模式的經典參考書；書目名稱和作者依出版社版次為準。

[E-07-16] van der Schaft, A. J., & Maschke, B. M. (2002). “Hamiltonian formulation of distributed-parameter systems with boundary energy flow.” Journal of Geometry and Physics, 42, 166–194. DOI: 10.1016/S0393-0440(01)00083-3。連續介質與邊界能流的分布參數形式。

C. 跨文獻對照的警戒

本篇不把 2026 年的 Stuart–Landau 網路研究誤寫成所有非線性耦合的完備答案。其方法依模型與拓撲有適用邊界。
本篇不把新的模式名稱或張力術語當作對已知 Floquet／MSF 分析的替代。
本篇採用的聲學與熱模型在前作 GMDIS-06 中已有能量對帳；任何加入新場的式子都應依所採模型重新推導，不能只做「更多物理項」的拼接。
對 GMDIS 真正的新增性，需要由可操作的反事實、模型可比較性、觀察充分性及可重複增量而非跨領域類比來判定。

附錄 A｜符號、單位與型別對照表

符號	定義	型別／單位	警示
$X(t)$	完整世界狀態	有限維／多場函數空間	取決於指定模型
$R$	模型規則與邊界	規則結構	不等於電阻或機械力
$H_m$	機械儲能	J	不等於 GRST 勢
$E_a,E_{th}$	聲學與熱儲能	J	不能同一份能量雙重計數
$P_{\mathrm{drive}}$	外部輸入機械功率	W	若源在容器內應改列內部流
$D_m,D_a$	指定模型之轉熱功率	W，非負	非所有「損失」都一定入同一熱節點
$q_j,\dot q_j$	位移及速度	m, m/s	若採擺角需改為 rad, rad/s
$m,c,k,\kappa$	質量、阻尼、個體剛度、互連剛度	kg, kg/s, N/m, N/m	須維持參數非負條件
$F,\Omega$	外部週期力、角頻率	N, s $^{-1}$	共用外力不是互相同步
$a,e$	平均運動／相對差異	m	$e\to0$ 不代表 $a\to0$
$\lambda,\omega,K$	Stuart–Landau 正規型參數	依無量綱時間定義	$K$ 非物理牛頓力，也不同於機械 $k$
$L,\mu_j$	圖拉普拉斯與特徵值	隨圖權重尺度	須明示圖連通、對稱條件
$\gamma_\perp$	橫向擾動指數	時間倒數	局部判據，非吸引域大小
$\Phi(t,s)$	線性化傳遞算子	映射／算子	必須有合法演化型別
$S$	事件 saltation matrix	線性映射	掠觸條件不能直接套用
$\Sigma_R$	GRST 全域規則勢	未來結構語義	非 J、N 或 Pa 的普遍單位
$\mathfrak F$	SIO 自由結構	容許狀態集合／纖維	非假定的物理纖維粒子
$G_{j,\Gamma}$	指定輸入／輸出增益	由尺度規約決定	不可未無量綱化直接比較
$\mathcal M_R$	不變／集體模式集合	狀態空間子集合	同時記錄穩定性與適用域
$\nu_0$	初態製備分布	機率測度	不等於宇宙固有成功率

附錄 B｜兩個解析模型的最小覆核

B.1 受迫雙振子分解

對方程 $m\ddot q_j+c\dot q_j+kq_j+\kappa(q_j-q_{3-j})=F\cos\Omega t$ ，令 $a=(q_1+q_2)/2$ 、 $e=q_1-q_2$ 。相加與相減可得：

m\ddot a+c\dot a+ka=F\cos\Omega t,

m\ddot e+c\dot e+(k+2\kappa)e=0.

相對模態的根由 $mr^2+cr+k+2\kappa=0$ 決定。 $m,c,k>0$ 、 $\kappa\ge0$ 保證根實部為負。平均模態受迫穩態之複幅滿足：

(k-m\Omega^2+i c\Omega)\widehat a=F.

總機械能在第 4.3 節定義，其時間導數中互連彈簧的雙向功率相消；剩下的是共同外力做功與兩個黏性阻尼耗散。若 $\kappa=0$ ，相對模態仍因個別阻尼衰減，所以此構造不提供自主耦合同步的專屬證據。

B.2 正耦合正規型線性化

旋轉座標為 $w_j=e^{-i\omega t}z_j$ 。令同向解 $w_j^\ast=\sqrt\lambda$ ，擾動 $\delta w_j=u_j+i v_j$ 。在實數微分下：

\frac{\partial}{\partial u} [(\lambda-|w|^2)w]\bigg|_{w=\sqrt\lambda} =-2\lambda,

相角方向的一階係數為零。因此：

\dot u_j=-2\lambda u_j+K(\bar u-u_j),

\dot v_j=K(\bar v-v_j).

寫 $u_j=\bar u+u_j^\perp$ 、 $\sum u_j^\perp=0$ ， $v$ 同理。均值徑向速率 $-2\lambda$ ，均值角向為零；各橫向徑向 $-2\lambda-K$ ，橫向角向 $-K$ 。對 $K>0$ ，相位差具有局部線性恢復能力；不是所有初態都已證明會進入該模式。

B.3 反相正規型線性化

對 $N=2$ 、 $K<0$ ，反相解在旋轉座標為 $(r,-r)$ 、 $r^2=\lambda-K$ 。在 $r$ 或 $-r$ 的線性化中，局部單體向量場的實向徑向係數為 $\lambda-3r^2$ ，虛向切向係數為 $\lambda-r^2$ 。對稱擾動不受平均場差分項影響，故：

\gamma_{\mathrm{sym},r}=\lambda-3r^2=-2\lambda+3K,

\gamma_{\mathrm{sym},\phi}=\lambda-r^2=K.

反對稱擾動額外獲得耦合係數 $-K$ ：

\gamma_{\mathrm{anti},r}=\lambda-3r^2-K=-2(\lambda-K),

\gamma_{\mathrm{anti},\phi}=\lambda-r^2-K=0.

後者為週期軌道的共同相位方向；其餘在 $K<0$ 、 $\lambda>0$ 下均為負。此證明只適用該正規型的局部週期解，不含大擾動、參數異質或聲學延遲。

B.4 圖拉普拉斯版本

以 $-KLz$ 替代平均場 $K(\bar z-z)$ ，對無向非負圖有 $L=L^\top\succeq0$ 。在同向振盪解的旋轉座標內，徑向與角向擾動分別滿足：

\dot u=(-2\lambda I-KL)u,

\dot v=-KLv.

在拉普拉斯本徵模態 $v_j$ 上的指數為 $-2\lambda-K\mu_j$ 和 $-K\mu_j$ 。對連通圖的 $\mu_2>0$ 、 $K>0$ ，除共同相位中性方向外所有橫向指數為負。若網路不連通，則 $\mu=0$ 的重數大於一，每個不連通分量可以保有自己的獨立相位；此時不可宣稱整個網路全域相位鎖定。

附錄 C｜可重現的 NumPy 數值一致性驗證

用途與界線： 下列程式可獨立重算第 4–5 節的兩種模型。它驗證數學推導與特定數值方法之間的一致性，不是新物理裝置實驗，也不代表上游多介質 PDE 模型已被實物量測。程式只需要 Python 3 與 NumPy，不存取網路，不讀取私人檔案，不產生不可追溯的隨機樣本。執行指令：

python -m pip install numpy
python gmdis07_validation.py

以下是與本文數值核對一致的完整最小程式：

"""GMDIS-07: deterministic analytical/numerical consistency checks.
Models: (A) a damped, commonly forced mechanical pair;
        (B) mean-field Stuart–Landau oscillators.
No network I/O. Requires numpy. Not physical experimental validation.
"""

import numpy as np
from math import cos


def mechanical_pair():
    # Dimensional interpretation if q has units m and time has units s:
    # m=kg, c=kg/s, k=kappa=N/m, F=N, Omega=1/s.
    m = 1.0
    c = 0.2
    k = 1.0
    kappa = 0.3
    F = 0.2
    Omega = 1.1
    dt = 0.01
    steps = 30000

    L2 = np.array([[1., -1.], [-1., 1.]])
    A = np.block([
        [np.zeros((2, 2)), np.eye(2)],
        [-(k * np.eye(2) + kappa * L2) / m, -c * np.eye(2) / m]
    ])
    I = np.eye(4)
    left_inv = np.linalg.inv(I - dt * A / 2.0)
    Q = left_inv @ (I + dt * A / 2.0)

    def H(x):
        q, v = x[:2], x[2:]
        return (0.5 * m * (v @ v)
                + 0.5 * k * (q @ q)
                + 0.5 * kappa * (q[0] - q[1]) ** 2)

    x = np.array([0.7, -0.4, 0.0, 0.2], dtype=float)
    max_ledger_error = 0.0
    accumulated_error = 0.0
    samples = {}
    for n in range(steps):
        t_mid = (n + 0.5) * dt
        f_mid = F * cos(Omega * t_mid)
        b_mid = np.array([0., 0., f_mid / m, f_mid / m])
        x_new = Q @ x + dt * left_inv @ b_mid
        v_mid = (x[2:] + x_new[2:]) / 2.0
        expected_delta_H = dt * (
            f_mid * np.sum(v_mid) - c * np.dot(v_mid, v_mid)
        )
        error = H(x_new) - H(x) - expected_delta_H
        max_ledger_error = max(max_ledger_error, abs(error))
        accumulated_error += error
        x = x_new
        if n + 1 in (1000, 5000, 10000, 30000):
            samples[(n + 1) * dt] = float(abs(x[0] - x[1]))

    # Analytic steady-state amplitude of the common coordinate a.
    amplitude = abs(F) / np.sqrt(
        (k - m * Omega ** 2) ** 2 + (c * Omega) ** 2
    )
    assert max_ledger_error < 1e-10, max_ledger_error
    assert samples[300.] < 1e-9, samples[300.]
    assert amplitude > 0
    print("Mechanical pair")
    print("  Analytic common-mode amplitude:", float(amplitude))
    print("  Relative-distance samples:", samples)
    print("  Max per-step energy-ledger error:", float(max_ledger_error))
    print("  Accumulated ledger roundoff:", float(accumulated_error))


def stuart_landau():
    lam = 1.0
    omega = 2.0
    dt = 0.01
    steps = 12000

    def rhs(z, K):
        # z is a NumPy complex array; average is mean over oscillators.
        return ((lam + 1j * omega - np.abs(z) ** 2) * z
                + K * (np.mean(z) - z))

    def rk4(z, K):
        a = rhs(z, K)
        b = rhs(z + 0.5 * dt * a, K)
        c = rhs(z + 0.5 * dt * b, K)
        d = rhs(z + dt * c, K)
        return z + dt * (a + 2*b + 2*c + d) / 6.0

    def run(K, initial, mode):
        z = initial.copy()
        samples = {}
        for n in range(steps):
            z = rk4(z, K)
            if n + 1 in (500, 2000, 5000, 12000):
                if mode == "sync":
                    distance = abs((z[0] - z[1]) / 2.0)
                else:
                    distance = abs((z[0] + z[1]) / 2.0)
                samples[(n + 1) * dt] = float(distance)
        return samples, z

    initial = np.array([
        (1.0 + 0.015) * np.exp(0.03j),
        (1.0 - 0.015) * np.exp(-0.03j)
    ], dtype=np.complex128)
    positive, z_pos = run(0.4, initial, "sync")
    zero, z_zero = run(0.0, initial, "sync")

    amp_anti = np.sqrt(lam + 0.4)
    initial_anti = np.array([
        amp_anti * (1.0 + 0.015) * np.exp(0.03j),
        -amp_anti * (1.0 - 0.015) * np.exp(-0.03j)
    ], dtype=np.complex128)
    negative, z_neg = run(-0.4, initial_anti, "anti")

    # Orbital directions: in-phase K>0 and anti-phase N=2 K<0.
    positive_exponents = [-2*lam, 0.0, -(2*lam + 0.4), -0.4]
    negative_exponents = [0.0, -2*(lam + 0.4), -0.4,
                          -2*lam + 3*(-0.4)]
    assert positive[50.] < 1e-8
    assert zero[50.] > 0.02
    assert negative[50.] < 1e-8
    assert abs(abs(z_neg[0]) - amp_anti) < 1e-6
    print("Stuart-Landau model")
    print("  K=+0.4, in-phase distance:", positive)
    print("  K= 0.0, phase disagreement persists:", zero)
    print("  K=-0.4, anti-phase distance:", negative)
    print("  Analytic in-phase exponents:", positive_exponents)
    print("  Analytic anti-phase exponents:", negative_exponents)
    print("  Final anti-phase amplitude:", float(abs(z_neg[0])))
    print("  All checks passed")


if __name__ == "__main__":
    mechanical_pair()
    stuart_landau()

C.1 已重算的參數與結果

受迫雙振子採 $m=1$ 、 $c=0.2$ 、 $k=1$ 、 $\kappa=0.3$ 、 $F=0.2$ 、 $\Omega=1.1$ 、 $\Delta t=0.01$ ，初態 $(q_1,q_2,\dot q_1,\dot q_2)=(0.7,-0.4,0,0.2)$ 。隱式中點法運行 30,000 步後，程式的局部能量帳本每步最大殘差約 $4.0\times10^{-16}$ ，累積浮點殘差約 $-8.4\times10^{-13}$ 。解析共同運動振幅約 $0.65759595$ ，相對位移絕對值如下：

| 時刻 $t$ | $|q_1-q_2|$ （同初態與程式） | |---:|---:| | 10 | 約 0.4031812302 | | 50 | 約 0.0071400567 | | 100 | 約 0.0000440637 | | 300 | 約 $4.5\times10^{-14}$ |

此收斂說明的是指定阻尼、線性耦合與共同週期力下的數值結果。結果沒有證明真實發條節拍器在同樣位置或聲學條件下也會達到相同恢復率。

Stuart–Landau 模型採 $\lambda=1$ 、 $\omega=2$ 、 $\Delta t=0.01$ ，各組 12,000 步，初態非常接近相應同向或反相週期解。若 $K=+0.4$ ，同向差異距離在 $t=5,20,50$ 約為 $0.0040611,1.01\times10^{-5},6.19\times10^{-11}$ ；若 $K=0$ ，差異約維持 $0.0299955$ ；若 $K=-0.4$ ，反相偏離距離在 $t=5,20,50$ 約為 $0.0048052,1.19\times10^{-5},7.32\times10^{-11}$ ，反相振幅趨近 $\sqrt{1.4}\approx1.183216$ 。

浮點算術與離散時間格式使最終顯示零誤差不等於精確數學零；本篇已提供獨立解析特徵指數及程式可執行重算，因此不依靠畫面感受作為證明。

C.2 誤差區分

應至少分為：

解析模型誤差： 正規型或線性彈簧本身與真實裝置的落差。
數值離散誤差： 時間步與演算法對連續模型的近似。
能量帳本殘差： 在本離散方案下，模型指定的能量恆等式是否被保存。
參數不確定性： 質量、摩擦、聲學邊界與擒縱供能難以準確量測。
模式判定誤差： 有限時間、採樣頻率與閾值可能混淆同步、反相與長暫態。

第三項接近機器精度，不代表前四項也已接近零。

附錄 D｜最小實驗登錄樣板

以下是建議的 YAML 格式，可納入後續 GMDIS-12／13 的模型註冊與資料驗證。所列內容是樣板，不是已執行的實物紀錄。

study_id: "GMDIS-07-E2-SAMPLE"
document_version: "v0.1"
evidence_level: "simulation_proposal"
physical_claim_status: "not_experimentally_verified"
world:
  device: "metronome_array"
  oscillator_count: 2
  model_class: "autonomous_escapement_with_movable_base"
  physical_dimensions: ["q", "qdot", "base_translation", "drive_state"]
  geometry:
    coordinates_m: [[-0.08, 0.0], [0.08, 0.0]]
    base_translation: true
    base_rotation: false
    base_flexure: false
  environment:
    acoustic_feedback: false
    material_thermal_model: false
    background_flow: false
  clocks:
    integration_step_s: 0.0001
    observation_interval_s: 0.01
    total_duration_s: 120.0
    convergence_warmup_s: 20.0
  rule_set_hash: "TO_BE_FILLED"
  initial_state_hash: "TO_BE_FILLED"
  solver_version: "TO_BE_FILLED"
  random_seed: null
  observation_contract:
    outputs: ["angle", "angular_velocity", "base_motion", "input_power"]
    units: ["rad", "rad/s", "m", "W"]
    mode_labels: ["in_phase", "anti_phase", "cluster", "unsynchronized", "unclassified"]
    phase_tolerance_rad: 0.05
    recovery_distance_tolerance: 0.05
    tolerance_pre_registered: true
  intervention_contract:
    type: "state"
    target_oscillator: 0
    angular_impulse_rad_s: 0.1
    injection_time_s: 60.0
    controls: ["no_intervention", "fixed_base", "no_coupling"]
  reported_outputs:
    energy_residual: null
    transverse_recovery_rate_per_s: null
    mode_and_confidence: null
    unclassified_fraction: null
    runtime_seconds: null
  provenance:
    raw_artifact_uri: null
    script_hash: null
    measurement_calibration_ref: null
    failure_log_uri: null

樣板需依實驗級別調整。 對真正的實物研究，應填入儀器型號、感測器校準、實測環境、操作者、原始時間戳、驅動供能測量與實驗批次；對數值研究則應填入求解器容差、空間網格、時間積分收斂測試、電腦精度與模型版本。若初始條件、規則與求解器已有版本化快照，需使用真正的雜湊值而不是上面佔位字串。

附錄 E｜審稿與後續交付自我檢查

物理語義

文章中使用的力、功率、能量、應力、阻尼及熱流均有合法型別與單位。
未將「總力量守恆」稱作標準物理守恆定律。
未將 GRST 全域勢或 SIO 資訊纖維等同新增物理粒子或焦耳量。
對聲學、流固與熱場沒有雙重計能或未追蹤的能量來源。

數學條件

共同受迫雙振子與自治 Stuart–Landau 正規型已分離。
正耦合、零耦合與負耦合的穩定性結論只在指定模型域成立。
局部軌道穩定與所有初態的全域收斂已明確區分。
對週期軌道保留一個共同相位中性方向。
對混合擒縱模型沒有省略事件與 saltation 的必要條件。
長期穩定指標與有限時間擾動增益分別報告。

實驗判敗

以事前固定的時間窗、初態分布、感測精度與模式標籤做比較。
有無共同外力、有無耦合的控制組能隔離不同同步機制。
剛體平移的幾何零效應基線已保留。
對聲學、底座旋轉、材料阻尼設消融測試。
與相位—振幅、Floquet、MSF、完整多體／FSI 等成熟強基線公平競爭。
無法分類的軌道、失敗資料和負面結果會被保存。

檔案與版本

UTF-8 編碼、換行、Unicode 與 LaTeX delimiter 配對已驗證。
完整原稿只使用 $...$ 或 $$...$$ 的數學分隔符。
公式內複雜命令保留為正式 Markdown 原始碼，不以聊天渲染器的缺陷破壞數學語義。
對外 DOI、內部前作、數值程式與符號表都有來源或可重現路徑。
本篇暫未做出的真實裝置實驗沒有被標為已完成。

文件狀態聲明

GMDIS-07 v0.1 是可提交後續形式檢查、數值重現及外部學術審查的研究原稿，而非已受同行審查的物理定理合集。 文中模型解析結果及數值一致性測試已明確標示條件；任何跨模型普遍性、GRST／SIO 額外預測能力，以及千台真實節拍器的大規模同步行為仍屬待驗證問題。後續 GMDIS-08 至 GMDIS-14 只有在繼續遵守這些邊界時，才能累積可比較的科學證據。

GMDIS-07｜持續作用與動態穩定結構：驅動、耗散、橫向穩定與集體模態選擇