【正文】 endfunction FileMenu_Callback(hObject, eventdata, handles)function OpenMenuItem_Callback(hObject, eventdata, handles)file = 。 case 4 plot(membrane)。 case 2 plot(sin(1::))。)。popup_sel_index = get(, 39。function pushbutton1_Callback(hObject, eventdata, handles)axes()。) plot(rand(5))。),39。if strcmp(get(hObject,39。endfunction xtjm_OpeningF(hObject, eventdata, handles, varargin) = hObject。endif nargout [varargout{1:nargout}] = gui_mainf(gui_State, varargin{:})。amp。, [])。, [] , ... 39。, xtjm_OutputF, ... 39。, xtjm_OpeningF, ... 39。, gui_Singleton, ... 39。, mfilename, ... 39。gui_State = struct(39。另外，我還要感謝郭豐倫老師給我提出的寶貴意見，對(duì)我后來系統(tǒng)界面的修改有很大的幫助。她有著一顆對(duì)學(xué)生積極負(fù)責(zé)的心，經(jīng)常主動(dòng)詢問我畢業(yè)設(shè)計(jì)的進(jìn)展以及遇到的問題，促使我更加認(rèn)真扎實(shí)地投入到畢業(yè)設(shè)計(jì)中。致 謝致 謝本畢業(yè)設(shè)計(jì)和論文的書寫是在張海玲老師的精心指導(dǎo)下完成的。但是缺乏更多的數(shù)據(jù)、更長(zhǎng)期的檢驗(yàn)來驗(yàn)證其穩(wěn)定性和精度，是該系統(tǒng)的不足之處。（即北京54坐標(biāo)、國(guó)家80坐標(biāo)、WGS84坐標(biāo)、2000國(guó)家坐標(biāo)以及地方獨(dú)立坐標(biāo)）⑶ 該系統(tǒng)還能將轉(zhuǎn)換后的坐標(biāo)和參數(shù)存儲(chǔ)到程序路徑下的word文件中。該系統(tǒng)滿足了目前國(guó)內(nèi)坐標(biāo)轉(zhuǎn)換的需求，利用參數(shù)模型（即三參數(shù)模型、四參數(shù)模型、七參數(shù)模型）和MATLAB平臺(tái)實(shí)現(xiàn)了以下功能：⑴ 在選擇參考橢球體以后，系統(tǒng)輸出長(zhǎng)半軸和扁率，能夠?qū)崿F(xiàn)相同坐標(biāo)系下不同坐標(biāo)之間的轉(zhuǎn)換。本文詳細(xì)介紹了大地測(cè)量系統(tǒng)，包括地球表面、大地水準(zhǔn)面、基準(zhǔn)、三種坐標(biāo)形式、常用坐標(biāo)系以及坐標(biāo)轉(zhuǎn)換模型等，論述坐標(biāo)轉(zhuǎn)換理論為系統(tǒng)的開發(fā)奠定基礎(chǔ)。 3 第六章 總結(jié)與展望本設(shè)計(jì)是基于MATLAB平臺(tái)的坐標(biāo)轉(zhuǎn)換系統(tǒng)的開發(fā)，實(shí)現(xiàn)了同一坐標(biāo)系下不同坐標(biāo)之間的轉(zhuǎn)換和不同坐標(biāo)系之間的相互轉(zhuǎn)換，能夠滿足不同工程測(cè)量對(duì)坐標(biāo)的需求?，F(xiàn)將轉(zhuǎn)換成果展示如下： 大地坐標(biāo)轉(zhuǎn)換為空間直角坐標(biāo)輸出的坐標(biāo)文件內(nèi)容如下：長(zhǎng)半軸= 短半軸= 扁 率= 第一偏心率= 轉(zhuǎn)換后的大地坐標(biāo)（空間直角坐標(biāo)） 點(diǎn)名 X(m) Y(m) Z(m) 1 2 3 4 5 不同坐標(biāo)系之間的相互轉(zhuǎn)換，選擇不同的模型參數(shù)不同，轉(zhuǎn)換精度不一樣。)。,39。),39。,39。else set(findobj(39。)。,39。fprintf(myy,39。w39。39。)。,39。),39。,39。)==1 set(findobj(39。),39。,39。) returns toggle state of radiobutton2% Executes on button press in radiobutton3.if get(findobj(39。 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。K39。ForegroundColor39。radiobutton239。tag39。fclose(myy)。139。%s39。)。,39。 myy=fopen(39。R39。ForegroundColor39。radiobutton239。tag39。Value39。radiobutton239。tag39。Value39。通常只能選擇其中一個(gè)，當(dāng)被選中時(shí)圓圈中心會(huì)變成黑點(diǎn)，其余未被選中的選項(xiàng)成了空白。在相應(yīng)控件下輸入要實(shí)現(xiàn)功能的語(yǔ)句或要調(diào)用的函數(shù)，即可完成該系統(tǒng)最終的設(shè)計(jì)目的。qtdf(:,4)=qtdf(:,3)+gcyc。for i=1:mm gcyc=gcyc+pt(i,3)ps(i,3)。cs(2)=css(2)。zxs(2)]。zxt(2)][cs(3) cs(4)。 %平移因子（x0,y0）rf0=atan(cs(4)/cs(3))。 %尺度因子py=[cs(1)。 cs(2)]。 %參數(shù)參數(shù)改正cs=cscz+csgz。*bb)*bb39。li]。bbi]。pt(i,2)]bbi*cscz。0 1 ps(i,2) ps(i,1)]。ll=[]。y2]。y1。1 0 xx2 yy2。cscz=inv([1 0 xx1 yy1。xx2=ps(2,1)。xx1=ps(1,1)。x2=pt(2,1)。endx1=pt(1,1)。 ps(i,1)=ps(i,1)zxs(1)。 %中心化for i=1:mmpt(i,1)=pt(i,1)zxt(1)。zxt=sum(pt)/mm。pt0=pt。根據(jù)以上平差過程，即可完成程序的編寫。此處以四參數(shù)模型為例進(jìn)行說明。z]。y]。x]。 z=(N(i)*(1e^2)+H(i))*sin(B(i))。for i=1:h x=(N(i)+H(i))*cos(B(i))*cos(L(i))。Y=[]。n]。for i=1:h n=cz/sqrt((1e^2*(sin(B(i)))^2))。 [h l]=size(pt)。大地坐標(biāo)轉(zhuǎn)換為空間直角坐標(biāo)的部分程序如下（完整程序見附錄）：dz=cz*(1bl)。弧度轉(zhuǎn)度分秒的程序見附錄。 jrad(myi,myj)=(degd+degm/+degs)*pi/。 degm=fix(degf)。 degd=fix(ddeg)。180程序如下：%function jrad=deg2rad(deg)[my1,my2]=size(deg)。 %關(guān)閉文件 角度和弧度的轉(zhuǎn)換大地坐標(biāo)一般是用角度值表示的，而在轉(zhuǎn)換計(jì)算過程中所有的值都是弧度值，所以需要先把角度和弧度進(jìn)行轉(zhuǎn)化。survey(mm).x survey(mm).y survey(mm).h]。,1)。 survey(mm).h=fscanf(fm,39。%f39。,1)。%survey(mm)為為結(jié)構(gòu)體，用fscanf方式讀取一組字符型數(shù)據(jù)存入survey(mm).qdm survey(mm).x=fscanf(fm,39。%s39。 %定義變量while (~feof(fm)) %循環(huán)體while mm=mm+1。 %打開文件mm=0。r39。 %將路徑和文件名賦給變量fnalfpath=P。請(qǐng)選擇公共點(diǎn)舊坐標(biāo)系的坐標(biāo)文件：39。*.dat39。文件讀取程序：[F,P]=uigetfile(39。 坐標(biāo)轉(zhuǎn)換系統(tǒng)程序 文件讀取進(jìn)行坐標(biāo)轉(zhuǎn)換時(shí)首先要進(jìn)行坐標(biāo)文件的讀取，包括公共文件的讀取和非公共點(diǎn)（即要轉(zhuǎn)換的坐標(biāo)文件）的讀取。程序內(nèi)容包括坐標(biāo)文件的讀取、三種參數(shù)模型的運(yùn)算、數(shù)據(jù)結(jié)果的輸出、角度和弧度的相互轉(zhuǎn)換以及系統(tǒng)界面的總程序。（以下示圖以WGS84坐標(biāo)系轉(zhuǎn)西安80坐標(biāo)系為例）圖47 選擇坐標(biāo)系和轉(zhuǎn)換模型圖 48 選擇公共點(diǎn)WGS84坐標(biāo)系文件圖 49 選擇公共點(diǎn)西安80坐標(biāo)系文件圖 410 選擇要轉(zhuǎn)換的非公共點(diǎn)WGS84坐標(biāo)系文件以上即完成了坐標(biāo)轉(zhuǎn)換系統(tǒng)的系統(tǒng)界面設(shè)計(jì)，通過點(diǎn)擊控件即可完成坐標(biāo)轉(zhuǎn)換，坐標(biāo)轉(zhuǎn)換系統(tǒng)所調(diào)用的程序詳細(xì)見第五章。同一坐標(biāo)系下不同坐標(biāo)之間進(jìn)行轉(zhuǎn)換：①首先選擇自左向右轉(zhuǎn)換的兩種坐標(biāo)；②選擇參考橢球體，這時(shí)會(huì)返回橢球體的長(zhǎng)半軸和扁率；③點(diǎn)擊“轉(zhuǎn)換”按鈕，按提示選擇文件即可完成轉(zhuǎn)換，轉(zhuǎn)換后的坐標(biāo)在程序路徑下的word文件里。如圖展示界面設(shè)計(jì)過程：圖43 新建GUIDE 在【Create New GUI】的頁(yè)面中可以選擇四種模版： Blank GUI (Default) ：空白樣版 GUI with Uicontrols ：含有UI控件物體的樣版 GUI with Axes and Menu ：含有圖軸和下拉選單的樣版 Modal Question Dialog ：含有對(duì)話框的樣版 在本文中坐標(biāo)轉(zhuǎn)換系統(tǒng)界面設(shè)計(jì)時(shí)采用的是第二種模版，即含有UI控件物體的樣版。用戶通過一定的方法（如鼠標(biāo)或鍵盤）選擇、激活這些圖形對(duì)象，使計(jì)算機(jī)產(chǎn)生某種動(dòng)作或變化，比如實(shí)現(xiàn)計(jì)算、繪圖等。 總體功能設(shè)計(jì)不同坐標(biāo)系之間的轉(zhuǎn)換 國(guó)家80坐標(biāo)系 北京54坐標(biāo)系 地方獨(dú)立坐標(biāo)系 2000國(guó)家坐標(biāo)系 WGS84坐標(biāo)系 三參數(shù)模型、四參數(shù)模型、 七參數(shù)模型 相同坐標(biāo)系下不同坐標(biāo)之間的轉(zhuǎn)換 空間直角坐標(biāo) （X,Y,Z）大地坐標(biāo)（L,B,H） 轉(zhuǎn)換 系統(tǒng)設(shè)計(jì)算法框圖 開始 坐標(biāo)轉(zhuǎn)換總程序流程圖 讀取已知數(shù)據(jù) 進(jìn)行參數(shù)計(jì)算 選擇不同參考橢球體 結(jié)束轉(zhuǎn)換結(jié)果輸出 進(jìn)行坐標(biāo)轉(zhuǎn)換 參數(shù)設(shè)置轉(zhuǎn)換 參數(shù)計(jì)算流程圖 開始 讀取公共點(diǎn)坐標(biāo) 組成系數(shù)矩陣B 組成 組成閉合差L 組成常數(shù)項(xiàng) 結(jié)束 計(jì)算 系統(tǒng)界面設(shè)計(jì) Matlab中系統(tǒng)界面設(shè)計(jì)是指：人和機(jī)器（程序）之間交互作用的工具和方法。圖 42 Matlab圖形用戶界面（GUIDE） 系統(tǒng)結(jié)構(gòu)與功能設(shè)計(jì) 根據(jù)工程測(cè)量坐標(biāo)轉(zhuǎn)換的需求，要實(shí)現(xiàn)坐標(biāo)方便快捷地轉(zhuǎn)換，利用Matlab設(shè)計(jì)坐標(biāo)轉(zhuǎn)換系統(tǒng)。通過GUIDE設(shè)計(jì)用戶界面，通過提供的各種控件我們能夠很快的設(shè)計(jì)我們想要的界面。Matlab作為坐標(biāo)轉(zhuǎn)換的平臺(tái)，能夠通過編寫程序讀取大批量的坐標(biāo)數(shù)據(jù)，根據(jù)轉(zhuǎn)換模型結(jié)合最小二乘原理，利用其矩陣運(yùn)算功能便能計(jì)算出模型參數(shù)，最終把參數(shù)代入要轉(zhuǎn)換的坐標(biāo)公式中就能實(shí)現(xiàn)坐標(biāo)轉(zhuǎn)換。無(wú)論是不同坐標(biāo)系之間的相互轉(zhuǎn)換，還是相同坐標(biāo)系下不同坐標(biāo)之間的轉(zhuǎn)換，都會(huì)涉及矩陣加減、矩陣乘除、矩陣求逆、三角函數(shù)運(yùn)算等等，而且數(shù)值位數(shù)多，單純依靠人工計(jì)算很容易出錯(cuò)，而且當(dāng)數(shù)據(jù)量很大時(shí)非常耗時(shí)。測(cè)繪工程專業(yè)特點(diǎn)之一就是計(jì)算量大，要求計(jì)算能力比較強(qiáng)，數(shù)學(xué)測(cè)繪高級(jí)人才的搖籃，Matlab是多功能的數(shù)學(xué)應(yīng)用軟件，成為測(cè)繪工作者和高校師生重要計(jì)算工具。其目的是方便其所設(shè)立的線性代數(shù)課程的教學(xué)。 大地坐標(biāo)轉(zhuǎn)換成空間直角坐標(biāo) 同一坐標(biāo)系下的大地坐標(biāo)(L,B,H)轉(zhuǎn)換為空間直角坐標(biāo)(X,Y,Z)的公式: (316) (317) (318)上式