1、 计算一棵树树叶的数量以及质量，根据叶镶嵌原理，我们将一棵树的树叶分成 n 层，并且假设树的轮廓近似为一个圆锥形，如图所示：圆锥侧面展开之后成为一个扇形，假设每一层的树叶相邻排列，并且每片树叶形状相同大小相等，那么树叶覆盖扇面的方式如下图所示：，圆锥表面与阳光强度最大时（即正午）直照光相垂直，此时阳光照射方向与水平方向的夹角为 （即太阳高度角） 。这种关系可以用下图来表示：假设每一层树叶与圆锥表面平行。这样每一层树叶形成了以 为底面半径， 为nRtanRh高的圆锥体。假设每一层的底面半径从外向内等差递减，即 （d 为常数） ，即1。那么每层圆锥体除底表面积为 。假设在0nRd22sinnnnS
2、h阳光最大位置时，照在外层圆锥除底表面的平均光照强度为 ，那么第一层的光照总强0sL度为 。根据假设光照到树叶上被完全吸收和反射，由于从上一层透过来的日光0osLS可以看做一个有一个均匀分布的光源群，根据无影灯原理，可以认为每一层的日光照射都是均匀分布的。从而我们知道由外向内的第二层总光照强度为 ，其100()sLSqab中， 第一层的叶子数；第二层总光照强度为 ，其中， 第0q 10()sSq1二层的叶子数；这样每一层的光照总强度为 ，每一层的平均光照100nnsLab强度为.那么我们假设只要有 且第 层的圆锥的体积1002sin()nssnnLSqabLSR0snL应该大于或者等于树木的体
3、积而第 层的圆锥的体积应该小于或者等于树木的体积。那么树叶就可以生长。下面计算每一层的树叶数:将每一层圆锥展开，成为一个扇形，扇形半径为 ，则每一层2sinnRrh的树叶数为 ，其中 。sin0 0211Rrc cn na an ckkaqbbsinkr那么一棵树的树叶的总质量为：sin0021Rcnna nckraMqmqb ， 其中 代表树的高度， 代表树冠底面距离地面的高度（即主0tanRh0lthl th干高） 。模型的关键是得到树叶的层数，这个要考虑体积的原因。这样就可以得到，树叶的总质量，与树的尺寸特征（高度、质量、体积）是有联系的。)女贞叶树叶成卵形或卵状披针形，可近似看作为标准
4、椭圆形，长 514cm，宽 3.56cm，在我们学校中，抽取 5 棵大致相同的大女贞树进行测量，得到它们的叶子大小、树高、主干高、胸径的平均值为：长 8cm、宽 4cm，树高 6 米，主干 2 米，胸径 40cm 的研究树：假设每一层的树叶的形状大小没有变化，只是质量（与厚度有关）发生了变化，每一片树叶的质量和平均光照强度成函数关系。在研究不同平均光照强度对女贞树的树叶质量的影响时，我们查找资料，得出了不同光照之下，叶子质量的数据，如下：光照强度叶子质量树木1树叶质量（g）树木2树叶质量(g)树木3树叶质量(g)树木4树叶质量(g)树木5树叶质量(g)叶子平均质量(g)40 0.42 0.47
5、 0.68 0.52 0.62 0.54280 0.53 0.52 0.82 0.62 0.65 0.628160 0.61 0.59 0.86 0.69 0.72 0.694320 0.72 0.68 0.96 0.81 0.78 0.79640 0.68 0.61 0.89 0.74 0.74 0.732、则可以拟合出叶子质量与平均光照强度的函数关系：0 100 200 300 400 500 600 700 8000.450.50.550.60.650.70.750.8The average illumination Themass of leafRaw dataThe fitting
6、curvey =0.5420 0.6280 0.6940 0.7900 0.7320p =0.0000 -0.0000 0.0019 0.4786 0.31.60nsnL0.19.4786nsnmL由于在光照强度达到 640 之后，叶子的质量开始变小，这样不利于研究，我们假设树木生活在合适的光照强度下，那么重新对数据进行拟合，得到叶子质量与平均光照强度的函数关系：0 50 100 150 200 250 300 350 4000.450.50.550.60.650.70.750.8The average illumination Themass of leafRaw dataThe fitti
7、ng curvex=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 quality )y =.5420 0.6280 0.6940 0.7900p =0.0000 0.0018 0.4853 其中平均光照强度 的单位为 ，树叶的质snL21
8、mols量 的单位为 。0mg。.180.453nsnL由于西安位于东经 1074010949和北纬 33393445之间，为了研究简单，我们假设我们学校位于东经 108和北纬 处。80对于太阳高度角即 的确定，首先计算处 2 月 11 日时太阳直射点的地理纬度，根据某日太阳直射点的地理纬度位置=2326N（某日6 月 22 日）*（2326*4/365）可以计算出 2 月 11 日太阳直射点的地理纬度 2364236()()5NTimeJunN位置为 N，则又由正午太阳高度计算的公式：H90（H 49087为正午太阳高度， 为当地纬度， 为直射点的纬度） 。计算出我们学校的正午时太阳高度： 32530参数的大小如下图所示：参数 ( )(m)a（m）b(cm)d(m)lh(m)0(0sL21mol)数值 30 0.08 0.04 0.1 6 2 320经过 matlab 编程之后即可计算出研究树的树叶总质量以及总叶子数:则女贞树的树叶总质量为; 0nMqm13.84kg树的树叶层次： 24总叶子数为：（片）40.7310nQq
