方程求解的实验 matlab,matlab 实验四 求微分方程的解

实际应用问题通过数学建模所归纳而得到的方程，绝大多数都是微分方程，真正能得到代数方程的机会很少．另一方面，能够求解的微分方程也是十分有限的，特别是高阶方程和偏微分方程(组)．这就要求我们必须研究微分方程(组)的解法，既要研究微分方程(组)的解析解法(精确解)，更要研究微分方程(组)的数值解法(近似解)．

对微分方程(组)的解析解法(精确解)，Matlab有专门的函数可以用，本实验将作一定的介绍．

本实验将主要研究微分方程(组)的数值解法(近似解)，重点介绍Euler折线法．

1．dsolve('equ1','equ2',…)：Matlab求微分方程的解析解．equ1、equ2、…为方程(或条件)．写方程(或条件)时用Dy表示y关于自变量的一阶导数，用用D2y表示y关于自变量的二阶导数，依此类推．

2．simplify(s)：对表达式s使用maple的化简规则进行化简．

例如：

syms x

simplify(sin(x)^2 + cos(x)^2)

ans=1

3．[r,how]=simple(s)：由于Matlab提供了多种化简规则，simple命令就是对表达式s用各种规则进行化简，然后用r返回最简形式，how返回形成这种形式所用的规则．

例如：

syms x

[r,how]=simple(cos(x)^2-sin(x)^2)

r = cos(2*x)

how = combine

4．[T,Y] = solver(odefun,tspan,y0)求微分方程的数值解．

说明：

(1)其中的solver为命令ode45、ode23、ode113、ode15s、ode23s、ode23t、ode23tb之一．

(2) odefun是显式常微分方程：

(3)在积分区间tspan=

上，从到，用初始条件求解．

(4)要获得问题在其他指定时间点上的解，则令tspan=

(要求是单调的)．

(5)因为没有一种算法可以有效地解决所有的ODE问题，为此，Matlab提供了多种求解器Solver，对于不同的ODE问题，采用不同的Solver．

求解器

Solver

ODE类型

特点

说明

ode45

非刚性