【正文】 換的4個參數(shù)：,，其中，；,P取單位矩陣。根據(jù)以上平差過程，即可完成程序的編寫。部分程序如下（完整的程序見附錄）：mm=mm1。pt0=pt。ps0=ps。zxt=sum(pt)/mm。zxs=sum(ps)/mm。 %中心化for i=1:mmpt(i,1)=pt(i,1)zxt(1)。pt(i,2)=pt(i,2)zxt(2)。 ps(i,1)=ps(i,1)zxs(1)。ps(i,2)=ps(i,2)zxs(2)。endx1=pt(1,1)。y1=pt(1,2)。x2=pt(2,1)。y2=pt(2,2)。xx1=ps(1,1)。yy1=ps(1,2)。xx2=ps(2,1)。yy2=ps(2,2)。cscz=inv([1 0 xx1 yy1。0 1 yy1 xx1。1 0 xx2 yy2。0 1 yy2 xx2])*[x1。y1。x2。y2]。bb=[]。ll=[]。mm=mm1for i=1:mm bbi=[1 0 ps(i,1) ps(i,2)。0 1 ps(i,2) ps(i,1)]。 li=[pt(i,1)。pt(i,2)]bbi*cscz。 bb=[bb。bbi]。ll=[ll。li]。endcsgz=inv(bb39。*bb)*bb39。*ll。 %參數(shù)參數(shù)改正cs=cscz+csgz。 %參數(shù)py=[cs(1)。 cs(2)]。lmd=sqrt(cs(3)*cs(3)+cs(4)*cs(4))。 %尺度因子py=[cs(1)。 cs(2)]/lmd。 %平移因子（x0,y0）rf0=atan(cs(4)/cs(3))。 %轉(zhuǎn)角因子css=py+[zxt(1)。zxt(2)][cs(3) cs(4)。cs(4) cs(3)]*[zxs(1)。zxs(2)]。cs(1)=css(1)。cs(2)=css(2)。%計算平均高程異常gcyc=0。for i=1:mm gcyc=gcyc+pt(i,3)ps(i,3)。Endgcyc=gcyc/mm。qtdf(:,4)=qtdf(:,3)+gcyc。 系統(tǒng)控件程序在利用guide命令完成了系統(tǒng)界面設(shè)計后，會自動生成包含已設(shè)計控件的程序。在相應控件下輸入要實現(xiàn)功能的語句或要調(diào)用的函數(shù)，即可完成該系統(tǒng)最終的設(shè)計目的。以radio button控件進行說明，該控件是單選按鈕，是一組帶有文字提示的選擇項。通常只能選擇其中一個，當被選中時圓圈中心會變成黑點，其余未被選中的選項成了空白。程序如下：function radiobutton2_Callback(hObject, eventdata, handles)% hObject handle to radiobutton2 (see GCBO)% eventdata reserved to be defined in a future version of MATLAB% handles structure with handles and user data (see GUIDATA)% Hint: get(hObject,39。Value39。) returns toggle state of radiobutton2% Executes on button press in radiobutton3.if get(findobj(39。tag39。,39。radiobutton239。),39。Value39。)==1 set(findobj(39。tag39。,39。radiobutton239。),39。ForegroundColor39。,39。R39。)。 myy=fopen(39。39。,39。w39。)。fprintf(myy,39。%s39。,39。139。)。fclose(myy)。else set(findobj(39。tag39。,39。radiobutton239。),39。ForegroundColor39。,39。K39。)。 end 2_Callback(hObject, eventdata, handles)% hObject handle to radiobutton2 (see GCBO)% eventdata reserved to be defined in a future version of MATLAB% handles structure with handles and user data (see GUIDATA)% Hint: get(hObject,39。Value39。) returns toggle state of radiobutton2% Executes on button press in radiobutton3.if get(findobj(39。tag39。,39。radiobutton239。),39。Value39。)==1 set(findobj(39。tag39。,39。radiobutton239。),39。ForegroundColor39。,39。R39。)。 myy=fopen(39。39。,39。w39。)。fprintf(myy,39。%s39。,39。139。)。fclose(myy)。else set(findobj(39。tag39。,39。radiobutton239。),39。ForegroundColor39。,39。K39。)。 end 實例檢測通過本坐標轉(zhuǎn)換系統(tǒng)可完成相同坐標系下不同坐標的轉(zhuǎn)換和不同坐標系之間的相互轉(zhuǎn)換?，F(xiàn)將轉(zhuǎn)換成果展示如下： 大地坐標轉(zhuǎn)換為空間直角坐標輸出的坐標文件內(nèi)容如下：長半軸= 短半軸= 扁 率= 第一偏心率= 轉(zhuǎn)換后的大地坐標（空間直角坐標） 點名 X(m) Y(m) Z(m) 1 2 3 4 5 不同坐標系之間的相互轉(zhuǎn)換，選擇不同的模型參數(shù)不同，轉(zhuǎn)換精度不一樣。七參數(shù)模型輸出的坐標文件內(nèi)容如下：七個轉(zhuǎn)換參數(shù)三個平移參數(shù),一個尺度因子,三個旋轉(zhuǎn)參數(shù)為：三個平移參數(shù)：dx= dy= dz= 尺度因子:m=三個旋轉(zhuǎn)參數(shù): a1= a2= a2= 點名 X(m) Y(m) Z(m) 1 2 3 4 5 因為篇幅有限，其它模型的輸出數(shù)據(jù)內(nèi)容沒有展示，可查看附錄程序了解輸出的內(nèi)容。 3 第六章 總結(jié)與展望本設(shè)計是基于MATLAB平臺的坐標轉(zhuǎn)換系統(tǒng)的開發(fā)，實現(xiàn)了同一坐標系下不同坐標之間的轉(zhuǎn)換和不同坐標系之間的相互轉(zhuǎn)換，能夠滿足不同工程測量對坐標的需求。本文中就坐標轉(zhuǎn)換的理論和系統(tǒng)平臺的操作以及程序的編寫進行了論述。本文詳細介紹了大地測量系統(tǒng)，包括地球表面、大地水準面、基準、三種坐標形式、常用坐標系以及坐標轉(zhuǎn)換模型等，論述坐標轉(zhuǎn)換理論為系統(tǒng)的開發(fā)奠定基礎(chǔ)?；贛ATLAB平臺的該系統(tǒng)，涉及了MATLAB軟件中的諸多功能，如文件數(shù)據(jù)的批量讀取、數(shù)據(jù)的批量計算（即強大的矩陣運算能力）、函數(shù)的應用、用戶界面的設(shè)計（即GUIDE命令）等。該系統(tǒng)滿足了目前國內(nèi)坐標轉(zhuǎn)換的需求，利用參數(shù)模型（即三參數(shù)模型、四參數(shù)模型、七參數(shù)模型）和MATLAB平臺實現(xiàn)了以下功能：⑴ 在選擇參考橢球體以后，系統(tǒng)輸出長半軸和扁率，能夠?qū)崿F(xiàn)相同坐標系下不同坐標之間的轉(zhuǎn)換。⑵ 選擇轉(zhuǎn)換參數(shù)模型后，通過讀取公共點坐標可以實現(xiàn)不同坐標系之間的相互轉(zhuǎn)換。（即北京54坐標、國家80坐標、WGS84坐標、2000國家坐標以及地方獨立坐標）⑶ 該系統(tǒng)還能將轉(zhuǎn)換后的坐標和參數(shù)存儲到程序路徑下的word文件中。本系統(tǒng)完成了初步檢測，功能都能得到實現(xiàn)，運算數(shù)據(jù)沒有錯誤。但是缺乏更多的數(shù)據(jù)、更長期的檢驗來驗證其穩(wěn)定性和精度，是該系統(tǒng)的不足之處。通過本系統(tǒng)可以真正實現(xiàn)坐標轉(zhuǎn)換的方便、快捷與簡潔，同時在以后測繪中新出現(xiàn)的轉(zhuǎn)換模型和新的坐標系都可以添加到該系統(tǒng)中去，使得該系統(tǒng)具有可升級性。致 謝致 謝本畢業(yè)設(shè)計和論文的書寫是在張海玲老師的精心指導下完成的。張老師有著嚴謹務實的治學態(tài)度，為我畢業(yè)設(shè)計提出許多寶貴的意見和建議，并指導我不斷的修改。她有著一顆對學生積極負責的心，經(jīng)常主動詢問我畢業(yè)設(shè)計的進展以及遇到的問題，促使我更加認真扎實地投入到畢業(yè)設(shè)計中。同時在此次畢業(yè)設(shè)計中，我對坐標轉(zhuǎn)換理論和MATLAB有了更加深入的認識，又學到了不少東西。另外，我還要感謝郭豐倫老師給我提出的寶貴意見，對我后來系統(tǒng)界面的修改有很大的幫助。最后我要感謝我的父母，謝謝他們這些年的支持讓我順利完成學業(yè)！ 35 參考文獻[1寧津生. .，~11[2]武漢大學測繪學院測量平差學科組,誤差理論與測量平差基礎(chǔ)[M]武漢,武漢大學出版社,2009.[3]金書林,馬金鈴,[M],北京,煤炭工業(yè)出版社,1989.[4]王正林,[M].北京,電子工業(yè)出版社,2006.[5][M],徐州,中國礦業(yè)大學出版社,2005.[6]姚吉利,[M]山東理工大學,2010.[7]張勤,[M]北京,科學出版社,2005.[8]張志永,[M]北京,北京航空航天大學出版社,2006.[9]李宜龍,徐衛(wèi)明. 基于轉(zhuǎn)換異常場的平差四參數(shù)法 GPS 坐標轉(zhuǎn)換技術(shù),~32[10]劉大杰,陶本藻. 實用測量數(shù)據(jù)處理方法[M].北京：測繪出版社,2000[11]李宜龍,徐衛(wèi)明. 基于轉(zhuǎn)換異常場的平差四參數(shù)法 GPS 坐標轉(zhuǎn)換技術(shù),~32[12]孫健. 北京 54 坐標系中的坐標計算. :77~82[13]孔祥元,郭際明,[M]武漢,武漢大學出版社,2010. 附 錄附 錄系統(tǒng)總程序（其中調(diào)用的子程序?qū)⒃诤竺嬉来胃缴希篺unction varargout = xtjm(varargin)gui_Singleton = 1。gui_State = struct(39。gui_Name39。, mfilename, ... 39。gui_Singleton39。, gui_Singleton, ... 39。gui_OpeningF39。, @xtjm_OpeningF, ... 39。gui_OutputF39。, @xtjm_OutputF, ... 39。gui_LayoutF39。, [] , ... 39。gui_Callback39。, [])。if nargin amp。amp。 ischar(varargin{1}) = str2func(varargin{1})。endif nargout [varargout{1:nargout}] = gui_mainf(gui_State, varargin{:})。else gui_mainf(gui_State, varargin{:})。endfunction xtjm_OpeningF(hObject, eventdata, handles, varargin) = hObject。guidata(hObject, handles)。if strcmp(get(hObject,39。Visible39。),39。off39。) plot(rand(5))。endfunction varargout = xtjm_OutputF(hObject, eventdata, handles)varargout{1} = 。function pushbutton1_Callback(hObject, eventdata, handles)axes()。cla。popup_sel_index = get(, 39。Value39。)。switch popup_sel_index case 1 plot(rand(5))。 case 2 plot(sin(1::))。 case 3 bar(1:.5:10)。 case 4 plot(membrane)。 case 5 surf(peaks)。endfunction FileMenu_Callback(hObject, eventdata, handles)function OpenMenuItem_Callback(hObject, eventdata, handles)file