Matlab使用ode45解微分方程组的基本问题,下面有一个简单的微分方程组求代码
Matlab使用ode45解微分方程组的基本问题,下面有一个简单的微分方程组求代码
日期:2016-04-13 11:27:57 人气:1
建立myfun.m文件,把以下代码复制进去
function dx=myfun(t,x)dx=zeros(2,1);dx(1)=4*x(1)+x(2).^2;dx(2)=8*x(1)+2*x(2).^2+5*x(2);运行以下代码
t=[0 0.1];
x0=[1 1];
[tt,xx]=ode45('myfun',t,x0);
plot(tt,xx(:,1),'-',tt,xx(:,2),'-.')%画出两者图像