| 查看: 741 | 回复: 2 | |||||
[交流]
求助1stopt解决简单地常微分方程参数拟合！有代码！跪谢！ 已有1人参与
|
|
Variables x,y; Parameters a; ODEFunction y'=(pi-y)/30/Ln(a^-1)*sin(y)*(2-sin(pi/2+y))*(cos(y)+0.217); data; 0 0.375594855 30 1.018748684 60 1.484227996 90 1.667487567 120 1.745852851 150 1.760339084 180 1.776919711 240 1.789486082 求问哪里错了？？结果总是显示： expression expected but ） found |
回复此楼» 猜你喜欢
070200求调剂，一志愿某211，288分
已经有20人回复
-
法国博士后职位
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有270人回复
-
重庆交大26年硕士生招生拟调剂通知已出！欢迎加入机器视觉与3D光学成像课题组。
已经有0人回复
-
**
已经有1人回复
-
0702一志愿吉大B区求调剂 本科期间发表一篇Sci
已经有3人回复
-
基底STO,薄膜SRO，XRD里面的振荡，是laue震荡还是kiessig振荡? 怎么判断？
已经有2人回复
-
B区学生调剂-兰州交通大学材料科学与工程学院
已经有8人回复
-
山西大同大学物理学专业还有调剂名额，欢迎调剂！
已经有9人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 求助Matlab大神帮忙求解复杂的参数方程组
已经有3人回复
- 求助帮忙用1stopt解一下这个方程组，方程组见附件，已经参考文件编写成1stopt的形式
已经有3人回复
- 1stOpt拟合代码的问题求助
已经有5人回复
- 进行线性拟合时取的试验点越多，拟合出的方程的方程越准确
已经有6人回复
- 使用MATLAB拟合动力学参数问题
已经有18人回复
- 看了下1stOpt的介绍，感觉略显浮夸
已经有14人回复
- 求助 1stopt低版本 方程的拟合（含若干未知参数）
已经有10人回复
- 求大哥大姐帮跑一下1stopt~新手啊，版本低啊
已经有6人回复
- 1stopt自变量超多时是否有简单语法？AND 自变量多 数据量大的话1stopt都会算的特慢么
已经有7人回复
- 1stopt数据拟合求助
已经有8人回复
- 微分方程组参数拟合的问题（多参数动力系统参数估计）
已经有11人回复
- 大家在使用1stOpt进行数据拟合的时候有没有遇到过这种情况？
已经有8人回复
- 1stopt4阶非线性常微分方程，帮小弟看一下，跪求。
已经有7人回复
- 一阶微分方程求解，并优化参数值
已经有14人回复
- matlab 拟合反应动力学参数结果很差。大家帮忙看一下
已经有14人回复
- matlab-常微分方程参数估计
已经有12人回复
月只蓝
主管区长 (职业作家)
-
专家经验: +1059 - 计算强帖: 8
- 应助: 1712 (讲师)
- 贵宾: 8.888
- 金币: 68129.7
- 散金: 1938
- 红花: 443
- 沙发: 4
- 帖子: 4373
- 在线: 3291.4小时
- 虫号: 1122189
- 注册: 2010-10-14
- 专业: 宇宙学
- 管辖: 计算模拟区
★ ★
小木虫: 金币+0.5, 给个红包，谢谢回帖
fegg7502: 金币+1, 鼓励交流 2015-06-29 08:37:19
小木虫: 金币+0.5, 给个红包，谢谢回帖
fegg7502: 金币+1, 鼓励交流 2015-06-29 08:37:19
|
ODEFunction y'=(pi-y)/30/Ln(a^-1)*sin(y)*(2-sin(pi/2+y))*(cos(y)+0.217); 改成 ODEFunction y'=(pi-y)/30/Ln(a^(-1))*sin(y)*(2-sin(pi/2+y))*(cos(y)+0.217); 计算结果： Root of Mean Square Error (RMSE): 0.0158516811713416 Residual Sum of Squares (SSE): 0.00201020636766293 Correlation Coef. (R): 0.999269586709141 R-Square: 0.998539706921858 Determination Coef. (DC): 0.997326463543269 Parameters Best Estimate -------------------- ------------- a 0.0585027054245842 ====== Output Results ===== File: Data file - 1 No Obs. y Cal. y 1 1.018748684 0.99044826487802 2 1.484227996 1.50307876691741 3 1.667487567 1.6894972796383 4 1.745852851 1.75335693506481 5 1.760339084 1.77620648682724 6 1.776919711 1.78458812214004 7 1.789486082 1.78786640578681 |
2楼2015-06-28 16:50:15
3楼2015-06-28 18:46:13
