matlab高斯列主元消去法
合集下载
相关主题
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
%创建高斯列主元素消去法函数
function X=liezhuyuan(A,b)
B=[A b'];
n=length(b);
X=zeros(n,1);
C=zeros(1,n+1);
for k=1:n-1
[Y,j]=max(abs(B(k:n,k)));
C=B(k,:);
B(k,:)=B(j+k-1,:);
B(j+k-1,:)=C;
for i=k+1:n
m=B(i,k)/B(k,k);
B(i,k:n+1)=B(i,k:n+1)-m*B(k,k:n+1);
end
end
b=B(1:n,n+1);
A=B(1:n,1:n);
X(n)=b(n)/A(n,n);
for k=n-1:-1:1
X(k)=(b(k)-sum(A(k,k+1:n)*X(k+1:n)))/A(k,k); End