| 查看: 524 | 回复: 2 | ||
[求助]
求1stopt高版本软件帮忙运行程序 已有1人参与
|
|
那位大虾帮我用高版本的1stopt运行一下下面两段的程序：微分方程参数拟合： Variable X,S,P,T; Parameter r=[0,1],m=[0,1],a=[0,1],b=[0,1],rr=[0,1]; ODEFunction X'=r*X*(1-X/m)-(c*S)/(S+d); S'=a*X*(0.008343*P+1.258)-b*(67.39*T/(T+13.53))*a*X*(0.008343*P+1.258)-rr*X; Data; 0.024 0.0322 108.2 1 0.026 0.0284 196.3 2 0.048 0.0233 168.3 3 0.052 0.0235 281.4 4 0.087 0.0224 144.1 5 0.174 0.018 119.7 6 0.13 0.017 173.6 7 0.164 0.0167 187.1 8 0.1632 0.0166 214.9 9 0.1542 0.0152 121 10 0.1952 0.0142 180.7 11 0.2289 0.0132 181.8 12 0.2585 0.016 201.8 13 0.2888 0.0263 219.9 14 0.3522 0.0198 215.4 15 0.287 0.0145 157.8 16 0.2222 0.0137 200.5 17 0.2323 0.0109 141.2 18 0.2138 0.0113 151.6 19 0.1972 0.0128 225.1 20 0.2105 0.01635 286.6 21 0.1794 0.0151 246.6 22 0.1661 0.01215 136.6 23 0.1822 0.01215 74.4 24 0.1738 0.01185 124.6 25 0.1904 0.0112 271.2 26 0.1612 0.01085 137.9 27 0.1497 0.01045 153 28 0.1422 0.0111 137 29 0.1233 0.0127 185.8 30 0.1026 0.0145 178.6 31 0.1408 0.0146 215.4 32 0.1349 0.01335 157.8 33 0.1368 0.0157 200.5 34 0.1322 0.01 141.2 35 0.12015 0.0103 151.6 36 0.11775 0.0135 225.1 37 0.1083 0.01875 286.6 38 0.1014 0.02045 246.6 39 0.10965 0.0163 136.6 40 0.1035 0.0134 74.4 41 0.09215 0.01595 124.6 42 0.0895 0.0137 271.2 43 0.0924 0.01365 137.9 44 0.09105 0.0173 153 45 0.09765 0.0193 137 46 0.0952 0.0202 185.8 47 0.09495 0.0157 178.6 48 0.09 0.0086 74.4 49 0.09 0.0083 124.6 50 0.091 0.0109 271.2 51 0.093 0.0159 137.9 52 0.093 0.0148 153 53 0.095 0.0143 137 54 0.093 0.0177 185.8 55 0.096 0.0109 178.6 56 ========================================================= 第二段程序 Variable X,S,P,T; Parameter r=[0,1],b=[0,1],rr=[0,1]; ODEFunction X'=r*X*(1-X/0.38)-(c*S)/(S+d); S'=0.08*X*(0.008343*P+1.258)-b*(67.39*T/(T+13.53))*0.08*X*(0.008343*P+1.258)-rr*X; Data; 0.024 0.0322 108.2 1 0.026 0.0284 196.3 2 0.048 0.0233 168.3 3 0.052 0.0235 281.4 4 0.087 0.0224 144.1 5 0.174 0.018 119.7 6 0.13 0.017 173.6 7 0.164 0.0167 187.1 8 0.1632 0.0166 214.9 9 0.1542 0.0152 121 10 0.1952 0.0142 180.7 11 0.2289 0.0132 181.8 12 0.2585 0.016 201.8 13 0.2888 0.0263 219.9 14 0.3522 0.0198 215.4 15 0.287 0.0145 157.8 16 0.2222 0.0137 200.5 17 0.2323 0.0109 141.2 18 0.2138 0.0113 151.6 19 0.1972 0.0128 225.1 20 0.2105 0.01635 286.6 21 0.1794 0.0151 246.6 22 0.1661 0.01215 136.6 23 0.1822 0.01215 74.4 24 0.1738 0.01185 124.6 25 0.1904 0.0112 271.2 26 0.1612 0.01085 137.9 27 0.1497 0.01045 153 28 0.1422 0.0111 137 29 0.1233 0.0127 185.8 30 0.1026 0.0145 178.6 31 0.1408 0.0146 215.4 32 0.1349 0.01335 157.8 33 0.1368 0.0157 200.5 34 0.1322 0.01 141.2 35 0.12015 0.0103 151.6 36 0.11775 0.0135 225.1 37 0.1083 0.01875 286.6 38 0.1014 0.02045 246.6 39 0.10965 0.0163 136.6 40 0.1035 0.0134 74.4 41 0.09215 0.01595 124.6 42 0.0895 0.0137 271.2 43 0.0924 0.01365 137.9 44 0.09105 0.0173 153 45 0.09765 0.0193 137 46 0.0952 0.0202 185.8 47 0.09495 0.0157 178.6 48 0.09 0.0086 74.4 49 0.09 0.0083 124.6 50 0.091 0.0109 271.2 51 0.093 0.0159 137.9 52 0.093 0.0148 153 53 0.095 0.0143 137 54 0.093 0.0177 185.8 55 0.096 0.0109 178.6 56 |
回复此楼» 猜你喜欢
【复旦大学】二维材料方向招收2026年博士研究生1名
已经有0人回复
-
北京纳米能源与系统研究所 王中林院士/曹南颖研究员课题组2026级硕/博/博后招生
已经有10人回复
-
物理学I论文润色/翻译怎么收费?
已经有180人回复
-
荷兰Utrecht University超快太赫兹光谱王海教授课题招收2026 CSC博士生
已经有23人回复
-
反铁磁体中的磁性切换：两种不同的机制已成功可视化
已经有0人回复
-
求标准粉末衍射卡号 ICDD 01-076-1802
已经有0人回复
-
新西兰Robinson研究所招收全奖PhD
已经有0人回复
-
石墨烯转移--二氧化硅衬底石墨烯
已经有0人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 求高版本1stOpT跑代码
已经有8人回复
- 看了下1stOpt的介绍，感觉略显浮夸
已经有14人回复
- 1stopt运行报错“参数虽被定义，但是没有在函数表示中出现”，顺便求高手帮忙跑下程序
已经有3人回复
- 求大哥大姐帮跑一下1stopt~新手啊，版本低啊
已经有6人回复
- 谁能帮忙用1Stopt运行下程序呀？谢谢，我的1.5版本的运行不了啊
已经有6人回复
- 求大神跑下1stOpt，一个sharemodel，谢谢！
已经有7人回复
- 求高手给做一道题，用MATLAB或1stOpt，请把结果直接给出来，我的软件就不行，谢谢
已经有5人回复
- 求高版本1stopt跑一个程序
已经有7人回复
- 用1stOpt进行多模型多变量非线性拟合得不到结果
已经有4人回复
- 1stopt能进行隐函数非线性拟合吗？谁能帮拟合一下新型液体的粘度混合计算公式
已经有5人回复
- 我的MS低版本打不开高版本，怎么办啊，急求助！！！
已经有9人回复
- 求前辈帮忙用1stOpt你个一下数据，写论文急用，先谢过
已经有14人回复
- 求大哥大姐帮跑一下1stopt~新手啊，版本低啊
已经有9人回复
- 哪位朋友有注册版1stOpt软件？急需
已经有5人回复
2楼2014-11-06 09:20:57
dingd
铁杆木虫 (职业作家)
- 计算强帖: 4
- 应助: 1641 (讲师)
- 金币: 15037.3
- 散金: 101
- 红花: 234
- 帖子: 3410
- 在线: 1223.5小时
- 虫号: 291104
- 注册: 2006-10-28
【答案】应助回帖
感谢参与，应助指数 +1
|
微分方程拟合计算比较费时间，下面是第一段代码的参考结果（有时间的话结果还会更好）： 均方差(RMSE):0.017640728970135 残差平方和(SSE):0.0342314850457537 相关系数(R): 0.602461362176299 相关系数之平方(R^2): 0.362959692915321 确定系数(DC): -3.36028539491201 F统计(F-Statistic): 69.8177377311799 参数 最佳估算 -------------------- ------------- c 0.00106428654192994 d -0.0349011012047928 r 3.60972530286868E-10 m 9.33157757082492E-10 a 0.00332753093676446 b 0.0278157876334486 rr 0.0014446672737485 ====== 结果输出 ====== 文件: 数据文件-1 No t 目标 x 计算 x 目标 s 计算 s 1 2 0.026 0.036819251843059 0.0284 0.0323928918894434 2 3 0.048 0.0504793431206788 0.0233 0.0326048751860358 3 4 0.052 0.0655576197566643 0.0235 0.0329407245701933 4 5 0.087 0.0825052174306297 0.0224 0.0331619435163049 5 6 0.174 0.100795442428006 0.018 0.033346624951942 6 7 0.13 0.120787975446028 0.017 0.0335716045574381 7 8 0.164 0.143535218303242 0.0167 0.0337889386727208 8 9 0.1632 0.170114630181354 0.0166 0.0339994979630296 9 10 0.1542 0.19858412512844 0.0152 0.0340509079602557 10 11 0.1952 0.225043607501893 0.0142 0.0340911826461687 11 12 0.2289 0.247476562408349 0.0132 0.0340519379474219 12 13 0.2585 0.262343553713506 0.016 0.033930360088174 13 14 0.2888 0.268735755210552 0.0263 0.0337193453012847 14 15 0.3522 0.267836891537942 0.0198 0.0334143145239424 15 16 0.287 0.261864960294848 0.0145 0.0330303388838698 16 17 0.2222 0.25297900847369 0.0137 0.0325839453327286 17 18 0.2323 0.242727942646882 0.0109 0.0321090741112015 18 19 0.2138 0.232090832828099 0.0113 0.0316021888326738 19 20 0.1972 0.221521024683804 0.0128 0.0310278125687688 20 21 0.2105 0.211226796686975 0.01635 0.0303814698534273 21 22 0.1794 0.201384721932077 0.0151 0.029737451164964 22 23 0.1661 0.192168163806734 0.01215 0.0291740022818871 23 24 0.1822 0.183664080011589 0.01215 0.0286729929308619 24 25 0.1738 0.175829539643054 0.01185 0.0281343466862331 25 26 0.1904 0.168535998576585 0.0112 0.0274510459221761 26 27 0.1612 0.161724726074676 0.01085 0.0268766834354867 27 28 0.1497 0.155410096276067 0.01045 0.02630701093047 28 29 0.1422 0.149543711382441 0.0111 0.0257545687657379 29 30 0.1233 0.144079761827116 0.0127 0.0251582642253836 30 31 0.1026 0.138970741825484 0.0145 0.0245667266219661 31 32 0.1408 0.134183341186696 0.0146 0.0239425583475144 32 33 0.1349 0.12969106566469 0.01335 0.0233783820402628 33 34 0.1368 0.12547546471084 0.0157 0.0227831318167776 34 35 0.1322 0.121511828246697 0.01 0.0222496977090231 35 36 0.12015 0.1177872212671 0.0103 0.021719627760067 36 37 0.11775 0.114272485486526 0.0135 0.0211174136925277 37 38 0.1083 0.110936188795761 0.01875 0.020447330994152 38 39 0.1014 0.107762174781625 0.02045 0.0198169763772046 39 40 0.10965 0.104752119949238 0.0163 0.0193123707600452 40 41 0.1035 0.101909009228637 0.0134 0.0188928965660205 41 42 0.09215 0.0992215745590295 0.01595 0.0184365715356579 42 43 0.0895 0.0966637264522639 0.0137 0.017831130642767 43 44 0.0924 0.0942207813135252 0.01365 0.017346513857017 44 45 0.09105 0.0918949943467024 0.0173 0.0168708559987182 45 46 0.09765 0.0896772250290444 0.0193 0.0164149611221576 46 47 0.0952 0.0875586327235844 0.0202 0.0159191627761857 47 48 0.09495 0.0855293652404816 0.0157 0.0154300542654557 48 49 0.09 0.0835880867476915 0.0086 0.015048237961479 49 50 0.09 0.0817331351387703 0.0083 0.0146369814626851 50 51 0.091 0.0799502666950306 0.0109 0.0140859026699785 51 52 0.093 0.0782316922833754 0.0159 0.0136476948874128 52 53 0.093 0.0765801833523736 0.0148 0.013217663923481 53 54 0.095 0.0749912461179193 0.0143 0.0128058267597263 54 55 0.093 0.0734603614707211 0.0177 0.0123566752817365 55 56 0.096 0.0719821261000118 0.0109 0.0119135640050082 |
3楼2014-11-06 10:55:29