24小时热门版块排行榜

返回列表

myw052

金虫 (正式写手)

博士

应助: 9 (幼儿园)
金币: 542
散金: 300
红花: 2
帖子: 409
在线: 164.1小时
虫号: 1320903
注册: 2011-06-11
性别: GG
专业: 固体力学

[求助] 散点图函数拟合

本人需要对下面的数据进行拟合，拟合函数为x=B1+B2*lg(y-B3),其中lg表示以10为底的log函数。需要拟合出B1,B2,B3的数值。

式中x为下面数据的第一列，y为下列数据的第二列。
能用matlab编制程序最好，不能也行，只要能拟合出来系数B1 B2 B3的值也行。

希望版主能出来帮下忙。。。
3.004  640.19
3.297  591.08
3.345  583.34
3.465  564.69
3.585  546.54
3.769  519.82
3.959  493.58
4.018  485.67
4.278  452.43
4.475  400.28
4.802  379.72
5.036  367.44
5.151  368.26
5.282  367.44
5.315  341.13
5.386  356.70
5.398  369.09
5.729  356.70
5.866  356.70
6.296  359.17
6.302  343.58

回复此楼

» 猜你喜欢

职称评审没过，求安慰已经有21人回复
垃圾破二本职称评审标准已经有15人回复
投稿Elsevier的Neoplasia杂志，到最后选publishing options时页面空白，不能完成投稿已经有20人回复
EST投稿状态问题已经有7人回复
毕业后当辅导员了，天天各种学生超烦已经有4人回复
聘U V热熔胶研究人员已经有10人回复
求助文献已经有3人回复
投稿返修后收到这样的回复，还有希望吗已经有8人回复
三无产品还有机会吗已经有6人回复

» 本主题相关价值贴推荐，对您同样有帮助:

关于拟合的思想的一些交流已经有8人回复
请问一个有理函数的曲线拟合过程如何做，给出数据，请求拟合，已经有14人回复
很没有规律的散点图用什么方法拟合？已经有5人回复
根据数据，如何得到拟合公式已经有7人回复
求助matlab隐函数参数拟合已经有16人回复
大家帮我看看这个图怎么做啊已经有9人回复
求助Origin8.0作图的问题已经有2人回复
请教如何用origin 对柱状图数据拟合已经有3人回复
【一个数据拟合的问题，比较麻烦】已经有17人回复
origin 作相关图时怎么加上斜线。已经有9人回复
一条曲线如何进行分段拟合已经有13人回复
【讨论】急！！请教大家Origin7.5的非线性拟合中的初始值设定！已经有2人回复
excel绘图中计算的r2与用回归分析工具计算的不一致已经有9人回复
请教关于origin中自定义函数非线性回归拟合问题！已经有10人回复
【求助】如何把origin中曲线的拟合线延长？已经有7人回复

天天快乐

1楼2012-04-15 11:19:33

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

dbb627

荣誉版主 (著名写手)

专家经验: +4
应助: 289 (大学生)
贵宾: 0.589
金币: 24640.4
散金: 551
红花: 61
沙发: 1
帖子: 1246
在线: 1794.8小时
虫号: 149791
注册: 2005-12-29
性别: GG
专业: 污染控制化学
管辖: 计算模拟

【答案】应助回帖

感谢参与，应助指数 +1

>> A=[3.004  640.19
3.297  591.08
3.345  583.34
3.465  564.69
3.585  546.54
3.769  519.82
3.959  493.58
4.018  485.67
4.278  452.43
4.475  400.28
4.802  379.72
5.036  367.44
5.151  368.26
5.282  367.44
5.315  341.13
5.386  356.70
5.398  369.09
5.729  356.70
5.866  356.70
6.296  359.17
6.302  343.58];
x=A(:,1);y=A(:,2);
st_ = [2 2 3] ;
ft_ = fittype('B1+B2*log10(y-B3)','dependent',{'x'},'independent',{'y'},'coefficients',{'B1', 'B2','B3'});
[cf_,good]= fit(y,x,ft_ ,'Startpoint',st_)
h_ = plot(cf_,'fit',0.95);
legend off;  % turn off legend from plot method call
set(h_(1),'Color',[1 0 0],...
   'LineStyle','-', 'LineWidth',2,...
   'Marker','none', 'MarkerSize',6);
hold on,plot(y,x,'*')

cf_ =

   General model:
   cf_(y) = B1+B2*log10(y-B3)
   Coefficients (with 95% confidence bounds):
   B1 =    11.65  (5.446, 17.86)
   B2 =    -3.371  (-5.779, -0.9636)
   B3 =    292.6  (210.3, 375)

good =

         sse: 1.7251
   rsquare: 0.9180
         dfe: 18
adjrsquare: 0.9089
      rmse: 0.3096