A1=[3.01,6.03,1.99;1.27,4.16,-1.23;0.987,-4.81,9.34;];
A2=[3.00,6.03,1.99;1.27,4.16,-1.23;0.990,-4.81,9.34;];
A=[A1,B];
[m,n]=size(A);
x=zeros(1,m);
for j=1:(m-1)
[a,b]=max(abs(A(j:m,j))) ;%找出最大值和位置
b=b+j-1;
if b~=j
A([j,b],:)=A([b,j],:);
end
for i=(j+1):m
A(i,:)=A(i,:)-A(i,j)/A(j,j)*A(j,:);
end
end
for i=(m-1):-1:1
for j=(i+1):m
A(i,:)=A(i,:)-A(i,j)/A(j,j)*A(j,:);
end
end
for i=1:m
x(i)=A(i,n)/A(i,i);
end
A1
B
detA=det(A1)
x
condA_inf=cond(A1,inf)
matlab实验-高斯列主元消去法(草稿)
最新推荐文章于 2021-12-10 20:58:43 发布