列主元高斯消去法matlab

列主元高斯消去法（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

给个栗子：