二元材积方程参数优化的研究.pdfVIP

福 建 林 业 科 技 20 ({)：71—75 1993 Jour 0j Fujian Forestry Sci ＆ Toch 二元材积方程参数优化的研究 董希 恩 林文靖 谢健平 阮桂琴 江希钿 — — (宁德地区林业局，福建林学院) —— 7f 一 摘要 车文以38株杉木的拌术资料，从材积残差平方和最小出发t分别用单纯形法和加救最小 二乘法对=元材积方程的参数优化问题进行丁分析研究。结果表明：这二种方{去均能有啦地降低 材积残差平方和 ，优化其参数值，是提高二元材积表精度的有效途径之一 关键词 二元材积方程 兰茎垡丝j单纯形法；苎尘三苎查 用二元材积方程V—aDH 编制二元材积表 ，在求解参数a、b、c时，通常的作法是对材 积方程取对数作线性化处理 ，而后用最小二乘法求解 。然而，按这一方法求解的参数，在某 些情况下 ，将失去使误差平方和达到最小的意义，致使所求的参数并非最佳 为此 ，本文从 材积残差平方和最小出发 ，分别用单纯形法和加权最小二乘法对二元材积方程的参数优化问 题进行 了研究。 1 基础数据 以闽东地区38株杉木的样术资料 每株样木渊胸径、树高、材积，原始数据如表 l所示 用通常的最小二乘法求得二元材积方程为： V 一 0 。0 H 。“。 材积残差平方和 Q一∑ (v一v)=O 2 加权最小二乘法求解参数 设材积与胸径D树高H满足如下关系： V=aD H (1) 现根据 N株样木资料估计方程参数a、b、c时，通常的作法是令新变量 y=logV，x=]ogD， x：~logH，a：10Ba，则 (1)式化为直线形式 ： y=a+bXl+cX2 (2) 然后用最小二乘法即可求得方程参数估计值 a、b、c，但这样求得的a、b、c仅仅是新变量y 的最小二乘逼近，它并不满足原变量材积的残差平方和最小，故所求参数并非最佳。下面我 们将推导出满足原变量材积的最小二乘逼近的a、b、c值。 考查新变量残差项 (y1一 )，在不同的 处 ，相 同的离差转换到原变量上，会产生不同 的原变量 v的离差 ，可以近似认为： ． () 故 5 福 建 林 业 科 技 第 20卷 ’ V一·=·(dV)(y一·) (3) 为此 ，我们用最小二乘法求解 (2)式回归系数对，用 ( )t值对每一新变量残差 (y一 )加 权 ，则残羞平方和为： Q一∑[()l(y．～)]=∑[(dV)-(yj--ax一cX2)) v 设w=·＼／d ~ dy ， ， 根据微积分中求撅值的方法，令詈一o，--~=0，警=0 则有 ： f一一22Cw(Yl--～bXI广cx]=-- ． {署=一2ECw(y--|a--bx一“c：)x-t]一o {警==一2P．Cw．(YJ--a--bx。l—c))‘x=0 faN∑WI+b∑wx‘I1+c∑w -2 ~WtYl {aEW)‘(ti+b∑W，Xl+‘c∑Wxlx=Ew 1YJ (5) a【EW ‘2+b∑WLXX“2+‘c∑wLx=Ew)(t2 V 解方程组即可求得 a、b、