zyj8119

木蟲 (著名寫手)

應(yīng)助: 65 (初中生)
貴賓: 0.003
金幣: 915.1
帖子: 2936
在線: 1329.4小時(shí)
蟲號(hào): 664177

[交流] 【求助】polyfit的擬合效果不好，應(yīng)該怎么修改？【已解決】

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 ]

回復(fù)此樓

» 猜你喜歡

環(huán)境領(lǐng)域全國重點(diǎn)實(shí)驗(yàn)室招收博士1-2名已經(jīng)有3人回復(fù)
材料專碩306英一數(shù)二已經(jīng)有10人回復(fù)
301求調(diào)劑已經(jīng)有6人回復(fù)
一志愿天津大學(xué)化學(xué)工藝專業(yè)（081702）315分求調(diào)劑已經(jīng)有7人回復(fù)
302求調(diào)劑已經(jīng)有6人回復(fù)
26博士申請(qǐng) 已經(jīng)有3人回復(fù)
268求調(diào)劑已經(jīng)有3人回復(fù)
311求調(diào)劑已經(jīng)有10人回復(fù)
被我言中：新模板不強(qiáng)調(diào)格式了，假專家開始管格式了已經(jīng)有4人回復(fù)
303求調(diào)劑已經(jīng)有5人回復(fù)

» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:

投了篇ANGEW，編輯送審后建議修改之后投下邊的子刊，新人求助已經(jīng)有11人回復(fù)
請(qǐng)問在manuscriptcentral系統(tǒng)里面已經(jīng)submitted的manuscript怎樣再修改啊已經(jīng)有12人回復(fù)
氫化鋁鋰還原氰基用THF作溶劑，怎樣使氫化鋁鋰溶于溶劑呢如果不溶效果會(huì)不會(huì)不好已經(jīng)有6人回復(fù)
剛上傳了修改稿，祈禱啊已經(jīng)有9人回復(fù)
大家看看這個(gè)反應(yīng)怎么做啊已經(jīng)有22人回復(fù)
晶體指標(biāo)化不好，一般怎么改進(jìn)才能得到好的晶體，新手求教啊已經(jīng)有5人回復(fù)
求助修改說明的寫法已經(jīng)有4人回復(fù)
C18柱子沖洗效果不理想怎么回事？？已經(jīng)有4人回復(fù)
高考狀元為啥“混”不好？已經(jīng)有38人回復(fù)
急問�。�！第一次投稿，意見是小修，但是有一個(gè)關(guān)于修改的問題向大家請(qǐng)教已經(jīng)有38人回復(fù)

» 搶金幣啦！回帖就可以得到:

查看全部散金貼

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

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

zyj8119

木蟲 (著名寫手)

應(yīng)助: 65 (初中生)
貴賓: 0.003
金幣: 915.1
帖子: 2936
在線: 1329.4小時(shí)
蟲號(hào): 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. ...

綠色曲線開始的部分效果很差。。。。粉紅色的也是的。。。。

贊一下

回復(fù)此樓

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

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

筱晴

銅蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 372
帖子: 26
在線: 36小時(shí)
蟲號(hào): 1117100

★
zyj8119(金幣+1):謝謝參與
zyj8119(金幣+5):這個(gè)貌似根本就沒啥擬合啊。 2010-11-30 13:04:20

這個(gè)是n為6 的時(shí)候的擬合結(jié)果

贊一下(1人)

回復(fù)此樓

3樓2010-11-30 13:01:12

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

筱晴

銅蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 372
帖子: 26
在線: 36小時(shí)
蟲號(hào): 1117100

zyj8119(金幣+5):6的擬合效果跟不擬合都差別不大，呵呵。 2010-11-30 13:38:59

我就是把你程序中的n設(shè)置為6，出來的圖像就是這樣的

贊一下

回復(fù)此樓

4樓2010-11-30 13:37:11

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

change0618

鐵桿木蟲 (著名寫手)

應(yīng)助: 44 (小學(xué)生)
金幣: 17724.5
帖子: 2413
在線: 546.7小時(shí)
蟲號(hào): 496517

★
zyj8119(金幣+1):謝謝參與

N的值越大出現(xiàn)的龍格現(xiàn)象越嚴(yán)重，你只是顯示了原數(shù)據(jù)點(diǎn)的圖，并沒有顯示全區(qū)間的曲線。下面的就是n=6時(shí)的擬合圖形

贊一下(1人)

回復(fù)此樓

5樓2010-11-30 16:02:30

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

zyj8119

木蟲 (著名寫手)

應(yīng)助: 65 (初中生)
貴賓: 0.003
金幣: 915.1
帖子: 2936
在線: 1329.4小時(shí)
蟲號(hào): 664177

引用回帖:

Originally posted by change0618 at 2010-11-30 16:02:30:
N的值越大出現(xiàn)的龍格現(xiàn)象越嚴(yán)重，你只是顯示了原數(shù)據(jù)點(diǎn)的圖，并沒有顯示全區(qū)間的曲線。下面的就是n=6時(shí)的擬合圖形

怎么上面的那個(gè)人n=6的時(shí)候，擬合與不擬合差別很�。�

贊一下

回復(fù)此樓

6樓2010-11-30 16:08:57

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

change0618

鐵桿木蟲 (著名寫手)

應(yīng)助: 44 (小學(xué)生)
金幣: 17724.5
帖子: 2413
在線: 546.7小時(shí)
蟲號(hào): 496517

★ ★
nono2009(金幣+2):鼓勵(lì)應(yīng)助 2010-12-01 08:18:59

因?yàn)樗划嬃肆鶄€(gè)擬合值的折線圖，n=6時(shí)，擬合值跟數(shù)據(jù)點(diǎn)基本重合，但是非數(shù)據(jù)點(diǎn)就不怎么好了。數(shù)值分析課本上說過多項(xiàng)式擬合N值越大，龍格現(xiàn)象越嚴(yán)重，所以采用三次分段插值擬合。

[ Last edited by change0618 on 2010-11-30 at 16:19 ]

贊一下(1人)

回復(fù)此樓

7樓2010-11-30 16:15:56

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

zyj8119

木蟲 (著名寫手)

應(yīng)助: 65 (初中生)
貴賓: 0.003
金幣: 915.1
帖子: 2936
在線: 1329.4小時(shí)
蟲號(hào): 664177

引用回帖:

Originally posted by change0618 at 2010-11-30 16:15:56:
因?yàn)樗划嬃肆鶄€(gè)擬合值的折線圖，n=6時(shí)，擬合值跟數(shù)據(jù)點(diǎn)基本重合，但是非數(shù)據(jù)點(diǎn)就不怎么好了

我在cfttools里面，用GUI中樣條函數(shù)擬合了，效果還行，但是不知道怎么把格式搞成與plot一樣的

legend里面多余的那幾條線的標(biāo)識(shí)不知道怎么去掉。。。

贊一下

回復(fù)此樓

8樓2010-11-30 16:20:51

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

change0618

鐵桿木蟲 (著名寫手)

應(yīng)助: 44 (小學(xué)生)
金幣: 17724.5
帖子: 2413
在線: 546.7小時(shí)
蟲號(hào): 496517

zyj8119(金幣+4):謝謝提示。。。 2010-11-30 16:23:44

可以分別試試 interp1 spline csaps spap2 spaps 等函數(shù)

贊一下

回復(fù)此樓

9樓2010-11-30 16:22:09

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

change0618

鐵桿木蟲 (著名寫手)

應(yīng)助: 44 (小學(xué)生)
金幣: 17724.5
帖子: 2413
在線: 546.7小時(shí)
蟲號(hào): 496517

zyj8119(金幣+4):非常感謝，此貼結(jié)貼。。。 2010-11-30 16:59:48

cftool用的是fit函數(shù)，'smoothingspline‘選項(xiàng)

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

贊一下

回復(fù)此樓

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

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

最具人氣熱帖推薦 [查看全部]	作者	回/看	最后發(fā)表

[基金申請(qǐng)] 被我言中：新模板不強(qiáng)調(diào)格式了，假專家開始管格式了 +4	beefly 2026-03-14	4/200	2026-03-17 22:04 by 黃鳥于飛Chao
[考研] 材料，紡織，生物（0856、0710），化學(xué)招生啦 +3	Eember. 2026-03-17	7/350	2026-03-17 20:20 by 花125533
[考研] 341求調(diào)劑 +5	搗蛋豬豬 2026-03-11	7/350	2026-03-17 19:09 by 搗蛋豬豬
[考研] 344求調(diào)劑 +4	knight344 2026-03-16	4/200	2026-03-17 17:27 by ruiyingmiao
[考研] 332求調(diào)劑 +6	Zz版 2026-03-13	6/300	2026-03-17 17:03 by ruiyingmiao
[考研] 290求調(diào)劑 +3	p asserby. 2026-03-15	4/200	2026-03-17 16:35 by wangkm
[考研] 一志愿蘇州大學(xué)材料工程（085601）專碩有科研經(jīng)歷三項(xiàng)國獎(jiǎng)兩個(gè)實(shí)用型專利一項(xiàng)省級(jí)立項(xiàng) +6	大火山小火山 2026-03-16	8/400	2026-03-17 15:05 by 無懈可擊111
[考研] 283求調(diào)劑 +3	聽風(fēng)就是雨； 2026-03-16	3/150	2026-03-17 07:41 by 熱情沙漠
[考研] 286求調(diào)劑 +3	lemonzzn 2026-03-16	5/250	2026-03-16 20:43 by lemonzzn
[考研] 環(huán)境工程調(diào)劑 +6	大可digkids 2026-03-16	6/300	2026-03-16 17:16 by barlinike
[考研] 材料與化工一志愿南昌大學(xué)327求調(diào)劑推薦 +7	Ncdx123456 2026-03-13	8/400	2026-03-16 12:15 by karry wen
[考研] 求老師收留調(diào)劑 +4	jiang姜66 2026-03-14	5/250	2026-03-15 20:11 by Winj1e
[考研] 26考研一志愿中國石油大學(xué)(華東)305分求調(diào)劑 +3	嘉年新程 2026-03-15	3/150	2026-03-15 13:58 by 哈哈哈哈嘿嘿嘿
[考研] 085601材料工程315分求調(diào)劑 +3	yang_0104 2026-03-15	3/150	2026-03-15 10:58 by peike
[考研] 本科南京大學(xué)一志愿川大藥學(xué)327 +3	麥田耕者 2026-03-14	3/150	2026-03-14 20:04 by 外星文明
[考研] 材料與化工 323 英一+數(shù)二+物化，一志愿：哈工大本人本科雙一流 +4	自由的_飛翔 2026-03-13	5/250	2026-03-14 19:39 by hmn_wj
[考研] 復(fù)試調(diào)劑 +3	呼呼？~+123456 2026-03-14	3/150	2026-03-14 16:53 by WTUChen
[考研] 304求調(diào)劑 +7	7712b 2026-03-13	7/350	2026-03-13 21:42 by peike
[考研] 307求調(diào)劑 +5	超級(jí)伊昂大王 2026-03-12	5/250	2026-03-13 15:56 by 棒棒球手
[考研] 工科材料085601 279求調(diào)劑 +8	困于星晨 2026-03-12	10/500	2026-03-13 15:42 by ms629