MATLAB实现列主元消去法求解任意线性方程组的解

%MATLAB实现列主元消去法求解任意线性方程组的解 clear all close all clc format longe disp('请输入参数'); K=input('维数n='); tic A=100*rand(K);% A元素是0-100 for i=1:K A(i,i)=sum(abs(A(i,:)))+25*rand(1); %对角占优的量为0~25 end b=zeros(K,1); for i=1:K; x=0; for r=1:K; x=x+A(i,r); end b(i,1)=x; end%产生b矩阵,b中的元素为A中对应行的和,目的是使方程解全为 1 jd=input('控制精度jd='); A=[A b]; %矩阵增广 for k=1:K-1 B=A(k:K,k); %先将第k列可能作为主元的元素取出方至矩阵B P=max(abs(B)); %%选主元P if(P