常微分方程课程设计课件.pptVIP

数学实验 求微分方程的解 问题背景和实验目的 Euler 折线法 初值问题的Euler折线法 Euler 折线法举例 Euler 折线法源程序 Euler折线法举例（续） Runge-Kutta 方法 Runge-Kutta 方法 四阶 R-K 方法源程序 Runge-Kutta 方法 Euler 法与 R-K法误差比较 Matlab 解初值问题 dsolve 求解析解 dsolve 的使用 dsolve 举例 dsolve 举例 dsolve 举例 Matlab函数数值求解 参数说明 数值求解举例 数值求解举例 Matlab 求解微分方程小结 ************************************ 想：无论是解析解还是数值解，都不如图形解直观明了。在MATLAB下用 什么命令可以用图形显示解析解或数值解？ 提示：如果微分方程（组）的解析解为y=f(x)，则可以用MATLAB函数fplot作出函数曲线y=f(x)的图形，具体用法 其中fun给出函数表达式；lims=[xmin xmax ymin ymax]限定坐标轴的大小．例如 如果得到微分方程（组）的数值解（x，y）．用MATLAB函数plot（x，y）直接作出图形，其中x和y为向量（或矩阵） *********************************** 上机练习（不是大作业） * * 自牛顿发明微积分以来，微分方程在描述事物运动规律上已发挥了重要的作用。实际应用问题通过数学建模所得到的方程，绝大多数是微分方程。 由于实际应用的需要，人们必须求解微分方程。然而能够求得解析解的微分方程十分有限，绝大多数微分方程需要利用数值方法来近似求解。 本实验主要研究如何用 Matlab 来计算微分方程（组）的数值解，并重点介绍一个求解微分方程的基本数值解法－－Euler折线法。 考虑一维经典初值问题 基本思想：用差商代替微商 根据 Talyor 公式，y(x) 在点 xk 处有 具体步骤： 等距剖分： 步长： 分割求解区间 差商代替微商 得方程组： 分割求解区间，差商代替微商，解代数方程 为分割点 k = 0, 1, 2, ..., n-1 yk 是 y (xk) 的近似 例：用 Euler 法解初值问题 取步长 h = (2 - 0)/n = 2/n，得差分方程 当 h=0.4，即 n=5 时，Matlab 源程序见 fulu1.m 解： clear f=sym(y+2*x/y^2); a=0; b=2; h=0.4; n=(b-a)/h+1; % n=(b-a)/h; x=0; y=1; szj=[x,y]; for i=1:n-1 % i=1:n y=y+h*subs(f,{‘x’,‘y’},{x,y});%第一次代初值 x=x+h; szj=[szj;x,y]; end szj plot(szj(:,1),szj(:,2),or-) 解析解： 解析解 近似解 为了减小误差，可采用以下方法： 让步长 h 取得更小一些； 改用具有较高精度的数值方法： 龙格-库塔方法 Runge-Kutta (龙格-库塔) 方法 是一类求解常微分方程的数值方法 有多种不同的迭代格式 用得较多的是 四阶R-K方法（教材第 79 页） 其中 clear; f=sym(y+2*x/y^2); a=0; b=2; h=0.4; n=(b-a)/h+1; % n=(b-a)/h; x=0; y=1; szj=[x,y]; for i=1:n-1 % i=1:n l1=subs(f,{x,y},{x,y}); l2=subs(f,{x,y},{x+h/2,y+l1*h/2}); l3=subs(f,{x,y},{x+h/2,y+l2*h/2}); l4=subs(f,{x,y},{x+h,y+l3*h}); y=y+h*(l1+2*l2+2*l3+l4)/6; x=x+h; szj=[szj;x,y]; end plot(szj(:,1),szj(:,2), dg-) 用 Maltab自带函数 解初值问题 求解析解：dsolve 求数值解：  ode45、ode23、 ode113、ode23t、ode15s、 ode23s、ode23tb dsolve 的使用 y=dsolve(eq1,eq2, ... ,cond1,cond2, ... ,v) 其中 y 为输出， eq1、eq2、...为微分方程，cond1、cond2、...为初值条件，v 为自变量。 例 1：求微分方程