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

2744557306

这是一个用于解线性方程组的 MATLAB 脚本，采用高斯消元法。下面是对代码的解释：
/ {8 ]! _" I1 |+ k+ X8 q5 c
  p" [: H- K+ \8 y/ C1.矩阵初始化：

   a1=[2,-2,5;2,3,1;-1,4,1];
: @& C1 ?: j- u7 t# ^) H   b=[6,13,3]';
% b( w1 B2 B/ P! P1 H   n=length(b); % 方程组大小为 n
$ M2 V& A, k$ M2 j; a# H   a=zeros(n,n+1);
1 U8 o7 @7 ?& S5 t# w& t   a(:,1:n)=a1;
( F* |0 @7 a" A5 {  I   a(1:n,n+1)=b; % 增广矩阵
) G% I+ t. a: {" d
9 ~' `5 H8 Y2 b这里给定了系数矩阵 a1 和右侧向量 b，然后初始化了增广矩阵 a。
9 z. B: [* X0 R) G2 r, L6 E
5 U; t: Y; h- T2.高斯消元过程：
) }& o) \7 A$ q
8 v( ~/ s( _  V" C   for i=1:n-1
7 r' C! i" Y$ v       % 确保主对角线元素不为0
       if(a(i:n,i)==0)
5 T" q7 g+ A! J: y; h          error('方程组没有唯一解');
       end
       % 选取主对角线元素不为0的行
# X) |+ O: m$ p       for p=i:n
- A# x; ~; x$ ~# P$ I3 U) V- \           if(a(p,i)~=0)
               p;
4 r! t, }) e) q4 ]! t2 B/ F$ b0 Q+ K, U               break;
) V  U* O; `1 X' z2 i" o$ a. G/ {           end
       end
* a' u  _8 f# M# w$ D3 h% n       % 如果选取的行不是当前行，则交换两行
( L+ v. c# `7 a       if(p~=i)
         t=a(i,;
         a(i,=a(p,;
         a(p,=t;
       end
       % 高斯消元
       for j=i+1:n
           a(j,i)=a(j,i)/a(i,i);
           a(j,=a(j,-a(j,i)*a(i,;
$ w3 U7 o& B. Q: M       end
6 a) v+ d  Q) _- R! l2 M/ m   end
# o3 Q" a# W2 N' @8 R7 G# w5 l   % 检查方程组是否有唯一解
   if (a(n,n)==0)
9 Y1 I( i, n8 N5 E. x        error('方程组没有唯一解');
   end

这个部分实现了高斯消元算法。它通过迭代将增广矩阵 a 转化为上三角形式。
  ]! A2 }' O8 D+ A
3.回代过程：

   % 回代
   x(n)=b(n)/a(n,n);
' F: V  L* G' {+ t$ A5 S$ o   for i=n-1:-1:1
       sum=0;
       for j=i+1:n
, ?4 o5 ^8 j) E) @( ~! n           sum=sum+a(i,j)*x(j);
       end
       x(i)=(b(i)-sum)/a(i,i);
. W* t9 D  P  C1 w1 C' }- |   end
. |: l# `9 z# x' E, `
; s4 `# k" F7 U, Y- B0 U) b这个部分实现了回代过程，得到方程组的解向量 x。

4.打印结果：
' O, h$ ?) g5 L/ n: I- R& Q' q
% |7 U# [! L- r% A, o4 ~8 n   jie=x'
& r; f: b7 ]4 m  v' O; i3 G! \
最后，打印求解得到的解向量。
在实际应用中，可以使用 MATLAB 提供的 linsolve 函数来更稳妥地解线性方程组。
0 l# J: F' ]3 g0 s" t( l+ c

# N8 X. b* C* D$ Y* H3 j