Znn3bq.jpeg
²é¿´: 714  |  »Ø¸´: 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µÄ»ØÌû

ÓÆÈ»Çå·ç

гæ (³õÈëÎÄ̳)

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µÄ»ØÌû

ÓÆÈ»Çå·ç

гæ (³õÈëÎÄ̳)

ÒýÓûØÌû:
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µÄ»ØÌû
Ïà¹Ø°æ¿éÌø×ª ÎÒÒª¶©ÔÄÂ¥Ö÷ ÓÆÈ»Çå·ç µÄÖ÷Ìâ¸üÐÂ
×î¾ßÈËÆøÈÈÌûÍÆ¼ö [²é¿´È«²¿] ×÷Õß »Ø/¿´ ×îºó·¢±í
[¿¼ÑÐ] Çóµ÷¼Á£¬262»úеר˶ +8 àÅyyl 2026-04-08 8/400 2026-04-12 02:31 by Çï¶¹²ËÑ¿
[¿¼ÑÐ] µ÷¼Á +6 ÔÂ@163.com 2026-04-11 6/300 2026-04-11 21:46 by lwk2004
[¿¼ÑÐ] 283Çóµ÷¼Á 086004¿¼Ó¢¶þÊý¶þ +17 ÄǸöàà×Ó 2026-04-10 18/900 2026-04-11 16:27 by Ã÷Ô´ËʱÓÐ
[¿¼ÑÐ] 308Çóµ÷¼Á +3 VvvvL 2026-04-10 3/150 2026-04-11 15:32 by xiangxu
[»ù½ðÉêÇë] ɽ¶«Ê¡»ù½ð2026 +4 jerry681 2026-04-08 5/250 2026-04-11 13:59 by laobibibi
[¿¼ÑÐ] 0858Çóµ÷¼Á 5+5 Gky09300550£¬ 2026-04-10 8/400 2026-04-10 19:13 by chemisry
[¿¼ÑÐ] 314Çóµ÷¼Á +23 wakeluofu 2026-04-09 24/1200 2026-04-10 15:31 by MOF_Catal
[¿¼ÑÐ] ³õÊÔ·Ö332£¬Ò»Ö¾Ô¸±¨¿¼Î÷±±¹¤Òµ´óѧ£¬ +11 ¹ÊÈË?? 2026-04-09 11/550 2026-04-09 21:54 by JineShine
[¿¼ÑÐ] ²ÄÁÏר˶³õÊÔ·Ö332Ò»Ö¾Ô¸Î÷±±¹¤Òµ´óѧ£¬ +12 ¹ÊÈË?? 2026-04-09 12/600 2026-04-09 18:34 by Ccclqqq
[¿¼ÑÐ] ÉúÎïѧ308·ÖÇóµ÷¼Á£¨Ò»Ö¾Ô¸»ª¶«Ê¦´ó£© +13 ÏàÐűػá¹ââÍòÕ 2026-04-06 16/800 2026-04-09 13:54 by ÐìÁ¼°×ü´óÏÀ
[¿¼ÑÐ] ²ÄÁϹ¤³Ì322 +18 ¹þ¹þ¹þºðºðºð¹þ 2026-04-07 19/950 2026-04-09 10:44 by cymywx
[¿¼ÑÐ] 328Çóµ÷¼Á +17 lftmya 2026-04-07 18/900 2026-04-09 08:05 by 5268321
[¿¼ÑÐ] 318Çóµ÷¼Á +13 ykyhsa 2026-04-05 15/750 2026-04-08 21:37 by wj165256
[¿¼ÑÐ] 326·Ö£¬Ò»Ö¾Ô¸»¦9£¬ÇóÉúÎïѧµ÷¼Á +4 Áõīī 2026-04-05 4/200 2026-04-08 06:22 by lijunpoly
[¿¼ÑÐ] 305Çóµ÷¼Á +4 77Qi 2026-04-06 4/200 2026-04-07 20:06 by shanqishi
[¿¼ÑÐ] 319·Ö085702°²È«¹¤³ÌÇóµ÷¼Á +6 rious 2026-04-05 6/300 2026-04-07 09:42 by jp9609
[¿¼ÑÐ] ²ÄÁϵ÷¼Á +5 СÁõͬѧ߹߹ 2026-04-06 5/250 2026-04-06 18:34 by sherry_1901
[¿¼ÑÐ] 285Çóµ÷¼Á +5 mapmath 2026-04-06 6/300 2026-04-06 17:18 by À¶ÔÆË¼Óê
[¿¼ÑÐ] 085405Èí¼þ¹¤³Ì301·ÖÇóµ÷¼Á£¬×¨Ë¶¿É¿çרҵ£¬ËÄÁù¼¶Òѹý +3 ¾²¾²ÏëÏë 2026-04-05 3/150 2026-04-06 15:23 by nepu_uu
[¿¼ÑÐ] ¿¼Ñе÷¼ÁÉúѰÕÒµ¼Ê¦ +3 ¹ËÕ°¿¼Ñа¡ 2026-04-05 3/150 2026-04-05 18:18 by à£à£à£0119
ÐÅÏ¢Ìáʾ
ÇëÌî´¦ÀíÒâ¼û