版塊導(dǎo)航: 正在加載中...

登錄注冊

應(yīng)《網(wǎng)絡(luò)安全法》要求，自2017年10月1日起，未進(jìn)行實(shí)名認(rèn)證將不得使用互聯(lián)網(wǎng)跟帖服務(wù)。為保障您的帳號能夠正常使用，請盡快對帳號進(jìn)行手機(jī)號驗(yàn)證，感謝您的理解與支持！

24小時(shí)熱門版塊排行榜

返回列表

17_李賞

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 935.1
帖子: 29
在線: 16.6小時(shí)
蟲號: 3162383
注冊: 2014-04-25
專業(yè): 工程熱力學(xué)

[求助] matlab曲線擬合及修正已有2人參與

本人新手，還未對matlab進(jìn)行深入學(xué)習(xí)，緊急求助~
已知：x=[8.9 11.7 13.2 13.4 14 14.2 14 10.7 21.671 23.487 20.902 16.507]；
y=[17.698 17.698 27.365 27.365 37.031 37.031 27.365 17.698 66.03 75.696 56.364 37.031]。
請問該如何用matlab擬合出多項(xiàng)式方程y=f（x），然后如何對方程進(jìn)行修正？

回復(fù)此樓

1樓 2014-06-27 11:38:58

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

回帖支持 ( 顯示支持度最高的前 50 名 )

cooooldog

鐵桿木蟲 (著名寫手)

ส็็็

應(yīng)助: 237 (大學(xué)生)
金幣: 6101.9
散金: 1114
紅花: 39
帖子: 1380
在線: 553.8小時(shí)
蟲號: 506699
注冊: 2008-02-18
專業(yè): 模式識別

【答案】應(yīng)助回帖

★
感謝參與，應(yīng)助指數(shù) +1
17_李賞: 金幣+1, ★有幫助 2014-06-29 14:03:46

Matlab不是很熟，我用gnuplot+CTeX擬合了一下，湊合著看看吧：

代碼：

CODE:

\documentclass[border=10pt]{standalone}

\usepackage{siunitx}

\usepackage{pgfplots}

\pgfplotsset{compat=newest}

\usepackage{pgfplotstable}

\usepackage{filecontents}

\begin{filecontents*}{myData.dat}

X     Y

8.900        17.698

11.700        17.698

13.200        27.365

13.400        27.365

14.000        37.031

14.200        37.031

14.000        27.365

10.700        17.698

21.671        66.030

23.487        75.696

20.902        56.364

16.507        37.031

\end{filecontents*}

\begin{document}

\begin{tikzpicture}

\pgfplotsset{

width=10cm,

legend style={font=\scriptsize}

}

\begin{axis}[%

xlabel=$x$,

ylabel=$y$,

ymin=0,

xmin=0,

scaled y ticks=base 10:0,

legend cell align = left,

legend pos = north west

]

\addplot[only marks] table {mydata.dat};

\addlegendentry{scatter plot};

\addplot+[no markers,red] table [y={create col/linear regression={y=Y}}]{myData.dat};

\addlegendentry{%

Linear regression:  $y=\pgfmathprintnumber{\pgfplotstableregressiona} x

\pgfmathprintnumber[print sign]{\pgfplotstableregressionb}$};

% polynomial fit

\addplot [no markers, blue] gnuplot [raw gnuplot] { % allows arbitrary gnuplot commands

f(x) = a*x**2+b*x+c;     % Define the function to fit

a=10;b=10; c=10;         % Set reasonable starting values here

fit f(x) 'myData.dat' u 1:2 via a,b,c; % Select the file, starts at col 1 and two variables

plot [x=0:28] f(x);    % Specify the range to plot

set print "parameters.dat";  % Open a file to save the parameters

print a, b, c;                  % Write the parameters to file

};

\addlegendentry{\pgfplotstableread{parameters.dat}\parameters % Open the file Gnuplot wrote

\pgfplotstablegetelem{0}{0}\of\parameters \pgfmathsetmacro\paramA{\pgfplotsretval} % Get first element, save into \paramA

\pgfplotstablegetelem{0}{1}\of\parameters \pgfmathsetmacro\paramB{\pgfplotsretval}

\pgfplotstablegetelem{0}{2}\of\parameters \pgfmathsetmacro\paramC{\pgfplotsretval}

polynomial:$y=\pgfmathprintnumber{\paramA} x^2 \pgfmathprintnumber[print sign]{\paramB} x \pgfmathprintnumber[print sign]{\paramC}$

}

\end{axis}

\end{tikzpicture}

\end{document}

2014-06-27_21-50-07.png

贊一下(2人)

回復(fù)此樓

ส็็็็็็็็็็็็็็็็็็็็

5樓2014-06-27 21:50:27

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

阿兵將心

木蟲 (正式寫手)

清蟲

應(yīng)助: 2 (幼兒園)
金幣: 2966.5
散金: 50
紅花: 2
帖子: 398
在線: 198.2小時(shí)
蟲號: 1762729
注冊: 2012-04-18
性別: GG
專業(yè): 同步輻射技術(shù)及其應(yīng)用

matlab里面有選項(xiàng)來擬合的，高斯擬合，多項(xiàng)式擬合什么的都有，你百度就可以看到的。

[ 發(fā)自小木蟲客戶端 ]

贊一下(1人)

回復(fù)此樓

學(xué)術(shù)離我遙遠(yuǎn)嗎

2樓2014-06-27 15:56:00

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

zhangzhiguang

至尊木蟲 (文壇精英)

應(yīng)助: 121 (高中生)
金幣: 53088.8
紅花: 24
沙發(fā): 222
帖子: 26519
在線: 881.4小時(shí)
蟲號: 492270
注冊: 2008-01-08
性別: GG
專業(yè): 基礎(chǔ)物理學(xué)

【答案】應(yīng)助回帖

如果x不能調(diào)整，那就是一個(gè)較為復(fù)雜的非線性擬合，其擬合次數(shù)一般較高，你可以從4次開始不斷試一下

贊一下(1人)

回復(fù)此樓

好好學(xué)習(xí),天天向上!

8樓2014-06-28 11:31:55

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

普通回帖

17_李賞

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 935.1
帖子: 29
在線: 16.6小時(shí)
蟲號: 3162383
注冊: 2014-04-25
專業(yè): 工程熱力學(xué)

引用回帖:

2樓: Originally posted by 阿兵將心 at 2014-06-27 15:56:00
matlab里面有選項(xiàng)來擬合的，高斯擬合，多項(xiàng)式擬合什么的都有，你百度就可以看到的。

您好！謝謝您關(guān)注我的帖子，擬合函數(shù)求出來以后，應(yīng)該怎么進(jìn)行修正，您知道嗎？

贊一下

回復(fù)此樓

3樓2014-06-27 17:13:33

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

阿兵將心

木蟲 (正式寫手)

清蟲

應(yīng)助: 2 (幼兒園)
金幣: 2966.5
散金: 50
紅花: 2
帖子: 398
在線: 198.2小時(shí)
蟲號: 1762729
注冊: 2012-04-18
性別: GG
專業(yè): 同步輻射技術(shù)及其應(yīng)用

引用回帖:

3樓: Originally posted by 17_李賞 at 2014-06-27 17:13:33
您好！謝謝您關(guān)注我的帖子，擬合函數(shù)求出來以后，應(yīng)該怎么進(jìn)行修正，您知道嗎？...

我還沒有仔細(xì)看過這方面的東西，有文獻(xiàn)資料嗎？為何要修正呢，擬合不就是修正嗎

[ 發(fā)自小木蟲客戶端 ]

贊一下

回復(fù)此樓

學(xué)術(shù)離我遙遠(yuǎn)嗎

4樓2014-06-27 20:36:36

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

zhangzhiguang

至尊木蟲 (文壇精英)

應(yīng)助: 121 (高中生)
金幣: 53088.8
紅花: 24
沙發(fā): 222
帖子: 26519
在線: 881.4小時(shí)
蟲號: 492270
注冊: 2008-01-08
性別: GG
專業(yè): 基礎(chǔ)物理學(xué)

【答案】應(yīng)助回帖

感謝參與，應(yīng)助指數(shù) +1

可以使用polyfit來進(jìn)行擬合。但是你給的數(shù)據(jù)比較雜亂，例如x有升有降，不知道是否可以調(diào)整順序。如果能夠調(diào)整的話可以這樣

贊一下

回復(fù)此樓

好好學(xué)習(xí),天天向上!

6樓2014-06-28 11:28:00

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

zhangzhiguang

至尊木蟲 (文壇精英)

應(yīng)助: 121 (高中生)
金幣: 53088.8
紅花: 24
沙發(fā): 222
帖子: 26519
在線: 881.4小時(shí)
蟲號: 492270
注冊: 2008-01-08
性別: GG
專業(yè): 基礎(chǔ)物理學(xué)

【答案】應(yīng)助回帖

★ ★ ★ ★
17_李賞: 金幣+4, ★★★★★最佳答案 2014-06-29 14:03:29

>> x=[8.9 11.7 13.2 13.4 14 14.2 14 10.7 21.671 23.487 20.902 16.507];
>> y=[17.698 17.698 27.365 27.365 37.031 37.031 27.365 17.698 66.03 75.696 56.364 37.031];
>> p=polyfit(x,y,2)
>> y1=polyval(p,x);
>> plot(x,y,'o',x,y1).其中 p=polyfit(x,y,2)中的2是擬合項(xiàng)的次數(shù)，可以修改，最高可到10次

贊一下

回復(fù)此樓

好好學(xué)習(xí),天天向上!

7樓2014-06-28 11:29:39

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

17_李賞

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 935.1
帖子: 29
在線: 16.6小時(shí)
蟲號: 3162383
注冊: 2014-04-25
專業(yè): 工程熱力學(xué)

引用回帖:

4樓: Originally posted by 阿兵將心 at 2014-06-27 20:36:36
我還沒有仔細(xì)看過這方面的東西，有文獻(xiàn)資料嗎？為何要修正呢，擬合不就是修正嗎
...

可能我的想法太過理想化了，我在想，擬合曲線生成以后，再編一段調(diào)試程序，調(diào)換各階系數(shù)，使更多的點(diǎn)落在曲線上或者更接近于曲線。我再找找相關(guān)資料吧，如果發(fā)現(xiàn)了，再來找您交流哈！

贊一下

回復(fù)此樓

9樓2014-06-29 14:09:00

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

17_李賞

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 935.1
帖子: 29
在線: 16.6小時(shí)
蟲號: 3162383
注冊: 2014-04-25
專業(yè): 工程熱力學(xué)

引用回帖:

8樓: Originally posted by zhangzhiguang at 2014-06-28 11:31:55
如果x不能調(diào)整，那就是一個(gè)較為復(fù)雜的非線性擬合，其擬合次數(shù)一般較高，你可以從4次開始不斷試一下

這個(gè)例子中的x不能單獨(dú)調(diào)整，x和y值是對應(yīng)的，x調(diào)整為從小到大的順序，則y值也跟著相應(yīng)的x值調(diào)整到相應(yīng)位置，即點(diǎn)的位置固定，驗(yàn)證得到的擬合曲線是同一條。謝謝您提供的思路！

贊一下

回復(fù)此樓

10樓2014-06-29 15:52:42

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

相關(guān)版塊跳轉(zhuǎn) 我要訂閱樓主 17_李賞的主題更新

返回列表

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

[考研] 一志愿蘇州大學(xué)材料工程（085601）專碩有科研經(jīng)歷三項(xiàng)國獎兩個(gè)實(shí)用型專利一項(xiàng)省級立項(xiàng) +6	大火山小火山 2026-03-16	8/400	2026-03-17 15:05 by 無懈可擊111
[考研] 302求調(diào)劑 +4	小賈同學(xué)123 2026-03-15	8/400	2026-03-17 10:33 by 小賈同學(xué)123
[考研] 化學(xué)調(diào)劑0703 +8	啊我我的 2026-03-11	8/400	2026-03-16 17:23 by 我的船我的海
[考研] 283求調(diào)劑 +10	小樓。 2026-03-12	14/700	2026-03-16 16:08 by 13811244083
[考研] 0703化學(xué)調(diào)劑 290分有科研經(jīng)歷，論文在投 +7	膩膩gk 2026-03-14	7/350	2026-03-16 10:12 by houyaoxu
[考研] 326求調(diào)劑 +3	mlpqaz03 2026-03-15	3/150	2026-03-16 07:33 by Iveryant
[考研] 085601材料工程315分求調(diào)劑 +3	yang_0104 2026-03-15	3/150	2026-03-15 10:58 by peike
[考研] 288求調(diào)劑 +4	奇點(diǎn)0314 2026-03-14	4/200	2026-03-14 23:04 by JourneyLucky
[考研] 080500，材料學(xué)碩302分求調(diào)劑學(xué)校 +4	初識可樂 2026-03-14	5/250	2026-03-14 21:08 by peike
[考研] 中科大材料專碩319求調(diào)劑 +3	孟鑫材料 2026-03-13	3/150	2026-03-14 18:10 by houyaoxu
[考研] 一志愿哈工大材料324分求調(diào)劑 +5	閆旭東 2026-03-14	5/250	2026-03-14 14:53 by 木瓜膏
[考研] 266求調(diào)劑 +4	學(xué)員97LZgn 2026-03-13	4/200	2026-03-14 08:37 by zhukairuo
[考研] 招收0805（材料）調(diào)劑 +3	18595523086 2026-03-13	3/150	2026-03-14 00:33 by 123%、
[考研] 0703化學(xué)調(diào)劑 +4	快樂的香蕉 2026-03-11	4/200	2026-03-13 22:41 by JourneyLucky
[考研] 304求調(diào)劑 +6	Mochaaaa 2026-03-12	7/350	2026-03-13 22:18 by 星空星月
[考研] 26調(diào)劑/材料科學(xué)與工程/總分295/求收留 +9	2026調(diào)劑俠 2026-03-12	9/450	2026-03-13 20:46 by 18595523086
[考研] 材料專碩350 求調(diào)劑 +4	王金科 2026-03-12	4/200	2026-03-13 16:02 by ruiyingmiao
[考研] 工科材料085601 279求調(diào)劑 +8	困于星晨 2026-03-12	10/500	2026-03-13 15:42 by ms629
[考研] 290求調(diào)劑 +7	ADT 2026-03-12	7/350	2026-03-13 15:17 by JourneyLucky
[考研] 085600材料與化工 309分請求調(diào)劑 +7	dtdxzxx 2026-03-12	8/400	2026-03-13 14:43 by jxchenghu

24小時(shí)熱門版塊排行榜

[求助] matlab曲線擬合及修正 已有2人參與

» 收錄本帖的淘帖專輯推薦

» 猜你喜歡

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

【答案】應(yīng)助回帖

【答案】應(yīng)助回帖

【答案】應(yīng)助回帖

【答案】應(yīng)助回帖

[求助] matlab曲線擬合及修正已有2人參與

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