森林面积资源的计算机仿真及其应用*张冰峰****（北京林业大学信息学院，100083）摘要本文论述了森林面积资源的计算机仿真理论、方法及其应用。森林面积资源仿真是进行森林面积资源现代化管理的理论依据和重要的系统方法之一，具有广泛的应用。关键词森林面积资源计算机仿真应用1引言仿真技术已成为研究各种复杂系统的重要工具，在生态环境和林业上的应用已成为研究的热点之一。森林面积资源的计算机仿真就是对描述森林面积资源的数学模型在计算机上进行各种模拟、试验和分析，以便描述森林面积资源的动态变化规律和运动轨迹。目前国内外林业上进行的森林和森林资源的计算机仿真所使用的数学模型不是统计模型就是经验模型，并且多为常微分方程，这些模型都不是根据森林资源的动态变化规律和生灭过程而推导出来的数理方程，因而它们不能精确的对森林资源系统进行描述。本文所使用的模型是根据森林面积资源的动态变化规律和生灭过程而推导出来的数理方程，是偏微分方程，它能对森林面积资源的动态变化规律和生灭过程做出精确的表达，能建立科学的森林面积资源的仿真理论。它无论是在理论上，还是在实际应用中都有重要的意义。2森林面积资源仿真的数学模型2.1仿真的连续型数学模型描述森林面积资源的方程一般有：2.1.1同龄纯林的林龄面积分布结构变化方程[1]：P(τ,t)：表示同龄纯林面积的林龄分布密度函数μ(τ,t)：表示同龄纯林的相对采消率函数*本文为国家自然科学基金资助项目。****张冰峰,硕士研究生,主要研究方向:森林资源系统管理的仿真技术{∂P(τ,t)∂t+¿∂P(τ,t)∂τl=−μ(τ,t)P(τ,t)−g(τ,t)¿{P(τ,0)=P0(τ)¿¿¿¿l:表示单位时间内林龄的变化率，即林龄与时间的换算量P(τ,0)：表示方程的初始条件φ(t)：表示t时刻，单位时间内造林更新并郁闭成林的面积β(t):表示t时刻单位时间内的造林更新率γ(t):表示t时刻单位时间内的郁闭成林率g(τ,t)：为干扰项2.1.2面积按树种混交系数·平均年龄分布结构变化方程[2]：A：表示林分的树种组成系数或比例τ：表示平均年龄P(A,τ,t)：表示t时刻树种混交的隶属度为A,平均年龄数为τ的分布密度函数连续型模型在理论分析中是比较方便的，但不适合实际计算和统计分析。为了与现实中采集的数据相结合和实现森林面积资源数学模型的计算机仿真，须对模型离散化。2.2离散型模型为了便于讨论，下面以同龄纯林的林龄面积分布结构变化方程为例，给出该模型离散化后以向量形式表示的代数方程和以矩阵形式表示的矩阵方程。（1）式离散化后的向量代数方程形式为[1]：{S(k+1)=H(k)S(k)+γ(k)β(k)Μ(k)S(k)−G(k)¿¿¿¿（3）S(k+1):表示k+1分期同龄纯林面积向量H(k):表示面积的林龄转移矩阵S(k):表示k分期同龄纯林面积向量M(k):表示消耗率矩阵G(k):表示同龄纯林因突然大面积的火灾、滥砍乱伐等干扰因素使森林减少的面积向量S(0):表示初始条件向量{∂P(A,τ,t)∂t+¿∂∂A[^ψ1(A,τ,t)P(A,τ,t)]+∂∂τ[^ψ2(A,τ,t)P(A,τ,t)]=−μ(A,τ,t)P(A,τ,t)¿{¿{P(A,τ,0)=P0(A,τ)¿¿¿¿（1）式离散化后的矩阵方程为[1]：[S1(k+1)¿][S2(k+1)¿][...¿]¿¿¿¿=[00...001−μ1(k)0...0001−μ2(k)...00.......0001−μm−1(k)0].[S1(k)¿][S2(k)¿][....¿]¿¿¿¿+γ(k)β(k)[μ1(k)...μm(k)0...0][S1(k)¿][S2(k)¿][....¿]¿¿¿¿-[g1(k)¿][g2(k)¿][...¿]¿¿¿¿（4）M(k)：为消耗率矩阵H(k)：为林龄面积转移矩阵Si(k+1)：为k+1分期第i龄级的面积Si(k)：为k分期第i龄级的面积β(k)：为k分期的造林更新率γ(k)：为k分期的郁闭成林率gi(k)：为k分期第i龄级的同龄纯林因突然的干扰因素使森林减少的面积离散化后上述矩阵方程可视为是同龄纯林面积资源仿真的仿真模型，便可根据分布参数的仿真语言编制仿真软件，从而在计算机上实现仿真计算。3森林面积资源仿真的计算设计3．1仿真的辅助计算设计要在计算机上实现森林面积资源的仿真计算，需要许多参数和数据，而这些参数和数据必须预先从许多图面资料、调查数据、统计资料经过提取、统计、分析或外推而得到。有了这些准确的参数和数据才能进行准确的计算机仿真。而为了获取这些参数和...
VIP
