| 4 | 1/1 | 返回列表 |
| 查看: 716 | 回復(fù): 3 | ||
champion3568木蟲 (職業(yè)作家)
|
[求助]
請各位大神幫我看看關(guān)于擬牛頓法Maltab代碼到底哪里有錯??
|
|
代碼我是依照一本教程做的,可運(yùn)行時就出現(xiàn)如下錯誤: syms x1 x2 >> x0=[1,2]'; >> [x,val,k]=sr1(fun,gfun,x0) ??? Undefined function or variable 'fun'. 這個是代碼: function [x,val,k] =sr1(fun,gfun,x0) %功能:用對稱秩1算法求解無約束問題:min f(x) %輸入:X0是初始點(diǎn),fun,gfun分別是目標(biāo)函數(shù)及其梯度 %輸出:x,val分別是近似最優(yōu)點(diǎn)和最優(yōu)值,k是迭代次數(shù). maxk=500; %給出最大迭代次數(shù) rho=0.55;sigma=0.4;epsilon=le-4; k=0; n=length(x0);Hk=eye(n); while(k<maxk) gk=feval(gfun,x0);%計算梯度 dk=-Hk*gk;%計算搜索方向 if(norm(gk)<epsilon),break;end %檢驗終止準(zhǔn)則 m=0;mk=0; while(m<20) %用Armijo搜索求步長 if(feval(fun,x0+rho^m*dk)<feval(fun,x0)+sigma*rho^m*gk'*dk) mk=m;break; end m=m+1; end x=x0+rho^mk*dk; sk=x-x0; yk=feval(gfun,x)-gk; Hk=Hk+(sk-Hk*yk)*(sk-Hk*yk)'/((sk-Hk*yk)'*yk); %秩1校正 k=k+1; x0=x; end val=feval(fun,x0); |
新蟲 (正式寫手)
木蟲 (小有名氣)
木蟲 (職業(yè)作家)
| 4 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 307求調(diào)劑 +6 | 冷笙123 2026-03-17 | 6/300 |
|
|---|---|---|---|---|
|
[考研] 085600材料與化工調(diào)劑 324分 +10 | llllkkkhh 2026-03-18 | 12/600 |
|
|
[考研] 求調(diào)劑,一志愿:南京航空航天大學(xué)大學(xué) ,080500材料科學(xué)與工程學(xué)碩,總分289分 +3 | @taotao 2026-03-19 | 3/150 |
|
|
[考研] 294求調(diào)劑材料與化工專碩 +11 | 陌の森林 2026-03-18 | 11/550 |
|
|
[考研] 一志愿中海洋材料工程專碩330分求調(diào)劑 +7 | 小材化本科 2026-03-18 | 7/350 |
|
|
[考研] 一志愿吉林大學(xué)材料學(xué)碩321求調(diào)劑 +6 | Ymlll 2026-03-18 | 9/450 |
|
|
[考研] 085700資源與環(huán)境308求調(diào)劑 +3 | 墨墨漠 2026-03-18 | 3/150 |
|
|
[考研] 303求調(diào)劑 +4 | 睿08 2026-03-17 | 6/300 |
|
|
[考研] 278求調(diào)劑 +5 | 煙火先于春 2026-03-17 | 5/250 |
|
|
[考研] 296求調(diào)劑 +5 | 大口吃飯 身體健 2026-03-13 | 5/250 |
|
|
[考研] 332求調(diào)劑 +6 | Zz版 2026-03-13 | 6/300 |
|
|
[考研] 材料與化工專碩調(diào)劑 +5 | heming3743 2026-03-16 | 5/250 |
|
|
[考研] 274求調(diào)劑 +5 | 時間點(diǎn) 2026-03-13 | 5/250 |
|
|
[考研]
|
zhouzhen654 2026-03-16 | 3/150 |
|
|
[考研] 本科南京大學(xué)一志愿川大藥學(xué)327 +3 | 麥田耕者 2026-03-14 | 3/150 |
|
|
[考研] 復(fù)試調(diào)劑 +4 | z1z2z3879 2026-03-14 | 5/250 |
|
|
[考研] 一志愿哈工大材料324分求調(diào)劑 +5 | 閆旭東 2026-03-14 | 5/250 |
|
|
[考研] 304求調(diào)劑 +7 | 7712b 2026-03-13 | 7/350 |
|
|
[考研] 308求調(diào)劑 +3 | 是Lupa啊 2026-03-12 | 3/150 |
|
|
[考研] 一志愿山大07化學(xué) 332分 四六級已過 本科山東雙非 求調(diào)劑! +3 | 不想理你 2026-03-12 | 3/150 |
|