第五章条件平差

第五章条件平差5.1条件平差原理间接平差－－对参数平差条件平差－－对观测值平差例：为测定三角形形状，分别测量了其三内角的值L1、L2、L3。因有多余观测及观测误差的存在，观测值间产生矛盾。此题中观测值个数：n=3必要观测数：t=2多余观测数：r=n-t=1条件方程代入上式＋＋平差值条件式：，消除矛盾。使平差值和为最或是值（平差值），角小二乘原则，计算出各条件平差目的：遵循最＋＋产生闭合差：－＋＋由于：0ˆ180ˆˆˆ180180180321321321321WVVVVLLLLLLLLWLLLiii条件方程矩阵形式：一、条件方程的个数及选择条件方程个数=多余观测数=n-t=r条件方程的选择：条件方程间应线性无关。即条件方程中的系数矩阵A阵应为行满秩阵。即。平差值条件方程：A0为已知的常数向量。1110rrnnrWVArARnr0ˆ0ALA•条件方程一般形式：rbannnnrrnnbbnnaann212121210221102211022110000其中：－可得矩阵形式：，，，二、法方程及改正数方程。的原则作用于条件方程将minPVVTkQAkAPVAkPVVkkkkTTTTTrbar11022,,得对新函数求导：条件方程联系数式中WAVkPVVTT2－组成新函数:---改正数方程得基础方程：由上式得法方程：得观测值的平差值为：2011WAVkAPVT。数）式，可得观测值改正代入（其中，设或VkWNkAAPNWNkWkAAPTT100111VLL＋ˆ向量函数求导规律dxdFGdxdGFdxFGddxGFdxgfmgggGfffFTTTTiiTmTm有的函数，均为、维向量函数，且均为设),,(),,(2121注意：当N为满秩方阵时，才有唯一存在,法方程才有唯一解。所以，r个条件方程必须线性独立。1N5.2精度评定一、单位权中误差NkkkQAPkQAPVVkNvppvvPVVPVVrPVVtnPVVTTTTTniiiTTTT求：和用）直接法：（的计算：)2(1ˆ120kWkPQAVPVVkWTTTT求：、用)3(AQNQAQAALNQAkQAVAQNQAQQANQALANQAIVLLTVVTTTLLTT1011ˆˆ011ˆ二、观测值平差值及观测值改正数的权逆阵MATLAB平差应用简介1、启动与退出双击MATLAB图标，进入MATLAB命令窗口。建立新文件：File－New－M-file进入编辑窗。法1、直接编程，再用编辑窗Debug菜单中的Run运行，然后到命令窗中看结果；或法2、在编辑窗中编制函数，再在命令窗或另一个文件中调用。建立文件：function[x,Qxx]=ABC(A,L,d,P);文件名：ABC.mx,Qxx：要求输出内容A,L,d,P—输入内容分号：该句内容不显示2、矩阵的输入A＝[123;456;789]同行的元素可用空格或逗号隔开，式中分号也可换成击回车键。可得：3、矩阵转置A＝A’4、矩阵求逆M＝inv(N)----M为N的逆阵987654321＝A5、矩阵的＋、－、×用符号＋、－、*表示，运算时注意矩阵的阶数。6、特殊矩阵产生A=eye(n)----产生n维单位阵A=zeros(n,m)-----产生n×m维0矩阵7、B=det(A)A的行列式，A必须是方阵B=rank(A)A的秩B=trace(A)A的迹，A阵主元之和8、输出格式formatshout短格式，显示5位。如：1.4142formatlong长格式，15位。如：1.41421356237310三、平差值函数的权逆阵VLLLffVfffVfVfVfVfFLLLLLfFiiLLiiTTTTnnnin~ˆ,ˆˆ~ˆˆˆˆˆˆ22211121其中：展开，得权函数式：的真值点在，，设有平差值的函数：fQffAQNQAQfAQfNQAfQffQQQQAQNQAQQQQffQQQQffQLLTTTTTTFFVVVVTVVLLVVVVTTFFˆˆ11ˆˆ1ˆˆ－－－求观测值的平差值函数的权逆阵过程可记忆为：fQffQfQFVLdLffLdfLdfLdfLdfFLLLLLfFLLTLdLdTFFLLiiTnniniiˆˆˆˆˆˆˆ221121ˆˆ,ˆˆˆˆˆˆ~ˆˆˆˆˆ的权逆阵为则其中：展开，得权函数式：的真值点在，，设有平差值的函数：补充：控制测量概念一、控制测量任务1。为测绘地形图建立控制网2。建立施工控制网3。建立变形监测网二、控制网的形式（按观测值分类）1。高程控制网－－水准网，三角高程网2。平面控制网－－三角网，导线网，GPS网测方向网测角网测角网三角网测边网测边长测边＋测方向边角网（导线网）测边＋测角三、三角网的布设－－从高级到低级逐级布设四、三角网平差的方法1。严密平差----遵守VTPV=min原则；2。近似平差5.3测角网条件平差独立网（经典自由网）－-－只有必要起算数据d。非独立网（附合网）－-－已知条件超过必要起算数据。（必要起算数据个数d－-－水准网：d=1，测角网：d=4，测边、边角网：d=3）一、独立测角网按角度坐标平差r=n-2p=n-tp：网中待定点数1、单三角形V1+V2+V3-W=0-W=(L1+L2+L3-180)－-－图形条件2、中点多边形测角网n=9,t=4,r=5个：边条件极条件个圆周角条件：个图形条件：1)(135rBACD123456789极条件1ˆsinˆsinˆsinˆsinˆsinˆsin852741LLLLLLDCDBDADADCDB＝以“D”点为极点，得平差值条件方程：线性化后的极条件：741852885522774411sinsinsinsinsinsin10cotcotcotcotcotcotLLLLLLWWVLVLVLVLVLVL式中闭合差为：•CD边的相对中误差：LdfSLdLLdLLdLLdLSLdLLdLLdLLdLLLLLABdSLLLLSTCDCDCDCDˆ)ˆcotˆcotˆcotˆcot)ˆcotˆcotˆcotˆcotsinsinsinsinˆsinˆsinˆsinˆsinABˆ887733228877332273827382＋－－（＝＋－－（＝全微分后得：＝得测边相对中误差为：3、大地四边形测角网n=8,t=4,r=4fQffDfSLdfSdSLLTLLTCDSTCDCDCDˆˆ20ˆˆ2ˆˆˆ13极条件：图形条件：fQfSLLTCDSCDˆˆ0ˆˆ＝4、混合测角网列条件式原则－－将混合网“化整为零”，拆分成各种基本图形（单三角形、中点多边形、大地四边形和扇形P84例5-3）5、测角网极条件个数判断极条件：其作用是使经过不同三角形（指推算路线）推算的同一条边长具有相同的长度。极条件个数：极条件仅存在于有角度观测值的控制网中。知点）网中所有点数（包括已－－测角网中的总边数极plplr32二、测角网按方向条件平差（独立网）必要观测数：t=2p+设站点数Mp－－待定点数1、单三角形n=6t=2+3=5r=12、中点多边形n=12t=2×2+4=8r=4测方向网无圆周角条件。测方向网的图形条件：13极条件：图形条件：0110115612Wvvvvvv3、大地四边形r=n-t=12-8=4三、附合测角网（非独立网）条件平差附合条件数＝多余已知数据个数条件式个数：r1－－独立网条件个数r2－－附合条件数（强制附合条件）13极图21rrtnr强制附合条件r2一般有三种形式：（1）边长条件（基线条件）条件方程个数=多余已知边数；（2）方位角条件（固定角条件）条件方程个数=多余已知方位角个数；（3）纵横坐标条件条件方程个数=多余已知点组个数的2倍已知点组：用已知边和已知方位角将已知点连接在一起时，称为一个已知点组。如附合导线就有两个已知点组。5.4边角网条件平差边角网条件方程的构成：必要起算数据：d=3,必要观测数据：t=P×2-1P－－网中待定点的关系描述观测边与观测角间余弦条件正弦条件测边图形条件与测角网相同极条件圆周角条件测角网图形条件654321常用正弦条件代替1、单三角形边角网：如图：观测了三角、三边。n=6，t=1+2=3注意：每个单三角形中只能有两正弦（或一正一余、或二余）条件，列三个正余弦条件必定相关。观测值：观测边－－用表示。观测角－－用表示。AcabCB112213余弦正弦或或余弦正弦图形rcba，，CBA，，正弦条件方程：0ˆsinˆˆsinˆˆsinˆˆsinˆAbBaBbAa即由上式得正弦条件的线性化条件方程：01cos1cossinsinWVBaVAbAVBVBAbaAbBaWsinsin式中闭合差为：取秒为单位。、单位，也应化成、为单位，则取、若BAbaVVcmbacmVV注意：•余弦条件方程：0ˆˆˆ2ˆˆˆcos222AcbacbarcAbcacbarcWbcaKWKAVKABVCVVAcba2cos0sinsincoscos222式中闭合差：其中设：由上式得余弦条件的线性化条件方程：2、中点多边形边角网（中点三边形为例）n=15,t=53、大地四边形n=14，t=5511361310正弦极条件圆周图形而不必正弦圆周图形r639正弦图形r4、测距仪的标称测距精度为：乘常数加常数basbammisi610Si－所测距离，以KM为单位。如索佳SET2100型全站仪：mmmmsiSppmSm22102265、边角网观测值的权通常设：6、边角网VTPV的计算有单位－－测边权无单位测角权则测角中误差22001sisipp角边pvvpvvPVVT5.5条件平差应用举例：一、单导线条件平差测角：n个，测边：n-1条，待定点个数：n-2，则r=n+n-1-(n-2)×2=31、方位角条件（设观测角为左角）1121横坐标条件纵坐标条件坐标条件方位角条件niCDABinnWWVVV1211800闭合差：2、纵、横坐标条件值。是由观测值计算的近似、、式中闭合差：00011001,11001,1111000111100001sin01cosiiinicnCiiBynicnCiiBxniniyiinSiininixiinSiiyxyyyyyWxxxxxWWVxxVWVyyV二、GPS平面网条件平差