拟合闽北杉木干形曲线的三次样条函数的节点选择-福建农林大学图书馆.pdfVIP

维普资讯 浙 江 林 学 院 学 报 2003，20 2 ：219 223 Joun~／ofZhejiangForestryCollege 文章编号 ：1000．5692 2003 02．0219．05 拟合 闽北杉木干形 曲线的三次样条函数 的节点选择 刘 健 ，郭建宏 ，郭进辉 ，郑淑娟 福建农林大学 林学院，福建 南平 353001 摘要 ：收集闺北杉木解析木材料 ，运用统计学的知识并结合计算机技术，用三次样条函数拟 合杉木干形 曲线 ，以样本相关 系数均值为主要评价指标 ，筛选 出拟合杉木干形 曲线的三次样 条函数的节点 ，即从树干基部算起 ，依 次为树 高 的 0，l0％，l5％一20％，50％一6o％， 80％一85％和 100％处。经3处方案的拟合 比较和实例验证 ，结果表明所选节点的似合效果 明显优 于用传统经验节点的拟合效果；所选节点整体拟合精度 高，适应性好 ，能够全面确切 反映杉木完整形状。表 4参8 关键词 ：样条 函数 ；杉木 ；干形 曲线；节点选择 中图分类号 ：$758 文献标识码 ：A 树干纵断面轮廓的对称 曲线称为干形曲线。该 曲线 自基部 向梢端变化大体呈凹曲线 、平行于 轴 的直线 、抛物线和相交于 轴 的直线这4种类型。为了探索能够表达整株树干的形状 ，大偶真一和 门捷列夫等提出过许多干形曲线式…。这些干形 曲线式一般只能表示出中上部形状 ，对基部适应性较 差 。样条 函数是光滑连接 的分段多项式，具有很强的适应数据微小变化 的能力 ，并具有整体光滑性， 因而成为曲线修匀的理想工具 】。美 国NCHS标准就是采用三次样条函数修匀 的∞。近年来 ，样条函 数在医学和物理学等领域得到广泛的研究和应用D 】，在林学和生态学方面 的研究还不多见 。薛亚琴 等 根据传统的经验数据作为节点，用样条函数拟合林木干形曲线，取得了较好的拟合效果 。因受 计算技术的限制，早期的研究没有进行节点的选择，且随着树种、气候等的差异，林木的干 曲线变化 较大n 。为此 ，本 研 究 运用 统 计学 的知识 ，结 合 计 算机 技 术 ，用 三 次样 条 函数 拟 合 杉木 Cunn／ngham／aa／nceoa／ta干形曲线 ，筛选 出整体拟合精度高，适应性好 ，能够全面确切反映杉木完整形 状 的三次样条 函数节点，并进行拟合计算及精度 比较 。 l 研 究方法 1．1 三次样条 函数简介H 所谓样条 函数 ，从数学上说 ，就是按一定光滑性要求 “对接”起来 的分段多项式 。具体地说 ，在 区间 [口，b]上，给定一个分划 ：口 。 ： ，… b，这样我们在区间 [ ，Ⅲ]上可以找 到近似 的三次多项式 P i 1，2，…，n一1 ，要求 ： P f ，Pi +。 f +。 ， 1，2，…，n一1 ， 1 且 P Pi一。 ，P P．。 ， 2；3，…，n—1 ， 2 收稿 日期 ：2002．09．28；修 回日期 ：20173．02．23 基金项 目：福建省教育厅科学基金资助项 目 JB0l165 作者简介 ：刘健 1963一 ，男，福建福州人 ，副教授 ，从事森林资源资产评估和森林可持续经营等研究 。 维普资讯 浙 江 林 学 院 学 报 2003年 6月 则称这样 的分段多项式函数为样条 函数 ，点 称为样条函数 的节点。为使表达式简单 ，记 ：h Ⅲ 一 ， P 。由定义知，每一个 P 都是 的线性函数，因此在区间 [ ，…] 上可 由区间端点的值 和 +。用线性插值按 Langrange公式直接写 出：即：P + 。 将上式积分二次，利用条件 1 得 出： fbi 高 扎。一+等 [每一 +[鲁一警 …一 3 式 中：i 1，2，…，n一1 。由 3 式可以看 出，只要求 出每个 待估参数 ，就可 由 3 式求 出所有的三次插值多项式。通常采用追赶法求解 ，树干顶部边界条件采用 自由边界，根径处边界采用 固定边界，一般取 一2．5左右为宜。 1．2 节点选择方法 用样条函数方法模拟生物体形状 ，应遵循客观规律 ，了解其生物背景 。根据林木干形 的生物学特 性，用样条函数拟合林木干形 曲线时，一般取 6个节点为宜 包括根径处和梢顶处 。用杉木解析材 料 ，从树干基部算起 ，把树高按 5％间隔区分 ，计算各相对高处断面高 i 0，1，2，…，20 ，利用插 值法求各断面高处的断面半径 Y i 0，1，2，…，20 ，然后以组合的方式 ，选取 6个点 P ，Y 作为三 次样条 函数一组节点，6个点 中要求必须选根径处 P 。，Y。 和梢顶处 P 加，Y加 。这样 的组合共有 c：