24小时热门版块排行榜    

查看: 1768  |  回复: 9

zyj8119

木虫 (著名写手)


[交流] 【求助】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 ]
回复此楼

» 猜你喜欢

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

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

查看全部散金贴

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

zyj8119

木虫 (著名写手)


引用回帖:
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的回帖

筱晴

铜虫 (初入文坛)



zyj8119(金币+1):谢谢参与
zyj8119(金币+5):这个貌似根本就没啥拟合啊。 2010-11-30 13:04:20

这个是n为6 的时候的拟合结果
3楼2010-11-30 13:01:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

筱晴

铜虫 (初入文坛)


zyj8119(金币+5):6的拟合效果跟不拟合都差别不大,呵呵。 2010-11-30 13:38:59
我就是把你程序中的n设置为6,出来的图像就是这样的
4楼2010-11-30 13:37:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

change0618

铁杆木虫 (著名写手)



zyj8119(金币+1):谢谢参与
N的值越大出现的龙格现象越严重,你只是显示了原数据点的图,并没有显示全区间的曲线。下面的就是n=6时的拟合图形

5楼2010-11-30 16:02:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zyj8119

木虫 (著名写手)


引用回帖:
Originally posted by change0618 at 2010-11-30 16:02:30:
N的值越大出现的龙格现象越严重,你只是显示了原数据点的图,并没有显示全区间的曲线。下面的就是n=6时的拟合图形


怎么上面的那个人n=6的时候,拟合与不拟合差别很小?
6楼2010-11-30 16:08:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

change0618

铁杆木虫 (著名写手)


★ ★
nono2009(金币+2):鼓励应助 2010-12-01 08:18:59
因为他只画了六个拟合值的折线图,n=6时,拟合值跟数据点基本重合,但是非数据点就不怎么好了。数值分析课本上说过多项式拟合N值越大,龙格现象越严重,所以采用三次分段插值拟合。

[ Last edited by change0618 on 2010-11-30 at 16:19 ]
7楼2010-11-30 16:15:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zyj8119

木虫 (著名写手)


引用回帖:
Originally posted by change0618 at 2010-11-30 16:15:56:
因为他只画了六个拟合值的折线图,n=6时,拟合值跟数据点基本重合,但是非数据点就不怎么好了

我在cfttools里面,用GUI中样条函数拟合了,效果还行,但是不知道怎么把格式搞成与plot一样的legend里面多余的那几条线的标识不知道怎么去掉。。。
8楼2010-11-30 16:20:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

change0618

铁杆木虫 (著名写手)


zyj8119(金币+4):谢谢提示。。。 2010-11-30 16:23:44
可以分别试试 interp1   spline  csaps  spap2  spaps 等函数
9楼2010-11-30 16:22:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

change0618

铁杆木虫 (著名写手)


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的回帖
相关版块跳转 我要订阅楼主 zyj8119 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 286求调剂 +3 lemonzzn 2026-03-16 5/250 2026-03-16 20:43 by lemonzzn
[考研] 328求调剂,英语六级551,有科研经历 +3 生物工程调剂 2026-03-16 4/200 2026-03-16 20:13 by Wangjingyue
[考研] 本科南京大学一志愿川大药学327 +3 麦田耕者 2026-03-14 3/150 2026-03-14 20:04 by 外星文明
[考研] 290求调剂 +4 @将就将就看 2026-03-10 8/400 2026-03-14 14:23 by 千千运气
[考研] 271求调剂 +10 生如夏花… 2026-03-11 10/500 2026-03-14 00:35 by 卖报员小雨
[考研] 311求调剂 +8 zchqwer 2026-03-10 8/400 2026-03-14 00:01 by JourneyLucky
[考研] 材料371求调剂 +9 鳄鱼? 2026-03-11 11/550 2026-03-13 22:53 by JourneyLucky
[考研] 332求调剂 +3 zjy101327 2026-03-11 6/300 2026-03-13 22:48 by JourneyLucky
[考研] 材料与化工求调剂一志愿 985 总分 295 +8 dream…… 2026-03-12 8/400 2026-03-13 22:17 by 星空星月
[考研] 0856材料与化工301求调剂 +5 奕束光 2026-03-13 5/250 2026-03-13 22:00 by 星空星月
[考研] (081700)化学工程与技术-298分求调剂 +12 11啦啦啦 2026-03-11 35/1750 2026-03-13 21:25 by JourneyLucky
[考研] 332求调剂 +3 Zz版 2026-03-13 3/150 2026-03-13 20:36 by 18595523086
[考研] 302求调剂 +6 负心者当诛 2026-03-11 6/300 2026-03-13 16:11 by JourneyLucky
[考研] 材料专硕350 求调剂 +4 王金科 2026-03-12 4/200 2026-03-13 16:02 by ruiyingmiao
[考研] 290求调剂 +7 ADT 2026-03-12 7/350 2026-03-13 15:17 by JourneyLucky
[考研] 328化工专硕求调剂 +4 。,。,。,。i 2026-03-12 4/200 2026-03-13 14:44 by JourneyLucky
[考研] 268求调剂 +4 好运连绵不绝 2026-03-12 4/200 2026-03-13 10:45 by hyswxzs
[考研] 工科0856专硕化学工程269能调剂吗 +10 我想读研11 2026-03-10 10/500 2026-03-13 10:14 by Yuyi.
[考研] 270求调剂 085600材料与化工专硕 +3 YXCT 2026-03-11 3/150 2026-03-13 10:13 by houyaoxu
[考研] 化工0817调剂 +8 灿若星晨 2026-03-10 8/400 2026-03-10 22:44 by 星空星月
信息提示
请填处理意见