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 ...

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];

3楼: Originally posted by 无尘茶 at 2012-07-10 16:28:18
你能帮我讲一下如何操作吗？我是一个菜鸟，谢谢...

4楼: Originally posted by dbb627 at 2012-07-10 16:32:23
上面的代码复制之后，打开editor，粘贴后保存为mulVNewton.m文件在当前目录
然后在命令窗口输入=mulVNewton()%x1-x3的三个初值

r为迭代后的值，如果收敛就是解
n迭代次数...