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

我用MATLAB7.4的ODE45解四元常微分方程组,程序如下:functionxdot=eqx(t,x);m=0.5;alf=0.15;K=[0.0265,0.1,0.... 我用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, ...
展开
日向淳正
2009-02-08 · TA获得超过2811个赞
知道大有可为答主
回答量:2770
采纳率:0%
帮助的人:1403万
展开全部
调用方法:[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));
本回答被提问者采纳
已赞过 已踩过<
你对这个回答的评价是?
评论 收起
富港检测技术(东莞)有限公司_
2024-04-02 广告
正弦振动多用于找出产品设计或包装设计的脆弱点。看在哪一个具体频率点响应最大(共振点);正弦振动在任一瞬间只包含一种频率的振动,而随机振动在任一瞬间包含频谱范围内的各种频率的振动。由于随机振动包含频谱内所有的频率,所以样品上的共振点会同时激发... 点击进入详情页
本回答由富港检测技术(东莞)有限公司_提供
greatdju
2009-02-07 · TA获得超过2.2万个赞
知道大有可为答主
回答量:2486
采纳率:50%
帮助的人:4187万
展开全部
你那个函数里K(4)是多少啊?没有定义
已赞过 已踩过<
你对这个回答的评价是?
评论 收起
收起 1条折叠回答
推荐律师服务: 若未解决您的问题,请您详细描述您的问题,通过百度律临进行免费专业咨询

为你推荐:

下载百度知道APP,抢鲜体验
使用百度知道APP,立即抢鲜体验。你的手机镜头里或许有别人想知道的答案。
扫描二维码下载
×

类别

我们会通过消息、邮箱等方式尽快将举报结果通知您。

说明

0/200

提交
取消

辅 助

模 式