 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 4 | 1/1 | 返回列表 |
| 查看: 883 | 回復: 3 | ||||
[求助]
哪位MATLAB高手能幫忙調試一下程序，不是遇到這問題就是那問題 已有1人參與
|
|
function [x,val,k]=frcg(fun,gfun,x0) % 功能: 用FR共軛梯度法求解無約束問題: min f(x) %輸入: x0是初始點, fun, gfun分別是目標函數(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)點和最優(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); %計算梯du itern=k-(n+1)*floor(k/(n+1)); itern=itern+1; %計算搜索方向 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 %檢驗終止條件 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'. 怎么回事呢 |
回復此樓» 猜你喜歡
289求調劑
已經有4人回復
-
292求調劑
已經有4人回復
-
080500求調劑
已經有3人回復
-
291求調劑
已經有3人回復
-
材料292調劑
已經有8人回復
-
一志愿211 初試270分 求調劑
已經有6人回復
-
求調劑
已經有3人回復
-
274求調劑
已經有4人回復
-
300求調劑，材料科學英一數(shù)二
已經有5人回復
-
招08考數(shù)學
已經有15人回復
» 本主題相關價值貼推薦，對您同樣有幫助:
- 高手指點：如何使用matlab使數(shù)據行與行間空一行
已經有6人回復
- matlab編程 ODE45
已經有4人回復
- linux 下 matlab 的調用問題
已經有3人回復
- matlab畫圖問題: 怎樣用小圖標或自定義圖形代替+ x . o等做圖
已經有10人回復
- 一個關于matlab的simulink的變頻空調PI控制的問題
已經有9人回復
- 請大俠指教matlab的這個錯誤提示到底是個神馬意思？
已經有12人回復
- 有個matlab程序，但是運行不下來，求高手
已經有23人回復
- matlab程序中調用m文件的問題
已經有13人回復
- fortran的debug調試
已經有10人回復
- 求高手解答Matlab求解微分方程的初始值和邊界值問題
已經有5人回復
- 請matlab高手看一下這個問題怎樣編程才好，附程序代碼，謝謝！
已經有9人回復
- MATLAB程序調試
已經有9人回復
- 求Matlab高手解決線性方程組的迭代求解問題
已經有17人回復
- Matlab中如何自己定義一個新函數(shù)，然后計算時可以調用？
已經有12人回復
- matlab調試問題
已經有3人回復
- matlab程序問題，總是出錯，麻煩大家看看問題在那里？
已經有6人回復
- matlab重復調用子程序的問題
已經有6人回復
- 求助大俠幫忙，有關matlab模擬數(shù)據
已經有29人回復
- 求助matlab編程
已經有13人回復
- VS2005下matlab調用
已經有8人回復
- 【求助】一個MATLAB程序，讀不懂，老師說是用隨機游走算pi的值【已完結】
已經有6人回復
wurongjun
專家顧問 (職業(yè)作家)
-
專家經驗: +831 - 數(shù)學EPI: 9
- 應助: 791 (博后)
- 貴賓: 0.308
- 金幣: 24609
- 散金: 310
- 紅花: 75
- 帖子: 3004
- 在線: 881.4小時
- 蟲號: 1368482
- 注冊: 2011-08-14
- 性別: GG
- 專業(yè): 計算數(shù)學與科學工程計算
- 管轄: 數(shù)學
【答案】應助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感謝參與，應助指數(shù) +1
guan0603: 金幣+10, ★★★很有幫助 2014-12-24 11:58:46
感謝參與，應助指數(shù) +1
guan0603: 金幣+10, ★★★很有幫助 2014-12-24 11:58:46
|
簡單修改就行,如下: function [x,val,k]=frcg(fun,gfun,x0) % 功能: 用FR共軛梯度法求解無約束問題: min f(x) %輸入: x0是初始點, fun, gfun分別是目標函數(shù)和梯度 %輸出: x, val分別是近似最優(yōu)點和最優(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); %計算梯du itern=k-(n+1)*floor(k/(n+1)); itern=itern+1; %計算搜索方向 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 %檢驗終止條件 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; 下面是運行結果: >> 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 |
» 本帖已獲得的紅花（最新10朵）
guan0603
2樓2014-12-24 11:13:16
3樓2014-12-24 11:58:06
4樓2014-12-24 14:26:39
相關版塊跳轉
我要訂閱樓主 guan0603 的主題更新
| 4 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 086003食品工程求調劑 +4 | 淼淼111 2026-03-24 | 4/200 |
|
|---|---|---|---|---|
|
[考研] 344求調劑 +3 | desto 2026-03-24 | 3/150 |
|
|
|
[考研] 求調劑 +7 | 十三加油 2026-03-21 | 7/350 |
|
|
|
[考研] 一志愿國科過程所081700，274求調劑 +3 | 三水研0水立方 2026-03-23 | 3/150 |
|
|
|
[考研] 269求調劑 +4 | 我想讀研11 2026-03-23 | 4/200 |
|
|
|
[考研] 333求調劑 +3 | ALULU4408 2026-03-23 | 3/150 |
|
|
|
[考研] 070300，一志愿北航320求調劑 +3 | Jerry0216 2026-03-22 | 5/250 |
|
|
|
[考研] 275求調劑 +6 | shansx 2026-03-22 | 8/400 |
|
|
|
[考研] 311求調劑 +3 | 26研0 2026-03-20 | 3/150 |
|
|
|
[考研] 求調劑 +5 | Zhangbod 2026-03-21 | 7/350 |
|
|
|
[考研] 0805 316求調劑 +3 | 大雪深藏 2026-03-18 | 3/150 |
|
|
|
[考研] 0703化學調劑 +4 | 妮妮ninicgb 2026-03-21 | 4/200 |
|
|
|
[考研] 265求調劑 +12 | 梁梁校校 2026-03-19 | 14/700 |
|
|
|
[考研] 330求調劑0854 +3 | assdll 2026-03-21 | 3/150 |
|
|
|
[考研] 085601調劑 358分 +3 | zzzzggh 2026-03-20 | 4/200 |
|
|
|
[考研] 華東師范大學-071000生物學-293分-求調劑 +3 | 研究生何瑤明 2026-03-18 | 3/150 |
|
|
|
[考研] 一志愿華中科技大學，080502，354分求調劑 +5 | 守候夕陽CF 2026-03-18 | 5/250 |
|
|
|
[考研] 274求調劑 +10 | S.H1 2026-03-18 | 10/500 |
|
|
|
[考研] 求調劑一志愿南京航空航天大學289分 +3 | @taotao 2026-03-19 | 3/150 |
|
|
|
[考研] 一志愿西安交通大學 學碩 354求調劑211或者雙一流 +3 | 我想要讀研究生 2026-03-20 | 3/150 |
|