²é¿´: 907  |  »Ø¸´: 5

jydtjm

гæ (³õÈëÎÄ̳)

[ÇóÖú] ¹ØÓÚNewton-pcg·¨½â·½³ÌµÄ±à³Ì˼· ÒÑÓÐ1È˲ÎÓë

ÏȰÝл£¬ÓÃNewton-pcg·¨½â·½³Ì£¬ÓÃMATLAB±à³Ì£¬±à³ÌÖ÷¿ò¼ÜÓ¦¸ÃÈçºÎд£¿£¿£¿ÓÐÄÄλ´óÉñ¿ÉÒÔ°ïæ½â´ð£¬Ð»Ð»£¬·Ç³£¸Ðл£¬ÒòΪĿǰֻÓÐÕâô¶à½ð±Ò£¬Çë¶®ÕâÒ»·½ÃæµÄ´óÉñ´Í½Ì£¬ÔٴθÐл£¡£¡£¡
»Ø¸´´ËÂ¥
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû

jydtjm

гæ (³õÈëÎÄ̳)

Çó½âÏßÐÔ·½³Ì×éµÄnewton-pcg²¢ÐÐËã·¨£¬¼´¶þ´Îµü´ú·¨£¬ÔÚnewton·¨µÄÿһ²½µü´úÖУ¬²ÉÓÃpcg·¨Çó½ânewton·½³Ì×飬ÓÃMATLAB±à³Ì£¬±à³ÌÖ÷¿ò¼ÜÓ¦¸ÃÈçºÎд£¿£¿£¿ÓÐÄÄλ´óÉñ¿ÉÒÔ°ïæ½â´ð£¬Ð»Ð»£¬·Ç³£¸Ðл
2Â¥2015-03-30 15:54:27
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû

wurongjun

ר¼Ò¹ËÎÊ (Ö°Òµ×÷¼Ò)

¡¾´ð°¸¡¿Ó¦Öú»ØÌû

¸Ðл²ÎÓ룬ӦÖúÖ¸Êý +1
ÒÔǰÓÐÈËÎʹýÕâ¸ö!
ÔÙ¸øÄãÌùÒ»ÏÂ:
function [x,val,k]=frcg(fun,gfun,x0)
% ¹¦ÄÜ: ÓÃFR¹²éîÌݶȷ¨Çó½âÎÞÔ¼ÊøÎÊÌâ: min f(x)
%ÊäÈë: x0Êdzõʼµã, fun, gfun·Ö±ðÊÇÄ¿±êº¯ÊýºÍÌݶÈ
%Êä³ö: x, val·Ö±ðÊǽüËÆ×îÓŵãºÍ×îÓÅÖµ, kÊǵü´ú´ÎÊý.
maxk=5000; %×î´óµü´ú´ÎÊý
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;
ÉÆ¶ñµ½Í·ÖÕÓб¨,È˼äÕýµÀÊDz×É£.
3Â¥2015-03-30 17:16:29
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû

jydtjm

гæ (³õÈëÎÄ̳)

ÒýÓûØÌû:
3Â¥: Originally posted by wurongjun at 2015-03-30 17:16:29
ÒÔǰÓÐÈËÎʹýÕâ¸ö!
ÔÙ¸øÄãÌùÒ»ÏÂ:
function =frcg(fun,gfun,x0)
% ¹¦ÄÜ: ÓÃFR¹²éîÌݶȷ¨Çó½âÎÞÔ¼ÊøÎÊÌâ: min f(x)
%ÊäÈë: x0Êdzõʼµã, fun, gfun·Ö±ðÊÇÄ¿±êº¯ÊýºÍÌݶÈ
%Êä³ö: x, val·Ö±ðÊǽüËÆ×îÓŵãºÍ×îÓÅÖµ, ...

ÄúºÃ£¬Ê×ÏȸÐлÄúµÄ»Ø¸´£¬µ«ÊÇË¡ÎÒÓÞÃÁ£¬ÊµÔÚ¿´²»³öÕâ»Ø´ðÓëÎÒÕâÎÊÌâÓкÎÁªÏµ£¬ÇëÖ¸½Ì£¡
4Â¥2015-03-30 21:21:45
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû

wurongjun

ר¼Ò¹ËÎÊ (Ö°Òµ×÷¼Ò)

ÒýÓûØÌû:
4Â¥: Originally posted by jydtjm at 2015-03-30 21:21:45
ÄúºÃ£¬Ê×ÏȸÐлÄúµÄ»Ø¸´£¬µ«ÊÇË¡ÎÒÓÞÃÁ£¬ÊµÔÚ¿´²»³öÕâ»Ø´ðÓëÎÒÕâÎÊÌâÓкÎÁªÏµ£¬ÇëÖ¸½Ì£¡...

Sorry!
Õ³´íÀ²!
ÉÆ¶ñµ½Í·ÖÕÓб¨,È˼äÕýµÀÊDz×É£.
5Â¥2015-03-31 08:26:17
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû

jydtjm

гæ (³õÈëÎÄ̳)

ÒýÓûØÌû:
5Â¥: Originally posted by wurongjun at 2015-03-31 08:26:17
Sorry!
Õ³´íÀ²!...

ûÊ£¬Ð»Ð»£¡ÇëÎÊÄúÖªµÀÕâÖֽⷨµÄÏà¹Ø³ÌÐòÂð£¿
6Â¥2015-03-31 14:57:49
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû
Ïà¹Ø°æ¿éÌø×ª ÎÒÒª¶©ÔÄÂ¥Ö÷ jydtjm µÄÖ÷Ìâ¸üÐÂ
×î¾ßÈËÆøÈÈÌûÍÆ¼ö [²é¿´È«²¿] ×÷Õß »Ø/¿´ ×îºó·¢±í
[¿¼ÑÐ] 0703»¯Ñ§µ÷¼Á£¬Çóµ¼Ê¦ÊÕ +6 ÌìÌìºÃÔËÀ´Éϰ¶° 2026-03-24 6/300 2026-03-24 17:19 by dick_runner
[¿¼ÑÐ] Çóµ÷¼Á +6 ÑÐÑУ¬½Óµç»° 2026-03-24 7/350 2026-03-24 17:01 by barlinike
[¿¼ÑÐ] Ò»Ö¾Ô¸211 ³õÊÔ270·Ö Çóµ÷¼Á +5 ¹ÈÓêÉϰ¶ 2026-03-23 6/300 2026-03-24 16:32 by laoshidan
[¿¼ÑÐ] 300Çóµ÷¼Á£¬²ÄÁÏ¿ÆÑ§Ó¢Ò»Êý¶þ +5 leaflight 2026-03-24 5/250 2026-03-24 16:25 by laoshidan
[¿¼ÑÐ] ²ÄÁÏר˶331Çóµ÷¼Á +4 Ïʵ±Å£ 2026-03-24 4/200 2026-03-24 15:58 by JourneyLucky
[¿¼ÑÐ] 307Çóµ÷¼Á +3 ÓàÒâÇä 2026-03-21 6/300 2026-03-24 15:03 by ÓàÒâÇä
[¿¼ÑÐ] ²ÄÁÏר˶ӢһÊý¶þ306 +8 z1z2z3879 2026-03-18 8/400 2026-03-23 20:49 by baobaoye
[¿¼ÑÐ] 336Çóµ÷¼Á +4 ÊÕµ½VS 2026-03-20 4/200 2026-03-23 19:02 by macy2011
[¿¼ÑÐ] 336»¯¹¤µ÷¼Á +4 Íõ´ó̹1 2026-03-23 5/250 2026-03-23 18:32 by allen-yin
[¿¼ÑÐ] 323Çóµ÷¼Á +6 ÍÝСͰ 2026-03-18 6/300 2026-03-23 00:29 by king123£¡
[¿¼ÑÐ] Ò»Ö¾Ô¸Î÷°²½»Í¨´óѧ²ÄÁϹ¤³Ìרҵ 282·ÖÇóµ÷¼Á +11 ·ãÇÅZL 2026-03-18 13/650 2026-03-22 20:26 by edmund7
[¿¼ÑÐ] ѰÕÒµ÷¼Á +4 ¾óǿâ? 2026-03-21 4/200 2026-03-22 16:14 by ľÍÐζ¶
[¿¼ÑÐ] 260Çóµ÷¼Á +3 ÖìÜÆÁÕ 2026-03-20 4/200 2026-03-22 15:12 by ÖìÜÆÁÕ
[¿¼ÑÐ] 303Çóµ÷¼Á +5 °²ÒäÁé 2026-03-22 6/300 2026-03-22 12:46 by ËØÑÕÇã³Ç1988
[¿¼ÑÐ] ²ÄÁÏ 271Çóµ÷¼Á +5 Õ¹ÐÅÔÃ_ 2026-03-21 5/250 2026-03-21 17:29 by ѧԱ8dgXkO
[¿¼ÑÐ] Çóµ÷¼Á +6 Mqqqqqq 2026-03-19 6/300 2026-03-21 08:04 by JourneyLucky
[¿¼ÑÐ] 324·Ö 085600²ÄÁÏ»¯¹¤Çóµ÷¼Á +4 llllkkkhh 2026-03-18 4/200 2026-03-21 01:24 by JourneyLucky
[¿¼ÑÐ] ²ÄÁÏѧ˶297ÒѹýËÄÁù¼¶Çóµ÷¼ÁÍÆ¼ö +11 adaie 2026-03-19 11/550 2026-03-20 21:30 by laoshidan
[¿¼ÑÐ] Çóµ÷¼Á +3 @taotao 2026-03-20 3/150 2026-03-20 19:35 by JourneyLucky
[¿¼ÑÐ] Ò»Ö¾Ô¸ÄÏÀí¹¤085701»·¾³302Çóµ÷¼ÁԺУ +3 ¿ûè÷ÎÀ¶Ó 2026-03-20 3/150 2026-03-20 19:28 by zhukairuo
ÐÅÏ¢Ìáʾ
ÇëÌî´¦ÀíÒâ¼û