2009现代测量平差跟半参数估计.ppt

现代测量平差与半参数估计 陶本藻教授 武汉大学测绘学院 地球空间环境与大地测量教育部重点实验室 4.3、最小二乘核估计与附加系统参数平差 附加系统参数估计 核估计（偏核光滑估计） 4.6、确定 的周江文方法 时间序列系统误差模型设为 例如 的协方差阵为 等价模型 (1) 假设检验 (2) （1）、（2）联合平差 拒绝域 7、时间序列分析与残差分析 7.1、AR（P）模型 求解： 平差求解可得 给定显著性水平（例如＝0.05，0.01）， 查F分布表可得临界值 若 ，则 是合适的（检验不显著）； 若 ，则 是不合适的（检验显著）； 7.2、半参数模型 模型：Y=X +S+A 式中， A= 7.3、两模型的比较 分别用AR(p)模型和半参数模型对香河台东西向长约700米的短水准 观测数据进行计算比较。 上面的图中虚线表示用AR(p)模型求得的平差值，实线表示半参数模型求得的平差值从图中可以看出用AR(p)模型求得的平差值的变化较大，而且快，曲线也不光滑；而用半参数模型求得的平差值的变化小，而且慢，曲线也非常光滑；由此可以看出当观测值中存在粗差或模型中存在模型误差时得出来的结果是不精确的或是错误的；而引入半参数后，就可以检查出一般模型不能检查出的系统误差或模型误差，从而提高结果的精度。 谢 谢！ F检验定阶 S= X= ， = AR(p)模型求得的平差值图 半参数模型求得的平差值图 分别用AR(p)模型和半参数模型求得的平差值比较图 * 函数模型是描述观测量与待求参数间的数学函数关系的模型，是确定客观实际的本质或特征的模型。 随机模型是描述平差问题中的随机量（如观测量）及其相互间统计相关性质的模型。 1、测量平差数学模型 经典平差模型 R（A）=U R（Q）=n X为非随机参数 经典平差公式 R（A）=tu d=u-t R(Q)=n X非随机 秩亏自由网平差 R（Q）=gn R(A)=u X非随机 具有奇异协方差的平差模型 配置（似合推估）模型 广义测量平差原理 最小二乘统一理论 2、广义高斯—马尔柯夫模型，最小二乘统一理论 Rao在文中提出的最小二小乘准则是 R（A）=t R(D)=g n≥gt u≥t X非随机 陶本藻、刘大杰[5]（[5]1990）从奇异正态分布的密度函数 3、平差系统的模型误差 模型误差分为函数模型误差和随机模型误差两类 最小二乘平差参数X的估值具有最优无偏性，单位权方差的估值具有无偏性和渐进最优性。这些良好的统计性质都是基于模型误差不显著的情况。 但在实际平差系统中，由于种种原因的建模近似，例如非线性观测方程的线性化；未顾及或近似考虑某种系统误差影响；观测值的先验协方差阵不尽合理等原因都会造成函数模型和随机模型产生误差。模型近似在回归拟合模型中则更为突出。 4、测量中的半参数回归模型及估计方法 4.1、测量中的半参数回归模型和平差估计准则（补偿最小二乘原理） 半参数回归模型： 平差估计准则： 是一个适当给定的矩阵，称为正规化矩阵； 是一个给定的纯量因子，在极小化过程中对 和起平衡作用，称为平滑因子 。 4.2、正规化矩阵正定时半参数模型的估计 由拉格朗日乘数法，构造函数 参数估计： 非参数估计： 4.4、最小二乘配置与补偿最小二乘 半参数： 配置： 4.5、正规化矩阵的一个选择 如果观测值是在时刻 得到的一个时间序列，若假设相邻时刻的模型 误差与 的差别不应太大。因此可令 ： 5、半参数估计的自然样条函数法 5.1半参数回归模型中的自然样条插值函数 设 为区间 上的自然样条插值函数， 为节点，且