liuliping3

銀蟲 (小有名氣)

應(yīng)助: 1 (幼兒園)
金幣: 60.8
散金: 50
帖子: 54
在線: 36.6小時
蟲號: 2125053
注冊: 2012-11-13
性別: MM
專業(yè): 分離過程

[求助] Matlab自定義擬合函數(shù) 已有3人參與

matlab里面做自定義擬合，小妹不會編碼，由于時間緊迫，想讓大家?guī)蛶兔�，不甚感�?br />
x=[0.5 1 1.5 2 2.5 3 3.5 4 4.5 5 5.5 6 6.5 7 8.7 10.4 12.1 13.8 15.5 17.2 18.9 20.6 22.3 24]
y=[0.96564 0.96454 0.95738 0.9543 0.82423 0.6315 0.51949 0.36035 0.2554 0.20826 0.14967 0.06784 0.05 0.01123 0.00749 0.000220264 0 0 0 0 0 0 0 0]

函數(shù)見附件，求取參數(shù)k值。其余數(shù)為常數(shù)

MATLAB.gif

回復(fù)此樓

» 猜你喜歡

哈爾濱理工大學(xué)物理系招收物理學(xué)考研調(diào)劑已經(jīng)有7人回復(fù)
0702一志愿吉大B區(qū)求調(diào)劑已經(jīng)有5人回復(fù)
物理學(xué)I論文潤色/翻譯怎么收費? 已經(jīng)有155人回復(fù)
求調(diào)劑已經(jīng)有0人回復(fù)
0702一志愿吉大B區(qū)求調(diào)劑有論文已經(jīng)有0人回復(fù)
請問還有沒有用Latex寫文章的小伙伴們？已經(jīng)有0人回復(fù)
光學(xué)工程學(xué)碩調(diào)劑信息已經(jīng)有31人回復(fù)
歡迎加入課題組已經(jīng)有0人回復(fù)
散金幣，求好運，祝面上順利！已經(jīng)有68人回復(fù)
華東師范大學(xué)芯片設(shè)計徐瓏真實水平咋樣已經(jīng)有1人回復(fù)
調(diào)劑已經(jīng)有0人回復(fù)

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

MATLAB自定義公式擬合曲線已經(jīng)有3人回復(fù)
Matlab自定義擬合函數(shù)編程問題已經(jīng)有9人回復(fù)
matlab自定義沉降曲線的擬合已經(jīng)有0人回復(fù)
Matlab曲線擬合求助已經(jīng)有1人回復(fù)
MATLAB程序設(shè)計與工程應(yīng)用已經(jīng)有144人回復(fù)
《MATLAB編程與工程應(yīng)用（第二版）》---Stormy Attaway 已經(jīng)有284人回復(fù)
懸賞金幣求助matlab自定義函數(shù)擬合曲線已經(jīng)有0人回復(fù)
matlab或 origin 自定義非線性擬合問題,求高手幫忙!! 已經(jīng)有1人回復(fù)
求助如何用軟件擬合自定義函數(shù) 已經(jīng)有5人回復(fù)
matlab高手幫忙做個自定義復(fù)雜函數(shù)擬合，含有復(fù)數(shù)，已經(jīng)有25人回復(fù)
origin8 自定義擬合函數(shù)設(shè)定邊界條件語句和初始化條件語句的含義已經(jīng)有2人回復(fù)
origin8 自定義擬合函數(shù)設(shè)定邊界條件語句和初始化條件語句的含義已經(jīng)有0人回復(fù)
求可以做曲面擬合的軟件已經(jīng)有10人回復(fù)
Matlab自定義擬合函數(shù) 已經(jīng)有0人回復(fù)
請問origin里面怎樣自定義歸一化的函數(shù)？已經(jīng)有3人回復(fù)
Matlab 的自定義函數(shù)擬合已經(jīng)有3人回復(fù)

我就是我

1樓 2015-11-25 14:54:48

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

月只藍(lán)

主管區(qū)長 (職業(yè)作家)

專家經(jīng)驗: +1059
計算強帖: 8
應(yīng)助: 1712 (講師)
貴賓: 8.888
金幣: 68127.7
散金: 1938
紅花: 443
沙發(fā): 4
帖子: 4373
在線: 3291.4小時
蟲號: 1122189
注冊: 2010-10-14
專業(yè): 宇宙學(xué)
管轄: 計算模擬區(qū)

公式里變量那么多，哪些是已知的，哪些是未知待擬合的，需要說明。

贊一下

回復(fù)此樓

MATLAB、MS小問題、普通問題請發(fā)帖求助！時間精力有限，恕不接受無償私信求助。

2樓2015-11-25 16:20:38

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

liuliping3

銀蟲 (小有名氣)

應(yīng)助: 1 (幼兒園)
金幣: 60.8
散金: 50
帖子: 54
在線: 36.6小時
蟲號: 2125053
注冊: 2012-11-13
性別: MM
專業(yè): 分離過程

引用回帖:

2樓: Originally posted by 月只藍(lán) at 2015-11-25 16:20:38
公式里變量那么多，哪些是已知的，哪些是未知待擬合的，需要說明。

x=[0.5 1 1.5 2 2.5 3 3.5 4 4.5 5 5.5 6 6.5 7 8.7 10.4 12.1 13.8 15.5 17.2 18.9 20.6 22.3 24];
y=[0.96564 0.96454 0.95738 0.9543 0.82423 0.6315 0.51949 0.36035 0.2554 0.20826 0.14967 0.06784 0.05 0.01123 0.00749 0.000220264 0 0 0 0 0 0 0 0];

用Matlab擬合方程的參數(shù)，方程如下：

其中：k為待求的參數(shù)，
P=101325, Y=0.0112, M=0.032, t=4/60, S=9.08, ρ=2857.14， R=8.314，T=973，u=815.29
是對應(yīng)于兩組x, y的已知值。不同P, Y, M, t, S, ρ， R，T，u下測的x, y值不同，x, y值都為已知�，F(xiàn)在希望利用這兩組數(shù)據(jù)求參數(shù)k。同時畫出原圖和擬合圖形。

1.gif

贊一下

回復(fù)此樓

我就是我

3樓2015-11-25 17:00:11

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

whyjackeyson

金蟲 (小有名氣)

應(yīng)助: 20 (小學(xué)生)
金幣: 574.1
散金: 2
紅花: 3
帖子: 83
在線: 32.3小時
蟲號: 3529482
注冊: 2014-11-10
性別: GG
專業(yè): 傳遞過程

明天晚上可以做

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

回復(fù)此樓

4樓2015-11-26 00:24:19

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

月只藍(lán)

主管區(qū)長 (職業(yè)作家)

專家經(jīng)驗: +1059
計算強帖: 8
應(yīng)助: 1712 (講師)
貴賓: 8.888
金幣: 68127.7
散金: 1938
紅花: 443
沙發(fā): 4
帖子: 4373
在線: 3291.4小時
蟲號: 1122189
注冊: 2010-10-14
專業(yè): 宇宙學(xué)
管轄: 計算模擬區(qū)

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

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感謝參與，應(yīng)助指數(shù) +1
liuliping3: 金幣+20, ★★★很有幫助, 謝謝你的幫助，還在努力裝1stOpt程序，沒用過。 2015-11-26 15:07:41

引用回帖:

3樓: Originally posted by liuliping3 at 2015-11-25 17:00:11
x=;
y=;

用Matlab擬合方程的參數(shù)，方程如下：

其中：k為待求的參數(shù)，
P=101325, Y=0.0112, M=0.032, t=4/60, S=9.08, ρ=2857.14， R=8.314，T=973，u=815.29
是對應(yīng)于兩組x, y的已知值。不同P, Y, M, ...

擬合結(jié)果不好：

相關(guān)系數(shù)(R): 0.894753969475065
相關(guān)系數(shù)之平方(R^2): 0.800584665891385
決定系數(shù)(DC): 0.528999050577239
卡方系數(shù)(Chi-Square): 2.00819821290935
F統(tǒng)計(F-Statistic): 88.3225090404399

參數(shù) 最佳估算
---------- -------------
k       6.91279108338133E-6

用1stopt軟件做的，代碼：

Constant PP=101325,YY=0.0112,M=0.032, t=4/60, S=9.08,p=2857.14, R=8.314,TT=973,u=815.29 ;
Variable x,y;
Function y=exp(k*PP*YY*M*t/S)/( exp(k*PP*YY*M*t/S)+  exp(k*p*R*TT*x/u)       );
Data;
0.5000 0.9656
1.0000 0.9645
1.5000 0.9574
2.0000 0.9543
2.5000 0.8242
3.0000 0.6315
3.5000 0.5195
4.0000 0.3604
4.5000 0.2554
5.0000 0.2083
5.5000 0.1497
6.0000 0.0678
6.5000 0.0500
7.0000 0.0112
8.7000 0.0075
10.4000 0.0002
12.1000       0
13.8000       0
15.5000       0
17.2000       0
18.9000       0
20.6000       0
22.3000       0
24.0000       0

如果要用MATLAB，程序編寫可參見：http://www.gaoyang168.com/bbs/viewthread.php?tid=6572963

贊一下(1人)

回復(fù)此樓

MATLAB、MS小問題、普通問題請發(fā)帖求助！時間精力有限，恕不接受無償私信求助。

5樓2015-11-26 14:00:54

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

nankailiang

鐵蟲 (正式寫手)

應(yīng)助: 24 (小學(xué)生)
金幣: 185.1
散金: 91
紅花: 5
帖子: 812
在線: 402.6小時
蟲號: 1294784
注冊: 2011-05-13
專業(yè): 無機材料化學(xué)

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

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

matlab不適合做這種事自定義的如果你不知道每個參數(shù)的大概范圍那很難得到滿意的擬合結(jié)果目前最好的軟件是1stopt

贊一下

回復(fù)此樓

6樓2015-11-26 21:51:45

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

liuliping3

銀蟲 (小有名氣)

應(yīng)助: 1 (幼兒園)
金幣: 60.8
散金: 50
帖子: 54
在線: 36.6小時
蟲號: 2125053
注冊: 2012-11-13
性別: MM
專業(yè): 分離過程

引用回帖:

6樓: Originally posted by nankailiang at 2015-11-26 21:51:45
matlab不適合做這種事自定義的如果你不知道每個參數(shù)的大概范圍那很難得到滿意的擬合結(jié)果目前最好的軟件是1stopt

我用1stopt擬合了，效果特別不好。只有一個未知參數(shù)k，其余都是已知的？

贊一下

回復(fù)此樓

我就是我

7樓2015-11-27 15:08:42

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

dingd

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

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

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

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

那很可能就是你公式和數(shù)據(jù)的問題了。

贊一下

回復(fù)此樓

8樓2015-11-27 17:26:55

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

dingd

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

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

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

1stOpt自動搜索，下面一個供參考：

y = 1/(p1+p2*x^2)^p3;

均方差(RMSE):2.9336044540013E-6
殘差平方和(SSE):3.4424140370145E-11
相關(guān)系數(shù)(R): 0.999989449208796
相關(guān)系數(shù)之平方(R^2): 0.99997889852891
確定系數(shù)(DC): 0.999978894811308
卡方系數(shù)(Chi-Square): 1.6179549837862E-8
F統(tǒng)計(F-Statistic): 23690.8566555606

參數(shù) 最佳估算
---------- -------------
p1 12797275.6381992
p2 -703427.377456069
p3 0.462780398685235

====== 結(jié)果輸出 ======

No. 實測值 y 計算值 y
1 0.002591 0.00259108278497855
2 0.000964 0.000960580435794965
3 0.001099 0.00110358427571721
4 0.001523 0.00152169298668189

贊一下

回復(fù)此樓

9樓2015-11-27 17:45:26

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

nankailiang

鐵蟲 (正式寫手)

應(yīng)助: 24 (小學(xué)生)
金幣: 185.1
散金: 91
紅花: 5
帖子: 812
在線: 402.6小時
蟲號: 1294784
注冊: 2011-05-13
專業(yè): 無機材料化學(xué)

引用回帖:

7樓: Originally posted by liuliping3 at 2015-11-27 15:08:42
我用1stopt擬合了，效果特別不好。只有一個未知參數(shù)k，其余都是已知的？...

原因是你的1stopt版本太低吧

贊一下

回復(fù)此樓

10樓2015-11-28 22:20:18

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

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

[考研] 0703化學(xué)調(diào)劑 +11	妮妮ninicgb 2026-03-15	15/750	2026-03-21 19:15 by ColorlessPI
[考研] 297求調(diào)劑 +3	喜歡還是不甘心 2026-03-20	3/150	2026-03-21 18:33 by 學(xué)員8dgXkO
[考研] 求調(diào)劑 +3	13341 2026-03-20	3/150	2026-03-21 18:28 by 學(xué)員8dgXkO
[考研] 311求調(diào)劑 +3	勇敢的小吳 2026-03-20	3/150	2026-03-21 17:40 by ColorlessPI
[考研] 317求調(diào)劑 +9	申子申申 2026-03-19	15/750	2026-03-21 17:31 by 學(xué)員8dgXkO
[考研] 材料 271求調(diào)劑 +5	展信悅_ 2026-03-21	5/250	2026-03-21 17:29 by 學(xué)員8dgXkO
[考研] 286分人工智能專業(yè)請求調(diào)劑愿意跨考！ +4	lemonzzn 2026-03-17	7/350	2026-03-21 16:09 by lemonzzn
[考研] 330求調(diào)劑0854 +3	assdll 2026-03-21	3/150	2026-03-21 13:01 by 搏擊518
[考研] 0856材料專碩353求調(diào)劑 +3	NIFFFfff 2026-03-20	3/150	2026-03-21 10:23 by luoyongfeng
[考研] 材料 336 求調(diào)劑 +3	An@. 2026-03-18	4/200	2026-03-21 01:39 by JourneyLucky
[考研] 311求調(diào)劑 +5	冬十三 2026-03-18	5/250	2026-03-21 00:16 by JourneyLucky
[考研] 296求調(diào)劑 +6	www_q 2026-03-18	10/500	2026-03-20 23:56 by JourneyLucky
[考研] 308求調(diào)劑 +3	阿姐阿姐家啊 2026-03-18	3/150	2026-03-20 23:24 by JourneyLucky
[考研] 一志愿武漢理工材料工程專碩調(diào)劑 +9	Doleres 2026-03-19	9/450	2026-03-20 22:36 by JourneyLucky
[考研] 085600材料與化工 +8	安全上岸！ 2026-03-16	8/400	2026-03-20 22:13 by luoyongfeng
[考研] 一志愿西安交通大學(xué) 學(xué)碩 354求調(diào)劑211或者雙一流 +3	我想要讀研究生 2026-03-20	3/150	2026-03-20 20:13 by JourneyLucky
[考研] 招收調(diào)劑碩士 +4	lidianxing 2026-03-19	12/600	2026-03-20 12:25 by lidianxing
[考研] 288求調(diào)劑，一志愿華南理工大學(xué)071005 +5	ioodiiij 2026-03-17	5/250	2026-03-19 18:22 by zcl123
[考研] 收復(fù)試調(diào)劑生 +4	雨后秋荷 2026-03-18	4/200	2026-03-18 14:16 by elevennnne
[考研] 材料專碩326求調(diào)劑 +6	墨煜姒莘 2026-03-15	7/350	2026-03-17 17:10 by ruiyingmiao