| 查看: 897 | 回复: 6 | ||
hua99dushu铁虫 (小有名气)
|
[求助]
求帮忙运行一下下面的1stopt的程序，微分方程-参数估计 已有1人参与
|
|
Parameters a1=0.001[0,1],a2=0.001[0,1],a3=0.001[0,1],b1=0.001[0,1],b2=0.001[0,1],b3=0.001[0,1],c1=0.001[0,1],c2=0.001[0,1],c3=0.001[0,1],M=[0,1],C=[0,1],E=[0,1],N=[0,1]; 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.098,1.226,1.426,1.709,2.102,2.442,2.906,3.438,3.796,4.532,5.189]; y=[1.0000,1.017,1.081,1.219,1.470,1.838,2.064,2.247,2.397,2.449,2.742,3.209]; z=[1.0000,1.106,1.240,1.455,1.754,2.174,2.542,3.042,3.622,4.028,4.843,5.741]; 最好能给出误差图和相图呢。。 |
回复此楼» 猜你喜欢
山西大同大学物理学硕士研究生招收调剂生
已经有0人回复
-
[调剂信息]211智能人工感知方向国家青年特聘专家课题组招收调剂研究生
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有86人回复
-
［调剂信息］211智能人工感知方向国家青年特聘专家课题组招收调剂研究生
已经有0人回复
-
[调剂信息]211智能人工感知方向国家青年特聘专家课题组招收调剂研究生
已经有0人回复
-
0702一志愿吉大B区求调剂 本科期间发表一篇Sci
已经有2人回复
-
0702一志愿吉大B区求调剂 本科期间发表一篇Sci
已经有0人回复
-
070200求调剂，一志愿某211，288分
已经有18人回复
-
法国博士后职位
已经有0人回复
-
重庆交大26年硕士生招生拟调剂通知已出！欢迎加入机器视觉与3D光学成像课题组。
已经有0人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 求版主帮忙跑下面的1stopt的程序，参数估计。
已经有6人回复
- 用1stopt软件做微分方程组的参数估计
已经有13人回复
- 怎么解这个复数方程
已经有3人回复
- 求用高版本的1stOpt拟合方程
已经有8人回复
- 求大神用matlab给我处理一组数据，有人没有啊？
已经有15人回复
- matlab如何求解含有bessel方程中的参数
已经有5人回复
- 1stopt拟合问题求助
已经有7人回复
- 1stopt 拟合计算
已经有6人回复
- 用Matlab线性拟合求解两个参数，十万火急！
已经有9人回复
- 1stopt中使用sharedmodel拟合曲线出现问题，求指示！！
已经有6人回复
- 重金请教-如何用MAtlab对力学模型进行参数识别
已经有12人回复
- matlab 拟合反应动力学参数结果很差。大家帮忙看一下
已经有14人回复
- 跪求1stOpt2.5以上版本，或者请高人帮忙拟合包含复数的函数的参数
已经有8人回复
- 需要用matlab拟合数据的外行人急需帮助
已经有10人回复
hua99dushu
铁虫 (小有名气)
- 应助: 1 (幼儿园)
- 金币: 74.1
- 帖子: 161
- 在线: 25.5小时
- 虫号: 2078718
- 注册: 2012-10-22
- 性别: MM
- 专业: 控制论中的数学方法
2楼2014-03-22 11:38:21
hua99dushu
铁虫 (小有名气)
- 应助: 1 (幼儿园)
- 金币: 74.1
- 帖子: 161
- 在线: 25.5小时
- 虫号: 2078718
- 注册: 2012-10-22
- 性别: MM
- 专业: 控制论中的数学方法
3楼2014-03-22 11:39:41
hua99dushu
铁虫 (小有名气)
- 应助: 1 (幼儿园)
- 金币: 74.1
- 帖子: 161
- 在线: 25.5小时
- 虫号: 2078718
- 注册: 2012-10-22
- 性别: MM
- 专业: 控制论中的数学方法
★
fegg7502: 金币+1, 鼓励交流 2014-03-24 08:33:10
fegg7502: 金币+1, 鼓励交流 2014-03-24 08:33:10
|
上面的错误挺多的，程序更正确定如下： Parameters a1=0.001[0,1],a2=0.001[0,1],a3=0.001[0,1],b1=0.001[0,1],b2=0.001[0,1],b3=0.001[0,1], c1=0.001[0,1],c2=0.001[0,1],c3=0.001[0,1],M=[0,1],C=[0,1],E=[0,1],N=[0,1]; 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： t=[0,1,2,3,4,5,6,7,8,9,10,11]; x=[1.0000,1.098,1.226,1.426,1.709,2.102,2.442,2.906,3.438,3.796,4.532,5.189]; y=[1.0000,1.017,1.081,1.219,1.470,1.838,2.064,2.247,2.397,2.449,2.742,3.209]; z=[1.0000,1.106,1.240,1.455,1.754,2.174,2.542,3.042,3.622,4.028,4.843,5.741]; 希望能够得到估计的参数。。最好能有相关的图 |
4楼2014-03-22 14:24:35
hua99dushu
铁虫 (小有名气)
- 应助: 1 (幼儿园)
- 金币: 74.1
- 帖子: 161
- 在线: 25.5小时
- 虫号: 2078718
- 注册: 2012-10-22
- 性别: MM
- 专业: 控制论中的数学方法
|
请版主看此楼。。。。。快点帮忙运行一下吧。。 Parameters a1=0.001[0,1],a2=0.001[0,1],a3=0.001[0,1],b1=0.001[0,1],b2=0.001[0,1],b3=0.001[0,1], c1=0.001[0,1],c2=0.001[0,1],c3=0.001[0,1],M=[0,5],C=[0,10],E=[0,10],N=[0,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： t=[0,1,2,3,4,5,6,7,8,9,10,11]; 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]; |
5楼2014-03-23 20:18:44
dingd
铁杆木虫 (职业作家)
- 计算强帖: 4
- 应助: 1641 (讲师)
- 金币: 15037.3
- 散金: 101
- 红花: 234
- 帖子: 3410
- 在线: 1223.7小时
- 虫号: 291104
- 注册: 2006-10-28
【答案】应助回帖
★ ★
感谢参与，应助指数 +1
fegg7502: 金币+2, 应助指数+1, 3ks 2014-03-24 08:33:25
感谢参与，应助指数 +1
fegg7502: 金币+2, 应助指数+1, 3ks 2014-03-24 08:33:25
|
只大概算了下，供参考： 均方差(RMSE):0.227821606236476 残差平方和(SSE):1.71278858084954 相关系数(R): 0.968347715467601 相关系数之平方(R^2): 0.937697298051321 确定系数(DC): 0.810472079088207 F统计(F-Statistic): -6.59531364009876 参数 最佳估算 -------------------- ------------- a1 3.45198672383651E-11 a2 5.83733964128734E-8 a3 0.138472304076658 b1 1.44029340345734E-6 b2 0.0518147236689287 b3 0.0406931861755383 c1 3.70541984860404E-6 c2 1.37133289578244E-7 c3 8.59791920071737E-10 m 3.16205888404511 c 2.08529659127593 e 3.70210463181138 n 6.20898311572885E-5 ====== 结果输出 ====== 文件: 数据文件-1 No t 目标 x 计算 x 目标 y 计算 y 目标 z 计算 z 1 1 1.0777 1.14300503145253 1.0597 1.05729662499709 1.092 1.06851434515071 2 2 1.1443 1.29685305954089 1.0982 1.11588828762999 1.2033 1.1572932896264 3 3 1.3497 1.46465905303486 1.2035 1.17578665477138 1.4031 1.2709880836245 4 4 1.5881 1.65027065718425 1.3566 1.23690130994621 1.692 1.41560520068449 5 5 1.9814 1.85849139764935 1.461 1.29894832240681 1.9382 1.59905260863551 6 6 2.1905 2.09539536339523 1.4771 1.36129294495917 2.2159 1.8319321926032 7 7 2.5799 2.36877609710511 1.4862 1.42267433005819 2.6186 2.12871429124705 8 8 2.9763 2.68879681981474 1.6052 1.48071370354723 2.9489 2.50951876009061 9 9 3.2445 3.06894995110712 1.6094 1.53101337771447 3.1536 3.00288297139041 10 10 3.5765 3.52750422345532 1.8892 1.56545537275185 3.5979 3.65018117220115 11 11 3.9315 4.08974278597234 2.0334 1.56886778417041 4.0233 4.51289544111373 |
6楼2014-03-23 21:52:53
hua99dushu
铁虫 (小有名气)
- 应助: 1 (幼儿园)
- 金币: 74.1
- 帖子: 161
- 在线: 25.5小时
- 虫号: 2078718
- 注册: 2012-10-22
- 性别: MM
- 专业: 控制论中的数学方法
7楼2014-03-24 20:49:49
