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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 计算模拟 » 程序代码 » 求高版本的1stopt跑一下以下代码，谢谢！
31/1返回列表
查看: 763  |  回复: 2
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

iloveyou130

金虫 (小有名气)

[求助] 求高版本的1stopt跑一下以下代码，谢谢！ 已有1人参与

第一次发帖写程序，也不知道有没有写对，悬赏金币也不知是否合适。在以下程序中，变量y8没有在微分方程中出现过，由于y12~y14不可测量，所以没有相关数据，不知道这样可以拟合k的值吗？
---------------------------------------------------------------------------
Parameter k(17);
k(17)>=0;
Variable t,y1,y2,y3,y4,y5,y6,y7,y8,y9,y10,y11;
ODEFunction y1'=-2*k1*y2*y1-2*k3*y2*y1-8*k16*y2*y1-k17*y7*y1+k17/2.2*y6*y9;
y2'=-k1*y2*y1-k3*y2*y1-k5*y2*y6-k7*y2-k16*y2*y1;
y3'=k2*y5*y6+k4*y4*y6+k4*y4*y6+k5*y2*y6-k6*y3;
y4'=k3*y2*y6-k4*y4*y6;
y5'=k1*y2-k2*y5*y6;
y6'=-2*k1*y2*y1-k2*y5*y6-k4*y4*y6-4*k5*y2*y6+k9*y14+3*k15-y11+10*k16*y2*y1+k17*y7*y1-k17/2.2*y6*y9;
y7'=2*k1*y2*y1+k10*y14+8*k16*y2*y1-k17*y7*y1+k17/2.2*y6*y9;
y8'=k4*y4*y6+k11*y14;
y9'=k3*y2*y1+k12*y14+k17*y7*y1-k17/2.2*y6*y9;
y10'=k5*y2*y6+k13*y14;
y11'=k2*y5*y6+k3*y2*y1+k5*y2*y6+k14*y14-2*k15*y11;
y12'=k15*y11;
y13'=k8*y14;
y14'=k6*y3+k7*y2-k8*y14;
data;
10        5.41        0.19        9.56E-05        0.000966        0.000549        0.00165        0.00119        0.00134        0.00162        0        0.00169
20        5.45        0.17        2.62E-05        0.0023        0.00112        0.00302        0.00221        0.00174        0.00338        0.000283        0.00452
30        5.46        0.15        0.000786        0.00427        2.90E-03        0.00403        0.00317        0.00209        0.000633        0.00071        0.00527
45        5.49        0.12        0.00109        0.00461        0.0026        0.00998        0.00599        0.00644        0.00282        0.000824        0.0052
60        5.51        0.11        0.00208        0.00633        0.00434        0.0211        0.0053        0.0184        0.00841        0.00221        0.0078
80        5.44        0.1        0.00361        0.00783        0.00651        0.0136        0.00622        0.0103        0.00308        0.00153        0.00844
回复此楼

» 猜你喜欢
1楼 2020-06-16 20:36:32
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

独孤神宇

版主 (知名作家)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
iloveyou130: 金币+20, ★★★很有帮助 2020-06-18 11:17:03
均方差(RMSE): 0.0199436505934012
残差平方和(SSR): 0.021876205944542
相关系数(R): 0.406206827247245
相关系数之平方(R^2): 0.165003986502274
修正R平方(Adj. R^2): 1.15044908256948
确定系数(DC): -0.68967324594123
F统计(F-Statistic): -0.66404334401412

参数                  最佳估算
--------------------        -------------
k1        4.60883336079867E-8
k2        1.54247889887596
k3        5.1150136876246E-16
k4        0.000387115599728942
k5        3.30143387048788
k6        9.90173114206637E-12
k7        1.41860538976745
k8        19771.887540242
k9        0.341602464209552
k10        1125.69386161916
k11        2369.80489888995
k12        693.150581014912
k13        33.3865189378566
k14        1095.01407663084
k15        0.0176173700760016
k16        4.28454221172775E-13
k17        5.96905321164476E-13
y12 初值         0
y13 初值         0
y14 初值         0.000805185581681898

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

文件: 数据文件-1
No        t        目标 y1        计算 y1        目标 y2        计算 y2        目标 y3        计算 y3        目标 y4        计算 y4        目标 y5        计算 y5        目标 y6        计算 y6        目标 y7        计算 y7        目标 y8        计算 y8        目标 y9        计算 y9        目标 y10        计算 y10        目标 y11        计算 y11
1        20        5.45        5.40999999188834        0.17        0.16673223415736        2.62E-5        0.000288890181847264        0.0023        0.000965999875847417        0.00112        0.00054871967692699        0.00302        0.0054791674073478        0.00221        0.00254889247969801        0.00174        0.00420071647181485        0.00338        0.00245673858528659        0.000283        0.000233311486798891        0.00452        0.0031972088649323
2        30        5.46        5.40999998463729        0.15        0.14576533273367        0.000786        0.000629754879635674        0.00427        0.00096599962755348        0.0029        0.000548158657804571        0.00403        0.00867615151602474        0.00317        0.00372337399212577        0.00209        0.00667321386445097        0.000633        0.00317992589617602        0.00071        0.000608447276701977        0.00527        0.00466761215373619
3        45        5.49        5.40999997830098        0.12        0.127319431959793        0.00109        0.00105189712260829        0.00461        0.000965999276436763        0.0026        0.000547365929143034        0.00998        0.0114145598811571        0.00599        0.00474963372253263        0.00644        0.00883367755620465        0.00282        0.00381184560421971        0.000824        0.00106023275575494        0.0052        0.00607038597326002
4        60        5.51        5.40999997276868        0.11        0.111119483171172        0.00208        0.00151543628490384        0.00633        0.000965998835359388        0.00434        0.000546371494207135        0.0211        0.0138612901082428        0.0053        0.0056457324530104        0.0184        0.0107201287662758        0.00841        0.00436361864140109        0.00221        0.00154935297132375        0.0078        0.00738348099987025
5        80        5.44        5.40999996794202        0.1        0.0969107562795792        0.00361        0.00199525344001146        0.00783        0.000965998312451843        0.00651        0.000545194749301669        0.0136        0.016125635375396        0.00622        0.00642747209680412        0.0103        0.012365833385529        0.00308        0.00484497503960053        0.00153        0.00205117706028891        0.00844        0.00859746204094071
赞一下(1人) 回复此楼
数值计算
2楼2020-06-16 21:59:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

iloveyou130

金虫 (小有名气)
引用回帖:
2楼: Originally posted by 独孤神宇 at 2020-06-16 21:59:52
均方差(RMSE): 0.0199436505934012
残差平方和(SSR): 0.021876205944542
相关系数(R): 0.406206827247245
相关系数之平方(R^2): 0.165003986502274
修正R平方(Adj. R^2): 1.15044908256948
确定系数(DC): -0.6 ...

谢谢！
回复此楼
3楼2020-06-18 11:09:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 iloveyou130 的主题更新
31/1返回列表
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭