[发明专利]一种基于迭代加权最小二乘估计法的木材密度测定方法

摘要

一种基于迭代加权最小二乘估计法的木材密度测定方法，涉及木材密度测定的估计方法。本发明为了传统的回归分析方法等方法对于批次木材的密度测定进行估计存在误差较大的问题。本方法以有限样本下具有无偏估计性质的期望平均估计作为初始均值估计值，以迭代加权最小二乘法做回归估计微调，权系数取样本密度的频数与误差距离倒数的乘积，可以任意收敛到设定的门限阈值内而不出现发散现象，能够较准确估计反映出总体样本的密度。本发明适用于批次木材的密度估计领域。

主权项

一种基于迭代加权最小二乘估计法的木材密度测定方法，其特征在于包括以下步骤：步骤1、针对一批木材随机抽取试样，计算出各试样密度ρi；步骤2、试样密度的期望平均估计：步骤2.1、试样密度的频数统计：从ρi中确定出极大值ρmax和极小值ρmin；取a＝[ρmin*10l‑1‑0.5]/10l‑1，b＝[ρmax*10l‑1+0.5]/10l‑1；l为ρi、ρmax或ρmin小数点之后的有效位数，[·]运算表示取整数；将区间[a,b]等分成m个小区间，计算落入各小区间的试样密度个数pj，此pj即为频数，总的取样数步骤2.2、期望平均估计：由式(2)可计算出试样密度的期望平均估计ρ^e=Σi=1npi·ρiΣi=1npi---(2)]]>每个ρi所对应的频数为自身落入小区间所对应的频数pi；n个ρi会对应n个pi；步骤3、密度的加权最小二乘迭代估计：各试样的体积vi构成矩阵V＝[v1 v2 … vn]T，质量gi构成矩阵G＝[g1 g2 … gn]T，密度ρi构成矩阵P＝[ρ1 ρ2 … ρn]T，将ρ1至ρn所对应的n个pi构成对角矩阵W；步骤3.1：a、将代入公式(3)中，得出初始误差距离E0：E0=|ρ^eV-G|---(3)]]>b、将E0扩成对角阵，则初始权系数由式(4)确定；Wwls0=WE0-1---(4)]]>c、用式(5)的加权最小二乘法计算初始密度估计ρ^wls0=(VTWwls0V)-1VTWwls0G---(5)]]>步骤3.2：d、将代入公式(6)中，得出密度估计的第k次误差距离Ek；Ek=|ρ^wlsk-1V-G|---(6)]]>其中，是第k‑1次估计的密度；e、将Ek扩成对角阵，则第k次权系数由式(7)确定；Wwlsk=WEk-1---(7)]]>f、用式(8)的加权最小二乘法计算第k次密度估计ρ^wlsk=(VTWwlskV)-1VTWwlskG---(8)]]>步骤3.3：g、用式(9)计算密度估计的第k次方差σk2=1n-1Σi=1n(ρi-ρ^wlsk)2---(9)]]>h、δ为设定的阈值，判断如不满足条件则转入步骤3.2继续加权最小二乘迭代估计计算，如满足条件则结束并输出密度估计值