| 网站首页 | 业界新闻 | 小组 | 威客 | 人才 | 下载频道 | 博客 | 代码贴 | 在线编程 | 编程论坛
欢迎加入我们,一同切磋技术
用户名:   
 
密 码:  
共有 2160 人关注过本帖
标题:求助:非线性拟合求解参数
只看楼主 加入收藏
florence
Rank: 1
等 级:新手上路
帖 子:16
专家分:0
注 册:2007-11-10
收藏
 问题点数:0 回复次数:6 
求助:非线性拟合求解参数
我有这样一个函数式:
y=a1*((a2*(a3*x1.^(-a4)+(1-a3)*x2.^(-a4)).^(a5/a4))+(1-a2)*x3.^(-a5)).^(-a6/a5)
现已知x1,x2,x3,y的值,
x1=4.7613 14.3462 21.9311 16.7593 15.7033 9.5024 8.6414 16.574 24.646 32.2199 38.6204
x2=496.25 451.23 415.66 383.29 350.89 309.21 273.36 220.11 177.65 150.9 151
x3=1.03 1.05 1.01 1.1 1.15 1.11 1.19 1.27 1.3 1.29 1.3
y=134.14 147.04 164.76 178.6 195.34 221.88 250.37 287.35 353.24 415.79 476.93

想求解其中的6个参数
我用nlinfit和lsqcurvefit初始值对结果影响很大,调整很多次,结果无法解释,望高手能指点迷津!

感激中.............
搜索更多相关主题的帖子: 非线性 参数 拟合 求解 
2007-11-12 10:44
ydgsl
Rank: 3Rank: 3
等 级:论坛游民
威 望:8
帖 子:107
专家分:20
注 册:2006-9-20
收藏
得分:0 
SSy = 1.3069e+005
a =[129.4777588, .4323169081e-5, .2617801050e-2, 202.3083741, -3.429133891, 3.412052175]
RSS = 2828.828306
MSe = 565.7657
R^2 = .978354
另:可能模型并不十分恰当。
2007-11-12 11:23
florence
Rank: 1
等 级:新手上路
帖 子:16
专家分:0
注 册:2007-11-10
收藏
得分:0 

谢谢大侠回复

这个函数形式是固定的,0<a2,a3<1, a4,a5>-1 , a6>0
我只能做以下处理:
y = ((a1*(a2*x1.^(-a3)+(1-a2)*x2.^(-a3)).^(a4/a3))+(1-a1)*x3.^(-a4)).^(-a5/a4)
其中0<a1,a2<1, a3,a4>-1 , a5>0
去掉原函数的a1,对a1另行处理,其他的没法处理

还有就是,我在使用 lsqcurvefit进行运算时经常出现这样的提示:Maximum number of function evaluations exceeded;
increase options.MaxFunEvals
我调整了几次options都不行,还是提示

2007-11-12 23:45
ydgsl
Rank: 3Rank: 3
等 级:论坛游民
威 望:8
帖 子:107
专家分:20
注 册:2006-9-20
收藏
得分:0 
第一个参数a6是必不可少的。
[a6 a1 a2 a3 a4 a5]=[.2735758395e-19, .1834991482e-4, 6.990171934, -9.670005208, 3.511673612, 129.4592106]
RSS=2014.392605
MSe=402.8785
R^2=.984586
这很可能是该问题的最优拟合。

2007-11-14 09:51
florence
Rank: 1
等 级:新手上路
帖 子:16
专家分:0
注 册:2007-11-10
收藏
得分:0 

谢谢ydgsl的回复
从结果上来看,这个组合是不能选择的
我做的这个拟合求解参数是为了选择x1,x2,x3之间的最佳组合
拟合优度是考虑的一方面,但不是全部,其参数值还要符合限制条件。
其它两组组和为:
(一)
x1=4.7613 14.3462 21.9311 16.7593 15.7033 9.5024 8.6414 16.574 24.646 32.2199 38.6204
x2=1.03 1.05 1.01 1.1 1.15 1.11 1.19 1.27 1.3 1.29 1.3
x3=496.25 451.23 415.66 383.29 350.89 309.21 273.36 220.11 177.65 150.9 151
(二)
x1=1.03 1.05 1.01 1.1 1.15 1.11 1.19 1.27 1.3 1.29 1.3
x2=496.25 451.23 415.66 383.29 350.89 309.21 273.36 220.11 177.65 150.9 151
x3=4.7613 14.3462 21.9311 16.7593 15.7033 9.5024 8.6414 16.574 24.646 32.2199 38.6204

2007-11-14 20:22
ydgsl
Rank: 3Rank: 3
等 级:论坛游民
威 望:8
帖 子:107
专家分:20
注 册:2006-9-20
收藏
得分:0 
若都符合的参数,可能的解为:
a =[114.3824422, .1817293910e-1, .6867645237e-59, 176.3446065, .1322030417, 3.429403762]
RSS = 4364.783850
MSe = 872.9568
R^2 =.966601
2007-11-15 18:55
florence
Rank: 1
等 级:新手上路
帖 子:16
专家分:0
注 册:2007-11-10
收藏
得分:0 

我不知道这样的结果是怎么得到的
我的意思是对于每一组数据对应一组参数值,然后从三组参数中选一组较为合理的。
6参数中a1结果过大,我还不能做出合理解释。

2007-11-16 22:39
快速回复:求助:非线性拟合求解参数
数据加载中...
 
   



关于我们 | 广告合作 | 编程中国 | 清除Cookies | TOP | 手机版

编程中国 版权所有,并保留所有权利。
Powered by Discuz, Processed in 0.030067 second(s), 9 queries.
Copyright©2004-2024, BCCN.NET, All Rights Reserved