列主元高斯消去法matlab_matlab 用高斯列主元法求解-CSDN博客

本文链接：https://blog.csdn.net/Grace_five/article/details/134755637

本文介绍了使用列主元策略改进的高斯消元法在Matlab中的实现，重点讨论了如何选择绝对值最大的元素作为主元，以及如何处理零元素和避免“大数吃小数”现象，提供了一个完整的Matlab函数示例。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

列主元高斯消去法（matlab）

1.相较于顺序消元法的优点
① 避免了顺序消元法一旦出现aii=0，不能进行计算的缺点；
②避免了“大数吃小数”的现象，数值稳定。

准备：
① 列主元的主要思想是在顺序消元法的基础上，实现在每次消元前，将第k行至n行绝对值 的最大元素选为主元akk，然后实施消元。
②matlab中的指令：abs（*）求绝对值；
max最值函数主要有以下指令：

M = max(A) %求矩阵每列的最大值；
M = max(A,[],dim) %dim=1时，求每列最大元；dim=2时，求每行最大元；
[M,I] = max(*)%M为列最大元；I为最大元位置；
-C = max(A,B) C为A和B中的最大元素.
③对m=index+k-2 的理解（把相对位置转换成绝对位置）（原因：max(abs(b(k-1:n,k-1)))求的时相对位置的index）
max(abs(b(k-1:n,k-1))）从第k-1行开始比较大小，因此返回的index是少了n-(k）行进行大小比较，index+k-1两个相加再减去重复的一行。（自己能看懂的解释）

以上准备就绪，开始写程序:

function X=gaosi2(A,cf)
%用高斯消元法求解方程组
%A表示方程组系数矩阵，b为右端列向量
%x为解向   问题：变量之间混乱 
%% 例子
A=[1 1 1;-1 3 1;2 -6 1];  cf=[6;4;-5];
%% 消元
n=length(A);a=[A,cf];b=a;
for k=2:n
    [~,index]=max(abs(b(k-1:n,k-1)))  %找出绝对值最大数的相对位置
   m=index+k-2 % 把相对位置转换成绝对位置！！！绝对=相对于任何（加到if前面！在else后面不行)
    if a(m,k-1)==a(k-1,k-1)
       b=b;
    else
       z=b(k-1,:)
       b(k-1,:)=b(m,:)
       b(m,:)=z  
   end
   for i=k:n
   b(i,k-1:n+1)=b(i,k-1:n+1)-b(k-1,k-1:n+1).*(b(i,k-1)/b(k-1,k-1)) ; 
   end
end   
%% 回代
X=zeros(n,1); 
X(n,1)=b(n,n+1)/b(n,n); %h放在这里的话会出错，要放在第一个循环后面
%此下求解分为两部分，需要写入一个求和
for i=(n-1):-1:1
    h=0;
    for j=i+1:n
      h=b(i,j)*X(j,1)+h;
      h;
    end
   X(i,1)=(b(i,n+1)-h)/b(i,i);
   X
end