施工控制网的优化设计_顾利亚

施工控制网的优化设计

顾利亚　岑敏仪

(西南交通大学　测量工程系　成都　610031)

【摘　要】　根据施工控制网的特点,提出了用解析法进行控制网优化设计的新方法,介绍了在平

均可靠率和精度的约束下使用0-1规划进行网形设计的算法。实例验证,精度函数增量的“A ”标准和“E ”

标准均可作为控制网图形设计的目标函数。【关键词】　优化设计;0-1规划;测量控制网【分类号】　T P 391.41;T U 198

根据作业的过程,通常将施工控制网的优化设计划分为四个阶段,即:零类设计、一类设计、二类设计和三类设计。零类设计是控制网参考系或基准的设计问题,它包括数据处理的方法和坐标系的选择,不同用途的控制网选择不同的数据处理方法。由于施工控制网要考虑相对点位的精度问题,因此零类设计通常采用传统的习惯做法。一类设计是控制网的网形设计问题,是在预定测量精度的前提下,确定最佳的点位概略坐标和联系方式。控制点的设计位置,主要受施工放样的需要及地形和设备条件的制约,有些因素目前还很难用数学的方式表示。而控制网的图形(即控制点之间的联系方式)对网的图形强度影响较大,它是一类设计的主要研究内容,亦是本文的核心内容。二类设计是控制网在图形固定的前提下,寻求最佳的精度配置,它是控制网优化设计的热点问题。三类设计则是对已有控制网的改善,它一般要包含零类、一类和二类设计。

施工控制网优化设计的作用,是使所求解的控制网的图形和观测纲要在高精度、高可靠性及低成本意义上为最优。本文针对施工控制网设计的特点,在其图形设计中建立求解模型,使求出的图形和观测纲要同时满足预先规定的优化设计指标。

1　优化设计指标

控制网的优化设计指标包括精度、可靠性和经济费用指标。精度指标一般通过精度约束函数来满足。可靠性分为内部可靠性和外部可靠性,常用的指标有:观测量的多余观测分量、可发现粗差的下界值、外部可靠性尺度等。这些指标均对某些特定的条件有显著作用。根据施工控制网的特点,其可靠性指标可用平均可靠率来表示[1]

r 0=r /n

(1)

式中,r 为多余观测数,n 为总观测数。

控制网的费用标准一般可用下式表示

收稿日期:1996-10-08 顾利亚:女,1956年生,讲师。

第32卷第2期1997年4月

西南交通大学学报

JOU RNAL OF SOU THWEST JIAOT ONG UNIVERSITY

Vo l.32N o.2A pr.　1997　

S=∑CP(2)式中C称为费用系数,它是目前优化设计中最难确定的参数,也是数学优化设计方法中有待进一步研究的课题。P为观测值的权。通常,用观测次数的多寡来表征同一量纲观测量精度的费用,可用线性函数表示,而补充或删除一个观测量与增减重复观测次数的费用之间的函数关系则较难确定。控制网最终的优化结果,是各个阶段优化设计的总和。因此,在各个阶段的优化设计上不必强求同时满足精度、可靠性和费用指标,而最后的优化设计结果中达到这三项指标便可。因此,首先利用控制网的完全观测图形,在一定的平均可靠率和精度约束下,解算出最佳的观测图形,然后在此图形设计的基础上求解满足精度约束条件、费用最省的观测方案,这样,分两步将控制网图形与观测纲要优化设计用解析法直接求解。施工控制网的观测纲要优化设计的具体方法文献[2]已有研究。为节省篇幅,下面主要介绍控制网图形优化设计的方法。

2　图形优化

设施工控制网初始设计方案的平差模型为

V n×1=A

n×t

X

t×1

-L

n×1

L～N(A X~, 20P-1)

(3)式中L为观测向量;V为残差向量;A为模型系数阵,且R(A)=t; 20为方差因子;P为观测量权的对角线阵。初始设计方案的参数协方差阵为

D X= 20N-1X= 20(A T P A)-1(4) 设删除第i个观测量,其模型系数为A i,权系数值为P i,则方案修改后的平差模型为

V1=A1X-L1

L1～N(A1X~, 20P1)-1

(5)且有

V=V1

V i

A=

A1

A i

L=

L1

L i

删除观测量L i后,其参数协方差阵的增量为

D X i= 20[(A T1P1A1)-1-(A T PA)-1]= 20[(A T PA-A T i P i A i)-1-(A T PA)-1]

=-D X A T i(- 20P-1i+A i D X A T i)-1A i D X(6)由于A i为1×t的行矩阵,故A i D X A T i为标量,则有

D X i=-D X A T

i A i D X

- 20P-1i+A i D X A T i

(7) 当控制网为两级网,采用固定数据平差的数据处理方法,则完整地考虑原始数据误差影响的控制网参数协方差阵为[2]

D X=

D Y D Y,Z

D Z,Y D Z

=

D Y-D Y L T ZY G T Z

-G Z L ZY D Y 20N-1Z+G Z L ZY D Y L T ZY G T Z

(8)

式中

D Y= 20Y N-1Y= 20Y(A T Y P Y A Y)-1161

第2期 顾利亚等:施工控制网的优化设计

N Z=A T Z P Z A Z G Z=N-1Z A T Z P Z L ZY=( l Z Y)

删除观测量L i,等于该观测量的权P i=0,故控制网的参数协方差阵为

D X1=D X Pi=0

协方差增量为

D X i=D X1-D X(9)若按控制网的精度要求建立的权函数为

=FX

则控制网的精度函数为

D =FD X F T(10)精度函数的增量为

D i=F D xi F T

在控制网的纯精度指标中,如果 1, 2,…, t为参数协方差阵D X相应的特征值,则:

(a)若t r(D X)= 1+ 2+…+ t=min,称A最优;

(b)若 max=m in,称E最优, max为矩阵的最大特征值。

根据文献[3]的研究,矩阵D X中的每个元素不但受到D X的所有特征值 的影响,而且还受到其相应特征向量的影响。假如单独比较特征值,它只能获得抽象的数学概念。考虑到施工控制网的精度指标就是使精度函数D 的主对角线元素满足对应的限差要求。以精度函数增量的“A”最优或“E”最优作为优化设计中的目标函数,控制网的图形优化设计的数学模型为:

(1)“A”最优

t r( D i)=m in i∈n(11)即以 D i的迹最小为目标函数来删除观测量L i;

(2)“E”最优

diag( D i)max=min i∈n(12)即以所有的 D i最大的主对角线元素中,挑选最小的 D i max作为目标函数来删除观测量L i。以精度函数和平均可靠率为约束条件

diag(D )≤G

∑sign(P)≤t/(1-r0)

sign(P)=1,P i>0

0,P i=0

(13)

式中,G为精度函数的限差列阵;t为必要观测量个数。

参照整数规划中的0-1规划法,对式(11)(或式(12))和式(13)的数学模型进行控制网的图形优化,具体步骤为:

(1)建立控制网的完全观测图形,并按量纲相同的观测量等精度观测的要求给权赋初值,使diag(D )≤G;

(2)计算 D i,i=1,2,…,n;

(3)比较目标函数t r( D i)=m in;

162西　南　交　通　大　学　学　报 第32卷