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

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

B、P、L、X0否是 否是 R(A)=t<m? 結(jié)束 實驗結(jié)果水準網(wǎng)水準網(wǎng)如圖所示，測得觀測高差為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)如圖所示的測邊自由網(wǎng)（假定各邊等權(quán)），觀測數(shù)據(jù)如表所示，各點近似坐標(biāo)給定如下。邊長觀測值邊號邊長S/m邊號邊長S/m122141.335436577.034227908.063520480.046320044.592629402.438近似坐標(biāo)點名X/mY/mA00B22141.3350C19187.33520265.887D-10068.38617332.434BXB C CAADYD部分程序截屏：程序運行結(jié)果：實驗總結(jié)秩虧自由網(wǎng)平差建立在經(jīng)典平差之上，由于平差中R（N）=t<m，法方程的逆無法求的，采用附加條件未知數(shù)改正和為最小，通過求解法方程的最小范數(shù)逆、偽逆等求解未知數(shù)解；在給定間接平差程序的基礎(chǔ)上添加求解偽逆的部分，只是在間接平差基礎(chǔ)上稍作修改；在程序編譯運行階段由于數(shù)據(jù)類型的錯誤，導(dǎo)致輸出結(jié)果有誤，解決時將每一部分的輸出結(jié)果分別輸出，查找錯誤源。由于沒有搞清楚函數(shù)形參傳遞未知數(shù)個數(shù)還是觀測值個數(shù)，導(dǎo)