利用常微分方程数值解函数 ode( ),可以得到 y(t)的数值解。其求解方法:
1、建立常微分方程降阶方程组函数,myodefun(t,y)
2、利用 ode45()函数求解。该函数计算比其他函数(如ode15s, ode23s)精度要高
[t,y] = ode45(@myodefun,tspan,y0,options)
3、利用plot()函数绘图。绘出t—y(t)曲线图,y(t)—y’(t)相平面图
4、利用find()函数计算y(0.5)的值。k = find(t==0.5);y=y(k)
5、运行结果,y(t)值
y(0.5)=1.7351