实验五一元材积表的编制.docVIP

实验五 一元材积表的编制 一、实验目的 了解编制材积表所需资料及收集方法；掌握数式法编制一元材积表的过程和方法。 二、实验内容和方法 （一）一元材积表的直接编制 1、编表资料的收集 编制一元材积表需在材积表的使用地区范围内随机抽取样木，伐倒后用中央断面区分求积式测定其材积，并精确测定其胸径，取得V-D的成对值。样木不少于100株，并要求分布在各种立地条件，各径阶都要有一定的数量，否则要加大样本数，以保证编表资料能反映材积表使用地区的材积平均水平。 各人利用老师提供的资料进行编制，每人的数据都不相同。 编表资料的整理 将原始资料按2cm径阶分组统计，求出各径阶的算术平均胸径和平均材积，统计于表5-1。 表5-1一元材积表原始资料汇总表 径阶 平均胸径 平均材积 株数 2 2.85 0.00165 4 4 4.20 0.004318 44 6 5.96 0.008219 43 8 7.85 0.01868 25 10 10.17 0.028187 23 12 12.00 0.0402 9 14 14.24 0.09032 5 16 16.23 0.15345 4 18 17.60 0.1767 4 20 19.75 0.237125 4 22 21.95 0.29455 2 24 24.55 0.4412 2  绘制描绘直径－材积相关曲线 以横坐标表示直径，以纵坐标表示材积，根据各编表样木的胸径与材积，在计算机上以胸径为横坐标，材积为纵坐标作散点图。 描绘平均直径－平均材积相关曲线 绘制描绘平均直径－平均材积相关曲线：以横坐标表示各个径阶的平均直径，以纵坐标表示各个径阶的平均材积根据各编表样木的胸径与材积，在计算机上以径阶为横坐标，材积为纵坐标作折线散点图。 选择回归模型 根据曲线趋势选择最能模拟此曲线的回归方程，这次实习我们统一用。 方程参数求解 方程是非线性回归方程，在进行方程参数求解时，可将其直线化。对两边取对数，可得,这样可化为()，即非线性回归方程线性化。以下采用一般最小二乘法求解。 设原始资料的总样本单元数为n，第i株树的直径和和材积分别为（Di，Vi），取对数后有（xi，yi）。用Excel求出： ，，，，，， 再求离差平方和： 则 ; 这样可求得回归方程中的a，b两个参数。求得幂函数方程后，将实际胸径代入，可得该株树的理论材积，与原始数据中各实际材积相比较，可求得回归剩余离差平方和。  数据计算结果如下表： 参 数 n 值 169 143.6865507 -315.6559466 130.264614 641.1612267 -248.570128 0.850216276 -1.867786666 由上表并结合公式可求得下表： 参数 a 值 8.099969916 51192.445156684 -3.946698677 0.000113058 求得冥函数方程：V=0.000113058D2.445156684 方程检验 标准差，n为样本数，m为自变量的个数 标准误差 注意这里的不是线性方程中的的平均值，而是幂函数的V，因为我们最终是对幂函数进行拟合，而不是对线性方程拟合，对线性方程的拟合只是拟合幂函数的一个步骤和手段。 相对误差，为平均材积实际值。 t0.05是自由度为n－1的t分布的0.05双侧分位数。 精度P＝1-E 注意Q要用幂函数的回归残差。而不是线性化后的Q 相关指数＝ 参数 Q S E P 相关指数 值 0 0.027070053 0.002082312 11.09% 88.91% 87.46% 利用SPSS拟合材积方程 在SPSS中，拟合方程，求出参数值与统计量，并与你计算的参数值和相应统计量进行比较。除此之外，再在课本129页表5-5中任选一方程在SPSS中拟合，也求出参数和统计量。 参数 方程 a b Residual a.R squared V=aDb 0.0000365 2.937 0.023 0.974 V=a+bd2 -0.019 0.001 0.058 0.934 由上表可得方程V=aDb较符合该胸径与材积的关系。 列出一元材积表 将各径阶中值代入最优立木材积式，即求出对应径阶的材积，即为各径阶相应单株平均材积，将其列成表即为一元材积表。 径阶 材积 2 0.00028 4 0.002141 6 0.007042 8 0.016393 10 0.031571 12 0.053932 14 0.084815 16 0.125543 18 0.17743 20 0.241779