论坛 › Matlab讨论区 › 这个高斯方程（不是列主元消去法）程序有什么问题？谢谢

cannor 发表于 2006-11-16 20:06

这个高斯方程（不是列主元消去法）程序有什么问题？谢谢

为什么一直算出来，跟列主元消去法的一样？



function x=nomal1(A,b)
%x=nomal1(A,b)
n=length(A);
a=;
for k=1:n-1
    if a(k,k)==0
      return;
    end
    for i=k+1:n
      l(i,k)=a(i,k)/a(k,k);
      a(i,k+1:n+1)=a(i,k+1:n+1)-l(i,k).*a(k,k+1:n+1);
    end
end

%回代
if a(n,n)==0
    return
end
x(n)=a(n,n+1)/a(n,n);
for i=n-1:-1:1
    x(i)=(a(i,n+1)-sum(a(i,i+1:n).*x(i+1:n)))/a(i,i);
end

branch 发表于 2006-11-16 21:56

第一列有做消去吗??

cannor 发表于 2006-11-16 23:22

根据算法第一列好像应该不用消去吧

cannor 发表于 2006-11-16 23:29

哦不好意思 我看出第一行了   非常谢谢！！！！！

查看完整版本: 这个高斯方程（不是列主元消去法）程序有什么问题？谢谢