灯湖矿区控制网设计.doc

灯湖矿区控制网设计一、工程控制网设计任务书灯湖矿区控制网设计任务通知书XX测绘队：灯湖矿区位于XX省XX煤田东北域。本矿区已列入国家经济建设计划，准备进行重点开发。根据建设计划，首先需施测80平方公里的1：2000地形图。图一为该矿区1：50000地形图，其中虚线范围为1：2000测图区。为了加快建井速度，需在竖井间进行对向掘进，最大贯通距离为3km，矿井最大单翼长度5km，井深一般不超过800m。为满足矿山设计、建井、生产三阶段测绘各种比例尺地形图、井巷贯通以及工业场地施工测量的需要，应在全矿区建立统一的具有足够精度密度的平面控制网与高程控制网。要求你队在任务通知书下达后一周内完成控制网技术设计，一月内完成建网任务。XX省煤炭管理局 年 月 日二、测区自然地理条件及已有测绘成果资料（一）自然地理条件1地理概况本测区中心位置为东经11713，北纬3412。测区地面高程为+30+244m。位于XX省XX县境内，属于华东丘陵地带，南部及中部为山地。其上林木繁茂，不利通视，其余为平原，密布村庄和厂矿，建筑群遍地，影响观测。灯湖位于测区东部，湖面开阔，中有大堤。2交通情况测区内有铁路支线通过，公路四通八达，村间大道可行汽车，交通方便。3气候情况测区地处华东近海地区，气候宜人。全年平均降雨量为1040mm，雨量集中于6、7、8三个月份。全年平均气温+15，夏季气温较高，一般为30左右，冬季有雪，但不寒冷，最低温度为-5，冻土线深度为0.1m。全年平均风力为23级，夏季略受台风影响。宜于野外作业时间为311月份，年平均作业时间利用率为21天/月。4居民及居民点测区内地少人多，人口稠密，多为汉族，少数为回族。测量作业所需人力、物力、财料及食宿均可就地解决。（二）已有测绘成果资料1三角网成果及其精度测区内及附近有国家二等网点三个：玉山、太山、广具山。系1959年由XX测绘队施测，作业所依据的规范为一、二、三、四等三角测量细则（1958年）。三点标石保存完好。坐标系统为1954年北京坐标系，三度分带，中央子午线为117。2高程网成果及其精度国家等水准路线由西向东横穿测区北部。根据二等水准路线图，本测区内及附近应有三个二等水准点，依点之记只找到其中两个：-2，-4，标石保存完好。该二等水准路线系国家测绘队于1979年施测。施测精度及埋石质量均符合规范要求。高程为1956年黄海高程系统。3地形图资料1984年XX航测队航测本地区1：50000地形图，成图质量良好，可供技术设计图上选点之用。二等三角点、水准点坐标和高程： 点 名 等级 x（m） y（m） H（m）玉 山 二 3793400 20506000 195.4广具山 二 3790725 20515900 135.0太 山 二 3786800 20510250 244.0-4 二 35.000-2 二 46.000 三、测区已有地形图（四）坐标系统的选择及处理方法的论证、起始数据的配置和处理工程平面控制网的坐标系统以与国家大地网的坐标系统取得一致为最理想，此时，工程控制网便能够成为国家大地网的组成部分。但是要使工程控制网完全采用国家统一坐标系统，必须具备下列条件之一： (1)工程区域位于高斯投影3带中央子午线附近，同时测区平均高程面接近平均海水面； (2)虽然不能符合上述条件，但根据测区的y坐标和高程H计算的相对长度综合变形小于1/40000。事实上很多工矿测区不能满足上列条件之一。此时应该结合测区已有国家高级点的数目、位置、精度等情况作下列选择： 1. 如果测区内国家大地控制点的精度符合现行规范对起算数据的精度要求，且具有必要的或多余的(四个或四个以上)起算数据时，可以选择测区局部坐标系统(或称测区独立坐标系统)，选择方案可为： (1) 选择抵偿高程面作为投影面，按高斯正形投影3带计算平面直角坐标； (2) 保持国家统一的椭球面作投影面，选择抵偿投影带，按高斯正形投影计算平面直角坐标； (3) 选择测区平均高程面作为投影面，通过测区中心的子午线作为中央子午线，按高斯正形投影计算平面直角坐标。2. 测区内虽有多个(两个及两个以上)国家大地控制点，但其精度不详或低于现行规范对起算数据的精度要求。此时可选择其中一个标石完好，位置适中的国家点作原点，利用其坐标(x,y)和一条边的坐标方位角T作为起算数据，另外重新测定一条或一条以上的起算边长，将测距结果化算至测区平均高程面。按这种方法建立的坐标系统虽属矿区局部坐标系统，但其原点的坐标和方位角仍保持了国家统一坐标系统。这种方法已为许多工矿区和城市所采用，它是弥补国家高级网某些缺陷的好方法。该矿区坐标系统选择过程：Ym=(y1+y2)/2;Hm=(H1+H2)/2经过计算可得：Ym=-11.5km； Hm=137m所以： =(0.00123*(-11.406)2-15.7*137)/100000=1/503461/40000所以直接采用3度带国家高斯投影。（五）平面控制网布设方案1）首级GPS网（E级）设计：首级控制网的等级与要求精度及测区面积有关。一个点能够控制多大的面积，或者某一面积范围内应有多少个点，这是点的密度问题。在城市测量规范和工程测量规范中已有各等级的边长系列，实际上已明确了点的密度。 考虑到工程地区大都采用GPS网作为首级控制网，且GPS网大多布设成类似三角网的形式，。结合本次任务测区内实际情况，选择平均边长为2.5km。按四等GPS网的主要技术要求布设，在布设过程中，也应注意一点：虽然GPS测量时不要求相邻点之间通视，布网方式也非常灵活，但考虑到加密导线测量的需要，每个GPS点应与相邻的12个点通视。综合考虑以上各种因素，测区面积约为75平方千米，每个控制点控制的面积约为5.6平方千米(按公式S=0.87D2粗略计算)，所以在测区内布设E级GPS控制点11个，与3个已知点组成由14个点构成的GPS平面控制网。设计草图见下页：1，设计简要说明：四台GPS同时观测，共8个时段，可获独立基线24条，具体可见下页图。2，精度评定结果：运用matlab软件编程估算其精度，具体代码见matlab7.0程序附件gps2.m，下表（一）是各控制点点位精度和误差椭圆参数。（A表示椭圆方向角，E,F分别表示极大值和极小值误差，Mx,My分别表示X方向和Y方向的误差，M表示点位误差）；表（二）是相对点位精度和误差椭圆参数。（T表示椭圆方向角，E,F分别表示极大值和极小值误差，s是基线长，Ms表示基线误差，Ms/S表示相对边长误差，单位为万分之一）协因素阵：由精度估算公式：可得：点位中误差：表（一） 点号 A(dms) E(mm) F(mm) Mx (mm) My (mm) M(mm)相对点位中误差：表（二） 基线号 T(dms) E(mm) F(mm) S(km) Ms (mm) Ms/S (1/10000)-1-1 1-22-3 -3-3-3-76-75-64-5-1-41-52-62-73-7-3-98-97-86-86-114-11-2-11-2-108-109-10-2-83，由表（一）可知最弱点中误差为E09号点，点位精度29.7mm(小于30mm),符合要求，相对点位精度最大方向误差值为24.6mm（9-10）,相对边长精度最大值为-3-7号基线(-3号点为广具山)，其相对精度为1/192382，均符合规范精度要求。2）一级加密网设计： 根据一级导线的规范要求，结合测区实际情况，布设成共8个点（包括两个E级已知点E004,E011）的附合导线网，平均边长约为0.58km，网型布设情况见前面首级GPS网型设计草图。已知点和布设的E级GPS控制点的图上粗略量测坐标如下表： （1954年西安坐标系、中央子午线经度117度）序号点号转换后X坐标转换后Y坐标1-1(玉山)3793400205060002-2(E04)3790590205065903F013792880205063304F023792400205065605F033791910205064506F043791500205066707F05379100020506760精度评定结果：运用matlab软件编程估算其精度，具体计算过程及代码实现见matlab7.0程序附件daoxian2.m，下表（一）是各控制点点位精度和误差椭圆参数。（A表示椭圆方向角，E,F分别表示极大值和极小值误差，Mx,My分别表示X方向和Y方向的误差，M表示点位误差）；表（二）是相对点位精度和误差椭圆参数。（T表示椭圆方向角，E,F分别表示极大值和极小值误差，s是基线长，Ms表示基线误差，Ms/S表示相对边长误差，单位为万分之一）协因素阵：由精度估算公式：可得：点位中误差：表（一） 点号 A(dms) E(mm) F(mm) Mx (mm) My (mm) M(mm)相对点位中误差：表（二） (-1代表E011点位, -2代表E004点位) 基线号 T(dms) E(mm) F(mm) S(km) Ms (mm) Ms/S (1/10000)-1-11-22-33-44-5-2-5由表（一）可知最弱点中误差为F03号点，点位精度约为13.9mm(小于20mm),符合要求，相对点位精度最大方向误差值约为9.2mm（1-2）,相对边长精度最大值为-1-1（即玉山-F01）号基线，其相对精度为1/157673，均符合规范精度要求。(六)高程控制网设计： 根据两个个已知的二等高程控制点II-2，II-4，根据工程实际要求，布设四等水准控制网。1，两个已知控制点，,水准高程具体如下表：点名H/m备注II-246.000高斯投影L0=1171954国家高程基准II-435.0002，水准网型布设草图见下页：点位精度评定： 1，根据布网，找出高程精度弱点和结点，用边长定权，估算出其点位高程精度，所依据依据公式：2，编写程序计算，具体见附录matlab7.0程序附件shuizhun2.m，精度评定结果：协因素阵：由高差中误差公式可得： 点号 点位精度3，由估算结果可知，最弱点为F19号点，点位精度约为19.1mm,小于限制误差20mm，符合规范精度要求。（七）技术依据及作业方法1. 测量规范及技术规程（1）、全球定位系统（GPS）测量规范，GB/T18314-2001。（2）、全球定位系统城市测量技术规范，CJJ73-97。（3）、工程测量规范，GB12897-91。（4）、城市测量规范，CJJ8-99。(5 )、国家三、四等水准测量规范，GB12898-91。(6 )、1：500、1：1000、1：2000地形图图式，GB/T7929-1995。2.标石图（1）A.2.5 四等导线/三角形网/GPS平面控制点标石埋设规格应符合图A.2.5规定。图A.2.5 四等导线/三角形网/GPS平面控制点标石埋设图（单位：mm）注：1土面；2素土；3冻土线； 4贫混凝土（2）A.2.6 一、二级导线/五等GPS平面控制点标石埋设规格应符合图A.2.6规定。图A.2.6 一、二级导线/五等GPS平面控制点标石埋设图（单位：mm）1 土面；2素土；3贫混凝土（3）A.3.3 四等水准点标石埋设规格应符合图A.3.3的规定。图A.3.3 四等水准点标石埋设图（单位：mm）1地面；2素土；3冻土线；4贫混凝土3，仪器的选择及检验项目要求（1）仪器：GPS单频接收机，全站仪，DS3水准仪（2）检验项目： 水准仪的光学测微器效用正确性的检验和分划值的确定；（3）视准轴与水准轴相互关系的检验；（4）水准尺检验：水准尺分划面弯曲差的检验；（5）水准标尺每千米分划间隔真长的测定；（6）一对水准标尺零点差及基，辅分划读数差常数的测定；4.观测方法及各项限差（1）E级GPS观测方法及各项限差要求：1）观测组应严格按调度表规定的时间进行作业，以保证同步观测同一卫星组。当情况有变化需修改调度计划时，应经作业队负责人同意，观测组不得擅自更改计划。2）接收机电源电缆和天线应连接无误，接收机预置状态应正确，然后方可启动接收机进行观测。3）各观测时段的前后各量取天线高一次，两次量高之差不大于3mm。取平均值作为最后天线高，记录在手簿。若互差超限，应查明原因，提出处理意见记入手簿备注栏中。天线高是指观测时天线平均相位中心至测站中心标志面的高度，分为上、下两段：上段是指相位中心至天线底面的高度，这是常数hc，由厂家给出；下段是从天线底面至测站中心标志面的高度，由观测员在现场采用倾斜测量方法直接量取。具体方法是：从三脚架三个空档（互成120）测量天线底盘下表面至测站中心标志面的距离，互差应小于3mm，取平均值为L，天线底盘半径为R，再利用厂家提供的hc，按天线高 求出。4）接收机开始记录数据后，作业人员可使用专用功能键选择菜单，查看测站信息、接收卫星数、卫星号、各通道信噪比、实时定位结果及存贮介质记录情况等。5）仪器工作正常后，作业员及时（每隔15min）逐项填写测量手簿中各项内容。6）一个时段观测过程中不得进行以下操作：关闭接收机以重新启动；进行自测试（发现故障除外）；改变卫星截止高度角；改变数据采样间隔；改变天线位置；按动关闭文件和删除文件等功能。7）观测员在作业期间不得擅自离开测站，并防止仪器受震动和被移动，防止人和其它物体靠近天线，遮挡卫星信号。8）接收机在观测过程中不应在接收机近旁使用对讲机和手机等通讯设备；雷雨过境时应关机停测，并卸下天线以防雷击。9）观测中应保证接收机工作正常，数据记录正确，每日观测结束后，应及时将数据下载到计算机硬、软盘上，确保观测数据不丢失。表1GPS网的主要技术要求级别平均距离(km)(mm)(110-6)最弱边相对中误差E级0.2510201/45000表2E级GPS测量作业的基本技术要求级别卫星截止高度角()有效观测卫星数平均重复设站数时段长度(min)数据采样间隔(s)PDOP值E级1541.6401510表3同步环坐标分量及环线全长相对闭合差的规定（110-6）等级限差类型E 级坐标分量相对闭合差6.0环线全长相对闭合差10.0（2）一级导线观测方法及各项限差要求：1，左角采用前-后-后-前，右角采用后-前-前-后，2，水平角采用测回法，竖直角采用中丝法，三丝法3，导线施测采用三联脚架全园观测法施测，水平角观测的技术要求按工程测量规范GB50026-932.3.12.3.10执行。4，测角导线转折角有左角和右角之分。在导线前进方向左侧的水平角称为左角，右侧的水平角称为右角。闭合导线一般测其内角，在公路测量中，附和导线一般测右角，注意全线应统一。各等级的导线测角要求，应满足规范。5，测边光电测距，光电测距导线边采用单向或往返观测，导线边长均观测2测回，每测回4次读数，一测回内读数较差应小于5mm,单程测回间较差应小于10mm。表一，光电测距导线的水平角技术要求等级仪器等级测回数左+右-360之差 /上下半测回之差/测角中误差/方位角闭合差/一级I级210105表二，光电测距导线的竖直角技术要求等级仪器等级测回数竖直角互差/指标差互差/一级I级21010表三，光电测距导线的测距技术要求 等级仪器等级导线长度/km测回数测距中误差/mm相对误差导线全长相对闭合差一级I级44151/30 0001/14 000（3）四等水准观测方法及各项限差要求：5，概算内容及平差方法：（1）导线测量概算得主要内容：作业测量观测数据的检查，将观测值划算到标石中心，将观测值划算到椭球面上和高斯投影平面上。（2）GPS控制网的概算：数据处理及基线解算，观测结果的外业的检核。（3）水准测量的概算得主要内容：观测高差各项改正数（包括水准标尺每千米长度改正数，水准面不平行改正数）的计算和水准点概略高程表的编制。导线控制网，水准网采用间接平差方法，GPS控制网用GPS基线向量网三维无约束平差。6，提供成果资料（1）GPS网概算结果及平差结果；（2）一级导线网概算结果及平差结果；（3）一级导线网概算结果及平差结果。.（八）工作量综合计算及工作进程计划表1，查找资料，考察测区：5天；2，控制网优化设计：一周7天；3，准备相关材料及租借仪器和组织人员：一周7天；4，实地施测：8天；5，数据处理及成果提交：3天；共计：5+7+7+8+3=30天，由于天气等因素影响，再加上弹性时间5天，合计35天。（九）需要设备及经费预算1，仪器：GPS单频接收机四台，全站仪，DS3水准仪，光电测距仪等。2，粗略经费预算：（1）、设计费和施测费：10万元 （2）、材料费用：2万 （3）、人员工资：9万 （4）、仪器租用费：2万 （5）、交通等其他费用：2万 工程预计费用合计：25万元 （十）上交的资料清单1，控制测量技术设计书；2，外业观测原始资料；3，内业计算成果表及控制网图。附录matlab程序代码：1, gps2.mB=zeros(48,22);JK=1,99;1,2;2,3;3,99;7,99;6,7;5,6;4,5;4,99;1,5;2,6;2,7;3,7;9,99;8,9;7,8;6,8;6,11;4,11;11,99;10,99;8,10;9,10;8,99; %基线向量i=1; while 1 %给系数阵B阵赋值 if i=23; break; endj=JK(i,1);k=JK(i,2); B(i-1)*2+1,(j-1)*2+1)=1; B(i-1)*2+2,(j-1)*2+2)=1; if k=99 a=1; else B(i-1)*2+1,(k-1)*2+1)=-1; B(i-1)*2+2,(k-1)*2+2)=-1; end i=i+1; end S=0.64*4.6,4.6,4.0,4.0,4.2,4.2,3.8,3.8,5.0,5.0,3.7,3.7,2.2,2.2,4.2,4.2,4.3,4.3,3.8,3.8,3.5,3.5,4.0,4.0,3.9,3.9,4.1,4.1,3.6,3.6,3.0,3.0,3.6,3.6,3.9,3.9,4.1,4.1,4.0,4.0,4.0,4.0,2.6,2.6,3.5,3.5,3.7,3.7; %0.64是从图上距离转换到实际距离的参数 MO=27; C=(S.*10).2.+100)./(MO2);%以2.5千米为单位权中误差(其基线向量中误差为27mm) Q=diag(C); P=inv(Q);%权阵 BTPB=B*P*B; QXX=inv(BTPB);%协因素阵 D=zeros(11,6);%点位精度阵 i=1; while 1 if i=12 disp(D); break; end H=sqrt(QXX(2*(i-1)+1,2*(i-1)+1)-QXX(2*(i-1)+2,2*(i-1)+2)2+4*QXX(2*(i-1)+1,2*(i-1)+2)2); A=90.0;%误差椭圆的最大值方位角 E=MO*sqrt(0.5*(QXX(2*(i-1)+1,2*(i-1)+1)+QXX(2*(i-1)+2,2*(i-1)+2)+H); F=MO*sqrt(0.5*(QXX(2*(i-1)+1,2*(i-1)+1)+QXX(2*(i-1)+2,2*(i-1)+2)-H); Mx=MO*sqrt(QXX(2*(i-1)+1,2*(i-1)+1); My=MO*sqrt(QXX(2*(i-1)+2,2*(i-1)+2); M=MO*sqrt(QXX(2*(i-1)+1,2*(i-1)+1)+QXX(2*(i-1)+2,2*(i-1)+2);%点位误差 D(i,:)=A,E,F,Mx,My,M;%误差椭圆及点位误差阵 i=i+1; end DX=ones(24,6); i=1; while 1 if i=25 disp(DX); break; end j=JK(i,1);k=JK(i,2); if k=99 Qax=QXX(2*(j-1)+1,2*(j-1)+1); Qay=QXX(2*(j-1)+2,2*(j-1)+2); Qaxy=QXX(2*(j-1)+1,2*(j-1)+2); else Qax=QXX(2*(j-1)+1,2*(j-1)+1)+QXX(2*(k-1)+1,2*(k-1)+1)-2*QXX(2*(j-1)+1,2*(k-1)+1); Qay=QXX(2*(j-1)+2,2*(j-1)+2)+QXX(2*(k-1)+2,2*(k-1)+2)-2*QXX(2*(j-1)+2,2*(k-1)+2); Qaxy=QXX(2*(j-1)+1,2*(j-1)+2)+QXX(2*(k-1)+1,2*(k-1)+2)-QXX(2*(j-1)+1,2*(k-1)+2)-QXX(2*(k-1)+1,2*(j-1)+2); end H=sqrt(Qax-Qay)2+4*Qaxy2); E=MO*sqrt(0.5*(Qax+Qay+H); F=MO*sqrt(0.5*(Qax+Qay-H); T=A; %由于Qaxy=0; s=S(1,2*(i-1)+1);%边长 Ms=sqrt(cos(T)2*Qax+(sin(T)2*Qay+sin(2*T)*Qaxy);%边长中误差 f=(100*Ms)/s;%相对边长中误差（单位为万分之一） DX(i,:)=T,E,F,s,Ms,f;%相对误差椭圆阵 i=i+1; end 2, daoxian2.mB=zeros(13,10); % 给误差方程系数阵的元素赋0,行数是观测值个数,列数是未知数个数.JKH=-2,-1,1;-1,1,2;1,2,3;2,3,4;3,4,5;4,5,-2;5,-2,-1;JK=1,-1;1,2;2,3;3,4;4,5;5,-2;XYJS=3792880,20506330;3792400,20506560;3791910,20506450;3791500,20506670;3791000,20506760;3793400,20506000;3790590,20506590;ro=2062.65;i=1; while 1 if i=8;%观测角的个数是7,因此i=9时中断 break; endj=JKH(i,1);k=JKH(i,2);h=JKH(i,3); if j=0 j=5+abs(j);%未知点个数是6个 end if k=0 k=5+abs(k);%未知点个数是6个 end if h=0 h=5+abs(h);%未知点个数是6个 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; j=JKH(i,1);k=JKH(i,2);h=JKH(i,3); if j=0 & k=0 A=0; end if j=0 B(i,k*2-1)=-ro*dyjk/sjk0;B(i,k*2)=ro*dxjk/sjk0; end if k=0 B(i,j*2-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 B(i,h*2-1)=ro*dyjh/sjh0;B(i,h*2)=-ro*dxjh/sjh0; end if j=0 & h=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=14;%观测角加上观测边的总数是13,因此i=14时中断 break; endj=JK(i-7,1);k=JK(i-7,2);%减去观测角的个数7 if j=0 j=5+abs(j);%未知点个数是5个 end if k=0 k=5+abs(k);%未知点个数是5个 end dxjk=XYJS(k,1)-XYJS(j,1);dyjk=XYJS(k,2)-XYJS(j,2); sjk0=sqrt(dxjk*dxjk+dyjk*dyjk); j=JK(i-7,1);k=JK(i-7,2); if j=0 & k=0 A=0; end if j=0 B(i,k*2-1)=dxjk/sjk0;B(i,k*2)=dyjk/sjk0; end if k=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; end K=1,1,1,1,1,1,1,0.25,0.25,0.25,0.25,0.25,0.25;%粗略定权，平均边长约为0.5km P=diag(K);%观测角和观测边的权阵 BTPB=B*P*B; QXX=inv(BTPB);%点位协因素阵 D=zeros(5,7);%点位精度阵 i=1; MO=5;%单位权中误差为5mm while 1 if i=6 disp(D); break; end H=sqrt(QXX(2*(i-1)+1,2*(i-1)+1)-QXX(2*(i-1)+2,2*(i-1)+2)2+4*QXX(2*(i-1)+1,2*(i-1)+2)2); A=0.5*atan(2*QXX(2*(i-1)+1,2*(i-1)+2)/(QXX(2*(i-1)+1,2*(i-1)+1)-QXX(2*(i-1)+2,2*(i-1)+2)+1.57;%误差椭圆方向角0- E=MO*sqrt(0.5*(QXX(2*(i-1)+1,2*(i-1)+1)+QXX(2*(i-1)+2,2*(i-1)+2)+H); F=MO*sqrt(0.5*(QXX(2*(i-1)+1,2*(i-1)+1)+QXX(2*(i-1)+2,2*(i-1)+2)-H); Mx=MO*sqrt(QXX(2*(i-1)+1,2*(i-1)+1); My=MO*sqrt(QXX(2*(i-1)+2,2*(i-1)+2); M=MO*sqrt(QXX(2*(i-1)+1,2*(i-1)+1)+QXX(2*(i-1)+2,2*(i-1)+2);%点位误差 D(i,:)=i,A,E,F,Mx,My,M;%误差椭圆及点位误差阵 i=i+1; end DX=ones(6,6); i=1; S=0.54,0.55,0.54,0.57,0.60,0.59,0.65;%边长 while 1 if i=7 disp(DX); break; end j=JK(i,1);k=JK(i,2); % if k=-1|k=-2|k=-3 if k=0 k=6+abs(k);%未知