秩虧自由網(wǎng)平差程序設(shè)計(jì)上機(jī)實(shí)習(xí)

地質(zhì)工程與測(cè)繪學(xué)院實(shí)習(xí)報(bào)告課程名稱：近代測(cè)量數(shù)據(jù)處理與應(yīng)用實(shí)習(xí)名稱：秩虧自由網(wǎng)平差程序設(shè)計(jì)上機(jī)實(shí)習(xí)班級(jí)：學(xué)號(hào)：姓名：指導(dǎo)教師：實(shí)習(xí)時(shí)間：2013年10月28日實(shí)習(xí)目的理解秩虧自由網(wǎng)平差的函數(shù)模型和隨機(jī)模型；理解廣義逆、最小范數(shù)逆、偽逆的概念學(xué)會(huì)如何求解最小范數(shù)逆和偽逆；理解秩虧自由網(wǎng)平差數(shù)據(jù)的處理的一般過(guò)程。加強(qiáng)程序的理解力、使用給定的程序處理數(shù)據(jù)；實(shí)驗(yàn)基本原理誤差方程由最小二乘原理，可得由于R(N)=R（A）=t<u，故N為奇異矩陣，其正則逆不存在，得不到唯一解。增加約束條件組成新的函數(shù)：求偏導(dǎo)數(shù)并令其等于零，得： 為最小范數(shù)逆 求解方法有：最小范數(shù)逆法 在法方程系數(shù)陣N的對(duì)角線元素加微小量，將秩虧問題轉(zhuǎn)化為經(jīng)典自由網(wǎng)平差采用降秩法求解最小范數(shù)逆?zhèn)文娼夥ú捎脻M秩分解法求解偽逆程序操作流程開始

B、P、L、X0否是 否是 R(A)=t<m? 結(jié)束 實(shí)驗(yàn)結(jié)果水準(zhǔn)網(wǎng)水準(zhǔn)網(wǎng)如圖所示，測(cè)得觀測(cè)高差為h1=12.345mh2=3.478mh3=-15.817m各段路線距離相等（等權(quán)），取H10=0H20=12.345mH30=15.823m1）誤差方程V=B*X-lB=,l=P=,X0=[012.34515.823]T平面網(wǎng)如圖所示的測(cè)邊自由網(wǎng)（假定各邊等權(quán)），觀測(cè)數(shù)據(jù)如表所示，各點(diǎn)近似坐標(biāo)給定如下。邊長(zhǎng)觀測(cè)值邊號(hào)邊長(zhǎng)S/m邊號(hào)邊長(zhǎng)S/m122141.335436577.034227908.063520480.046320044.592629402.438近似坐標(biāo)點(diǎn)名X/mY/mA00B22141.3350C19187.33520265.887D-10068.38617332.434BXB C CAADYD部分程序截屏：程序運(yùn)行結(jié)果：實(shí)驗(yàn)總結(jié)秩虧自由網(wǎng)平差建立在經(jīng)典平差之上，由于平差中R（N）=t<m，法方程的逆無(wú)法求的，采用附加條件未知數(shù)改正和為最小，通過(guò)求解法方程的最小范數(shù)逆、偽逆等求解未知數(shù)解；在給定間接平差程序的基礎(chǔ)上添加求解偽逆的部分，只是在間接平差基礎(chǔ)上稍作修改；在程序編譯運(yùn)行階段由于數(shù)據(jù)類型的錯(cuò)誤，導(dǎo)致輸出結(jié)果有誤，解決時(shí)將每一部分的輸出結(jié)果分別輸出，查找錯(cuò)誤源。由于沒有搞清楚函數(shù)形參傳遞未知數(shù)個(gè)數(shù)還是觀測(cè)值個(gè)數(shù)，導(dǎo)