应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 仿真模拟 » MATLAB » 【求助】polyfit的拟合效果不好，应该怎么修改？【已解决】
北京石油化工学院2026年研究生招生接收调剂公告
101/1返回列表
查看: 1810  |  回复: 9
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

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 ]
回复此楼

» 猜你喜欢

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

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

查看全部散金贴
1楼 2010-11-29 12:38:24
已阅   回复此楼   关注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 的时候的拟合结果
赞一下(1人) 回复此楼
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时的拟合图形
赞一下(1人) 回复此楼
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 ]
赞一下(1人) 回复此楼
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 的主题更新
101/1返回列表
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 一志愿苏州大学材料工程（085601）专硕有科研经历三项国奖两个实用型专利一项省级立项 +10 大火山小火山 2026-04-05 10/500 2026-04-05 19:58 by 小蜗牛小蜗牛
[考研] 一志愿北交大材料工程总分358求调剂 +6 cs0106 2026-04-05 6/300 2026-04-05 16:34 by imissbao
[考研] 0703化学 +12 妮妮ninicgb 2026-04-04 13/650 2026-04-05 10:46 by 啊俊！
[考研] 290求调剂085701 +10 1314捧花 2026-04-02 10/500 2026-04-05 10:19 by Sealedwind
[考研] 288求调剂，一志愿华南理工大学071005 +6 ioodiiij 2026-04-04 6/300 2026-04-05 10:09 by guoweigw
[考研] 求调剂 +3 小沢 2026-04-03 3/150 2026-04-05 09:10 by sihailian3
[考研] 材料295 +13 小英11 2026-04-03 14/700 2026-04-04 09:02 by 来看流星雨10
[考研] 292求调剂 +21 是妍子也是研子 2026-03-30 22/1100 2026-04-03 21:44 by qlm5820
[考研] 求调剂 +4 压力??大 2026-04-03 4/200 2026-04-03 21:36 by 啵啵啵0119
[考研] 274求调剂 +9 顺理成张 2026-04-03 10/500 2026-04-03 15:10 by 啊俊！
[考研] 求调剂 +3 usbdndj 2026-04-03 3/150 2026-04-03 14:10 by dxiaoxin
[考研] 262求调剂 +6 励志一定发文章 2026-04-02 7/350 2026-04-03 09:54 by linyelide
[考研] 生物学308分求调剂（一志愿华东师大） +6 相信必会光芒万� 2026-03-31 7/350 2026-04-02 23:16 by JourneyLucky
[考研] 一志愿复旦材料，英一专硕，总分357调剂 +4 1050389037 2026-04-02 5/250 2026-04-02 21:40 by dongzh2009
[考研] 初试301，代码085701环境工程，本硕一致，四六级已过，有二区一作，共发表5篇论文 +6 axibli 2026-04-01 6/300 2026-04-02 13:42 by Ecowxq666！
[考研] 279求调剂 +6 学而思兮知 2026-04-01 6/300 2026-04-02 09:16 by vgtyfty
[考研] 085602化学工程268分蹲调剂 +8 月照花林。 2026-04-01 8/400 2026-04-01 22:08 by 无际的草原
[考研] 298求调剂 +4 什么是胖头鱼 2026-03-30 6/300 2026-04-01 22:06 by 客尔美德
[考研] 085600 一志愿9 总分351 求调剂学校 +7 czhcz 2026-03-31 9/450 2026-04-01 19:24 by 无际的草原
[考研] 349求调剂 +6 吃的不少 2026-04-01 6/300 2026-04-01 17:55 by JYD2011
信息提示
关闭
请填处理意见
关闭