杉木马尾松几个模型的预报有效性的计算机模拟试验

《杉木马尾松几个模型的预报有效性的计算机模拟试验》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、维普资讯http://www.cqvip.comi瞽一ts浙江椿学院学报1993，1O(2)l184~188Journalo，ZhefiangForestryCollege杉木马尾松几个模型的预报有效性的计算机模拟试验林思祖黄青峥杉／吴旺民琳7(7／．。f●摘要用计算机模拟试验的方法，对杉木和马尾松平均树高，平均胸经厦最大密●度的数学模型的预报有效性进行检验。结果表明；各模型的模椒检验精度均达到99·O以上，井且P{>)>85．0．可认为各模型的总体相关系数厦其预报效果是理想的。、／，U／关键词垫堡垫，模拟兰坚数学

2、填型，兰查，曼丝中田分类号$711。在森林经营活动中，许多指标都需要以平均胸径，平均树高及最大密度模型为基础。这些模型不仅是制定产量预测的依据和经营数表及林分收获表的基础，而且是计算机辅助造林设计系统中的重要组成部分。这些数学模型即“理想的模型是来源于“现实系统中的调查观察，而这些模型的可靠性，实用性和精确性及稳定性决定了计算机辅助造林设计系统能否准确可靠稳定地反映“现实系统。一般人们判断数学模型优劣的指标是相关系数的大小，这是正确的，但仅仅视相关系数的大小还是不够的。建立这些模型的资料是从该总体中随机抽取的1个相

3、对于总体单元数来说是很小容量的样本它们是否有足够的代表性和可靠性反映林分生长过程，是值得我们进一步探讨的。在建立数学模型的实践中可以看到来源于同一地区的样本，各个学者用同一数学方程所建立的数学模型的相关系数是有差异的。理论上研究也证实相关系数是有差异的[n】。如果用同一总体的部分资料建立的数学模型的相关系数很大，但重演时相关系数若有很大变化的话，那么应用该模型预澳l同一总体的另一部分资料必定是不够准确的。所以必须对上述几个模型的预报效果的可靠性和稳定性进行统计检验。以往人们的一般做法是把收集样本分成两部分；一部分用

4、来建立模型，另一部分用来检验模型预报效果的有效性。这是困难的。因此，为了克服这种困难，笔者提出利用统计模拟方法用计算机伪造调查现溟l数据中的预报量来检验上述几个模型预报的有效性。啦藕H期t1992—03—15丰建宙林业厅翼助谭时子谋腰之一j术术现在寿宁县林业局工作，书书术现在尤撬县中韩韭靖维普资讯http://www.cqvip.com2期林思祖等：杉术马尾松几十模型的预报有效性的计算机模拟试验185i模型和计算机模拟的原理从总体中抽取若干组实测数据(t，，，，⋯，，{)f=1，2，⋯，．假定预测量满足线性回归方程

5、：．l(x．)+￡i(1)式中，为自变量，玑为预报量，￡是相互独立并且遵从N(0，)分布的随机变量。应用回归分析法得到回归方程的估计式，并可求得相关系数的估计值：，———————～’一㈣’∑(i—)式中；表示回归方程中给出的预报值即理论值，=1罩，是预报量的平均值。为了检验模型(1)和随机预报差异是否显著，由计算机产生遵从N(0，)的伪随机数，从而伪造了几个因变量，计算伪造因变量；的回归方程并求得相关系数R的抽样值”。重复上述Monte-Carlo模型，得到R的平均值，并寻找R的经验分布，从而判断模型的优劣。为了得

6、到正态分布的伪随机数，应用下列近似正态抽样方法抽样p]。假定}(=1，2，⋯，Ⅲ)是相互独立并且遵从N(0，1)上的均匀分布的随机变量，而}一号就是相互独立遵从在(一i1，号)上的均匀分布的随机变量。由(Ⅱ且HyH0B)中心极限定理推知，当m充分大时，随机变量=[(一+(一号)+⋯+(一号)]·~／詈=[(⋯卜刳‘√等(3)是近似遵扶Ⅳ(0，1)，所以⋯={⋯一‘√等(4)是遵从N(o，口)的随机变量，Shreider(1964)：~出Ⅲ=6时，(4)式已足够近似正态分布￡”。此处Ⅲ=6，因此产生均值为0，方差为的

7、正态分布随机变量公式为￡】t‘玎=．～，i亭+言+言a+++芋。一3(5)2数据来源和计算机模拟笔者于1984年在福建省闽北地区8个县市收集了339块样地。样地进行每木检尺，分别径阶测树高，调查大于和等于平均术的林术冠幅，实测2664株树冠，建立了计算机辅助造林设计系统￡。】。其中平均高，平均胸径及最大密度模型分别为}口=+日优(§)维普资讯http://www.cqvip.com浙江林学院学报l0卷D=口十坍(7)Ⅳ皿“=口十b／D(8)裹1树种平均高数学接受·Table1育一日dm。de1s。fCunni”。i

8、ⅡlanceoandPnd∞¨fdnd●裹3树种最大密度敛学横型”Table3Maximumdensity-瓦modelsofCunninghamiaf口月0fD，dandPin，rE站0ni。”D树种最大舒度模型相关系数回归剩采标准差杉木Ⅳm二；2．6892+234B．995／DO．09098O．67i4马尾橙Jvmx=120997+1410．9zs／D-