劉志平,張書畢,卞和方
(中國(guó)礦業(yè)大學(xué) 國(guó)土環(huán)境與災(zāi)害監(jiān)測(cè)國(guó)家測(cè)繪地理信息局重點(diǎn)實(shí)驗(yàn)室,江蘇 徐州221116)
《誤差理論與測(cè)量平差基礎(chǔ)》是測(cè)繪專業(yè)的八大公共專業(yè)基礎(chǔ)課之一,該課程教學(xué)內(nèi)容受到高度重視[1-2]。經(jīng)典平差函數(shù)模型已發(fā)展為條件平差、附有參數(shù)的條件平差、間接平差、附有限制條件的間接平差以及附有限制條件的條件平差模型五種基本形式,但在過(guò)去較長(zhǎng)時(shí)期內(nèi)沿用平差方法介紹孤立的體系[3]。對(duì)此,於宗儔[3]教授提出以附有限制條件的條件平差模型統(tǒng)一各種經(jīng)典平差模型(為敘述方便,簡(jiǎn)稱於氏概括模型),對(duì)引導(dǎo)學(xué)生理解各種平差函數(shù)模型的區(qū)別與聯(lián)系起到了極大的作用。此外,王新洲[4]教授指出概括平差模型只是在形式上統(tǒng)一各種平差模型,但沒有真正揭示各種經(jīng)典平差模型之間的內(nèi)在聯(lián)系,并證明附有限制條件的間接平差即為各種經(jīng)典平差模型的概括模型(簡(jiǎn)稱王氏概括模型)。
在證明各種經(jīng)典平差模型的數(shù)學(xué)等價(jià)性方面,王氏概括模型優(yōu)于於氏概括模型;而在指導(dǎo)測(cè)量平差應(yīng)用中的參數(shù)選取等方面,於氏概括模型優(yōu)于王氏概括模型。例如,文獻(xiàn)[4]在證明條件平差模型是王氏概括模型的特例時(shí),設(shè)參數(shù)個(gè)數(shù)與觀測(cè)數(shù)相等。這本身只是一個(gè)數(shù)學(xué)推導(dǎo)過(guò)程,易使初學(xué)者對(duì)剛掌握的“觀測(cè)數(shù)”、“多余觀測(cè)”、“參數(shù)”等重要概念產(chǎn)生混淆[5]。鑒此,為進(jìn)一步在教學(xué)過(guò)程中突出測(cè)量平差理論的核心、反映經(jīng)典測(cè)量平差的研究對(duì)象-偶然誤差或帶有偶然誤差的觀測(cè)值,并在此基礎(chǔ)上提出另外一種概括模型形式-等價(jià)條件平差模型。
對(duì)于某一平差問(wèn)題,設(shè)觀測(cè)數(shù)為n,必要觀測(cè)數(shù)為t,則多余觀測(cè)數(shù)r=n-t。r個(gè)多余觀測(cè)必然產(chǎn)生r個(gè)關(guān)于觀測(cè)值的一般條件方程。若選用u個(gè)參數(shù),且參數(shù)之間可以形成s個(gè)函數(shù)關(guān)系式,則必然產(chǎn)生u-s個(gè)關(guān)于觀測(cè)值與參數(shù)的一般條件方程、s個(gè)關(guān)于參數(shù)的限制條件方程。因此,一般條件方程個(gè)數(shù)為c=r+u-s,限制條件方程個(gè)數(shù)為s。因此,根據(jù)上述平差問(wèn)題,可形成於氏概括平差模型:
式中:A,B,C為已知的系數(shù)矩陣;W=-(AL+Bx0+A0)為具有參數(shù)的條件方程閉合差;Z=-(Cx0+C0)為限制條件方程閉合差;DL為觀測(cè)值L的方差陣;V,x分別為待求的殘差與參數(shù)改正向量。
為保證平差函數(shù)模型解估計(jì)的存在性,應(yīng)滿足
因此,平差模型中的矩陣應(yīng)滿足:A為行滿秩矩陣,B為列滿秩矩陣或?yàn)榱袧M秩矩陣
C為行滿秩矩陣?;?,結(jié)合三大類參數(shù)選取方式,可以獲得五種基本平差模型形式,見表1。
表1 完備的平差函數(shù)模型
分析三大類參數(shù)選取及其對(duì)應(yīng)概括模型中的矩陣性質(zhì):
1)當(dāng)不選取任何參數(shù)時(shí),其參數(shù)數(shù)目u=0,可建立條件平差模型,即B=0,C=0,A為行滿秩矩陣的概括模型。
2)選取獨(dú)立參數(shù),其數(shù)目理論上滿足0<u≤t。若0<u<t,則可建立附有參數(shù)的條件平差模型,即相當(dāng)于C=0、A為行滿秩矩陣、B為列滿秩矩陣的概括模型。若u=t則可列立間接平差模型,即相當(dāng)于A=-I,C=0。
3)當(dāng)選取的參數(shù)具有函數(shù)相關(guān)的性質(zhì),其數(shù)目理論上滿足u>0。若其中參數(shù)獨(dú)立數(shù)不超過(guò)必要參數(shù)個(gè)數(shù)t,則參數(shù)數(shù)目u>0,可列立附有限制條件的條件平差模型,相當(dāng)于A為行滿秩矩陣,ˉB為列滿秩矩陣,C為行滿秩矩陣的概括模型;若其中參數(shù)獨(dú)立數(shù)恰好為必要參數(shù)個(gè)數(shù)t,則參數(shù)數(shù)目u>t,可列立附限制條件的間接平差模型,即相當(dāng)于A=-I,B為列滿秩矩陣,C為行滿秩矩陣的概括模型。
綜上所述,從參數(shù)數(shù)目u取值范圍的角度看,表1描述的經(jīng)典平差模型是完備的。然而現(xiàn)有文獻(xiàn)[6]介紹的附有限制條件的條件平差,對(duì)參數(shù)數(shù)目的要求(0<u<t)是不完備的,現(xiàn)有教材對(duì)u≥t且其中獨(dú)立參數(shù)少于t的情況沒有考慮。對(duì)此,曾在文獻(xiàn)[7]針對(duì)參數(shù)u=t且其獨(dú)立參數(shù)少于t的附有限制條件的條件平差解法進(jìn)行完善,其思路可以推廣至參數(shù)取值范圍完備的於氏概括平差。
討論參數(shù)完備條件下的於氏概括模型,為平差模型應(yīng)用中的參數(shù)選取提供依據(jù),同時(shí)也建立必要觀測(cè)、多余觀測(cè)、參數(shù)等基本概念與測(cè)量平差模型的聯(lián)系。
1)附有參數(shù)的條件平差模型、間接平差模型。首先,對(duì)任意列滿秩矩陣B,其零空間投影算子[8]HB可以仿照平差因子[9-12]得到
其中,零空間投影算子滿足HBB=0。
進(jìn)而,對(duì)附有參數(shù)的條件平差模型、間接平差模型分別左乘上述零空間投影算子HB,便可消去其中的獨(dú)立參數(shù),由此可得相同形式的條件平差模型
經(jīng)零空間投影算子HB處理后,上述兩種基本平差模型均具有“條件平差模型”的形式。分析零空間投影算子可知,零空間投影算子HB的維數(shù)為c×c,由于附有參數(shù)的條件平差模型與間接平差模型的限制條件方程個(gè)數(shù)s=0,則方程總數(shù)為c=r+u。另外,顧及到零空間投影算子的秩rank(HB)=c-u,表明其僅r個(gè)行向量相互獨(dú)立。因此該“條件平差模型”具有r個(gè)有效的條件方程。
2)附限制條件的條件平差模型、附有限制條件的間接平差模型。與式(2)同理,令ˉBT=[BC]可得列滿秩矩陣B的零空間投影算子 HˉB。
與式(3)、式(4)同理,對(duì)附有限制條件的條件平差模型、附有限制條件的間接平差模型分別左乘上述零空間投影算子HˉB,便可消去其中的相關(guān)參數(shù),得到相同形式的條件平差模型
經(jīng)零空間投影算子HˉB處理之后,上述兩種基本平差模型均具有“條件平差”的形式。分別分析兩類零空間投影算子可知,HˉB的維數(shù)為(c+s )×(c+s ) ,由于附有限制條件的條件平差模型、附有限制條件的間接平差模型具有s個(gè)限制條件,則方程總數(shù)c+s=r+u。另外,考慮零空間投影算子的秩rank (HˉB)=c+s-u,表明其僅r個(gè)行向量相互獨(dú)立。因此該“條件平差模型”具有r個(gè)有效的條件方程。
綜上所述,式(2)、式(5)實(shí)現(xiàn)了四種基本平差模型至條件平差模型的轉(zhuǎn)換,證明了五種基本平差模型的等價(jià)性。因此,可將式(3)、式(4)、式(6)和式(7)以等價(jià)條件平差模型概括為
由前述分析可知,式(3)、式(4)、式(6)和式(7)的方程個(gè)數(shù)不少于r,但僅含有r個(gè)函數(shù)獨(dú)立的條件方程,即等價(jià)條件平差模型。此外,現(xiàn)有的於氏概括模型、王氏概括模型和本文概括形式——等價(jià)條件平差模型,從三個(gè)不同的角度均對(duì)各種經(jīng)典平差模型進(jìn)行正確的概括,這本身也證明了各種經(jīng)典平差模型的等價(jià)性。本文僅對(duì)包含殘差向量V的等價(jià)條件平差模型進(jìn)行概括,該特點(diǎn)便于在教學(xué)過(guò)程中使學(xué)生加深理解平差準(zhǔn)則這一核心內(nèi)容。
實(shí)際上,在教學(xué)中要強(qiáng)調(diào)平差參數(shù)的引入僅是源于測(cè)量工作的實(shí)際需要,所需參數(shù)可能是函數(shù)相關(guān)、也可能是函數(shù)獨(dú)立的,而且參數(shù)個(gè)數(shù)也可因需要確定,但都不影響其平差實(shí)質(zhì)——對(duì)偶然誤差或帶有偶然誤差的觀測(cè)值按最小二乘準(zhǔn)則進(jìn)行調(diào)整。換言之,若不考慮實(shí)際工作,首先,不選任何參數(shù)而建立條件平差模型,之后根據(jù)所需參數(shù)與觀測(cè)值的函數(shù)關(guān)系進(jìn)行參數(shù)估計(jì)、以及利用協(xié)方差傳播律進(jìn)行參數(shù)精度評(píng)定。由此可知,本文概括模型凸顯了條件平差模型及與其他各種經(jīng)典平差模型的關(guān)系,也加強(qiáng)條件獨(dú)立與參數(shù)獨(dú)立、平差準(zhǔn)則與平差模型等基礎(chǔ)概念的聯(lián)系,易使學(xué)生更加深刻理解“同一平差問(wèn)題的平差結(jié)果不因平差函數(shù)模型的不同而不同”。
當(dāng)前,許多學(xué)者對(duì)各種經(jīng)典平差模型的等價(jià)性及其內(nèi)在聯(lián)系作了較為廣泛的探討。本文對(duì)比分析於氏概括模型、王氏概括模型的特點(diǎn),同時(shí)從參數(shù)數(shù)目取值范圍的角度,完善於氏概括模型。在此基礎(chǔ)上,利用零空間投影算子(或稱平差因子)將各種經(jīng)典平差函數(shù)模型轉(zhuǎn)換為等價(jià)條件平差模型,由此得出經(jīng)典平差模型的另一種概括形式。同時(shí)指出,等價(jià)條件平差函數(shù)模型僅包含殘差向量,便于加深學(xué)生對(duì)測(cè)量平差研究對(duì)象(偶然誤差或帶有偶然誤差的觀測(cè)值)、經(jīng)典平差理論核心(最小二乘平差準(zhǔn)則)的理解。此外,推導(dǎo)過(guò)程中零空間投影算子的應(yīng)用,可讓學(xué)生對(duì)矩陣運(yùn)算的重要性留下深刻印象,為學(xué)習(xí)高等平差做鋪墊。
[1] 鄧興升.提高測(cè)量平差課程教學(xué)質(zhì)量的措施[J].測(cè)繪工程,2010,19(3):74-76.
[2] 邱衛(wèi)寧,陶本藻,姚宜斌.誤差理論與測(cè)量平差基礎(chǔ)精品課程的建設(shè)與實(shí)踐[J].測(cè)繪工程,2011,20(1):77-81
[3] 武漢測(cè)繪科技大學(xué)測(cè)量平差教研室.測(cè)量平差基礎(chǔ)[M].3版.北京:測(cè)繪出版社,2007.
[4] 王新洲.論經(jīng)典測(cè)量平差模型的內(nèi)在聯(lián)系[J].測(cè)繪通報(bào),2004,33(2):1-4.
[5] 姚宜斌,邱衛(wèi)寧.測(cè)量平差問(wèn)題中必要觀測(cè)數(shù)的確定[J].2007,36(3):14-15.
[6] 邱衛(wèi)寧,王新洲.測(cè)量平差教學(xué)體系的設(shè)計(jì)與研究[J].測(cè)繪通報(bào),2006,35(2):67-69.
[7] 張書畢,秦永寬,陸翠,等.概括平差模型解法的完善[J].測(cè)繪通報(bào),2007,36(11):7-9.
[8] 劉志平.等價(jià)條件閉合差的方差-協(xié)方差分量估計(jì)解析法[J].測(cè)繪學(xué)報(bào),2013,42(5):648-653.
[9] 陶本澡.平差因子與平差結(jié)構(gòu)[J].大地測(cè)量與地球動(dòng)力學(xué),2002,22(3):6-10.
[10]劉志平,張書畢.方差協(xié)方差分量估計(jì)的概括平差因子法[J].武漢大學(xué)學(xué)報(bào):信息科學(xué)版,2013,38(8):925-929.
[11]張海霞,陳向陽(yáng),趙文普.基于自由網(wǎng)平差的一體化DLG生產(chǎn)試驗(yàn)[J].測(cè)繪工程,2015,24(1):69-72.
[12]黃志偉,劉曉華,龔率,等.基于三維平差的基坑監(jiān)測(cè)數(shù)據(jù)處理方案探討[J].測(cè)繪工程,2015,24(2):51-54.