MATLAB测量平差程序实习报告

I三」测量平差程序设计课程设计任务书专业班级：指导教师：小组成员:TOC\o"1-5"\h\z设计题目2\o"CurrentDocument"设计资料:3\o"CurrentDocument"一、课程设计的目的3\o"CurrentDocument"二、课程设计的任务和内容4\o"CurrentDocument"三、课程设计阶段4\o"CurrentDocument"四、组织方式进度安排5\o"CurrentDocument"五、考核与成绩评定6\o"CurrentDocument"六、参考文献:6\o"CurrentDocument"七、实习报告:5设计题目

边角三角网平差程序设计设计资料:表1已知点坐标JH1=1.点号X坐标Y坐标备注Bj011149.5525529.757已知点Bj021351.0826907.547已知点Bj031466.8006231.200待定点Bj041548.2005403.600待定点表2观测值精度边长观测值固定误差(mm)1比例误差(ppm)1方向观测值(')1表3边长观测值起点终点h、八、、边长(m)Bj03Bj04832.388Bj01Bj03770.170Bj03Bj02685.482Bj02Bj011392.450Bj02Bj041516.951Bj01Bj04418.530表4方向观测值起点终点h、八、、方向观测值(°'〃)起点终点h、八、、方向观测值(°'〃)Bj01Bj030Bj03Bj010Bj01Bj0215.57244Bj03Bj0429.57251Bj01Bj04276.43169Bj03Bj02213.5653Bj02Bj010Bj04Bj030Bj02Bj0415.48098Bj04Bj021.48093Bj02Bj0317.59294Bj04Bj0166.45528、课程设计的目的学生在学习完误差理论与测量平差基础、测量平差程序设计基础等课程的基础上，设计一个完整的测量数据处理程序，培养学生综合应用量数据处理与计算机应用能力，培养学生主动学习，创新设计能力。二、课程设计的任务和内容课程设计任务：在两周的时间内应用者Matlab程序设计语言编制一个完整的边角网严密平差程序，要求有简易的界面，数据输入采用文本输入，采用间接平差模型完成平差的基本计算，能够画出控制网图，输出基本的计算结果，并根据设计过程完成设计报告。程序设计主要内容包括：系统功能设计界面设计流程设计代码书写程序调试三、课程设计阶段准备阶段研究设计任务书，分析设计题目，熟悉原始数据，明确设计内容和要求；制定课程设计计划和进度。熟悉算法模型阅读误差理论与测量平差基础教材，掌握平面控制网数据处理的数学模型，这里主要是指方向观测量、角度观测量、边长观测量的观测方程和误差方程的构成，研究平面观测数据的组织方法，设计Matlab算法，实现计算的自动表达。功能设计阶段设计程序要实现的功能平差程序的基本功能包括数据的输入，平差计算，精度评定、成果输出等；流程和界面设计阶段根据平差计算的过程和程序功能，画出流程图，设计简易界面实现数据的输入和平差计算和成果输出。在此基础上，根据功能要求，设计简便的界面。代码书写和调试阶段按照计算流程图和界面设计，根据方向观测值，边长观测值的误差方程的组成，设计Matlab算法，实现误差方程的自动构成，分阶段书写代码，调试实现各个阶段的功能。设计报告撰写阶段设计报告是对整个设计过程进行综合总结提高，内容包括课设的目的意义、程序设计的内容、算法设计、设计心得等根据设计过程和对测量数据处理以及程序设计的理解进行独立撰写。四、组织方式进度安排以小组为单位，每小组5-6人，分工合作共同完成程序设计任务，时间两周，进度安排如下：第一周周一熟悉设计书内容，制定计划熟悉平差模型周二功能设计流程设计绘制流程图观测值数据组织周三熟悉观测方程误差方程的构成，列出计算公式，设计Matlab算法周四界面设计按照分工进行分模块代码设计完成数据读入模块周五分模块进行调试完成误差方程构成模块第二周周一代码书写调试完成平差计算模块周二代码书写调试完成输出模块周三代码书写调试完成整体调试周四撰写设计报告周五综合考核五、考核与成绩评定考核内容根据考勤、课程设计报告、程序结构、程序代码进行综合评定。成绩评定方法平时成绩20%，设计报告30%；程序代码50%（程序20%，运行正确30%），总评成绩分为优、良、中、差四个等级。设计报告内容完整，图表清晰，代码书写规范，计算结果正确评为优秀。六、参考文献：1.武汉大学测绘学院误差理论与测量平差基础20092.北京建筑工程学院测量平差计算自编资料20113.同济大学Matlab控制网测量平差2006七、实习报告：目的与任务：目的：基于学过的测量平差理论的基础知识，在matlab中编写一个相对完整的数据平差处理程序，旨在让同学们能够将学到的平差理论知识和计算机的MATLAB软件程序编写有机地结合在起来，让同学在实践过程中更加深刻地理解理论知识并尝试去应用，由此提高每个人学生的思考和动手能力。任务：主要任务就是在实习的两周时间内，分组完成一套完整的关于边角三角网的平差数据处理的matlab的程序代码编写与试调任务，其中主要包括：数据组织与输入、平差计算、精度评定，成果可视化、平差成果输出等。然后小组上交一份实习总结和小组成员各上交一份个人总结（设计书及总结）。算法与模型：算法：间接平差方法；模型：平面控制网中的边角三角网本次的平差数据处理是采用间接平差的方法（数学模型）而进行的，就三个待定点xy坐标改正数及涉及到的五个观测站的方位角改正数为未知数（11个），进行误差方程的编写和平差计算，精度评定等。而实际应用是采用所提供的平面控制网中的边角三角网的平差数据模型，分别是就着边长观测量和方向角观测量的观测数据进行平差程序的编写。程序功能设计：通过手动地选择数据文件，而后自动地对数据文件中涉及的变量进行赋值和所获得的数据的输出，而后通过实际情况分析得出需要进行平差计算和精度评定等的必要步骤，如本次数据平差的必要步骤就是待定点近似坐标的计算，紧接着就是误差方程的编写，包括系数矩阵和

观测权阵的求得，而后就是利用最小二乘法的平差计算和精度评定，并且计算误差椭圆的参数和该椭圆的自动绘制，最后平差成果的输出。分模块讲解：数据组织输入和变量赋值已知点坐标信息：包括已知点数n0，待定点数nd，点号名pn和已知点坐标xy0观测值精度信息：方向观测值精度md，边长观测值精度ma,mb边长观测值信息：包括起点bfl,终点bf2，边长观测值L2方向观测值：包括测站点dfl,照准点df2,方向观测值L1具体代码截图：宵槌a。meciclfrHZ■Ji—L:EC■占dfile'IEETfill);仆ununhTECfLI白tvsMidt1-fopentin);ID]2]31ID]2]31HTEl]a£(:'EtnFutFil«arPathirnotcorrect'>'yTainini1□'Tam?rrtura;etid书一—冏的岫惟且寿兀湖mdJitt-fscanftf；LtJ,'?]);的卖H![a]=fscanf[-fit1,Wd*r3),ni上我:1，；£已坷点耒妊直⑵；lidl二边1，\,hp=n&+nd,iori,-];npbpe口诃刑i]=pntiJ如d[ffl]=ijEanf(iitL'^frB2tnJjj;如it1=11mQByajfltijL]=a(2*i-Li;co£iE-XY(](ij2]j,EdaytUliD'yfl'];Cal,seiiptCalE.11QIUtItrt5sHKkrojiis：ni>>v0寻・1>■■•：■：•&漩♦■孔□LUr?EJ]AbdayfJ-riU'yD"];[a]=fsc：anffiH1»'阳'?3*i,35:n2=fi(n3fin^a(3)3F冒-c?adIlin?dirtsrwobEerrtnnn.邮isid=O.3B|i]=f5csrf^fitLB，«iJ2：l:ob5_na=a[lj:obfl_nb=-i[0J;4zdi3t=rsc==cnf(titl3'iifftHr',[3何，43dizs-dist'4dbfkdisi;,1]:■15bf2J:46L2=di2>：!.31.47oL2-dla=[tfLhf2LZ]；4Bdis-obEadis;■IPobEnDb=[b£lb£2L2].5D以h邳b，：L「l=|dfL±E2iL]Edtibfl=MLtJ2];战LML,11]=Ckkdat(tip,削hfl\-欢邑点票点号起点号5Jmr务13]=dikdat帖缶Pm仃舛；*&4nb-tbflbf2：l;55forL?l;h2;/皆重部iiChfiUl=hf2Ci)>；57i3=L;E5Ewtfid的end6D列—4readdirretinnbb/rrwartim61idit=U;V方向喊Sltl■|EcriftLxLCdL?-\SfrXrutFilffirkEllaFilllert互iHKkrojiirniwrm□■■〔。■■•:八&漩♦■孔□LUr?ES3ldiE=O；诘向御[fi~TiKhEwanfIfit]j?>f'jD;63击％rijta=nd.&4lor日ml鸵t】：Li4t#=aidK：fitL•心蹄/3)SOtJTIi：:L)=tUjna(])6F6B(i)=dRf_rad(tlljiaf3*);融wd.花df：i=df：lr;71df2=df2r;TZLJ=LLa；73nsD-tdfld€2：l;74cbs-dir^Edfldf2LJ].75dir=cb5*diE°?6[dflBil]m独日雕让皿gdf1);77[df：,1Z]=c-JtktoWaPmdi2);?Brrrtdfldf2]TP[fileriw.^filejathl='jipu.tfil*Cieru.ltdata.txt'>TEST£i】/),8DitipnunSstr(iilciLsnej;®i也妃二f■叩5〔如§2if[fid2==-l]SJFfc：｛bi：it[[rCnpuTFileorPatJiisnot'■:mnf);&4ie+uxn.35end.$6Ssoperia.file+□vriteariprintf(fi犯.，T&HFttK花宿息■W5；--r(....■：■：.i./：;-..n.-.L□'I：'：；-mrnr~ni.*,住曰京.节:■■.：).3Piprintf^fidS,'W];iprintf'*5膏IIIIOdKIiM'ui、nrt3戒,np』nL;£Uiprirrtf(fid2d'戚J;92fpiiritf<fid2/^);toriplLxLCdL待定点近似坐标的计算通过所提供的已知点坐标计算得出的已知方位角，在就方向观测值，求出其每一次观测的夹角，而后找出每个测站所拥有的已知点（一个观测方向的两个点）确定其所在观测的行数与其他待定点所在观测方向（必须有一个点是已知点）的关系，求出其他待定点所在观测方向的方位角（注意方位角是否合理问题），而后在边长观测值中找到对应的边长，利用极坐标法最终求得待定点的近似坐标°x0y0.同时还可以将近似坐标以文本的文件输出，计算代码截图如下：

agLE.11QIUtItrt5sHKkrojiis：ni>>vD■■|■■■•：，g漩牛*取LJ□匕河匚agL]陀nUbflT]fZ)CZ72岫na^dfldf£)L]3:nflP=fLeU]';AnjdE-trt>+litxplr;5iiifatdii-ChaLIzetcas:(ni.,,1)叮亡roaljiL.l)xtxosC&L.1]]:0cuMdfi1)-7dy-iyOIjiUPtZ],DLrflPLL?)8加];gtnd=fimM血aJ；i■1)=3iJflP(1))4.(na(i2)==»UP(2))];LOtnitzfinrf((naJj.2]=nOP[I]〕d1；i==nOP:L]ihfafdjr(tril：,4)=ta]fi.LZtats=t：=iir：=r+pi.13if+af-3>2*pi,L4tifa=tsfa-2«pL,Liend留ihfLUtx4)-1fifa;XEnri-fcnJ]}]);rpL9trd=find£fnaf;v1)==1)h(na1-j2)=D、;20afarinfDiditU+1Lftnij-Ll(21;2]ilifQt：d3JTa；22toj=eLfid(M:.,1)==ijL[gW)=L):.Z3网iitafa>2+pi,tafa=tafa.~2*pL,霍end27ikfotduf(ttijr4!-1fZBAhd界艾1kQ=EkO(]^xwDsdjL]irrpsdj1]];3】yfl=[yCi[B却顷EC-rofiOjL]zeros1]]:lor1=1rd;V,serijst|LaLCalI9Editiu-GN曲\gTH.・■"固「XDlililtgk互iiKkrojirniwHelpr>mO寻・1>■■•：■：•3漩♦■七□LUr?Ei新59ilU住134.4■IS燹H4B■195062M55牖5?煲6]yO=[yn(]JyD(2)Earoa(]jL)EdTDfl(]j1)];■fori-];rd；ti=ndPti];lrid.=imd((!iribCieJ==d)&[itfa(ij21=ti^|<lnh(EjL]=tis!*^bCh2J=1J]);tS=±is(lEid^S);L)==l.)l(rra(i?Z)==ti]J;alfa-nirfoidie[taLfaidjfl).a£i(tiJ=JcO[l)i-t5>coEihHlfi>,?Cifti)=70(l)H;dt=]■forj-]：nd;ti=n.dFtj];iori-2it-i-lt：m=ELtidi：(if1)—ill2)=tlJ);t：rij=t：Lfid(g1：!.,Ij==i)h2)==L)iaf-srinfDidir(trik,4l?+£L(tnda1~L](tnk).infotdiTtrii=findj；aJ；H1)=tij&(naIjW；i==d3)taf*=af+pi；iftifa>3vpi;t：=ifa=taf：=r2*pl；?ndinfDidiiCtnjP=ta£jend.t-t4-L:ift—up-]brvikend.etid7Q|Ecrip-ILxLCdL山功IdtI«rtSba.«krojli:ni->vHelpkm□Li^■IKJ•・rf旨漩**虹□LUr?Elafa^tac£s-2»pi;&4tifdEH5irifaidiittiij；,^)=t：aCa;&bend57t=+tl,&S£Sift—up-100bt皿blend成eni亮2thewftd堂.JLfilJ.：随比虬糠tt=3x£l'yiT]：scS5'n=[pxian'河1；67^hn'yTi'].[fileriMi^filejath]=ijdpu.tfile-CxQyO_d3ta.txt'jTEETfile'),@9itipnui25trfiLiensnel;70折±1=£咖wGn'irt’n?]it(fit：!-=-]]'mrhoH：「[npmFiLflnrPa+hisnMcocrflcV,"aTrim/jj'mrn「)；ie+uxn;end也吐ntf〔旺t，’=-一近帆坐桥一-=赋)；tjsrlinrCLVJ/序号点号E：n)TOWlh");fiprintf(fit1,Yr');fnri=j-lwiethLxOjfprurtf(fitJji'阳源】1dK2LSft2Q.Hin、%叩[i]BxOCi).yU(ij);tlid田trprtmt[txt,_\n、\S3EclpsrLfit]?:|EcrifiLxLCdL误差方程式组成本次平面控制网平差模型中观测量为方向观测量和边长观测量，未知数分为待定点坐标和定点角两类未知数，边角网误差方程的形式如下：方向观测量:L^+v^=&-Zj误差方程：（系数为：b1,c1）4△AAAAAv——Z+a—L——z+ax+by—ax—by—ljkjjkjkjjk狄jkjjkkjkkjk式中常数项：I*—L*+Zo-a°k变长观测量：L.—L+vi=&Xk-乂,）2+（「Yj）2误差方程：（系数为：b2,c2）AXoAAYoAAXoAAYoAv=Kx.-——Ky+监x+——总yiSojSojSokSok其中l—L-So=t'（Xo-X）2+（Yo-Yo）2iijkkjokj误差方程的编写要注意两个问题：系数的位置确定和是否待定点的判断。对于方向观测值，未知数不仅包括各待定点坐标（6个），同样包括各测站的定向角，在未知数排序中，定向角未知数在前，待定点坐标在后，则未知数个数为nz+2*nd=11,对于方向观测值，每一列与定向角对应的系数均为-1，与坐标对应的系数即每一列有5个系数（最多），其他系数均为0，及对应的系数矩阵的数值的确定（注意正负号）。而对于边长观测值，每一列包括边的两个点的对应的改正数（如果两个点都是待定点的话）四个参数，其他均为0，同样还要注意系数的位置问题。由于这里分别对方向观测量和边长观测量设置两个数组，分别存储系数阵常数项变量，这里设置blb2clc2四个数组。根据观测值序列和未知数序列确定系数矩阵中各元素的位置和大小；这里将误差方程构成设计成一个function函数文件，输入量为已知点信息，近似坐标，边长观测量，方向观测量，等输出量为误差方程系数，常数项。代码编写截图：

3Editi>r-mEEffXE.11QIUtItrt5sHKkrojiis：ni>>vD寻.I出J•・rc3漩♦■兀□LUr?|E21bl<ijdir(iP2);-a22b1(i,diiZ3ard网?n.db2=.[itco!：(n%],蹈fari=l:n227血二邛WiHi.U.幻.1);ZEdy=iy(dis(ij1),-uy>：di2(Ij舞sGj.11-sQc-t(da^dz+<lv*d?r>;30ifdiH姑L)侦顶"uNi疽MM扪b2<iJdir(ielh2)=-dr^B(i,L);32(□.,dir(ir1J*2+1]=-dy/b(i.1):30elstifdis[i,LKi皿Hl2(in2)hiO孰L2dtfi(i?2)»2)=d^/a(I.,l：l;36b2知由工(i?的*2+1)■面&&]:';36elsflifdis[i,L)<adUdiFti9<hd37t2<iji]=iJ；签el^t翌liSdj(ijL);40bid.L5])*2+1)・-此，剥Lp1);4」bH土dLs&勤蛇扯心此虬1.〕;dir(ie卫A家1)&1>;必;end<4end15tar1=]?jiz4后brrarfind(infqrdi工L]"iSinfor±LJC;?5i'--n?,47ZLLEi-j^*Cni-L)+l;Cnz^1)1]=in.fordixCln?u;皓endcL-L1+21-1^0^11(2^1;51nn.d52LxLCdL5.平差计算和精度评定本次平差的特殊就是结合两个误差方程的数据进行同时平差，具体步骤在得到误差方程的系数之后，根据方向观测值和边长观测值的初步计算观测值权阵(plp2)，严格的精确计算可以采用方差分量估计方法精确计算方向和边长观测值的权。在得到误差方程之后，可以直接利用矩阵进行方程组的解算(最小二乘法)，误差方程分为方向观测误差方程，边长观测误差方程，两者可分别计算，然后再相加。由于近似坐标不精确，在求取参数时，程序设置了迭代计算并确定终止条件，在改正数满足条件后，即可进行平差值计算和精度统计。程序截图：13*Editor-^:\SJr\jincdu.■|•X23a667891013121314ISie171819202122Z32425272S29303132*4;Vidf-find<LnfordiiC:731—01;ro=36COtl80/pi.nx”・S<j>FJ)y;dxk=xeios(ra.l>;pl=X9XO2《H，nl);p2=2«ro2<riZ,n2>;yl-3oxo5<n]/IJ,v2-zoio5<n2/1).scnO=cbs.rbta.wobaaa;ib=ob5.址：fori=3:nlpltl.0=1;oiviMlfori=l：n2abi=dir(i,5)*rb/LCOOi5i2-ria'2*nbi'2p2(i,i)=a^no'2/>3i26ivdM2forkk=1:8lb，cL.^2,c2]^plcb^e^CnO.nl.LLdix,r>2.L2,dis.xy,inf&vdir)atpa=br»pl*b!4bZ«p2«b2i»tp1=br»pl»cl-H)Z«p2«c2<^X=<inv<axpa)>»atpl.v.)-bl*dxk-Gl.v_2=b2<dxk-c2.Lpj=LL4v_L：Lp2=L24v_2.Mt=[rifidxk]bz—1;fork=nOtt:np,j^z*-<k-n0>*2;xyli.0=x7l>.l>4dxk<j-1)/300.teripl田卬曰5|口F11S*Fditor-jincdu.■仁~l耳又2S29303132335435为3738¥40414243444546474842W5)52S3S4WS657SS69Mtt=[Ziff心2bz—1,k=nOtt:nj>.g:Hk-nD)・2;g,L5=x7l>.lMdxkO-1)/300;町(k,2)=x^2)-htok<j)/IOO.oivd^kfork=nsM:nxif&bffldxkl»)<O.l;bx=l:en偷ifbx>0breake*mfendikkCoti=l:nlLp](iJ-rad_d>i5llpl(i)Jendgy&IJ;Yp=xy<!.2);lfYp=xy.txrp!:rd;obad”p=hufriaOLlCv_lLpJ];vt=*bl>dxk-cl.v2=b2>dtX-<2;HvI.WJ;pw=u「<pL<vl*v2・»p2»92；nOJ-IA'v/<nI»nd-rra>;nQ=冲tM2),<|xx=iriv(atpal;Dxx=rD®i|xx.certpl田卬曰5|口Col误差椭圆参数计算参考书上的误差椭圆参数计算模型由于误差椭圆描述了点位精度及其在各个方向的分布，表示了网中点或点之间的误差分布情况，常用误差椭圆对布网方案做精度分析，误差椭圆三个参数为椭圆的长半轴E，短半轴F和主轴方向本次参数计算是通过单位权中误差，协因素阵，待定点的协方差阵等上一步平差得出的结果或推算结果，来进行误差椭圆确定。代码截图：

"9Editiu-ucijum.u._酒XE.ilaFilllert互iHKkrojiirniwFlMD■□LUr?E]]节==========薰7Q2帷善棉E3朝计算QxaQn(np+JI部山.r(r+Li«ndl;50fari=]2sid7iimDH;3gOTfcHXL2*i);LOirndEm-LW蚯J;L]LZRswqy；13yy"砂，L4II=OK+Yjy,Liipc'grt(ez);留if[alsarisJ<1d-]0]L7fl[i^fiigriiDO.Ojqay).18bIspL9fiti)=staniyy/xK)*57»2058,20end2]if[n>=ILLSyy>=|].0)22fi[iJ=fi<£>/2-0：Z3slseifiis>=DnWyy<=Dr0：1网elspif(xx-40.0)霍fii；iJ=i：fiCi^+L$fluCi^2.y；27errdZBai(1^=nDsi£qin:(-,阁hi,Li?-nO*^qrtC£qK+T.v-kX^2-.呻«艾1abitil=aiCjJ~td.U];in1iJ^D^jqrt；»/li)=nfl>scjkrt(qyl.,serijst|LaLCalI控制网图形绘制和误差椭圆绘制边角三角网控制网图形绘制主要考虑绘制控制网点位分布图，点位标注，边长观测量和方向观测量的绘制等。基本思路是首先根据坐标绘制控制网散点图，并进行标注，然后根据距离观测值，方向观测值连线，构成完整的控制网图。

至于误差椭圆的绘制则是在控制网图的基础上再加上待定点对应的误差椭圆的绘制就完成了。（利用Plot命令）其中有一个特别注意的地方就是参数的设置，要根据实际情况而定。绘制代码截图：平面控制网成果输出以生成txt格式文档的形式进行成果输出。其内容主要包括五个方面内容：平差计算后的坐标值，待定点点位中误差，误差椭圆参数计算等。信息具体内容如下：平差计算成果成果总输出；坐标平差成果输出点位误差椭圆参数误差椭圆的绘制方向观测平差成果：方向观测值，平差值，改正数；边长观测平差成果：边长观测值，边长平差值，改正数；代码截图（其中的几个）：3Editiu-■G-;^SfT\Finvl.ai■X□

平差主界面的设计此界面为本组实习成果的最大创新之一，就是我们在原本要求的平差程序界面（边角三角网平差程序）的基础上增加了另一个界面，就是小组信息的界面（包括有小组成员照），而两个界面之间利用一个总的界面（小组平差程序设计）进行选择性打开代码截图：主界面一小组平差程序设计□LUr?EF-V总的设计流程图:代码书写与调试：由于上面已经展示过几乎所有的代码（根据每一个模块），真正运行整一个程序的时候却是按照步骤一步步地一个模块地运行，但是问题一定会出现的，而且还不是很少，关键就是怎么样通过一步步地发现问题所在，并根据实际情况去试调程序，合理地按照平差理论去运行程序。一下列举了一些出现的主要问题。方位角由于所提供的数据时边角三角网的数据，这样就会涉及到方位角问题，而且还不止一次的出现方位角的运算（如待定点的近似坐标的求解过程中出现过），这就应该惊醒一个重要的细节：就是所求的方位角是否合理？怎么样去判断呢？就是通过if（elseif）语句去判断其与0和360度或者2*pi的关系，在去改正方位角。原则上一旦出现方位角的求解出都是要用这个语句去判断，这样才能保证精密。和提高正确性。单位问题单位问题出现在两处，就是观测值的单位权的求解处，二是误差方程的系数求解处。由于这两者是在后面的法方程求解时都要用上所以，在此之前要保证两者的单位相对统一，这个统一要分两处，边长观测量方程系数和方向观测方程系数，特别注意的是原始数据提供的观测精度（两种）的单位，和误差方程的系数的单位，而在实际中，我们组应用的是方向观测用秒为单位，而变成观测以米为单位，这样出来的结果最后的平差数据处理，如精度评定，误差椭圆计算等，也要注意单位问题。比例参数根据实际求出的误差椭圆参数（EFfi）的大小和所绘制的三角网的大小，这两者之间通过比例参数的调整才能更好地在显示误差椭圆。弧度转角度对于弧度转角度这个小问题，我不得不提，这是一个关键问题。为什么？因为我们组把其他问题都解决了后，组的坐标改正数，误差椭圆参数等都很不正常的显示在电脑屏幕上，一次次地打击我们的信心。最后，我在仔细地对照书本，细心查询，观察异同发现，一个最不关紧要的差异，就是弧度转角度的小问题，我们原来采用的是MATLAB系统自带的弧度转角度的函数，结果我们发现这个问题直接就是导致待定点的近似坐标差异特别大，尤其在待定点3的x,y坐标上面，差了好几米（跟提供的平差结果数据相比），这样后面的就影响到误差方程系数的不正常，和平差计算得出的改正数，单位全中误差，点位误差和误差椭圆参数的异常。而这样我不得不自己再编两个弧度和角度相互转换的函数文件。于是重新采用转换函数程序之后，我们的数据就非常正常的出来的。平差程序设计成果展示进入界面13叵］区I导入数据文件HFigur^2=平面两平差计算程序读人数据文件近仍坐标计耸平差计食误差椭圆计耸篷制网囹显示计耸结果退出欢迎使用13叵］区I导入数据文件平面网平差计算程序小组成员：王子强，缪云龙，王萝琦，闪永川，王瑞,

宁子铭数据文件的查询查询结果：Prvffiil1dDtQ.txt-Q回区1fJr299,5®一臀方向观测值（卜邳11WBuSaMC2150.?7®5fl用14JlrSZT?®1|21e.a»Me520u275S1623731n.99Mey34e.E2£B5y23.73416193fl.OiBMf11W1211.1AE2A-巽侦•穿％'近似坐标的计算计算结果：控制网点为分布6907.5>i66231.8185403.503—迁船坐标一XO(n))114y控制网点为分布6907.5>i66231.8185403.503—迁船坐标一XO(n))显示平差计算成果平差计算成果总输出点位误差椭圆参数中站站

位大均

兽率mmP-tuoyuan_data,txt-记事本文件（X）编辑©椿式①查亏⑦越助叫中站站

位大均

兽率mm--误差椭圆参教一pn(i)mxmymnabFi32.3666.0772.0036.2BJ|105.37842.2632.3403.2552.6531.88647.905左=1上25£19差mm-3.25E(cn)实习总结:两周的实习历经百般波折磨难，终于结束了，这两周的时间里，大家经历的是前所未有的实习体验，让我们不单单收货了专业能力，更是给我们的人生上了重大的一课，在课堂上我们都是学习的书本上的理论知识，这次我们把理论知识和计算机编程结合在一起，这才是真正学习到了未来生活工作中所实用的东西，提高了我们真正面对测量平差工作还有matlab软件时的操作能力。我们在实习中发现自己的问题并及时的解决，让自己在以后的学习工作中不犯同样错误。在刚拿到实习任务的第一天我们就进行了分组，分完组后，为了能够更快更有效率的完成此次实习，也为了能让组内的每个成员参与进来，首先，我们在组内对各自的任务进行了分配。大家各自负责一部分的程序代码编写内容，各自针对自己所主要面对的内容查书并向其他各组询问讨论，然后根据书上的代码提示，调通自己负责的这段程序代码，每个人把自己的部分完成后，大家在一起编程查错，经过反复修改弄出整个的成果输出，还有简单的界面设计。我们的组长是王子强，他和王瑞负责的是主要的数据计算部分还有整体修正还有精度评定部分，这几点可以说是这次程序代码中的核心所以需要他们两个人分担重任，宁子铭负责的是误差椭圆参数的计算和图形的绘制内容，缪云龙负责的则是数据文件的输入部分的程序编写，闪永川负责的是最后的成果输出部分，王梦琦则是负责最后的界面部分程序代码的设计，这些分工只是各自主要负责的项目，但是因为大家能力都有限，所以出现问题后还是集合在一起解决问题再继续进行之后的实习。一开始，先是将已知点信息录入，这块代码是由我们组的缪云龙主要负责的，这部分的代码看似简单，但还是给我们带来了不小的麻烦，需要在matlab中单独建立一个以'chkdat’命名的‘function’文件。这些看似简单的问题却也给我们带来过不少的烦恼。接下来就是由我们组长王子强和王瑞同学负责的近似坐标计算和相关精度评定的代码编写，这部分非常的复杂，它包含了方位角计算，编写误差方式的程序代码等或大或小的问题。遇到的首要问题就是我们看不懂书上的近似坐标计算的代码，而书上的代码又存在问题，缺少nk的定义，而想要修正这一点需要读懂整段代码并将nk通过书写给予定义。还有就是误差方程的系数矩阵的编写程序。这里有好几个因素要考虑：方向观测边的方位角计算，未知数的系数位置安放问题，还有就是点号是否是待定点的判断等等诸多问题。在不断地学习中向其他组成员请教后我们总算输入了正确的代码修正了这部分工作。算完上面这些计算之后，就是由宁子铭负责的误差椭圆参数的计算和图形的绘制了，这部分的代码编与之前的计算相比就显得简单了一点，但由于我们只能参考教材，而教材的程序又多少有些不适用，程序还是出现了问题，在大家共同的努力下，在同学的帮助下我们还是完成了程序，但是，由于我们数据的单位不统一，导致我们的椭圆参数一直较大，调节了很久才找到问题的所在，纠正后总算是可以了。而图形的绘制相对就简单些，主要都是将计算与程序相连，再把x，y给予定义就可以了。出图也完成之后就是由闪永川负责的成果的输出，这部分不难，我们很顺利的完成了，当然这也源自闪永川之前的准备比较充足，让全组的效率得到了提高。最后的就是，由王梦琦