| 7 | 1/1 | 返回列表 |
| 查看: 1138 | 回復(fù): 6 | |||
[求助]
MATLAB非線性擬合 已有1人參與
|
|
數(shù)據(jù)是每一個x對應(yīng)一個y x=10 20 30 40 50 60 70 80 90 100 110 120 130 140 150 160 170 180 190 200 y=417.308 390.957 359.313 347.39 315.77 309.942 290.582 275.24 251.62 200.9 180.76 134.65 149.45 108.823 138.25 103.431 93.75 104.87 83.871 108.475 擬合要用的方程 y=(700+x)*(2*x/(700+x))*(n1*h1*k1*L/(1+k1*L)+n2*k2*h2*L/(1+k2*L)),其中L很復(fù)雜,用下面的方程單獨(dú)表示出來: L=2*(-1*Y)^0.5*cos((θ)/3)-(B)/3/(A); θ=arccos((-1*(x))/(((-1*(y))^3)^0.5)); A=k1*k2; B=k1+k2+k1*k2*(2*x/(700+x)*(n1+n2)-7000/(700+x)); C=1+2*x/(700+x)*(n1*k1+n2*k2)-7000/(700+x)*(k1+k2); D=-1*7000/(700+X); Y=C/3/(A)-((B)/3/(A))^2; X=((B)/3/(A))^3-(B)*(C)/6/((A)^2)+(D)/2/(A);注意這里的X,Y是大寫的 自變量、因變量是小寫的x,y n1、n2、K1、K2、H1、H2是需要擬合的參數(shù)。 請大牛幫我把matlab的程序?qū)懗鰜恚救耸浅鯇W(xué)者,僅僅看了matlab的教程而已,求指教。 |

主管區(qū)長 (職業(yè)作家)
![]() |
專家經(jīng)驗(yàn): +1059 |

捐助貴賓 (小有名氣)
|
你的問題可參照下面的帖子進(jìn)行,用lsqnonlin函數(shù)。 http://www.gaoyang168.com/html/201407/7603645.html |
木蟲 (正式寫手)

鐵桿木蟲 (著名寫手)
ส็็็

鐵桿木蟲 (著名寫手)
ส็็็
|
這問題實(shí)際是 參數(shù)個數(shù)(6)維 的非線性優(yōu)化,所以,在參數(shù)的較為準(zhǔn)確的初值無法知道的情況下,用普通迭代方法lsqcurvefit, lsqnonlin之類都不太可能給出好結(jié)果。 從實(shí)際使用看,1stOpt的differential evolution的確是不錯的選項(xiàng)。 用matlab也可以,但是要很擅長隨機(jī)優(yōu)化才行。自帶的ParticalSwarmOptimization函數(shù)優(yōu)化能力很菜。Mathematica內(nèi)置了Differential Evolution,但是某些情況下(很可能大多數(shù)情況下)不及1stOpt |


| 7 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 0703化學(xué)調(diào)劑 ,六級已過,有科研經(jīng)歷 +11 | 曦熙兮 2026-03-15 | 11/550 |
|
|---|---|---|---|---|
|
[考研] 286求調(diào)劑 +6 | lemonzzn 2026-03-16 | 10/500 |
|
|
[考研] 085600材料與化工求調(diào)劑 +6 | 緒幸與子 2026-03-17 | 6/300 |
|
|
[考研] 346求調(diào)劑[0856] +3 | WayneLim327 2026-03-16 | 6/300 |
|
|
[考研] 材料專業(yè)求調(diào)劑 +5 | hanamiko 2026-03-18 | 5/250 |
|
|
[考研] 344求調(diào)劑 +6 | knight344 2026-03-16 | 7/350 |
|
|
[考研] 26調(diào)劑/材料/英一數(shù)二/總分289/已過A區(qū)線 +7 | 步川酷紫123 2026-03-13 | 7/350 |
|
|
[考研] 297求調(diào)劑 +8 | 戲精丹丹丹 2026-03-17 | 8/400 |
|
|
[考研] 280求調(diào)劑 +6 | 咕嚕曉曉 2026-03-18 | 7/350 |
|
|
[考研] 環(huán)境工程調(diào)劑 +8 | 大可digkids 2026-03-16 | 8/400 |
|
|
[考研] 268求調(diào)劑 +6 | 簡單點(diǎn)0 2026-03-17 | 6/300 |
|
|
[論文投稿] 有沒有大佬發(fā)小論文能帶我個二作 +3 | 增銳漏人 2026-03-17 | 4/200 |
|
|
[考研] 考研調(diào)劑 +3 | 淇ya_~ 2026-03-17 | 5/250 |
|
|
[考研] 東南大學(xué)364求調(diào)劑 +5 | JasonYuiui 2026-03-15 | 5/250 |
|
|
[考研] 機(jī)械專碩325,尋找調(diào)劑院校 +3 | y9999 2026-03-15 | 5/250 |
|
|
[考研] 0854控制工程 359求調(diào)劑 可跨專業(yè) +3 | 626776879 2026-03-14 | 9/450 |
|
|
[考研] 一志愿華中師范071000,325求調(diào)劑 +6 | RuitingC 2026-03-12 | 6/300 |
|
|
[考研] 0703 物理化學(xué)調(diào)劑 +3 | 我可以上岸的對?/a> 2026-03-13 | 5/250 |
|
|
[考研] 0703化學(xué)調(diào)劑 290分有科研經(jīng)歷,論文在投 +7 | 膩膩gk 2026-03-14 | 7/350 |
|
|
[考研] 22408總分284求調(diào)劑 +3 | InAspic 2026-03-13 | 3/150 |
|