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

登錄注冊

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

24小時熱門版塊排行榜

返回列表

luckyman0125

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

應(yīng)助: 19 (小學(xué)生)
金幣: 12063.2
散金: 1838
紅花: 73
帖子: 4445
在線: 1211小時
蟲號: 1093343
注冊: 2010-09-08
性別: GG
專業(yè): 化學(xué)工程及工業(yè)化學(xué)

[求助] 求助，matlab/ 1stOpt 自定義微分方程擬合實驗數(shù)據(jù)來求方程中的參數(shù)

各位大神，請問如何用matlab/1stOpt 自定義微分方程來擬合實驗數(shù)據(jù)來求方程中的參數(shù)。

本人實驗數(shù)據(jù)如下
x = [0 600 1200 1800 2400 3600 5400 7200 14400 18000 43200 57600 86400 136800 151200]

y =[0.026 0.022 0.02 0.019 0.018 0.0175 0.0165 0.0158 0.015 0.0145 0.013 0.0125 0.0118 0.0105 0.0103]

需要用 dy/dx=r^2/(0.1424*y) * (588*y - 2.16*10^(-3)*cos(t)/r)擬合來求出 r 和 t。

如果用公式 dy/dx=(1.58*10^(-7))/y * (588*y - 14.4*cos(t)), 擬合 t，又是多少

謝謝

[ Last edited by luckyman0125 on 2013-6-13 at 00:35 ]

回復(fù)此樓

» 猜你喜歡

323求調(diào)劑已經(jīng)有6人回復(fù)
一志愿北京化工大學(xué) 070300 學(xué)碩 336分求調(diào)劑已經(jīng)有4人回復(fù)
352求調(diào)劑已經(jīng)有3人回復(fù)
一志愿東華大學(xué)化學(xué)070300，求調(diào)劑已經(jīng)有8人回復(fù)
277材料科學(xué)與工程080500求調(diào)劑已經(jīng)有7人回復(fù)
317求調(diào)劑已經(jīng)有18人回復(fù)
293求調(diào)劑已經(jīng)有5人回復(fù)
280分求調(diào)劑一志愿085802 已經(jīng)有7人回復(fù)
0854電子信息求調(diào)劑已經(jīng)有3人回復(fù)
263求調(diào)劑已經(jīng)有4人回復(fù)

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

求助，matlab/ 1stOpt 自定義微分方程擬合實驗數(shù)據(jù)來求方程中的參數(shù) 已經(jīng)有4人回復(fù)
origin 非線性自定義擬合問題~~方程組編寫及運行~~ 已經(jīng)有8人回復(fù)
求助，matlab自定義微分方程擬合實驗數(shù)據(jù)來求方程中的參數(shù) 已經(jīng)有4人回復(fù)
微分方程組求參數(shù)問題，求高人指教，感謝已經(jīng)有11人回復(fù)
微分方程組參數(shù)擬合的問題（多參數(shù)動力系統(tǒng)參數(shù)估計）已經(jīng)有11人回復(fù)
用1stopt擬合非線性方程結(jié)果與其他軟件擬合結(jié)果差異大已經(jīng)有5人回復(fù)
matlab高手幫忙做個自定義復(fù)雜函數(shù)擬合，含有復(fù)數(shù)，已經(jīng)有25人回復(fù)
用MATLAB或者ORIGIN怎么擬合復(fù)數(shù)數(shù)據(jù) 已經(jīng)有10人回復(fù)
matlab求解方程中的參數(shù) 已經(jīng)有21人回復(fù)
請問matlab擬合二元方程的問題已經(jīng)有3人回復(fù)
1stopt或matlab如何求解以下的非線性方程并擬合出相應(yīng)曲線？已經(jīng)有13人回復(fù)
origin8.0擬合方程中參數(shù)怎么求已經(jīng)有5人回復(fù)
一階微分方程求解，并優(yōu)化參數(shù)值已經(jīng)有14人回復(fù)
matlab 擬合反應(yīng)動力學(xué)參數(shù)結(jié)果很差。大家?guī)兔匆幌?/a> 已經(jīng)有14人回復(fù)

matlab微分方程組參數(shù)擬合，以周為單位求解，匯總后以年為單位進(jìn)行數(shù)值比較已經(jīng)有4人回復(fù)
matlab擬合方程參數(shù)時初值的選擇已經(jīng)有15人回復(fù)
學(xué)1stopt要多久？已經(jīng)有12人回復(fù)
高金求助matlab解微分方程組已經(jīng)有12人回復(fù)
求教origin擬合簡單的滲流方程（又稱逾滲方程），做復(fù)合材料的應(yīng)該會遇到吧已經(jīng)有20人回復(fù)
怎么由測試數(shù)據(jù)得到擬合方程已經(jīng)有11人回復(fù)
【求助】用matlab最優(yōu)化方法進(jìn)行參數(shù)擬合已經(jīng)有17人回復(fù)
用1stopt進(jìn)行數(shù)據(jù)擬合時，無法運行sharedmodel ,求助高手已經(jīng)有7人回復(fù)
求助Origin8.0擬合方程已經(jīng)有6人回復(fù)

有夢想才有追求

1樓 2013-06-12 06:30:09

已閱回復(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
luckyman0125: 金幣+20, ★★★很有幫助 2013-06-13 09:29:45

1stOpt求解：

1：dy/dx=r^2/(0.1424*y) * (588*y - 2.16*10^(-3)*cos(t)/r)
均方差(RMSE):0.00756626515276929
殘差平方和(SSE):0.000801477157068152
相關(guān)系數(shù)(R): 0.792802636782712
相關(guān)系數(shù)之平方(R^2): 0.628536020889622
確定系數(shù)(DC): -3.76235514275764
F統(tǒng)計(F-Statistic): 2.51976168099286

參數(shù) 最佳估算
-------------------- -------------
r -1.38058887495168E-7
t -3.14159293861944

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

文件: 數(shù)據(jù)文件-1
No x 目標(biāo) y 計算 y
1 600 0.022 0.0259516756931473
2 1200 0.02 0.0259032611459708
3 1800 0.019 0.0258547558513642
4 2400 0.018 0.0258061592974539
5 3600 0.0175 0.0257086903400129
6 5400 0.0165 0.0255617893810344
7 7200 0.0158 0.0254140384823027
8 14400 0.015 0.024814230435684
9 18000 0.0145 0.0245088170910798
10 43200 0.013 0.0222537567477591
11 57600 0.0125 0.0208558624037211
12 86400 0.0118 0.0177322184557397
13 136800 0.0105 0.0101712955998023
14 151200 0.0103 0.00656991795119879

2： dy/dx=(1.58*10^(-7))/y * (588*y - 14.4*cos(t))
周期性函數(shù)會導(dǎo)致無窮多個t，效果不好！

贊一下

回復(fù)此樓

2樓2013-06-13 09:02:58

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

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