²é¿´: 709  |  »Ø¸´: 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µÄ»ØÌû
×î¾ßÈËÆøÈÈÌûÍÆ¼ö [²é¿´È«²¿] ×÷Õß »Ø/¿´ ×îºó·¢±í
[¿¼ÑÐ] Ò»Ö¾Ô¸ÖØÇì´óѧ085700×ÊÔ´Óë»·¾³×Ü·Ö308Çóµ÷¼Á +6 īīĮ 2026-03-20 6/300 2026-03-21 09:23 by ÌìÇàµÈÑÌÓê
[¿¼ÑÐ] 268Çóµ÷¼Á +7 ¼òµ¥µã0 2026-03-17 7/350 2026-03-21 08:32 by JourneyLucky
[¿¼ÑÐ] Çóµ÷¼Á +6 Mqqqqqq 2026-03-19 6/300 2026-03-21 08:04 by JourneyLucky
[¿¼ÑÐ] 301Çóµ÷¼Á +10 yyÒªÉϰ¶Ñ½ 2026-03-17 10/500 2026-03-21 03:14 by JourneyLucky
[¿¼ÑÐ] 303Çóµ÷¼Á +5 î£08 2026-03-17 7/350 2026-03-21 03:11 by JourneyLucky
[¿¼ÑÐ] Ò»Ö¾Ô¸ÎäÀí²ÄÁÏ305·ÖÇóµ÷¼Á +6 ÏëÉϰ¶µÄÀðÓã 2026-03-18 7/350 2026-03-21 01:03 by JourneyLucky
[¿¼ÑÐ] Ò»Ö¾Ô¸ÄϾ©Àí¹¤´óѧ085701×ÊÔ´Óë»·¾³302·ÖÇóµ÷¼Á +4 ¿ûè÷ÎÀ¶Ó 2026-03-18 6/300 2026-03-20 23:02 by JourneyLucky
[¿¼ÑÐ] ¿¼Ñе÷¼ÁÇóÑ§Ð£ÍÆ¼ö +3 ²®ÀÖ29 2026-03-18 5/250 2026-03-20 22:59 by JourneyLucky
[¿¼ÑÐ] 324Çóµ÷¼Á +5 luckyѽѽѽѼ 2026-03-20 5/250 2026-03-20 22:30 by ´ÙÌì³É
[¿¼ÑÐ] 317Çóµ÷¼Á +5 Éê×ÓÉêÉê 2026-03-19 9/450 2026-03-20 22:26 by JourneyLucky
[¿¼ÑÐ] Ò»Ö¾Ô¸ËÕÖÝ´óѧ²ÄÁÏÇóµ÷¼Á£¬×Ü·Ö315£¨Ó¢Ò»£© +5 sbdksD 2026-03-19 5/250 2026-03-20 22:10 by luoyongfeng
[¿¼ÑÐ] Ò»Ö¾Ô¸ÖÐÄÏ»¯Ñ§£¨0703£©×Ü·Ö337Çóµ÷¼Á +8 niko- 2026-03-19 9/450 2026-03-20 21:57 by luoyongfeng
[¿¼ÑÐ] ÕÐÊÕµ÷¼Á˶ʿ +4 lidianxing 2026-03-19 12/600 2026-03-20 12:25 by lidianxing
[¿¼²©] ÕÐÊÕ²©Ê¿1-2ÈË +3 QGZDSYS 2026-03-18 3/150 2026-03-20 11:58 by ßÉßÉßÉßɽÐ
[¿¼ÑÐ] 085410È˹¤ÖÇÄÜר˶317Çóµ÷¼Á£¨0854¶¼¿ÉÒÔ£© +4 xbxudjdn 2026-03-18 4/200 2026-03-20 09:07 by ²»168
[¿¼ÑÐ] ²ÄÁÏר˶326Çóµ÷¼Á +6 Ä«ìÏæ¦Ý· 2026-03-15 7/350 2026-03-17 17:10 by ruiyingmiao
[¿¼ÑÐ] Ò»Ö¾Ô¸£¬¸£ÖÝ´óѧ²ÄÁÏר˶339·ÖÇóµ÷¼Á +3 ľ×ÓmomoÇàÕù 2026-03-15 3/150 2026-03-17 07:52 by laoshidan
[¿¼ÑÐ] 326Çóµ÷¼Á +4 ŵ±´¶û»¯Ñ§½±êéê 2026-03-15 7/350 2026-03-16 17:11 by ŵ±´¶û»¯Ñ§½±êéê
[¿¼ÑÐ] 070305Çóµ÷¼Á +3 mlpqaz03 2026-03-14 4/200 2026-03-15 11:04 by peike
[¿¼ÑÐ] ¸´ÊÔµ÷¼Á +3 ºôºô£¿~+123456 2026-03-14 3/150 2026-03-14 16:53 by WTUChen
ÐÅÏ¢Ìáʾ
ÇëÌî´¦ÀíÒâ¼û