《MatLab教学资料》06实验六微分方程.pptVIP

实验六 微分方程和微分方程组的求解 一、实验目的 熟悉MATLAB软件中关于求解微分方程和微分方程组的各种命令，掌握利用MATLAB软件进行微分方程和微分方程组的求解。 二、相关知识 在MATLAB中，由函数dsolve()解决微分方程（组）的求解问题，其具体格式如下： X=dsolve(‘eqn1’,’eqn2’,…) 函数dsolve用来解符号微分方程、方程组，如果没有初始条件，则求出通解，如果有初始条件，则求出特解。 例1：求解微分方程 的MATLAB程序为：dsolve(Dy=1/(x+y),x)，注意，系统缺省的自变量为t，因此这里要把自变量写明。 二、相关知识 结果为：-lambertw(-C1*exp(-x-1))-x-1 其中：Y=lambertw(X)表示函数关系Y*exp(Y)=X。 例2：求解微分方程 。 解：先把方程改写为 例1：求解微分方程 的MATLAB程序为： dsolve(Dy=1/(x+y),x)，注意，系统缺省的自变量为t，因此这里要把自变量写明。 结果为：-lambertw(-C1*exp(-x-1))-x-1 其中：Y=lambertw(X)表示函数关系Y*exp(Y)=X。 例2：求解微分方程 。 解：先把方程改写为 然后编写MATLAB程序： Y2=dsolve(D2y-Dy^2/y=0, x) 结果为： Y2 = exp(C1*x)*C2 例3：求微分方程组 通解的MATLAB程序为： [X,Y]=dsolve(Dx+5*x+y=exp(t),Dy-x-3*y=exp(2*t),t) 然后编写MATLAB程序： Y2=dsolve(D2y-Dy^2/y=0, x) 结果为： Y2 = exp(C1*x)*C2 例3：求微分方程组 通解的MATLAB程序为： [X,Y]=dsolve(Dx+5*x+y=exp(t),Dy-x-3*y=exp(2*t),t) 例4：求微分方程组 特解的MATLAB程序为： [X,Y]=dsolve(Dx+2*x-Dy=10*cos(t), Dx+Dy+2*y=4*exp(-2*t),x(0)=2,y(0)=0) 以上这些都是微分方程的精确解法，也称为微分方程的符号解。我们知道，有大量的微分方程虽然从理论上讲，其解是存在的，但我们却无法求出其解析解，此时，我们需要寻求方程的数值解，在求微分方程数值解方面，MATLAB具有丰富的函数，我们将其统称为solver，其一般格式为： 例4：求微分方程组 特解的MATLAB程序为： [X,Y]=dsolve(Dx+2*x-Dy=10*cos(t), Dx+Dy+2*y=4*exp(-2*t),x(0)=2,y(0)=0) 以上这些都是微分方程的精确解法，也称为微分方程的符号解。我们知道，有大量的微分方程虽然从理论上讲，其解是存在的，但我们却无法求出其解析解，此时，我们需要寻求方程的数值解，在求微分方程数值解方面，MATLAB具有丰富的函数，我们将其统称为solver，其一般格式为： [T,Y]=solver(odefun,tspan,y0) 该函数表示在区间tspan=[t0,tf]上，用初始条件y0求解显式常微分方程 。 solver为命令ode45，ode23，ode113，ode15s，ode23s，ode23t，ode23tb之一，这些命令各有特点。我们列表说明如下： [T,Y]=solver(odefun,tspan,y0) 该函数表示在区间tspan=[t0,tf]上，用初始条件y0求解显式常微分方程 。 solver为命令ode45，ode23，ode113，ode15s，ode23s，ode23t，ode23tb之一，这些命令各有特点。我们列表说明如下： 求解器 ODE类型 特点 说明 ode45 非刚性 一步算法，4,5阶Runge-Kutta方法累积截断误差 大部分场合的首选算法 ode23 非刚性 一步算法，2,3阶Runge-Kutta方法累积截断误差 使用于精度较低的情形 求解器 ODE类型 特点 说明 ode45 非刚性 一步算法，4,5阶Runge-Kutta方法累积截断误差 大部分场合的首选算法 ode23 非刚性 一步算法，2,3阶Runge-Kutta方法累积截断误差 使用于精度较低的情形 ode113 非刚性 多步法，Adams算法，高低精度均可达到 计算时间比ode45短 ode23t 适度刚性 采用梯形算法 适度刚性情形 o