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

24小时热门版块排行榜    

CyRhmU.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 计算模拟 » 求版主帮忙跑下面的1stopt的程序，参数估计。
南方科技大学公共卫生及应急管理学院2026级博士研究生招生报考通知（长期有效）
71/1返回列表
查看: 554  |  回复: 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的回帖

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
赞一下(1人) 回复此楼
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 的主题更新
71/1返回列表
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭