林木生长规律的数学模型模拟.docxVIP

林木生长规律的数学模型模拟 1 辽宁内的其他县 本根直径材料的数据基本上覆盖了辽宁省的各个地区，树种全面。主要数据采集地区:辽宁东部的抚顺县、新宾县、顺城区、本溪县、桓仁县、宽甸县、西丰县、岫岩县;辽宁南部的金州区、普兰店、旅顺口、庄河县;辽宁西部的北票县、凌海县、北宁市、连山区、绥中县、兴城市、朝阳县。共31 859个样本,其中红松采集了3 179组数据,落叶松采集了5 359组数据,油松采集了1 282组数据,刺槐采集了4 184组数据,柞树采集了5 235组数据,杨柳采集了4 372组数据,硬阔采集了5 146组数据,软阔采集了3 102组数据。 2 处理数据 2.1 相同根径值的样本胸径值计算 在收集的样本数据中,剔除数据异常值,对相同根径值的样本胸径值取平均数,即:胸径值=∑样本胸径该径阶样本数量=∑样本胸径该径阶样本数量。 2.2 数据分类 按各树种组根径值由小到大次序排列,形成根径—胸径有序数据对,以便建立数学模型。 3 各树种组相关指标的比较 3.1各树种组样本数据散点图绘制(略) 3.2模拟方程类型的选择与建立 通过对各树种组胸径—根径散点图分析,为各树种组建立了5种模拟方程,包括:线性模拟方程、对数模拟方程、三阶二项式模拟方程、指数模拟方程和乘幂模拟方程。各树种组模拟方程如下: 线性模拟方程:Y=a·X-b; 对数模拟方程:Y=a·Ln(X)-b; 乘幂模拟方程:Y=a·Xb; 指数型模拟方程:Y=a·ebx; 三阶二项式模拟方程:Y=a·X3+b·X2+ c·X+d; 其中:Y—胸径模拟值,X—根径值,a,b,c,d—参数。 3.3模拟方程优化筛选 各树种组相关指标比较以红松为例,见表1,其它树种组略去。 根据上述各模拟方程指标比较,三阶二项式模拟方程各项指标最优。因此,各树种组最优模拟方程选定为: 1)红松树种组: Y=-03+02+ 0+62)落叶松树种组: Y=-03+02+ 0+03)油松树种组: Y=-03+02+ 0+24)刺槐组: Y=-03+02+ 0+05)硬阔组: Y=03-02+ 0+06)软阔组: Y=-03+02+ 0+17)柞树组: Y=03-02+ 0-08)杨柳组: Y=-03+02+ 0-03.4单株材积计算 采用辽宁省地方标准DB21-778-814-94林业经营数表一元立木材积式: V=a1×[f(D)]a2×{b1+b2/[f(D)+k]}a3 其中:f(D)=c1+c2D围 计算立木积,材积式中的参数a1,a2,a3,b1,b2,k,c1,c2已设定,具体数值略去。 4 瓦房店采集样本采集 检验数据采集于辽宁省东部地区的本溪县、新宾县、抚顺县、辽西地区的北票县、中部地区的新民县、辽南地区的瓦房店,共采集样本1 029株,其中:辽西杨柳130株、辽西刺槐27株、辽西油松48株、辽东刺槐71株、辽东红松110株、辽东柞树144株、辽东硬阔121株、辽东软阔77株、辽东落叶松72株、辽东油松116株、辽南刺槐113株。检验样本占原始数据总量的55.35%。 4.1 计算径生长空间 Kd(i)=(D0?D1)/D0×100%Kv(i)=(V0?V1)/V0×100%Kdˉˉˉˉˉ=∑i=1nKd(i);Kvˉˉˉˉˉ=∑i=1nKv(i)Κd(i)=(D0-D1)/D0×100%Κv(i)=(V0-V1)/V0×100%Κdˉ=∑i=1nΚd(i);Κvˉ=∑i=1nΚv(i) 式中:D0—胸径观测值;D1—胸径模拟值;Kd(i)—模拟胸径精度;Kv(i)—模拟材积精度;V0—观测值一元材积表材积值;V1—模拟材积值;KdˉˉˉˉˉΚdˉ—胸径精度平均数;KvˉˉˉˉˉΚvˉ—材积精度平均数;n—观测值样本数量;I—第i个样本。 4.2 材积表材积与模拟材积绝对误差计数模型试验 P(i)=V0(i)-V1(i);S=∑I=1NP(i)S=∑Ι=1ΝΡ(i) 式中:P(i)—一元材积表材积与模拟材积的绝对误差;S—一元材积表材积与模拟材积绝对误