24小时热门版块排行榜

返回列表

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

zyj8119

木虫 (著名写手)

应助: 65 (初中生)
贵宾: 0.003
金币: 915.1
帖子: 2936
在线: 1329.4小时
虫号: 664177

[交流] 【求助】polyfit的拟合效果不好，应该怎么修改？【已解决】

CODE:

function dpolydengyaxuanzexing(n)

clc;

x=[238,258,278,298,318,338,358];

y1=[2.409,4.219,5.553,9.241,23.425,23.75,25.142];

y2=[3.060,5.270,9.421,13.423,25.909,31.125,31.833];

y3=[3.393,5.452,9.077,11.5,25.454,29.315,29.333];

y4=[3.803,5.381,7.078,9.667,25.3,28.5,28.666];

plot(x,y1,'gx')

hold on

plot(x,y2,'b*')

hold on

plot(x,y3,'r+')

hold on

plot(x,y4,'mo')

xlabel('temperature(K)');ylabel('adsorption capacity(mmol/g)');

p1=polyfit(x,y1,n);

p2=polyfit(x,y2,n);

p3=polyfit(x,y3,n);

p4=polyfit(x,y4,n);

y1=polyval(p1,x);

y2=polyval(p2,x);

y3=polyval(p3,x);

y4=polyval(p4,x);

plot(x,y1,'g-');

hold on

plot(x,y2,'b-')

hold on

plot(x,y3,'r-')

hold on

plot(x,y4,'m-')

legend('100kpa in MCM-41 without APTS','100kpa in MCM-41 with 15APTS',...

    '100kpa in MCM-41 with 30APTS','100kpa in MCM-41 with 45APTS');

hold off

end

得到的图画是这样的：

[ Last edited by nono2009 on 2010-12-1 at 08:19 ]

回复此楼

» 猜你喜欢

拟解决的关键科学问题还要不要写已经有11人回复
救命帖已经有5人回复
限项规定已经有5人回复
为什么nbs上溴没有产物点出现呢已经有9人回复
招博士已经有3人回复
存款400万可以在学校里躺平吗已经有35人回复
最失望的一年已经有18人回复
求推荐博导已经有4人回复
求推荐英文EI期刊已经有5人回复
疑惑？已经有5人回复

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

投了篇ANGEW，编辑送审后建议修改之后投下边的子刊，新人求助已经有11人回复
请问在manuscriptcentral系统里面已经submitted的manuscript怎样再修改啊已经有12人回复
氢化铝锂还原氰基用THF作溶剂，怎样使氢化铝锂溶于溶剂呢如果不溶效果会不会不好已经有6人回复
刚上传了修改稿，祈祷啊已经有9人回复
大家看看这个反应怎么做啊已经有22人回复
晶体指标化不好，一般怎么改进才能得到好的晶体，新手求教啊已经有5人回复
求助修改说明的写法已经有4人回复
C18柱子冲洗效果不理想怎么回事？？已经有4人回复
高考状元为啥“混”不好？已经有38人回复
急问！！！第一次投稿，意见是小修，但是有一个关于修改的问题向大家请教已经有38人回复

» 抢金币啦！回帖就可以得到:

查看全部散金贴

1楼 2010-11-29 12:38:24

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

change0618

铁杆木虫 (著名写手)

应助: 44 (小学生)
金币: 17724.5
帖子: 2413
在线: 546.7小时
虫号: 496517

zyj8119(金币+4):非常感谢，此贴结贴。。。 2010-11-30 16:59:48

cftool用的是fit函数，'smoothingspline‘选项

function dpolydengyaxuanzexing
clc;
x=[238,258,278,298,318,338,358];
y1=[2.409,4.219,5.553,9.241,23.425,23.75,25.142];
y2=[3.060,5.270,9.421,13.423,25.909,31.125,31.833];
y3=[3.393,5.452,9.077,11.5,25.454,29.315,29.333];
y4=[3.803,5.381,7.078,9.667,25.3,28.5,28.666];
plot(x,y1,'gx')
hold on
plot(x,y2,'b*')
hold on
plot(x,y3,'r+')
hold on
plot(x,y4,'mo')
xlabel('temperature(K)');ylabel('adsorption capacity(mmol/g)');
yy1=fit(x',y1','smoothingspline');
yy2=fit(x',y2','smoothingspline');
yy3=fit(x',y3','smoothingspline');
yy4=fit(x',y4','smoothingspline');
plot(yy1,'g-');
hold on
plot(yy2,'b-')
hold on
plot(yy3,'r-')
hold on
plot(yy4,'m-')
legend('100kpa in MCM-41 without APTS','100kpa in MCM-41 with 15APTS',...
'100kpa in MCM-41 with 30APTS','100kpa in MCM-41 with 45APTS');
hold off
end

赞一下

回复此楼

10楼2010-11-30 16:55:46

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

查看全部 10 个回答

zyj8119

木虫 (著名写手)

应助: 65 (初中生)
贵宾: 0.003
金币: 915.1
帖子: 2936
在线: 1329.4小时
虫号: 664177

引用回帖:

Originally posted by zyj8119 at 2010-11-29 12:38:24:
[code]function dpolydengyaxuanzexing(n)
clc;
x=[238,258,278,298,318,338,358];
y1=[2.409,4.219,5.553,9.241,23.425,23.75,25.142];
y2=[3.060,5.270,9.421,13.423,25.909,31.125,31.833];
y3=[3.393,5. ...

绿色曲线开始的部分效果很差。。。。粉红色的也是的。。。。

赞一下

回复此楼

2楼2010-11-29 12:40:18

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