采用高斯消元法解线性方程组

2744557306

这是一个用于解线性方程组的 MATLAB 脚本，采用高斯消元法。下面是对代码的解释：
5 o* x# m/ G# J, _( w( E* n- W
1.矩阵初始化：

/ |8 ~9 h" U! s* g0 C, L! I. N   a1=[2,-2,5;2,3,1;-1,4,1];
( \: A8 m9 L" n: A   b=[6,13,3]';
   n=length(b); % 方程组大小为 n
   a=zeros(n,n+1);
   a(:,1:n)=a1;
   a(1:n,n+1)=b; % 增广矩阵
6 M1 X$ l0 x- u1 o
这里给定了系数矩阵 a1 和右侧向量 b，然后初始化了增广矩阵 a。
+ K, m5 N' ^  O$ X/ A2 ?- E
2.高斯消元过程：

   for i=1:n-1
       % 确保主对角线元素不为0
* g7 L% _( j$ H% m' I       if(a(i:n,i)==0)
          error('方程组没有唯一解');
       end
5 S6 p& e1 \( W0 {) e8 v: d       % 选取主对角线元素不为0的行
- t' }0 U& D  k2 Z& ?+ J! a- B0 N       for p=i:n
. ^" o; F& W* M) D: H! T6 U/ }; d           if(a(p,i)~=0)
7 o1 v! p' Y: k  X( u               p;
               break;
; Z" M: c# s5 u! b+ Y3 J5 W           end
- D& N: x* {( D$ t! U$ P- V+ W       end
5 O- v/ G1 G1 R7 h$ S       % 如果选取的行不是当前行，则交换两行
       if(p~=i)
  |  }5 Z4 ]/ B" h: e' @         t=a(i,;
8 [' E& O0 M: t* x( Z) j. m         a(i,=a(p,;
: i) M$ u( w& Z         a(p,=t;
       end
       % 高斯消元
       for j=i+1:n
- f4 C  m5 c; Y1 u/ c5 A9 Y           a(j,i)=a(j,i)/a(i,i);
           a(j,=a(j,-a(j,i)*a(i,;
       end
) a3 [( a6 \& }# g. y& L% Q   end
" E1 M3 Q# Y3 J2 S$ b# V   % 检查方程组是否有唯一解
   if (a(n,n)==0)
        error('方程组没有唯一解');
1 y3 h$ e% H) b  {8 \2 }   end

这个部分实现了高斯消元算法。它通过迭代将增广矩阵 a 转化为上三角形式。
6 k+ l# p% _: K% t) a2 G3 H1 n3 a
3.回代过程：
9 d/ G7 L2 |, \$ |% p
! S, y. s( [5 v   % 回代
% M  ?$ ^( H* I7 g' g& Q2 J   x(n)=b(n)/a(n,n);
7 O( E' D& @" `   for i=n-1:-1:1
       sum=0;
  Q- X) {. n: v9 Y6 m& k       for j=i+1:n
) l# o! U4 W& O# C8 }           sum=sum+a(i,j)*x(j);
       end
       x(i)=(b(i)-sum)/a(i,i);
   end

7 z  i4 ^) i" d3 j5 C# ^这个部分实现了回代过程，得到方程组的解向量 x。
4 M/ Y5 Q% o- l$ U% `5 D
4.打印结果：
$ M4 l/ r  k3 P& i/ K
) z: _8 H, z) [1 B  W( d9 k  c   jie=x'

最后，打印求解得到的解向量。
在实际应用中，可以使用 MATLAB 提供的 linsolve 函数来更稳妥地解线性方程组。