當(dāng)前只顯示滿(mǎn)足指定條件的回帖，點(diǎn)擊這里查看本話(huà)題的所有回帖

guan0603

應(yīng)助: 0 (幼兒園)
金幣: 67.1
帖子: 8
在線(xiàn): 9.8小時(shí)
蟲(chóng)號(hào): 3599995
注冊(cè): 2014-12-17
專(zhuān)業(yè): 運(yùn)籌學(xué)

[求助] 哪位MATLAB高手能幫忙調(diào)試一下程序，不是遇到這問(wèn)題就是那問(wèn)題已有1人參與

function [x,val,k]=frcg(fun,gfun,x0)
% 功能: 用FR共軛梯度法求解無(wú)約束問(wèn)題: min f(x)
%輸入: x0是初始點(diǎn), fun, gfun分別是目標(biāo)函數(shù)和梯度
%function g=gfun(x) g=[400*x(1)*(x(1)^2-x(2))+2*(x(1)-1),-200*(x(1)^2-x(2))]';
%function f=fun(x) f=100*(x(1)^2-x(2))^2+(x(1)-1)^2;
%輸出: x, val分別是近似最優(yōu)點(diǎn)和最優(yōu)值, k是迭代次數(shù).
maxk=5000; %最大迭代次數(shù)
rho=0.6;sigma=0.4;
k=0;epsilon=1e-4;
n=length(x0);
while(k<maxk)
g=feval(gfun,x0); %計(jì)算梯du
itern=k-(n+1)*floor(k/(n+1));
itern=itern+1;
%計(jì)算搜索方向
if(itern==1)
d=-g;
else
beta=(g'*g)/(g0'*g0);
d=-g+beta*d0; gd=g'*d;
if(gd>=0.0)
d=-g;
end
end
if(norm(g)<epsilon),break;end %檢驗(yàn)終止條件
m=0;mk=0;
while(m<20) %Armijo搜索
if(feval(fun,x0+rho^m*d)<feval(fun,x0)+sigma*rho^m*g'*d)
mk=m;break;
end
m=m+1;
end
x0=x0+rho^mk*d;
val=feval(fun,x0);
g0=g;d0=d;
k=k+1;
end
x=x0;
val=feval(fun,x);
%輸入x0=[-1.2 1]';
%[x,val,k]=frcg('fun','gfun',x0)
在命令窗口輸入
x0=[-1.2 1]';
>> [x,val,k]=frcg('fun','gfun',x0)
??? Undefined command/function 'frcg'.
怎么回事呢

回復(fù)此樓

» 猜你喜歡

303求調(diào)劑已經(jīng)有5人回復(fù)
一志愿武理085500機(jī)械專(zhuān)業(yè)總分300求調(diào)劑已經(jīng)有7人回復(fù)
材料學(xué)求調(diào)劑已經(jīng)有5人回復(fù)
考研調(diào)劑已經(jīng)有4人回復(fù)
281求調(diào)劑已經(jīng)有4人回復(fù)
0805 316求調(diào)劑已經(jīng)有6人回復(fù)
085601求調(diào)劑總分293英一數(shù)二已經(jīng)有3人回復(fù)
08工學(xué)調(diào)劑已經(jīng)有17人回復(fù)
340求調(diào)劑已經(jīng)有4人回復(fù)
311求調(diào)劑已經(jīng)有3人回復(fù)

» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:

高手指點(diǎn)：如何使用matlab使數(shù)據(jù)行與行間空一行已經(jīng)有6人回復(fù)
matlab編程 ODE45 已經(jīng)有4人回復(fù)
linux 下 matlab 的調(diào)用問(wèn)題已經(jīng)有3人回復(fù)
matlab畫(huà)圖問(wèn)題: 怎樣用小圖標(biāo)或自定義圖形代替+ x . o等做圖已經(jīng)有10人回復(fù)
一個(gè)關(guān)于matlab的simulink的變頻空調(diào)PI控制的問(wèn)題已經(jīng)有9人回復(fù)
請(qǐng)大俠指教matlab的這個(gè)錯(cuò)誤提示到底是個(gè)神馬意思？已經(jīng)有12人回復(fù)
有個(gè)matlab程序，但是運(yùn)行不下來(lái)，求高手已經(jīng)有23人回復(fù)
matlab程序中調(diào)用m文件的問(wèn)題已經(jīng)有13人回復(fù)
fortran的debug調(diào)試已經(jīng)有10人回復(fù)
求高手解答Matlab求解微分方程的初始值和邊界值問(wèn)題已經(jīng)有5人回復(fù)
請(qǐng)matlab高手看一下這個(gè)問(wèn)題怎樣編程才好，附程序代碼，謝謝！已經(jīng)有9人回復(fù)
MATLAB程序調(diào)試已經(jīng)有9人回復(fù)
求Matlab高手解決線(xiàn)性方程組的迭代求解問(wèn)題已經(jīng)有17人回復(fù)
Matlab中如何自己定義一個(gè)新函數(shù)，然后計(jì)算時(shí)可以調(diào)用？已經(jīng)有12人回復(fù)
matlab調(diào)試問(wèn)題已經(jīng)有3人回復(fù)
matlab程序問(wèn)題，總是出錯(cuò)，麻煩大家看看問(wèn)題在那里？已經(jīng)有6人回復(fù)
matlab重復(fù)調(diào)用子程序的問(wèn)題已經(jīng)有6人回復(fù)
求助大俠幫忙，有關(guān)matlab模擬數(shù)據(jù) 已經(jīng)有29人回復(fù)
求助matlab編程已經(jīng)有13人回復(fù)
VS2005下matlab調(diào)用已經(jīng)有8人回復(fù)
【求助】一個(gè)MATLAB程序，讀不懂，老師說(shuō)是用隨機(jī)游走算pi的值【已完結(jié)】已經(jīng)有6人回復(fù)

1樓 2014-12-24 09:18:39

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

guan0603

新蟲(chóng) (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 67.1
帖子: 8
在線(xiàn): 9.8小時(shí)
蟲(chóng)號(hào): 3599995
注冊(cè): 2014-12-17
專(zhuān)業(yè): 運(yùn)籌學(xué)

引用回帖:

2樓: Originally posted by wurongjun at 2014-12-24 11:13:16
簡(jiǎn)單修改就行,如下:
function =frcg(fun,gfun,x0)
% 功能: 用FR共軛梯度法求解無(wú)約束問(wèn)題: min f(x)
%輸入: x0是初始點(diǎn), fun, gfun分別是目標(biāo)函數(shù)和梯度
%輸出: x, val分別是近似最優(yōu)點(diǎn)和最優(yōu)值, k是迭代次數(shù).
...

>> x0=[-1.2 2]'

x0 =

-1.2000
2.0000

>> [x,val,k]=frcg('fun','gfun',x0)
??? Error using ==> frcg
Too many input arguments.
這是怎么回事呢，我剛接觸MATLAB，編程基礎(chǔ)很差，謝謝幫助

贊一下

回復(fù)此樓

3樓2014-12-24 11:58:06

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

查看全部 4 個(gè)回答

wurongjun

專(zhuān)家顧問(wèn) (職業(yè)作家)

專(zhuān)家經(jīng)驗(yàn): +831
數(shù)學(xué)EPI: 9
應(yīng)助: 791 (博后)
貴賓: 0.308
金幣: 24609
散金: 310
紅花: 75
帖子: 3004
在線(xiàn): 881.4小時(shí)
蟲(chóng)號(hào): 1368482
注冊(cè): 2011-08-14
性別: GG
專(zhuān)業(yè): 計(jì)算數(shù)學(xué)與科學(xué)工程計(jì)算
管轄: 數(shù)學(xué)

【答案】應(yīng)助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感謝參與，應(yīng)助指數(shù) +1
guan0603: 金幣+10, ★★★很有幫助 2014-12-24 11:58:46

簡(jiǎn)單修改就行,如下:
function [x,val,k]=frcg(fun,gfun,x0)
% 功能: 用FR共軛梯度法求解無(wú)約束問(wèn)題: min f(x)
%輸入: x0是初始點(diǎn), fun, gfun分別是目標(biāo)函數(shù)和梯度
%輸出: x, val分別是近似最優(yōu)點(diǎn)和最優(yōu)值, k是迭代次數(shù).
maxk=5000; %最大迭代次數(shù)
rho=0.6;sigma=0.4;
k=0;epsilon=1e-4;
n=length(x0);
while(k<maxk)
g=feval(gfun,x0); %計(jì)算梯du
itern=k-(n+1)*floor(k/(n+1));
itern=itern+1;
%計(jì)算搜索方向
if(itern==1)
d=-g;
else
beta=(g'*g)/(g0'*g0);
d=-g+beta*d0; gd=g'*d;
if(gd>=0.0)
d=-g;
end
end
if(norm(g)<epsilon),break;end %檢驗(yàn)終止條件
m=0;mk=0;
while(m<20) %Armijo搜索
if(feval(fun,x0+rho^m*d)<feval(fun,x0)+sigma*rho^m*g'*d)
mk=m;break;
end
m=m+1;
end
x0=x0+rho^mk*d;
val=feval(fun,x0);
g0=g;d0=d;
k=k+1;
end
x=x0;
val=feval(fun,x);
function g=gfun(x)
g=[400*x(1)*(x(1)^2-x(2))+2*(x(1)-1),-200*(x(1)^2-x(2))]';
function f=fun(x)
f=100*(x(1)^2-x(2))^2+(x(1)-1)^2;

下面是運(yùn)行結(jié)果:
>> x0=[-1.2 1]'
x0 =
-1.2000
1.0000
>> [x,val,k]=frcg('fun','gfun',x0)
x =
0.9999
0.9999
val =
2.9396e-009
k =
44

贊一下(1人)

回復(fù)此樓

» 本帖已獲得的紅花（最新10朵）

guan0603

善惡到頭終有報(bào),人間正道是滄桑.

2樓2014-12-24 11:13:16

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

guan0603

新蟲(chóng) (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 67.1
帖子: 8
在線(xiàn): 9.8小時(shí)
蟲(chóng)號(hào): 3599995
注冊(cè): 2014-12-17
專(zhuān)業(yè): 運(yùn)籌學(xué)

送紅花一朵

引用回帖:

出來(lái)了，謝謝了

回復(fù)此樓

4樓2014-12-24 14:26:39

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

最具人氣熱帖推薦 [查看全部]	作者	回/看	最后發(fā)表

[考研] 求調(diào)劑323材料與化工 +3	1124361 2026-03-24	3/150	2026-03-24 20:24 by peike
[考研] 求調(diào)劑一志愿武漢理工大學(xué)材料工程（085601） +5	WW.' 2026-03-23	7/350	2026-03-24 14:50 by sprinining
[考研] 293求調(diào)劑 +6	加一一九 2026-03-24	6/300	2026-03-24 14:29 by JourneyLucky
[考研] 求材料，環(huán)境專(zhuān)業(yè)調(diào)劑 +3	18567500178 2026-03-18	3/150	2026-03-23 23:50 by 熱情沙漠
[考研] 335求調(diào)劑 +4	yuyu宇 2026-03-23	5/250	2026-03-23 23:49 by Txy@872106
[考研] 上海電力大學(xué)材料防護(hù)與新材料重點(diǎn)實(shí)驗(yàn)室招收調(diào)劑研究生（材料、化學(xué)、電化學(xué)，環(huán)境） +3	我愛(ài)學(xué)電池 2026-03-23	3/150	2026-03-23 17:16 by AZMK
[考研] 0854電子信息求調(diào)劑 324 +3	Promise-jyl 2026-03-23	3/150	2026-03-23 13:43 by wangkm
[考研] 求調(diào)劑材料學(xué)碩080500，總分289分 5+3	@taotao 2026-03-19	21/1050	2026-03-23 10:17 by 冠c哥
[考研] 070300，一志愿北航320求調(diào)劑 +3	Jerry0216 2026-03-22	5/250	2026-03-23 09:16 by 。。堂堂
[考研] 287求調(diào)劑 +8	晨昏線(xiàn)與星海 2026-03-19	9/450	2026-03-22 17:01 by i_cooler
[考研] 306求調(diào)劑 +5	來(lái)好運(yùn)來(lái)來(lái)來(lái) 2026-03-22	5/250	2026-03-22 16:17 by BruceLiu320
[考研] 275求調(diào)劑 +6	shansx 2026-03-22	8/400	2026-03-22 15:27 by barlinike
[考研] 085600材料與化工306 +4	z1z2z3879 2026-03-21	4/200	2026-03-21 23:44 by ms629
[考研] 一志愿南大，0703化學(xué)，分?jǐn)?shù)336，求調(diào)劑 +3	收到VS 2026-03-21	3/150	2026-03-21 18:42 by 學(xué)員8dgXkO
[考研] 308求調(diào)劑 +3	阿姐阿姐家啊 2026-03-18	3/150	2026-03-20 23:24 by JourneyLucky
[考研] 一志愿武漢理工材料工程專(zhuān)碩調(diào)劑 +9	Doleres 2026-03-19	9/450	2026-03-20 22:36 by JourneyLucky
[考研] 材料與化工 322求調(diào)劑 +4	然11 2026-03-19	4/200	2026-03-20 22:12 by luoyongfeng
[考研] A區(qū)線(xiàn)材料學(xué)調(diào)劑 +5	周周無(wú)極 2026-03-20	5/250	2026-03-20 21:33 by laoshidan
[考研] 一志愿南京航空航天大學(xué)大學(xué) ，080500材料科學(xué)與工程學(xué)碩 +5	@taotao 2026-03-20	5/250	2026-03-20 20:16 by JourneyLucky
[考研] 320求調(diào)劑0856 +3	不想起名字112 2026-03-19	3/150	2026-03-19 22:53 by 學(xué)員8dgXkO