| 5 | 1/1 | 返回列表 |
| 查看: 1544 | 回復: 3 | ||
| 當前只顯示滿足指定條件的回帖,點擊這里查看本話題的所有回帖 | ||
rita128新蟲 (初入文壇)
|
[求助]
1stopt做最小二乘法回歸參數(shù),請求大家?guī)兔χ笇?
|
|
|
第一次用這個做最小二乘法回歸參數(shù)。一共六個參數(shù)。有八組數(shù)據(jù),目標方程一個。應該是最簡單的那種吧。代碼如下,運行不出來,不知道怎么定義常量或者輸入數(shù)據(jù)。求助懂的人幫忙指點一下。謝謝 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è)作家)
鐵桿木蟲 (職業(yè)作家)
|
幾點錯誤: 1:“l(fā)nx6”應該是“l(fā)n(x6)”吧 2:1stOpt不區(qū)分大小寫 3:“Sum”的寫法有問題 正確代碼如下: 結(jié)果: 目標函數(shù)值(最小): 0.000360743527320632 x1: -0.140392701334565 x2: 0.216320938632417 x3: 2.29529578782098 x4: -0.370075186018573 x5: 0.00657741865607194 x6: 1.24296439011043 |
新蟲 (初入文壇)
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 341求調(diào)劑(一志愿湖南大學070300) +5 | 番茄頭--- 2026-03-22 | 6/300 |
|
|---|---|---|---|---|
|
[考研] 361求調(diào)劑 +3 | Glack 2026-03-22 | 3/150 |
|
|
[考研] 070300化學求調(diào)劑 +8 | 苑豆豆 2026-03-20 | 8/400 |
|
|
[考研] 一志愿重慶大學085700資源與環(huán)境,總分308求調(diào)劑 +7 | 墨墨漠 2026-03-23 | 8/400 |
|
|
[考研] 328求調(diào)劑,英語六級551,有科研經(jīng)歷 +7 | 生物工程調(diào)劑 2026-03-17 | 12/600 |
|
|
[考研] 上海電力大學材料防護與新材料重點實驗室招收調(diào)劑研究生(材料、化學、電化學,環(huán)境) +3 | 我愛學電池 2026-03-23 | 3/150 |
|
|
[論文投稿] 急發(fā)核心期刊論文 +3 | 賢達問津 2026-03-23 | 5/250 |
|
|
[考研] 一志愿華中農(nóng)業(yè)071010,總分320求調(diào)劑 +5 | 困困困困坤坤 2026-03-20 | 6/300 |
|
|
[考研] 289求調(diào)劑 +7 | 懷瑾握瑜l 2026-03-20 | 7/350 |
|
|
[考研] 354求調(diào)劑 +7 | Tyoumou 2026-03-18 | 10/500 |
|
|
[考研] 求調(diào)劑 +4 | 要好好無聊 2026-03-21 | 4/200 |
|
|
[考研] 297求調(diào)劑 +3 | 喜歡還是不甘心 2026-03-20 | 3/150 |
|
|
[考研] 材料學碩333求調(diào)劑 +3 | 北道巷 2026-03-18 | 3/150 |
|
|
[考研] 二本跨考鄭大材料306英一數(shù)二 +3 | z1z2z3879 2026-03-17 | 3/150 |
|
|
[考研] 295求調(diào)劑 +4 | 一志愿京區(qū)211 2026-03-18 | 6/300 |
|
|
[考研] 317求調(diào)劑 +5 | 申子申申 2026-03-19 | 9/450 |
|
|
[考研] 求調(diào)劑一志愿南京航空航天大學289分 +3 | @taotao 2026-03-19 | 3/150 |
|
|
[考研]
|
不想起名字112 2026-03-19 | 3/150 |
|
|
[考研] 085600材料與化工求調(diào)劑 +6 | 緒幸與子 2026-03-17 | 6/300 |
|
|
[考研] 考研求調(diào)劑 +3 | 橘頌. 2026-03-17 | 4/200 |
|