微分方程数值解.pdfVIP

数学模型与实验（三） 微分方程数值解 常微分方程求解 dx/dt=rx(1-x/k) 1. 求表达式 （符号运算） S=dsolve(‘Dx=r*x*(1-x/k)’,’x(0)=0.01’) %转变为函数值 r=0.3;k=8; s=subs(S) t=0:40; ss=subs(s,’t’,t); plot(t,ss),grid 2. 求数值解 [tout,xout]=ode23(fun,[to,tfinal],x0): %调用ode23(二、三阶龙格库塔法)解微分方程 输入参数说明： fun:，被积分的函数名, 用函数文件fun.m定义： function y=fun(x,t) y=0.3*x*(1-x/8); [to, tfinal]=[t的初始值, t的终值] 。 xo ：x 的初始值。 输出参数说明： tout: 离散的自变量值，xout: 离散的函数值。 [tout,xout]=ode23(@fun,[0, 40],0.01) plot(tout,xout) 3. 其他的求解语句 ode45 ode113 ode 23s ode23t