±±¾©Ê¯ÓÍ»¯¹¤Ñ§Ôº2026ÄêÑо¿ÉúÕÐÉú½ÓÊÕµ÷¼Á¹«¸æ
²é¿´: 2349  |  »Ø¸´: 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

гæ (СÓÐÃûÆø)

±ÈÈç 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µÄ»ØÌû
²é¿´È«²¿ 11 ¸ö»Ø´ð

buhely

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

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

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

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µÄ»ØÌû
×î¾ßÈËÆøÈÈÌûÍÆ¼ö [²é¿´È«²¿] ×÷Õß »Ø/¿´ ×îºó·¢±í
[¿¼ÑÐ] »¯Ñ§308·Öµ÷¼Á +14 ÄãºÃÃ÷ÌìÄãºÃ 2026-03-30 15/750 2026-04-02 13:39 by ÎûÎûЦ22
[»ù½ðÉêÇë] ÇëÎʹ²Í¬Í¨Ñ¶ºÍ¹²Í¬Ò»×÷µÄÈϿɶÈÎÊÌâ 10+3 psa1234 2026-04-01 4/200 2026-04-02 13:24 by psa1234
[¿¼ÑÐ] 311Çóµ÷¼Á +14 À¶ÔÂÁÁÁÁ 2026-03-30 14/700 2026-04-02 12:18 by 1753564080
[¿¼ÑÐ] 366Çóµ÷¼ÁÒ»Ö¾Ô¸¶«±±´óѧ +6 ÔËÆøÀ´µÃÈôÓÐËÆÎ 2026-04-02 6/300 2026-04-02 10:31 by ²»³Ôô~µÄ؈
[¿¼ÑÐ] 274Çóµ÷¼Á +8 Ѧ¶¨Ú̵Ļ¢¡£ 2026-04-01 8/400 2026-04-02 09:21 by Çç¿Õ210210
[¿¼ÑÐ] Ò»Ö¾Ô¸°²»Õ´óѧ¼ÆËã»ú¿ÆÑ§Óë¼¼Êõѧ˶£¬331·ÖÇóµ÷¼Á +5 ½¯²ýÅôqtj 2026-04-01 5/250 2026-04-02 08:10 by fxue1114
[¿¼ÑÐ] 285Çóµ÷¼Á +7 AZMK 2026-03-30 13/650 2026-04-01 17:00 by Æß¶È²»ÐÅÈÎ
[¿¼ÑÐ] 0703Çóµ÷¼Á +4 zizimo 2026-03-31 4/200 2026-04-01 16:04 by yanflower7133
[¿¼ÑÐ] 285Çóµ÷¼Á +5 FZAC123 2026-03-30 5/250 2026-04-01 15:50 by º«Ó꺭
[¿¼ÑÐ] 358Çóµ÷¼Á +3 ÍõÏòÑô»¨ 2026-03-31 3/150 2026-04-01 09:56 by zzchen2000
[¿¼ÑÐ] 0856µ÷¼Á +7 ÇúÌýóÞ 2026-03-30 7/350 2026-04-01 08:51 by JourneyLucky
[¿¼ÑÐ] Ò»Ö¾Ô¸aÇø211£¬085601-307·ÖÇóµ÷¼Á +10 µ³¼ÎºÀ 2026-03-31 23/1150 2026-03-31 18:44 by JourneyLucky
[¿¼ÑÐ] ºÏ·ÊÇøÓòÐÔÖØµãÒ»±¾ÕÐÊÕµ÷¼Á +4 6266jl 2026-03-30 8/400 2026-03-31 18:43 by 6266jl
[¿¼ÑÐ] ±¾2Ò»Ö¾Ô¸C9-333·Ö£¬²ÄÁÏ¿ÆÑ§Ó빤³Ì£¬Çóµ÷¼Á +9 ÉýÉý²»½µ 2026-03-31 9/450 2026-03-31 18:01 by Î޼ʵIJÝÔ­
[¿¼ÑÐ] 329Çóµ÷¼Á£¬Ò»Ö¾Ô¸Î÷±±¹¤Òµ´óѧ£¬²ÄÁϹ¤³Ì£¨085601£© +6 СС»úÁ鳿 2026-03-29 12/600 2026-03-31 16:58 by ¼Çʱ¾2026
[¿¼ÑÐ] 286Çóµ÷¼Á +6 Faune 2026-03-30 6/300 2026-03-31 14:37 by jp9609
[¿¼ÑÐ] 22408 359·Öµ÷¼Á +4 Qshers 2026-03-27 8/400 2026-03-31 08:53 by Qshers
[¿¼ÑÐ] 08¹¤¿ÆÇóµ÷¼Á286 +5 tgs_001 2026-03-28 5/250 2026-03-31 08:18 by Ò»Ö»ºÃ¹û×Ó?
[¿¼ÑÐ] 085701»·¾³¹¤³ÌÇóµ÷¼Á +11 ¶à¾ÃÉϿΠ2026-03-27 12/600 2026-03-30 21:21 by Ñо¿É®µ¼µ¼
[¿¼ÑÐ] 340Çóµ÷¼Á +5 jhx777 2026-03-27 5/250 2026-03-28 04:18 by fmesaito
ÐÅÏ¢Ìáʾ
ÇëÌî´¦ÀíÒâ¼û