Znn3bq.jpeg
²é¿´: 715  |  »Ø¸´: 3
µ±Ç°Ö»ÏÔʾÂú×ãÖ¸¶¨Ìõ¼þµÄ»ØÌû£¬µã»÷ÕâÀï²é¿´±¾»°ÌâµÄËùÓлØÌû

ÓÆÈ»Çå·ç

гæ (³õÈëÎÄ̳)

[ÇóÖú] matlab¸ßάº¯ÊýÈçºÎµ÷Óõ½Ö÷º¯ÊýÖУ¿ ÒÑÓÐ1È˲ÎÓë

%×îËÙϽµ·¨
%×îËÙϽµ·¨
function [x,k,tol,fmin]=minzs(fun,gfun,x0,eps)
format long;
tic
if nargin==3
    eps=1.0e-6;
end
tol=1;
syms r;
k=0;
while tol>eps
   g0=feval(gfun,x0);
    tol=norm(g0);
    m=x0+r*g0;
    yf=feval(fun,m);
    r1=minnd(yf,0);
    x1=x0+r1*g0;
    x0=x1;
    k=k+1;
end
x=vpa(x0,5);
fmin=feval(fun,x0);
fmin=vpa(fmin,5);
format short;


fun1.mº¯Êý
function f=fun1(x)
y=ones(1,20);
var=vpa(ones(1,20));
for i=1:20;
    x(i)=vpa([\'x\',num2str(i)]);
    var(i)=x(i);
    y(i)=1.0;
end

m=0;
for i=1:10
    m=m+(x(2*i-1))^2+100*(x(2*i))^2;
end

    f=0.5*m;

gfun1º¯Êý
function gf=gfun1(x)
var=vpa(ones(1,20));
for i=1:20;
   x(i)=vpa([\'x\',num2str(i)]);
    var(i)=x(i);
    if mod(i,2)==0
        x(i)=100*x(i);
    end
end
gf=x;



clear
x0=ones(1,20);
[x,k,tol,fmin]=minzs(\'fun\',\'gfun\',x0)
Ϊʲô²»ÄÜÔËÐÐ
The following error occurred converting from sym to double:
Error using mupadmex
Error in MuPAD command: DOUBLE cannot convert the input expression
into a double array.

If the input expression contains a symbolic variable, use the VPA
function instead.

Error in gfun1 (line 4)
   x(i)=vpa([\'x\',num2str(i)]);

Error in minzs (line 12)
   g0=feval(gfun,x0);
»Ø¸´´ËÂ¥

» ²ÂÄãϲ»¶

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

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

ÓÆÈ»Çå·ç

гæ (³õÈëÎÄ̳)

ÒýÓûØÌû:
3Â¥: Originally posted by wurongjun at 2016-03-25 17:52:35
ÄãµÄº¯Êýfun1д·¨²»ÕýÈ·: x(i)=vpa();Õâ¿éÓÐÎÊÌâ!

лл£¬ÎÒÒѾ­½â¾öÁË¡£
4Â¥2016-03-26 10:36:51
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû
²é¿´È«²¿ 4 ¸ö»Ø´ð

ÓÆÈ»Çå·ç

гæ (³õÈëÎÄ̳)

function [x,fmin]=minnd(f,x0,eps)
format long;
if nargin==2
  eps=1.0e-6;
end
g1=diff(f);
g2=diff(g1);
k=0;
t=1;
while t>eps
h1=subs(g1,findsym(g1),x0);
h2=subs(g2,findsym(g2),x0);
t=abs(h1);
x0=x0-h1/h2;
k=k+1;
end
x=x0;
fmin=subs(f,findsym(f),x);
format short;
2Â¥2016-03-25 16:09:53
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû

wurongjun

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

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

¡ï ¡ï ¡ï ¡ï ¡ï
¸Ðл²ÎÓ룬ӦÖúÖ¸Êý +1
ÓÆÈ»Çå·ç: ½ð±Ò+5, »¹ÊÇлл£¬ÎÒÒѾ­½â¾öÁË 2016-03-26 18:12:32
ÄãµÄº¯Êýfun1д·¨²»ÕýÈ·: x(i)=vpa([\'x\',num2str(i)]);Õâ¿éÓÐÎÊÌâ!
ÉÆ¶ñµ½Í·ÖÕÓб¨,È˼äÕýµÀÊDz×É£.
3Â¥2016-03-25 17:52:35
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû
×î¾ßÈËÆøÈÈÌûÍÆ¼ö [²é¿´È«²¿] ×÷Õß »Ø/¿´ ×îºó·¢±í
[¿¼ÑÐ] 085404 293Çóµ÷¼Á +9 ÓÂÔ¶¿â°®314 2026-04-08 9/450 2026-04-12 02:24 by Çï¶¹²ËÑ¿
[¿¼ÑÐ] 270Çóµ÷¼Á +14 ÑîÀÖ369 2026-04-11 14/700 2026-04-11 20:16 by À¶ÔÆË¼Óê
[¿¼ÑÐ] Çóµ÷¼Á +11 ôæôæÒ»ÊéÉú 2026-04-09 11/550 2026-04-11 19:57 by ÄæË®³Ë·ç
[¿¼ÑÐ] 295Çóµ÷¼Á +6 ?ÒªÉϰ¶? 2026-04-05 7/350 2026-04-11 19:02 by laoshidan
[¿¼ÑÐ] ũѧ0904 312Çóµ÷¼Á +3 Say Never 2026-04-11 3/150 2026-04-11 17:22 by daydayup2005
[¿¼ÑÐ] Çóµ÷¼Á +13 ÕÅ·¬ÇѲ»³´µ° 2026-04-10 14/700 2026-04-11 16:02 by vgtyfty
[¿¼ÑÐ] »¯Ñ§308·ÖÇóµ÷¼Á +22 ÄãºÃÃ÷ÌìÄãºÃ 2026-04-07 24/1200 2026-04-11 11:14 by ChemPharm
[¿¼ÑÐ] 085600²ÄÁÏÓ뻯¹¤329·ÖÇóµ÷¼Á +16 Ò¶zilin 2026-04-10 16/800 2026-04-11 11:04 by may_ÐÂÓî
[¿¼ÑÐ] Çóµ÷¼Á +6 archer.. 2026-04-09 8/400 2026-04-11 10:55 by zhq0425
[¿¼ÑÐ] 300·ÖÇóµ÷¼Á £¨085501»úеר˶£¬±¾¿ÆÑï´ó£© +8 xu@841019 2026-04-11 8/400 2026-04-11 10:46 by qingpingzhu
[¿¼ÑÐ] ¼ÆËã»ú11408£¬286·ÖÇóµ÷¼Á +10 ľ×ÓÄî•„ 2026-04-05 10/500 2026-04-11 10:44 by zhq0425
[¿¼ÑÐ] Ò»Ö¾Ô¸985»úеѧ˶380Çóµ÷¼Á +5 ¹Ø¹ØöÂð¯10 2026-04-11 5/250 2026-04-11 10:10 by ÖªÄî¡£A
[¿¼ÑÐ] ¸´ÊÔµ÷¼Á +9 »ý¼«ÏòÉÏ£» 2026-04-10 11/550 2026-04-11 09:25 by Öí»á·É
[¿¼ÑÐ] 336²ÄÁÏÓ뻯¹¤085600Çóµ÷¼Á +21 Ë®ÐǼÇinfp 2026-04-05 24/1200 2026-04-10 15:28 by luoyongfeng
[¿¼ÑÐ] Çóµ÷¼Á +11 ôæôæÒ»ÊéÉú 2026-04-09 13/650 2026-04-10 10:27 by liuhuiying09
[¿¼ÑÐ] Çóµ÷¼Á +8 ³Ô¿Ú±ù¼¤Áè 2026-04-07 8/400 2026-04-09 08:03 by 5268321
[¿¼ÑÐ] ²ÄÁϹ¤³ÌרҵÈÕÓïÉúÇóµ÷¼Á +9 111623 2026-04-07 9/450 2026-04-07 23:31 by Ò»Ö»ºÃ¹û×Ó?
[¿¼ÑÐ] ²ÄÁϵ÷¼Á +5 СÁõͬѧ߹߹ 2026-04-06 5/250 2026-04-06 18:34 by sherry_1901
[¿¼ÑÐ] ¹¤¿Æ370Çóµ÷¼Á +3 äçÐļ弦µ° 2026-04-05 3/150 2026-04-06 10:55 by ÕâÊÇÒ»¸öÎÞÁĵÄê
[¿¼ÑÐ] 0857´óÀà»·¾³¹¤³ÌBÇøÇóµ÷¼Á +3 ¹¨ÓíÃú 2026-04-05 3/150 2026-04-06 10:22 by À¶ÔÆË¼Óê
ÐÅÏ¢Ìáʾ
ÇëÌî´¦ÀíÒâ¼û