第7章代数方程与微分方程的 计算机求解.ppt

MATLAB语言与应用 第 7 章代数方程和微分方程的 计算机求解 现代设计与分析研究所 王 雷 主要内容 代数方程的求解 微分方程的求解 7.1 代数方程的求解 代数方程的图解法 多项式型方程的准解析解法 一般非线性方程数值解 7.1.1 代数方程的图解法7.1.1.1 一元方程的图解法 7.1.1.2 二元方程的图解法 7.1.2 多项式型方程的准解析解法 【例7-4】 【例7-5】 【例7-6】 【例7-7】 7.1.3 一般非线性方程数值解 【例7-8】 【例7-9】 【例7-10】 7.2 微分方程的求解 常系数线性微分方程的解析解方法 微分方程问题的数值解法 7.2.1 常系数线性微分方程的解析解方法 线性常系数微分方程解析解的数学描述 微分方程的解析解方法 7.2.1.1 线性常系数微分方程解析解的数学描述 7.2.1.2 微分方程的解析解方法 【例7-11】 【例7-12】 7.2.1.3 特殊非线性微分方程的解析解 7.2.2 微分方程问题的数值解法 微分方程问题算法概述 一阶微分方程组的数值解 微分方程转换 7.2.2.1 微分方程问题算法概述 7.2.2.2 一阶微分方程组的数值解四阶五级Runge-Kutta-Felhberg算法 【例7-14】 【例7-15】 附加参数的微分方程求解 7.2.2.3 微分方程转换单个高阶常微分方程处理方法 【例7-17】 高阶常微分方程组的变换方法 【例7-18】 数学方程求解是科学与工程研究中经常遇到的问题，本章先介绍了简单方程的图解法，给出了方程求解的基本概念，并介绍了基于符号运算工具箱中 solve( ) 函数的多项式类方程的准解析求解算法，还介绍了fsolve( ) 函数求取一般非线性方程的数值解法。 介绍了基于符号运算工具箱 dsolve()函数的线性微分方程的解析解方法。引入了数值解的概念，并以最简单的一阶微分方程的 Euler 算法为例，介绍了一般数值解法的思路并介绍了变步长求解的概念，还介绍了 MATLAB 下的微分方程数值求解函数 ode45( )，通过例子演示了该函数的使用方法。 对一般非线性微分方程来说，解析解是不存在的，只能依赖数值解的方法对其进行研究。 微分方程初值函数能直接求解的方程是一阶显式微分方程组，若给出的方程不是这类函数，则需要通过本书介绍的方法选择一组状态变量，将原方程变换成一阶显式微分方程组，以便用给定的求解函数直接求解。 求解函数 【例7-16】 本章内容简介 微分方程的控制参数 odeset（） 四阶五阶Runge-Kutta-Fahlberg算法 ode45（） 常微分方程的解析解 dsolve（） 控制参数 optimset（） 方程的数值解 fsolve（） 方程的解析解，尤适于多项式方程 solve（） 函数功能 函数名 本章内容简介 本章内容简介 * * 【例7-1】 t=3.5203; vpa(exp(-3*t)*sin(4*t+2)+… 4*exp(-0.5*t)*cos(2*t)-0.5) 【例7-2】 【例7-3】 y0=lambertw(xx); plot(xx, y0) 处理系数后方程的解： 【例7-13】