MATLAB函数微分方程的数值解

微分方程的数值解
求解微积分方程你的数值解常用的MATLAB函数调用如下：
[t, x] = ode23(‘xprime’, t0, tf, x0, tol, trace)
[t, x] = ode45(‘xprime’, t0, tf, x0, tol, trace)
或
[t, x] = ode23(‘xprime’, [t0, tf], x0, tol, trace)
[t, x] = ode45(‘xprime’, [t0, tf], x0, tol, trace)
说明：
（1）两个指令的调用格式相同，均为Runge-Kutta法。
（2）该指令是针对一阶常微分设计的。因此，假如待解的是高阶微分方程，那么它必修先演化为形如 的一阶微分方程组，即“状态方程”。
（3）’xprime’是定义 的函数名。该函数文件必须以 为一个列向量输出，以 为输入参量（注意输入变量之间的 关系，先“时间变量”后“状态变量”）。
（4）输入参量 和 分别是积分的起始值和终止值。
（5）输入参量 为初始状态列向量。
（6）输出参量 和 分别给出“时间”向量和相应的状态向量。
（7）tol控制解的精度，可缺省。缺省时，ode23默认tol=1.e-3；ode23默认tol=1.e-6。
（8）输入参量trace控制求解的中间结果是否显示，可缺省，缺省时，默认为tol=0，不显示中间结果。
（9）一般地，两者分别采用自适应变步长࿰