| 查看: 882 | 回复: 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]; 最好能给出误差图和相图呢。。 |
回复此楼» 猜你喜欢
帮我的英语口语老师找学生
已经有4人回复
-
计算三阶力常数用thirdorder扩胞时参数如何选择？临近原子数怎么取值
已经有7人回复
-
物理学I论文润色/翻译怎么收费?
已经有109人回复
-
2026年循环经济功能材料国际会议（ICFMCE 2026）
已经有0人回复
-
2026年第五届电气、电子与信息工程国际会议（ISEEIE 2026）
已经有0人回复
-
哈尔滨理工大学物理系招收物理学考研调剂
已经有4人回复
-
0702一志愿吉大B区求调剂
已经有5人回复
-
求调剂
已经有0人回复
-
0702一志愿吉大B区求调剂有论文
已经有0人回复
-
请问还有没有用Latex写文章的小伙伴们？
已经有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
- 专业: 控制论中的数学方法
3楼2014-03-22 11:39:41
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
- 专业: 控制论中的数学方法
★
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
20