应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

CyRhmU.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 计算模拟 » 求帮忙运行一下下面的1stopt的程序，微分方程-参数估计
51/1返回列表
查看: 782  |  回复: 6
只看楼主@他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

hua99dushu

铁虫 (小有名气)

[求助] 求帮忙运行一下下面的1stopt的程序，微分方程-参数估计已有1人参与

Parameters a1=0.001[0,1],a2=0.001[0,1],a3=0.001[0,1],b1=0.001[0,1],b2=0.001[0,1],b3=0.001[0,1],c1=0.001[0,1],c2=0.001[0,1],c3=0.001[0,1],M=[0,1],C=[0,1],E=[0,1],N=[0,1];
Sharedmodel;
Variable t, x,y,z;
Function x'=a1*x*(y/M-1)-a2*y+a3*z;
y'=-b1*x+b2*y*(1-y/C)+b3*z*(1-z/E);
z'=c1*x*(x/N-1)-c2*y-c3*z;
RowData：
x=[1.0000,1.098,1.226,1.426,1.709,2.102,2.442,2.906,3.438,3.796,4.532,5.189];
y=[1.0000,1.017,1.081,1.219,1.470,1.838,2.064,2.247,2.397,2.449,2.742,3.209];
z=[1.0000,1.106,1.240,1.455,1.754,2.174,2.542,3.042,3.622,4.028,4.843,5.741];
最好能给出误差图和相图呢。。
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼2014-03-22 11:36:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hua99dushu

铁虫 (小有名气)
★
fegg7502: 金币+1, 鼓励交流 2014-03-24 08:33:10
上面的错误挺多的，程序更正确定如下：
Parameters a1=0.001[0,1],a2=0.001[0,1],a3=0.001[0,1],b1=0.001[0,1],b2=0.001[0,1],b3=0.001[0,1],
c1=0.001[0,1],c2=0.001[0,1],c3=0.001[0,1],M=[0,1],C=[0,1],E=[0,1],N=[0,1];
Sharedmodel;
Variable t,x,y,z;
Function x'=a1*x*(y/M-1)-a2*y+a3*z;
y'=-b1*x+b2*y*(1-y/C)+b3*z*(1-z/E);
z'=c1*x*(x/N-1)-c2*y-c3*z;
RowData：
t=[0,1,2,3,4,5,6,7,8,9,10,11];
x=[1.0000,1.098,1.226,1.426,1.709,2.102,2.442,2.906,3.438,3.796,4.532,5.189];
y=[1.0000,1.017,1.081,1.219,1.470,1.838,2.064,2.247,2.397,2.449,2.742,3.209];
z=[1.0000,1.106,1.240,1.455,1.754,2.174,2.542,3.042,3.622,4.028,4.843,5.741];
希望能够得到估计的参数。。最好能有相关的图
赞一下 回复此楼
4楼2014-03-22 14:24:35
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 7 个回答

hua99dushu

铁虫 (小有名气)
★
fegg7502: 金币+1, 3ks 2014-03-24 08:32:54
fegg7502: 鼓励交流 2014-03-24 08:33:13
a1=0.001[0,1]这种格式应该是表示初值是0.001，范围是0~1之间哟
我也是从教程里看到的，如果不对，请帮忙改正呢。
赞一下 回复此楼
2楼2014-03-22 11:38:21
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hua99dushu

铁虫 (小有名气)
额，忘了t了
t=[0,1,2,3,4,5,6,7,8,9,10,11];
回复此楼
3楼2014-03-22 11:39:41
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hua99dushu

铁虫 (小有名气)
请版主看此楼。。。。。快点帮忙运行一下吧。。
Parameters a1=0.001[0,1],a2=0.001[0,1],a3=0.001[0,1],b1=0.001[0,1],b2=0.001[0,1],b3=0.001[0,1],
c1=0.001[0,1],c2=0.001[0,1],c3=0.001[0,1],M=[0,5],C=[0,10],E=[0,10],N=[0,5];
Sharedmodel;
Variable t,x,y,z;
Function x'=a1*x*(y/M-1)-a2*y+a3*z;
y'=-b1*x+b2*y*(1-y/C)+b3*z*(1-z/E);
z'=c1*x*(x/N-1)-c2*y-c3*z;
RowData：
t=[0,1,2,3,4,5,6,7,8,9,10,11];
x=[1.0000,1.0777,1.1443,1.3497,1.5881,1.9814,2.1905,2.5799,2.9763,3.2445,3.5765,3.9315];
y=[1.0000,1.0597,1.0982,1.2035,1.3566,1.4610,1.4771,1.4862,1.6052,1.6094,1.8892,2.0334];
z=[1.0000,1.0920,1.2033,1.4031,1.6920,1.9382,2.2159,2.6186,2.9489,3.1536,3.5979,4.0233];
赞一下 回复此楼
5楼2014-03-23 20:18:44
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 7 个回答
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭