| 查看: 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论文润色/翻译怎么收费?
已经有58人回复
-
法国博士后职位
已经有0人回复
-
重庆交大26年硕士生招生拟调剂通知已出！欢迎加入机器视觉与3D光学成像课题组。
已经有0人回复
-
**
已经有1人回复
-
0702一志愿吉大B区求调剂 本科期间发表一篇Sci
已经有3人回复
-
基底STO,薄膜SRO，XRD里面的振荡，是laue震荡还是kiessig振荡? 怎么判断？
已经有0人回复
-
B区学生调剂-兰州交通大学材料科学与工程学院
已经有2人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 寻求方程组的解法，急急急！！
已经有7人回复
- origin曲线拟合求参数 出错
已经有7人回复
- 请教一个二元回归问题，求高手解决，谢谢。
已经有4人回复
- 求助，matlab/ 1stOpt 自定义微分方程拟合实验数据来求方程中的参数
已经有4人回复
- 两个自变量怎么拟合求参数值呢？
已经有7人回复
- 解带参数方程，参数在变
已经有6人回复
- 微分方程组参数拟合的问题（多参数动力系统参数估计）
已经有11人回复
- 求大哥大姐帮跑一下1stopt~新手啊，版本低啊
已经有9人回复
- 关于matlab微分、及非线性拟合问题
已经有5人回复
- MATLAB微分方程参数拟合问题，求大神
已经有7人回复
- 1stopt4阶非线性常微分方程，帮小弟看一下，跪求。
已经有7人回复
- matlab求解方程中的参数
已经有21人回复
- 形如1-(1-x)^a=by的方程怎么做拟合?
已经有11人回复
- 请问matlab拟合二元方程的问题
已经有3人回复
- 一阶微分方程求解，并优化参数值
已经有14人回复
- matlab微分方程组参数拟合，以周为单位求解，汇总后以年为单位进行数值比较
已经有4人回复
- 高金求助matlab解微分方程组
已经有12人回复
- 求回归方程参数估计的最优化算法matlab代码
已经有9人回复
- 【紧急求助】 带有微分平方项的ODE求解
已经有11人回复
- 重金加人民币求拟合wilson方程和NRTL方程参数的程序
已经有17人回复
- 【求助】求助一个一阶常微分方程的初值问题
已经有16人回复
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
- 专业: 控制论中的数学方法
2楼2014-02-10 13:54:34
hua99dushu
铁虫 (小有名气)
- 应助: 1 (幼儿园)
- 金币: 74.1
- 帖子: 161
- 在线: 25.5小时
- 虫号: 2078718
- 注册: 2012-10-22
- 性别: MM
- 专业: 控制论中的数学方法
4楼2014-02-11 10:18:42
hua99dushu
铁虫 (小有名气)
- 应助: 1 (幼儿园)
- 金币: 74.1
- 帖子: 161
- 在线: 25.5小时
- 虫号: 2078718
- 注册: 2012-10-22
- 性别: MM
- 专业: 控制论中的数学方法
★
fegg7502: 金币+1, 鼓励交流 2014-02-17 12:03:30
fegg7502: 金币+1, 鼓励交流 2014-02-17 12:03:30
|
Title "test 1"; Parameters a1[0,],a2[0,],a3[0,],b1[0,],b2[0,],b3[0,],c1[0,],c2[0,],c3[0,],M[8000,],C[y,],E[z,],N[300,]; 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; |
5楼2014-02-11 10:31:18
5