程序运行后全是NAN，问题出在哪

s20071041

function [B,w]=wu(M,K,eps)
D=inv(K)*M;      
n=rank(D);
A0=ones(n,1);
A=D*A0;
c=0;
while c<=n
    if norm(A-A0)>=eps
    A0=A;
    A=D*A0;
    else
    c=c+1;
    w(c,1)=1/A(n,1);
    A0=A*w(c,1);
    B(:,c)=A0(:,1);
    end
end

调用后
>> M=[2 0 0;0 1.5 0; 0 0 1];
>> K=[5 -2 0;-2 3 -1;0 -1 1];
>> eps=1e-6;
>> wu(M,K,eps)

ans =

     0     0     0     0
   NaN   NaN   NaN   NaN
   NaN   NaN   NaN   NaN

[ 本帖最后由 eight 于 2007-12-4 22:22 编辑 ]

sigma665

function [B,w]=wu(M,K,eps)

这个是什么

s20071041

前面是编写的一个函数呀  后面调用他

jinqiankun

NAN 的意思是非数 即0/0 或无穷大/无穷大
你 看看你的程序 是不是 有这种可能？

eight

原帖由 s20071041 于 2007-12-4 21:25 发表

                               
登录/注册后可看大图

function =wu(M,K,eps)
D=inv(K)*M;      
n=rank(D);
A0=ones(n,1);
A=D*A0;
c=0;
while c=eps
    A0=A;
    A=D*A0;
    else
    c=c+1;
    w(c,1)=1/A(n,1);
    A0=A*w(c,1);
    B(:,c)= ...

像这种问题自己设置断点调试一下即可

s20071041

应该有结果  我是书上例子   我用程序解他  答案没出来  应该程序出问题了

eight

原帖由 s20071041 于 2007-12-4 22:35 发表

                               
登录/注册后可看大图

应该有结果  我是书上例子   我用程序解他  答案没出来  应该程序出问题了

估计课本没有问题，你先看看是否打字的时候输入有误

s20071041

课本是手算的结果  我们要求编程序计算  还好编出来了   是程序问题  谢谢你了

花如月

希望可以分享一下你的程序