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

　　

【正文】 循環(huán)，那么其執(zhí)行速度將會(huì)接近于純粹用 C語言來寫的程序代碼。 MATLAB是一個(gè)完整的科學(xué)計(jì)算與可視化環(huán)境，但在多數(shù)情況下，如 MATLAB編程的代碼執(zhí)行效率不佳、你希望隱藏你的程序代碼以保護(hù)產(chǎn)權(quán)、或您想要進(jìn)行任何只有在某些高級(jí)語言程序才能做到的事，這時(shí)就需要用 MATLAB 所提供的應(yīng)用程序接口。 東 華理工大學(xué)畢業(yè)設(shè)計(jì) 致謝 30 致 謝 本論文是在導(dǎo)師吳良才教授的熱忱關(guān)懷和細(xì)心指導(dǎo)下完成的。從論文選題到最后的成文，無不凝結(jié)著導(dǎo)師的心血。導(dǎo)師開拓創(chuàng)新的精神、嚴(yán)謹(jǐn)踏實(shí)的治學(xué)態(tài)度以及寬厚的為人都使我深受熏陶，并將使我終身受益。在此謹(jǐn)向?qū)熤乱灾孕牡母兄x和誠摯的敬意。此外，也感謝測(cè)繪專業(yè)的全體老師，正因?yàn)橛辛怂麄兊奶撔慕虒?dǎo)，使我的知識(shí)得到一點(diǎn)一滴的積累，才從容的完成今天的畢業(yè)設(shè)計(jì)。 東華理工大學(xué)畢業(yè)設(shè)計(jì) 參考文獻(xiàn) 31 參考文獻(xiàn) [1] 崔明理 . 控制測(cè)量手冊(cè) . 山西 : 山西科學(xué)技術(shù)出版社 , . [2] 吳俊昶 . 控制網(wǎng)測(cè)量平差 2 版 . 北京 : 北京測(cè)繪出版社 , . [3] 譚 輝 . 測(cè)量程序與新型全站儀的應(yīng)用 . 北京 : 機(jī)械工業(yè)出版社 , 2021. [4] 孔祥元 , 梅是義 . 控制測(cè)量學(xué)下冊(cè) 2 版 . 武漢 : 武漢大學(xué)出版社 , . [5] 武漢大學(xué)測(cè)繪學(xué)院測(cè)量平差學(xué)科組 . 誤差理論與測(cè)量平差基礎(chǔ) . 武漢 : 武漢大學(xué)出版社 , . [6] 張智星 . MATLAB 程序設(shè) 計(jì)與應(yīng)用 . 北京 : 清華大學(xué)出版社 , . [7] 郭九訓(xùn) . 控制網(wǎng)平差程序設(shè)計(jì) . 北京 : 原子能出版社 , . [8] 姚連壁 , 周小平 . 基于 MATLAB 的控制網(wǎng)平差程序設(shè)計(jì) . 上海 : 同濟(jì)大學(xué)出版社 , . [9] 清源計(jì)算機(jī)工作室 . MATLAB 基礎(chǔ)及其應(yīng)用 [M]. 北京 : 機(jī)械工業(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) 東華理工大學(xué)畢業(yè)設(shè)計(jì) 附錄一 32 附錄一 水準(zhǔn)網(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)。 %輸出水準(zhǔn)網(wǎng)計(jì)算結(jié)果 return function[dh,h,V,L,uw0,uwh,uw1]=calculatelevel(ed,dd,sd,pn,h0,be,en,hd,distance)。 % 水準(zhǔn)平差網(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ù)項(xiàng) 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)。 %待定點(diǎn)高程近似值 h=h0+dh。 %待定點(diǎn)高程平差值 V=A*dh1。 %高差觀測(cè)值改正數(shù) L=hd+V。 %高差觀測(cè)值平差值 %精度評(píng)定 uw0=sqrt(V39。*p*V/(gddd))。 %單位權(quán)中誤差 Qxx=inv(A39。*p*A)。 uwh=uw0*sqrt(diag(Qxx))。 %待定點(diǎn)高程平差值中誤差 uwh(dd+1:ed+dd)=。 東華理工大學(xué)畢業(yè)設(shè)計(jì) 附錄一 33 Qff=A*Qxx*A39。 uwh=uw0*sqrt(diag(Qff))。 %高差平差值中誤差 h=[h。h00]。 %所有點(diǎn)高程 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)。 %已知點(diǎn)個(gè)數(shù) dd=fscanf(fid1,39。%f39。,1)。 %未知點(diǎn)個(gè)數(shù) sd=ed+dd。 %總點(diǎn)數(shù) gd=fscanf(fid1,39。%f39。,1)。 %觀測(cè)點(diǎn)個(gè)數(shù) pn=fscanf(fid1,39。%f39。,sd)。 %點(diǎn)號(hào) %known data h0=fscanf(fid1,39。%f39。,ed)。 %已知點(diǎn)高程 h0(dd+1:ed+dd)=h0(1:ed) heightdiff=fscanf(fid1,39。%f39。,[4,gd])。 heightdiff=heightdiff39。 k1=heightdiff(:,1)。 %起點(diǎn) k2=heightdiff(:,2)。 %終點(diǎn) 東華理工大學(xué)畢業(yè)設(shè)計(jì) 附錄一 34 k11=heightdiff(:,1)。 %起點(diǎn) k12=heightdiff(:,2)。 %終點(diǎn) h1=heightdiff(:,3)。 %高差 s=heightdiff(:,4)。 %距離 fclose(39。all39。)。 %點(diǎn)號(hào)轉(zhuǎn)換 [k1,k01]=chkdat(sd,pn,k1)。 [k2,k02]=chkdat(sd,pn,k2)。 h0(1:dd)=20210.。 ie=0。 while(1)%計(jì)算近似高程 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 東華理工大學(xué)畢業(yè)設(shè)計(jì) 附錄一 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 東華理工大學(xué)畢業(yè)設(shè)計(jì) 附錄二 36 附錄二 測(cè)邊網(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