用Matlab编程计算二阶常微分方程组的数值解的方法有很多种,最常用有ode45函数。
ode45函数的使用格式:
【t,y】=ode45(odefun,tspan,y0)
例如:
>> odefun=@(t,y)[y(1)+2*y(2);3*y(1)+2*y(2)]; %自定义函数
>> tspan=[0 20]; %时间
>> y0=[1 1]; %初值
>> [t,y]= ode45(odefun,tspan,y0);
>> A=[t,y];
>> A
运行结果