| 4 | 1/1 | 返回列表 |
| 查看: 714 | 回復(fù): 3 | ||
champion3568木蟲 (職業(yè)作家)
|
[求助]
請各位大神幫我看看關(guān)于擬牛頓法Maltab代碼到底哪里有錯(cuò)??
|
|
代碼我是依照一本教程做的,可運(yùn)行時(shí)就出現(xiàn)如下錯(cuò)誤: syms x1 x2 >> x0=[1,2]'; >> [x,val,k]=sr1(fun,gfun,x0) ??? Undefined function or variable 'fun'. 這個(gè)是代碼: 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);%計(jì)算梯度 dk=-Hk*gk;%計(jì)算搜索方向 if(norm(gk)<epsilon),break;end %檢驗(yàn)終止準(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ā)表 | |
|---|---|---|---|---|
|
[考研] 287求調(diào)劑 +3 | 晨昏線與星海 2026-03-19 | 4/200 |
|
|---|---|---|---|---|
|
[考研] 311求調(diào)劑 +4 | 冬十三 2026-03-18 | 4/200 |
|
|
[考研] 一志愿武理材料305分求調(diào)劑 +5 | 想上岸的鯉魚 2026-03-18 | 6/300 |
|
|
[考研] 085600材料與化工 +5 | 安全上岸! 2026-03-16 | 5/250 |
|
|
[考研] 286求調(diào)劑 +6 | lemonzzn 2026-03-16 | 9/450 |
|
|
[考研] 266求調(diào)劑 +5 | 陽陽哇塞 2026-03-14 | 9/450 |
|
|
[考研] 材料專碩274一志愿陜西師范大學(xué)求調(diào)劑 +6 | 薛云鵬 2026-03-13 | 6/300 |
|
|
[考研] 280求調(diào)劑 +6 | 咕嚕曉曉 2026-03-18 | 7/350 |
|
|
[考研] 材料,紡織,生物(0856、0710),化學(xué)招生啦 +3 | Eember. 2026-03-17 | 9/450 |
|
|
[考研] 環(huán)境工程調(diào)劑 +8 | 大可digkids 2026-03-16 | 8/400 |
|
|
[考研] 考研求調(diào)劑 +3 | 橘頌. 2026-03-17 | 4/200 |
|
|
[考研] 268求調(diào)劑 +7 | 好運(yùn)連綿不絕 2026-03-12 | 8/400 |
|
|
[考研] 277調(diào)劑 +5 | 自由煎餅果子 2026-03-16 | 6/300 |
|
|
[考研] 278求調(diào)劑 +3 | Yy7400 2026-03-13 | 3/150 |
|
|
[考研] [導(dǎo)師推薦]西南科技大學(xué)國防/材料導(dǎo)師推薦 +3 | 尖角小荷 2026-03-16 | 6/300 |
|
|
[考研] 機(jī)械專碩325,尋找調(diào)劑院校 +3 | y9999 2026-03-15 | 5/250 |
|
|
[考研] 318求調(diào)劑 +3 | Yanyali 2026-03-15 | 3/150 |
|
|
[考研] 070300化學(xué)學(xué)碩求調(diào)劑 +6 | 太想進(jìn)步了0608 2026-03-16 | 6/300 |
|
|
[考研] 327求調(diào)劑 +6 | 拾光任染 2026-03-15 | 11/550 |
|
|
[考研] 工科材料085601 279求調(diào)劑 +8 | 困于星晨 2026-03-12 | 10/500 |
|