计算一棵树树叶的数量以及质量

计算一棵树树叶的数量以及质量，根据叶镶嵌原理，我们将一棵树的树叶分成n层，

并且假设树的轮廓近似为一个圆锥形，如图所示：

圆锥侧面展开之后成为一个扇形，假设每一层的树叶相邻排列，并且每片树叶形状相同大小相等，那么树叶覆盖扇面的方式如下图所示：

，圆锥表面与阳光强度最大时（即正午）直照光相垂直，此时阳光照射方向与水平方向的夹角为（即太阳高度角）。这种关系可以用下图来表示：

θ

假设每一层树叶与圆锥表面平行。这样每一层树叶形成了以Rn为底面半径，hn

Rntan

为

高的圆锥体。假设每一层的底面半径从外向内等差递减，即RnRn1d（d为常数），即

RnR0n

d。那么每层圆锥体除底表面积为SnRn

Rn

2

sin

。假设在

阳光最大位置时，照在外层圆锥除底表面的平均光照强度为Ls0，那么第一层的光照总强度

为LoLs0S0。根据假设光照到树叶上被完全吸收和反射，由于从上一层透过来的日光可以看做一个有一个均匀分布的光源群，根据无影灯原理，可以认为每一层的日光照射都是均匀分布的。从而我们知道由外向内的第二层总光照强度为L1Ls0(S0q0ab)，其中，q0第一层的叶子数；第二层总光照强度为L1Ls0[S0(q0q1)ab]，其中，q1第二层的叶

n1

子数；这样每一层的光照总强度为LnLs0(S0qnab)，每一层的平均光照强度为

n0

n1

Lsn

LnSn

sin[Ls0(S0qnab)]

Rn

n0

2

.那么我们假设只要有Lsn0且第n层的圆锥的体积

应该大于或者等于树木的体积而第n1层的圆锥的体积应该小于或者等于树木的体积。那么树叶就可以生长。

下面计算每一层的树叶数:

将每一层圆锥展开，成为一个扇形，扇形半径为rn

ra



Rnsin

，则每一层的

c

树叶数为qn



c0

2c1

2krna

2



b

c

Rnasin



c0

2c1

2krna

R2

，其中knsin。

rb

Rnasin

nn

n

c

那么一棵树的树叶的总质量为：M

q

n0

mn



n0

c0

2c1

2krna

2

qn

b

R0h0tan，h0hlht其中hl代表树的高度，ht代表树冠底面距离地面的高度（即主

干高）。

模型的关键是得到树叶的层数，这个要考虑体积的原因。

这样就可以得到，树叶的总质量，与树的尺寸特征（高度、质量、体积）是有联系的。女贞叶树叶成卵形或卵状披针形，可近似看作为标准椭圆形，长5～14cm，宽3.5～6cm，

在我们学校中，抽取5棵大致相同的大女贞树进行测量，得到它们的叶子大小、树高、

主干高、胸径的平均值为：长8cm、宽4cm，树高6米，主干2米，胸径40cm的研究树：

假设每一层的树叶的形状大小没有变化，只是质量（与厚度有关）发生了变化，每一片树叶的质量和平均光照强度成函数关系。

在研究不同平均光照强度对女贞树的树叶质量的影响时，我们查找资料，得出了不同光照之下，叶子质量的数据，如下：

光照强度 叶子质量 40 80 160 320 640

树木1树叶质量（g） 0.42 0.53 0.61 0.72 0.68

树木2树叶质量(g) 0.47 0.52 0.59 0.68 0.61

树木3树叶质量(g) 0.68 0.82 0.86 0.96 0.89

树木4树叶质量(g) 0.52 0.62 0.69 0.81 0.74

树木5树叶质量(g) 0.62 0.65 0.72 0.78 0.74

叶子平均质量(g) 0.542 0.628 0.694 0.79 0.732

、则可以拟合出叶子质量与平均光照强度的函数关系：

0.8

0.75

0.7The mass of leaf

0.65

0.6

0.55

0.5

0.450

100200

300400500The average illumination

600700800

y =0.5420 0.6280 0.6940 0.7900 0.7320

p =0.0000 -0.0000 0.0019 0.4786mn0.0031Lsn1.3603

mn0.0019Lsn0.4786

由于在光照强度达到640之后，叶子的质量开始变小，这样不利于研究，我们假设树木生活在合适的光照强度下，那么重新对数据进行拟合，得到叶子质量与平均光照强度的函数关系：

0.8

0.75

0.7The mass of leaf

0.65

0.6

0.55

0.5

0.450

50100

150200250The average illumination

300350400

x=[40 80 160 320];

y=[0.542 0.628 0.694 0.79] n=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.7900

p =0.0000 0.0018 0.4853其中平均光照强度Lsn的单位为molm量m0的单位为g。

mn0.0018Lsn0.4853。

2

s

1

，树叶的质

由于西安位于东经107°40′～109°49′和北纬33°39′～34°45′之间，为了研究简单，我们假设我们学校位于东经108°和北纬3380处。

对于太阳高度角即的确定，首先计算处2月11日时太阳直射点的地理纬度，根据某日太阳直射点的地理纬度位置=23°26′N—（某日—6月22日）*（23°26′*4/365）



'

2326N(TimeJune22)(

'

23264365

'

)N可以计算出2月11日太阳直射点的地理纬度位

置为2326'4908'2572'N，则又由正午太阳高度计算的公式：H＝90°－φ＋δ（H为正午太阳高度，φ为当地纬度，δ为直射点的纬度）。计算出我们学校的正午时太阳高度：

'''9033802572302830

参数的大小如下图所示：

经过matlab编程之后即可计算出研究树的树叶总质量以及总叶子数: 则女贞树的树叶总质量为;

n

M

q

n0

n

mn13.384kg

树的树叶层次： n24

总叶子数为：

24

Q

q

n0

n

2.734410（片）

4