| 4 | 1/1 | 返回列表 |
| 查看: 1534 | 回復(fù): 3 | ||
rita128新蟲 (初入文壇)
|
[求助]
1stopt做最小二乘法回歸參數(shù),請(qǐng)求大家?guī)兔χ笇?dǎo)下
|
|
第一次用這個(gè)做最小二乘法回歸參數(shù)。一共六個(gè)參數(shù)。有八組數(shù)據(jù),目標(biāo)方程一個(gè)。應(yīng)該是最簡單的那種吧。代碼如下,運(yùn)行不出來,不知道怎么定義常量或者輸入數(shù)據(jù)。求助懂的人幫忙指點(diǎn)一下。謝謝 Title "Type your title here"i; Parameters x1,x2,x3,x4,x5,x6 ; Constant a(1:8)= [0.593492611, 0.615383905,0.605448133,0.692841215,0.655361944,0.757020028,0.671381587,0.735878492]; Constant b(1:8)=[ 5.1612, 4.6650 ,3.7857,6.7947,2.3167 ,1.4693,7.0426 ,10.0069]; Conststr I=3*a; Conststr f1=-0.76*(I^0.5/(1+1.2*I^0.5)+2/1.2*ln(1+1.2*I^0.5));//f'(i) Conststr B1=-0.00804+1.4354/(2*I)*(1-(1+2*I^0.5)*EXP((-2)*I^0.5));//B Conststr B2=1.4354/(2*I^2)*(-1+(1+2*I^0.5+2*I)*exp((-2)*I^0.5));//B' Conststr F1=1/(2*I)*(1-(1+2*I^0.5)*exp((-2)*I^0.5));//F(i) Conststr F2=1/(2*I^2)*(-1+(1+2*I^0.5+2*I)*EXP((-2)*I^0.5));//F'()i MinFunction F=sum(i=1:8)(((f1+8/3*a*(B1+2*a*0.001838)+2*a^2*(2*B2+4/3*0.001838)+b*x3/3+4/3*b*a*(x1+F1*x2)+8/3*a*b*(x1+F1*x2)+4*a^2*b*F2*x2)+ln(a)+1.5*ln(b)+1.5*(a*x3+4*a^2*(x1+F1*x2)+2*b*x4+3*b^2*x5)-lnx6)^2); |
鐵桿木蟲 (職業(yè)作家)
|
幾點(diǎn)錯(cuò)誤: 1:“l(fā)nx6”應(yīng)該是“l(fā)n(x6)”吧 2:1stOpt不區(qū)分大小寫 3:“Sum”的寫法有問題 正確代碼如下: 結(jié)果: 目標(biāo)函數(shù)值(最小): 0.000360743527320632 x1: -0.140392701334565 x2: 0.216320938632417 x3: 2.29529578782098 x4: -0.370075186018573 x5: 0.00657741865607194 x6: 1.24296439011043 |
新蟲 (初入文壇)
鐵桿木蟲 (職業(yè)作家)
| 4 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 求材料調(diào)劑 +8 | 隔壁陳先生 2026-03-12 | 8/400 |
|
|---|---|---|---|---|
|
[考研] 321求調(diào)劑 +3 | 何潤采123 2026-03-18 | 3/150 |
|
|
[考研] 274求調(diào)劑 +5 | S.H1 2026-03-18 | 5/250 |
|
|
[考研] 281求調(diào)劑(0805) +6 | 煙汐憶海 2026-03-16 | 15/750 |
|
|
[考研] 材料專碩英一數(shù)二306 +4 | z1z2z3879 2026-03-18 | 4/200 |
|
|
[考研] 0703化學(xué)調(diào)劑,求各位老師收留 +9 | 秋有木北 2026-03-14 | 9/450 |
|
|
[考研] 【同濟(jì)軟件】軟件(085405)考研求調(diào)劑 +3 | 2026eternal 2026-03-18 | 3/150 |
|
|
[考研] 085601材料工程專碩求調(diào)劑 +6 | 慕寒mio 2026-03-16 | 6/300 |
|
|
[考研] 304求調(diào)劑 +12 | 小熊joy 2026-03-14 | 13/650 |
|
|
[考研]
|
胡辣湯放糖 2026-03-15 | 6/300 |
|
|
[考研] 307求調(diào)劑 +3 | 冷笙123 2026-03-17 | 3/150 |
|
|
[考研] 334求調(diào)劑 +3 | 志存高遠(yuǎn)意在機(jī)?/a> 2026-03-16 | 3/150 |
|
|
[考研] 26考研求調(diào)劑 +6 | 丶宏Sir 2026-03-13 | 6/300 |
|
|
[考研] 一志愿蘇州大學(xué)材料工程(085601)專碩有科研經(jīng)歷三項(xiàng)國獎(jiǎng)兩個(gè)實(shí)用型專利一項(xiàng)省級(jí)立項(xiàng) +6 | 大火山小火山 2026-03-16 | 8/400 |
|
|
[考研] 285化工學(xué)碩求調(diào)劑(081700) +9 | 柴郡貓_ 2026-03-12 | 9/450 |
|
|
[考研] 304求調(diào)劑 +5 | 素年祭語 2026-03-15 | 5/250 |
|
|
[考研] 0703化學(xué)調(diào)劑 290分有科研經(jīng)歷,論文在投 +7 | 膩膩gk 2026-03-14 | 7/350 |
|
|
[考研] 0703化學(xué)求調(diào)劑 +7 | 綠豆芹菜湯 2026-03-12 | 7/350 |
|
|
[考研] 工科材料085601 279求調(diào)劑 +8 | 困于星晨 2026-03-12 | 10/500 |
|
|
[考博] 26讀博 +4 | Rui135246 2026-03-12 | 10/500 |
|