萬應(yīng)玲（甘肅建筑職業(yè)技術(shù)學(xué)院，甘肅 蘭州 730050）

雙向消元法在解法方程的應(yīng)用

萬應(yīng)玲（甘肅建筑職業(yè)技術(shù)學(xué)院，甘肅 蘭州 730050）

測量平差中的條件平差法和間接平差法需要解算法方程。法方程是一些多元一次方程構(gòu)成的方程組，教學(xué)過程中學(xué)生手算求解過程一般采用高斯約化法，由于缺乏對結(jié)果的校核而容易出現(xiàn)計(jì)算錯(cuò)誤。本文作者結(jié)合自己的教學(xué)經(jīng)驗(yàn)，通過具體實(shí)例介紹了雙向消元的解析過程，優(yōu)化了其計(jì)算方法并經(jīng)過步步校核消元確保了結(jié)果的正確性。

高斯約化法；雙向消元法；檢核計(jì)算

在測量平差中, 各種觀測量與已知條件以及各種圖形幾何條件之間可組成條件方程組或觀測量的誤差方程組。進(jìn)行平差解算時(shí),首先對方程組線性化, 形成滿足最小二乘(或滿足最小范數(shù))條件的線性化平差函數(shù)模型[1]。平差問題的求解最后都?xì)w結(jié)為線性方程組的求解過程。因而, 解算線性方程組是解算平差問題的關(guān)鍵[2]。在傳統(tǒng)的解算過程中是將各種平差方程組形成法方程, 通過高斯約化來進(jìn)行求解[3]。高斯約化法是按照方程的順序，從上向下逐個(gè)消元，最后使方程組變?yōu)橐粋€(gè)三角形方程組[4]。高斯約化法對消元過程提出了具體規(guī)定和要求，編制了一定的格式，具有較強(qiáng)的規(guī)律性。但高斯約化法原理的推導(dǎo)符號多，過程繁[5]。有些數(shù)據(jù)的得出需多個(gè)數(shù)據(jù)參與計(jì)算，若有差錯(cuò)，難以快速找到原因[6]。本文采用雙向消元法，最后將每個(gè)方程簡化為一元一次方程并經(jīng)過步步校核消元確保了結(jié)果的準(zhǔn)確性。

1 計(jì)算實(shí)例

設(shè)某法方程中有r個(gè)多元一次方程，采用雙向消元法解算過程如下：1）將第一個(gè)方程的第一個(gè)未知數(shù)系數(shù)變?yōu)?,用其消去其它方程的第一個(gè)未知數(shù)；2）在第一步所得的方程組中，將第二個(gè)方程的第二個(gè)未知數(shù)系數(shù)變?yōu)?（第一個(gè)未知數(shù)已消去），用其向上、向下消去各方程的第二個(gè)未知數(shù)。3）這樣依次處理，到第r步時(shí)，第r個(gè)方程就變?yōu)橐辉淮畏匠?，各方程未知?shù)的序號與該方程序號相同，且系數(shù)為1，這樣各未知數(shù)的解隨之可得。為了進(jìn)一步明確雙向消元法的具體解算過程，可以通過下列表中的數(shù)據(jù)實(shí)例進(jìn)行計(jì)算。

表1列出了某一條件平差計(jì)算過程的法方程系數(shù)。為了敘述簡便，將表1、表2中的某幾個(gè)方程編上帶圓圈的序號。將表1中方程(1)的各系數(shù)（包括∑，∑表示該行左邊各數(shù)之和）除以K1的系數(shù)2.7552，得表2中的方程(1)。在表2方程（1）中，S是由表1中方程（1）的∑除以2.7552而來，其∑仍為左邊各數(shù)之和，當(dāng)∑和S相等時(shí)，說明方程（1）從表1到表2的推算沒有錯(cuò)誤。

表1 法方程系數(shù)表

表2 法方程系數(shù)推算表

將方程①乘以方程②K1系數(shù)的相反數(shù)0.0178并與方程②相加，則消去了方程②的未知數(shù)K1得方程③,該過程是逐項(xiàng)進(jìn)行的，即：

將∑的結(jié)果填到方程③的5欄中，它是推算的方程③中各項(xiàng)系數(shù)應(yīng)有的和。然后再對方程③各項(xiàng)系數(shù)實(shí)際相加，其和也是7.5124。這說明由方程①、②推算方程③是正確的。

從方程①、②、③來看，其中∑是本方程內(nèi)左邊各系數(shù)之和，稱為實(shí)際和。方程①、③是由相關(guān)方程推算來的，其中S是由這些相關(guān)方程中的∑推算的，稱為推算和。若個(gè)別相關(guān)方程中的∑和S不一致時(shí)，仍要用∑推算S,這樣可以減少誤差的積累。當(dāng)∑和S相等或末位數(shù)相差不超過3時(shí)，說明推算過程正確，若∑和S相差較多，說明推算過程有誤，這時(shí)應(yīng)先檢查S是否正確，當(dāng)S無誤時(shí)，再對各系數(shù)的計(jì)算進(jìn)行檢查。

用同樣方法消去表1中方程（3）、（4）的未知數(shù)K1，得表2中的方程（3）、（4）。表2中的方程（1）是用來消去其它方程未知數(shù)K1的，是消元的基礎(chǔ)，把它稱為“基方程”。為了計(jì)算中減少差錯(cuò)，將基方程的系數(shù)下畫一橫線以與其它方程相區(qū)別，而S下不需要畫線。因S是消元中計(jì)算的結(jié)果，沒有參與任何計(jì)算，其作用就是與∑相互校核，檢查消元過程中有無錯(cuò)誤。消去K1后，將表2中方程（2）的K1系數(shù)變?yōu)?，用其向上、向下消去該表中其它方程中的K2，得到表3的結(jié)果。

表3 法方程系數(shù)推算表

表4 法方程系數(shù)推算表

表5 法方程系數(shù)推算表

用同樣的方法消去K3、K4，其結(jié)果如表4和表5所示。為使計(jì)算表格簡潔，可將消去未知數(shù)的欄進(jìn)行合并。如表4的K1～K2欄，表5的K1～K3欄。該欄內(nèi)的數(shù)字1，表示的只是和相應(yīng)方程同序號的未知數(shù)系數(shù)為1，其余未知數(shù)系數(shù)為0；而該欄內(nèi)的0，表示合寫的未知數(shù)都已消去，系數(shù)都為0。

由表5可直接得出未知數(shù)K的值。

2 分析討論

若采用高斯約化法利用消元回帶過程求解方程，如果中途有一個(gè)未知數(shù)求解有誤，則后面所有未知數(shù)的解都是錯(cuò)的，且很難快速找到錯(cuò)誤的原因。高斯約化法充分利用了方程的對稱性, 對于解非線性對稱方程組問題很難實(shí)現(xiàn)[7]，而在雙向消元法中只有向下消元時(shí)才利用對稱性。若有時(shí)計(jì)算的值并不完全對稱，這時(shí)以計(jì)算值為準(zhǔn)，保留其“不對稱性”，以保證該方程的∑和S不會相差太多。

3 結(jié) 語

高斯約化法和雙向消元法的基本思想都通過將一個(gè)方程乘以或除以某個(gè)常數(shù)，以及將兩個(gè)方程相加減這兩種手續(xù)，逐步減少方程中的變元的數(shù)目，最終使每個(gè)方程僅含一個(gè)變元，從而得出所求的解。高斯約化法的求解過程分為消元過程和回代過程兩個(gè)環(huán)節(jié)。這種方法缺點(diǎn)是一旦計(jì)算過程出現(xiàn)錯(cuò)誤從而導(dǎo)致整個(gè)計(jì)算錯(cuò)誤，且很難發(fā)現(xiàn)錯(cuò)誤的原因。而采用雙向消元法可同時(shí)求出法方程中各個(gè)未知數(shù)的解，不需要回代，并且這些一元一次方程都是經(jīng)過一步步校核消元而得出的，因此這些未知數(shù)的解都是正確可靠的。

[1]陶本藻．非線性與線性平差偏差的分布特征[J]．測繪工程,1998(4):7-12．

[2]張東明,呂翠華,李明．點(diǎn)松弛法在測量平差方程組解算中的應(yīng)用[J]．科學(xué)技術(shù)與工程,2010,10(24):5981-5984．

[3]張東明,李明,肖建虹．利用迭代方法求解測量平差中的線性方程組[J]．地礦測繪,2003,19(4):15-17．

[4]趙艷宇 ．Gauss消去法[J]．電大理工,2013(2):51-52．

[5]劉成軍．基于消息傳遞接口的線性方程組并行計(jì)算研究—以改進(jìn)的高斯消元法為例[J]．軟件,2013,34(1):119-120．

[6]顏平．測量平差[M]．北京：中國建筑工業(yè)出版社,2002:42～ 45．

[7]袁功林,李向榮．解非線性對稱方程組問題的具有下降方向的近似高斯-牛頓基礎(chǔ)的BFGS方法[J]．運(yùn)籌學(xué)學(xué)報(bào),2004,8(4):10-26．

Two-way elimination method applied in solution normal equation

Conditional adjustment method and adjustment of observation equations needs solution normal equation．The normal equation is an equation of some multivariate equations．In the process of teaching,the calculation process of the students is generally adopted by Gaussian elimination,Due to the lack of verification of the results and prone to error．In this paper, the author combined with his own teaching experience,Through concrete examples this paper introduces the analytical process of two-way elimination method．Its calculation method was optimized,And after the check elimination step by step to ensure the validity of the results．

Gaussian elimination；two-way elimination method；check calculation

TB302文獻(xiàn)辨識碼：B

1003-8965(2017)01-0050-02