| 查看: 565 | 回复: 2 | ||||
坤少pop新虫 (正式写手)
|
[求助]
1stopt运行不了，望大家赐教！ 已有2人参与
|
|
Title "非线性拟合";Parameters b(1:10);SharedModel;Variable x,y1,y2,y3;Function y1 = b7-((1/b8+b5/b9)*((b1*b2*exp(b3*x)/(b2-b1+b1*exp(b3*x)))-b1))-(b2/b3)*(b6/b9+b10)*log((b2-b1+b1*exp(b3*x))/b2);y2 = b4+b5*((b1*b2*exp(b3*x)/(b2-b1+b1*exp(b3*x)))-b1)+(b6*b2/b3)*log((b2-b1+b1*exp(b3*x))/b2);y3 = b1*b2*exp(b3*x)/(b2-b1+b1*exp(b3*x));Data;//x,y1,y2,y30 28.21 0.39 1.156 27.65 0.92 2.212 26.96 1.97 4.118 22.88 2.41 7.2524 18.04 2.92 9.130 12.81 3.2 10.0536 8.71 3.42 10.6542 6.59 3.42 10.748 6.22 3.68 10.7554 5.86 3.43 1160 5.74 3.44 10.766 5.7 3.52 10.672 5.43 3.62 10.4578 5.26 3.52 10.184 4.14 3.41 10 @月只蓝 发自小木虫Android客户端 |
回复此楼» 猜你喜欢
浙大物院表面物理课题组长期招收普博及直博学生
已经有259人回复
-
【征稿】Springer Nature热电专题 Discover Materials / Discover Chemistry
已经有47人回复
-
物理学I论文润色/翻译怎么收费?
已经有263人回复
-
25届硕士申请，材料计算方向
已经有2人回复
-
课题组招1-2名硕士
已经有1人回复
-
机时越来越贵之后，我把一部分 DFT 流程搬到了 AI 计算平台上
已经有0人回复
-
斯坦福大学与Elsevier世界前2%学者新排名查询
已经有1人回复
-
27申博推荐：南京航空航天大学国际前沿院光学方向招收博士生！
已经有19人回复
-
华南师范大学+双一流学科物理学+广东省杰青团队2027年博士招生+机器学习驱动材料设计
已经有1人回复
-
apl投稿状态
已经有2人回复
-
华南师范大学+双一流学科物理学+广东省杰青团队2027年博士招生+机器学习驱动材料设计
已经有0人回复
月只蓝
主管区长 (职业作家)
-
专家经验: +1059 - 计算强帖: 8
- 应助: 1712 (讲师)
- 贵宾: 8.888
- 金币: 68182.7
- 散金: 1938
- 红花: 443
- 沙发: 4
- 帖子: 4373
- 在线: 3292.4小时
- 虫号: 1122189
- 注册: 2010-10-14
- 专业: 宇宙学
- 管辖: 计算模拟区
2楼2017-04-28 10:51:23
lipenggg
专家顾问 (正式写手)
-
专家经验: +65 - 应助: 122 (高中生)
- 金币: 4841.3
- 红花: 43
- 帖子: 722
- 在线: 435.7小时
- 虫号: 3929201
- 注册: 2015-06-17
- 性别: GG
- 专业: 流体力学
- 管辖: 计算模拟
【答案】应助回帖
|
仅供参考 迭代数: 21 计算用时(时:分:秒:微秒): 00:00:08:165 优化算法: 麦夸特法(Levenberg-Marquardt) + 通用全局优化法 计算结束原因: 达到收敛判定标准 均方差(RMSE): 161.57197020968 残差平方和(SSE): 939798.056067761 相关系数(R): 0.616840194121187 相关系数之平方(R^2): 0.380491825083464 决定系数(DC): 0.380504389127482 F统计(F-Statistic): 6.141837032816 参数 最佳估算 -------------------- ------------- b1 2.33716162201893E-19 b2 17.3958605321308 b3 8.46776051292736 b4 -2.0582990250142 b5 21.1186527230704 b6 -1.39858295194152 b7 10.890468507053 b8 -34615.4838100747 b9 -279078.556769959 b10 0.0290034992138533 ====== 输出结果 ===== No 实测y1 计算y1 实测y2 计算y2 实测y3 计算y3 1 28.21 10.890468507053 0.39 -2.0582990250142 1.156 2.33716162201893E-19 2 0.92 -0.333867754925933 2.212 -175.926090470188 26.96 17.3958605321308 3 4.12 10.890468507053 22.88 -2.05829902492817 2.41 4.10550916516905E-12 4 18.04 9.96054354576693 2.92 320.396750581853 9.13 17.3958577116967 5 3.2 7.15481184556644 10.0536 185.124490512736 8.71 17.3958605321308 6 10.65 10.8904685041219 6.59 -2.05828053229706 3.42 8.82560361355933E-7 7 6.22 8.19434553673991 3.68 235.243371700392 10.7554 17.3958605321308 8 3.43 10.6607023622293 1160 346.60961417123 5.74 17.0385627076581 9 10.77 10.890468503581 5.7 -2.05827711963598 3.52 1.04542876899124E-6 10 5.43 8.23471577717384 3.62 237.189736018359 10.4578 17.3958605321308 11 3.52 10.8754135740444 10.184 81.1938854301014 4.14 3.97743539181013 12 10 10.8904685043599 0 -2.05828203374925 0 8.10903926099291E-7  y1.png  y2.png  y3.png |
3楼2017-05-02 10:56:20
50