高鐵CPⅢ平面程序(matlab版)
1、matlab版CPⅢ平面網(wǎng),用in2文件處理。
2、添加helmert-WF估計(jì)法,模擬CPⅢ數(shù)據(jù)的時(shí)候遇到負(fù)方差問題看文獻(xiàn)了解到的。個(gè)人認(rèn)為要比helmert估計(jì)法原理與實(shí)際效果更好,有興趣可以繼續(xù)研究。
3、推導(dǎo)了helmert估計(jì)法下面的參數(shù)加權(quán)平差,起因是一個(gè)隧道約束點(diǎn)超限問題。因?yàn)榭床欢S加納老師的論文,所以結(jié)合廣義平差推導(dǎo)了helmert估計(jì)法下面的參數(shù)加權(quán)平差,可以幫我看下推導(dǎo)對不對。有興趣可以繼續(xù)研究下,里面有很多有意思的地方。參數(shù)加權(quán)平差,如果方差設(shè)置0.00001,理論上得出的結(jié)果與約束平差要一致。如果不推導(dǎo),直接參數(shù)平差,4公里數(shù)據(jù)可能看不出問題,處理29公里高鐵數(shù)據(jù)會(huì)有大致1.5mm的坐標(biāo)波動(dòng)。同樣如果方差設(shè)置從100~1000的過程屬于趨近自由網(wǎng)平差,但是坐標(biāo)來回?cái)[動(dòng),暫時(shí)沒有找到原因。數(shù)學(xué)好一點(diǎn)的可以用分塊矩陣分析下,會(huì)有一些有意思的結(jié)果。
4、遇到了估計(jì)的負(fù)方差問題,手動(dòng)取巧的方式解決,網(wǎng)上有很多方法解決,暫時(shí)無精力編程,有興趣可以了解下不等式約束平差。
5、優(yōu)化問題。29公里CPⅢ計(jì)算大約3分鐘,很多地方可以優(yōu)化,比如最大的概率坐標(biāo)計(jì)算,想到了一個(gè)方法,沒空去重新梳理,整體算下來估計(jì)可以到1.5分鐘內(nèi)。
6、里面很多地方稍微改動(dòng)代碼(matlab好處在于把精力在想法上)就可以得到其他想要的結(jié)果。
以上僅供參考,不作為其他任何用途。由于工作原因,時(shí)間不夠,如果有哪位達(dá)人有時(shí)間,希望能夠提供一個(gè)界面版本的軟件供大家學(xué)習(xí),避免后面學(xué)習(xí)的人走彎路。
鏈接:https://pan.baidu.com/s/1DpWT4HMKJqX-zB0013PDdQ?pwd=pjtb
提取碼:pjtb