关于ode45函数的说明.pdfVIP

关于ode45函数的说明 求解器solver 在求解常微分方程方面，matlab具有丰富的函 数，将其统称为solver ，其一般格式为： [T,Y]=solver(odefun,tspan,y0) 该函数表示在区间tspan=[t0,t1]上，用初始条件y0 求解显式常微分方程y’=f(t,y) 。 odefun为显式常微分方程中的f(t,y)函数， tspan为求解区间。要获得问题在制定时刻的解， 可以令tspan=t0,t1,t2… （单调递增或递减），y0为 初始条件。 非刚性ODE求解命令 求解器solver 功能 说明 ode45 一步算法；4 ，5阶龙格库塔方程； 大部分尝试的首 累计截断误差(Δx)3 选算法 ode23 一步算法；2，3阶龙哥库塔方程； 适用于精度较低 累计截断误差(Δx)3 的情形 ode113 多步算法；Adams算法；高低精 计算时间比 度均可到10-6~10-3 ode45短 刚性ODE求解命令 求解器solver 功能 说明 ode23t 梯形算法 适度刚性情形 ode15s 多步法；Gear’s反向数值微分； 若ode45失效时，可 精度中等 尝试使用 ode23s 一步法；2阶Rosebrock算法； 当精度较低时，计算 精度低 时间比ode15s短 ode23tb 梯形算法；精度低 当精度较低时，计算 时间比ode15s短 [T,Y]=ode45(‘odefun’,tspan,y0) tspan为求解区间 y0为初始条件 （1）根据问题所属学科中的规律、定律、公式， 用微分方程与初始条件进行描述。 n     微分方程 F y , y , y ,..., y ,t  0 n1 初始条件 y 0 c , y  0 c ,..., y   0 c   0   1   n1 （2 ）若微分方程的阶数（n阶）大于1，则作变量 替换： n1     y1 y , y 2