n层,并且假设树的轮廓近似为一个圆锥形,如图所示:计算一棵树树叶的数量以与质量,根据叶镶嵌原理,我们将一棵树的树叶分成圆锥侧面展开之后成为一个扇形,假设每一层的树叶相邻排列,并且每片树叶形状相同大小相等,那么树叶覆盖扇面的方式如下图所示:,圆锥表面与阳光强度最大时〔即正午〕直照光相垂直,此时阳光照射方向与水平方向的夹假设每一层树叶与圆锥表面平行这样每一层树叶形成了以R为底面半径,hn—为tan高的圆锥体假设每一层的底面半径从外向内等差递减,即RnRn1d〔d为常数〕,即RnRond那么每层圆锥体除底表面积为SnRn,h2Rn2假设在sin阳光最大位置时,照在外层圆锥除底表面的平均光照强度为Lso,那么第一层的光照总强度为LoLsoSo根据假设光照到树叶上被完全吸收和反射,由于从上一层透过来的日光可以看做一个有一个均匀分布的光源群,根据无影灯原理,可以认为每一层的日光照射都是均匀分布的从而我们知道由外向内的第二层总光照强度为LiLso(Soqoab),其中,q0第一层的叶子数;第二层总光照强度为L1Lso[So(q0q1)ab],其中,q第二层的叶子数;这样每一层的光照总强度为l-nLs0(S01qnab),每一层的平均光照强度为0Lsin[Lso(So§n1qnn0Rn"ab)]——.那么我们假设只要有Lsn0且第n层的圆锥的体积应该大于或者等于树木的体积而第么树叶就可以生长。
下面计算每一层的树叶数:1层的圆锥的体积应该小于或者等于树木的体积那将每一层圆锥展开,成为一个扇形,扇形半径为Rnsin,则每一层的树叶数为qnr2c1c2krnaa2basin2krn2c1a2b,其中$sinr那么一棵树的树叶的总质量为:q.mnn0Rohotan,h°hht其中h代表树的高度,Rasin2krn2c1a2bqnht代表树冠底面距离地面的高度〔即主咼〕模型的关键是得到树叶的层数,这个要考虑体积的原因这样就可以得到,树叶的总质量,与树的尺寸特征〔高度、质量、体积〕是有联系的女贞叶树叶成卵形或卵状披针形,可近似看作为标准椭圆形,长5〜14cm,宽3.5〜6cm,在我们学校中,抽取5棵大致相同的大女贞树进行测量,得到它们的叶子大小、树高、主干高、胸径的平均值为:长8cm、宽4cm,树高6米,主干2米,胸径40cm的研究树:假设每一层的树叶的形状大小没有变化,只是质量〔与厚度有关〕发生了变化,每一片树叶的质量和平均光照强度成函数关系在研究不同平均光照强度对女贞树的树叶质量的影响时,我们查找资料,得出了不同光照之下,叶子质量的数据,如下:度叶子质量树木1树叶质量〔g〕树木2树叶质量(g)树木3树叶质量(g)树木4树叶质量(g)树木5树叶质量(g)叶子平均质量(g)400.420.470.680.520.620.542800.530.520.820.620.650.6281600.610.590.860.690.720.6943200.720.680.960.810.780.796400.680.610.890.740.740.732、则可以拟合出叶子质量与平均光照强度的函数关系:faprfoQwsameh0.50.450100200300400500600700800Theaverageiliumination0.4786mn0.0031Lsn1.3603mn0.0019Lsn0.4786由于在光照强度达到640之后,叶子的质量开始变小,这样不利于研究,我们假设树木生活在合适的光照强度下,那么重新对数据进行拟合,得到叶子质量与平均光照强度的函数关系:7565aGa5oofaprfoQwsameo—RawdataThefittingcurveTheaverageiliuminationx=[4080160320];y=[0.5420.6280.6940.79]n=2;p=polyfit(x,y,n)xi=linspace(0,400,10000);z=polyval(p,xi);%多项式求值其中平均光照强度Lsn的单位为molms,树叶的质mn0.0018Lsn0.4853。
由于XX位于东经107°4~'109°4和北纬33°3434°4之间,为了研究简单,我们假设我们学校位于东经108°和北纬3380'处对于太阳高度角即的确定,首先计算处2月11日时太阳直射点的地理纬度,根据某日太阳直射点的地理纬度位置=23°26'N—■〔某日一6月22日〕*〔23°26'*4/36〕<■'232642326N(TimeJune22)()N可以计算出2月11日太阳直射点的地理纬度位365置为2326'49:08'2572'N,则又由正午太阳高度计算的公式:H=90—时HH为正午太阳高度,0为当地纬度,S为直射点的纬度〕计算出我们学校的正午时太阳高度:■-'••'■■-'■■90・3380257230^2830”参数的大小如下图所示:参数()a(m)b〔m〕d(cm)hl(m)h)(m)Lse(molm2s1)数值300.080.040.162320经过matlab编程之后即可计算出研究树的树叶总质量以与总叶子数则女贞树的树叶总质量为;nMqng13.384kgn0树的树叶层次:n24总叶子数为:24Qqn2.7344104〔片〕n0。