liuliping3

銀蟲 (小有名氣)

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

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

matlab里面做自定義擬合，小妹不會編碼，由于時(shí)間緊迫，想讓大家?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論文潤色/翻譯怎么收費(fèi)? 已經(jīng)有222人回復(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ù)
散金幣，求好運(yùn)，祝面上順利！已經(jīng)有68人回復(fù)
華東師范大學(xué)芯片設(shè)計(jì)徐瓏真實(shí)水平咋樣已經(jīng)有1人回復(fù)
調(diào)劑已經(jīng)有0人回復(fù)

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

MATLAB自定義公式擬合曲線已經(jīng)有3人回復(fù)
Matlab自定義擬合函數(shù)編程問題已經(jīng)有9人回復(fù)
matlab自定義沉降曲線的擬合已經(jīng)有0人回復(fù)
Matlab曲線擬合求助已經(jīng)有1人回復(fù)
MATLAB程序設(shè)計(jì)與工程應(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高手幫忙做個(gè)自定義復(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)驗(yàn): +1059
計(jì)算強(qiáng)帖: 8
應(yīng)助: 1712 (講師)
貴賓: 8.888
金幣: 68127.7
散金: 1938
紅花: 443
沙發(fā): 4
帖子: 4373
在線: 3291.4小時(shí)
蟲號: 1122189
注冊: 2010-10-14
專業(yè): 宇宙學(xué)
管轄: 計(jì)算模擬區(qū)

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

贊一下

回復(fù)此樓

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

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小時(shí)
蟲號: 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。同時(shí)畫出原圖和擬合圖形。

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小時(shí)
蟲號: 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)驗(yàn): +1059
計(jì)算強(qiáng)帖: 8
應(yīng)助: 1712 (講師)
貴賓: 8.888
金幣: 68127.7
散金: 1938
紅花: 443
沙發(fā): 4
帖子: 4373
在線: 3291.4小時(shí)
蟲號: 1122189
注冊: 2010-10-14
專業(yè): 宇宙學(xué)
管轄: 計(jì)算模擬區(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)計(jì)(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ā)帖求助！時(shí)間精力有限，恕不接受無償私信求助。

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小時(shí)
蟲號: 1294784
注冊: 2011-05-13
專業(yè): 無機(jī)材料化學(xué)

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

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

matlab不適合做這種事自定義的如果你不知道每個(gè)參數(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小時(shí)
蟲號: 2125053
注冊: 2012-11-13
性別: MM
專業(yè): 分離過程

引用回帖:

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

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

贊一下

回復(fù)此樓

我就是我

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

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

dingd

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

計(jì)算強(qiáng)帖: 4
應(yīng)助: 1641 (講師)
金幣: 15037.3
散金: 101
紅花: 234
帖子: 3410
在線: 1223.7小時(shí)
蟲號: 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è)作家)

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

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

1stOpt自動搜索，下面一個(gè)供參考：

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)計(jì)(F-Statistic): 23690.8566555606

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

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

No. 實(shí)測值 y 計(jì)算值 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小時(shí)
蟲號: 1294784
注冊: 2011-05-13
專業(yè): 無機(jī)材料化學(xué)

引用回帖:

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

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

贊一下

回復(fù)此樓

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

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

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

[考研] 生物學(xué)調(diào)劑 +3	Surekei 2026-03-21	3/150	2026-03-21 18:31 by 學(xué)員8dgXkO
[考研] 材料 271求調(diào)劑 +5	展信悅_ 2026-03-21	5/250	2026-03-21 17:29 by 學(xué)員8dgXkO
[考研] 070300化學(xué)319求調(diào)劑 +7	錦鯉0909 2026-03-17	7/350	2026-03-21 03:46 by JourneyLucky
[考研] 294求調(diào)劑材料與化工專碩 +15	陌の森林 2026-03-18	15/750	2026-03-20 23:28 by JourneyLucky
[考研] 321求調(diào)劑 +9	何潤采123 2026-03-18	11/550	2026-03-20 23:19 by JourneyLucky
[考研] 330求調(diào)劑 +4	小材化本科 2026-03-18	4/200	2026-03-20 23:13 by JourneyLucky
[考研] 287求調(diào)劑 +7	晨昏線與星海 2026-03-19	8/400	2026-03-20 22:19 by JourneyLucky
[考研] 085600材料與化工 +8	安全上岸！ 2026-03-16	8/400	2026-03-20 22:13 by luoyongfeng
[考研] 一志愿南京航空航天大學(xué)大學(xué) ，080500材料科學(xué)與工程學(xué)碩 +5	@taotao 2026-03-20	5/250	2026-03-20 20:16 by JourneyLucky
[考研] 0703化學(xué)調(diào)劑，六級已過，有科研經(jīng)歷 +13	曦熙兮 2026-03-15	13/650	2026-03-20 19:35 by Dream007008
[考研] 求調(diào)劑 +3	@taotao 2026-03-20	3/150	2026-03-20 19:35 by JourneyLucky
[考研] 工科材料085601 279求調(diào)劑 +7	困于星晨 2026-03-17	9/450	2026-03-20 17:38 by 無懈可擊111
[考研] 08工學(xué)調(diào)劑 +5	用戶573181 2026-03-20	5/250	2026-03-20 15:47 by xia_2003
[考研] 281求調(diào)劑（0805） +14	煙汐憶海 2026-03-16	25/1250	2026-03-20 15:47 by yuncha
[考博] 招收博士1-2人 +3	QGZDSYS 2026-03-18	3/150	2026-03-20 11:58 by 呱呱呱呱叫
[考研] 生物學(xué)調(diào)劑招人�。�！ +3	山海天嵐 2026-03-17	4/200	2026-03-19 21:34 by 怎么釋懷
[考研] 收復(fù)試調(diào)劑生 +4	雨后秋荷 2026-03-18	4/200	2026-03-18 14:16 by elevennnne
[考研] 考研化學(xué)學(xué)碩調(diào)劑，一志愿985 +4	張vvvv 2026-03-15	6/300	2026-03-17 17:15 by ruiyingmiao
[考研] 326求調(diào)劑 +4	諾貝爾化學(xué)獎覬?/a> 2026-03-15	7/350	2026-03-16 17:11 by 諾貝爾化學(xué)獎覬?/a>
[考研] 0856專碩279求調(diào)劑 +5	加油加油！? 2026-03-15	5/250	2026-03-15 11:58 by 2020015