| 查看: 1063 | 回复: 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]; 最好能给出误差图和相图呢。。 |
回复此楼» 猜你喜欢
澳门大学IAPME李海峰教授常年招收博士研究生
已经有0人回复
-
深圳大学应用技术学院招聘凝聚态物理博士后
已经有123人回复
-
物理学I论文润色/翻译怎么收费?
已经有96人回复
-
求助中医院临床中医中标的青年标书学习，谢谢
已经有5人回复
-
浙大物院表面物理课题组长期招收普博及直博学生
已经有259人回复
-
【征稿】Springer Nature热电专题 Discover Materials / Discover Chemistry
已经有47人回复
-
25届硕士申请，材料计算方向
已经有2人回复
-
课题组招1-2名硕士
已经有1人回复
-
机时越来越贵之后，我把一部分 DFT 流程搬到了 AI 计算平台上
已经有0人回复
-
斯坦福大学与Elsevier世界前2%学者新排名查询
已经有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
- 专业: 控制论中的数学方法
|
请版主看此楼。。。。。快点帮忙运行一下吧。。 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
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
