±±¾©Ê¯ÓÍ»¯¹¤Ñ§Ôº2026ÄêÑо¿ÉúÕÐÉú½ÓÊÕµ÷¼Á¹«¸æ
²é¿´: 2347  |  »Ø¸´: 10
µ±Ç°Ö»ÏÔʾÂú×ãÖ¸¶¨Ìõ¼þµÄ»ØÌû£¬µã»÷ÕâÀï²é¿´±¾»°ÌâµÄËùÓлØÌû

jjliu2012

гæ (СÓÐÃûÆø)

[ÇóÖú] ÐÂÊÖmatlabÇóÖú£¬Ð»Ð»¸÷λ

×îÓÅ¿ØÖÆÖеÄ×îËÙϽµ·¨£¬ÈôÊÇÇó²»³öº¯Êý¼¸ºõº¯ÊýÌݶȵĽâÎö±í´ïʽ£¬×îËÙϽµ·¨µÄ³ÌÐò¸ÃÔõôÑù´¦ÀíÄØ£¿
ÒÔÏÂÊÇÎÒÃÇÆ½Ê±ÓõÄ×îËÙϽµ·¨³ÌÐò

function [x,val,k]=grad(fun,gfun,x0)
%¹¦ÄÜ: ÓÃ×îËÙϽµ·¨Çó½âÎÞÔ¼ÊøÎÊÌâ: min f(x)
%ÊäÈë: x0Êdzõʼµã, fun, gfun·Ö±ðÊÇÄ¿±êº¯ÊýºÍÌݶÈ
%Êä³ö: x, val·Ö±ðÊǽüËÆ×îÓŵãºÍ×îÓÅÖµ, kÊǵü´ú´ÎÊý.
maxk=5000; %×î´óµü´ú´ÎÊý
rho=0.5;sigma=0.4;
k=0; epsilon=1e-5;
while(k¡maxk)
g=feval(gfun,x0); %¼ÆËãÌݶÈ
d=-g; %¼ÆËãËÑË÷·½Ïò
if(norm(d)¡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;
k=k+1;
end
x=x0;
val=feval(fun,x0);

%ÕâÀïµÄfunÊǺ¯Êý£¬fÊǺ¯Êý£¬gÊÇfº¯ÊýµÄÌݶȣ¬ÕâÀïµÄfºÍg¶¼ÓнâÎöµÄ±í´ïʽ
function f=fun(x)
f=100*(x(1)^2-x(2))^2+(x(1)-1)^2;
function g=gfun(x)
g=[400*x(1)*(x(1)^2-x(2))+2*(x(1)-1), -200*(x(1)^2-x(2))]¡¯;
»Ø¸´´ËÂ¥

» ²ÂÄãϲ»¶

» ±¾Ö÷ÌâÏà¹Ø¼ÛÖµÌùÍÆ¼ö£¬¶ÔÄúͬÑùÓаïÖú:

ÒÑÔÄ   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû

jjliu2012

гæ (СÓÐÃûÆø)

ÒýÓûØÌû:
2Â¥: Originally posted by buhely at 2013-09-26 10:00:19
Èç¹ûÄãµÄfunÊÇÓÃÃ÷È·µÄÊýѧÓïÑÔÃèÊöµÄ£¬µ«²»ÄÜÖ±½ÓËã³ögfun£¬¿ÉÒÔ¿¼ÂÇÓÃmatlabµÄ·ûºÅ±äÁ¿±à³Ì£¬È»ºóËã³ögfunµÄ±í´ïʽ£»
ÔÙ²»ÐÐÃ²ËÆ¾ÍÖ»ÄÜÔÙÌí¼Ó¸÷ÖÖÆô·¢Ê½ËÑËØËã·¨ÁË

ÕâÀïµÄxÂú×ãµÄÊÇÒ»¸ö·ÇÏßÐÔµÄ΢·Ö·½³Ì£¬x¾ÍûÓнâÎö½â£¬funÊǹØÓÚxµÄ£¬Ò²Ã»ÓнâÎö±í´ïʽ°¡
3Â¥2013-09-26 10:23:53
ÒÑÔÄ   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû
²é¿´È«²¿ 11 ¸ö»Ø´ð

buhely

Òø³æ (СÓÐÃûÆø)

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

¸Ðл²ÎÓ룬ӦÖúÖ¸Êý +1
Èç¹ûÄãµÄfunÊÇÓÃÃ÷È·µÄÊýѧÓïÑÔÃèÊöµÄ£¬µ«²»ÄÜÖ±½ÓËã³ögfun£¬¿ÉÒÔ¿¼ÂÇÓÃmatlabµÄ·ûºÅ±äÁ¿±à³Ì£¬È»ºóËã³ögfunµÄ±í´ïʽ£»
ÔÙ²»ÐÐÃ²ËÆ¾ÍÖ»ÄÜÔÙÌí¼Ó¸÷ÖÖÆô·¢Ê½ËÑËØËã·¨ÁË
2Â¥2013-09-26 10:00:19
ÒÑÔÄ   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû

virtualzx

ľ³æ (ÖøÃûдÊÖ)

ÒýÓûØÌû:
3Â¥: Originally posted by jjliu2012 at 2013-09-25 21:23:53
ÕâÀïµÄxÂú×ãµÄÊÇÒ»¸ö·ÇÏßÐÔµÄ΢·Ö·½³Ì£¬x¾ÍûÓнâÎö½â£¬funÊǹØÓÚxµÄ£¬Ò²Ã»ÓнâÎö±í´ïʽ°¡...

ÄãÄܲ»ÄܰÑÄãµÃÎÊÌâ˵µÃÔÙÃ÷°×µã£¬ÌùµÄ³ÌÐòÃ²ËÆºÍÄãµÄÎÊÌâÒ²²»Ïà¹Ø
4Â¥2013-09-26 11:24:25
ÒÑÔÄ   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû

jjliu2012

гæ (СÓÐÃûÆø)

±ÈÈç dx/dt=f(x(t)),ÕâÀïµÄfÊÇÒ»¸ö·ÇÏßÐԵĺ¯Êý¡£¶øfun=int£¨x^2,a,b),ÕâÖÖÇé¿öx£¨t£©Ã»ÓнâÎö±í´ïʽ
5Â¥2013-09-26 15:19:50
ÒÑÔÄ   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû
×î¾ßÈËÆøÈÈÌûÍÆ¼ö [²é¿´È«²¿] ×÷Õß »Ø/¿´ ×îºó·¢±í
[¿¼ÑÐ] 265Çóµ÷¼Á +3 ÁºÁºÐ£Ð£ 2026-04-01 3/150 2026-04-01 22:36 by barlinike
[¿¼ÑÐ] 0856£¬269·ÖÇóµ÷¼Á +8 ÓÐѧÉϾÍÐÐÇóÇóÁ 2026-03-30 11/550 2026-04-01 22:33 by 2026²ÄÁϵ÷¼Á
[¿¼ÑÐ] 286Çóµ÷¼Á +16 PolarBear11 2026-03-26 16/800 2026-04-01 21:31 by Æß¶È²»ÐÅÈÎ
[¿¼ÑÐ] 805600ר˶²ÄÁÏÓ뻯¹¤348·ÖÇóµ÷¼Á +4 ÉÏѧÀ²£¡ 2026-04-01 4/200 2026-04-01 21:10 by lijunpoly
[¿¼ÑÐ] µ÷¼Á +3 ºÃºÃ¶ÁÊé¡£ 2026-04-01 5/250 2026-04-01 15:39 by ºÃºÃ¶ÁÊé¡£
[¿¼ÑÐ] Ò»Ö¾Ô¸ ÄϾ©º½¿Õº½Ìì´óѧ £¬080500²ÄÁÏ¿ÆÑ§Ó빤³Ìѧ˶ +7 @taotao 2026-03-30 7/350 2026-04-01 14:30 by chenqifeng666
[˶²©¼ÒÔ°] ²©Ò»±»ËͳöÁªÅà¸Ð¾õ²»ÊÊÓ¦Ôõô°ì +3 È«´åµÄ¹· 2026-03-31 3/150 2026-04-01 10:44 by 328838485
[¿¼ÑÐ] 318Ò»Ö¾Ô¸¼ªÁÖ´óѧÉúÎïÓëÒ½Ò© Çóµ÷¼Á +6 óÆÐÐÖÂÔ¶. 2026-03-28 6/300 2026-04-01 09:28 by oooqiao
[¿¼ÑÐ] ±¾2Ò»Ö¾Ô¸C9-333·Ö£¬²ÄÁÏ¿ÆÑ§Ó빤³Ì£¬Çóµ÷¼Á +9 ÉýÉý²»½µ 2026-03-31 9/450 2026-03-31 18:01 by Î޼ʵIJÝÔ­
[¿¼ÑÐ] 336²ÄÁÏÇóµ÷¼Á +10 ³ÂäÞÓ¨ 2026-03-26 12/600 2026-03-31 17:59 by jp9609
[¿¼ÑÐ] Çó»¯Ñ§µ÷¼Á +12 wulanna 2026-03-28 12/600 2026-03-31 16:38 by 690616278
[¿¼ÑÐ] 353Çóµ÷¼Á +3 ½­ÉÏ·ã_26 2026-03-28 3/150 2026-03-31 15:53 by jp9609
[¿¼ÑÐ] 354Çóµ÷¼Á +3 lxb598 2026-03-31 4/200 2026-03-31 13:42 by sophie2180
[¿¼ÑÐ] 22408 359·Öµ÷¼Á +4 Qshers 2026-03-27 8/400 2026-03-31 08:53 by Qshers
[¿¼ÑÐ] ¸£½¨Àí¹¤´óѧ²ÄÁÏѧԺÏȽøºÏ½ðÍŶÓÕÐÊÕ¿¼Ñе÷¼ÁѧÉú +3 ´ó»ª½ðÉ̶¼ 2026-03-30 4/200 2026-03-31 01:04 by ·½Ó¢¿¡602
[¿¼ÑÐ] Çóµ÷¼Á +7 Çà´º²ÃΪÈý½Ø 2026-03-29 7/350 2026-03-30 13:14 by laoshidan
[¿¼ÑÐ] 348Çóµ÷¼Á +6 СÀÁ³æ²»ÀÁÁË 2026-03-28 6/300 2026-03-30 10:29 by Evan_Liu
[˶²©¼ÒÔ°] ÕÐÊÕÉúÎïѧ/ϸ°ûÉúÎïѧµ÷¼Á +4 IceGuo 2026-03-26 5/250 2026-03-29 01:25 by griffith2014
[¿¼ÑÐ] 081200-314 +3 LILIQQ 2026-03-27 4/200 2026-03-28 09:41 by ±£»¤µØÇòÄãÎÒ×öÆ
[¿¼ÑÐ] 081200-11408-276ѧ˶Çóµ÷¼Á +4 ´Þwj 2026-03-26 4/200 2026-03-27 08:04 by chemisry
ÐÅÏ¢Ìáʾ
ÇëÌî´¦ÀíÒâ¼û