24小时热门版块排行榜    

查看: 778  |  回复: 6

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];

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

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

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的回帖

dingd

铁杆木虫 (职业作家)

参考下:

均方差(RMSE):0.0743428646356062
残差平方和(SSE):0.182386430233526
相关系数(R): 0.990490604242565
相关系数之平方(R^2): 0.981071637092801
确定系数(DC): 0.97907392333002
F统计(F-Statistic): -17.029556181244

参数                  最佳估算
--------------------        -------------
a1        1.30190551784766
a2        4.06003441093189
a3        0.869490127930172
b1        9.42946312355037
b2        1.57560369388824
b3        5.9857438663248
c1        0.0249019284686679
c2        13.4132204832955
c3        9.74396616982886E-6
m        332.144120934994
c        22.8549320352576
e        45695.8913468096
n        1.78908285577578

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

文件: 数据文件-1
No        t        目标 x        计算 x        目标 y        计算 y        目标 z        计算 z
1        1        1.0777        1.13953501723234        1.0597        1.05984607578657        1.092        1.13757057269921
2        2        1.1443        1.30021884916371        1.0982        1.12347532732008        1.2033        1.29904743776794
3        3        1.3497        1.48448399130672        1.2035        1.19141949444282        1.4031        1.48688285600025
4        4        1.5881        1.69476430026936        1.3566        1.26432646601749        1.692        1.70323733352619
5        5        1.9814        1.93339002993677        1.461        1.34297964122961        1.9382        1.94974242279121
6        6        2.1905        2.20245049650347        1.4771        1.42831743928452        2.2159        2.22719556655783
7        7        2.5799        2.50362263043615        1.4862        1.52145197784229        2.6186        2.5351804050362
8        8        2.9763        2.8379668528597        1.6052        1.62368427128187        2.9489        2.87161074307044
9        9        3.2445        3.20569653533687        1.6094        1.73651368452624        3.1536        3.23220275868481
10        10        3.5765        3.60593368625635        1.8892        1.8616395305816        3.5979        3.60988868394394
11        11        3.9315        4.03647149994331        2.0334        2.00094968015341        4.0233        3.99419780445337
5楼2014-02-24 11:52:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hua99dushu

铁虫 (小有名气)

引用回帖:
5楼: Originally posted by dingd at 2014-02-24 11:52:03
参考下:

均方差(RMSE):0.0743428646356062
残差平方和(SSE):0.182386430233526
相关系数(R): 0.990490604242565
相关系数之平方(R^2): 0.981071637092801
确定系数(DC): 0.97907392333002
F统计(F-Statist ...

真的很感谢你如此迅速的回复,多次认真的解答我的问题。
根据实际的条件,已经确定了一些参数。
[/CODE]
Parameters a1=[0,],a2=[0,],a3=[0,],b1=[0,],b2=[0,],b3=[0,],c1=[0,],c2=[0,],c3=[0,];
Variable t,x,y,z;
ODEFunction  x'=a1*x*(y/33.214-1)-a2*y+a3*z;
y'=-b1*x+b2*y*(1-y/22.854)+b3*z*(1-z/45.695);
z'=c1*x*(x/1.789-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];
求帮忙再跑一下上面的代码,还有就是是采用的什么算法啊?
6楼2014-02-24 14:46:18
已阅   回复此楼   关注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的回帖
相关版块跳转 我要订阅楼主 hua99dushu 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 只有每年这种时候来逛逛小木虫 +8 yaoyewhu2008 2026-08-20 8/400 2026-08-20 08:46 by since—2010
[基金申请] 今天基金会出结果吗?20260819 +14 kkkl_v 2026-08-19 15/750 2026-08-20 08:21 by sunzitan
[基金申请] 欢迎发来filecode的Mz6后的代码验证其规律 +55 医学老男孩 2026-08-13 123/6150 2026-08-19 23:44 by stan_et
[基金申请] filecode,4个jtjc了 +11 ziyangfang 2026-08-19 13/650 2026-08-19 21:28 by aasahr
[教师之家] 为什么余额宝的年化利率越来越低?主要原因有哪些? +5 瞬息宇宙 2026-08-15 5/250 2026-08-19 20:23 by super2002521
[基金申请] 今天放榜没戏了吧 +4 yuleib84 2026-08-19 5/250 2026-08-19 19:52 by hhs666
[基金申请] 今天放榜吗? +14 布布和一二 2026-08-19 15/750 2026-08-19 18:07 by gltch
[基金申请] filecode=后面第一个是大写字母 +10 wangze12014 2026-08-14 12/600 2026-08-19 16:56 by 苦难博士
[基金申请] 2027广东省杰青 +4 奶牛小黑 2026-08-15 10/500 2026-08-19 11:02 by wanfengnew
[基金申请] 朋友圈看到的 +6 wangzilk 2026-08-18 8/400 2026-08-19 10:55 by Haru815
[基金申请] 快农历七夕节了,轻松一下,男人悄悄话,女施主请不要进来。 +6 Tide man 2026-08-14 7/350 2026-08-19 09:56 by ZJTJZ
[基金申请] 什么时候开奖? +6 CrisMessi 2026-08-18 6/300 2026-08-19 08:06 by Equinoxhua
[基金申请] 重要消息,中午系统在维护 +10 yuleib84 2026-08-18 11/550 2026-08-19 07:30 by 羽毛枫f
[基金申请] 明天放榜? +5 Shxjjxjkx 2026-08-18 5/250 2026-08-18 18:14 by -大大大大大-
[论文投稿] 投稿咨询 +4 wwm09 2026-08-17 6/300 2026-08-18 15:36 by wwm09
[基金申请] 哪位老哥知道今年的国自然具体哪一天放榜? +13 Ldrop2023 2026-08-13 16/800 2026-08-18 12:25 by 淀粉搬运工
[基金申请] 93BebMhtakh前后11位开头都是大写 +4 且听虎啸 2026-08-17 5/250 2026-08-18 00:49 by 蔡棒棒菂
[基金申请] 感觉是下周放榜了 +6 angus9576 2026-08-17 11/550 2026-08-17 23:57 by angus9576
[精细化工] 招聘 金属平磨液,抛光液研发工程师 +3 小天0311 2026-08-14 3/150 2026-08-16 07:31 by H9PLUS
[基金申请] 各位道友,我要去昆明玩几天,回来见。 +7 Tide man 2026-08-14 8/400 2026-08-15 01:11 by arzu_hma
信息提示
请填处理意见