在matlab中求解一个常微分方程组！

来源：百度知道 编辑：UC知道 时间：2024/05/24 16:00:57
w'(t)=-k*f(t)
f'(t)=g*f(t)*[m*w(t)-f(t)]

用maple或者mathematic也行
要有代码哦。。

谢谢～！！！！

建立一个文件test.m

function dx=test(t,x)
dx=zeros(2,1);

dx(1)=k*x(2);
dx(2)=g*x(2)(m*x(1)-x(2));

保存,

在matlab里输入:

[t,x]=ode45('test',[0 15],[25 2]);

maple 的
restart:with(PDEtools):
alias(w=w(t),f=f(t)):
eq1:=diff(w,t)=k*f;
eq2:=diff(f,t)=g*f*(m*w-f);
dsolve({w,f});
程序是这样的，但是方程我试过了，不能直接给解。不是maple的问题，因为方程本身可能不可积，是数学问题。建议先分析一下可积性。
鉴定完毕