liuliping3

銀蟲 (小有名氣)

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

[求助] Matlab自定義擬合函數(shù)編程問題已有2人參與

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];

用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ù)此樓

» 猜你喜歡

288求調(diào)劑，一志愿華南理工大學(xué)071005 已經(jīng)有5人回復(fù)
一志愿中國海洋大學(xué)，生物學(xué)，301分，求調(diào)劑已經(jīng)有5人回復(fù)
294求調(diào)劑材料與化工專碩已經(jīng)有12人回復(fù)
0856調(diào)劑，是學(xué)校就去已經(jīng)有5人回復(fù)
0817 化學(xué)工程 299分求調(diào)劑有科研經(jīng)歷有二區(qū)文章已經(jīng)有13人回復(fù)
梁成偉老師課題組歡迎你的加入已經(jīng)有11人回復(fù)
復(fù)試調(diào)劑已經(jīng)有6人回復(fù)
0703化學(xué)調(diào)劑，六級已過，有科研經(jīng)歷已經(jīng)有11人回復(fù)
321求調(diào)劑已經(jīng)有10人回復(fù)
東華理工大學(xué)化材專業(yè)26屆碩士博士申請已經(jīng)有8人回復(fù)

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

Origin自定義模擬與matlab中cftool工具箱擬合結(jié)果比較已經(jīng)有4人回復(fù)
通過Matlab軟件自行計算了一個自由曲面擬合成非球面多項式或者zernike多項式的問題已經(jīng)有3人回復(fù)
關(guān)于使用origin數(shù)據(jù)擬合的問題請教。已經(jīng)有6人回復(fù)
【求助】請各位高手告訴我：這種想法能不能用matlab實現(xiàn)？200金幣求助已經(jīng)有16人回復(fù)

我就是我

1樓 2015-11-25 19:35:56

已閱回復(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

兩組值在哪兒？公式最好先給出文本格式

贊一下

回復(fù)此樓

2樓2015-11-25 21:00:47

已閱回復(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 dingd at 2015-11-25 21:00:47
兩組值在哪兒？公式最好先給出文本格式

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];

y=exp(k*P*Y*M*t/S)/(exp(k*P*Y*M*t/S)+exp(k*ρ*R*T*x/u)-1);

P=101325, Y=0.0112, M=0.032, t=4/60, S=9.08, ρ=2857.14， R=8.314，T=973，u=815.29

贊一下

回復(fù)此樓

我就是我

3樓2015-11-25 21:56:49

已閱回復(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 dingd at 2015-11-25 21:00:47
兩組值在哪兒？公式最好先給出文本格式

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];

y= exp(k*P*Y*M*t/S)/( exp(k*P*Y*M*t/S)+exp(k*ρ*R*T*x/u)-1);

P=101325, Y=0.0112, M=0.032, t=4/60, S=9.08, ρ=2857.14， R=8.314，T=973，u=815.29

贊一下

回復(fù)此樓

我就是我

4樓2015-11-25 21:59:44

已閱回復(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)助回帖

不是兩組x, y數(shù)值嗎？對應(yīng)不同的P, Y, M, t, S, ρ， R，T，u值

如果只有這一組，1stOpt試試：

CODE:

Constant P=101325, Y1=0.0112, M=0.032, t1=4/60, S=9.08, p1=2857.14, R=8.314,T=973,u=815.29;

Function y= exp(k*P*Y1*M*t1/S)/( exp(k*P*Y1*M*t1/S)+exp(k*p1*R*T*x/u)-1);

Data;

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];

均方差(RMSE):0.13879469945569
殘差平方和(SSE):0.462335246327889
相關(guān)系數(shù)(R): 0.958172545145573
相關(guān)系數(shù)之平方(R^2): 0.918094626270745
確定系數(shù)(DC): 0.860178141633635
卡方系數(shù)(Chi-Square): 1.00096945974885

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

贊一下(1人)

回復(fù)此樓

5樓2015-11-25 22:46: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è): 分離過程

引用回帖:

5樓: Originally posted by dingd at 2015-11-25 22:46:38
不是兩組x, y數(shù)值嗎？對應(yīng)不同的P, Y, M, t, S, ρ， R，T，u值

如果只有這一組，1stOpt試試：

Constant P=101325, Y1=0.0112, M=0.032, t1=4/60, S=9.08, p1=2857.14, R=8.314,T=973,u=815.29;
Function y= ...

謝謝你的回復(fù)�？赡芪冶磉_(dá)的不清楚，我現(xiàn)在差不多有22組x,y值，P，R，p1是始終不變的，其余的都變，我本來想的是讓大家?guī)臀揖幰粋€只有一組x,y的程序，然后計算的時候我自己改變數(shù)。我沒學(xué)過MATLAB，不懂，所以是可以將全部數(shù)據(jù)都上傳嗎？另外1stOpt也沒用過，沒有相關(guān)的軟件。

贊一下

回復(fù)此樓

我就是我

6樓2015-11-26 09:09:52

已閱回復(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 liuliping3 at 2015-11-26 09:09:52
謝謝你的回復(fù)。可能我表達(dá)的不清楚，我現(xiàn)在差不多有22組x,y值，P，R，p1是始終不變的，其余的都變，我本來想的是讓大家?guī)臀揖幰粋€只有一組x,y的程序，然后計算的時候我自己改變數(shù)。我沒學(xué)過MATLAB，不懂，所以是可 ...

再麻煩您一下，那個1stOpt軟件擬合的時候可以畫出原圖和擬合圖嗎？另外您有1stOpt的程序沒？

贊一下

回復(fù)此樓

我就是我

7樓2015-11-26 09:23:17

已閱回復(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)助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
liuliping3: 金幣+10, ★有幫助, 這個軟件擬合的不準(zhǔn)確，偏差特別大，不過還是謝謝了。 2015-11-26 17:55:44

自己用1stOpt試試吧，完全可以滿足的原問題的要求。至于軟件自己解決。

贊一下

回復(fù)此樓

8樓2015-11-26 11:12:15

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

享受記憶

新蟲 (小有名氣)

應(yīng)助: 2 (幼兒園)
金幣: 1579.2
帖子: 181
在線: 9.9小時
蟲號: 3742115
注冊: 2015-03-16
專業(yè): 人工智能與知識工程

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

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

什么時候用，這一兩天有空可以幫你編編

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

贊一下

回復(fù)此樓

9樓2015-11-26 14:28:57

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

liuliping3

銀蟲 (小有名氣)

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

引用回帖:

9樓: Originally posted by 享受記憶 at 2015-11-26 14:28:57
什么時候用，這一兩天有空可以幫你編編

恩，有點著急想用，麻煩您有空幫忙編輯一下了，非常感謝。

贊一下

回復(fù)此樓

我就是我

10樓2015-11-26 15:01:50

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

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

[考研] 0817 化學(xué)工程 299分求調(diào)劑有科研經(jīng)歷有二區(qū)文章 +13	rare12345 2026-03-18	13/650	2026-03-19 17:30 by jean5056
[考研] 286求調(diào)劑 +6	lemonzzn 2026-03-16	10/500	2026-03-19 14:31 by lemonzzn
[考研] 一志愿西北大學(xué) ，070300化學(xué)學(xué)碩，總分287，雙非一本，求調(diào)劑。 +3	晨昏線與星海 2026-03-19	3/150	2026-03-19 13:36 by houyaoxu
[考研] 287求調(diào)劑 +3	晨昏線與星海 2026-03-19	4/200	2026-03-19 12:32 by peike
[考研] 材料080500調(diào)劑求收留 +4	一顆meteor 2026-03-13	4/200	2026-03-19 10:32 by 30660438
[考研] 0854，計算機類招收調(diào)劑 +3	胡辣湯放糖 2026-03-15	6/300	2026-03-18 12:09 by 上岸上岸……..
[考博] 26博士申請 +3	1042136743 2026-03-17	3/150	2026-03-17 23:30 by 輕松不少隨
[基金申請] 被我言中：新模板不強調(diào)格式了，假專家開始管格式了 +4	beefly 2026-03-14	4/200	2026-03-17 22:04 by 黃鳥于飛Chao
[考研] 考研求調(diào)劑 +3	橘頌. 2026-03-17	4/200	2026-03-17 21:43 by 有只貍奴
[考研] 一志愿，福州大學(xué)材料專碩339分求調(diào)劑 +3	木子momo青爭 2026-03-15	3/150	2026-03-17 07:52 by laoshidan
[考研] 藥學(xué)383 求調(diào)劑 +3	藥學(xué)chy 2026-03-15	4/200	2026-03-16 20:51 by 元子^0^
[考研] 070300化學(xué)學(xué)碩求調(diào)劑 +6	太想進(jìn)步了0608 2026-03-16	6/300	2026-03-16 16:13 by kykm678
[考研] 085600材料與化工求調(diào)劑 +13	enenenhui 2026-03-13	14/700	2026-03-16 15:19 by 了了了了。。
[考研] 0856專碩279求調(diào)劑 +5	加油加油！? 2026-03-15	5/250	2026-03-15 11:58 by 2020015
[考研] 22408總分284求調(diào)劑 +3	InAspic 2026-03-13	3/150	2026-03-15 11:10 by zhq0425
[考研] 289求調(diào)劑 +4	這么名字咋樣 2026-03-14	6/300	2026-03-14 18:58 by userper
[考研] 297一志愿上交085600求調(diào)劑 +5	指尖八千里 2026-03-14	5/250	2026-03-14 17:26 by a不易
[基金申請] 現(xiàn)在如何回避去年的某一個專家，不知道名字 +3	zk200107 2026-03-12	6/300	2026-03-14 17:13 by zk200107
[考研] 333求調(diào)劑 +3	152697 2026-03-12	4/200	2026-03-13 07:08 by Iveryant
[考研] 081200-11408-276學(xué)碩求調(diào)劑 +3	崔wj 2026-03-12	4/200	2026-03-12 19:33 by 求調(diào)劑zz