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

24小时热门版块排行榜    

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

hua99dushu

铁虫 (小有名气)

[求助] 求版主帮忙跑下面的1stopt的程序，参数估计。 已有1人参与

[/code]Title "test 1";
Parameters a1[0,],a2[0,],a3[0,],b1[0,],b2[0,],b3[0,],c1[0,],c2[0,],c3[0,],M[1.5,],C[y,],E[z,],N[1.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：
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];

感谢小木虫论坛，感谢斑竹们的耐心帮忙。
回复此楼

» 猜你喜欢

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

hua99dushu

铁虫 (小有名气)
有一行错了，正确的是：ODEFunction  x'=a1*x*(y/13.214-1)-a2*y+a3*z;
y/M，m=13.214.
赞一下 回复此楼
7楼2014-02-24 15:15:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 7 个回答

dingd

铁杆木虫 (职业作家)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
hua99dushu: 金币+30, ★★★很有帮助, 麻烦斑竹啦，版主好人。。 2014-02-24 10:06:49
参考：
CODE:
Parameters a1=[0,],a2=[0,],a3=[0,],b1=[0,],b2=[0,],b3=[0,],c1=[0,],c2=[0,],c3=[0,],M=[0,],C=[0,],E=[0,],N=[0,];
Variable t,x,y,z;
ODEFunction  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;
Data;
t=[0:1:12];
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];

均方差(RMSE):0.054186898245058
残差平方和(SSE):0.096895258066869
相关系数(R): 0.990468602959854
相关系数之平方(R^2): 0.981028053449245
确定系数(DC): 0.981119298577706
F统计(F-Statistic): -68.497450670183

参数                  最佳估算
--------------------        -------------
a1        0.934922929556056
a2        43.1779982925144
a3        0.634247035099968
b1        0.885424066820235
b2        1.72323467449167E-14
b3        1.64537787581792
c1        7.17385064510645E-11
c2        47.0215125361918
c3        0.0444107007235367
m        0.4861243914792
c        0.00874309356262975
e        201.823147971463
n        5.99589073558176

====== 结果输出 ======

文件: 数据文件-1
No        t        目标 x        计算 x        目标 y        计算 y        目标 z        计算 z
1        1        1.0777        1.03073467525285        1.0597        1.06691160587989        1.092        1.06234022501214
2        2        1.1443        1.18178677578727        1.0982        1.1387423298674        1.2033        1.22168768708318
3        3        1.3497        1.39065370562373        1.2035        1.21472123201098        1.4031        1.43111557008028
4        4        1.5881        1.63651472709221        1.3566        1.29454185100693        1.692        1.67502754880883
5        5        1.9814        1.91230103994829        1.461        1.37802854578025        1.9382        1.94813144944932
6        6        2.1905        2.21483623814851        1.4771        1.46501620348821        2.2159        2.24783814736362
7        7        2.5799        2.54118886530888        1.4862        1.5553142562366        2.6186        2.57123967133884
8        8        2.9763        2.88710686907347        1.6052        1.64870518739119        2.9489        2.91357690174086
9        9        3.2445        3.2461270702916        1.6094        1.74495815962043        3.1536        3.26712262655401
10        10        3.5765        3.60891777176678        1.8892        1.84385184177429        3.5979        3.62014395661707
11        11        3.9315        3.96272367218441        2.0334        1.9451949237521        4.0233        3.95587454727496
赞一下 回复此楼
2楼2014-02-20 20:40:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hua99dushu

铁虫 (小有名气)
引用回帖:
2楼: Originally posted by dingd at 2014-02-20 20:40:16
参考：

Parameters a1=,a2=,a3=,b1=,b2=,b3=,c1=,c2=,c3=,M=,C=,E=,N=;
Variable t,x,y,z;
ODEFunction  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;
Dat ...

有几个参数是这样设置的：
M[1.5,],C[1,],E[1,],N[1.5,]。
求帮着再跑一下啊。。
赞一下 回复此楼
3楼2014-02-21 14:54:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hua99dushu

铁虫 (小有名气)
引用回帖:
2楼: Originally posted by dingd at 2014-02-20 20:40:16
参考：

Parameters a1=,a2=,a3=,b1=,b2=,b3=,c1=,c2=,c3=,M=,C=,E=,N=;
Variable t,x,y,z;
ODEFunction  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;
Dat ...

之前程序里参数的设置我没有讲清楚。
是这样的：M[1.5,],C[1,],E[1,],N[1.5,].
帮忙再来一遍。。急着对比结果啊
赞一下 回复此楼
4楼2014-02-24 10:09:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 7 个回答
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[论文投稿] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +3 QTy3jDtz1uLt 2026-08-21 8/400 2026-08-22 23:51 by OEbVnUOu01ol
[硕博家园] 售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急 +4 QTy3jDtz1uLt 2026-08-21 10/500 2026-08-22 23:42 by OEbVnUOu01ol
[论文投稿] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +3 2JOx3r2CYEgw 2026-08-22 3/150 2026-08-22 22:03 by OEbVnUOu01ol
[基金申请] 93BebMhtakh前后11位开头都是大写 +7 且听虎啸 2026-08-17 8/400 2026-08-22 21:55 by 医学老男孩
[硕博家园] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +3 2JOx3r2CYEgw 2026-08-22 5/250 2026-08-22 21:42 by OEbVnUOu01ol
[教师之家] 跳槽后在研项目怎么办？ +3 简单化xn 2026-08-22 7/350 2026-08-22 19:46 by 简单化xn
[基金申请] 什么时候开奖？ +7 CrisMessi 2026-08-18 8/400 2026-08-22 18:29 by 淀粉搬运工
[基金申请] 只有每年这种时候来逛逛小木虫 +24 yaoyewhu2008 2026-08-20 26/1300 2026-08-22 17:43 by kammury
[基金申请] filecode，4个jtjc了 +13 ziyangfang 2026-08-19 16/800 2026-08-22 17:08 by WH3796
[基金申请] 放榜前的不淡定 40+4 snowwithsea 2026-08-19 14/700 2026-08-21 23:51 by cratir
[基金申请] 时间戳又变了 +13 wuchongjun 2026-08-20 19/950 2026-08-21 17:21 by 紫杉醇
[基金申请] 让我中一个面上吧！ +11 大萍1987 2026-08-20 13/650 2026-08-21 14:48 by 20120902066
[论文投稿] 投稿咨询 +5 wwm09 2026-08-17 7/350 2026-08-21 10:11 by 期刊论文帮手
[基金申请] 今天放榜没戏了吧 +9 yuleib84 2026-08-19 11/550 2026-08-21 10:06 by gltch
[基金申请] 估计是周四 +3 archvillain 2026-08-18 3/150 2026-08-21 01:48 by jnhyjjm
[基金申请] 今天系统多次维护，明天很可能放榜！ +10 zju2000 2026-08-16 11/550 2026-08-20 20:02 by cl479861084
[基金申请] 时间戳变了，能看出什么问题？ +18 基诺咪客 2026-08-17 23/1150 2026-08-20 17:19 by Godzela
[基金申请] 今天放榜吗？ +14 布布和一二 2026-08-19 15/750 2026-08-19 18:07 by gltch
[基金申请] 朋友圈看到的 +6 wangzilk 2026-08-18 8/400 2026-08-19 10:55 by Haru815
[基金申请] 明天放榜？ +5 Shxjjxjkx 2026-08-18 5/250 2026-08-18 18:14 by -大大大大大-
信息提示
关闭
请填处理意见
关闭