4Â¥: Originally posted by dbb627 at 2012-07-10 16:32:23
ÉÏÃæµÄ´úÂ븴ÖÆÖ®ºó£¬´ò¿ªeditor£¬Õ³Ìùºó±£´æΪmulVNewton.mÎļþÔÚµ±Ç°Ä¿Â¼
È»ºóÔÚÃüÁî´°¿ÚÊäÈë=mulVNewton()%x1-x3µÄÈý¸ö³õÖµ

rΪµü´úºóµÄÖµ£¬Èç¹ûÊÕÁ²¾ÍÊǽâ
nµü´ú´ÎÊý...

function [r,n]=mulVNewton(x0,A,eps)
if nargin==1
    A=eye(length(x0));
    eps=1.0e-4;
else
    if nargin==2
        eps=1.0e-4;
    end
end

r=x0-myf(x0)/A;
n=1;
tol=1;
while tol>eps
    x0=r;
    r=x0-myf(x0)/A;
    y=r-x0;
    z=myf(r)-myf(x0);
    A1=A+(z-y*A)'*y/norm(y);
    A=A1;
    n=n+1;
    if(n>100000)
        disp(' µü´ú²½ÊýÌ«¶à£¬¿ÉÄܲ»ÊÕÁ²£¡');
        return;
    end
    tol=norm(r-x0);
end                                       
%------------------
function f=myf(x)

x1=x(1);x2=x(2);x3=x(3);
f=[x1-(x2/81.54-1)*x3*log(1-81.54/x2)+8.80/81.54,    x1-(x2/84.46-1)*x3*log(1-84.46/x2)+9.55/84.46,    x1-(x2/86.76-1)*x3*log(1-86.76/x2)+9.90/86.76];

2Â¥: Originally posted by dbb627 at 2012-07-10 16:09:14
Å£¶Ù·¨

function =mulVNewton(x0,A,eps)
if nargin==1
    A=eye(length(x0));
    eps=1.0e-4;
else
    if nargin==2
        eps=1.0e-4;
    end
end

r=x0-myf(x0)/A;
n=1;
tol=1;
while to ...

3Â¥: Originally posted by ÎÞ³¾²è at 2012-07-10 16:28:18
ÄãÄÜ°ïÎÒ½²Ò»ÏÂÈçºÎ²Ù×÷Âð£¿ÎÒÊÇÒ»¸ö²ËÄñ£¬Ð»Ð»...