测量平差实验报告3

1、实验三秩亏自由网平差一、实验目的和要求1. 掌握秩亏自由网平差的函数模型及原理;2. 提高编制程序、使用相关软件的能力；3. 熟练使用秩亏自由网准则处理测量数据。二、实验时间及地点三、实验内容设有模拟水准网如图所示，测得观测高差为h】=12.345m h2=3.478mh3=-15.817m各段路线距离相等（等权），取h|°=0h2°= 12.345mh3°= 15.823m1）误差方程v=b*x-1-1 1 00100b=0 -1 1,1=0知：p=010,x0=012.34515.823 t1 0 -160012）组成方程 法方程系数：n二b'*p*b

2、_ 2 -1-1"得：n= - 12- 1-1-12w= w=b'*p*1w =60-63)计算(nn) 和 nm二n (nn) n二 b 怯 pbn =2 1-11 2-1-1 12n*nans =6-33363-3-36pinv(n*n)ans =0.22220.111100.11110.22220000a4)求xxa = (2 0-2)5)平差结果h1=h. o+xt =0.002mv l=-2mm hal= hl+vl= 12.343h2= h2(>+xa2= 12.345m v2=-2mm ha2= h2+vl= 3.476h3= h.o+xa3 =15.82

3、1m v3=-2mm ha3= h3+vl= -15.819a6) x的协因数qaxax=ans =2/901/9-1/92/9-1/9-1/9-1/92/9四、结果与分析此例闭合差w若按经典的条件平差或间接平差也可求得v ,两种方法所得改正数相 同，这是由于最小二乘的准则是相同的，v的值也是相同的唯一的，这是秩亏自由网平差的 一个重要性质。五、总结与体会这次实习的内容是做秩亏自由网平差的计算，这种是比较常见的内容，这样在以后实际应用 中就能够更好的使用，这次实习可以较好的使用软件來进行我们项进行的计算。以后可以进 行较为简单的数据独立处理。六、程序代码（计算清单，可另附纸张，可以打印后粘贴）

4、j lan.ablil edit vjwjindov helpworkspacenyp 1 b h 711】u八-ni】【ii】 二二 rhu【u -二二12 double array24 double array72 dcubl« wray72 double array|si”bjrus clb =» 1=0;0;6command xmndov/006? current directory: | cwtlabgpnwork-11oo-1110-1x3x1x3x33 3 3 3launch pad workspacereadyto get started, select p

5、atlab help* fron the kelp nenu» p=(l 0 0:0 1 0:0 0 1leellt1d q爲 g 46 g c7c«jrrent directory: p2:watla06p1wotk二|flit edit viev w«b vindov kelp» f=r *p»blaunch padworkspacesizebytes» r stjj» if=r *p»1> pinv(n»x)ans =247224727232花s1 p brm rm【n】ni 】【二-rm r

6、mdoublewraydoublearraydoublearraydoublearraydoublearraydoablovrrvydoublearrayclass0.0741 0.0370-0.0370-0.03700.0741-0 0370-0.0370-0.03700 0741ready丿 iatlabpfile edit vi tvjindov help? current directory: |cjmatl/aep讥workworjdrpaceh e71 sze 一 y pianx 1 g binlmniy 二! 1rm【ii】【71】 rm rm rm rm rm rm n -m

7、rm【tn rm 5【5 二 f】rm m rml iclass double d«ubl« double double double double double double doublearrayarryarrayorr»yarrayarraywryarrayarraydoublearrayreadylaunch pad workspace» p=l 0 0；0 10；0 0 1p =1 000 100 01j» m=r* *p«bh =2 1-1-1 2-1-1 -12» 5 *p«1w =606»

8、 w*mans =6-3-3-36-33-36|>|/ xd q7c«jrrent directory: p2:watla06p1wotk二|flit edit viev w«b vindov kelplaunch padworkspacesizebytes g|stack247224727232花s1 p brm rm【n】ni 】【二-rm rmdoublewraydoublearraydoublearraydoublearraydoublearraydoablovrrvydoublearrayclassreadyj iat1.abflit edit vievvi

9、ndov kelpworkspsca>> pinv(n*n)ans =0.0741-0.0370-0.03700.07410.0370 0.0370> pinv(ans)onz =6.0000-3.00)0-3 00006.00003.0000 -3.0000» fl= 8 -3.-3 6m =» invlm)ans =0.22220. 11110.11110.2222-o 0370-0.03700.0741-3 oooo -3 00006 oooo» p:nv(h»k)=0.2222 0. 1111 0：0 1111 0.2222 0

10、:0 0 0 <1艮 ? cerentdirectory： |cwatla06p1 workd2$jzlsizebytesclass“ sistackdoublearraydoubley f 5 x 1 g y hixnh rm rm【111ni 】【i = rm rm rm rm【i!launch pad workspacereadydoubledoubledoubledoubledouble doublearraytrr ayarraywrayarraywrdoubledoublevrxvfarray-36-3-3-36>> pinv(m»x)4bs =0.0741-0.0370-0.0370-0.03700.0741-0.0370-0.0370-0.03700 0741>> pinv(oxis)ws =6.00003.000