新能源开发之麻风树问题10

新能源开发问题摘要本文以麻疯树施肥立意，以判断施肥对样地有无明显改善效果为目标建立模型。对于问题1，我们从非萌生枝长、非萌生枝径、非萌生生殖次数、非萌生枝数量、萌生枝长、萌生枝径、萌生生殖次数、萌生枝数量和萌生率九个方面来考虑施肥对样地有无明显的改善。首先用Excel对个体构件调查表的数据进行处理，画出每个方面在施肥和不施肥两种情况下的折线图，浅层次的判断出施肥的效果。再运用方差分析，对九个方面分别进行显著性检验。最后得出施肥对样地的具体影响结果，见下表：非萌生枝枝长非萌生枝径萌生枝生殖次数萌生率效果抑制抑制促进抑制显著性显著显著显著显著然后对非萌生枝长、非萌生枝径、萌生枝长和萌生枝径进行区间估计，求出在置信度为95%的置信区间，求得施肥对枝长或枝径增长值的定量数据结果，见下表：非萌生枝枝长非萌生枝径萌生枝长萌生枝径置信区间【-23.848，-3.554】【-0.5256，-0.1664】【-13.2552,27.8052】【-0.4123,0.1463】问题2，我们根据对硼肥作用和麻疯树果实发育阶段划分与性状变化动态的研究，选择坐果率和出油量作为衡量果实产量的指标，另外由于不同方位麻疯树光合作用不同，我们分别从东、南、西、北、中五个方位研究坐果率。首先对数据分析画出在施肥和不施肥两种情况下，五个方位坐果率和每株麻疯树平均出油量的折线图，然后分别求出施肥样地和不施肥样地的坐果率和出油量，比较大小得出施肥能起到减少坐果率、提高出油量的效果，影响效率分别为36.65%和35.49%。关键词：浅层次判断方差分析显著性检验区间估计11.问题重述生物柴油作为生物质液体燃料，受到了世界主要石油进口国的高度重视，已成为世界各国国家能源战略资源的重要组成部分。开发生物柴油是我国石油安全战略的重要组成部分，是实现能源、经济和生态效益的良性循环共同发展的战略途径之一。目前，四川已成为了全国规模最大的麻疯树资源分布地区，但是麻疯树丰产栽培与资源高效培育尚未在生产中大面积应用，急需开展丰产栽培技术试验与应用，尽快研究出高产和稳产的定向培育技术，满足大面积高产原料培育的生产需要，推进麻疯树原料资源规范化、规模化培育，为生物柴油产业化发展打好资源基础。其中对麻疯树配方施肥试验的研究，是麻疯树果实（种子）产量高产、稳产问题中一个重要的技术问题。现以四川攀西地区的某地的硼肥喷施实验结果（见附表）为基础：问题1：利用个体构件调查表的数据，通过对施肥样地和未施肥样地的数据的对比分析，判断对样地的施肥的效果有无明显的改善。若有改善，给出施肥对枝长或者枝径增长值的定量数据结果。问题2：利用开花座果调查表的数据，分析施肥对麻疯树果实产量有无明显的影响？若有，试建立相应的数学模型并给出定量的数据结果。2.问题假设1.植株同一方向上所受光合作用相同；2.施肥对花的总数没有影响；3.只考虑枝径与枝长之间的相互作用；4.植株生长过程中相互独立，互不影响；5.不考虑施肥以外的因素对植株的影响；6.对所有植株施肥剂量相同。7.第2、3茬全为皱果或圆果。3.符号说明d出油量；1ijn皱果数；2ijn圆果数；33ijn黄果数；4ijn干果数；1,2i=表示施肥或未施肥；1,2ijn=…在同因素下不同的观测值；21,2,3,4,5k=五个方位；m所有数据的总均值；im施肥组或未施肥组的均值；ia施肥组或未施肥组均值与总体均值之差ijx随机误差项，且2(0,)ijNxs；ikja坐果率；ikjb雌花数目；ikjc总结果数；ijd出油量；h施肥的影响效率。4.问题分析问题1，首先我们查找得到了硼肥的作用：（1）促进作用，硼能促进碳水化合物的运转，植物体内含硼量适宜，能改善作物各器官的有机物供应，使作物生长正常，提高结实率和坐果率；（2）特殊作用，硼对受精过程有特殊作用，它在花粉中的量，以柱头和子房含量最多，能刺激花粉的萌发和花粉管的伸长，使授粉能顺利进行。通过对个体构件调查表的数据进行分析，我们得到每株麻疯树的枝数各不相同，每条枝的枝长、枝径和生殖次数也都不相同，而且两种条件下的样地的萌生枝比率也不同，所以我们可以从枝数、枝长、枝径、生殖次数和萌生率五个方面进行对比分析。第一步，运用Excel求出施肥样地和未施肥样地每株麻疯树的平均枝长、枝径和生殖次数。在施肥样地和未施肥样地两种情况下，分别画出平均枝长、枝径、生殖次数和枝数的折线图，从折线图可以看出施肥对平均枝长、枝径、生殖次数和枝数具有改善作用。第二步，运用Matlab软件分别对平均枝长、枝径、生殖次数和枝数进行方差分析，这样可以得出施肥分别对于四者的显著性情况。第三步，对于萌生枝比率，可以通过比较大小的方法求得。这样就能判断出施肥对于样地有无显著性改变。对于施肥对枝长或枝径增长量的定量描述，可以采用区间估计的方法。以施肥与未施肥的枝长（或枝径）长度均值之差来表示增长值，求出95%置信度下的估计区间就是施肥对枝长或枝径增长量的定量数据结果。问题2，通过对开花座果调查表的数据分析，得到植株在不同方位的开数目不同，使得果数不同，基于此，我们分别对不同方位的开花座果情况进行研究。要研究施肥对对麻疯树果实产量的影响，必须先找到衡量果实产量的指标。第一步，题目中已经假设施肥不对开花数量起作用，而且了解到硼肥主要增加植株的坐果率，所以我们可以以坐3果率（=雌花数目坐果率结果数）来衡量施肥对对麻疯树果实产量的影响。第二步，通过对麻疯树果实发育阶段划分与性状变化动态的研究，得知麻疯树在皱果期、圆果期、黄果期和干果期的水分含量和颜色有明显的差别，这导致了麻疯树的出油率在四个阶段不同。经过研究我们得到皱果、圆果、黄果和干果的出油参数分别为0.9、0.9、0.5和0.2。这样我们可以吧出油量（12340.90.90.50.2dnnnn=????）作为衡量施肥对对麻疯树果实产量的影响的指标。第三步，分别求出施肥样地和没施肥样地每株麻疯树的坐果率和出油量，画出折线图，就能浅层次的判断出施肥对麻疯树果实产量的影响。第四步，运用Matlab软件分别对坐果率和出油量进行方差分析，这样可以得出施肥分别对于二者的显著性情况。第五步，对于施肥对坐果率和出油量的定量描述，同样可以采用区间估计的方法。以施肥与未施肥的坐果率（或出油量）均值之差来表示增长值，求出95%置信度下的估计区间就是施肥对坐果率（或出油量）增长量的定量数据结果。5.模型建立模型1：通过对个体构件调查表的数据的分析，我们从枝数、枝长、枝径、生殖次数和萌生率五个方面来判断施肥对样地有无明显的改善。首先进行浅层次的初步判断，然后进行方差分析判断改善的显著性。萌发枝为萌发新枝的树枝，其各项数据都与非萌发枝有一定的区别，在这里，我们对萌发枝和非萌发枝分开进行分析。下面以非萌发枝枝长在施肥和不施肥两种情况下的数据分析为例。第一步：首先运用Excel求出施肥样地和未施肥样地每株麻疯树的平均枝长，这样得到了两组数据，再画出折线图。通过折线图可以初步判断出施肥对枝长的改善。第二步：对非萌发枝枝长进行方差分析，对施肥和未施肥非萌发枝麻疯树的平均枝长进行统计，检验得出数据服从正态分布。假设检验：012:Hmm=112:Hmm¹方差分析公式：ijiijcmaz=++则全体数据对x的偏差平方和2211()inTijijSxx===-邋各组均值对总方差的偏差平方和，即组间平方和221()AiiiSnxx·==-å4各组内的数据对均值方差偏差平方和得总和,即组内平方和2211()inEijiijSxx·===-邋TAESSS=+注意到21()inijijxx·=-å是总体2(,)iNms的样本方差的1in-倍，于是有2221()/(1)inijiijxxnsc·=--å由2c分布的可加性知2221/(1)EiISnsc=轾犏-犏臌å即22/(2)ESnsc-2(2)ESEns=-对AS作进一步分析得2221(21)ASiiIEnsa==-+å当0H成立2ASEs=可知若0H成立，AS只反映随机波动，若0H不成立，它就反映了A的不同水平的效应ia。但从数值上看/(21)1/(2)AESSn-»-该值服从自由度11n=，22nn=-的F分布。/(21)(1,2)/(2)AESFFnSn-=?-为检验0H，给定显著性水平a，记F分布的1a-分位数为1(1,2)Fna--。检验规则为1(1,2)FFna--时，接受0H；否则拒绝。用Matlab即可求得枝长的显著性水平P，若0.05P,则为极显著；若0.050.1P,则为显著；若0.1P,则为不显著。这样就能精确的判断出施肥对非萌生枝长有无明显的改善。5第三步：用第一步和第二步的方法判断出施肥分别对非萌生枝径、非萌生生殖次数、非萌生枝数量、萌生枝长、萌生枝径、萌生枝生殖次数和萌生枝数量有无明显的改善。第四步：考虑枝长与枝径之间的相互关系。首先对施肥后的平均枝长、平均枝茎进行回归分析，选择拟和度最大的线，进行参数估计最终得到平均枝茎关于平均枝长的关系式11()yfx=。同理，我们对未施肥树的平均枝长、平均枝茎进行回归分析，平均枝茎关于平均枝长的关系式22()yfx=。则321()()yfxfx=-表示没有施肥的与施过肥的平均枝茎之差，用Matlab作出图形观察判断施肥对枝径的改善效果。对11()yfx=和22()yfx=函数进行求导得：11'()xfy=和22'()xfy=，则321xxx=-表示没有施肥的与施过肥的平均枝长之差，用Matlab作出图形观察判断施肥对枝长的改善效果。第四步：判断施肥对萌生率有无明显的改善。由=萌生枝数萌生率总枝数分别算出施肥样地和未施肥样地的萌生率，比较二者大小即可得出施肥对萌生率的影响。综合上面四步可判断出施肥对样地有无明显的改变。下面给出若有改变的情况下，施肥对枝长或枝径增长值的定量数据结果。我们采用区间估计的方法来计算施肥对枝长或枝径增长值的定量数据结果，以非萌生枝枝径为例。分析可知，施肥与未施肥的非萌生枝径均服从正态分布，记为：211,N，222,N。设置信度为1，1210,,,XXX为施肥的枝径总体的样本，1210,,,YYY为未施肥的枝径总体的样本，,XY和2212,SS分别为两个样本的均值和方差。,XY分别为12,的无偏估计，可以得出XY是12的无偏估计，由概率论知识可知211,10XN222,10YN221212,1010XYN因此12的一个置信度为1的近似置信区间为：6221221010XYz可以得到121111()101niiSXX222211()101niiSXX因此12的一个置信度为1的近似置信区间可写为：221221010SSXYz其中，0.05a=，21.96Z对每株非萌生枝麻疯树求平均枝径，得到两个样本的值，用Matlab即可得到施肥与未施肥的非萌生枝枝径的平均长度之差的95%的置信区间，即得到了施肥对非萌生枝枝径增长值的定量数据结果。用同样方法可得到施肥对非萌生枝枝长、萌生枝枝径和萌生枝枝长增长值的定量数据结果。模型2：要求施肥对麻疯树果实产量有无明显改善，必须先找出衡量疯树果实产量的指标。题目中已经假设施肥不对开花数量起作用，而且了解到硼肥主要增加植株的坐果率，所以我们可以以坐果率（=雌花数目坐果率结果数）作为衡量施肥对麻疯树果实产量影响的指标。通过对麻疯树果实发育阶段划分与性状变化动态的研究，得知麻疯树在皱果期、圆果期、黄果期和干果期的水分含量和颜色有明显的差别，这导致了麻疯树的出油率在四个阶段不同。经过研究我们得到皱果、圆果、黄果和干果的出油参数分别为0.9、0.9、0.5和0.2。这样我们可以吧出油量（12340.90.90.50.2dnnnn=????）作为衡量施肥对对麻疯树果实产量的影响的指标。这样问题2就可以转化为“分析施肥对麻疯树坐果率和出油量有无明显的影响，并给出定量的数据结果”。从开花座果调查表的数据表得到，不同方位的麻疯树雌花