24小时热门版块排行榜    

查看: 1787  |  回复: 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 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 085700资源与环境308求调剂 +6 墨墨漠 2026-03-18 6/300 2026-03-20 00:02 by 23Postgrad
[论文投稿] 申请回稿延期一个月,编辑同意了。但系统上的时间没变,给编辑又写邮件了,没回复 10+3 wangf9518 2026-03-17 4/200 2026-03-19 23:55 by babero
[考研] 304求调剂 +3 曼殊2266 2026-03-18 3/150 2026-03-19 14:42 by peike
[考研] 085600材料与化工调剂 324分 +10 llllkkkhh 2026-03-18 12/600 2026-03-19 14:33 by llllkkkhh
[考研] 085410人工智能专硕317求调剂(0854都可以) +3 xbxudjdn 2026-03-18 3/150 2026-03-18 22:14 by zhq0425
[考研] 298-一志愿中国农业大学-求调剂 +7 手机用户 2026-03-17 7/350 2026-03-18 14:34 by vgtyfty
[考研] 297求调剂 +8 戏精丹丹丹 2026-03-17 8/400 2026-03-18 14:30 by laoshidan
[考研] 303求调剂 +4 睿08 2026-03-17 6/300 2026-03-18 11:01 by Iveryant
[考研] 0703化学调剂 +3 妮妮ninicgb 2026-03-17 3/150 2026-03-18 10:29 by macy2011
[考研] 278求调剂 +5 烟火先于春 2026-03-17 5/250 2026-03-18 08:43 by 星空星月
[考研] 考研求调剂 +3 橘颂. 2026-03-17 4/200 2026-03-17 21:43 by 有只狸奴
[考研] 308求调剂 +4 是Lupa啊 2026-03-16 4/200 2026-03-17 17:12 by ruiyingmiao
[考研] 332求调剂 +6 Zz版 2026-03-13 6/300 2026-03-17 17:03 by ruiyingmiao
[考研] 有没有道铁/土木的想调剂南林,给自己招师弟中~ +3 TqlXswl 2026-03-16 7/350 2026-03-17 15:23 by TqlXswl
[考博] 26申博 +4 八6八68 2026-03-16 4/200 2026-03-17 13:00 by 轻松不少随
[考研] 333求调剂 +3 文思客 2026-03-16 7/350 2026-03-16 18:21 by 文思客
[考研] 321求调剂 +5 大米饭! 2026-03-15 5/250 2026-03-16 16:33 by houyaoxu
[考研] 070300化学学硕求调剂 +6 太想进步了0608 2026-03-16 6/300 2026-03-16 16:13 by kykm678
[考研] 070303 总分349求调剂 +3 LJY9966 2026-03-15 5/250 2026-03-16 14:24 by xwxstudy
[考研] 289求调剂 +4 这么名字咋样 2026-03-14 6/300 2026-03-14 18:58 by userper
信息提示
请填处理意见