森林资源调查数据的稳健估计及分析-中国科技论文在线.PDFVIP

森林资源调查数据的稳健估计及分析 梁长秀 冯仲科 郎南军 黄华国 周科亮 北京林业大学资源与环境学院 Email:Fengzhongke@126.com 摘 要:该文指出 : 以往森林经理调查中 ,航片判读小班蓄积与地面实测蓄积之间的统计关 系多建立在经典的最小二乘原理之基础上 . 当调查数据中不含有粗差时 ,最小二乘估计可以 得到较好的结果 . 然而 , 由于各种因素的影响 ,观测成果中总会含有粗差 .利用稳健估计 , 可以使含有粗差的观测成果得到较为理想的估计值 .该文通过稳健估计的原理介绍和算法 推证 , 以及三个实例分析 ,证明了稳健估计原理可使森林调查得到更好的结果 . 关键词: 最小二乘估计 ,森林调查 ,稳健估计 建立在∑ｎｉ= 1 ＰｉＶ2 ｉ= ｍｉｎ基础上的参数估计是基于经典的最小二乘原理 , 已在森林资源调查中得到了广泛的应用 .理论和实践证明 ,观测 (调查 )结果中总是不可避 免的要产生误差和粗差 .对于误差 ,可采用最小二乘原理平差解决 ;对于粗差 ,则需要外业 剔除 ,外业不能发现的粗差则需要内业通过一定的理论或算法剔除 ,稳健估计就是一种处理 粗差的现代理论[1～ 1 5]. 1 稳健估计原理及算法 设有某森林调查统计回归问题为 : Ｚｉ = ａ0 + ａ1 ｘ1 + ａ2 ｘ2 +… + ａｍｘｍ ( ｉ=1 ,2 ,3 ,… , ｎ) (1 ) 当有ｎ个点时 ,可表示为矩阵形式 : Ｖ= ＢＸ- Ｚ (2 ) (2 )式中 : Ｖ =( Ｖ1 Ｖ2 … Ｖｎ) Ｔ Ｘ =( ａ0 ａ1 ａ2 … ａｍ) Ｔ Ｚ =( Ｚ1 Ｚ2 … Ｚｎ) Ｔ Ｂ = 1 ｘ1 1 ｘ1 2 …ｘ1ｍ 1 ｘ2 1 ｘ2 2 …ｘ2 ｍ 　　　　　 1 ｘｎ1 ｘｎ2 …ｘｎｍ 依最小二乘原理 , 即ＶＴＶ = ｍｉｎ ,令 2 ＶＴｄＶ ｄＸ =0 ,则 ＶＴＢ =0 或ＢＴＶ=0 ＢＴ( ＢＸ- Ｚ) =0 从而有Ｘ =( ＢＴＢ) -1 ＢＴＺ (3 ) (3 )式为最小二乘解 .最小二乘解最大的缺点在于对于粗差的不敏感性 , Ｖ2 ｉ 成几何倍数 增加 ,有一个观测异常点 ,就会破坏整个趋势面 . 因此 ,必须选择对于粗差 (异常点 )不敏 感的稳健估计理论来建立森林调查统计的空间模型 .如选择最小一乘原理 ( ∑｜Ｖｉ｜= ｍｉｎ ,残差绝对值之和为最小 ,又称一次范数最小 ) ,使异常点 (通常含有 粗差 ,| Ｖｉ|较大 ) 占有很小的权引入权阵 ( ｃ为小常数 ) .定权时 ,为了避免因ｖ =0 而出 现的计算问题 ,可取 ｐｉ =1 ｜ｖｉ｜+ ｃ (4) 综上所述 ,一次范数最小估计的步骤是 : 转载 中国科技论文在线