第三章 线性方程组的解算.docVIP

第三章 线性方程组的解算 测绘观测数据处理主要采用最小二乘法进行解算，平差计算中的法方程式为线性方程组。本章主要介绍线性方程组解算及相应的程序设计。 §3-1 高斯（GAUSS）消去法 3.1.1 非对称方程组求解 设有非对称方程Ax = b，其中，A是 m×m阶非奇异实矩阵，b是给定的m×1阶向量。为了说明简便，下面表3-1以4阶方程组为例介绍消去法的具体做法。（设矩阵A的主元素不等于零） 表3-1 消去法解非对称线性方程 a11 a12 a13 a14 a21 a22 a23 a24 a31 a32 a33 a34 a41 a42 a43 a44 b1 b2 b3 b4 消 去 第 一 列 a11 a12 a13 a14 0 a′22 a′23 a′24 0 a′32 a′33 a′34 0 a′42 a′43 a′44 b1 b′2 b′3 b′4 a′ij = aij－ai1 a1j / a11 b′i = bi －ai1 b1 / a11 i , j = 2,3,4 消 去 第 二 列 a11 a12 a13 a14 0 a′22 a′23 a′24 0 0 a33 a34 0 0 a43 a44 b1 b′2 b3 b4 aij=a′ij－a′i2a′2j / a′22 bi=b′i－a′i2 b′2 / a′22 i , j = 3,4 ……………………… … …………………… 根据上表的消去过程，归纳如下： 消去第k个未知数是从k+1行开始，使第k列所有元素为零。（k = 1，2，…， m－1） 从第k+1列开始，各元素为 (3-1) (3-2) 用赋值语句写为 (3-3) (3-4) k = 1，2，…， m－1；i = k+1， … ，m； j = k+1， … ，m。 将上述消去过程用VB语言写成程序如下： Sub GS1(m as Integer) Dim c as Double For k=1 To m-1 For i=k+1 To m c=a(I,k)/a(k,k) For j=k+1 to m a(I,J)=a(I,J)-c *a(k,J) Next j b(i)=b(i)-b(k)*c Next i Next k End Sub 3.1.2 列主元消去法 为了避免主元素等于零而使消去过程失败，有时由于主元素的有效数字严重损失而使舍入误差的相对误差增大，舍入误差扩散使得最终答解极不准确。为克服这些困难，选主元消去法是有效的。选主元消去法有列主元消去法、行主元消去法、全选主元消去法几种。列主元消去法是从第k列（k = 1,2, …, m）的各元素寻出最大元作为主元，其具体做法如下： 从第1列开始，寻找本列绝对值最大的元素，记下该元素所在的行数（i0），并把该行（i0）与第1行（k行）的所有元素互换（即把第i0行移至第k行，使主元素是本列中最大元）。 按式（3-1）和（3-2）进行消元。 将上述消去过程用VB语言写成程序如下： Sub GS2(m As Integer) Dim t As Double Dim i0 As Integer Dim c As Double For k= 1 To m - 1 c = 0# i0 = k For i = k To m If Abs(A(i, k)) Abs(c) Then c = a(i, k) i0 = i End If Next i If i0 =k Then GoTo 30 For i =k To m t = A(k, j) a(k, j) = a(i0, j) a(i0,j) = t Next j t = b(k) b(k) = b(i0) b(i0) = t 30: For i = k + 1 To m c = a(i, k) / a(k, k) For j= k + 1 To m a(i, j) = a(i, j) - c * a(k, j) Next j b(i) =