1、 计算一棵树树叶的数量以及质量，根据叶镶嵌原理，我们将一棵树的树叶分成n层，并且假设树的轮廓近似为一个圆锥形，如图所示：圆锥侧面展开之后成为一个扇形，假设每一层的树叶相邻排列，并且每片树叶形状相同大小相等，那么树叶覆盖扇面的方式如下图所示：，圆锥表面与阳光强度最大时（即正午）直照光相垂直，此时阳光照射方向与水平方向的夹角为（即太阳高度角）。这种关系可以用下图来表示：假设每一层树叶与圆锥表面平行。这样每一层树叶形成了以为底面半径，为高的圆锥体。假设每一层的底面半径从外向内等差递减，即（d为常数），即。那么每层圆锥体除底表面积为。假设在阳光最大位置时，照在外层圆锥除底表面的平均光照强度为，那么第一
2、层的光照总强度为。根据假设光照到树叶上被完全吸收和反射，由于从上一层透过来的日光可以看做一个有一个均匀分布的光源群，根据无影灯原理，可以认为每一层的日光照射都是均匀分布的。从而我们知道由外向内的第二层总光照强度为，其中，第一层的叶子数；第二层总光照强度为，其中，第二层的叶子数；这样每一层的光照总强度为，每一层的平均光照强度为.那么我们假设只要有且第层的圆锥的体积应该大于或者等于树木的体积而第层的圆锥的体积应该小于或者等于树木的体积。那么树叶就可以生长。下面计算每一层的树叶数:将每一层圆锥展开，成为一个扇形，扇形半径为，则每一层的树叶数为，其中。那么一棵树的树叶的总质量为：，其中代表树的高度，代
3、表树冠底面距离地面的高度（即主干高）。 模型的关键是得到树叶的层数，这个要考虑体积的原因。 这样就可以得到，树叶的总质量，与树的尺寸特征（高度、质量、体积）是有联系的。)女贞叶树叶成卵形或卵状披针形，可近似看作为标准椭圆形，长514cm，宽3.56cm，在我们学校中，抽取5棵大致相同的大女贞树进行测量，得到它们的叶子大小、树高、主干高、胸径的平均值为：长8cm、宽4cm，树高6米，主干2米，胸径40cm的研究树：假设每一层的树叶的形状大小没有变化，只是质量（与厚度有关）发生了变化，每一片树叶的质量和平均光照强度成函数关系。在研究不同平均光照强度对女贞树的树叶质量的影响时，我们查找资料，得出了不
4、同光照之下，叶子质量的数据，如下：光照强度叶子质量树木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、则可以拟合出叶子质量与平均光照强度的函数关系：y =0.5420 0.6280 0.6940 0.7900 0.7320p =0.0000 -0.0000
5、0.0019 0.4786由于在光照强度达到640之后，叶子的质量开始变小，这样不利于研究，我们假设树木生活在合适的光照强度下，那么重新对数据进行拟合，得到叶子质量与平均光照强度的函数关系：x=40 80 160 320;y=0.542 0.628 0.694 0.79n=2;p=polyfit(x,y,n)xi=linspace(0,400,10000);z=polyval(p,xi); %多项式求值plot(x,y,o:,xi,z,r,x,y, b )legend(Raw data,The fitting curve )xlabel(Light intensity )ylabel(Leaf
6、 quality )y =.5420 0.6280 0.6940 0.7900p =0.0000 0.0018 0.4853其中平均光照强度的单位为，树叶的质量的单位为。 由于西安位于东经1074010949和北纬33393445之间，为了研究简单，我们假设我们学校位于东经108和北纬处。对于太阳高度角即的确定，首先计算处2月11日时太阳直射点的地理纬度，根据某日太阳直射点的地理纬度位置=2326N（某日6月22日）*（2326*4/365）可以计算出2月11日太阳直射点的地理纬度位置为N，则又由正午太阳高度计算的公式：H90（H为正午太阳高度，为当地纬度，为直射点的纬度）。计算出我们学校的正午时太阳高度：参数的大小如下图所示：参数()(m)（m）(cm)(m)(m)()数值300.080.040.162320经过matlab编程之后即可计算出研究树的树叶总质量以及总叶子数:则女贞树的树叶总质量为;树的树叶层次：总叶子数为：（片）
