【導讀】計算、科學試驗處理等任務(wù)。以往的測量數(shù)據(jù)處理方法需要編制特定的處理矩陣運算。程序，而且程度復雜，難度大。言相比，具有編程簡單，運算速度快的特點。文中分別闡述了水準網(wǎng)和測邊網(wǎng)程序的。理論基礎(chǔ)、實現(xiàn)步驟和運行結(jié)果。通過實例的分析，總結(jié)出利用MATLAB對測量數(shù)。據(jù)處理有很大的應(yīng)用價值，它縮短了編程的時間，提高工作效率。Keywords：MATLAB；levelingwork；measuringedgework；programming

　　

【正文】 循環(huán)，那么其執(zhí)行速度將會接近于純粹用 C語言來寫的程序代碼。 MATLAB是一個完整的科學計算與可視化環(huán)境，但在多數(shù)情況下，如 MATLAB編程的代碼執(zhí)行效率不佳、你希望隱藏你的程序代碼以保護產(chǎn)權(quán)、或您想要進行任何只有在某些高級語言程序才能做到的事，這時就需要用 MATLAB 所提供的應(yīng)用程序接口。 東 華理工大學畢業(yè)設(shè)計 致謝 30 致 謝 本論文是在導師吳良才教授的熱忱關(guān)懷和細心指導下完成的。從論文選題到最后的成文，無不凝結(jié)著導師的心血。導師開拓創(chuàng)新的精神、嚴謹踏實的治學態(tài)度以及寬厚的為人都使我深受熏陶，并將使我終身受益。在此謹向?qū)熤乱灾孕牡母兄x和誠摯的敬意。此外，也感謝測繪專業(yè)的全體老師，正因為有了他們的虛心教導，使我的知識得到一點一滴的積累，才從容的完成今天的畢業(yè)設(shè)計。 東華理工大學畢業(yè)設(shè)計 參考文獻 31 參考文獻 [1] 崔明理 . 控制測量手冊 . 山西 : 山西科學技術(shù)出版社 , . [2] 吳俊昶 . 控制網(wǎng)測量平差 2 版 . 北京 : 北京測繪出版社 , . [3] 譚 輝 . 測量程序與新型全站儀的應(yīng)用 . 北京 : 機械工業(yè)出版社 , 2021. [4] 孔祥元 , 梅是義 . 控制測量學下冊 2 版 . 武漢 : 武漢大學出版社 , . [5] 武漢大學測繪學院測量平差學科組 . 誤差理論與測量平差基礎(chǔ) . 武漢 : 武漢大學出版社 , . [6] 張智星 . MATLAB 程序設(shè) 計與應(yīng)用 . 北京 : 清華大學出版社 , . [7] 郭九訓 . 控制網(wǎng)平差程序設(shè)計 . 北京 : 原子能出版社 , . [8] 姚連壁 , 周小平 . 基于 MATLAB 的控制網(wǎng)平差程序設(shè)計 . 上海 : 同濟大學出版社 , . [9] 清源計算機工作室 . MATLAB 基礎(chǔ)及其應(yīng)用 [M]. 北京 : 機械工業(yè)出版社 , 2021. [10] 程衛(wèi)國 ,馮 峰 , 姚 東 , 徐 昕 , 等 . MATLAB 應(yīng)用指南 [M ]. 北京 : 人民郵電出版社 , 1999. [11] CHAI Yan2jun , et al . Quasi2accurate Detection Methodsof Gross Errors with MATLAB Language [J ] . Journal of Shandong University of Science and Technology (NatureScience) , 2021 , (3) [12] NIE Gui2gen. Application of MATLAB to Surveying Da2ta Processing [ J ] . Bulletin of Surveying and Mapping ,2021 , (2) 東華理工大學畢業(yè)設(shè)計 附錄一 32 附錄一 水準網(wǎng)平差程序代碼： function level3 [ed,dd,sd,gd,pn,h0,k1,k2,h1,s]=readleveldata。 global pathname _name s_datafile a1_datafile。 global ed dd sd pn gd h0 k1 k2 h1 s dh。 [dh,h,V,L,uw0,uwh,uw1]=calculatelevel(ed,dd,sd,pn,gd,h0,k1,k2,h1,s)。 writeleveldata(pn,k1,k2,h1,V,L,h0,dh,h,uwh,uw0)。 %輸出水準網(wǎng)計算結(jié)果 return function[dh,h,V,L,uw0,uwh,uw1]=calculatelevel(ed,dd,sd,pn,h0,be,en,hd,distance)。 % 水準平差網(wǎng) A=sparse(zeros(sd,gd))。 %求解系數(shù)陣 b=(0:(gd1)*sd)。 A(be39。+b)=1。 A(en39。+b)=1。 A=A39。 A=A(:,1:dd)。 l=zeros(gd,1)。 %求解常數(shù)項 l=h0(be)h0(en)+hd。 p=diag(1./distance)。 %權(quán)陣 dh=inv(A39。*p*A)*A39。*p*1。 %高程改正數(shù) h00=h0(dd+1:sd)。 h0=h0(1:dd)。 %待定點高程近似值 h=h0+dh。 %待定點高程平差值 V=A*dh1。 %高差觀測值改正數(shù) L=hd+V。 %高差觀測值平差值 %精度評定 uw0=sqrt(V39。*p*V/(gddd))。 %單位權(quán)中誤差 Qxx=inv(A39。*p*A)。 uwh=uw0*sqrt(diag(Qxx))。 %待定點高程平差值中誤差 uwh(dd+1:ed+dd)=。 東華理工大學畢業(yè)設(shè)計 附錄一 33 Qff=A*Qxx*A39。 uwh=uw0*sqrt(diag(Qff))。 %高差平差值中誤差 h=[h。h00]。 %所有點高程 h0=[h0。h00]。 dh=[dh。zeros(ed,1)]。 return function [ed,dd,sd,gd,pn,h0,k1,k2,h1,s]=readleveldata。 global pathname _name s_datafile b_datafile。 global ed dd sd pn gd h0 k1 k2 h1 s k11 k12。 k1=[]。k2=[]。h=[]。s=[]。 if(isempty(pathname)|isempty(_name)) [filename,pathname]=uigetfile(39。*.txt39。,39。Input filename39。)。 i=find(39。.39。==filename)。 _name=filename(1:i1)。 end fid1=fopen(strcat(pathname,_name,s_datafile),39。rt39。)。 if(fid1==1) msgbox(39。Input File or Path is not correct39。,39。Warning39。,39。warn39。)。 return。 end %open afile to read %open afile to read ed=fscanf(fid1,39。%f39。,1)。 %已知點個數(shù) dd=fscanf(fid1,39。%f39。,1)。 %未知點個數(shù) sd=ed+dd。 %總點數(shù) gd=fscanf(fid1,39。%f39。,1)。 %觀測點個數(shù) pn=fscanf(fid1,39。%f39。,sd)。 %點號 %known data h0=fscanf(fid1,39。%f39。,ed)。 %已知點高程 h0(dd+1:ed+dd)=h0(1:ed) heightdiff=fscanf(fid1,39。%f39。,[4,gd])。 heightdiff=heightdiff39。 k1=heightdiff(:,1)。 %起點 k2=heightdiff(:,2)。 %終點 東華理工大學畢業(yè)設(shè)計 附錄一 34 k11=heightdiff(:,1)。 %起點 k12=heightdiff(:,2)。 %終點 h1=heightdiff(:,3)。 %高差 s=heightdiff(:,4)。 %距離 fclose(39。all39。)。 %點號轉(zhuǎn)換 [k1,k01]=chkdat(sd,pn,k1)。 [k2,k02]=chkdat(sd,pn,k2)。 h0(1:dd)=20210.。 ie=0。 while(1)%計算近似高程 for k=1:gd i=k1(k)。 j=k2(k)。 if(h0(i)le4amp。h0(j)le4) h0(j)=h0(i)+h1(k)。 ie=ie+1。 end if(h0(i)le4amp。h0(j)le4) h0(i)=h0(j)h1(k)。 ie=ie+1。 end end if(ie==dd) break。 end end h0=reshape(h0,length(h0),1)。 return function[n1,k]=chkdat(sd,pn,n1) n=length(n1)。 k=0。 for i=1:n 東華理工大學畢業(yè)設(shè)計 附錄一 35 i1=0。 for j=1:sd if(n1(i)==pn(j)) i1=1。 n1(i)=j。 break。 end end if(i1==0) % fprintf(fit2,39。%5d%5d\n39。,n1(i))。 k=1。 end end return 東華理工大學畢業(yè)設(shè)計 附錄二 36 附錄二 測邊網(wǎng)平差程序代碼： global ed dd sd dd1 pn x0 y0 m1 m2 m3 ms pp e d sid md g f aa bb cc rt rr tt ql w global pathname _name s_datafile b_datafile filename。 x0=[]。y0=[]。e=[]。d=[]。sid=[]。g=[]。f=[]。pn=[]。 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %%%%%%%數(shù)據(jù)的讀入 if(isempty(pathname)|isempty(_name)) [filename,pathname]=uigetfile(39。*.txt39。,39。Input filename39。)。 i=find(39。.39。==filename)。 _name=filename(1:i1)。 end fit1=fopen(strcat(pathname,filename),39。rt39。)。 if(fit1==1) msgbox(39。Input File or Path is not correct39。,39。Warning39。,39。warn39。)。 return。 end %open afile to read b_datafile=39。39。 fit2=fopen(strcat(pathname,_name,b_datafile),39。wt39。)。 if(fit1==1) msgbox(39。Input File or Path is not correct39。,39。Warning39。,39。warn39。)。 return。 end