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

2744557306

这是一个用于解线性方程组的 MATLAB 脚本，采用高斯消元法。下面是对代码的解释：

1.矩阵初始化：
' w+ d$ r6 |+ q8 T
9 u! X+ s  [! j- q+ j) g# v   a1=[2,-2,5;2,3,1;-1,4,1];
   b=[6,13,3]';
   n=length(b); % 方程组大小为 n
   a=zeros(n,n+1);
' r3 J/ P/ f, z5 y, E( V8 H   a(:,1:n)=a1;
( ]) a; X9 E7 s, t( s1 l' Y   a(1:n,n+1)=b; % 增广矩阵

这里给定了系数矩阵 a1 和右侧向量 b，然后初始化了增广矩阵 a。
# J6 |! ^) R% d
2.高斯消元过程：

   for i=1:n-1
       % 确保主对角线元素不为0
( e. h3 S; k; A/ d0 }6 P: M       if(a(i:n,i)==0)
% y, a: W& \  i1 E          error('方程组没有唯一解');
  K' |  n, W8 F: g; P. W       end
       % 选取主对角线元素不为0的行
- z. h* L9 b& j9 |, ]       for p=i:n
! R/ n1 j! N1 F+ P           if(a(p,i)~=0)
$ _2 m: ^8 X1 Z" x               p;
  u$ F0 w% N& _  j. v               break;
           end
. G; B7 [1 l3 _/ X. V  ?! D8 ^- G  Y       end
$ b5 h9 H4 r) P0 M0 y       % 如果选取的行不是当前行，则交换两行
/ n# y! s; L, _       if(p~=i)
         t=a(i,;
         a(i,=a(p,;
0 h1 B  X4 T: ~# b         a(p,=t;
" T' `8 q3 E2 b& r$ n# [! X       end
5 }" e# J2 f% G/ o0 ?# O: b9 J! r       % 高斯消元
       for j=i+1:n
           a(j,i)=a(j,i)/a(i,i);
) \8 S6 b  K4 S           a(j,=a(j,-a(j,i)*a(i,;
       end
  m4 e' B9 O% j$ b2 F5 E   end
   % 检查方程组是否有唯一解
! r1 J3 o9 `. P% ]. {- a   if (a(n,n)==0)
$ `! E5 [( c" H0 c" ?        error('方程组没有唯一解');
   end
/ m+ [' A& B7 x
2 }. U$ y& b3 d/ O) N' p这个部分实现了高斯消元算法。它通过迭代将增广矩阵 a 转化为上三角形式。
* i/ q3 Y! t' }! \/ n* G$ A0 x6 e
3.回代过程：
) k% X5 L  M  m+ i
   % 回代
   x(n)=b(n)/a(n,n);
& W  c7 X! s& R. j! b   for i=n-1:-1:1
       sum=0;
. x6 W+ h4 [4 r       for j=i+1:n
           sum=sum+a(i,j)*x(j);
( k1 Z( V0 u" N       end
       x(i)=(b(i)-sum)/a(i,i);
   end
0 w5 ^7 S9 d( Q9 \9 w& f6 A
这个部分实现了回代过程，得到方程组的解向量 x。
) Q5 v' ?; ~2 X) ^
4.打印结果：

6 K: f8 f4 H% ~   jie=x'

最后，打印求解得到的解向量。
+ {$ _3 `; Q6 l4 X' o在实际应用中，可以使用 MATLAB 提供的 linsolve 函数来更稳妥地解线性方程组。
" j$ E3 }+ d3 K2 ^1 g5 b! d