平差课程设计.doc

课程编号： 课程性质：必修误差理论与测量平差基础课程设计报告学院： 测绘学院 专业： 测绘工程 班级： 2011级学号： 姓名： 2013-2-25 至 2013-3-8（一）图1为一水准网，A、B、C为已知高程点，为9个待定高程点，第一次观测了图中115条水准路线的高差，各水准路线的观测高差、距离及已知点高程均列于表1，试求（1）1公里高差平差值中误差；（2）各待定点高程平差值及中误差；如果准备加测5段高差1620（图中用虚线表示），（3）试估算平差后各待定点的精度；（4）精度最弱的点发生改变没有；（5）哪一点精度提高得最多。 图1 表1 水准网观测数据线路号观测高差(m)距离(km)线路号观测高差(m)距离(km)11.49570.4112.04160.821.78920.5121.53550.735.20620.7131.74500.842.08530.9142.74340.453.41750.4151.37850.561.20801.0160.676.03730.6170.684.33070.6181.893.65100.6190.8102.51200.3201.2已知点高程：、数学模型： 间接平差： 函数模型： V = B x - l , l=L - L0; 基础方程：V = B x - l BT P V = 0 解向量： x=NBB-1 W 随机模型：D=02 Q Qxx=NBB-1计算过程：(1).列误差方程 易得到矩阵B和矩阵l：B=000100000100000000-100000000010-1000000-1000000000100000000-100000000001000000001-1000000001-100000000-10100010000-10001000000000100-1000000001-1 l= （单位：0.1mm）0-29030-17000200406(2).组成法方程取1公里的观测高差为单位权观测，即按 P=1/S定权，得观测值的权阵：P=2 1/200000000000000020000000000000001 3/70000000000000001 1/90000000000000002 1/200000000000000010000000000000001 2/30000000000000001 2/30000000000000001 2/30000000000000003 1/30000000000000001 1/40000000000000001 3/70000000000000001 1/40000000000000002 1/20000000000000002组成法方程： BTPBx - BTPl=0令：NBB=BTPB W=BTPl 得： x =NBB-1 W解法方程 x =-16.9170.93821-6.3750.049181.4821.072-1.13312.7434-1.6798改正数为： v =（单位：0.1mm）0.0491812.08316.917-2.111-0.938210.6256.3751.4820.410070.20503-0.54681.729-3.9508-1.2614-1.5767解：(1). 1公里高差平差值中误差： = 1.2422 mm(2). 各待定点高程平差值及中误差： 平差值：X = X0 + x 144.6687148.2427144.8917146.0689148.7982145.1062142.3737146.1807144.3605 因为Qxx=NBB-1 ，所以Dxx=Qxx DXX=0.450060000000000.4562400.094040.0298780.0402210.0453920.0429020.059182000.5786500000000.0940400.305630.0971050.130720.147520.139430.1923400.02987800.0971050.662530.571390.525820.547760.4042900.04022100.130720.571391.16090.99270.559320.5442400.04539200.147520.525820.99271.22610.56510.6142100.04290200.139430.547760.559320.56510.905220.5805200.05918200.192340.404290.544240.614210.580520.80081则各点中误差为： 0.67090.67540.76070.55280.81391.07741.10730.95140.8949(3).最弱点为P7 ，其精度为1.1073(4).加测后的B=000100000100000000-100000000010-1000000-1000000000100000000-100000000001000000001-1000000001-100000000-10100010000-10001000000000100-1000000001-1-11000000001-100000000-1100000000-110000000000-110加测后的权P= diag(5/2,2,10/7,10/9,5/2,1,5/3,5/3,5/3,10/3,5/4,10/7,5/4,5/2,2,5/3,5/3,5/9,5/4,5/6)所以加测后的Qxx=(NBB)-1 Dxx=QxxDXX=0.332530.090760.0330810.0188310.0103340.0116510.012310.0120960.014210.090760.277470.101130.0575680.0315920.035620.0376340.0369790.0434410.0330810.101130.355670.0489570.0268660.0302920.0320040.0314480.0369430.0188310.0575680.0489570.258120.141650.159710.168740.16580.194780.0103340.0315920.0268660.141650.49550.440630.413190.422110.334090.0116510.035620.0302920.159710.440630.986430.796410.510980.479980.012310.0376340.0320040.168740.413190.796410.988010.555410.552920.0120960.0369790.0314480.16580.422110.510980.555410.772430.529220.014210.0434410.0369430.194780.334090.479980.552920.529220.76325精度为： 0.57670.52680.59640.50800.70390.99320.99400.87890.8736(5).精度最弱的点没有改变，仍是P7。(6).比较可知，P3点的精度提高的最多。结果检核：h1= HP4 - HB =146.0640m - 144.5684m=1.4956m h2= HP1 - HB =146.3587m - 144.5684m=1.7903mh3= HA - HP1 = 151.5664 - 146.3587 = 5.2077m h4= HP2 - HP4 = 148.1590m - 146.0640m = 2.0850mh5= HA - HP2 = 151.5664m - 148.1590m = 3.4175mh6= HP3 - HC =145.5286m-144.3194m = 1.2092mh7= HA - HP3 = 151.5664m - 145.5286m = 6.0378mh8= HP5 - HC = 148.6501m - 144.3194m = 4.3307mh9= HP5 - HP6 = 148.6501m - 144.9991m = 3.6510mh10= HP6 - HP7 = 142.9991m - 142.4869m = 2.5122mh11= HP9 - HP7 = 144.5283m - 142.4869m = 2.0414mh12= HP4 - HP9 = 146.0640m - 144.5283m = 1.5357mh13= HP4 - HC = 146.0640m - 144.3194m = 1.7446mh14= HP5 - HP8 = 148.6501m - 145.9067m = 2.7434mh15= HP8 - HP9 = 145.9067m-144.5283m = 1.3784m经检核，由每个高程点平差后的高程得到的各点间的高差与已知高差平差后的结果相同。(二)（题目详见课程设计辅助系统）数学模型(1) 函数模型：测角网函数模型：测边网函数模型：（2）随机模型：（3）平差准则:VPV=min计算过程： 根据题意可知：n=46 t=30 r=16 现在选取待定点坐标平差值作为参数:X=X1 Y1 X2 Y2 X3 Y3 X4 Y4 X5 Y5 X6 Y6 X7 Y7 X8 Y8 X10 Y10 X11 Y11 X12 Y12 X13 Y13 X14 Y14 X15 Y15 T 设单位权中误差为=6 , 则角度观测的权为，导线边的权为 PSi = 02/s2 = 36/(k*1000*s)2，其k为边精度1/2000。利用Excel表格计算已知方位角和各个待定点的近似坐标如下：点号观测角坐标方位角()距离(m)坐标分量坐标0x(m)y(m)X/mY/m17126.6353573.431181.148230.5015 （17 5 6）300.148527246.783877849.091-19.35169801-45.11582943137.33289.4266 ( 5 6 7 ）182.246222249.030134.896-12.48849147-32.58478781117.978302244.31017067 （ 6 7 8）189.94822258.978322245.167-8.635044177-44.33389111105.4898105211.72538288 （ 7 8 9）229.510638308.488961169.70243.38000465-54.5577125796.85476635167.39149169 （8 9 14）288.50711156.9960722256.65230.8581476447.51024972140.234771112.833779114（9 14 15）180.10897257.1050444448.98226.6021499641.12860247171.0929186160.344028815（14 15 17）242.520027119.625072233.447-16.533598129.07476471197.6950686201.472631317（15 17 5）187.005611126.630683300181.1614705230.5473961457.1050555648.98226.6021419841.12860764171.0929437160.344035315（14 15 16）182.19211159.2971666751.65526.3742907944.41436491197.6950857201.472642916(15 16 2)171.97077751.2679444553.00633.1647511341.3489458224.0693765245.88700782 (16 2 3)291.549000162.816944429.068-27.770572268.587429205257.2341276287.23595363 (2 3 4)203.371277186.188222252.942-52.63351456-5.706882617229.4635553295.82338284 (3 4 5)174.754694180.942916739.543-39.53764535-0.650729796176.8300408290.11650025 (4 5 17)305.683472306.6263889137.2923954289.46577048308.488972269.70243.38001521-54.5577041796.85479002167.39148389 （8 9 10）213.0911111341.580083342.7340.54078105-13.50177662140.2348052112.833779610(9 10 11)212.165111113.7451944230.91630.030627657.34577831180.775586399.3320030111(10 11 12)183.897166717.6423610940.09838.2120648712.15268292210.8062139106.677781312(11 12 13)204.361944442.0043055352.58539.0756264435.1891694249.034118.77613(12 13 0)230.360027892.3643333156.115-2.31495084556.06722953288.1096264153.96516940 (13 0 1 )177.282527889.6468610945.5750.2808968545.57413435285.7946756210.03239891 ( 0 1 2 )222.7866111132.433472242.798-28.8772519631.58754697286.0755724255.60653332 ( 1 2 3 )210.3978333162.8313055287.1940803然后由测角网坐标平差的误差方程和测边的误差方程计算出系数矩阵B，根据 x=(BTPB)-1 BTPl 和 v=Bx-l计算出x v。 xv-7.5578-15.46336.827-16.092-13.063-5.788330.767-7.1603-3.47990.0104221.09395.6127-5.47760.354991.1964-10.778-1.071511.3930.769860.018061-1.719110.3290.0195080.2327-7.8670.56991-12.877-9.2728-13.11410.55-32.462-18.957-19.612-10.434-27.469-18.29110.475-17.426-40.792-16.85913.37-1.0674-45.344-5.04493.06592.955448.835-0.70401-9.27465.1665-9.1863-5.6353-1.3859-6.09364.5967-28.3670.155271.93896.4971-4.3346-1.05860.6597414.24320.228-7.951832.7561.7303-0.4776634.958-11.55820.96215.8582.4208-6.4896-4.6810.42 (1) 由 求得 13.62 (2)各个待定点的坐标平差值根据=X+v =Y+v求得，详见表格。 各个待定点的坐标中误差由 ，最后由求 得各个点的中误差。坐标答题报表 坐标编号 X平差值 Y平差值 X改正数 Y改正数 极大值方向 极大值 极小值0 285.786 210.069 -8 37 66 4.45 3.061 286.062 255.637 -13 31 66 3.49 2.342 257.231 287.236 -3 1 25 2.81 2.393 229.459 295.824 -5 1 5 2.48 2.054 176.829 290.117 -1 1 2 3.06 0.935 137.33 289.426 0 0 0 0 06 117.976 244.31 -2 0 56 2.1 17 105.476 211.712 -8 -13 68 3.26 1.718 96.841 167.359 -13 -32 90 2.53 2.49 140.214 112.813 -20 -27 32 2.9 2.2410 180.785 99.291 10 -41 14 3.36 1.9111 210.819 106.632 13 -45 17 3.03 1.2712 249.034 118.776 0 0 0 0 013 288.112 154.014 3 49 46 3.86 1.6714 171.083 160.335 -9 -9 63 3.25 2.0715 197.694 201.477 -1 5 299 2.18 0.5716 224.169 245.893 0 6 60 3.15 2.1317 181.148 230.501 0 0 0 0 0各点的点位中误差为：010.80111.1427.3338.0546.3864.7177.3388.5097.33107.70116.53138.38147.70154.30167.57（3）各观测值的平差值=观测值+改正数 边长答题报表 边长编号 边长观测值 边长平差值 边长改正数 0 - 1 45.575 45.569 -6 1 - 2 42.798 42.77 -28 2 - 3 29.068 29.07 2 3 - 4 52.942 52.938 -4 4 - 5 39.543 39.542 -1 5 - 6 49.091 49.092 1 6 - 7 34.896 34.91 14 7 - 8 45.167 45.187 20 8 - 9 69.702 69.694 -8 9 - 10 42.73 42.763 33 10 - 11 30.916 30.918 2 11 - 12 40.098 40.098 0 12 - 13 52.585 52.62 35 13 - 0 56.115 56.103 -12 9 - 14 56.652 56.673 21 14 - 15 48.982 48.998 16 15 - 16 51.655 51.657 2 16 - 2 53.006 53 -6 15 - 17 33.447 33.442 -5 17 - 5 73.431 73.431 0 测角答题报表 测角编号 测角观测值 测角平差值 测角改正数0 1 2 2224711.8 2224656.2 -15.51 0 13 182432.9 1824246.7 -16.11 2 3 2102352.2 2102346.4 -5.81 2 16 2785050.5 2785043.2 -7.22 3 4 2032216.6 2032216.6 03 2 16 68273.6 68279.1 5.63 4 5 1744516.9 1744517.2 0.44 5 6 2454956.7 2454945.9 -10.84 5 17 305410.5 3054111.8 11.45 6 7 1821446.4 1821446.4 06 5 17 59515.3 595115.3 106 7 8 1895653.6 1895653.5 07 8 9 2293038.3 2293039.2 18 9 10 2130528 2130518.9 -98 9 14 2883025.6 2883036.5 119 10 11 2120954.4 2120935.4 -1910 9 14 752449.9 752439.9 -1010 11 12 1835349.8 1835331.7 -1811 12 13 2042143 2042126 -1712 13 0 2302136.1 2302119.1 -179 14 15 1800632.3 1800631.2 -114 15 16 1821131.6 1821126.6 -514 15 17 2423112.1 2423115.1 315 16 2 1715814.8 1715813.8 -116 15 17 601937.7 601942.7 515 17 5 1870020.2 1870014.1 -6 （4）由 线性化的权函数式为 所得的系数阵为F，数值详见表格。 由以及计算得到各个边的中误差。 最弱边为最大的边长，由计算所得的结果，最弱边为1213 两点间的边长相对中误差为 1213/s1213= 1/1308 （5） 结合E F知最弱的误差椭圆为0号点。 实习心得上学期我们学习了误差理论与测量平差基础这门重要的专业基础课程。当时由于课程安排紧凑（每个星期有6个学时），课后也没有很多时间来复习和预习，加上这门课本身的难度，导致很多东西都没有及时消化，许多概念也是模棱两可，最后的期末考试成绩自然也不是很理想。但是通过这次误差理论与测量平差实习，我找到了解决实际问题的办法，锻炼了我分析问题、解决实际问题的能力，特别是计算能力。可以说这次的课程设计是很成功的。参加了实习动员大会后，我对这次实习一点底都没有，因为这次实习要编写程序。虽然C语言学的还可以，但也仅仅限于书上的基本理论知识，对于实际编程解决问题能力还是有所欠缺。后来老师说如果编程能力差点的同学可以用matlab或者excel来解决问题，这时我才增加了些许的信心。经过了一个多月的寒假，我对这门课已经忘得差不多了。拿到课程设计任务书后，看到虽然只有两题，但是我还是不敢怠慢。我做的第一件事就是把课本重新温习了一遍，尤其是间接平差部分，因为这次的课程设计主要是利用间接平差模型来解决问题。这样我不但把遗忘的理论的知识又重新记住了，而且还把很多上学期没有学的章节弄懂了。 第一题是水准网的平差，所以我选择间接平差来做。第一步是求未知点的近似高程，然后是列各观测值的误差方差，系数矩阵B和l可以很快地列出来，再一步步解出矩阵BTPl和NBB以及NBB的逆矩阵，求出改正数v和x，最后对所测的结果进行平差，平差后再进行检核。通过检核，得出平差无误的结论。因为第一题是书上例题的变型，所以没有遇到什么难处，矩阵的计算用matlab也很好的解决了。 第二题才是这次实习的重点吧，第一题与它比起来完全只是热身。它的计算量大的惊人，都是46*30的大矩阵。第二题是导线网的平差，它包括了测边网和测角网两种平差，难度可想而知。而当时上课的时候因为太复杂老师也没有将其作为重点来讲，所以我不得不再次回归书本，将那几节好好的看一遍，把书上例题的表格每一项是怎么来的都研究清楚了才开始做题。大体的思路和第一题差