實驗三常微分方程式(I).docVIP

實驗三 常微分方程式(I) 目的： 本實驗之主要目的在訓練同學們學習使用Euler，Euler-Cromer及Verlet等方法來求解常微分方程式。我們將以二維拋體運動，單擺運動為例，來實作上述方法之應用。 原理： 設想我們要描述一質點的運動狀態，我們將引入位置及速度等狀態的變數，藉由描述與的時間相關性，便能掌握期運動的軌跡。由於牛頓的貢獻，使我們瞭解到，速度的變化率( 即加速度 ) 如何因物體所受之力而改變，即牛頓第二運動定律： ( 3-1 ) ( 3-2 ) 只要我們能測得物體所受的力，即及初使( t=0 ) 的運動狀態，則可完全掌握其軌跡。本實驗及企圖使用數值計算的方法來獲取其運動之軌跡及性質，廣義來說，就是需求取常微分方程式之數值解，不僅限於物理中的力學問題。 ( 一 ) Euler Method 我們首先考慮一個最簡單的例子，以上一個實驗中學得的差商的觀念應用到常微分方程式之求解，即，其中是變數t ( 可視為時間 ) 的增量。利用此一近似，我們可以將牛頓運動方程式寫為 ( 3-3 )