| 查看: 1100 | 回复: 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 |
回复此楼» 猜你喜欢
8月时间戳变的，举个手。玩一下，释放压力
已经有12人回复
-
我的国基提前知道中了，可是同事的操作让我实在接受不了，怎么会有这样的人
已经有8人回复
-
据悉今年马上要出结果了
已经有8人回复
-
关于代码变化问题，想知道的进来
已经有23人回复
-
FileCode能看出啥？
已经有6人回复
-
综述论文作为代表作会不会影响评审专家的印象分？
已经有9人回复
-
应该是93bebmhtak前后十一个字符比较关键
已经有20人回复
-
2026国自然放榜时间
已经有9人回复
-
这样的filecode谁见过
已经有22人回复
-
面上项目filecode邪修
已经有7人回复
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
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
4楼2016-03-22 23:23:24
|
多谢，我查了一下，1stopt可以给参数加约束，能不能帮忙再计算一下，如果方便的话可不可以传一下破解的1stopt，可以用钱买。 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,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,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,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.1,0.2,0.3,0.4,0.5,0.6,0.7,0.8,0.9,1,1.1,1.2,1.3,1.4,1.5,1.6,1.7,1.8]'; 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,4.5,5,5.7,6,6.3,6.5,7.5,7.2,6.8,6.2,5,4.2,3.5,3.4,1.1,0,0,0,0]'; a0=[8.758,0.073,-1.571,0.069,0.412,13.128,1.195,0.566]'; fun=@(a,x)a(1)+a(2)*x1+a(3)*x2+a(4)*exp(a(5)*x3)+a(6)/a(7)/sqrt(pi./2)*exp(-2.*(x4-a(8)).^2/a(8)/a(8)); [a,r,j]=nlinfit(x,S,fun,a0) y1=fun(a,x) a(2)>0;a(3)<0;a(4)>0;a(5)>0;a(6)>0;0.8<a(7)<1.2;0.5<a(8)<0.8; |
5楼2016-03-22 23:34:53
20