程序设计实验报告(matlab)

1、登录界面:83登咗界面登陆界面水准平差界面:导入点要素:FileNa*e PathNaiie = uigetfile(r *. txtJ Txt FiLr txt T / +r xlsJ,J Txt Fipu (*. xls)|l，v* idwoisc fileJ):Llengt h (F i leName 1if L<5err or dUC导入文件失政'J文件J ):retijrTi;elseALSgboX C导入文件咸?；文件J )；旳dtest = F£leNanie(lyL-3 江)：switch t estcase " txtJstr «

2、LPathJIajie FileName;fin = f open (st r," F );stx = f£etl(fin).count er=l;vhtle feof(fin) 0舊文件猪针星否到未匡/str = fg«t 1 (fin):la b c d二 srtrjreadtslri %f Sf afr '' dLiait erB jR j" ) zpllabfcouiiter, :) = a b e d；counter = 匚口urrter+I:endcase " xls'str = PathNane Fileu

3、e:fin - misread(str):ptlabfin：endsave pt Tab pt Tab;set (handles, uitble 1 " datJ > ptT：ab);导入边要素:FileManie PathHane = uigetfile (f *. tKt* f J Txt File (*. txt)r / *. xlsJ J Txt File (Kxlsrij* dwose fileJ): L= 1 enft h. (Fi 1-eNa>e):if L.<5errordl( C导入文件共OT /文件J ；return;elsensgboi V导丸

4、文件成功文件巧：endtest - Fild?3M/】,L一引L:switch, testcase r,txt*str - PaihNane FilaNajie;fin = fispen(istr? ? r" ) i5r - fgetlifin);counter=l.讪nW feci'fin) = 0 $交件指针是否到耒產。stt - f grtl(£in) a b c di e = stiread(srtra' if %f Mdelimiter/ /):曲匚 aunt er j.;)占a b c de;counter = counter+l:endends

5、ave ImeTb linelab;set (handlesilit able23 ' dat aa j lin自lab"水准平差代码:iakptTablineTab = gtHQ(ptlabjline!al1: if not (isok)计莫概略离程眞败J;return;*ndP, S3 L, ptlabj lineTab二 <etVptTab(, lineTab) X = inv(B,tP*B)*B,*P*L;V = B*S+L;rovj col二 size (E i :mo = sqrt(矿 *P*V)/(row-col):M = moP 2inv CSJ *P*

6、B);rowj col = size (ptlab);for i=2:rowpt Tab (ia 2) - ptTab(if 2)+1(!-!, 1):enddat1 = ptTab,for i=l:rowlfor j=l:3|j) =ptlab(il j);%datal (ij j) - data ：end-endset (handles, uit able4j J dat aJ j dat a3):lineTabl = line!ab:global -_l'；glob&l d j- 2 ;flobal M;datalcet (handles, uitablelj '

7、 data' ' ; dat a2=get J handles, uitable2j ' dat s' ; dat 2 二 data2 :rawlj coll = slz亡(ti包七比£)"rows, CO12 = size <data2).if datal (row!., I) >0 I I dat alf- rovl-1, I) >0 pt lab = dat al :lineTab = dat a2:fori-lsrowHlsfor j=l:cq!1adata 1 ij j;ptXah (ij j) = dat al

8、i5 j :-endendfori=l: rowZlRfor j=l:co!2lineTab (ij j ) = dat ：a2 ij endtndtndCFileNaiiie FathNanie = uisetfilt (" *B tst" f * Txt File/ =*-k15?/ 1st File 'rKls)pchMe file*):1= length (FileNaune * ;if L<5MrardliCA文件失ST*文件：E EtUEJI ；els*nsebca C导入交萍成功"'t'文件'；test = Fi

9、lcNameChL-iL.switch testcase ?txVstr « PathNane FiltNanefin = fopentstTj'r');str = fget 1 (fin).cnuntei=l;while feof (fin) Q霽文件指甘是否到末昆&sti - fgetKfin;(a b c d = strreadtstrj' M Kf Mf Sf4 s' delimiter* s P / ) zptI謝bi；軸at亡巧 t) = abed；counter = caunter+|;endcase ' . xls-st

10、r = EPathName FiUNamel;fin = misread (str) rpt TabfixL：endsave pt lab ptlab;set handlesi-mtiableli "chth, pt Tab);global 11：global E；r V = sizfldinTabl):for i=1:rline lab I (ij 6) = l/lin-elabl (ij 5);p end11 = linelab1:m = dat a3 :输出平差报告:fprint-f (fid, *«*#*-4*4*#9#*a*«*««*

11、4!*#*«#*«*#Xn ； fpirdiTtf (f idj headl " n");c21 = size(dat 2);i = 1 : r2for j=l:c2£printf (f id, J Sosend.(£xd " n')Lr2'far篝披揩走格式依次写入毎一行勒掘'j xiuji.2st r ( d.t 2 (ij j ):-end fprintf (fid, fprintf (fid, £print f <f id,Lr c lyl £oz s昨軍章*审尊准*

12、來尊囂*冲:冷：*+：#*平碁结果卓事岸木聲准电术草*事字*車車*:戕*卓算；nJ); fJ畔卓鼻字斗卓甲岸卓尊軍字專字鼻車#字字双坝q精度卓字总家專*本朮畫电字曹車席蕈鼻字家耳学.“3 )； head : nJ );sxze(11):策扌畫指定格式依次写入鸟一行数攥for j = 1:c fprintf <f id, " *ios end£pr int f (f ' n'end.f口工f/fzLd, £'*电宙幺陪*量审丘時*水#世沟点位丰青度亍*审：* *半审碎审*审*电fl")；f print-f <f idj

13、head3 ' n/ )：fd u3= si蛊盘(dztl);head 二 °边号起点终点高差弘路线辰headl ='边号起点终点高差/m路线长head?-*点号高程5点标志未知点序号:head3 -尸点号点标志未知点序号点位精度'fid = fopenC 平差抿告 tict', J vt") global global glob al global观测精JT ;a = input dig C请输入现蒯者：平差抿告信息J； b二input dlECW输入双测曰期:/平差抿告信恵T); fprintf (fid, fprmtf (fid fpr

14、intf tf idj fprintf (fidj fprintf (fidj fprintf (f id, fprintf (fid,V观测者:F： tad 'm : £现刻日期匸: m叭；平差报皆nJ):當耳羞凍爲倉案耳倉*耳案當*事索重事覃拳电 点菱袁*柑*“和*1*林*：*心* nJ： Chead2 'if :rl cl = size (bt 1): for i = I i rlfor j=l:cl fprintf (fid, J %6s endfprintf (fid,nJ) hend舅拽指走格式依坎写入每一行数据,nwi2st r (dat (i, j)；

15、fprintf (fid* * 点位精度* 事*累*半*#*、"')； fprintf ffid? heads ' ri* J);r3 c3 = size(dat1): for i = 1 i r3熾指定格式依次写入每一行数据Ifor j=l：c3fprintf (fid,'那 EsnuniiSt r (dat I (i, j);endif i-l>0fprintf(fid,S6selsefprintf (fid, ''0endfprintf (fi 4 ' nJ); nun2str i.v(i-lj j-11):fclose(f

16、id): usgboxC告已输出')登陆界面:登底界面1侧绘2022210301登陆界面严超转换界面:转换代码:鳴和用坐标増里探取方f立窟子程序刪立于第一象限诵 斛立于第二象限话 和立于第三象眼角 紀垃于第四象眼用-1funct i on f> = getfurdet dety)i-f= 0if dety>0Tw - pi/2;retum:elseE-v = 3*pi/2 ;return.:endelsekx 芦 a.bs ( at an ( del： y/ det x ):«zidi£ det x>0 && dety>sO

17、fir n xx :ertdif detx<0 && dety>=0 fw = piM盘：endif detx<0 SlSl dety<=0 fw = pi-M<x :endif detx>0 && dety<=0 fw 三 2*pi-xx: end'： i while fw>2*pi fw = fw-2*pi: endIwhile fw<0fw = fw+2*pj.;endload xx:row; col = size (source) ; Si礬取吕gul亡巴恚中的要秦牛知 jsq = 0;巩+数

18、器養零z for i = 1: rowif source(ij 4) 1 邸艮制条件，从公共点中选取 jsq = jsq+1:舄计数器加！if jsq 乂 I廉计数器为1，说朋第一次巖取到公共点sxx 二 source(ij 21; sxy = source(ij 3' ：1和謎y为厳蜷点坐标sex = source(2'1 :scy = source(i, 3);加他和乩y为参照点坐标break；endendend觀确辭霍以上代码为获範源坐标系中讎挎点和蔘照点坐标 三血垃=近：*讦算號坐掃系中巻照点和掘菲克坐愫差sdety = scy-sxy;sf兰cetfv(sd«

19、;tK,sdety) ; *讦篇兼塑転系中盖解刖目对于施转点的方位角sd = sqrt (sdeli"2+sdety'2);筋计茸源坐标系中螯照点相对于證特点的距离tdetx 二 tes-tMK;tdety - tcy-txy;+f> = c«tf>(+ deti>tdety):务计耳目転坐桁採中掺赔点相对于龍轶点的方应角detfv =衣计為两坐标義彳目腐于戯转至的矗转甬度td三讯珂(tdrfK+tdety):斛中尊目标坐标栗冲番駅点I目对于诞砖点的距禹k = td/sd-i;帥计算比例系数程倉情况下対0fot i = l：r*sbaSc(i,I

20、= tsoutCe(i, 2'-Szk)*(+k).sbase (ij. 2) = (sciuice (i3 J)-siy) * 11+k):-end.戦纳这段慎玛目的是飛H源坐标系中各点相对于徒辂直的呈纸向里，弁将羽辛大门喩)晤后存储到前包丸矩陆中rovl, collj = siz«(target):中赛需个敎' 以豪暗症憎坏:女数jsq " 0;斷计tftjjl!星fflr i - I:rawlif t吉邙寶仏4)= I 、；.:戢t直：g： ：*亦i拦吕*旳案整忙址土jsq =加 Iif jsq = I病一次霰取独共百|txx = target fij

21、 2);txy = target fi, 3); 命灯和：和対瞌转占土怖«lsetci = tatgtt (ij 2).tcy = target <it 3>十：*it> "丁帯搭也坐厅break;#ndend+nd诵*那谄段代码吕的是兼刖慎坐标希中昔点相対于薛蜡点的基线向軍并椅其护兀H+k)梧后徉储到命豳亡葩陡中 XZ tl* I j = cps<d¥tfw':IZ tlj 2) = -sin'detfir):II(2t1)= sin(fittf*) ,|1Z (2, 2) = cos (detf¥):領融计耳龍

22、無葩阵tbase = IJi*5b35eJ :她先对命吃e中基绒向塑业行直特吏插庐生丽3丸矩阵，艮吕标坐标乗中各点掴対于就转点的基銀冋星匚月驭要注意的是，搭线表中设育点名信息，tltbas中基线末点悟息喪中次陳是一肘应的 jsq = 0； for i = lirpwsname = sourceUj l):占義取基飪末方点名found - 0;辛该点点宕在中存在，那么为“否那么尙0for j = l:rwl Kzai target表中査找当甫抠线耒卓是不是存在if snane tareet (j, 1)*星綾耒点是T?:get表中的点found = 1;乌该和志设为1tflgld町-lirKtn

23、S* (1, U :隔在当前第5刃和第6別填入持頼肩坐标tarettdjei = tsy+tbasf (if 2. *主要杲为了牡査转换是否正碣endendif found =0*14芯自前野主琳亢不呈：2rgZ丟中的产jSQ = jsq+J计数器idtarsctUcwl+jsq, L) = sna»e：MK在target表中未屋增加一行描入無虫 target 'rowL+jsq 5) = txx+tbased, i > ;tar set rcyvL+jsq, $) = txy+tbase (ij 2 ：endend醴唧6将tbd樂申基线向里除回标坐桶帝中的施转点坐桶阴

24、可荻译溥坐标系各点在目标坐标系中坐桶登陆界面:转换界面:导入公共点:Executes on button press in ptishbuttonlBfunction 3 uahbut t on 1C allb ack hObjtctj eventdat色 handle苦' I hObjectN event dat a=K handleshandle to pushbutt dr1 (s&e GCBO)reserved - to be defined m a future version of MATLAB structure with handles and useE data (ssee GUIDAIAjFileName Pathliame = uigetfile(C *a xls' f Encel File b *筈1% 丨"fho节百亡 file1 = length(FileNa».e),if 1<5errordlgC导入文件出错'文件): return：endstr = PathName FileNarn