MATLAB求解常微分方程数值解.docxVIP

PAGE \* MERGEFORMAT3 利用MATLAB求解常微分方程数值解 目录 TOC \o 1-3 \h \z \u HYPERLINK \l _Toc391567541 1. 内容简介 PAGEREF _Toc391567541 \h 1 HYPERLINK \l _Toc391567542 2. Euler Method(欧拉法)求解 PAGEREF _Toc391567542 \h 1 HYPERLINK \l _Toc391567543 2.1. 显式Euler法和隐式Euler法 PAGEREF _Toc391567543 \h 2 HYPERLINK \l _Toc391567544 2.2. 梯形公式和改进Euler法 PAGEREF _Toc391567544 \h 3 HYPERLINK \l _Toc391567545 2.3. Euler法实用性 PAGEREF _Toc391567545 \h 5 HYPERLINK \l _Toc391567546 3. Runge-Kutta Method(龙格库塔法)求解 PAGEREF _Toc391567546 \h 6 HYPERLINK \l _Toc391567547 3.1. Runge-Kutta基本原理 PAGEREF _Toc391567547 \h 6 HYPERLINK \l _Toc391567548 3.2. MATLAB中使用Runge-Kutta法的函数 PAGEREF _Toc391567548 \h 8 HYPERLINK \l _Toc391567549 4. 使用MATLAB求解常微分方程 PAGEREF _Toc391567549 \h 8 HYPERLINK \l _Toc391567550 4.1. 使用ode45函数求解非刚性常微分方程 PAGEREF _Toc391567550 \h 8 HYPERLINK \l _Toc391567551 4.2. 刚性常微分方程 PAGEREF _Toc391567551 \h 9 HYPERLINK \l _Toc391567552 5. 总结 PAGEREF _Toc391567552 \h 10 HYPERLINK \l _Toc391567553 参考文献 PAGEREF _Toc391567553 \h 11 HYPERLINK \l _Toc391567554 附录 PAGEREF _Toc391567554 \h 12 HYPERLINK \l _Toc391567555 1. 显式Euler法数值求解 PAGEREF _Toc391567555 \h 12 HYPERLINK \l _Toc391567556 2. 改进Euler法数值求解 PAGEREF _Toc391567556 \h 12 HYPERLINK \l _Toc391567557 3. 四阶四级Runge-Kutta法数值求解 PAGEREF _Toc391567557 \h 13 HYPERLINK \l _Toc391567558 4. 使用ode45求解 PAGEREF _Toc391567558 \h 14 PAGE \* MERGEFORMAT14 PAGE \* MERGEFORMAT12 内容简介 把《高等工程数学》看了一遍，增加对数学内容的了解，对其中数值解法比较感兴趣，这大概是因为在其它各方面的学习和研究中经常会遇到数值解法的问题。理解模型然后列出微分方程，却对着方程无从下手，无法得出精确结果实在是让人难受的一件事情。 实际问题中更多遇到的是利用数值法求解偏微分方程问题，但考虑到先从常微分方程下手更为简单有效率，所以本文只研究常微分方程的数值解法。把一个工程实际问题弄出精确结果远比弄清楚各种细枝末节更有意思，因此文章中不追求非常严格地证明，而是偏向如何利用工具实际求解出常微分方程的数值解，力求将课程上所学的知识真正地运用到实际方程的求解中去，在以后遇到微分方程的时候能够熟练运用MATLAB得到能够在工程上运用的结果。 文中求解过程中用到MATLAB进行数值求解，主要目的是弄清楚各个函数本质上是如何对常微分方程进行求解的，对各种方法进行MATLAB编程求解，并将求得的数值解与精确解对比，其中源程序在附录中。最后考察MATLAB中各个函数的适用范围 ，当遇到实际工程问题时能够正确地得到问题的数值解。 Euler Method(欧拉法)求解 Euler法求解常微分方程主要包括3种形式，即显式Euler法、隐式Euler法、梯形公式法，本节内容分别介绍这3种方法的具体内容，并在最后