林木脑径与材积的关系——Logistic衍生模型.doc

1、林木胸径与材积的关系Logistic 衍生模型张连翔，林 阳（辽宁省干旱地区造林研究所，辽宁 朝阳，122000）摘 要： 基于著名的 Logistic 模型，以时间因子作媒介，推导出一个能较准确描述林木 关系的一元VD材积模型Logi stic 衍生模型。这是一个既有理论依据，又有实用价值的新模型。实践证明， Logistic衍生模型具有较强的数据拟合能力和广泛的适用性，除林木 关系外，尚可应用于林木 关系和VDH植病 关系等项研究，具有理论和实用上的双重意义。SI-关键词： Logistic 衍生模型；林木 关系； 关系；一元材积式VDH分类号： S7571林木胸径与材积的关系（ 关系）是

2、测树学中的一个重要概念，作为一元材积式，广泛应用于森林调查中。特别是随着电子计算机的普及应用，不仅使材积表的准确度得到提高，更使复杂模型的求解以及多个材积式的选优与检验成为可能。尽管迄今为止各国林学家相继提出了几十个类型各异的一元材积模型 1，然而令人遗憾的是，有关方法和模型大都属经验式的，缺乏理论依据，同时模型的拟合精度也不很尽如人意。因此，有必要从理论上研究胸径与材积的关系，探讨建立既有理论依据，又有较高拟合精度的一元材积回归模型。基于 Logistic 模型，笔者推导出一个能准确描述林木胸径与材积关系的新模型Logistic 衍生模型，经实例验证，效果良好，具有理论和实用上的双重意义。1

3、 Logistic 衍生模型的导出众所周知，林木生长包括 3 个基本过程，即细胞分裂、细胞延长和细胞分化。虽然从理论上讲细胞和组织的生长潜力是无限的，可按指数式增长，但因受“环境阻力”的影响，使整个生长曲线呈现慢快慢停的“S”形生长节律。显示同样生长进程的调查因子有直径、树高、材积和质量等。因此可见，在林木生长过程中，时间始终是一个重要因子，利用时间动态方程描述林木生长进程，已有不少学者进行过较为深入的探讨 1。本文的目的，是利用林木生长进程中的时间因子作媒介，推导出一个新的 关系模型。VD在林木生长进程中，胸径和材积的生长大体可视为 2 个同步过程，时间是它们的共同因子。对于这 2 个生长过

4、程，可分别用方程（1） 、 （2） 表示 2,3：（1）tbaKD11exp（2）V22整理方程（1） 、 （2） ，解出 ，可得： t（3）111ln（4）222lbVKat显然，在林木生长进程中，对应于每一个 ，林木的 2 个表观生长过程均有相应的it和 ，这些 数值的关系及其变化，便构成了林木的 关系。iDiV VD联立式（3） 、 （4）并整理得：（5）bKa12式中： ； 。1212expbaa12为明了起见，用 ，进一步变换式（5） ，容易得到：KVD代，代（6）b式中： 材积； 胸径； 环境条件所允许的最大材积； 环境条件所允许VV DK的最大胸径； 待定参数。且有 0， 0，

5、（0， ） ， （0， ） 。ba,aVV在此，把式（6）称之为 Logistic 衍生模型。 2 Logistic 衍生模型的性质Logistic 衍生模型具有以下几点重要性质： 因 0 和 0，当 0 时， 的极限为零，当 时， ，这DDV种特征符合一般 关系。V 在开区间（0， ）内，模型有连续的导数，形式为：K（7）221dbDVab式（7）左边的值恒为正数，表明方程（6）在此区间内单调上升，没有极值和拐点，即随的增加， 也逐渐增大，但相对速率则由于 的不同而异，这种特征也符合一般DV关系。3 Logistic 衍生模型的求解Logistic 衍生模型参数的求解并不很复杂，首先按求解

6、Logistic 模型参数 的方法 4,5K确定 和 ，然后把式（6）还原为下面的形式：VKD（8）bDVKa11这样，即可用对数线性最小二乘法求得参数 。亦即：,（9）lnlln最后，再将式（9）变换为式（6）的形式。但是，笔者已经注意到，上述方法求解模型参数并不是最优的，在可能的情况下，可采用非线性方法最优拟合具有 4 个参数（ ）的 Logistic 衍生模型（式（6） ） ，baKDV,诸如 Marquardt 方法 6、单纯形法 7等，这时，上面采用对数线性化方法求得的参数可用于确定迭代初值。Logistic 衍生模型的非线性最优拟合方法及其有关问题（稳定性、收敛速度和收敛精度等）尚

7、待进一步探讨。不过，从实际应用来看，对于一组 关系数据，既VD使用目测法粗估 和 ，也可获得令人满意的效果，显示出 Logistic 衍生模型较强的数VKD据拟合能力。4 例证研究应用凉水林场落叶松人工林资料 1，分别拟合下列一元材积模型：（10）ba（11）（12）2V（13）D（14）1bKa用剩余平方和 和决定系数 。其中 为总平方Q)(2 SQR)(2和； 为回归平方和）两项指标评价模型的优劣 5。显然， 愈小， 愈大，说)(S 2R明模型的拟合效果愈好（见表 1） 。表 1 凉水林场落叶松人工林资料对不同模型的拟合结果 *一元材积模型 剩余平方和 决定系数)(Q)(2R0.006 0

8、.940594DV015.92.0.00049958 0.995054342980.00056622 0.99439428.0.00046643 0.99538260.00011128 0.99889810.57.06035.1*表中 Logistic 衍生模型的参数 和 的取值系按目测法给出的。VKD从表 1 清晰可见，尽管采用目测法粗估 （0.30）和 （25.0）两参数，但与其它VDK4 个常用的一元材积模型相比，本文提出的 Logistic 衍生模型还是以其剩余平方和 最小)(Q和决定系数（ ）最大而无可争辩地成为最优者。这一结果甚至比采用 Marquardt 方法最2R优拟合 模型获

9、得的结果（ ）还要好得多。因此，baDV 952.0,04826.RQ有理由相信，Logistic 衍生模型有希望成为最受欢迎的一元材积模型。5 结论与讨论以 Logistic 模型为基础，以时间因子为媒介，成功地推导出一个既有理论依据、又有实用价值，能准确描述林木胸径与材积关系的一元材积模型Logistic 衍生模型，实际应用效果良好，可在森林调查中推广应用。Logistic 衍生模型具有较强的数据拟合能力和广泛的适用性。特别是对试验数据不很完整，即在有关调查因子的生长末达平衡位置便终止试验时，该模型依然好用，这是很有意义的。除林木 关系外，Logistic 衍生模型尚可应用于林木 关系以及

10、生物量预VD HD估等更广泛的研究领域。推而广之，任一生物个体或其种群，若有经济意义的 2 个相关因子 的生长节律BA,或发展进程同呈“S”形关系，则可应用本文提出的 Logistic 衍生模型由其中一个易测因子实现对另一因子 的精确估计，反之亦然：BA（15）1bBKa需要说明的是，对于前期生长较慢的树种来说，在林木达到可测胸径之前，已有少量的材积积累，可能会出现 关系中时间因子的非同步现象，对此，Logistic 衍生模型VD未予考虑。但这并不对模型的准确程度构成较大影响。这也正是前面提及的胸径和材积的生长“大体可视为两个同步过程”的基本理由之所在。参 考 文 献：1 孟宪宇主编测树学M

11、北京：中国林业出版社，1996，132-135；193-2002 张连翔，刘学增逻辑斯谛曲线上两个重要特征点的分析及其应用J河北林学院学报，1992，7(2)：154-1583 张连翔，刘学增静电处理林木种子对苗木年高生长节律的影响J西北林学院学报，1991，6(4)：54-584 徐汝梅昆虫种群生态学M 北京：北京师范大学出版社，1987，111-1185 王振中，林孔勋逻辑斯谛曲线 K 值的四点式平均值估计法J生态学报，1987，7(3)：193-1976 张连翔，郑红旗用 Marquardt 方法最优拟合 Gerrand-Chiang 模型J 辽宁林业科技，1996，(2) ：34-35

12、7 董德元，杨 节，苏敏文，等试验研究的数理统计方法M北京：中国计量出版社，1987，309-317Research on the Logistic Extension Model About the Relationship Between the Diameter at Breast Height and Volume of TreesZHANG Lian-Xiang, LIN Yang(Research Institute of Afforestation in Arid Zones of Liaoning Province, Chaoyang, 122000) Abstract: Ba

13、sed on the famous Logistic model and time factor as medium, an one-way tree volume model-Logistic extension model, which can accurately describe the relationship between the diameter at breast height and volume of trees has been deduced. It is a new model both has theoretical basis and practical val

14、ue. It has been found that the Logistic extension model has strong data fitting ability and widespread applicability, except the relationship between the diameter at breast height and volume. It can still be used in the research on the relationship between the diameter at breast height and height and the relationship of plant disease(I-S). It is of both theoretical and practical significance.Key words: Logistic extension model, the relationship between the diameter at breast height and volume of trees, the diameter at breast height and height, one-way tree volume model