大地测量学课程设计

1、应用大地测量课程设计灯湖矿区控制网设计目录一、目的要求及任务范围2二、测区的自然地理条件3三、测区有关测绘资料 3四、测区已有地形图4五、平面坐标系统和高程系统 5六、与线网的建立 6七、高程控制测量 8八、埋标与经费预管级测绘彳2二i11九、工作量综合计算及工作进程计划表12十、控制网相关嫡令平H辑I晅 13十一、水准高程锲中布粉12282s 40十二、上交资料清单 462014-6-5应用大地测量学课程设计目的要求及任务范围1 目的：1 总结和检验大地测量学基础基本知识的学习情况。工程控制网分测图控制网施工控制网变形监测控制网等。通过设计控制网的技术实践，深化已有知识，拓宽新的知识，掌握控

2、制网设计的方法。2将大地测量学基础课程中涉及到的大量的零散的独立的观点和资料，按照设计任务通知书的要求，进行分析归纳综合，完成技术设计任务，达到培养和提高学生的逻辑思维和创造性思维能力的目的。1 技术设计说明书是对工程设计进行解释与说明的书面材料，是一种技术性文件。设计者通过技术设计，编写技术设计书，是进行科技写作的锻炼，培养科技 写作的实际能力。2 目的要求：1设计的项目和内容应该齐全并符合大纲和规范的要求。设计分为几个步骤：学习领会技术任务书图上设计构网做出精度估算制定观测方案绘制控制网图编写技术设计说明书。2设计的论点应该正确，明确表达设计者的主张意见和看法。论据力求做到真实充分新颖。公

3、式推导正确，推理符合逻辑。3认真编写技术设计说明书。在使用专业词语布局谋篇及至行款格式等方面，都要加强训练。3 任务要求：1阅读领会任务通知书2熟悉测区地理环境及原有测绘成果等情况；对原有控制测量成果进行分析和评定，确定其利用程度。3根据任务书要求，选择平面坐标系统和高程系统，拟定起始数据的配置和处理方案。4根据任务书要求和测区面积，确定首级控制网的等级和布网层次，确定作业的技术依据。5在地形图上进行首级平面控制网的选点设计，作出精度估算；平面控制网一般应设计两个方案，以便用人机对话方式进行优化，选出最佳方案。6在地形图上进行高程控制网的选线设计，确定水准网或测距三角高程网布设方案，并进行待定

4、点高程精度估算；注意平面网点的高程联测方案。7在首级控制网的基础上，再作出加密控制网布设方案。8确定控制点标石规格。9制定观测纲要；编制建网进度计划表。10编制经费预算表。11 编写技术设计说明书。12 .绘制所设计的平面和高程控制网图（A4图纸，着墨成图）。:、测区的自然地理条件1地理概况本测区中心位置为东经117° 13',北纬34° 12'。测区地面高程为+30+244m位于XX省XX县境内，属于华东丘陵地带，南部及中部为山地。其上林木繁茂，不利通视，其余为平原，密布村庄和厂矿，建筑群遍地，影响观测。灯湖位于测区东部，湖面开阔，中有大堤。2交通情况测区

5、内有铁路支线通过，公路四通八达，村间大道可行汽车，交通方便。3气候情况测区地处华东近海地区，气候宜人。全年平均降雨量为1040mm雨量集中于6、7、 8 三个月份。全年平均气温+15，夏季气温较高，一般为30左右，冬季有雪，但不寒冷，最低温度为-5C,冻土线深度为0.1m。全年平均风力为23级,夏季略受台风影响。宜于野外作业时间为311 月份，年平均作业时间利用率为21 天 / 月。4居民及居民点测区内地少人多，人口稠密，多为汉族，少数为回族。测量作业所需人力、物力、财料及食宿均可就地解决。三、测区有关测绘资料1三角网成果及其精度测区内及附近有国家二等网点三个：玉山、太山、广具山。系 1959

6、年由XX测绘队施测，作业所依据的规范为一、二、三、四等三角测量细则（1958 年）。三点标石保存完好。坐标系统为1954 年北京坐标系，三度分带，中央子午线为117°。2高程网成果及其精度国家n等水准路线由西向东横穿测区北部。 根据二等水准路线图，本测区内及附 近应有三个二等水准点，依点之记只找到其中两个：n-2, n-4,标石保存完好。该二等水准路线系国家测绘队于1979 年施测。施测精度及埋石质量均符合规范要求。高程为1956年黄海高程系统。3地形图资料1984年XX航测队航测本地区1: 50000地形图，成图质量良好，可供技术设计图上选点之用。二等三角点、水准点坐标和高程：点

7、名 等级x （ m）y （ m）H （ m）II-4二35.000n-2二46.000四、测区已有地形图通过对已有资料的分析及论证，我们决定平面坐标系统采用1954 年北京坐标系，而高程系统则采用1985国家高程基准。GPS3观测与数据处理采用基线解算、无约束平差、约 束平差，经坐标转换最后得到 GPSK的1954年北京坐标系坐标。GPS点的水准高程通过已知水准高程公共点采用多项式拟合法求得。在图上设计时我 们则暂时取图纸的西南角为坐标原点，定出测区内两个假定控制点的坐标，并以此为依据 进行水准及高程控制网的设计，及精度估计。待对该测区进行实测和进行工程建设时，只 需确定原点在国家80西安坐标

8、系中坐标，即可延展出所需点坐标六、导线网的建立1、1选点(1)相邻点之间应通视良好，其视线距障碍物的距离不宜小于1.5m(2)远离大功率无线电发射源(如电视台、微波站等)，其距离不小于400米，远离高 压输电线200米以上；(3)为避免多路径效应，点位附近不应有强烈干扰卫星信号接收的物体，并尽量避开 大面积水域；(4)交通方便，有利于其它测量手段扩展和联测；基础稳固，便于点的保存。1 、2埋石与点之记三角点标石应用混凝土灌制，或用相同规格的花岗石、青石等坚硬石料代替1 、3?使用仪器；拓普康GTS-225全站仪设备1套，脚架3个，反射棱镜2套，喷漆两瓶，钢钉若干。DS3自动安平水准仪1台，脚架

9、1个，黑红双面尺1对1、4?线路布设路线图x M犷JK垃犬捕t E蝇”红色为首级e级GPS)g,绿色为水准网(四等)，蓝色为一级附和导线。2、布网依据的原则(1)全球定位系统(GPS测量规范，GB/T 18314-2009。(2)全球定位系统城市测量技术规程，CJJ73-97。(3)工程测量规范，GB50026-2007(4)城市测量规范，CJJ8-99。(5)国家一、二等水准 测量规范，GB12897-91。(6)国家三、四等水准 测量规范。GB12898-91(7)1:500 1:1 000 1:2 000 地形图图式GB/T20257.1-2007 ;3、外业观测A、控制点布设要求(1)

10、点位选在质地坚硬、稳固可靠、便于保存的地方，视野应相对开阔，便于加密、 扩展和寻找，并按照要求绘制点之记；(2)相邻点之间应通视良好，其视线距障碍物的距离，三、四等不宜小于1.5m；四等以下宜保证便于观测，以不受旁折光的影响为原则；(3)当采用电磁波测距时，相邻点之间视线应避开烟囱、散热塔、散热池等发热体 及强电磁场；(4)相邻两点之间的视线倾角不宜太大；(5)充分利用旧有控制点。B 、观测的基本技术要求(1)导线测量技术要求(2)光 电测距 的技术 要求：C、观测 过程等级导线长平均边长每边测距中误差测角中水平角测回数导线全长闭合差导线全长相对闭合差四等9km1.5Km± 18mm

11、±2.5J265； 1/35000观测要求往返测较差测回数一测回读数较差测回较差一测回读数次数往、返15mm210mm15mm3(1)、观测路线：玉山-广具山-太山-玉山(2)人员安排：两人负责棱镜与配套脚架，一人为观测员，一人为记录员，一人为协 调员。(3)观测方法:在组织组员对测区进行了勘察情况下选定了控制点，在确定已知点坐标和正北方向后，按照三联脚架法对导线进行观测，具体步骤如下：1 ）将全站仪安置整平对中在 BM2点，棱镜分别安置在后视点BM侪口前视点A的基座中， 进行导线测量。2 ）采用两个测回观测数据，分别置为 0° 0' 0和90° 0

12、9; 0，读取水平角B ,和距离 S。3）迁站时，BM2点和A点上的脚架和基座不移动，将全站仪安置在A点的基座上，BM2 点上则安置棱镜，再将 BM1点的仪器迁到B点，随后再如前一站进行观测，直到导线测量 完毕。4、内业计算使用控制测量优化与平差软件进行数据计算，最后得出网形与精度统计、导线闭合差、控制点成果表、方向平差成果表、边长成果表、高差平差成果表、点位误差、点间误差。七、高程控制测量1、水准观测（一）、水准测量的主要技术要求，应符合有关的规定。（二）、水准测量所使用的仪器及水准尺，应符合下列规定：（1）、水准仪视准轴与水准管轴的夹角， DS1型不应超过15 ；，DS3型不应超过20 ；

13、（ 2） 、 水准尺上的米间隔平均长与名义长之差，对于因瓦水准尺，不应超过0.15mm，对于双面水准尺，不应超过0.5mm；（ 3）、水准点应选在土质坚硬、便于长期保存和使用方便的地点。点位应便于寻找、保存和引测。2、水准控制的一般规定（一）测区的高程系统，采用1985 国家高程基准。在已有高程控制网的地区进行测量时，可沿用原高程系统；当小测区联测有困难时，亦可采用假定高程系统。（二） 高程控制测量，可采用水准测量和电磁波测距三角高程测量。高程控制测量等级的划分，应依次为二、三、四等。(三)首级网应布设成环形网。当加密时，宣布设成附合路线或结点网。(四) 高程控制点间的距离，一般地区应为13k

14、m, 了业厂区、城镇建筑区宜小 于lkm。但一个测区及周围至少应有3个高程控制点。应布设成环形网。当加密时，宜 布设成附合路线或结点网。3、布网形式和要求(1)本测区以国家高等级水准点作为高程控制起算点，布设三等水准网，作为高程控制，以满足测区高程控制发展的需要。(2)等外水准、测距高程导线，自三等水准联测点起发展不得超过2次。(3)三等水准观测采用DS3型号的水准仪，中丝法读数，各测段测站为偶数。作业 前须对水准仪和水准尺进行检校。(4)控制网与水准网的连测由相应水准点处采用三角测量的方法传递。4 、实测方法及过程(一)仪器：三等水准测量中使用 DS3自动安平水准仪1台，脚架1个，黑红双面尺

15、1对，记录板1个。(二)三等水准测量技术要求前后视距较不前后视距累计差K+黑一红黑红面高差之差视线最长距离往返较差、闭合差<5m10m3mm5mm75m12/L mm注：表中L为闭合导线长度，以km为单位(三)、仪器检校水准仪i角检验(1)在比较平坦的地面上选取距离为 30米的A、B两点，然后在BA延长线上选取一 点1,使1到A的距离等于30米，最后在AB延长线上选取2点，使2点至ij B点的距离也 等于20.6米。(2)在A、B两点放上尺垫并在尺垫上竖立水准尺，在 1点架设水准仪分别读取在 A、 B两点所竖立的水准尺上的中丝读数 a1、b1 (为了提高a、b1的精度，可在A、B两点水准

16、 尺上读中丝读数3至4次，然后取其平均值)；(3)将水准仪搬至2点，以同样的方法分别获取A、B两点所竖立的水准尺上的中丝读数 a1、b1;（ 4）根据下式即可计算出水准轴与视准轴之间的夹角i 角：如下图所示：（ 四 ） 、作业方法（ 1）观测路线：韩愈路- 民主路- 神州路- 迎宾大道（ 2）人员安排：两人为跑尺员，一人为观测员，一人为记录员，一人为计算员。（ 3）观测方法：在选取了合适的路线中点之后进行观测，具体步骤如下：1 ）安置整平仪器，照准后尺黑面，依次读取上、下丝及中丝读数。2 ）转动仪器，照准前尺黑面，依次读取上、下丝及中丝读数。3 ）前尺变红面朝向仪器，读取中丝读数。4）后尺变红

17、面，仪器照准后尺红面，读取中丝读数。（五）注意事项：1 ）前后视距相等（ 在限差内） 。2）从后视转为前视时，望远镜不得从新调焦。3）读数前应消除视差。4 ）当用水准仪瞄准、读数时，水准尺必须立直。5）每站观测完毕后，必须及时进行计算，检核满足限差要求后才能搬站。注：当水准线路较长时，可设置转站，水准尺立于尺垫之上，并保证最后测站总数为偶数。八、埋标与经费预算标石埋设要求：（ 1）、盐碱地区埋设混凝土标石，须加涂沥青，以防腐蚀。（ 2）、在泥土松软、地下水位较高的地区或沼泽地区埋设标石时，除应尽量选择好埋石地点以外，应在盘石下边浇灌混凝土底层。（ 3） 、 埋石时， 须使各层标石的标志中心严格

18、在同一铅垂线上，其偏差不大于3mm。并用钢尺量取各层标石面间的垂直距离，填记于点之记的标石断面图中，结果取至厘米。装备：运送人员车辆一辆、运送材料卡车一辆；仪器：选才¥ 3台GPStl频接收机（配带脚架、卷尺），2台D3水准仪（配带脚架、标尺、尺垫）、全站仪一台（配带脚架、棱镜、气温气压计、卷尺）；器材：草帽、对讲机、标石、混凝土、重力计等经费预算：1 .工资：工程师：3000元/人？月*1人技术员：2000元/人？月*1人测工：1200元/人？月*3人临时工：800元/人？月*2人2 .仪器费用：单频接收机：每台15000元*3台全站仪每台：20000元*1台DS3水准仪：每台20

19、00元*1标石：每座100元*73座3 .生活补贴费：每人/每天100元*30天*7人4 .交通住宿费：每人/每天100元*30天*7人总费用=126000元九、工作量综合计算及工作进程计划表外业计算：用全站仪进行导线网观测时，应及时计算出每测回的2c值、2c互差值（一级导线小于13），多测回方向角（方向角互差小于 9）；水准仪观测时应及时计算出前后视距差（小于1米）、累计视距差（小于3米）；内业计算：概算和平差（见前述） 工作进程计划表：时间任务备注1天踏勘找点（已知点）查看标后情况7天选位和造标埋石GPSK、水准点、导线点需待标石凝固符 合要求后方可观 测，导线点需要通 视1天GPSR观测

20、每时段观测时间 不少于45分钟1天GP制数据解算及结果精度分析确保精度达标7天水准仪检校 水准网施测水准路线需要往 返观测，且路线要 至少经过4个GPS点3天水准网图程计算和GPSSJtlb程拟合 （包括错误测站重测时间和休息时问）确保精度达标3天全站仪检校 导线网施测及时计算发现改 正错误1天导线网计算及错误测站复测确保精度2天休息1天GPSR、水准网、导线网全面检查查看标后情况2天:整理资料及结果、上交图纸1天单位检验验收十、控制网相关参数与平差结果第一种方案【GPS网】的边长由MAPINFO估算得到,mapinfo测量数据起点终点距离km-113.904-143.904212.11025

21、2.110323.112363.376-232.321-263.270621.952652.532413.061452.638843.482853.270-381.108-352.954-363.904-392.638962.2169101.477763.8517-22.4791062.4791072.374平均边长2.771750KM基线向量的观测中误差：所以该网的观测值中误差为:29.4662( mn) = c 0将未知点坐标设为未知数，列(计算过程由 Matlab完成：):一 =a2b2 = - 102 (10 10-6 2771750)2 =从每个观测时段选取两条独立基线，最好各时段选

22、取的独立基线不重复，出用未知数表示观测值的方程组，从而可以得出未知数的系数矩阵，即B矩阵B=zeros(48,20);JK=1,99;4,99;1,2;2,5;3,2;3,6;3,99;6,99;6,2;6,5;4,1;4,5;8,4;8,5;8,99;5,99;6,99 ;9,99;9,6;9,10;7,6;7,99;10,6;10,7;i=1;while 1if i=25; break ;endj=JK(i,1);k=JK(i,2);if j=99;B(i-1)*2+1,(k-1)*2+1)=-1;B(i-1)*2+2,(k-1)*2+2)=-1;elseB(i-1)*2+1,(j-1)*

23、2+1)=1;B(i-1)*2+2,(j-1)*2+2)=1;end if k=99;a=1;elsewhile j=99B(i-1)*2+1,(k-1)*2+1)=-1;B(i-1)*2+2,(k-1)*2+2)=-1; break ;endendi=i+1;endi=1;s=0;D=blkdiag(3.904,3.904,2.110,2.110,3.112,3.376,2.321,3.270,1.952,2.532,3.061,2.638,3.482,3.270,1.108,2.954,3.90426382216,1.477,3.851,2.479,2.479237);while 1if

24、i=25; break ;ends=s+D(i,i);i=i+1;endq=100+(10+10*s/24)A2;C=zeros(48,48);i=1;while 1if i=25; break ;endC(i-1)*2+1,(i-1)*2+1)=100+(10+10*D(i,i)A2;C(i-1)*2+2,(i-1)*2+2)=100+(10+10*D(i,i)A2;i=i+1;endQ=C/q;P=inv(Q);BTPB=B'*P*B;QXX=inv(BTPB);晦阵Q:Qxx:由该表知1号点和2号点为最弱点，二 GPS =2.085cm < 5 cm【一级附和导线】符合导线

25、：600. 26204450. 49481478. 5028784531.1449316492. 4368942435. 2679316459.862018506.1750143409. 8531646384. G0077B2375. 2485302395. 3859258488. 9500987453. 6739984448. 9715423493. 6350001407. 0597704445. 6624879376. 2052594443. 555928442. 2854091453. 1100547389. 0355256起始点玉山，终点太山中间布设21个点，相邻点之间距离平均为450

26、.49481 m,且该导线经过主井石红。详细情况见图B:QXX:16号点的点位中误差最大。Qx=2.5640最弱点的点位中误差为：=2.66233cm 根据测量规范及测量任务书，此误差符合限差要求，小于 5cm, matlab代码如下： clear;yzds=2;wzds=21; gcjs=21;gcbs=22;cjzwc=5;cbzwc=1.45;cbblwc=0.0;JKH=1,-1,2;2,1,3;3,2,4;4,3,5;5,3,6;6,5,77,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,16;1

27、6,15,17;17,16,18;18,17,19;19,18,20;20,19, 21;21,20,-2;JK=-1,1;1,2;2,3;3,4;4,5;5,6;6,7;7,8;8,9;9,10;10,11;11,12;12,13;13,14;14,15;15 ,16;16,17;17,18;18,19;19,20;20,21;21,-2;XYJS=3793081.17875,6222.93767;3792631.83417,6164.64960;3792222.69855,6332.654 33;3791809.45075,6493.79958;3791435.20284,6532.125

28、80;3791005.64895,6650.85239; 3790618.03974,6775.18550;3790142.75376,6641.84246;3789694.40540,6665.46884;378 9240.96620,6650.85239;3788762.50062,6751.55912;3788368.50591,6718.41810;378800 9.65563,6828.13476;3787643.21735,6944.93473;3787248.32295,7054.65139;3786757.2 7726,7177.48358;3786582.73781,7602

29、.93406;3786790.96467,7985.16496;3786822.7430 0,8476.57491;3786692.45447,8991.49275;3786713.02634,9469.55456;3793400,6000; 3786800,10250;B=zeros(gcjs+gcbs,wzds*2);i=1;while 1if i=gcjs+1;break ;endj=JKH(i,1);k=JKH(i,2);h=JKH(i,3);if j<=0j=wzds+abs(j);endif k<=0k=wzds+abs(k);endif h<=0h=wzds+a

30、bs(h);enddxjk=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;j=JKH(i,1);k=JKH(i,2);h=JKH(i,3);if j<=0 & k<=0A=0;endif j<=0 & k>=0B(i,k*2-1)=+ro*dyjk/sjk0;B(i,k*2)=-ro*dxjk/sjk0;endif

31、 k<=0 & j>=0B(i,j*2-1)=-ro*dyjk/sjk0;B(i,j*2)=+ro*dxjk/sjk0;endif k>=0 & j>=0B(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;endif j<=0 & h<=0A=0;endif j<=0 & h>=0B(i,h*2-1)=-ro*dyjh/sjh0;B(i,h*2)=+ro*dxjh/sjh0;

32、endif j>=0 & h<=0B(i,j*2-1)=B(i,j*2-1)+ro*dyjh/sjh0;B(i,j*2)=B(i,j*2)-ro*dxjh/sjh0; endif j>=0 & h>=0B(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;endi=i+1;endwhile 2if i=gcjs+gcbs+1;break ;endj=JK(i-gcjs,1);k=J

33、K(i-gcjs,2);if j<=0j=wzds+abs(j);endif k<=0k=wzds+abs(k);enddxjk=XYJS(k,1)-XYJS(j,1);dyjk=XYJS(k,2)-XYJS(j,2);sjk0=sqrt(dxjk*dxjk+dyjk*dyjk);gcbc(i-gcjs)=sjk0;j=JK(i-gcjs,1);k=JK(i-gcjs,2);if j<=0 & k<=0A=0;endif j<=0 & k>=0B(i,k*2-1)=dxjk/sjk0;B(i,k*2)=dyjk/sjk0;endif k<

34、;=0 & j>=0B(i,j*2-1)=-dxjk/sjk0;B(i,j*2)=-dyjk/sjk0;endif k>=0 & j>=0B(i,j*2-1)=-dxjk/sjk0;B(i,j*2)=-dyjk/sjk0;B(i,k*2-1)=dxjk/sjk0;B(i,k*2)=dyjk/s jk0;endi=i+1;endfor i=1:gcjs;P(i,i)=cjzwc*cjzwc/(cjzwc*cjzwc);endfor i=gcjs+1:gcjs+gcbs;cbfc=(cbzwc+cbblwc*gcbc(i-gcjs)/1000)*(cbzwc+cb

35、blwc*gcbc(i-gcjs)/1000); P(i,i)=cjzwc*cjzwc/(cbfc);endBTPB=B'*P*B;QXX=inv(BTPB);while 3if j=wzds+1;break ;endsgmPQ)=cjzwc*sqrt(QXX(2*j-1,2*j-1)+QXX(2*j,2*j);j=j+1；end第二种方案：【GP轲】3460.11874-2837. 483792845. 0140852374. 7770613103. 95909530Q3. 3020812694. 3209373588. 2694412482. 7078743171. 6801732

36、006. 6492193399. 4279172933. 7609022559. 1619832472. 5293632530. 816042420. 5024283111. 9240673933. 5794433102. 3937942515. 468578280L 4S85712225. 53395423G8. 974129KM30.0854( mr)l =；二0将未知点坐标设为未知数，列(计算过程由 Matlab完成：)平均边长基线向量的观测中误差： 所以该网的观测值中误差为：-a2b2 = ；,；102( 10 1 0-62 8 3 7 4 84)2 二从每个观测时段选取两条独立基线，

37、最好各时段选取的独立基线不重复，出用未知数表示观测值的方程组，从而可以得出未知数的系数矩阵，即B矩阵B=zeros(46,18);JK=1,99;1,2;2,3;4,99;1,4;1,6;2,6;2,7;3,7;3,99;7,99;6,7;4,6;4,5;5,99;6,5;6,99;7,99;8,99;7,8;7,9;8,9;9,99;i=1;while 1if i=24;break ;endj=JK(i,1);k=JK(i,2);if j=99;B(i-1)*2+1,(k-1)*2+1)=-1;B(i-1)*2+2,(k-1)*2+2)=-1;elseB(i-1)*2+1,(j-1)*2+1

38、)=1;B(i-1)*2+2,(j-1)*2+2)=1; endif k=99;a=1;elsewhile j=99B(i-1)*2+1,(k-1)*2+1)=-1;B(i-1)*2+2,(k-1)*2+2)=-1; break ; end endi=i+1;endi=1;s=0;D=blkdiag(3460.11874,2845.014085,2374.777061,3103.659095,3063.362681,2684.3208 37,3588.269441,2482.707874,3171.686173,2006.649219,3399.427917,2983.760902,255 9

39、.161983,2472.529368,2536.81664,2420.502428,3111.924067,3983.579443,3102.39379 4,2515.468578,2801.488971,2225.533854,2368.974129);while 1if i=23; break ; ends=s+D(i,i);i=i+1;endq=100+(10+10*s/23)A2;C=zeros(46,46);i=1;while 1if i=24; break ;endC(i-1)*2+1,(i-1)*2+1)=100+(10+10*D(i,i)A2;C(i-1)*2+2,(i-1)

40、*2+2)=100+(10+10*D(i,i)A2;i=i+1;endQ=C/q;P=inv(Q);BTPB=B'*P*B;QXX=inv(BTPB);BQ:QXX由该表知1号点和2号点为最弱点，二 GPSI= 1.9902 cm < 5 cm根据测量规范及测量任务书，此误差符合限差要求，小于5cmi【一级附和导线】符合导线：起始点玉山，终点太山，中间他设16个点，相邻点之间距离平均为F 347.55891486. 38785435. 6936781458. 182330G440. 2375429460 4037217525. 0064879431. 0251532464. 33

41、43621561. 0021245529. 4506245565. 3465329451. 41956452也 9761193631. 1733585516. 0870496431. 530027过主井石红。i¥细情况见图：487. 1658802486.38785 m,且该导线经B:QXX:11号点的点位中误差最大。Qx=1.3860最弱点的点位中误差为：二一 =2.1195cm根据测量规范及测量任务书，此误差符合限差要求，小于5cmiMatlab计算程序：clear;yzds=2;wzds=16;gcjs=16;gcbs=17;cjzwc=5;cbzwc=1.49;cbblwc=

42、0.0;JKH=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,16;16,15,-2;JK=-1,1;1,2;2,3;3,4;4,5;5,6;6,7;7,8;8,9;9,10;10,11;11,12;12,13;13,14;14,15;15 ,16;16,-2;XYJS=3786970.60373,10706.31871;3787228.00838,11052.67465;3787194.26909,11567. 65845;37

43、87117.51580,12194.12589;3787084.33125,12720.05486;3787459.65822,12970. 86498;3787852.55604,13377.37570;3787742.63690,13895.28902;3788258.74888,14115. 17717;3788496.17672,14514.21451;3788438.62734,14941.38143;3788879.99306,15225. 68775;3789129.12649,15619.98177;3789555.41456,15729.92585;3789986.81992

44、,15884. 25519;3790381.09525,16069.67479;3786800,10250;3790725,15900;B=zeros(gcjs+gcbs,wzds*2);i=1;while 1if i=gcjs+1;break ;endj=JKH(i,1);k=JKH(i,2);h=JKH(i,3);if j<=0j=wzds+abs(j);endif k<=0k=wzds+abs(k);endif h<=0h=wzds+abs(h);enddxjk=XYJS(k,1)-XYJS(j,1);dyjk=XYJS(k,2)-XYJS(j,2);dxjh=XYJS

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

46、0;endif k>=0 & j>=0B(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;endif j<=0 & h<=0A=0;endif j<=0 & h>=0B(i,h*2-1)=-ro*dyjh/sjh0;B(i,h*2)=+ro*dxjh/sjh0;endif j>=0 & h<=0B(i,j*2-1)=B(i,j*2-1)+ro*dyjh/sjh0;B(i,j*

47、2)=B(i,j*2)-ro*dxjh/sjh0; endif j>=0 & h>=0B(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;endi=i+1;endwhile 2if i=gcjs+gcbs+1;break ;endj=JK(i-gcjs,1);k=JK(i-gcjs,2);if j<=0j=wzds+abs(j);endif k<=0k=wzds+abs(k);enddx

48、jk=XYJS(k,1)-XYJS(j,1);dyjk=XYJS(k,2)-XYJS(j,2);sjk0=sqrt(dxjk*dxjk+dyjk*dyjk);gcbc(i-gcjs)=sjk0;j=JK(i-gcjs,1);k=JK(i-gcjs,2);if j<=0 & k<=0A=0;endif j<=0 & k>=0B(i,k*2-1)=dxjk/sjk0;B(i,k*2)=dyjk/sjk0;endif k<=0 & j>=0B(i,j*2-1)=-dxjk/sjk0;B(i,j*2)=-dyjk/sjk0;endif k&g

49、t;=0 & j>=0B(i,j*2-1)=-dxjk/sjk0;B(i,j*2)=-dyjk/sjk0;B(i,k*2-1)=dxjk/sjk0;B(i,k*2)=dyjk/s jk0;endi=i+1;endfor i=1:gcjs;P(i,i)=cjzwc*cjzwc/(cjzwc*cjzwc);endfor i=gcjs+1:gcjs+gcbs;cbfc=(cbzwc+cbblwc*gcbc(i-gcjs)/1000)*(cbzwc+cbblwc*gcbc(i-gcjs)/1000);P(i,i尸cjzwc*cjzwc/(cbfc);endBTPB=B'*P*B;

50、QXX=inv(BTPB);while 3if j=wzds+1;break ;endsgmPQ)=cjzwc*sqrt(QXX(2*j-1,2*j-1)+QXX(2*j,2*j);j=j+1;end总结：方案二的精度更高，所以选择方案二为最终方案。十一、水准高程控制网布设方案1.水准网采取四等符合水准网平均12km一个水准点。从已知水准点II2到II4布设符合水准路线，水准点大多布 设于公路道路旁，详细情况见图。1.网中最弱点精度估算及结果附和水准路线：以km定为单位权，M h = M /L各水准测段权值如下(取对角线)：1219.266521539. 5478612827.375181700. 97921952.364351522. 474761110. 554162372. 4222010. 8081963. 685162440. 280041865,888252622.18387662.952021412. 266231761.926241761.92844607,16237828.59