²é¿´: 702  |  »Ø¸´: 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µÄ»ØÌû
Ïà¹Ø°æ¿éÌø×ª ÎÒÒª¶©ÔÄÂ¥Ö÷ ÓÆÈ»Çå·ç µÄÖ÷Ìâ¸üÐÂ
×î¾ßÈËÆøÈÈÌûÍÆ¼ö [²é¿´È«²¿] ×÷Õß »Ø/¿´ ×îºó·¢±í
[¿¼ÑÐ] 278Çóµ÷¼Á +6 ÑÌ»ðÏÈÓÚ´º 2026-03-17 6/300 2026-03-21 01:57 by JourneyLucky
[¿¼ÑÐ] ²ÄÁÏ 336 Çóµ÷¼Á +3 An@. 2026-03-18 4/200 2026-03-21 01:39 by JourneyLucky
[¿¼ÑÐ] Ò»Ö¾Ô¸»ªÖпƼ¼´óѧ£¬080502£¬354·ÖÇóµ÷¼Á +5 ÊØºòϦÑôCF 2026-03-18 5/250 2026-03-21 01:06 by JourneyLucky
[¿¼ÑÐ] 330Çóµ÷¼Á +4 С²Ä»¯±¾¿Æ 2026-03-18 4/200 2026-03-20 23:13 by JourneyLucky
[¿¼ÑÐ] 287Çóµ÷¼Á +7 ³¿»èÏßÓëÐǺ£ 2026-03-19 8/400 2026-03-20 22:19 by JourneyLucky
[¿¼ÑÐ] Ò»Ö¾Ô¸ËÕÖÝ´óѧ²ÄÁÏÇóµ÷¼Á£¬×Ü·Ö315£¨Ó¢Ò»£© +5 sbdksD 2026-03-19 5/250 2026-03-20 22:10 by luoyongfeng
[¿¼ÑÐ] 295¸´ÊÔµ÷¼Á +8 ¼òľChuFront 2026-03-19 8/400 2026-03-20 20:44 by zhukairuo
[¿¼ÑÐ] 0856µ÷¼Á£¬ÊÇѧУ¾ÍÈ¥ +8 sllhht 2026-03-19 9/450 2026-03-20 14:25 by ÎÞи¿É»÷111
[¿¼ÑÐ] ²ÄÁÏѧ˶318Çóµ÷¼Á +5 February_Feb 2026-03-19 5/250 2026-03-19 23:51 by 23Postgrad
[¿¼ÑÐ] Ò»Ö¾Ô¸Öйúº£Ñó´óѧ£¬ÉúÎïѧ£¬301·Ö£¬Çóµ÷¼Á +5 1ËïÎò¿Õ 2026-03-17 6/300 2026-03-19 23:46 by zcl123
[¿¼ÑÐ] ±¾¿ÆÖ£ÖÝ´óѧÎïÀíѧԺ£¬Ò»Ö¾Ô¸»ª¿Æ070200ѧ˶£¬346Çóµ÷¼Á +4 ÎÒ²»ÊÇÒ»¸ù´Ð 2026-03-18 4/200 2026-03-19 09:11 by ¸¡ÔÆ166
[¿¼ÑÐ] ²ÄÁϹ¤³Ìר˶µ÷¼Á +5 204818@lcx 2026-03-17 6/300 2026-03-18 22:55 by 204818@lcx
[¿¼ÑÐ] »¯Ñ§¹¤³Ì321·ÖÇóµ÷¼Á +15 ´óÃ×·¹£¡ 2026-03-15 18/900 2026-03-18 14:52 by haxia
[¿¼ÑÐ] 311Çóµ÷¼Á +6 26ÑÐ0 2026-03-15 6/300 2026-03-18 14:43 by haxia
[¿¼ÑÐ] 312Çóµ÷¼Á +8 İå·Ï£ 2026-03-16 9/450 2026-03-18 12:39 by Linda Hu
[¿¼ÑÐ] 299Çóµ÷¼Á +5 ¡÷С͸Ã÷* 2026-03-17 5/250 2026-03-18 11:49 by ¾¡Ë´Ò¢1
[¿¼ÑÐ] ÉúÎïѧ071000 329·ÖÇóµ÷¼Á +3 ÎÒ°®ÉúÎïÉúÎﰮΠ2026-03-17 3/150 2026-03-18 10:12 by macy2011
[¿¼ÑÐ] 304Çóµ÷¼Á +4 ahbd 2026-03-14 4/200 2026-03-16 16:48 by ÎҵĴ¬Îҵĺ£
[¿¼ÑÐ] ÇóÀÏʦÊÕÁôµ÷¼Á +4 jiang½ª66 2026-03-14 5/250 2026-03-15 20:11 by Winj1e
[¿¼ÑÐ] 085601²ÄÁϹ¤³Ì315·ÖÇóµ÷¼Á +3 yang_0104 2026-03-15 3/150 2026-03-15 10:58 by peike
ÐÅÏ¢Ìáʾ
ÇëÌî´¦ÀíÒâ¼û