| 查看: 23689 | 回复: 73 | ||||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||||
dshuanger木虫 (小有名气)
|
[求助]
matlab非线性拟合模型参数，大神你在哪！
|
|||
|
本人非计算模拟专业，小白一枚，但最近需要模拟参数，不知如何运用matlab求参数，顾求助各位大神。 模型表达式： y=a/E1 +a/E2[1-exp(-x/t)]+a*x/E3 其中：x=[2，5，10，15，30，60，120，200，300，550，700]； y=[8.66667，10，10.66667，12.66667，13.33333，15，15.33333， 16.66667，17.33333，18.33333，20]； a=2 未知参数：E1，E2，t，E3 不知道能否通过matlab拟合参数 请各位大神能鼎力相助，最好能给出代码或者计算方法，谢谢大家！ [ Last edited by dshuanger on 2013-11-5 at 11:06 ] |
回复此楼» 猜你喜欢
深圳大学2026年秋博士招生-物理学-活性胶体方向-高永祥课题组
已经有18人回复
-
论物质与能量的统一模型及物理现象解释
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有205人回复
-
基于基元I统一理论的数学相关应用推导
已经有0人回复
-
基元I统一理论：宇宙本质、层级演化与修炼文明的本源规律
已经有1人回复
-
基元I理论下三大核心空间现象精准推导与细节解析
已经有0人回复
-
基于基元 I 统一理论的反重力理论推导
已经有0人回复
-
基于基元I统一理论的量子力学本源推导
已经有1人回复
-
推荐一款可以AI辅助写作的Latex编辑器SmartLatexEditor，超级好用，AI润色，全免费
已经有20人回复
-
【EI|Scopus 双检索】第六届智能机器人系统国际会议（ISoIRS 2026）
已经有0人回复
-
2026年第四届电动车与车辆工程国际会议（CEVVE 2026）
已经有0人回复
|
我自己实验数据的拟合代码： Parameters b0,b1,b2,b3,b4,b5,b6,b7,b8,b9; Variable X1,X2,X3,Y; Function Y=b0+b1*X1+b2*X2+b3*X3+b4*X1*X2+b5*X1*X3+b6*X2*X3+b7*X1^2+b8*X2^2+b9*X3^2; Data; 0.0000 8.0000 7.0000 63.9 0.0000 24.0000 11.5000 56.9 350.0000 40.0000 5.5000 95.4 350.0000 0.0000 11.5000 54.3 700.0000 16.0000 5.5000 92.7 700.0000 32.0000 10.000 92.9000 1050.0000 48.0000 4.0 95.80 1050.0000 0.0000 10.000 87.8 1400.0000 16.0000 4.0000 92.6000 1400.0000 32.0000 8.5000 98.700 1750.0000 48.0000 2.5 71.2000 1750.0000 8.0000 8.5000 95.1000 2100.0000 24.0000 2.5000 64.8000 2100.0000 40.0000 7.0000 67.2000 拟合数据结果： 均方差(RMSE): 1.54658329060161 残差平方和(SSE): 33.4868782467533 相关系数(R): 0.995247203455277 相关系数之平方(R^2): 0.990516995985549 决定系数(DC): 0.990516995985549 卡方系数(Chi-Square): 0.228909500020559 F统计(F-Statistic): 1253.42179901157 参数 最佳估算 ---------- ------------- b0 37.0459960300393 b1 0.0487926069483363 b2 1.50487196705523 b3 5.5454028094221 b4 -0.000558451632571275 b5 0.00153209146092655 b6 -0.0281388598195167 b7 -2.10419568056527E-5 b8 -0.0122539097071058 b9 -0.492832403832336 ====== 结果输出 ===== No 实测值y 计算值y 1 63.9 61.3939773 2 56.9 56.9333928 3 95.4 96.6663312 4 54.3 56.3074838 5 92.7 94.5904221 6 92.9 91.8798377 7 95.8 96.2627759 8 87.8 87.3372240 9 92.6 93.6201624 10 98.7 96.8095779 11 71.2 69.1925163 12 95.1 93.8336688 13 64.8 64.7666071 14 67.2 69.7060227 和文献里的不一致。文献的二次型多项式见21楼图1： |
25楼2015-01-21 14:15:38
jerkwin
专家顾问 (正式写手)
-
专家经验: +14 - 计算强帖: 1
- 应助: 454 (硕士)
- 金币: 20699.1
- 散金: 148
- 红花: 81
- 帖子: 813
- 在线: 2648.3小时
- 虫号: 1023452
- 注册: 2010-05-19
- 专业: 理论和计算化学
- 管辖: 分子模拟
2楼2013-11-05 11:23:25
feixiaolin
荣誉版主 (文坛精英)
-
专家经验: +518 - 应助: 942 (博后)
- 贵宾: 1.275
- 金币: 3430
- 散金: 58785
- 红花: 532
- 沙发: 11
- 帖子: 24215
- 在线: 2601.8小时
- 虫号: 2139575
- 注册: 2012-11-21
- 专业: 光学信息获取与处理
- 管辖: 数学
★
dingd: 金币+1, 谢谢分享。 2013-11-05 16:37:58
dingd: 金币+1, 谢谢分享。 2013-11-05 16:37:58
|
y=2/e1 +2/e2/(1-exp(-x/t))+a*x/3; // a=2 x=[2，5，10，15，30，60，120，200，300，550，700]； y=[8.66667，10，10.66667，12.66667，13.33333，15，15.33333， 16.66667，17.33333，18.33333，20]； 相关系数(R): 1 相关系数之平方(R^2): 1 决定系数(DC): 0.999999999996653 卡方系数(Chi-Square): 1.64243431014082E-11 参数 最佳估算 ---------- -------------- e1 =1.2090897803241 e2 =3.11591327892529 t =0.00209072546506921 ====== 输出结果 ===== No. Observed y Calculated y 1 8.66667 8.66666365993411 2 10 9.99999365993411 3 10.66667 10.6666636599341 4 12.66667 12.6666636599341 5 13.33333 13.3333236599341 6 15 14.9999936599341 7 15.33333 15.3333236599341 8 16.66667 16.6666636599341 9 17.33333 17.3333236599341 10 18.33333 18.3333236599341 11 20 19.9999936599341 |
3楼2013-11-05 11:54:42
feixiaolin
荣誉版主 (文坛精英)
-
专家经验: +518 - 应助: 942 (博后)
- 贵宾: 1.275
- 金币: 3430
- 散金: 58785
- 红花: 532
- 沙发: 11
- 帖子: 24215
- 在线: 2601.8小时
- 虫号: 2139575
- 注册: 2012-11-21
- 专业: 光学信息获取与处理
- 管辖: 数学
4楼2013-11-05 11:56:53