【正文】 常微分方程的高精度求解方法 安徽大學江淮學院 07 計算機 (1)班 安徽大學江淮學院 本科畢業(yè)論文（設計） 題 目： 常 微分方程求解的高階方法 學生姓名： 圣近 學號： JB074219 院（系）： 計算機科學與技術 專業(yè)： 計算機科學與技術 入學時間： 2021 年 9 月 導師姓 名： 汪繼文 職稱 /學位： 教授 導師所在單位： 安徽大學計算機科學與技術學院 常微分方程的高精度求解方法 安徽大學江淮學院 07 計算機 (1)班 I 常微分方程的高精度求解方法 摘 要 本文主要討論了常微分方程的高精度求解方法的相關解法問題。文章首先 案例引入微分方程概念，然后 給出了微分方程的基本概念?？茖W和工程中建立數(shù)學模型時常用到微分方程。由于它們通常沒有已知的解析解，因而需要求其數(shù)值近似解。先從常微分方程解析解法出發(fā)，分析解析解法在實際運用中的局限性，引入常微分方程的數(shù)值解法，呈現(xiàn)常微分方程數(shù)值求解的三個步驟：將問題離散化，建立或尋找一個遞推格式，按步進方式計算。再從對精度需求出發(fā)從低階數(shù)值方法到高階數(shù)值方法進行逐步的探討，分析 各種方法的數(shù)學原理，闡述其推導方法，比較不同方法的優(yōu)缺點，重點介紹實用的龍格 — 庫塔方法、歐拉方法、休恩方法、泰勒級數(shù)法和預報 — 校正方法，并以編寫相應程序作總結。最后，再討論高階常微分方程和一階常微分方程組：一般的高階常微分方程都可以通過相應的變量代換轉化為一階常微分方程組，一階常微分方程的初值問題求數(shù)值解與一階常微分方程的初值問題求數(shù)值解的方法基本相同。 關鍵詞：龍格 — 庫塔方法；歐拉方法；休恩方法；泰勒級數(shù)法；預報 — 校正方法； 常微分方程的高精度求解方法 安徽大學江淮學院 07 計算機 (1)班 II High accuracy method for solving ordinary differential equations Abstract This paper discusses the accuracy method for solving ordinary differential equations related to solution problems. The article first case to introduce the concept of differential equations, and then gives the basic concepts of differential equations. Science and engineering often use a mathematical model equations. As they often do not have known analytic solution, and thus demand for its numerical approximate solution. Start with the analytical solution of ordinary differential equations, analyzes the analytical solution of the limitations in the practical application, the introduction of numerical solution of ordinary differential equations, numerical solution of ordinary differential equations presented in three steps: discretization of the problem, create or find a recursive format is calculated by step. Starting from the demand for accuracy and then from low to high numerical method of step by step numerical method to analyze various methods of mathematical theory to explain their derivation, pare the advantages and disadvantages of different methods, focusing on practical Runge Kutta method, Euler method, Bethune method, Taylor series method and prediction correction methods, and procedures for the preparation of the corresponding summary. Finally, discuss the higher order ordinary differential equations and first order ordinary differential equations: general higher order ordinary differential equations can be substituted by the corresponding variable into a firstorder ordinary differential equations, first order initial value problems for ordinary differential equations numerical solution with an initial value problem of differential equation numerical solution method is basically the same. Keywords: Runge Kutta methods。 Euler method。 Huon method。 Taylor series method。 prediction correction。 常微分方程的高精度求解方法 安徽大學江淮學院 07 計算機 (1)班 III 目 錄 第一章 前 言 .................................................. 1 案例引入微分方程概念 .................................... 1 微分方程的基本概念 ...................................... 1 微分方程及微分方程的階 ............................. 1 微分方程的解、通解與特解 ............................ 1 微分方程的初值條件及其提法 .......................... 2 微分方程的解的幾何意義 .............................. 2 從解析方法到數(shù)值方法概述 ................................. 3 常溫分方程的離散化 ...................................... 4 第二章 數(shù)值解法公共程序模塊分析 ................................. 5 第三章 歐拉（ Euler）方法 ....................................... 7 Euler 方法思想 .......................................... 7 Euler 方法的誤差估計 .................................... 8 改進的 Euler 方法 ........................................ 8 梯形公式 .......................................... 8 改進 Euler 法 ....................................... 9 第四章 休恩方法 .............................................. 10 休恩方法思想 .......................................... 10 休恩方法的步長和誤差 ................................... 10 第五章 泰勒級數(shù)法 ............................................ 11 泰勒定理 .............................................. 11 N 次泰勒方法 ........................................... 12 第六章 龍格 庫塔（ Runge— Kutta 法） ............................ 13 龍格 庫塔（ Runge— Kutta）方法基本思想 .................... 13 階龍格 庫塔（ Runge— Kutta）方法公式 ..................... 14 第七章 預報 校正方法 .......................................... 15 MilneSimpon 方法 ...................................... 16 誤差估計于校正 ......................................... 16 正確的步長 ............................................ 17 第八章 一階微分方程組與高階微分方程的數(shù)值解法 ................... 17 一階微分方程組的數(shù)值解法 ............................... 17 高階微分方程的數(shù)值解法 ................................ 18 第九章 常微分方程模型數(shù)值解法在數(shù)學建模中的應用 ................. 19 常微分方程的高精度求解方法 安徽大學江淮學院 07 計算機 (1)班 IV 耐用消費新產品的銷售規(guī)律模型 ............................ 19 問題的提出 ....................................... 19 模型的構建 ....................................... 19 模型的求解 ....................................... 20 司機飲酒駕車防避模型的數(shù)值解法 .......................... 21 模型假設 ......................................... 22 模型建立 ......................................... 22 模型求解 ......................................... 24 模型評價 ......................................... 25 誠懇建議 ......................................... 25 模型推廣 ......................................... 26 主要參考文獻 ................................................. 26 致 謝 ..................................................... 27 常微分方程的高精度求解方法 安徽大學江淮學院 07 計算機 (1)班 1 第一章 前 言 案例引入微分方程概念 在科技、工程、經濟管理、生態(tài)、生態(tài)、刑偵等各個領域微分方程有著廣泛的應用。我們看一實例。 案例：一次謀殺案，在某天下午四點發(fā)現(xiàn)尸體，尸體的體溫為 30℃，假設當時屋內空間的溫度保護 20℃不變，現(xiàn)判斷謀殺是何時發(fā)生的？ 解決此問題首先必須要從尸體溫度的變化尋求關系式，這就需要知道物理學中的加熱與冷卻規(guī)律。物理學家牛頓（ Newton） 曾提出，一塊熱的物體，其溫度下降的速度是與它自身溫度的差值成正比。同樣，一塊冷的物體，其溫度