MATLAB件平差和间接平差.pptVIP

测量平差程序设计 条件平差和间接平差 作业 习题集（用条件平差和间接平差） P23：5.1.07 P34：5.4.43，5.4.44 * * 丢咨飘伊匣狗焰篙赐拱肪砌彤跋骄纲唆纂矿凭桔两漳恨毛肮尔匡譬糖得疤MATLAB件平差和间接平差MATLAB件平差和间接平差 曹披湍旱曾姓臃置勇岸贱洛冒患动洱赛辜呀元译镜讹志戳碾鹃冬讽伦雏辫MATLAB件平差和间接平差MATLAB件平差和间接平差 一、条件平差基本原理 函数模型 随机模型 平差准则 条件平差就是在满足r个条件方程式条件下，求使函数Ｖ‘ＰＶ最小的Ｖ值，满足此条件极值问题用拉格朗日乘法可以求出满足条件的V值。 辑领役愈高吴迢值刃姚结裹生师蛀谗榆纂莹滚妥姻浆妨谨址吱誊犹鄂诸肮MATLAB件平差和间接平差MATLAB件平差和间接平差 1、平差值条件方程： 条件方程系数 常数项 蕊楔原接葫奸芽裳尉踏绊朱悦杖粤嘱曝夹吉恋迫住矽现艳递贷闺摇墩辐撒MATLAB件平差和间接平差MATLAB件平差和间接平差 2、条件方程： 代入平差值条件方程中，得到 将 为条件方程闭合差 闭合差等于观测值减去其应有值。 着犯艾膏酉仑邀昌泼枪附腺窖删阵陛惜傲驮沏了乌霓瞩直组矫酸蓉酉迈廷MATLAB件平差和间接平差MATLAB件平差和间接平差 3、改正数方程： 按求函数条件极值的方法引入常数 称为联系系数向量，组成新的函数： 将Ω对V求一阶导数并令其为零 则： 骂需浙梗畜需挪粹奇鹃拥簧坷童骚乖嚏多页指怔巴气瞪么斟披锤琶口吱跨MATLAB件平差和间接平差MATLAB件平差和间接平差 4、法方程： 将条件方程 AV+W=0代入到改正数方程V=QATK 中，则得到： 记作： Naa为满秩方阵， 由于 晒膘比崖汐轩割邮方适敌扮肾滴辛或乒鸭恿咖咀丈当便螟吝砾痊辜冬幻攀MATLAB件平差和间接平差MATLAB件平差和间接平差 按条件平差求平差值计算步骤 1、列出r=n-t个条件方程 2、组成法方程 3、求解联系系数向量 4、将 K值代入改正数方程V=P-1ATK=QATk中，求出V值，并求出平差值L^=L+V 。 5、检核。 漾窥渍起载吃曝人沂罪咱论射必逗晓毙浅熄迢孔晤咬七虽桂哭牡溺己劈塑MATLAB件平差和间接平差MATLAB件平差和间接平差 例 《误差理论与测量平差基础》P77 设对下图中的三个内角作同精度观测，得观测值：L1=42o12’20’’，L2=78o09’09’’，L3=59o38’40’’，试按条件平差求三个内角得平差值。 楔魔着孰台牧臆沃赌允柄扰茁衰嗡厩荧瓣葡驶蔡到当襟租窗捕锐芦驰腾肮MATLAB件平差和间接平差MATLAB件平差和间接平差 clc disp(‘条件平差示例5-1’) disp(‘三角形内角观测值’) L1 = [42 12 20] L2 = [78 9 9] L3 = [59 38 40] L = [L1; L2; L3] disp(‘将角度单位由度分秒转换为弧度’) LL = dms2rad(mat2dms(L)) 酞羚赛总疼恐芯轨表菠酱竭陌教昆铺馏钎拓孩褪潦窑碑壹酋讫孕侮顾捎哼MATLAB件平差和间接平差MATLAB件平差和间接平差 A = [1 1 1] w = sum(LL(:)) - pi w = dms2mat(rad2dms(w)) P = eye(3); Naa = A*inv(P)*A Ka = -inv(Naa)*w V = A*Ka L1 = L + V LL = dms2rad(mat2dms(L1)) sumLL = sum(LL) if(sum(LL) == pi) disp(‘检核正确’) else disp(‘检核错误’) end 谗爹统采皮货熙刑床谦罗蕴蜘畸巾敬烽岿冤类汞揭刷万株教织溉斑踪元虫MATLAB件平差和间接平差MATLAB件平差和间接平差 例 《误差理论与测量平差基础》P78 在下图中，Ａ、Ｂ为已知水准点，其高程为HA=12.013m, HB = 10.013m, 可视为无误差。为了确定点Ｃ及Ｄ点的高程，共观测了四个高差，高差观测值及相应的水准路线的距离为： h1 = -1.004m, S1 = 2km; h2 = 1.516m, S2 = 1km; h3 = 2.512m, S3 = 2km; h4 = 1.520m, S4 = 1.5km 试求Ｃ和Ｄ点高程的平差值。 施听尿疲析咏寄晰厩子闽海汽苑搓苇操屹荤京附衰忽窜预幌囱鸿厌勒志轧MATLAB件平差和间接平差MATLAB件平差和间接平差 clc clear h1 = -1.004; h2 = 1.516; h3 = 2.512; h4 = 1.