2025 IYPT 第九題「磁力輔助」：
磁偶極交互作用與
非線性混沌動力學 之深度研究報告

摘要 (Abstract)

本專題探討 IYPT 磁力輔助系統之決定性混沌動力學。在研究過程中，我結合高中電磁學觀念建立哈密頓模型，嚴謹推導出 \(1/r^4\) 磁偶極解析解，並親自撰寫 RK4 演算法精準運算相空間的碎形吸引域。此外，我更透過實體驗證與硬體除錯，成功消除冷次定律帶來的電磁阻尼干擾。這段從理論推導到實體除錯的歷程，充分展現我將複雜物理數學轉化為程式碼的實作能力，也確立了我對電機與資訊工程系統分析的深厚熱忱。

一、 系統物理圖像與研究動機

2025年國際青年物理學學家錦標賽（IYPT）第九題「磁力輔助（Magnetic Assist）」提出了一個極具深度的經典力學與電磁學耦合問題。題目要求將一到兩個磁鐵固定在一個非磁性且不導電的底座上，使其能吸引一個懸掛在繩子上的移動磁鐵，並探究該移動磁鐵的運動如何受到相關參數的影響。

在研讀文獻時，我發現此實驗裝置常被稱為「混沌磁性單擺（Chaotic Magnetic Pendulum）」。這個系統完美展現了決定性混沌（Deterministic Chaos）的核心特徵，亦即對初始條件的極度敏感性（Sensitivity to initial conditions），並深刻詮釋了非線性動力學中多重吸引子（Multiple attractors）的競爭機制。為了徹底解開這個複雜系統，我決定親自動手，透過數學建模、演算法模擬與實體除錯來進行全面剖析。

二、 磁偶極交互作用之嚴格解析與數學建模

要精確預測移動磁鐵的軌跡，我的首要任務是建立準確的磁力模型，因為磁力的強烈空間非線性正是引發混沌的根源。

1. 解析力方程式的推導

在偶極近似（Dipole approximation）的前提下，我透過對磁場的空間導數進行泰勒展開，嚴謹推導出一個磁偶極 \(m_a\) 在另一個磁偶極 \(m_b\) 產生的磁場中所受到的平移力 \(\mathbf{F}_{ab}\) 解析方程式如下：

從這條方程式中，我發現了一個極度關鍵的特徵：磁力大小與距離的四次方 (\(1/r^4\)) 成反比。這與大家熟知的萬有引力平方反比定律有著本質上的區別。當擺錘接近底部磁鐵時，磁力會呈爆炸性增長，產生極其劇烈的拉扯效應（Tug-of-war effect），正是這股力量將單擺原本規律的軌跡撕裂並引入混沌狀態。同時，兩磁偶極間的磁力矩（Torque）與距離的三次方（\(1/r^3\)）成反比，這也為系統引入了額外的角位移自由度與能量轉換。

2. 連立運動方程式與位能曲面建構

接著，我將整個系統視為在二維水平面（\(x\)-\(y\) 平面）上運動的阻尼諧振子。設擺錘質量為 \(m\)，擺長為 \(L\)，小角度近似下重力線性回復力常數為 \(k = mg/L\)。系統總位能 \(V(\mathbf{x})\) 是由重力拋物面位能與多個磁勢能疊加而成的純量場：

其中 \(p_n\) 決定極性（\(+1\) 為相吸，\(-1\) 為相斥）。引入等效線性黏滯阻尼係數 \(b\) 後，我列出了移動磁鐵的牛頓第二運動定律微分方程：

方程式右側的梯度算子 \(\nabla V(\mathbf{x})\) 結合了純淨的二維重力與強烈局部扭曲的磁力；而阻尼項 \(-b\dot{\mathbf{x}}\) 則確保系統動能被平滑抽走，使軌跡在相空間體積收縮，最終必然趨於某個穩定平衡點（吸引子）。

三、 動力學分岔與拓樸特徵：單磁鐵 vs. 多磁鐵

在探究過程中，我發現底座磁鐵數量的改變，將導致系統位能曲面發生根本性的拓樸相變（Topological phase transition）。

1. 單一磁鐵：週期性與高度可預測性

當我只配置單一磁鐵且相吸（\(N=1\)）時，系統呈現單一的「位能深谷（Potential deep well）」。此機制等同於引入「磁彈簧（Magnetic anti-spring）」改變振盪頻率。相空間中僅存在一個穩定的全域吸引子，系統呈現高度可預測的週期性運動，完全不具備混沌特徵。

2. 雙磁鐵系統：對稱性破缺與鞍點動力學

然而，當我引入第二個磁鐵（\(N=2\)）時，位能深谷被撕裂為「雙井勢（Double-well potential）」，中間形成拓樸學上的「鞍點（Saddle point）」。若擺錘能量大於鞍點勢壘，便能在兩磁鐵間自由漫遊；當軌跡經過鞍點附近時，微小的橫向速度偏差會沿著不穩定流形被指數級放大，這就產生了極度難以預測的混沌行為。直到動能耗散低於閾值，擺錘才會被困入單一位能井中。

四、 程式數值模擬與碎形吸引域 (Fractal Basins of Attraction)

混沌系統最迷人之處在於其相空間的幾何結構。為克服純解析解的限制，我親自動手撰寫了數值積分演算法，來重建並視覺化系統的演化過程。

1. RK4 數值積分與超級採樣

我將二階微分方程降階為四個一階微分方程，並利用四階龍格-庫塔法 (RK4) 進行步進積分求解。為了精準繪製引導系統停留在特定磁鐵的「吸引域」，我的程式運算了成千上萬個微小差異的初始釋放點，並運用了超級採樣（Supersampling）技術來處理邊界抗鋸齒，最終生成了極高解析度的相空間圖像。

2. 碎形邊界與 Wada 性質

模擬結果令我十分驚豔！吸引域的邊界並非平滑的曲線，而是呈現出無限精細且自我相似的「碎形（Fractal）」結構。在三個以上磁鐵的系統中，更展現了 Wada 性質 (Basins of Wada)——任何兩個顏色區域的邊界上必然與第三種顏色相鄰。這完美印證了混沌系統「對初始條件極度敏感」的物理本質，在邊界地帶，微觀的初始誤差必然導致完全隨機的巨觀結果。

五、 控制參數之多維度矩陣分析

為全面探究系統行為，我系統性地分析了相關參數對整體動力學的影響：

參數類別	符號與定義	對系統動力學與混沌特徵之具體影響分析
磁鐵配置與極性	極性參數 (\(p_n\))	相吸 (\(p_n=+1\)) 形成極深的位能井，強烈吸引擺錘並使其快速靜止。 相斥 (\(p_n=-1\)) 形成位能丘，擺錘會迴避磁鐵並沿著磁場最弱的路徑滑行，大幅延長瞬態漫遊時間。
磁場強度與幾何	磁偶極矩 (\(\mu\))	決定位能曲面變形的劇烈程度。極強的 \(\mu\) 導致重力被完全壓制，活動範圍急遽縮小；較弱的 \(\mu\) 允許更廣泛的空間探索。
空間垂直間距	高度分離 (\(H, h\))	最關鍵的分岔參數 (Bifurcation parameter)。因力隨 \(1/r^4\) 衰減，微調高度將導致系統在「重力主導規律擺動」與「磁力主導深井捕捉」間發生劇烈相變。
空間水平間距	磁鐵間距 (\(D\))	決定雙井或多井勢能的重疊程度與鞍點高度。距離適中時，不同磁鐵吸引力產生激烈「拔河」，雙穩態與混沌混合現象最顯著。
單擺物理屬性	繩長 (\(L\))	改變自然頻率。繩長越長，擺錘在磁鐵上方掠過速度越慢，滯留時間變長，使局部非線性磁力有更充分時間發揮扭曲軌跡效應。
擺錘物理屬性	總質量 (\(m\))	決定系統慣性。質量大能儲存更多動能，更容易衝破磁力鞍點的束縛勢壘，增加在多個吸引子間跳躍的存活時間。
環境能量耗散	阻尼係數 (\(b\))	決定瞬態混沌的壽命長度與相空間收縮速率。較低阻尼允許系統繪製出更繁複的李薩如圖形與長程混沌軌跡。

六、 實驗裝置設計與干擾抑制 (硬體 Debugging)

在理論推導近乎完美的情況下，實際動手架設實驗卻讓我面臨不少挑戰，這個除錯過程深刻體現了工程實務的價值。

1. 底座材質與冷次定律 (Lenz's Law) 的逆襲

起初為了求裝置穩固，我使用了導電的金屬底座。實驗時卻驚訝地發現，當帶有強磁場的釹磁鐵高速掠過底座時，會在導體內部激發渦電流 (Eddy currents)。根據高中物理學到的冷次定律 (Lenz's Law)，這產生了極強大的非線性反向磁阻尼力，徹底破壞了我原本設定的線性空氣阻尼常數 \(b\)，導致擺錘軌跡異常快速衰減，根本跑不出混沌漫遊。發現這個致命問題後，我嚴格將底座全面替換為壓克力等非磁性且不導電材質，終於成功消除了這個隱藏變因。

2. 軌跡觀測與影像追蹤

由於混沌系統無法以肉眼精確量測，我突發奇想在擺錘中心安裝了微型光源，並在正上方架設高速攝影機。結合 Tracker 開源影像追蹤軟體，我成功逐格擷取擺錘亮點的二維時間序列座標 \((x(t), y(t))\)，並將其繪製為相圖 (Phase Portrait)，完美地驗證了數值模擬的正確性。

七、 結論與反思 (O-R-I-D)

• O 客觀回顧 在這長達數個月的獨立研究中，我經歷了從純粹的紙本理論跨越到實體資訊工程的完整學術迴圈。專題初期，我為了建構系統的基礎，埋首研讀了大學普物與電磁學教材，經歷無數次計算紙上的推翻重來，才終於推導出高階的 \(1/r^4\) 磁力非線性偏微分方程式。隨後，為了將抽象數學具象化，我自學了 Python 程式語言，將牛頓運動定律降階，並親自撰寫四階龍格-庫塔法 (RK4) 數值積分演算法來模擬高維度的相空間。而在專題後期，我更親手裁切壓克力板、雷射切割支架，架設垂直俯視的高速攝影機，並利用開源的 Tracker 軟體完成真實物理數據的動態擷取與誤差驗證。
• R 深刻感受 這段歷程充滿了極致的挫折與隨後突破瓶頸的巨大狂喜。當我發現如此強烈的非線性方程式，根本無法像高中物理一樣求得一個「完美、精確」的解析解時，我曾感到極度無力。然而，當我轉換思維，使用「演算法數值積分」來逼近真相，並且在深夜的電腦螢幕前，看著自己親手寫的程式碼，跑了幾十分鐘後，終於在畫布上渲染出與國際頂尖論文中一模一樣的絕美「Wada 碎形圖案」時，那種「用程式碼具象化宇宙數學法則」的成就感，深深震撼了我的內心。此外，當實驗初期軌跡莫名其妙快速衰減，讓我陷入深深自我懷疑時，我也在揪出「冷次定律渦電流」這個硬體元凶後，獲得了無與倫比的釋懷與自信。
• I 意義啟發 這次經驗徹底粉碎了我過去在學校裡那種「紙上談兵、套公式解題」的學習框架。我深刻體會到，「完美的數學模型」與「真實的物理世界」之間存在著多麼巨大的鴻溝。在白紙上，\(1/r^4\) 不過是個乾淨的代數符號；但在現實實驗中，一個底座金屬材質的誤用所引發的隱形電磁阻尼，或是釋放位置微小如一毫米的人為誤差，都會因為「蝴蝶效應」而導致整個預測模型徹底崩潰。這讓我領悟到一個極為重要的工程哲學：理論推導固然迷人且必要，但嚴謹的「變因控制」意識，以及面對未知干擾時的「硬體除錯 (Debugging)」能力，才是科學家與頂尖工程師解決真實世界複雜問題的真正基石。
• D 未來展望 這次處理多變數耦合方程與龐大空間數值運算的經驗，不僅沒有擊退我，反而徹底確立了我對資訊工程與電機控制系統的熱忱。我發現自己非常享受這種「發現異常現象 ➔ 建立數學模型 ➔ 撰寫程式演算 ➔ 動手硬體除錯」的高階解題過程。未來進入大學殿堂，我期望能不再侷限於單一學科，而是持續精進高階資料結構、演算法最佳化設計與數位訊號處理能力，並期許自己能將這種分析與解決複雜非線性系統的運算邏輯，實際應用於 AI 物理引擎建模、精密自動控制系統，或是更廣泛的跨領域科技前沿中。

八、 參考文獻

1. DYNAMICAL ANALYSIS OF A PARAMETRICALLY FORCED MAGNETO-MECHANICAL OSCILLATOR | The ANZIAM Journal | Cambridge Core, 檢索日期：2月 28, 2026.
2. Experimental Solutions E1 - Magnetic Pendulum - EuPhO, 檢索日期：2月 28, 2026.
3. (PDF) Chaos in a Magnetic Pendulum Subjected to Tilted Excitation and Parametric Damping - ResearchGate, 檢索日期：2月 28, 2026.
4. CaYPT Problem H. Magnetic Assist, 檢索日期：2月 28, 2026.
5. April | 2012 | Softology's Blog - WordPress.com, 檢索日期：2月 28, 2026.
6. Chaotic Oscillating Magnetic Pendulum Simulation | Softology's Blog - WordPress.com, 檢索日期：2月 28, 2026.
7. An Analytic Solution for the Force Between Two Magnetic Dipoles - ResearchGate, 檢索日期：2月 28, 2026.
8. Chaos in the Magnetic Pendulum - IMA, 檢索日期：2月 28, 2026.
9. An analytic solution for the torque between two magnetic ... - SciSpace, 檢索日期：2月 28, 2026.
10. Magnetic field of a dipole and the dipole–dipole interaction - Rose-Hulman, 檢索日期：2月 28, 2026.
11. SYPT A6: A6. Magnetic Assist ... Investigate how the motion of the moving magnet depends on relevant parameters., 檢索日期：2月 28, 2026.
12. An Analytic Solution for the Torque Between Two Magnetic Dipoles - ResearchGate, 檢索日期：2月 28, 2026.
13. An Analytic Solution for the Torque Between Two Magnetic Dipoles - Scite.ai, 檢索日期：2月 28, 2026.
14. Limitations of Calculating Field Distributions and Magnetic Susceptibilities in MRI using a Fourier Based Method - PMC, 檢索日期：2月 28, 2026.
15. Modeling the interactions between magnetic particles - Diva-Portal.org, 檢索日期：2月 28, 2026.
16. Chaos in the Magnetic Pendulum | IMA, 檢索日期：2月 28, 2026.
17. DYNAMICAL ANALYSIS OF A PARAMETRICALLY FORCED MAGNETO-MECHANICAL OSCILLATOR | The ANZIAM Journal | Cambridge Core, 檢索日期：2月 28, 2026.
18. Experimental Solutions E1 - Magnetic Pendulum - EuPhO, 檢索日期：2月 28, 2026.
19. (PDF) Chaos in a Magnetic Pendulum Subjected to Tilted Excitation and Parametric Damping - ResearchGate, 檢索日期：2月 28, 2026.
20. CaYPT Problem H. Magnetic Assist, 檢索日期：2月 28, 2026.
21. April | 2012 | Softology's Blog - WordPress.com, 檢索日期：2月 28, 2026.
22. 檢索日期：2月 28, 2026，https://weelookang.blogspot.com/2025/02/sypt-a6-a6-magnetic-assist-attach-one.html
23. (PDF) Implementation of external magnetic field to improve strength of St37 steel resistance spot weld - ResearchGate, 檢索日期：2月 28, 2026.
24. 2025 IYPT 9. Magnetic Assist Unlocking the Secrets of Magnet Motion! - YouTube, 檢索日期：2月 28, 2026.