| 查看: 790 | 回复: 2 | ||
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 |
回复此楼» 猜你喜欢
［调剂信息］211智能人工感知方向国家青年特聘专家课题组招收调剂研究生
已经有0人回复
-
[调剂信息]211智能人工感知方向国家青年特聘专家课题组招收调剂研究生
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有80人回复
-
0702一志愿吉大B区求调剂 本科期间发表一篇Sci
已经有2人回复
-
070200求调剂，一志愿某211，288分
已经有18人回复
-
法国博士后职位
已经有0人回复
-
重庆交大26年硕士生招生拟调剂通知已出！欢迎加入机器视觉与3D光学成像课题组。
已经有0人回复
-
广州大学光电信息工程专业调剂，招收物理学专业学生
已经有1人回复
-
0702一志愿吉大B区求调剂 本科期间发表一篇Sci
已经有0人回复
独孤神宇
版主 (知名作家)
- 应助: 490 (硕士)
- 贵宾: 0.008
- 金币: 31016.3
- 散金: 802
- 红花: 122
- 沙发: 1
- 帖子: 5600
- 在线: 856.5小时
- 虫号: 3522474
- 注册: 2014-11-06
- 性别: GG
- 专业: 机械动力学
- 管辖: 计算模拟
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
iloveyou130: 金币+20, ★★★很有帮助 2020-06-18 11:17:03
感谢参与，应助指数 +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 |
2楼2020-06-16 21:59:52
iloveyou130
金虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 1189.4
- 散金: 20
- 帖子: 80
- 在线: 18.8小时
- 虫号: 2835371
- 注册: 2013-11-29
- 专业: 多相流热物理学
3楼2020-06-18 11:09:00
