矿大测绘专业实习报告

1、大地测量基础与卫星导航与定位技术教学实习报告书 学院 环境与测绘学院 班级 测绘11-3班 学号 07113021 姓名 田孟浩 2013-07-23一、实习目的与要求1、进一步巩固、加深大地测量学、GPS的有关理论知识。2、熟练掌握大地测量学、GPS所使用的仪器及观测方法，提高控制测量实践能力。3、培养和提高解决实际问题及组织测绘生产的管理能力。4、进行热爱测绘事业的教育。二、实习内容与任务1、进行选点、埋石的实践训练。2、合作共同完成若干GPS点的静态相对定位观测。3、全组完成全站仪、一台水准仪及一对水准标尺的检验和记录。4、全组共同完成一条导线的水平角、垂直角的观测工作。每人独立完成至少

2、2个测站点的水平角和垂直角的观测工作，每人独立完成至少2个测站点的全部记录工作。5、全组共同完成一条导线的测距工作。（如用全站仪可同时进行角度与距离的观测）6、全组共同完成一条导线上若干测段的二等水准观测工作。每人都有不少于10站的观测和记录工作。7、每人独立完成一份平面控制网和高程控制网的概算及平差计算。8、每人编写一份实习报告书。三. 实习的组织和安排：测绘11-3班共分为五组，每组共6人，然后每组出一人参加由全专业共同组成的G-S组测量任务。每组设正副组长各一人，正组长负责分工安排全组的实习工作，由副组长兼任仪器保管员及时检查督促安全工作。实习宏观上由教师统一指挥，班干部积极协助指导教师

3、做好各项工作。我是在测绘11-3班实习三组，两个组长，四个组员，共六个人，分别是田孟浩、王名亿、蒋路、张学万、张宝科、张梦虹。四、执行规范名称（1）全球定位系统（GPS）测量规范GB/T 18314-2009。（2）全球定位系统城市测量技术规程CJJ73-97。（3）工程测量规范GB50026-2007。（4）城市测量规范CJJ8-99。（5）国家一、二等水准测量规范GB/T 12897-91。（6）国家三、四等水准测量规范GB12898-91。（7）国家三角测量规范GB/T17942-2000。五、实习地区一般情况介绍1. 实习地区所在地理位置实习地区为矿大南湖校区教学与生产实习基地五号区域

4、。GPS点有两个，分别为I24和I25点，基本上是绕着图书馆和镜湖区周围的闭合环。2. 交通情况我们的测区有三边是公路行走方便，但车辆较多，对我们的观测产生了一定影响。由于我们的测区四周与其他测区相连，导致我们的控制点经常被其他组使用。不利于我们的持续观测3. 地形地物（包括植被）概况测区有两面有着明显的地形起伏，高差较大，同时植被茂密，经常出现不通视的问题。并且测区内有一条河，河岸上的树有时会对测量造成一些不便。4. 已知控制点利用情况已知点有两个I24和I25，这两个点通过GPS组的观测，用Compass软件将基线进行解算可以得到这两个点的XY坐标。六、控制网的布设、作业计划和作业组织1.

5、 控制网（GPS网、导线网、水准网）示意图（1）GPS网示意图：2、测区分布示意图3、起算数据点号XYH（56高程）备注泉山（QSGY）3788109.4670514667.667755.950054坐标E0883787813.059513642.03439.966454坐标3787766.762513589.52039.96980坐标S0036.607概略高程七、参加选点在实地确定控制点的具体位置，力争做到点位安全可靠，便于观测和保存，点间相互通视，控制网的图形结构良好。在实地设立测量标志，埋设标石。其几何中心表示点位。八.仪器检验1、全站仪检验： （1）照准部水准轴应垂直于竖轴的检验和校正

6、检验时先将仪器大致整平,转动照准部使其水准管与任意两个脚螺旋的连线平行,调整脚螺旋使气泡居中,然后将照准部旋转180度,若气泡仍然居中则说明条件满足,否则应进行校正。校正的目的是使水准管轴垂直于竖轴即用校正针拨动水准管一端的校正螺钉，使气泡向正中间位置退回一半为使竖轴竖直，再用脚螺旋使气泡居中即可此项检验与校正必须反复进行，直到满足条件为止。 （2）十字丝竖丝应垂直于横轴的检验和校正检验时用十字丝竖丝瞄准一清晰小点，使望远镜绕横轴上下转动，如果小点始终在竖丝上移动则条件满足否则需要进行校正。校正时松开四个压环螺钉（装有十字丝环的目镜用压环和四个压环螺钉与望远镜筒相连接。转动目镜筒使小点始终在十

7、字丝竖丝上移动，校好后将压环螺钉旋紧。2、水准仪及水准尺检验：（1） 光学测微器效用正确性的检验和分划值的测定步骤：架设水准仪，将仪器视准轴瞄准水准尺，使用倾斜螺旋使水准气泡影像精密符合，在一测回内要严格保持倾斜螺旋的位置不变。旋转测微器螺旋，使测微器读数略大于零，指挥扶尺者将特制分划尺上下移动，直到特制分划尺上某一分划线被楔形丝平分。然后将特制分划尺固定，并在一测回中保持此位置不变。往测旋进测微螺旋，用楔形丝分别夹准特制分划尺上两相邻分划线，读取并记录分划线编号和相应的测微器读数。返测按相反的次序，旋出测微螺旋，使楔形丝分别夹准特制分划尺上两相邻分划线，读取并记录分划线编号和相应的测微器读数

8、。至此算作一个测回。每测完两个测回后，应略微移动分划尺或变更仪器高度，使每两个测回分别观测分划尺上不同的分划间隔。（2）、视准轴与水准管轴相互关系的检验校正：、i角检验：设水准管轴不平行于视准轴，它们之间的交角为i。当水准管气泡居中时，视准轴不在水平线上而倾斜了i角，水准仪至水准尺的距离越远，由此引起的读数偏差也越大。检验时，在平坦地面的一直线上选定J1，A，B，J2四点，点间距均为20.6m， 具体检验结果如下：i角误差的检验对于精密水准仪，要求i角小于20。超出限差时，必须对其进行校正。算出i=12<20,所以检验合格。九、GPS静态相对定位测量与计算成果：1、外业观测：（1) 调度

9、安排，确定每台接收机观测的测站、开机时间、搬站交通布置等。 （2) 观测组按调度表规定的时间进行作业，保证同步观测同一卫星组。（3) 每时段开机前，作业员量取天线高，并及时记录（或输入）测站名、年月日、时段号、天线高等信息。关机后再量取一次天线高作校核，两次量天线高互差不得大于3mm，取平均值作为最后结果。若互差超限，应查明原因，提出处理意见记入测量手簿备注栏中。（4) 仪器工作过程后，作业员应对照指示灯工作状态说明，判断仪器工作是否正常。 （5) 一个时段观测过程中不得进行以下操作：关闭接收机又重新启动；进行自测试(发现故障除外)；改变卫星高度角；改变数据采样间隔；改变天线位置；按动关闭文件

10、和删除文件等功能键。（6) 观测员在作业期间不得擅自离开测站，并应防止仪器受震动和被移动，防止人和其他物体靠近天线，遮挡卫星信号。（7) 接收机在观测过程中不应在接收机近旁使用对讲机；雷雨过境时应关机停测，并卸下天线以防雷击。（8) 应记录雨、晴、阴、云等天气状况。2、内业处理与平差计算：（1）数据传输；（2）数据预处理；（3）基线解算；（4）闭合差检验；（5）平差（自由网平差和约束平差）；（6）网平差成果报告。3、成果（1）平面坐标站点北向(x)/中误差(m)东向(y)/中误差(m)中误差(m)误差椭圆E(m)F(m)ET(度:分:秒)E0883787813.0590*513642.0340

11、*I0083787764.86480.0025513107.40210.00270.00370.00290.0023121:49:06I0093787637.20490.0047513006.96260.00510.00690.00530.0045121:55:46I0243787894.72700.0062512721.41030.00570.00850.00700.0047141:41:51I0253787992.32660.0070512587.96850.00610.00930.00750.0055148:05:22I0433787313.99260.0092512404.42800.

12、01150.01470.01150.009294:10:01I0453787466.54370.0102512392.40700.02070.02300.02220.0063112:07:59I0683786818.10160.0050512808.93240.00620.00790.00660.0044117:16:56I0743787778.56090.0090513570.51580.00740.01170.00900.0074176:07:29I0753787935.43510.0094513689.11150.00870.01280.00940.00878:08:46I6737866

13、57.22130.0053512678.78920.00610.00800.00670.0045124:14:35泉山点3788109.4670*514667.6677*2 基线向量及改正数 基线起点->终点.时段DX/改正数(m)DY/改正数(m)DZ/改正数(m)距离/改正数(m)中误差/相对误差(m)I025->I075.191B-955.5505-434.5615-93.10521053.84490.03240.0088-0.00420.0001-0.00631: 32550I67->I068.191H-74.4072-140.6239132.3218206.9314

14、0.0134-0.00120.0009-0.0005-0.00051: 15429I025->泉山点.191B-1788.7431-954.855959.34162028.51500.0171-0.00160.00190.00450.00061: 118946I67->I045.191H459.2072-270.4475673.0827858.51820.03860.0084-0.0091-0.00090.00671: 22243I025->E088.191B-945.3791-350.8361-193.50631026.77760.0174-0.00150.0012-0.0

15、0580.00211: 59125I67->I043.191G409.7618-199.8845546.6210711.79560.02650.0034-0.0129-0.00700.00021: 26888I025->I008.191B-482.2044-82.2523-232.0504541.41840.01730.0025-0.00660.0108-0.00581: 31253I67->I008.191G-97.0160-753.6495912.79751187.68640.01320.02810.0102-0.0257-0.02861: 90314I025->I

16、009.191B-425.472327.6369-337.3635543.69530.01890.0004-0.0032-0.01820.01081: 28816I67->I009.191G-40.2839-643.7602807.48441033.48000.01460.00230.00180.00280.00101: 70927I025->I024.191C-142.2979-14.6311-82.9673165.36730.0201-0.00010.00330.0036-0.00201: 8212I67->泉山点.191H-1403.5547-1626.25301204

17、.18942462.66860.01260.0015-0.00720.00910.00841: 195312I025->I074.191B-890.5598-300.8946-221.8900965.85170.0256-0.0018-0.0014-0.00320.00291: 37732I67->E088.191H-560.1907-1022.2333951.34151504.60130.01240.0009-0.00420.00510.00571: 121175I075->泉山点.191B-833.1926-520.2943152.4467994.06040.0296-0

18、.01240.01440.01780.00561: 33585I068->I045.191H533.6144-129.8236540.7609770.72740.0390-0.0008-0.00330.01040.00731: 19752I075->E088.191B10.171483.7254-100.4012131.12510.0297-0.01330.01080.00460.00231: 4417I068->I043.191G484.1689-59.2605414.2991639.98060.02680.0072-0.01080.00140.00731: 23838I0

19、75->I008.191B473.3461352.3092-138.9452606.20460.02960.0016-0.00330.0121-0.00351: 20497I068->I008.191G-22.6088-613.0255780.4756992.70020.01340.02070.0099-0.0240-0.02551: 74289I075->I009.191B530.0782462.1985-244.2583744.49480.0306-0.00280.0201-0.00660.01261: 24355I068->I009.191G34.1233-503

20、.1363675.1625842.70690.01500.00380.00260.00410.00191: 56141I075->I024.191C813.2526419.930410.1379915.32740.0313-0.00510.0021-0.0005-0.00361: 29233I068->泉山点.191H-1329.1475-1485.62911071.86762263.32210.01280.0026-0.00740.00930.00771: 177442I075->I074.191B64.9907133.6669-128.7848196.66260.0348

21、0.0035-0.01050.0003-0.00621: 5646I068->E088.191H-485.7835-881.6093819.01971297.69560.01260.0030-0.00440.00670.00611: 103286泉山点->E088.191B843.3640604.0197-252.84791067.72410.0086-0.0016-0.0016-0.01170.00061: 124423I045->I043.191G-49.445570.5630-126.4617153.02470.0444-0.0041-0.0032-0.00510.00

22、401: 3450泉山点->I008.191B1306.5387872.6035-291.39191597.93290.00800.0042-0.00760.0070-0.00201: 199783泉山点->I009.191B1363.2708982.4928-396.70501726.60770.01270.00180.0019-0.01600.00621: 135882I045->I009.191G-499.4911-373.3127134.4017637.90090.0389-0.00650.01240.0063-0.00091: 16400泉山点->I024.1

23、91C1646.4452940.2248-142.30881901.33010.01600.0013-0.0010-0.00260.00091: 118969I045->泉山点.191H-1862.7619-1355.8055531.10672364.35290.03850.0016-0.00490.00030.00161: 61463泉山点->I074.191B898.1833653.9613-281.23161146.07580.0226-0.0001-0.0029-0.00810.00021: 50678I045->E088.191H-1019.3979-751.785

24、8278.25881296.83530.03820.0024-0.00320.00680.00151: 33925E088->I008.191B463.1747268.5838-38.5440536.79950.00860.0054-0.00410.01880.00131: 62715I043->I008.191G-506.7778-553.7650366.1765835.20330.02620.0172-0.0065-0.0316-0.02001: 31935E088->I009.191B519.9068378.4731-143.8571658.96870.01230.00

25、560.0017-0.00790.00721: 53507I043->I009.191G-450.0456-443.8757260.8634683.82490.0265-0.0044-0.0152-0.00090.01251: 25821E088->I024.191C803.0812336.2050110.5391877.60590.01630.0027-0.00140.00510.00261: 53793I043->泉山点.191H-1813.3164-1426.3685657.56842398.96640.0260-0.0050-0.00410.00720.00821:

26、92260E088->I074.191B54.819349.9415-28.383779.40370.0229-0.00190.00350.0079-0.00191: 3474I043->E088.191H-969.9524-822.3488404.72051334.49010.0257-0.00900.01170.01810.00481: 51831I008->I009.191B56.7321109.8893-105.3131162.43470.0128-0.00570.0133-0.02070.02051: 12668I008->I009.191G56.732110

27、9.8893-105.3131162.43470.0128-0.0150-0.01310.0328-0.03541: 12668I008->I024.191C339.906567.6213149.0831377.27290.0163-0.00720.0208-0.0039-0.00431: 23094I008->泉山点.191H-1306.5387-872.6035291.39191597.93290.0080-0.0207-0.02150.03120.03431: 199783I008->I074.191B-408.3554-218.642310.1603463.31610

28、.02280.00190.0012-0.0186-0.00261: 20325I008->E088.191H-463.1747-268.583838.5440536.79950.0086-0.0273-0.01260.03330.03221: 62715I009->I024.191C283.1744-42.2680254.3962383.00360.0181-0.0051-0.00050.02100.01031: 21111I009->I074.191B-465.0875-328.5316115.4735581.01080.02420.0057-0.01020.00820.0

29、0281: 23991I009->E088.191H-519.9068-378.4731143.8571658.96870.0123-0.0011-0.00760.00190.00561: 53507I024->I074.191B-748.2619-286.2635-138.9228813.10650.0249-0.00160.00520.0015-0.00061: 32602泉山点->E088.191H843.3640604.0197-252.84791067.72410.00860.00020.0023-0.00280.00211: 1244233、高程拟合成果站点高程(

30、m)高程(m)E08839.9664I00841.14070.0346I00941.38220.0349I02446.51460.0359I02550.26550.0361I04351.59080.0385I04552.11880.0424I06844.90180.0349I07440.03500.0380I07538.62500.0414I6745.93510.0348泉山点55.9500十、水平角观测、内业计算与成果1.观测方法。水平角用DJ2经纬仪观测2个测回（注意瞄准觇板上标志中心或地面导线点中心）。每测回度盘位置分别为：0°00，90°30。每测回各方向2C互差1

31、3。角度测回互差限差为9。2.导线边改正后的斜距化算为高斯平面边长；测角中误差4；测距中误差12mm；计算导线方位角闭合差（限差为10n），合限并分配角度闭合差；计算导线坐标闭合差（限差为1/15000），合限后转入平差计算。提供各控制点的54坐标。精度评定：=+20秒<； =0.069m；K=0.0667/1555.9987=1/23328<3.水平角原始数据站点觇点读数左-右(2C)左减右(2C) °半侧回方向 °盘左°盘右°方向值I24I250 00 00179 59 53-7268 10 47 268 10 50I23 268 10

32、4788 10 45+2268 10 52I24I2590 30 01270 30 05-4268 10 41268 10 44I23358 40 42178 40 52-10268 10 47I17I180 00 00179 59 50+10268 36 36268 36 41I16268 36 3688 36 35+1268 36 45I09I1890 30 00270 29 47+13268 36 40268 36 37 I64359 06 40179 06 2020268 36 334.水平角观测成果导线计算代码：clear;%清空内存yzds=2;wzds=15;%yzds是已知点数

33、，wzds是未知点数gcjs=17;gcbs=17;cjzwc=1.5;cbzwc=2.0;cbblwc=0.0;%观测角数=12；观测边数=6；测角中误差=1.5秒；（测边固定误差、比例误差，以cm为单位）JKH=-1,-2,1;1,-1,2;2,1,3;3,2,4;4,3,5;5,4,6;6,5,7;7,6,8;8,7,9;9,8,10;10,9,11;11,10,12;12,11,13;13,12,14;14,13,15;15,14,-1;-1,15,1;%观测角的端点编号,顺序是测站，后视，前视。JK=-2,-1;-1,1;1,2;2,3;3,4;4,5;5,6;6,7;7,8;8,9

34、;9,10;10,11;11,12;12,13;13,14;14,15;15,-1;%观测边的端点编号YZDZB=3787894.7270,512721.4103;3787992.3266,512587.9685;%已知点坐标赋值JSWZD=-1,-2,1;1,-1,2;2,1,3;3,2,4;4,3,5;5,4,6;6,5,7;7,6,8;8,7,9;9,8,10;10,9,11;11,10,12;12,11,13;13,12,14;14,13,15;15,14,-1;-1,15,1;%未知点近似坐标推算路线，即测站、后视、前视；前视点就是需要推算的未知点，此时测站和后视点坐标已知或已经推出

35、。gcjd=268.17972,190.37583,205.09,107.91889,91.81528,258.61972,100.19028,171.92222,141.25639,179.33583,170.96639,100.06417,174.96361,268.61,110.55972,108.42306,139.89083;%观测角度gcbc=165.3265,191.6509,113.7752,62.6204,79.7123,174.0970,102.8376,120.5101,54.2983,66.2271,97.5891,36.6804,119.3408,171.8028,7

36、2.4477,50.5438,41.8650;%观测边长ro=2062.65;pi=3.141592654;%角度化弧度的常数for i=1:yzds; XYJS(i+wzds,1)=YZDZB(i,1);XYJS(i+wzds,2)=YZDZB(i,2);%将已知点坐标值赋给另外一个数组，此数组的前面将存贮未知点坐标.endfor i=1:wzds;%开始计算未知点坐标 j=JSWZD(i,1);k=JSWZD(i,2);h=JSWZD(i,3);%从未知点推算路线数组中读取推算信息 for jj=1:gcjs; if j=JKH(jj,1) & k=JKH(jj,2) &

37、h=JKH(jj,3) jds=jj; if j<=0 j=abs(j);cjjsjx=YZDZB(j,1);cjjsjy=YZDZB(j,2); else cjjsjx=XYJS(j,1);cjjsjy=XYJS(j,2); end if k<=0 k=abs(k);cjjskx=YZDZB(k,1);cjjsky=YZDZB(k,2); else cjjskx=XYJS(k,1);cjjsky=XYJS(k,2); end dxjk=cjjskx-cjjsjx;dyjk=cjjsky-cjjsjy;%以上为检索出J、k点的端点坐标if abs(dxjk) >= abs(d

38、yjk) afjkt=atan(dyjk/dxjk); if dxjk>=0 & dyjk>=0 afjk=afjkt; end if dxjk<0 afjk=pi+afjkt; end if dxjk>=0 & dyjk<0 afjk=2*pi+afjkt; end end if abs(dxjk) < abs(dyjk) afjkt=atan(dxjk/dyjk); if dyjk>=0 afjk=0.5*pi-afjkt; end if dyjk<0 afjk=1.5*pi-afjkt; end end c=gcjd(1,j

39、ds); chd=(fix(c)+fix(c-fix(c)*100)/60+(c-fix(c)-fix(c-fix(c)*100)/100)*100/36)*pi/180; afjh0=afjk+chd; end j=JSWZD(i,1);h=JSWZD(i,3); for hh=1:gcbs;%此循环是检索j、h对应的边长序号，以便调用该边长观测值 if j=JK(hh,1) & h=JK(hh,2) bcs=hh; end if j=JK(hh,2) & h=JK(hh,1) bcs=hh; end end end XYJS(h,1)=cjjsjx+gcbc(1,bcs)*

40、cos(afjh0);XYJS(h,2)=cjjsjy+gcbc(1,bcs)*sin(afjh0);%计算未知点近似坐标B=zeros(gcjs+gcbs,wzds*2); % 给误差方程系数阵的元素赋0,行数是观测值个数,列数是未知数个数.i=1; while 1 if i=gcjs+1;%观测角的个数是gcjs,因此i=gcjs+1时中断 break; endj=JKH(i,1);k=JKH(i,2);h=JKH(i,3);%读取角度信息 if j<=0 j=wzds+abs(j);%未知点个数是wzds个 end if k<=0 k=wzds+abs(k);%未知点个数是w

41、zds个 end if h<=0 h=wzds+abs(h);%未知点个数是wzds个 end dxjk=XYJS(k,1)-XYJS(j,1);dyjk=XYJS(k,2)-XYJS(j,2); dxjh=XYJS(h,1)-XYJS(j,1);dyjh=XYJS(h,2)-XYJS(j,2); sjk0=dxjk*dxjk+dyjk*dyjk; sjh0=dxjh*dxjh+dyjh*dyjh; if abs(dxjk) >= abs(dyjk) afjkt=atan(dyjk/dxjk); if dxjk>=0 & dyjk>=0 afjk=afjkt;

42、end if dxjk<0 afjk=pi+afjkt; end if dxjk>=0 & dyjk<0 afjk=2*pi+afjkt; end end if abs(dxjk) < abs(dyjk) afjkt=atan(dxjk/dyjk); if dyjk>=0 afjk=0.5*pi-afjkt; end if dyjk<0 afjk=1.5*pi-afjkt; end endif abs(dxjh) >= abs(dyjh) afjht=atan(dyjh/dxjh); if dxjh>=0 & dyjh>=0

43、 afjh=afjht; end if dxjh<0 afjh=pi+afjht; end if dxjh>=0 & dyjh<0 afjh=2*pi+afjht; end end if abs(dxjh) < abs(dyjh) afjht=atan(dxjh/dyjh); if dyjh>=0 afjh=0.5*pi-afjht; end if dyjh<0 afjh=1.5*pi-afjht; end endc=gcjd(1,i); chd=(fix(c)+fix(c-fix(c)*100)/60+(c-fix(c)-fix(c-fix(c)*

44、100)/100)*100/36)*pi/180; fwjhc=afjh-afjk; if fwjhc<0 fwjhc=2*pi+fwjhc; end l(i)=(chd-fwjhc)*206265;%此句为计算角度观测值的自由项，单位是秒j=JKH(i,1);k=JKH(i,2);h=JKH(i,3);if j<=0 & k<=0A=0; end if j<=0 & k>=0 B(i,k*2-1)=+ro*dyjk/sjk0;B(i,k*2)=-ro*dxjk/sjk0; end if k<=0 & j>=0 B(i,j*2-

45、1)=-ro*dyjk/sjk0;B(i,j*2)=+ro*dxjk/sjk0; end if k>=0 & j>=0 B(i,j*2-1)=-ro*dyjk/sjk0;B(i,j*2)=+ro*dxjk/sjk0;B(i,k*2-1)=+ro*dyjk/sjk0;B(i,k*2)=-ro*dxjk/sjk0; end if j<=0 & h<=0 A=0; end if j<=0 & h>=0 B(i,h*2-1)=-ro*dyjh/sjh0;B(i,h*2)=+ro*dxjh/sjh0; end if j>=0 &

46、h<=0 B(i,j*2-1)=B(i,j*2-1)+ro*dyjh/sjh0;B(i,j*2)=B(i,j*2)-ro*dxjh/sjh0; end if j>=0 & h>=0 B(i,j*2-1)=B(i,j*2-1)+ro*dyjh/sjh0;B(i,j*2)=B(i,j*2)-ro*dxjh/sjh0;B(i,h*2-1)=-ro*dyjh/sjh0;B(i,h*2)=+ro*dxjh/sjh0; end i=i+1; end while 2 if i=gcjs+gcbs+1;%观测角加上观测边的总数是gcjs+gcbs,因此i=gcjs+gcbs+1时中断

47、 break; endj=JK(i-gcjs,1);k=JK(i-gcjs,2);%减去观测边的个数,gcjs if j<=0 j=wzds+abs(j);%未知点个数是wzds个 end if k<=0 k=wzds+abs(k);%未知点个数是wzds个 end dxjk=XYJS(k,1)-XYJS(j,1);dyjk=XYJS(k,2)-XYJS(j,2); sjk0=sqrt(dxjk*dxjk+dyjk*dyjk); l(i)=(gcbc(1,i-gcjs)-sjk0)*100;%此句为计算边长观测值的自由项，单位是厘米。 j=JK(i-gcjs,1);k=JK(i-g

48、cjs,2);%减去观测边的个数,gcjs if j<=0 & k<=0 A=0; end if j<=0 & k>=0 B(i,k*2-1)=dxjk/sjk0;B(i,k*2)=dyjk/sjk0; end if k<=0 & j>=0 B(i,j*2-1)=-dxjk/sjk0;B(i,j*2)=-dyjk/sjk0; end if k>=0 & j>=0 B(i,j*2-1)=-dxjk/sjk0;B(i,j*2)=-dyjk/sjk0;B(i,k*2-1)=dxjk/sjk0;B(i,k*2)=dyjk/sjk0; end i=i+1; e