| 查看: 2948 | 回复: 13 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
hua99dushu铁虫 (小有名气)
|
[求助]
用1stopt软件做微分方程组的参数估计 已有1人参与
|
||
|
近一个月里才知道1stopt这个软件，可以做微分方程组的参数估计问题。 之前一直想用matlab的神经网络来做，由于水平有限，也做不出来。 网上下载的1stopt软件大多是无法写代码或者未注册，根据一些实例写了我的代码，可能有错误。 请帮忙改正后，跑一下软件。 Title "test 1"; Parameters a1,a2,a3,b1,b2,b3,c1,c2,c3,M,C,E,N[0,]; 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; 0,1,2,3,4,5,6,7,8,9,10; 265,291,325,378,453,557,647,770,911,1006,1201; 7217,7344,7800,8796,10607,13266,14893,16219,17303,17677,19786; 8554,9457,10607,12443,15004,18599,21742,26018,30981,34457,41425; |
回复此楼» 猜你喜欢
0702一志愿吉大B区求调剂 本科期间发表一篇Sci
已经有2人回复
-
070200求调剂，一志愿某211，288分
已经有18人回复
-
物理学I论文润色/翻译怎么收费?
已经有254人回复
-
法国博士后职位
已经有0人回复
-
重庆交大26年硕士生招生拟调剂通知已出！欢迎加入机器视觉与3D光学成像课题组。
已经有0人回复
-
**
已经有1人回复
-
0702一志愿吉大B区求调剂 本科期间发表一篇Sci
已经有3人回复
-
基底STO,薄膜SRO，XRD里面的振荡，是laue震荡还是kiessig振荡? 怎么判断？
已经有0人回复
-
B区学生调剂-兰州交通大学材料科学与工程学院
已经有2人回复
hua99dushu
铁虫 (小有名气)
- 应助: 1 (幼儿园)
- 金币: 74.1
- 帖子: 161
- 在线: 25.5小时
- 虫号: 2078718
- 注册: 2012-10-22
- 性别: MM
- 专业: 控制论中的数学方法
送红花一朵 |
之前你帮忙做出的结果已经很好了，但有个别数据太小，所以我做了替代数据，请帮忙再跑一下，万分感谢。。 下面是程序： 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]; |
11楼2014-02-18 18:48:11
hua99dushu
铁虫 (小有名气)
- 应助: 1 (幼儿园)
- 金币: 74.1
- 帖子: 161
- 在线: 25.5小时
- 虫号: 2078718
- 注册: 2012-10-22
- 性别: MM
- 专业: 控制论中的数学方法
2楼2014-02-10 13:54:34
dingd
铁杆木虫 (职业作家)
- 计算强帖: 4
- 应助: 1641 (讲师)
- 金币: 15037.3
- 散金: 101
- 红花: 234
- 帖子: 3410
- 在线: 1223.7小时
- 虫号: 291104
- 注册: 2006-10-28
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
hua99dushu: 金币+30, ★★★很有帮助, 辛苦了，能否帮忙再做一下？付报酬。 2014-02-13 12:37:55
fegg7502: 金币+1, 鼓励交流 2014-02-17 12:03:19
感谢参与，应助指数 +1
hua99dushu: 金币+30, ★★★很有帮助, 辛苦了，能否帮忙再做一下？付报酬。 2014-02-13 12:37:55
fegg7502: 金币+1, 鼓励交流 2014-02-17 12:03:19
|
感觉你的公式有些问题，似乎可以简化，下面结果参考下： 均方差(RMSE):666.824975358025 残差平方和(SSE):13339666.4328369 相关系数(R): 0.99281958012936 相关系数之平方(R^2): 0.985690718688239 确定系数(DC): 0.752302595710559 F统计(F-Statistic): -0.355256071121647 参数 最佳估算 -------------------- ------------- a1 0.063244391232729 a2 26.6537700533949 a3 425.865501005945 b1 3256273550.99063 b2 1385.92982355728 b3 190853579.985187 c1 1967.65266459696 c2 7557.60987652871 c3 174107.388897955 m 24291.4733546288 c 0.10182112015089 e 1.27923766610631E22 n 0.0972369125973638 文件: 数据文件-1 No t 目标 x 计算 x 目标 y 计算 y 目标 z 计算 z 1 1 291 295.919384351006 7344 7409.69953342258 9457 8964.92474439791 2 2 325 329.82579568465 7800 8490.90903548783 10607 10769.8722203566 3 3 378 365.520171155061 8796 9774.32316004032 12443 13050.8936628585 4 4 453 403.287476232365 10607 11128.0651594216 15004 15713.5200315771 5 5 557 443.400579346621 13266 12544.2720271355 18599 18789.0601075391 6 6 647 486.105229396296 14893 14028.6589510091 21742 22331.0481756248 7 7 770 531.633932176283 16219 15588.082085985 26018 26401.9444748202 8 8 911 580.212516031404 17303 17228.8956295004 30981 31071.3625784165 9 9 1006 632.066439410749 17677 18957.0942374736 34457 36416.4169412114 10 10 1201 687.422431267362 19786 20778.230707146 41425 42522.166618056 |
3楼2014-02-10 22:03:32
hua99dushu
铁虫 (小有名气)
- 应助: 1 (幼儿园)
- 金币: 74.1
- 帖子: 161
- 在线: 25.5小时
- 虫号: 2078718
- 注册: 2012-10-22
- 性别: MM
- 专业: 控制论中的数学方法
4楼2014-02-11 10:18:42
