[发明专利]一种毛竹胸径年龄联合分布动态模型的构建方法

权利要求书

1.一种毛竹胸径年龄联合分布动态模型的构建方法，其特征在于，包括以下步骤：

步骤一，收集研究区域以5年为一个复查周期的3期毛竹连续清查样地数据；

步骤二，将研究区域每期毛竹样地的连续清查数据按株数汇总成毛竹胸径年龄二维统计数据；

步骤三，根据步骤二汇总的各期毛竹林胸径年龄二维统计数据，计算各期数据毛竹胸径的边缘分布值u与年龄的边缘分布值v；

步骤四，以各期数据毛竹胸径与年龄边缘分布值为自变量，采用极大似然法，分别拟合若干种二元Copula密度函数，得每期数据的各二元Copula密度函数的参数向量β_i,j，其中，i为期数，取值为1、2、3，j为第j种二元Copula密度函数；

步骤五，根据AIC信息准则与确定系数得毛竹各期连续清查数据的胸径年龄最优二元Copula密度函数均为二元Gumbel Copula密度函数；

步骤六，根据步骤四与步骤五的计算结果，建立二元Gumbel Copula密度函数参数β与期数t的线性回归关系，关系式为：

β＝1.0926-0.0177t；

步骤七，由最优二元Copula密度函数及步骤六的参数β与期数t的线性回归关系，得出毛竹胸径年龄联合分布动态模型为：

式中t为毛竹连续清查数据期数，t＝1表示第1期，依次类推，β为最优二元Copula分布函数参数，u，v分别为毛竹胸径与年龄的边缘分布函数值。

2.根据权利要求1所述的一种毛竹胸径年龄联合分布动态模型的构建方法，其特征在于：所述若干种二元Copula函数包括二元正态Copula函数、二元t-Copula函数、二元Clayton Copula函数、二元Frank Copula函数及二元Gumbel Copula函数。

3.根据权利要求1或2所述的一种毛竹胸径年龄联合分布动态模型的构建方法，其特征在于：所述步骤五中AIC信息准则计算AIC值的计算式为：

式中：D(u_i,v_i；β)为所选二元Copula密度函数；β为所选二元Copula函数参数向量；m为所选二元Copula密度函数的参数个数；L为所选二元Copula密度函数的极大似然函数；N为要拟合的原始数据的样本量，当所选模型为二元Gaussian Copula、Frank Copula、GumbelCopula、Clayton Copula函数时，m等于1，所选模型为二元t Copula函数时，m等于2，u_i,v_i(i＝1，2，…N)分别为原始样本的第i个数据统计得到的胸径与年龄边缘分布值。

4.根据权利要求1或2所述的一种毛竹胸径年龄联合分布动态模型的构建方法，其特征在于：所述步骤五中确定系数公式为：

式中y_i为因变量的第i个实测值，为因变量的第i个估计值，N为样本个数，为因变量的平均值。