用matlab程序法编出高斯消元法
2017-06-20 分类:百科
TIPS:本文共有 429 个字,阅读大概需要 1 分钟。
function[x,XA]=GaussXQByOrder(A,b)%高斯顺序消元法N=size(A)n=N(1)
fori=1:(n-1)forj=(i+1):nif(A(i,i)==0)disp('对角元素为0!')
%防止对角元素为0returnendl=A(j,i)m=A(i,i)A(j,1:n)=A(j,1:n)-l*A(i,1:n)/m%消元方程b(j)=b(j)-l*b(i)/mendendx=SolveUpTriangle(A,b)
%通用的求上三角系数矩阵线性方程组的函数XA=A%消元后的系数矩阵functionx=SolveUpTriangle(A,b)N=size(A)
n=N(1)
fori=n:-1:1if(i<n)s=A(i,(i+1):n)*x((i+1):n,1)
elses=0endx(i,1)=(b(i)-s)/A(i,i)
end程序如上,自己算吧。这个东西应该尝试自己编程 参考资料: matlab常用算法
小编精心整理的这篇内容:用matlab程序法编出高斯消元法,如果你看到此处请一定要收藏哦!
阅读剩余内容