yqfan

鐵蟲 (小有名氣)

應(yīng)助: 3 (幼兒園)
金幣: 164.3
紅花: 1
帖子: 85
在線: 24.7小時
蟲號: 989990
注冊: 2010-04-06
專業(yè): 橡膠及彈性體

[求助] 如何用已知函數(shù)去擬合實(shí)驗(yàn)數(shù)據(jù)，得出參數(shù)

已知實(shí)驗(yàn)數(shù)據(jù)
t=[1 2.1 5.0 10 15 20 30 40 50 60 70 80 90 100];
y=[6.0 6.08 6.2 6.4 6.47 6.58 6.8 6.85 6.9 6.93 6.95 7.0 7.05 7.1];
如何用函數(shù)去最佳擬合實(shí)驗(yàn)數(shù)據(jù)并求出參數(shù)。本人剛學(xué)matlab，用了很長時間編程都錯的。求給編程（最好有解釋說明）和擬合曲線，謝謝

1.jpg

回復(fù)此樓

» 猜你喜歡

0703本科鄭州大學(xué)求調(diào)劑已經(jīng)有3人回復(fù)
311求調(diào)劑已經(jīng)有4人回復(fù)
一志愿華理，數(shù)一英一285求A區(qū)調(diào)劑已經(jīng)有11人回復(fù)
調(diào)劑310 已經(jīng)有5人回復(fù)
321求調(diào)劑已經(jīng)有6人回復(fù)
求調(diào)劑已經(jīng)有7人回復(fù)
349求調(diào)劑已經(jīng)有6人回復(fù)
本科新能源科學(xué)與工程，一志愿華理能動285求調(diào)劑已經(jīng)有4人回復(fù)
材料求調(diào)劑已經(jīng)有9人回復(fù)
299求調(diào)劑已經(jīng)有7人回復(fù)

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

怎樣擬合數(shù)據(jù)求出x=f(A,B,C),y=g(A,B,C,D)中的A,B,C,D參數(shù) 已經(jīng)有18人回復(fù)
在用origin進(jìn)行非線性擬合怎樣去選擇擬合函數(shù)？已經(jīng)有6人回復(fù)
通過實(shí)驗(yàn)數(shù)據(jù)擬合，求解公式中的參數(shù) 已經(jīng)有14人回復(fù)
求助一個邏輯斯蒂函數(shù)的擬合已經(jīng)有8人回復(fù)
用matlab里的nlinfit如何擬合分段函數(shù) 已經(jīng)有8人回復(fù)
如何用得到的吸附數(shù)據(jù)進(jìn)行Langmuir-Freundlich方程擬合，而得出相關(guān)參數(shù)的具體值已經(jīng)有5人回復(fù)
如何擬合得到動力學(xué)模型參數(shù)？？已經(jīng)有3人回復(fù)
origin7.5進(jìn)行自定義函數(shù)的擬合，參數(shù)初始值怎么確定已經(jīng)有5人回復(fù)
matlab高手幫忙做個自定義復(fù)雜函數(shù)擬合，含有復(fù)數(shù)，已經(jīng)有25人回復(fù)
非線性最小二乘擬合（初值問題）lsqcurvefit 函數(shù)的設(shè)定已經(jīng)有14人回復(fù)
matlab的殘差大小跟什么因素有關(guān)？擬合的函數(shù)選擇？還是給定值的范圍？已經(jīng)有4人回復(fù)
origin進(jìn)行非線性擬合如何選擇函數(shù)或者自定義合適的函數(shù) 已經(jīng)有19人回復(fù)
matlab的簡單函數(shù)擬合已經(jīng)有6人回復(fù)
origin進(jìn)行自定義函數(shù)的擬合已經(jīng)有3人回復(fù)
【求助】如何同時擬合兩條曲線或者兩個函數(shù)！【已解決】已經(jīng)有7人回復(fù)
【求助】如何由透射光譜擬合得到參數(shù) 已經(jīng)有17人回復(fù)
【討論】怎么對數(shù)據(jù)進(jìn)行特定函數(shù)的擬合已經(jīng)有4人回復(fù)
求助高手關(guān)于origin 8.0 用自定義函數(shù)進(jìn)行擬合已經(jīng)有14人回復(fù)

1樓 2013-11-04 17:25:52

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

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

dingd

鐵桿木蟲 (職業(yè)作家)

應(yīng)助: 1641 (講師)
金幣: 15037.3
散金: 101
紅花: 234
帖子: 3410
在線: 1223.7小時
蟲號: 291104
注冊: 2006-10-28

可以試試1stOpt：

CODE:

Parameter [m,n]=[0,1];

Function y=a*t^n*(cos(pi*n/2)+b*t^m*cos(pi*(m-n)/2))/(1+2*b*t^m*cos(pi*m/2)+b^2*t^(2*m));

Data;

t=[1 2.1 5.0 10 15 20 30 40 50 60 70 80 90 100];

y=[6.0 6.08 6.2 6.4 6.47 6.58 6.8 6.85 6.9 6.93 6.95 7.0 7.05 7.1];

均方差(RMSE):0.0329522283231971
殘差平方和(SSE):0.0152018909204975
相關(guān)系數(shù)(R): 0.995826351407069
相關(guān)系數(shù)之平方(R^2): 0.991670122156715
確定系數(shù)(DC): 0.991669968535852
卡方系數(shù)(Chi-Square): 0.00112001672548189
F統(tǒng)計(jì)(F-Statistic): 400.158552543249

參數(shù) 最佳估算
---------- -------------
m 0.559176569921799
n 0.0104584004879844
a 5.88620912838722
b -0.0275911544240939

如何用已知函數(shù)去擬合實(shí)驗(yàn)數(shù)據(jù)，得出參數(shù)-1

c12.jpg

贊一下(1人)

回復(fù)此樓

10樓2013-11-04 22:23:53

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

普通回帖

zixuan203344

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 12
紅花: 1
帖子: 5
在線: 2小時
蟲號: 1829189
注冊: 2012-05-22
專業(yè): 動力學(xué)與控制

★
jjdg: 金幣+1, 感謝參與 2013-11-05 01:44:09

用擬合工具cftool吧，手頭沒裝這軟件，不能幫你算了

贊一下

回復(fù)此樓

2樓2013-11-04 18:06:57

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

yqfan

鐵蟲 (小有名氣)

應(yīng)助: 3 (幼兒園)
金幣: 164.3
紅花: 1
帖子: 85
在線: 24.7小時
蟲號: 989990
注冊: 2010-04-06
專業(yè): 橡膠及彈性體

引用回帖:

2樓: Originally posted by zixuan203344 at 2013-11-04 18:06:57
用擬合工具cftool吧，手頭沒裝這軟件，不能幫你算了

謝謝！我只有matlab，還是希望用matlab擬合出來，求助一次，弄懂了以后我就可以自己擬合了

贊一下

回復(fù)此樓

3樓2013-11-04 20:43:05

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

dingd

鐵桿木蟲 (職業(yè)作家)

應(yīng)助: 1641 (講師)
金幣: 15037.3
散金: 101
紅花: 234
帖子: 3410
在線: 1223.7小時
蟲號: 291104
注冊: 2006-10-28

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

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

cftool就是Matlab的工具箱啊。

贊一下

回復(fù)此樓

4樓2013-11-04 21:27:56

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

yqfan

鐵蟲 (小有名氣)

應(yīng)助: 3 (幼兒園)
金幣: 164.3
紅花: 1
帖子: 85
在線: 24.7小時
蟲號: 989990
注冊: 2010-04-06
專業(yè): 橡膠及彈性體

引用回帖:

4樓: Originally posted by dingd at 2013-11-04 21:27:56
cftool就是Matlab的工具箱啊。

不好意思，初學(xué)matlab才兩天完全不清楚它的強(qiáng)大功能

贊一下

回復(fù)此樓

5樓2013-11-04 21:33:00

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

yqfan

鐵蟲 (小有名氣)

應(yīng)助: 3 (幼兒園)
金幣: 164.3
紅花: 1
帖子: 85
在線: 24.7小時
蟲號: 989990
注冊: 2010-04-06
專業(yè): 橡膠及彈性體

引用回帖:

4樓: Originally posted by dingd at 2013-11-04 21:27:56
cftool就是Matlab的工具箱啊。

呵呵，剛學(xué)不懂

回復(fù)此樓

6樓2013-11-04 21:33:41

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

baobiao007

木蟲 (職業(yè)作家)

中國特色

應(yīng)助: 201 (大學(xué)生)
金幣: 6482.7
散金: 557
紅花: 40
帖子: 3050
在線: 1009.9小時
蟲號: 505962
注冊: 2008-02-18
專業(yè): 應(yīng)用地球物理學(xué)

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

★
感謝參與，應(yīng)助指數(shù) +1
jjdg: 金幣+1, 感謝參與 2013-11-05 01:44:23

現(xiàn)在關(guān)于matlab數(shù)值計(jì)算的書多如牛毛，找本好的，看看里面講曲線擬合的部分，照著弄弄不就會了嘛

贊一下

回復(fù)此樓

我同意叔本華的觀點(diǎn)，人們投身藝術(shù)和科學(xué)領(lǐng)域的強(qiáng)烈愿望之一就是逃離痛苦、殘酷和枯燥無味的現(xiàn)實(shí)生活，逃離自己飄忽不定的七情六欲的桎梏。--愛因斯坦

7樓2013-11-04 21:35:47

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

yqfan

鐵蟲 (小有名氣)

應(yīng)助: 3 (幼兒園)
金幣: 164.3
紅花: 1
帖子: 85
在線: 24.7小時
蟲號: 989990
注冊: 2010-04-06
專業(yè): 橡膠及彈性體

問題補(bǔ)充：
M-file文件：function f = Example11(c,tdata,ydata)
tdata =[1 2.1 5.0 10 15 20 30 40 50 60 70 80 90 100];
ydata =[6.0 6.08 6.2 6.4 6.47 6.58 6.8 6.85 6.9 6.93 6.95 7.0 7.05 7.1];
f=c(1)*tdata^c(4)*(cos(pi*c(4)/2)+c(2)*tdata^c(3)*cos(pi*(c(3)-c(4))/2))/(1+2*c(2)*tdata^c(3)*cos(pi*c(3)/2)+c(2)^2*tdata^(2*c(3)));

然后在命令窗口運(yùn)行
c0=[1 1 1 1];
for i = 1: 50
c=lsqnonlin('Example',c0);c=c0;
end
c
錯誤提示：
Undefined function or variable 't'.
Error in Example (line 5)
f=c(1)*t^c(4)*(cos(pi*c(4)/2)+c(2)*t^c(3)*cos(pi*(c(3)-c(4))/2))/(1+2*c(2)*t^c(3)*cos(pi*c(3)/2)+c(2)^2*t^(2*c(3)));

Error in lsqnonlin (line 197)
initVals.F = feval(funfcn{3},xCurrent,varargin{:});

Caused by:
Failure in initial user-supplied objective function evaluation. LSQNONLIN cannot
continue.

贊一下

回復(fù)此樓

8樓2013-11-04 21:46:38

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

yqfan

鐵蟲 (小有名氣)

應(yīng)助: 3 (幼兒園)
金幣: 164.3
紅花: 1
帖子: 85
在線: 24.7小時
蟲號: 989990
注冊: 2010-04-06
專業(yè): 橡膠及彈性體

引用回帖:

7樓: Originally posted by baobiao007 at 2013-11-04 21:35:47
現(xiàn)在關(guān)于matlab數(shù)值計(jì)算的書多如牛毛，找本好的，看看里面講曲線擬合的部分，照著弄弄不就會了嘛

弄了，出錯了很多次，一再修改

贊一下

回復(fù)此樓

9樓2013-11-04 21:48:19

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

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

[考研] 266分，求材料冶金能源化工等調(diào)劑 +7	哇呼哼呼哼 2026-03-27	9/450	2026-03-28 12:22 by zllcz
[考研] 317分一志愿南理工材料工程本科湖工大求調(diào)劑 +4	芋泥小鈴鐺 2026-03-28	4/200	2026-03-28 12:04 by wxiongid
[考研] 藥學(xué)105500求調(diào)劑 +3	Ssun。。 2026-03-28	3/150	2026-03-28 11:24 by lxf170613
[考研] 張芳銘-中國農(nóng)業(yè)大學(xué)-環(huán)境工程專碩-298 +4	手機(jī)用戶 2026-03-26	4/200	2026-03-28 07:17 by mmm just
[考研] 085701環(huán)境工程求調(diào)劑 +9	多久上課 2026-03-27	9/450	2026-03-28 03:58 by fmesaito
[考研] 085602 化工專碩 338分求調(diào)劑 +10	路癡小琪 2026-03-27	10/500	2026-03-28 03:36 by fmesaito
[考研] 考研化學(xué)308分求調(diào)劑 +10	你好明天你好 2026-03-23	12/600	2026-03-27 14:43 by shangxh
[考研] 化學(xué)308分求調(diào)劑 +8	你好明天你好 2026-03-23	9/450	2026-03-27 14:01 by 楊光于青云
[考研] 348求調(diào)劑 +4	小懶蟲不懶了 2026-03-27	5/250	2026-03-27 12:47 by 果果媽咪
[考研] 一志愿陜師大生物學(xué)071000，298分，求調(diào)劑 +5	SYA！ 2026-03-23	5/250	2026-03-27 09:29 by 不吃魚的貓
[考研] 333求調(diào)劑 +6	wfh030413@ 2026-03-23	6/300	2026-03-26 22:45 by 學(xué)員8dgXkO
[考研] 325求調(diào)劑 +3	Aoyijiang 2026-03-23	3/150	2026-03-26 20:46 by 不吃魚的貓
[考研] 材料與化工328分調(diào)劑 +6	。，。，。，。i 2026-03-23	6/300	2026-03-25 22:30 by 418490947
[考研] 機(jī)械學(xué)碩總分317求調(diào)劑�。。。� +4	Acaciad 2026-03-25	4/200	2026-03-25 19:59 by hanserlol
[考研] 26考研-291分-廈門大學(xué)（085601）-柔性電子學(xué)院材料工程專業(yè)求調(diào)劑 +3	min3 2026-03-24	4/200	2026-03-25 18:22 by xcjcqu
[考研] 296求調(diào)劑 +4	汪！？！ 2026-03-25	7/350	2026-03-25 16:41 by 汪�。�！
[考研] 【2026考研調(diào)劑】制藥工程 284分求相關(guān)專業(yè)調(diào)劑名額 +4	袁奐奐 2026-03-25	8/400	2026-03-25 14:32 by lbsjt
[考研] 0854人工智能方向招收調(diào)劑 +4	章小魚567 2026-03-24	4/200	2026-03-25 13:29 by 2177681040
[考研] 0703化學(xué)調(diào)劑，求導(dǎo)師收 +7	天天好運(yùn)來上岸?/a> 2026-03-24	7/350	2026-03-24 20:26 by peike
[考研] 一志愿南大，0703化學(xué)，分?jǐn)?shù)336，求調(diào)劑 +3	收到VS 2026-03-21	3/150	2026-03-21 18:42 by 學(xué)員8dgXkO