±±¾©Ê¯ÓÍ»¯¹¤Ñ§Ôº2026ÄêÑо¿ÉúÕÐÉú½ÓÊÕµ÷¼Á¹«¸æ
²é¿´: 1511  |  »Ø¸´: 3

xuetb

Ìú³æ (³õÈëÎÄ̳)

[ÇóÖú] matlab×îС¶þ³ËÎÊÌâÇóÖú Ôõô¶¨Òå¶à±äÁ¿

×îС¶þ³ËÎÊÌ⣺
·ÇÏßÐÔ×îС¶þ³ËÎÊÌ⣺minf^2 ÆäÖУ¬ f=[x1+10*x2;sqrt(5)*(x3-x4);(x2-2*x3)^2;sqrt(10)*(x1-x4)^2];x0=[3,-1,0,1]
³ÌÐòÈçÏ£¬¿ÉÒÔÇóµ¥±äÁ¿£¬Ôõô¶¨Òåx=[x1,x2,x3,x4],Çó½âÄØ£¿  Çë¸ßÊÖÖ¸½Ì                        
function fv = Funval( f,varec,varval )
%UNTITLED Summary of this function goes here
%   Detailed explanation goes here
tic;
var = findsym(f);
varc = findsym(varec);
s1=length(var);
s2=length(varc);
m=floor((s1-1)/3+1);
varv=zeros(1,m);
if s1~=s2
    for i=0(s1-1)/3)
        k=findstr(varc,var(3*i+1));
        index=(k-1)/3;
        varv(i+1)=varval(index+1);
    end
    fv=subs(f,var,varv);
else
    fv=subs(f,varec,varval);
end
end
x=[x1,x2,x3,x4];
function [x,minf]= minGN(f,x0,var,eps)%Ä¿±êº¯Êý;f;³õʼµã:x0£¬×Ô±äÁ¿ÏòÁ¿;var,
  %¾«¶È£ºeps£¬Ä¿±êº¯Êý×îСֵʱµÄxÒÔ¼°Ä¿±êº¯Êý×îСֵ
  format long;
  if nargin==4
      eps=1.0e-5;
  end
  S=transpose(f)*f;
  k=length(f);
  n=length(x0);
  x0=transpose(x0);
  tol=1;
  A=jacobian(f,var);
  while tol>eps
      Fx=zeros(k,1);
      for i=1:k
          Fx(i,1)=Funval(f(i),var,x0);
      end
      Sx=Funval(S,var,x0);
      Ax=Funval(A,var,x0);
      gSx=transpose(Ax)*Fx;
      dx=-transpose(Ax)*Fx;
      x0=x0+dx;
      tol=norm(dx);
  end
  x=x0;
  minf=Funval(S,var,x);
  format short;
  toc;
  Iterations;
end
»Ø¸´´ËÂ¥

» ²ÂÄãϲ»¶

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

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

xuetb

Ìú³æ (³õÈëÎÄ̳)

²¹³ä£ºÅ¼ÔÚmatlabÖÐÊäÈ룺 syms x=[x1 x2 x3 x4];f=[x1+10*x2;sqrt(5)*(x3-x4);(x2-2*x3)^2;sqrt(10)*(x1-x4)^2];x0=[3,-1,0,1];[x,minf]=minGN(f,x0,[x],1.0e-5),S,A,k
´íÎó£º Error using ==> syms at 61
Not a valid variable name.
2Â¥2011-12-05 15:54:01
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû

zhfzh

ľ³æ (ÕýʽдÊÖ)

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

¡ï
ÓàÔó³É(½ð±Ò+1): лл²ÎÓëÓ¦Öú£¡»¶Ó­³£À´³ÌÐòÓïÑ԰棡 2011-12-08 14:06:03
Äã¿ÉÒÔÓÃstrcat¡¢  evalº¯Êýͨ¹ýÑ­»·À´ÊµÏÖ¶à¸ö±äÁ¿µÄ¶¨Òå
3Â¥2011-12-05 17:11:33
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû

xuetb

Ìú³æ (³õÈëÎÄ̳)

ÒýÓûØÌû:
3Â¥: Originally posted by zhfzh at 2011-12-05 17:11:33:
Äã¿ÉÒÔÓÃstrcat¡¢  evalº¯Êýͨ¹ýÑ­»·À´ÊµÏÖ¶à¸ö±äÁ¿µÄ¶¨Òå

ÄܰïÖú±àдÏÂô£¿ ÈçºÎʵÏÖ£¿Ð»Ð»À²
4Â¥2011-12-06 00:54:46
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû
Ïà¹Ø°æ¿éÌø×ª ÎÒÒª¶©ÔÄÂ¥Ö÷ xuetb µÄÖ÷Ìâ¸üÐÂ
×î¾ßÈËÆøÈÈÌûÍÆ¼ö [²é¿´È«²¿] ×÷Õß »Ø/¿´ ×îºó·¢±í
[¿¼ÑÐ] 080500Çóµ÷¼Á +6 »ÆÓ 2026-04-06 6/300 2026-04-06 17:41 by ÍÁľ˶ʿÕÐÉú
[¿¼ÑÐ] Çóµ÷¼Á +4 ÐܶþÏëÉϰ¶ 2026-04-06 4/200 2026-04-06 17:17 by dongzh2009
[¿¼ÑÐ] 277Çóµ÷¼Á ÊýÒ»104·Ö +9 Æ¿×ÓPZ 2026-04-05 11/550 2026-04-06 17:07 by À¶ÔÆË¼Óê
[¿¼ÑÐ] 308Çóµ÷¼Á +13 ÌÈÈôÆð·çÁËÄØ 2026-04-05 13/650 2026-04-06 14:20 by ½¯ð©Óí
[¿¼ÑÐ] ר˶0854³õÊÔ¿¼²Ä¿Æ»ù£¬Çóµ÷¼Á +6 3220548044 2026-04-06 9/450 2026-04-06 10:26 by barlinike
[¿¼ÑÐ] ÉúÎïÓëÒ½Ò©273Çóµ÷¼Á +7 ÀóÌâÄÏǽ 2026-04-05 7/350 2026-04-06 09:26 by 286640313
[¿¼ÑÐ] Ò»Ö¾Ô¸Î÷°²½»´ó²ÄÁÏѧ˶£¨Ó¢Ò»Êý¶þ£©347£¬Çóµ÷¼Áµ½¸ß·Ö×Ó/²ÄÁÏÏà¹Ø×¨Òµ +8 zju51 2026-03-31 10/500 2026-04-06 09:03 by ×íÎÌwl
[¿¼ÑÐ] 296Çóµ÷¼Á +3 Íô£¡£¿£¡ 2026-04-05 5/250 2026-04-05 17:38 by À¶ÔÆË¼Óê
[¿¼ÑÐ] Ò»Ö¾Ô¸µç×ӿƼ¼´óѧ085600²ÄÁÏÓ뻯¹¤ 329·ÖÇóµ÷¼Á +10 Naiko 2026-04-04 10/500 2026-04-05 09:40 by sam3303
[¿¼ÑÐ] 085602µ÷¼Á ³õÊÔ×Ü·Ö335 +12 19123253302 2026-04-04 12/600 2026-04-05 08:08 by 544594351
[¿¼ÑÐ] ±¾¿Æ211 ·ÖÊý293ÇëÇóµ÷¼Á +4 Á«²Ë¾ÍÊÇź°É 2026-04-01 4/200 2026-04-04 22:32 by hemengdong
[¿¼ÑÐ] Ò»Ö¾Ô¸¶«±±´óѧ085901ÍÁľר˶345Çóµ÷¼Á +3 zxt11111 2026-04-04 3/150 2026-04-04 14:21 by ÍÁľ˶ʿÕÐÉú
[¿¼ÑÐ] 312Çóµ÷¼Á +4 ÉÞÔÂÉ« 2026-04-02 5/250 2026-04-03 08:21 by fangshan711
[¿¼ÑÐ] 312Çóµ÷¼Á +6 ССī123 2026-04-02 7/350 2026-04-03 07:32 by jsw79
[¿¼ÑÐ] 285Çóµ÷¼Á +8 AZMK 2026-04-02 11/550 2026-04-02 20:16 by yulian1987
[¿¼ÑÐ] »·¾³¹¤³Ì297·ÖÇóµ÷¼ÁÒ»Ö¾Ô¸º¼¸ßÔº +15 GENJIOW 2026-03-31 16/800 2026-04-02 17:56 by cyh¡ª315
[¿¼ÑÐ] 348Çóµ÷¼Á +11 zzzzyk123 2026-04-01 11/550 2026-04-02 16:52 by Wang200018
[¿¼ÑÐ] Ò»Ö¾Ô¸ÉϺ£º£Ñó´óѧ083200ʳƷѧ˶£¬Çóµ÷¼Á£¬½ÓÊÜÆäËûרҵ +6 whatÕÅ 2026-04-01 7/350 2026-04-02 16:48 by zzsw+
[¿¼ÑÐ] Ò»Ö¾Ô¸9³õÊÔ366 ±¾Ë«·ÇÇóµ÷¼Á +4 ÔËÆøÀ´µÃÈôÓÐËÆÎ 2026-04-02 4/200 2026-04-02 09:56 by guanxin1001
[¿¼ÑÐ] ¸´ÊÔµ÷¼Á +7 Ë«ÂíβƦÀϰå2 2026-03-31 7/350 2026-03-31 19:49 by Dyhoer
ÐÅÏ¢Ìáʾ
ÇëÌî´¦ÀíÒâ¼û