2009-现代测量平差与半参数估计

《2009-现代测量平差与半参数估计》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、现代测量平差与半参数估计陶本藻教授武汉大学测绘学院地球空间环境与大地测量教育部重点实验室函数模型是描述观测量与待求参数间的数学函数关系的模型，是确定客观实际的本质或特征的模型。随机模型是描述平差问题中的随机量（如观测量）及其相互间统计相关性质的模型。1、测量平差数学模型经典平差模型R（A）=UR（Q）=nX为非随机参数经典平差公式R（A）=t

2、义高斯—马尔柯夫模型，最小二乘统一理论Rao在文中提出的最小二小乘准则是R（A）=tR(D)=gn≥g>tu≥tX非随机陶本藻、刘大杰[5]（[5]1990）从奇异正态分布的密度函数3、平差系统的模型误差模型误差分为函数模型误差和随机模型误差两类最小二乘平差参数X的估值具有最优无偏性，单位权方差的估值具有无偏性和渐进最优性。这些良好的统计性质都是基于模型误差不显著的情况。但在实际平差系统中，由于种种原因的建模近似，例如非线性观测方程的线性化；未顾及或近似考虑某种系统误差影响；观测值的先验协方差阵不尽合理等

3、原因都会造成函数模型和随机模型产生误差。模型近似在回归拟合模型中则更为突出。4、测量中的半参数回归模型及估计方法4.1、测量中的半参数回归模型和平差估计准则（补偿最小二乘原理）半参数回归模型：平差估计准则：是一个适当给定的矩阵，称为正规化矩阵；是一个给定的纯量因子，在极小化过程中对和起平衡作用，称为平滑因子。4.2、正规化矩阵正定时半参数模型的估计由拉格朗日乘数法，构造函数参数估计：非参数估计：4.3、最小二乘核估计与附加系统参数平差附加系统参数估计核估计（偏核光滑估计）4.4、最小二乘配置与补偿最小二乘

4、半参数：配置：4.5、正规化矩阵的一个选择如果观测值是在时刻得到的一个时间序列，若假设相邻时刻的模型误差与的差别不应太大。因此可令：4.6、确定的周江文方法时间序列系统误差模型设为例如的协方差阵为5、半参数估计的自然样条函数法5.1半参数回归模型中的自然样条插值函数设为区间上的自然样条插值函数，为节点，且，满足插值条件：满足上述条件的插值函数中，自然样条函数是最光滑的。(Green,Silverman1987)5.2补偿最小二乘原理及其解(Fessler,1991)G矩阵为严格的对角占优矩阵，即正定矩阵。

5、设：，是半正定矩阵补偿最小二乘准则项表示为：按照求条件极值的拉格朗日乘数法，构造函数得到法方程为：法方程系数阵满足：时可逆。5.3平滑参数的确定（交叉核实法）去除时刻的观测值，也就是以作为观测值，对于给定的平滑参数得到观测值估计和系统误差估计。通过三次样条函数，可以内插出时刻的系统误差估计，记为，从而得到时刻的观测值估计如果平滑参数选取的比较合适的话，那么若选取的值，使(Green.Silverman1994)等价模型6、半参数模型假设检验(1)假设检验(2)（1）、（2）联合平差拒绝域7、时间序列分析与

6、残差分析7.1、AR（P）模型求解：平差求解可得F检验定阶给定显著性水平（例如＝0.05，0.01），查F分布表可得临界值若，则是合适的（检验不显著）；若，则是不合适的（检验显著）；7.2、半参数模型模型：Y=X+S+A式中，A=S=X=，=7.3、两模型的比较分别用AR(p)模型和半参数模型对香河台东西向长约700米的短水准观测数据进行计算比较。AR(p)模型求得的平差值图半参数模型求得的平差值图分别用AR(p)模型和半参数模型求得的平差值比较图上面的图中虚线表示用AR(p)模型求得的平差值，实线表示半

7、参数模型求得的平差值从图中可以看出用AR(p)模型求得的平差值的变化较大，而且快，曲线也不光滑；而用半参数模型求得的平差值的变化小，而且慢，曲线也非常光滑；由此可以看出当观测值中存在粗差或模型中存在模型误差时得出来的结果是不精确的或是错误的；而引入半参数后，就可以检查出一般模型不能检查出的系统误差或模型误差，从而提高结果的精度。谢谢！