参数变换在参数估计中的应用.docxVIP

参数变换在参数估计中的应用 1 实验与模型的转换 现阶段，中国大部分测量结果如北京54坐标系和西安80地形坐标系。ccs2000标准将于2007年7月1日开始。因此，在不久的将来，将有一个问题改变北京54和ccs2000、西安80和ccs2000的坐标。对于平面坐标系的转换,常规的做法在理论上存在一些缺陷,即在求解过程中,往往把某一坐标系下的坐标当作是不存在误差的,仅对另一坐标系下的坐标进行改正。而实际上,测量获得的坐标成果都是有误差的,那么此时的最小二乘估计结果从统计的观点看就不再是最优的,而是有偏的。为了解决这个问题,本文推导了利用整体最小二乘求取坐标转换参数的公式,并用实测数据验证了算法的可行性,从而解决了传统理论上存在的缺陷。另外,相关研究表明在不同地区采用不同的坐标转换模型对转换结果也会有影响,为了适应不同工程的要求,也给出了平面四参数模型和仿射变换模型的参数估计公式。 2 算法导出 2.1 u3000定义 设有两个坐标系A系与B系,一组点在两套坐标系下的坐标分别为(XA,YA)i、(XB,YB)i,各点对应的平差坐标为(X′A,Y′A)i、(X′B,Y′B)i。 四参数转换模型为: [XY]B=m[cosαsinα-sinαcosα][XY]A+[ΔxΔy](1)[XY]B=m[cosα?sinαsinαcosα][XY]A+[ΔxΔy](1) 式中,[XY]A[XY]A为点在A坐标系的平面坐标,[XY]B[XY]B为点在B坐标系的平面坐标,m为由坐标系A到坐标系B转换的尺度参数,α为由坐标系A到坐标系B转换的旋转参数,[ΔxΔy][ΔxΔy]为由坐标系A到坐标系B转换的平移参数。 要求转换参数在满足式(1)条件下是残差平方和∑v2i2i=min的最优解,即: ∑v2i2i=∑(X′A-XA)2+∑(Y′A-YA)2+∑(X′B-XB)2+∑(Y′B-YB)2=min (2) 另,由式(1)可得: {X′B=X′Amcosα+Y′Amsinα+ΔxY′B=-X′Amsinα+Y′Amcosα+Δy(3){X′B=X′Amcosα+Y′Amsinα+ΔxY′B=?X′Amsinα+Y′Amcosα+Δy(3) 将式(3)代入式(2),有: ∑v2i2i=∑(X′A-XA)2+∑(Y′A-YA)2+∑(X′Amcosα+Y′Amsinα+Δx-XB)2+∑(-X′Amsinα+Y′Amcosα+Δy-YB)2=min (4) 将式(4) 设为F,对于具有n个公共点的观测数据来说,要求的参数是2n个点的坐标平差值以及4个坐标转换参数。即求得(X′A,Y′A)i、(X′B,Y′B)i和α\,m\,Δx\,Δy使得F取得极小值。 这样就转换成了直接求极值的问题,这与传统的用拉格朗日条件求极值的方法是等价的。下面就是要寻找参数的估值来使得上面的目标函数达到最小。对于式(4),此时的待求参数减少为n个点的坐标平差值以及4个坐标转换参数。把F分别对未知参数求导并令其等于零,有: ?F?XAi=(a2+b2+1)X′Ai+aΔX-bΔY-aXBi+bXBi-XAi=0?F?YAi=(a2+b2+1)Y′Ai+bΔX+aΔY-bXBi-aYBi-YAi=0?F?Δx=a∑X′Ai+b∑Y′Ai+nΔX-∑XBi=0?F?Δy=-b∑X′Ai+a∑Y′Ai+nΔY-∑YBi=0?F?a=(∑(X′Ai)2+∑(Y′Ai)2)a+ΔX∑X′Ai+ΔY∑Y′Ai-∑XBiX′Ai-∑YBiY′Ai=0?F?b=(∑(X′Ai)2+∑(Y′Ai)2)b+ΔX∑Y′Ai-ΔY∑X′Ai+∑YBiX′Ai-∑XBiY′Ai=0(5) 根据文献的结论,式(5)已经对坐标转换参数进行了整合,即:a=mcosα,b=msinα。式(5)就是待估参数所满足的方程组,从式(5)求解式中的(2n+4)个参数。 2.2 u3000计算公式 由于将式(5)线性化后会面临方程不收敛问题,所以本文采用参数变换方法求解,把(2n+4)个方程分成两部分,即对坐标求导的2n个公式,以及对转换参数求导的4个公式: {X′A1=-(aΔX-bΔY-aXB1+bXB1-XA1)a2+b2+1Y′A1=-(bΔX+aΔY-bXB1-aYB1-YA1)a2+b2+1?X′An=-(aΔX-bΔY-aXBn+bXBn-XAn)a2+b2+1Y′An=-(bΔX+aΔY-bXBn-aYBn-YAn)a2+b2+1(6) 从式(6)可以看出,有了原始的坐标观测值以及坐标转换参数的平差值就可以求得A系坐标的平差值,进而由式(3)求得B系坐标平差值。根据式(5)中对转换参数求导的公式整理得: [∑(X′A1)2+∑(Y′A1)20∑X′A1∑Y′A10∑(X′A