【正文】 Runge— Kutta）方法公式 ..................... 14 第七章 預(yù)報(bào) 校正方法 .......................................... 15 MilneSimpon 方法 ...................................... 16 誤差估計(jì)于校正 ......................................... 16 正確的步長(zhǎng) ............................................ 17 第八章 一階微分方程組與高階微分方程的數(shù)值解法 ................... 17 一階微分方程組的數(shù)值解法 ............................... 17 高階微分方程的數(shù)值解法 ................................ 18 第九章 常微分方程模型數(shù)值解法在數(shù)學(xué)建模中的應(yīng)用 ................. 19 常微分方程的高精度求解方法 安徽大學(xué)江淮學(xué)院 07 計(jì)算機(jī) (1)班 IV 耐用消費(fèi)新產(chǎn)品的銷售規(guī)律模型 ............................ 19 問(wèn)題的提出 ....................................... 19 模型的構(gòu)建 ....................................... 19 模型的求解 ....................................... 20 司機(jī)飲酒駕車防避模型的數(shù)值解法 .......................... 21 模型假設(shè) ......................................... 22 模型建立 ......................................... 22 模型求解 ......................................... 24 模型評(píng)價(jià) ......................................... 25 誠(chéng)懇建議 ......................................... 25 模型推廣 ......................................... 26 主要參考文獻(xiàn) ................................................. 26 致 謝 ..................................................... 27 常微分方程的高精度求解方法 安徽大學(xué)江淮學(xué)院 07 計(jì)算機(jī) (1)班 1 第一章 前 言 案例引入微分方程概念 在科技、工程、經(jīng)濟(jì)管理、生態(tài)、生態(tài)、刑偵等各個(gè)領(lǐng)域微分方程有著廣泛的應(yīng)用。 Taylor series method。 Euler method。最后，再討論高階常微分方程和一階常微分方程組：一般的高階常微分方程都可以通過(guò)相應(yīng)的變量代換轉(zhuǎn)化為一階常微分方程組，一階常微分方程的初值問(wèn)題求數(shù)值解與一階常微分方程的初值問(wèn)題求數(shù)值解的方法基本相同。先從常微分方程解析解法出發(fā)，分析解析解法在實(shí)際運(yùn)用中的局限性，引入常微分方程的數(shù)值解法，呈現(xiàn)常微分方程數(shù)值求解的三個(gè)步驟：將問(wèn)題離散化，建立或?qū)ふ乙粋€(gè)遞推格式，按步進(jìn)方式計(jì)算?？茖W(xué)和工程中建立數(shù)學(xué)模型時(shí)常用到微分方程。常微分方程的高精度求解方法 安徽大學(xué)江淮學(xué)院 07 計(jì)算機(jī) (1)班 安徽大學(xué)江淮學(xué)院 本科畢業(yè)論文（設(shè)計(jì)） 題 目： 常 微分方程求解的高階方法 學(xué)生姓名： 圣近 學(xué)號(hào)： JB074219 院（系）： 計(jì)算機(jī)科學(xué)與技術(shù) 專業(yè)： 計(jì)算機(jī)科學(xué)與技術(shù) 入學(xué)時(shí)間： 2021 年 9 月 導(dǎo)師姓 名： 汪繼文 職稱 /學(xué)位： 教授 導(dǎo)師所在單位： 安徽大學(xué)計(jì)算機(jī)科學(xué)與技術(shù)學(xué)院 常微分方程的高精度求解方法 安徽大學(xué)江淮學(xué)院 07 計(jì)算機(jī) (1)班 I 常微分方程的高精度求解方法 摘 要 本文主要討論了常微分方程的高精度求解方法的相關(guān)解法問(wèn)題。文章首先 案例引入微分方程概念，然后 給出了微分方程的基本概念。由于它們通常沒(méi)有已知的解析解，因而需要求其數(shù)值近似解。再?gòu)膶?duì)精度需求出發(fā)從低階數(shù)值方法到高階數(shù)值方法進(jìn)行逐步的探討，分析 各種方法的數(shù)學(xué)原理，闡述其推導(dǎo)方法，比較不同方法的優(yōu)缺點(diǎn)，重點(diǎn)介紹實(shí)用的龍格 — 庫(kù)塔方法、歐拉方法、休恩方法、泰勒級(jí)數(shù)法和預(yù)報(bào) — 校正方法，并以編寫相應(yīng)程序作總結(jié)。 關(guān)鍵詞：龍格 — 庫(kù)塔方法；歐拉方法；休恩方法；泰勒級(jí)數(shù)法；預(yù)報(bào) — 校正方法； 常微分方程的高精度求解方法 安徽大學(xué)江淮學(xué)院 07 計(jì)算機(jī) (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。 Huon method。 prediction correction。我們看一實(shí)例。物理學(xué)家牛頓（ Newton） 曾提出，一塊熱的物體，其溫度下降的速度是與它自身溫度的差值成正比。據(jù)此我們可找到溫度與時(shí)間之間的函數(shù)關(guān)系式，這事實(shí)上就是一個(gè)微分方程的建立問(wèn)題。通過(guò)求解微分方程，可以得到所需求的函數(shù)。?????????驗(yàn) 證 函 數(shù) ， 為 任 意 常 數(shù) 是二 階 微 分 方 程的 通 解 ， 并 求 此 微 分 方 程 滿 足 初 值 條 件 ：).的 特 解 3 Cxy ??例如 .d3d 2 的通解是 xxy ?.dd 22 的通解是 gt s ? 21 212 CtCgts ???又如 13 ?? xy例如 .d3d 2 的特解是 xxy ?221 gts ?又如 .dd 22 的特解是 gts ?.,)( ,| )( ,| 00000000000都是已知值其中或，或記作y39。xy39。y39。xy39。 CCy 22212221e2e2 ee ???? ?? ，xx CCy 2221 e4e4 ??? 2 2 2 21 2 1 2( 1 . 3 ) 4 4 e 4 e 4 e 4 e 0x x x xy y C C C C??? ? ? ? ? ?把 它 們 代 入 微 分 方 程 的 左 端 ， 得 所以函數(shù) xx CCy 2221 ee ??? 是所給微分方程 ()的解 .又因?yàn)檫@個(gè)解中含有兩個(gè)獨(dú)立的任意常數(shù)，任意常數(shù)的個(gè)數(shù)與微分方程 ()的階數(shù)相同，所以它是該方程的通解 . 0022122212( ) : 0 39。 2 e 2 exxxxxxyyy C Cy C C??????????把 式 中 的 條 件 “ ” 及 “ ” 分 別代 入及中 ， 得 .41411220212121??????????CCCCCC，解得， ).ee(41 22 xxy ???初值條件的特解為于是所求微分方程滿足 從解 析方法到數(shù)值方法概述 求解常微分方程的解析方法很多，像變量分離法，積分因子法，遺憾的是實(shí)際上得到的大部分常微分方程都不能使用這些理論上的方法。 常微分方程的高精度求解方法 安徽大學(xué)江淮學(xué)院 07 計(jì)算機(jī) (1)班 4 常溫分方程的離散化 下面主要討論 一階常微分方程的初值問(wèn)題，其一般形式是 0( , ) ( )()dy f x y a x bdxy a y? ? ? ???? ?? 在下面的討論 我們總假定函數(shù) f (x, y) 連續(xù)，且關(guān)于 y 滿足李普希茲(Lipschitz)條 件，即存在常數(shù) L ，使得 ( , ) ( , )f x y f x y L y y? ? ? 這樣，由常微分方程理論知，初值問(wèn)題 (1)的解必定存在唯一。今后如無(wú)特別說(shuō)明，我們總?cè)〔介L(zhǎng)為常量 h 。式（ 1,7）是個(gè)離散化的問(wèn)題，稱為差分方程初值問(wèn)題。 常微分方程的高精度求解方法 安徽大學(xué)江淮學(xué)院 07 計(jì)算機(jī) (1)班 5 （ ii）用數(shù)值積分方法 將問(wèn)題（ ）的解表成積分形式，用數(shù)值積分方法離散化。 （ iii） Taylor 多項(xiàng)式近似 將函數(shù) ()yx 在 nx 處展開，取一次 Taylor 多項(xiàng)式近似，則得