作业帮 > 综合 > 作业

急求MATLAB解常微分方程组数值解的解决方案

来源：学生作业帮 编辑：搜搜考试网作业帮 分类：综合作业 时间：2024/06/23 11:53:30
急求MATLAB解常微分方程组数值解的解决方案
我用MATLAB7.4的ODE45解四元常微分方程组,程序如下：
function xdot = eqx(t,x);
m = 0.5;
alf = 0.15;
K=[0.0265,0.1,0.15];
Q=20;
a=0.8;
xdot = zeros(4,1);
xdot(1) = Q - alf * x(1) * x(2) + a*m*(x(1)+x(2)+x(3)+x(4));
xdot(2) = x(2) * (-m + K(2) * alf * x(1) - alf * x(3));
xdot(3) = x(3) * (-m + K(3) * alf * x(2) - alf * x(4));
xdot(4) = x(4) * (-m + K(4) * alf * x(3));
执行
>> [t,x] = ode45('eqx',[0 0.2] ,[0 125] ,[0 12.5e+3] ,[0 0.2e+3]);
后得到如下出错信息,
Error using ==> funfun\private\odearguments at 44
Correct syntax is ode45(EQX,tspan,y0,options).
Error in ==> ode45 at 173
[neq,tspan,ntspan,next,t0,tfinal,tdir,y0,f0,odeArgs,odeFcn,...
急求MATLAB解常微分方程组数值解的解决方案
调用方法：[t,x] = ode45(@eqx,[0 0.2] ,[0 1 1.1 1.25])x的初始值写在一个矩阵里．
另外你的函数定义时K少一维K(4),我给你加了一个：
function xdot = eqx(t,x);
m = 0.5;
alf = 0.15;
K=[0.0265,0.1,0.15,0.12];
Q=20;
a=0.8;
xdot = zeros(4,1);
xdot(1) = Q - alf * x(1) * x(2) + a*m*(x(1)+x(2)+x(3)+x(4));
xdot(2) = x(2) * (-m + K(2) * alf * x(1) - alf * x(3));
xdot(3) = x(3) * (-m + K(3) * alf * x(2) - alf * x(4));
xdot(4) = x(4) * (-m + K(4) * alf * x(3));
急求MATLAB解常微分方程组数值解的解决方案 matlab求微分方程组数值解 编写Matlab程序求常微分方程组的解 利用MATLAB求常微分方程组 的特解 matlab解常微分方程组 常微分方程组的求解,数值解怎么求出来, matlab数值解微分方程组问题 怎样用matlab解常微分的数值解 如何用matlab求解微分方程组的数值解 用MATLAB解常微分方程组 用Matlab能解 常微分方程组 吗?我是要计算其数值解 matlab求解常微分方程组数值解,总出现这个,急死了