| 查看: 1112 | 回复: 8 | ||||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||||
[求助]
MATLAB拟合 已有2人参与
|
||||
|
我用下面的数据和公式进行拟合，出现错误，麻烦大神解答 clear all x1=[10,20,30,40,50,40,40,40,40,40,40,40,40,40,40,40,40,40,40,40,40,40,40,40]'; x2=[8,8,8,8,8,4,6,8,10,8,8,8,8,8,8,8,8,8,8,8,8,8,8,8]'; x3=[0.5,0.5,0.5,0.5,0.5,0.5,0.5,0.5,0.5,0.5,1,2,0.5,0.5,0.5,0.5,0.5,0.5,0.5,0.5,0.5,0.5,0.5,0.5]'; x4=[0,0,0,0,0,0,0,0,0,0,0,0,0,0.25,0.5,0.75,0.8,0.9,1,1.2,1.3,1.4,1.5,2]'; x=[x1,x2,x3,x4]; S=[2,2.8,3.5,4.5,4.8,3,3.3,4.5,6.8,4.5,3,2.6,4.5,6,6.5,7.6,6.8,6.2,5,3.5,3.4,1.1,0,0]'; a0=[3.208,0.073,3.265,-0.431,0.069,0.417,11,1.076,0.565]'; fun=inline('a(1)+a(2)*x1+a(3)*x2.^a(4)+a(5)*exp(a(6)*x3)+a(7)/a(8)/sqrt(pi./2)*exp(-2*(x4-a(9)).^2/a(8)/a(8))','a','x'); [a,r,j]=nlinfit(x,S,fun,a0)  错误.png |
回复此楼» 猜你喜欢
帮忙看看fileCode
已经有8人回复
-
小木虫上这么多卖论文的，真有人买论文么？感觉没必要啊
已经有3人回复
-
FileCode能看出啥？
已经有15人回复
-
filecode
已经有8人回复
-
确定了，国自然21号放榜
已经有7人回复
-
综述论文作为代表作会不会影响评审专家的印象分？
已经有11人回复
-
关于代码变化问题，想知道的进来
已经有24人回复
-
为什么网上很多人说本周 12号出结果
已经有3人回复
-
静等基金结果
已经有16人回复
-
我的国基提前知道中了，可是同事的操作让我实在接受不了，怎么会有这样的人
已经有11人回复
dingd
铁杆木虫 (职业作家)
- 应助: 1641 (讲师)
- 金币: 15037.3
- 散金: 101
- 红花: 234
- 帖子: 3410
- 在线: 1223.7小时
- 虫号: 291104
- 注册: 2006-10-28
|
均方差(RMSE): 0.786006701903222 残差平方和(SSE): 19.1520025985402 相关系数(R): 0.934123874070345 相关系数之平方(R^2): 0.872587412108189 确定系数(DC): 0.872587412108189 卡方系数(Chi-Square): 2.82861233469515 F统计(F-Statistic): 25.7879879851766 参数 最佳估算 ---------- ------------- a2 0.0644163727838743 a3 -1.91118384935733E-16 a4 1.22839069041116E-12 a5 6.51673289707249E-11 a6 13.1414143445168 a7 1.19543637850916 a8 0.566009585348225 a1 -4.18607760247816 |
9楼2016-03-23 17:53:20
wurongjun
专家顾问 (职业作家)
-
专家经验: +831 - 数学EPI: 9
- 应助: 791 (博后)
- 贵宾: 0.308
- 金币: 24609
- 散金: 310
- 红花: 75
- 帖子: 3004
- 在线: 881.4小时
- 虫号: 1368482
- 注册: 2011-08-14
- 性别: GG
- 专业: 计算数学与科学工程计算
- 管辖: 数学
【答案】应助回帖
★ ★
感谢参与，应助指数 +1
zhw921116: 金币+2, ★★★很有帮助 2016-03-22 23:36:41
感谢参与，应助指数 +1
zhw921116: 金币+2, ★★★很有帮助 2016-03-22 23:36:41
|
把fun改一下就行啦! fun=inline('a(1)+a(2)*x(1)+a(3)*x(2).^a(4)+a(5)*exp(a(6)*x(3))+a(7)/a(8)/sqrt(pi./2)*exp(-2*(x(4)-a(9)).^2/a(8)/a(8))','a','x'); 运行结果: >> a = 3.2080 0.0730 3.2650 -0.4310 -0.0000 0.3309 11.0000 1.0760 0.5650 r = -1.9958 -1.1958 -0.4958 0.5042 0.8042 -0.9958 -0.6958 0.5042 2.8042 0.5042 -0.9958 -1.3958 0.5042 2.0042 2.5042 3.6042 2.8042 2.2042 1.0042 -0.4958 -0.5958 -2.8958 -3.9958 -3.9958 j = 1.0e+004 * Columns 1 through 7 0.0001 0.0010 0.0000 0.0003 2.0452 -0.0025 0 0.0001 0.0010 0.0000 0.0003 2.0452 -0.0025 0 0.0001 0.0010 0.0000 0.0003 2.0452 -0.0025 0 0.0001 0.0010 0.0000 0.0003 2.0452 -0.0025 0 0.0001 0.0010 0.0000 0.0003 2.0452 -0.0025 0 0.0001 0.0010 0.0000 0.0003 2.0452 -0.0025 0 0.0001 0.0010 0.0000 0.0003 2.0452 -0.0025 0 0.0001 0.0010 0.0000 0.0003 2.0452 -0.0025 0 0.0001 0.0010 0.0000 0.0003 2.0452 -0.0025 0 0.0001 0.0010 0.0000 0.0003 2.0452 -0.0025 0 0.0001 0.0010 0.0000 0.0003 2.0452 -0.0025 0 0.0001 0.0010 0.0000 0.0003 2.0452 -0.0025 0 0.0001 0.0010 0.0000 0.0003 2.0452 -0.0025 0 0.0001 0.0010 0.0000 0.0003 2.0452 -0.0025 0 0.0001 0.0010 0.0000 0.0003 2.0452 -0.0025 0 0.0001 0.0010 0.0000 0.0003 2.0452 -0.0025 0 0.0001 0.0010 0.0000 0.0003 2.0452 -0.0025 0 0.0001 0.0010 0.0000 0.0003 2.0452 -0.0025 0 0.0001 0.0010 0.0000 0.0003 2.0452 -0.0025 0 0.0001 0.0010 0.0000 0.0003 2.0452 -0.0025 0 0.0001 0.0010 0.0000 0.0003 2.0452 -0.0025 0 0.0001 0.0010 0.0000 0.0003 2.0452 -0.0025 0 0.0001 0.0010 0.0000 0.0003 2.0452 -0.0025 0 0.0001 0.0010 0.0000 0.0003 2.0452 -0.0025 0 Columns 8 through 9 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 |
2楼2016-03-22 11:08:52
dingd
铁杆木虫 (职业作家)
- 应助: 1641 (讲师)
- 金币: 15037.3
- 散金: 101
- 红花: 234
- 帖子: 3410
- 在线: 1223.7小时
- 虫号: 291104
- 注册: 2006-10-28
【答案】应助回帖
★ ★ ★
感谢参与，应助指数 +1
zhw921116: 金币+3, ★★★很有帮助 2016-03-22 23:37:08
感谢参与，应助指数 +1
zhw921116: 金币+3, ★★★很有帮助 2016-03-22 23:37:08
|
1stOpt拟合结果： 均方差(RMSE): 0.43390091962119 残差平方和(SSE): 4.51848019315474 相关系数(R): 0.976560633341951 相关系数之平方(R^2): 0.953670670593233 确定系数(DC): 0.953670670593233 卡方系数(Chi-Square): 1.05109921892536 F统计(F-Statistic): 40.4711232389677 参数 最佳估算 ---------- ------------- a1 -6.72233963538864 a2 0.0736026176846468 a3 5.69457619273679E-5 a4 4.82855796950308 a5 6.728128434041 a6 -2.75700141028568 a7 11.4573837792192 a8 1.11137672190465 a9 0.549488798892685 |
3楼2016-03-22 11:36:08
4楼2016-03-22 23:23:24
