24小时热门版块排行榜

返回列表

当前主题已经存档。

【有奖交流】积极回复本帖子，参与交流，就有机会分得作者 feiliao8170 的 3 个金币

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

feiliao8170

木虫 (小有名气)

应助: 0 (幼儿园)
金币: 2496.9
帖子: 155
在线: 214.5小时
虫号: 579548
注册: 2008-07-11
专业: 机械摩擦学与表面技术

[交流] 【求助】Matlab数据拟合

数据如下：
x             y1          y2          y3             y4
308       0.1          0.07       0.01       0.006
313       0.11       0.068       0.0095       0.0065
323       0.115       0.045       0.009       0.007
333       0.098       0.025       0.0089       0.0075
343       0.04       0.015       0.0092       0.0078
353       0.02       0.01       0.009       0.008
363       0.015       0.0095    0.0095       0.009
373       0.01       0.012       0.012       0.01

P3          123          1233       2837       4933

拟合方程为：
y=0.006+0.114*(2838*P3/(2838*P3+exp(-P1/(8.314472*x))/P2))
其中x自变量，y为因变量，P1,P2为待优化参数，对应四组Y数据P3为固定值；

对于单列数据，采用origin很容易就可以得到P1和P2的优化值，但是如何对四列数据进行优化得到唯一的P1和P2值？

希望高手指点，谢谢。

回复此楼

» 猜你喜欢

青椒八年已不青，大家都被折磨成啥样了？已经有11人回复
限项规定已经有9人回复
免疫学博士有名额，速联系已经有4人回复
交叉科学部支持青年基金，对三无青椒是个机会吗？已经有5人回复
国家基金申请书模板内插入图片不可调整大小? 已经有6人回复
国家级人才课题组招收2026年入学博士已经有5人回复
Fe3O4@SiO2合成已经有6人回复
青年基金C终止已经有4人回复
26申博求博导推荐-遥感图像处理方向已经有4人回复
西南交通大学国家级人才团队2026年博士研究生招生（考核制）—机械、材料、力学方向已经有3人回复

1楼 2009-10-11 15:52:28

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

change0618

铁杆木虫 (著名写手)

方丈大师

应助: 44 (小学生)
金币: 17724.5
红花: 17
帖子: 2413
在线: 546.7小时
虫号: 496517
注册: 2008-01-19
专业: 化学反应工程

★ ★ ★ ★
feiliao8170(金币+2,VIP+0):谢谢 10-13 19:53
nono2009(金币+2,VIP+0):谢谢热心帮助。欢迎常来交流。 10-21 13:26

将x,P3看作是自变量，[y1;y2;y3;y4]看作是的因变量，做非线性拟合嘛

但是好像拟合不好

%% x y1 y2 y3 y4
data=[308 0.10 0.07 0.01 0.006
313 0.11 0.068 0.0095 0.0065
323 0.115 0.045 0.009 0.007
333 0.098 0.025 0.0089 0.0075
343 0.04 0.015 0.0092 0.0078
353 0.02 0.01 0.009 0.008
363 0.015 0.0095 0.0095 0.009
373 0.01 0.012 0.012 0.01];

xdata=data(:,1);
xdata=repmat(xdata,size(data,2)-1,1);
y1=data(:,2);y2=data(:,3);y3=data(:,4);y4=data(:,5);
y=[y1;y2;y3;y4];
P3 =ones(size(data,1),1)*[123 1233 2837 4933];
P3=P3(: );
fx=@(P,x)0.006+0.114*(2838*x(:,2)./(2838*x(:,2)+exp(-P(1)./(8.314472*x(:,1)))./P(2)));

XX=[xdata P3]
P=[10 20];
[P,r] = nlinfit(XX,y,fx,P);