²é¿´: 720  |  »Ø¸´: 3

champion3568

ľ³æ (Ö°Òµ×÷¼Ò)

[ÇóÖú] Çë¸÷λ´óÉñ°ïÎÒ¿´¿´¹ØÓÚÄâÅ£¶Ù·¨Maltab´úÂëµ½µ×ÄÄÀïÓÐ´í£¿£¿

´úÂëÎÒÊÇÒÀÕÕÒ»±¾½Ì³Ì×öµÄ£¬¿ÉÔËÐÐʱ¾Í³öÏÖÈçÏ´íÎó£º
syms x1 x2
>> x0=[1,2]';
>> [x,val,k]=sr1(fun,gfun,x0)
??? Undefined function or variable 'fun'.
Õâ¸öÊÇ´úÂ룺
function [x,val,k] =sr1(fun,gfun,x0)
%¹¦ÄÜ£ºÓöԳÆÖÈ1Ëã·¨Çó½âÎÞÔ¼ÊøÎÊÌ⣺min f(x)
%ÊäÈ룺X0Êdzõʼµã£¬fun,gfun·Ö±ðÊÇÄ¿±êº¯Êý¼°ÆäÌݶȠ  
%Êä³ö£ºx£¬val·Ö±ðÊǽüËÆ×îÓŵãºÍ×îÓÅÖµ£¬kÊǵü´ú´ÎÊý.
maxk=500; %¸ø³ö×î´óµü´ú´ÎÊý
rho=0.55;sigma=0.4;epsilon=le-4;
k=0; n=length(x0);Hk=eye(n);
while(k<maxk)
    gk=feval(gfun,x0);%¼ÆËãÌݶÈ
    dk=-Hk*gk;%¼ÆËãËÑË÷·½Ïò
    if(norm(gk)<epsilon),break;end  %¼ìÑéÖÕÖ¹×¼Ôò
    m=0;mk=0;
    while(m<20)  %ÓÃArmijoËÑË÷Çó²½³¤
      if(feval(fun,x0+rho^m*dk)<feval(fun,x0)+sigma*rho^m*gk'*dk)   
         mk=m;break;
      end
    m=m+1;
    end
    x=x0+rho^mk*dk;
    sk=x-x0; yk=feval(gfun,x)-gk;
    Hk=Hk+(sk-Hk*yk)*(sk-Hk*yk)'/((sk-Hk*yk)'*yk); %ÖÈ1УÕý
    k=k+1;  x0=x;
end
val=feval(fun,x0);
»Ø¸´´ËÂ¥
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû

yanshuiting

гæ (ÕýʽдÊÖ)

2Â¥2015-11-20 22:28:12
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû

heaven0419

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

Ö÷³ÌÐòÀïÃæÃ»ÓиøfunµÄÖµ£¿

·¢×ÔСľ³æIOS¿Í»§¶Ë
3Â¥2015-11-20 22:31:20
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû

champion3568

ľ³æ (Ö°Òµ×÷¼Ò)

ÒýÓûØÌû:
2Â¥: Originally posted by yanshuiting at 2015-11-20 22:28:12
ÄãµÄÊäÈë²»¶Ô

ÄÇôÈçºÎ½â¾öÄØ
4Â¥2015-11-20 23:45:39
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû
Ïà¹Ø°æ¿éÌø×ª ÎÒÒª¶©ÔÄÂ¥Ö÷ champion3568 µÄÖ÷Ìâ¸üÐÂ
×î¾ßÈËÆøÈÈÌûÍÆ¼ö [²é¿´È«²¿] ×÷Õß »Ø/¿´ ×îºó·¢±í
[¿¼ÑÐ] 22408 344·Ö Çóµ÷¼Á Ò»Ö¾Ô¸ »ªµç¼ÆËã»ú¼¼Êõ +4 solanXXX 2026-03-20 4/200 2026-03-20 23:49 by alg094825
[¿¼ÑÐ] Äϲý´óѧ²ÄÁÏר˶311·ÖÇóµ÷¼Á +5 77chaselx 2026-03-20 5/250 2026-03-20 23:42 by lovewei0727
[¿¼ÑÐ] ±±¿Æ281ѧ˶²ÄÁÏÇóµ÷¼Á +5 tcxiaoxx 2026-03-20 5/250 2026-03-20 21:35 by laoshidan
[¿¼ÑÐ] ²ÄÁÏѧ˶297ÒѹýËÄÁù¼¶Çóµ÷¼ÁÍÆ¼ö +11 adaie 2026-03-19 11/550 2026-03-20 21:30 by laoshidan
[¿¼ÑÐ] 316Çóµ÷¼Á +5 ÁºÜçö© 2026-03-19 5/250 2026-03-20 21:26 by ÄãºÃÄãºÃ666
[¿¼ÑÐ] Ò»Ö¾Ô¸Î÷ÄϽ»Í¨ ר˶ ²ÄÁÏ355 ±¾¿ÆË«·Ç Çóµ÷¼Á +5 Î÷ÄϽ»Í¨×¨²Ä355 2026-03-19 5/250 2026-03-20 21:10 by JourneyLucky
[¿¼ÑÐ] Ò»Ö¾Ô¸¼ªÁÖ´óѧ²ÄÁÏѧ˶321Çóµ÷¼Á +11 Ymlll 2026-03-18 15/750 2026-03-20 19:40 by ¶¡¶¡*
[¿¼ÑÐ] 0703»¯Ñ§µ÷¼Á £¬Áù¼¶Òѹý£¬ÓпÆÑо­Àú +13 êØÎõÙâ 2026-03-15 13/650 2026-03-20 19:35 by Dream007008
[¿¼ÑÐ] 261ÇóBÇøµ÷¼Á£¬¿ÆÑо­Àú·á¸» +3 Å£Ä̺Üæ 2026-03-20 4/200 2026-03-20 19:34 by JourneyLucky
[¿¼ÑÐ] Çóµ÷¼Á +3 eation27 2026-03-20 3/150 2026-03-20 19:32 by JourneyLucky
[¿¼ÑÐ] 086500 325 Çóµ÷¼Á +3 Áì´øÐ¡ÐÜ 2026-03-19 3/150 2026-03-20 18:38 by ¾¡Ë´Ò¢1
[¿¼ÑÐ] 281Çóµ÷¼Á£¨0805£© +14 ÑÌÏ«Ò亣 2026-03-16 25/1250 2026-03-20 15:47 by yuncha
[¿¼ÑÐ] 265Çóµ÷¼Á +8 ÁºÁºÐ£Ð£ 2026-03-17 8/400 2026-03-20 14:40 by 27µÀ¿ÆÌØ
[¿¼ÑÐ] 085600²ÄÁÏÓ뻯¹¤Çóµ÷¼Á +6 Ð÷ÐÒÓë×Ó 2026-03-17 6/300 2026-03-19 13:27 by houyaoxu
[¿¼ÑÐ] 344Çóµ÷¼Á +6 knight344 2026-03-16 7/350 2026-03-18 20:13 by walc
[¿¼ÑÐ] 085601ר˶£¬×Ü·Ö342Çóµ÷¼Á£¬µØÇø²»ÏÞ +5 share_joy 2026-03-16 5/250 2026-03-18 14:48 by haxia
[¿¼ÑÐ] 275Çóµ÷¼Á +4 Ì«Ñô»¨ÌìÌ쿪ÐÄ 2026-03-16 4/200 2026-03-17 10:53 by ¹¦·ò·è¿ñ
[¿¼ÑÐ] 326Çóµ÷¼Á +4 ŵ±´¶û»¯Ñ§½±êéê 2026-03-15 7/350 2026-03-16 17:11 by ŵ±´¶û»¯Ñ§½±êéê
[¿¼ÑÐ] 0856Çóµ÷¼Á +3 ÁõÃÎ΢ 2026-03-15 3/150 2026-03-16 10:00 by houyaoxu
[¿¼ÑÐ] 080500£¬²ÄÁÏѧ˶302·ÖÇóµ÷¼ÁѧУ +4 ³õʶ¿ÉÀÖ 2026-03-14 5/250 2026-03-14 21:08 by peike
ÐÅÏ¢Ìáʾ
ÇëÌî´¦ÀíÒâ¼û