Richards增长曲线的参数估计.pdf

第31卷 哈尔滨师范大学自然科学学报 NATURALSCIENCES OF 第5期 JOURNALHARBINNORMALUNIVERSITY Richards增长曲线的参数估计木 乔钰，徐文科一 (东北林业大学) 【摘要】介绍了Richards增长曲线方程，而后对Richards增长曲线方程的参 数进行研究．运用线性最小二乘法、双向差分误差平方和最小估计法和双向差分加 权最小二乘估计法对其参数进行估计，从定量的角度研究了Richards增长曲线方 程，为进一步的研究奠定了基础，最后将模型应用到传染病疫情实例中． 【关键词】Richards增长曲线；参数估计；最小二乘法；双向差分法 中图分类号：0241 文献标识码：A 文章编号：1000—5617(2015)05—0023—04 线拐点的位置．该增长曲线方程的图形是拐点随 1 Richards增长曲线方程 参数取值的变化而变化的“S”曲线‘4|．该文利用 Richards模型的微分形式‘13为 线性最小二乘估计法‘51，双向差分误差平方和最 ，，、 小估计法M1和双向差分加权最／J、-乘法口。81估 警=叼矿叫 一 计Richards增长曲线方程的参数． 这里警表示生物增长率，叼，y和m是异速生长常 2 线性最小二乘法估计Richards 数．(1)式为Bernoulli型微分方程，积分可得通 增长曲线参数 解，得到Richards增长曲线方程．在文献[1]中 令m一1=6，并将(3)式中的因变量K用 初值t=0，V=Vo，在文献[2]中选取初始条件 y(t)表示，(3)式可表为 为t=t。，K。=0，更一般的情况是在文献[3]中 选取初值t=t。，V=vo，得到通解结构 y(￡)=A(1一Be一缸)一百 (4) 对(4)式取对数求导整理可得 E=[旦一(卫一圪一“)e州1训(…。’]击(2) ’y ‘y (5) 方程(2)可以简写为 掣=(一争’e-XyS+1(t) 令一BA一=C，并对(5)式取对数，可得 K=A(1一Be啦)而 (3) T，l—m 这里A(卫)志，B=[1一号]e‰，K=7(1一 jY 1177 (6) m)．通常把(3)式叫做Richards增长曲线方程， 这里A，曰，K，m都是待估计参数，A是生物生长的 极限值，B是与生长初始条件相关的参数，K是描 譬≯，并令风=lniKc，卢。=6+1，卢：=一K设 述生长速率的参数，m决定着曲线的形状以及曲 误差为e(f)，所以由(6)式可得 收稿日期：2015—05—11 +黑龙江省自然科学基金项目(C201408) }}通讯作者：xuwenke2004@sina．coin 万方数据 24