《变形监测网非线性二类动态优化设计中观测权的解算方法》由会员分享，可在线阅读，更多相关《变形监测网非线性二类动态优化设计中观测权的解算方法（10页珍藏版）》请在金锄头文库上搜索。
1、测绘信息网 http:/测绘信息网测绘专业门户网站.变形监测网非线性二类动态优化设计中观测权的解算方法胡召玲摘要根据变形监测网非线性二类动态优化设计的数学模型,本文从一个简单的网形入手,导出测边网非线性二类动态优化设计数学模型的具体形式,并提出了采用非线性优化设计技术的方法投影梯度算法解算这种模型,从而将这种解算方法推广到求解任意网形的监测网非线性二类动态优化设计中的观测权。http:/ 测绘信息网关键词变形监测网; 非线性; 优化设计; 投影梯度中国图书资料分类号P221Solution Method of Observational Weight to the Nonlinear Seco
2、nd Dynamic Optimal Design of Deformation Monitoring NetworkHu Zhaoling(China University of Mining and Technology, Xuzhou, 221008)AbstractAccording to the mathematical model for the nonlinear second order dynamic optimal design of deformation monitoring network, a specific mathematical model of nonli
3、near second order dynamic optimal design of distance network is derived from a simple network. And gradient projection algorithm that is a kind of nonlinear optimal design technique method is put forward using solving the type of model. Thereby the solution method is popularized to solve observation
4、 weight of nonlinear second order dynamic optimal design of any deformation monitoring network.Key wordsDeformation monitoring network; Nonlinear; Optimal design; Gradient projection http:/ 测绘信息网0引言变形监测的主要目的,是对变形参数的探测。近年来用于变形监测的仪器和监测方法都有很大发展,对变形体的变形研究也越来越精细,因此对变形监测网的设计、观测和数据处理提出了更高的要求,变形监测网具有观测周期长、费
5、用高的特点,研究监测网的优化设计具有更重要的意义。而监测网的二类设计就是观测精度的选择和各类观测值权的确测绘信息网 http:/测绘信息网测绘专业门户网站.定问题。由于变形观测的精度要求高及其特殊要求,监测网的二类设计显得更有实际意义。国内外学者对变形监测网的优化设计已作过一些研究,但共同的作法都是采用线性规划设计方法。众所周知,在变形监测系统的优化设计模型中,客观上存在的绝大多数函模型关系都是非线性的,传统的作法是将非线性函模型按泰勒级数展开,略去二次及二次以上高次项,化为线性函数模型,用线性规划方法解算。但这是建立在近似值与其真值充分接近的基础上,研究表明,该条件难以满足,因而线性化必然会
6、影响到非线性函数模型的真实性,得不到满意的结果。科学发展到今天,有必要研究和应用非线性优化设计,尤其是非线性动态优化设计的函数模型及其解算方法。衡量变形监测网的质量准则,一般包括:精度、灵敏度、可靠性和观测费用。文献1根据精度要求直接构造了以变形参数精度为质量准则的非线性二类动态优化设计模型。针对这种模型,本文从一个简单的网形入手,导出二类优化设计数学模型的具体形式,并提出了采用非线性优化设计技术投影梯度算法解算这种模型,从而将这种解算方法推广到求解任意网形的监测网非线性二类动态优化设计。1变形监测网非线性二类动态优化设计的数学模型设有 n 个独立观测值 L 为某周期观测值,其改正数为 V,观
7、测值对角型权阵为 P,以相对于前周期的位移量 d 作为未知参数。将每个观测值的平差值 都表达成未知参数 d 的非线性函,即(1)用矢量表示: http:/ 测绘信息网V=f(d)-L(2)1.1目标函数以位移量 d 的协方差 Kdd作精度准则。观测值 L 为随机向量,其协方差矩阵为 KLL=P-1,在间接平差中,未知参数 d 是根据观测值 L 来计算的,则得未知参数 d 的协方差矩阵 Kdd为(3)为了计算 d/ L,需引用 d=( L)函数式,为此,由最小二乘原理 VTPV 最小的必要条件为 (VTPV)/ d=0 http:/ 测绘信息网测绘信息网 http:/测绘信息网测绘专业门户网站.
8、令 http:/ 测绘信息网(4)对(4)式求微分得:(5)为了根据(4)式计算 ,将 h(d)写成 ,其中(6)y 是 n 维向量, YT是 kkn 的矩阵, 是 nk 的矩阵, fs是 nk 维向量,令雅克比矩阵 ,则代入(5),(3)式得 http:/ 测绘信息网(7)取目标函数为:变形参数位移的协方差最小,即:(8)因为变形参数位移的协方差最小相当于它的协因数最小,也相当于它的权矩阵的最大,所以目标函也可变为: http:/ 测绘信息网测绘信息网 http:/测绘信息网测绘专业门户网站.(9)1.2约束条件由文献1,灵敏度约束为:(10)其中 0是单位权方差, 0是非中心化参数,可由显
9、著水平 ,检验功效 从有关诺谟图中查取, g 是给定的方向向量, d0是设计时预先要求或预先给定在某一方向 g 上要监测的最小位移量。(10)式也可变为:(11)令 Ag=K,则(11)式变为: http:/ 测绘信息网(12)采用矩阵的 Hadamard 积*,并将矩阵按对角线拉值后得:(KT*KT)p W(13)令 B=(KT*KT),得监测网二类设计的灵敏度约束:Bp W(14)根据文献2,监测网二类设计的可靠性约束:b(1) Dp b(2)(15)其中, D 为系数矩阵, b(1),b(2)为常数向量,可根据实际网形和单位权以及可发现的观测值粗差的最小值确定。监测网二类设计的费用约束为
10、: http:/ 测绘信息网Sp w(16)其中, S=(1 1 1),w 为一给定的限值。由式(9),(14),(15),(16)得变形监测网非线性二类动态优化设计的数学模型。测绘信息网 http:/测绘信息网测绘专业门户网站.2测边网数学模型 http:/ 测绘信息网如图所示,在固定点 A、 B、 C 的三角形插入一点 T,使待定点 T 的协方差阵最小(相当于点位误差最小), A、 B、 C 的坐标值及 T 的经一类优化设计后的坐标值列于表中。采用边长观测法,共测得三条边,观测值分别为 S1,S2,S3,求监测网非线性二类动态优化设计的结果 P1,P2,P3。图 1测边监测网表 1各点的坐
11、标点 X 坐标 Y 坐标A 1400m 100mB 1400m 1600mC 100m 840mT 950m 850m令 A、 B、 C、 T 点的坐标分别为( xA,yA),(xB,yB),(xC,yC),(x1,y1),则得平差值方程为: 雅克比矩阵 A 为: http:/ 测绘信息网因此 http:/ 测绘信息网测绘信息网 http:/测绘信息网测绘专业门户网站.(17)令 ,可以求得令 ,可以求得(18) http:/ 测绘信息网令(19)其中, 测绘信息网 http:/测绘信息网测绘专业门户网站.由(17)、(19)式得目标函数为: http:/ 测绘信息网F=-tr HTN-1H=
12、min(20)其中 H,N 均为 22 阶的函数矩阵, Hij,Nij为 p1,p2,p3的函数。因为约束条件均属线性约束,令其为:GTp b(21)其中, p=(p1,p2,p3)T,G,b 分别为约束系数矩阵和常数项。当观测值个数为 n,坐标未知数个数为 k 时,相应地,式(20)中的 H,N均为 kk 阶的函数矩阵, Hij,Nij为 p1,p2,pn的函数。3算法 http:/ 测绘信息网根据变形监测网纯量精度准则 F 最优,要使 QF=-fTHTN-1Hf=min,其中f 是变形参数 d 的线性化函数的系数向量。因此令 F=QF,变形监测网非线性二类动态优化设计的数学模型为目标函数
13、F=QF=-fTHTN-1Hf=min(22)约束条件 GTp b(23)这是仅带有线性约束条件的最优化问题。令 x=(x1,x2,xn)T,f(x)=F,利用投影梯度法进行迭代求解 x。投影梯度法是求解线性约束的非线性优化问题最有效的方法之一,它的基本思想是产生一个可行点列 x(k),满足 f(x(k+1) f(x(k),使得 x(k)收敛于约束问题的极小点(或 KT 点),或者使得 x(k)的任意聚点为一 KT 点。因此在 x(k)处的搜索方向不仅要是一个下降方向,而且要是一个可行方向。而对于线性约束问题,过边界点的任一方向在边界上的投影,都是可行方向,而负梯度方向的投影,将是一个下降方向
14、。根据文献3,当 f(x)连续可微,且任意可行点处有效约束的系数向量线性无关,则算法或者终止于一个 KT 点,或者生成一个无穷序列 x(k),其任意聚点为一 KT 点。目标函数 f(x)的梯度为:其中, ,因为 NN-1=I,I 是单位矩阵,将此式的两边同时对 xi求偏导得 http:/ 测绘信息网测绘信息网 http:/测绘信息网测绘专业门户网站.于是有利用投影梯度法求解优化模型的迭代步骤如下:(1) 给出初始可行点 x(1),令 k=1;(2) 确定 x(k)处 q 个有效约束的法向量 gi所构成的 nq 矩阵 Gq由式 Pq=I-Gq(GTqGq)-1GTq计算投影矩阵 Pq,若 q=0,则令 Pq=I;(3) 令 p(k)=-Pq f(x(k),若 p(k)=0 则转(5),否则转(4);(4) 计算 ,求 k使令 x(k+1)=x(k)+ kp(k),转(6);(5) 计算 =( GTqGq)-1GTq f(x(k),若 0,则 x(k)为 KT 点,停。否则,求 i=min i0,从 Gq中取掉对应于 i的向量 gi,得 Gq-1,由式Pq-1=I-Gq-1(GTq-
