基于集合经验模态分解与Elman神经网络的线椒株高预测

第 31 卷 第 18 期 农 业 工 程 学 报 Vol.31 No.18 2015 年 9月 Transactions of the Chinese Society of Agricultural Engineering Sep. 2015 169 基于集合经验模态分解与 Elman 神经网络的线椒株高预测张 瑜，汪小旵※，孙国祥，李永博 （ 1. 南京农业大学工学院，南京 210031； 2. 江苏省现代设施农业技术与装备工程实验室，南京 210031） 摘 要 为提高温室环境控制系统的有效性，针对作物生长量的变化与环境因子的变化存在时间尺度不统一的问题，该文基于集合经验模态分解（ ensemble empirical mode decomposition， EEMD）与 Elman 神经网络建模，提出一种线椒株高生长量预测方法。以 8819 线椒为试验对象，分别对线椒株高及其环境因子进行 EEMD 分解，对各尺度下的时间序列建立 EEMD-Elman 预测模型。结果表明应用 EEMD-Elman 神经网络建立线椒株高生长量预测模型，模型预测值与实测值的平均绝对误差为 1.69 cm，相关决定系数为 0.996，标准误为 1.104，模型预测结果与实测值呈极显著性相关。研究结果可以解决作物生长变化与环境变化时间尺度不统一的问题，为温室环境控制系统的控制目标的优化提供有效参数。 关键词 神经网络；模型；优化；集合经验模态分解； Elman 神经网络；线椒株高；预测 doi 10.11975/j.issn.1002-6819.2015.18.024 中图分类号 S625.5 文献标志码 A 文章编号 1002-68192015-18-0169-06 张 瑜，汪小旵，孙国祥，李永博. 基于集合经验模态分解与 Elman 神经网络的线椒株高预测[J]. 农业工程学报，2015，3118169－174. doi 10.11975/j.issn.1002-6819.2015.18.024 http//www.tcsae.org Zhang Yu, Wang Xiaochan, Sun Guoxiang, Li Yongbo. Prediction of cayenne pepper plant height based on ensemble empirical mode decomposition and Elman neural network[J]. Transactions of the Chinese Society of Agricultural Engineering Transactions of the CSAE, 2015, 3118 169－ 174. in Chinese with English abstract doi 10.11975/j.issn.1002-6819.2015.18.024 http//www.tcsae.org 0 引 言作物生长量预测为温室环境调控目标优化提供重要依据[1-2]。目前，温室作物生长量预测多通过建立作物外观形态模拟模型的方法实现[3-6]。 如， calvalho 等[7]， Larsen等[8]， Heuvelink 等[9]， Nothnagl 等[10]分别建立了菊花生长模拟模型，对菊花的株高、花朵数、花朵大小等生长量进行预测。上述模型研究的最直接意义是及时掌握温室环境控制系统对作物所产生的作用，为温室环境优化控制提供理论依据。然而，作物生长模拟模型无法真实地反映作物的实际生长情况，使得环境控制系统无法根据作物实际生长情况做出有效的响应，降低了温室环境控制系统的有效性。所以，根据作物实际生长情况对作物生长量进行预测，可以更真实地了解环境控制的效果[11]，提高温室环境控制系统的有效性[12]。 根据作物实际生长情况对环境进行调控是提高温室环境控制系统有效性的重要方法。但是作物生长量的变化与环境因子的变化存在时间尺度不统一的特点，作物生长量的变化是缓慢的，而环境因子的变化是瞬时的。收稿日期 2015-06-04 修订日期 2015-08-24 基金项目国家自然科学基金项目（ 61273227） 、江苏省自然科学基金（ BK20150686） 作者简介张瑜，女（汉） ，吉林长春人，讲师，博士，农业电气化与自动化、检测技术与自动化装置。南京 南京农业大学工学院， 210031。 Email ※通信作者汪小旵，男（汉） ，教授，博士，从事农业生物环境工程、设施农业温室环境模拟与调控方面研究。 南京 南京农业大学工学院， 210031。Email 单纯采用平均值对两个时间尺度不统一的时间序列进行建模、预测，无法清晰地反映两种时间序列在不同时间尺度下的波动情况，造成预测结果不准确。因此，需要对两者进行多时间尺度分解，利用不同时间尺度下的时间序列特征进行预测与重构，实现作物实际生长量的预测，为温室环境控制系统提供实际有效的控制依据。 线椒是温室中常见种植作物，株高是反映线椒长势的一项重要生长参数，能够直观地反映出在温室环境控制系统对作物生长的作用效果。因此，本文基于集合经验模态分解（ ensemble empirical mode decomposition，EEMD） 与 Elman 神经网络， 对线椒株高及其环境参数时间序列进行分解与重构，分别对分解后的时间序列建立非线性 Elman 神经网络，根据种植期间，线椒株高当前的实际生长量，重构并预测线椒株高在同一生长期内未来的生长量，为温室作物生长预测提供精确、实时的预测方法，为温室环境控制系统控制目标的优化提供实时有效的控制依据，实现提高温室环境控制系统有效性的目的。 1 数据与方法 1.1 试验数据采集 本文试验对象为 2014 年种植的 8819 线椒 定植后生长期内， 2014 年 5 月 28 日至 2014 年 7 月 8 日 及对应的温室内环境因子 （温室内温度、 相对湿度、 室内总辐射） 。试验于江苏省南京市南京农业大学工学院，博远楼顶venlo 型玻璃温室内进行。温室环境因子由 Eco-Watch 生态网络监测系统每 15 min 采集一次，取日平均值用于模农业工程学报（ http//www.tcsae.org） 2015 年 170 型建立。由于温室内温度湿度等环境因子分布的不均匀性，所以本文选取在温室内东面密集种植的 16 株线椒作为试验样本，以降低温室内不同小气候条件对作物生长形态的影响。作物均生长在适宜的环境条件下，未发生环境胁迫与营养胁迫情况。其中 10 株为测试集，用于建立模型，另外 6 株为验证集，用于模型验证。考虑到作物的个体差异，均取其平均值。网络的输入量为温室内环境因子的日平均值，输出量为测试集中 10 株线椒株高的每日测量值的平均值，测量方法采用文献 [13]中线性激光扫描法。线椒采用椰糠基质，独立种植于栽培袋中。灌溉营养液配方采用山东农业大学番茄辣椒配方， Ec 值为 2.0， pH 值为 6.5。 1.2 分解方法 在解决时间序列多尺度分解的方法中，小波分解和奇异谱分解等方法均需要依赖于预先给定的基函数，无法根据数据自身的时间尺度特征进行分解，可能会由此分解出虚假或无效的波动分量[14-20]。集合经验模态分解方法利用高斯白噪声具有频率均匀分布的统计特性，无需预先设定任何基函数，自适应地将数据分解为有限个本征模函数（ intrinsic mode function， IMF），所分解出来的各 IMF 分量包含了原始数据的不同时间尺度的局部特征信号[21-26]。所以，本文采用 EEMD 分解方法对试验数据（温室环境因子与线椒株高生长量）进行多尺度分解，分解结果可以由公式（ 1）求得，重构结果由公式（ 2）求得。 11 NjijiCt CtN∑（ 1） 1 nj njX tCtr∑（ 2） 式中 Cijt表示第 i次加入白噪声后分解所得的第 j个 IMF分量， t 表示时间； Cjt为分解得到的第 j 个 IMF 分量；N 为添加白噪声序列的数目； Xt为原始信号； rn为分解后的残留余项， n 表示分解次数，当 rn为单调信号时，分解终止。 1.3 预测方法 Elman 神经网络是一种带有反馈的两层 BP 网络结构，其反馈连接是从隐含层的输出到其输入端。这种反馈方式使得 Elman 网络能够探测和识别时变模式。其隐含层又称为反馈层，这种特殊的两层网络可以任意精度逼近任意函数，唯一的要求是其隐含层必须具有足够的神经元数目。隐含层神经元数越多，则逼近复杂函数的精度就越高[27-28]。 Elman 神经网络的非线性表达原理[29]如式（ 3）所示。 213 1 1 −⎧⎪⎨⎪⎩ccxk xkxk f w uk wx kyk gwxk（ 3） 式中 y、 x、 u、 xc分别表示 m 维输出节点向量、 n 维中间层节点单元向量、 l 维输入向量和 n 维反馈状态向量； w3、w2、 w1分别表示中间层到输出层、输入层到中间层、承接层到中间层的连接权值。 g 为输出神经元的传递函数，是中间层输出的线性组合。 f 为中间层神经元的 传递函数。 2 结果与分析 2.1 EEMD 分解与重构 图 1 为原始时间序列。根据式（ 1）本文分别对线椒株高及其环境因子（温室温度、相对湿度、室内总辐射）进行 EEMD 分解，分别得到了 5 个本征模函数，分别为IMF1～ IMF5，如图 2 所示。本征模函数表征了数据的内在的振动模式。图 2 中， IMF1～ IMF4 反应了线椒株高和各个环境因子在不同频率下振荡情况， IMF5 为变化趋势项，反应了线椒株高在生长期内的整体变化趋势。 EEMD分解重构后结果的相关参数如表 1 所示。由表 1 可知，各个环境因子和线椒株高时间序列经过 EEMD 分解重构后，与原始时间序列相关系数在有效数字范围内近似为1，标准误接近于 0，平均绝对误差接近于 0，表明重构结果十分有效。 图 1 原始时间序列 Fig.1 Original time series 2.2 Elman 神经网络生长量预测与分析 本文 Elman 网络的输入层为温室内环境因子，包括温室温度、相对湿度和室内总辐射，输入节点数为 3；输出层为线椒株高生长量的预测值，节点数为 1。本文采用双层反馈结构，第一层节点数为 10，第二层节点数为 3。反馈层神经元的传递函数为 tansig，线性输出层的传递函数为 purelin。网络训练样本组按照 0.7∶ 0.15∶ 0.15的比例划分为训练集、验证集和测试集。 为避免神经元饱和，在输入层对输入数据进行归一化处理，将各个数值换算至 [0 1]区间内，在输出层对得到的预测结果进行反归一化。根据 Elman 神经网络的非线性空间表达式（公式（ 3）），本文网络模型的输出结果可以表示为 1121122 2132 1 1 1 −−⎧⎪⎨⎪⎩xk fwxk w ukxk fwxk w xkyk gwx k（ 4） 式中 uk−1为输入层的输入； x1k为第 1 层反馈的输出；x2k为第 2 层网络的输出， yk为输出结果； wp1（ p1,2；p 为网络层数） 为第 1 层网络和第 2 层网络的承接层到中间层的连接权值； wp2为第 1 层网络和第 2 层网络的输入层到中间层的连接权值； w3为中间层到输出层的连接权值，本文中间层神经元传递函数采用双曲正切传递函数第 18 期 张 瑜等基于集合经验模态分解与 Elman 神经网络的线椒株高预测 171 tansig，输出层神经元采用线性传递函数 purelin。由双曲正切函数定义，可得到本文预测网络模型中间层传递函数表达式为 21212[ 1 1]2[ 1 1]11− −− −−ppppwuk wxkwuk wxkefe（ 5） 式中 p1,2；表示网络层数。 a. IFM1 b. IMF2 c. IMF3 d. IMF4 e. IMF5 注 IMF1～ IMF4 反应了线椒株高和各个环境因子在不同频率下的振荡情况， IMF5 反应了线椒株高和各个环境因子的变化趋势。 Note IMF1-IMF4 reflect the oscillation of variate, including plant height and every environment elements under different frequencies. IMF5 reflects the trend of plant height and every environment elements. 图 2 线椒株高与环境因子 EEMD 分解后各模态分量波动情况 Fig.2 Oscillation of EEMD results of plant height and environment elements 表 1 EEMD 分解重构结果分析 Table 1 Summary of EEMD and reconstruction 时间序列 Time series 决定系数 R2标准误 Standard error 平均绝对误差 Mean absolute error线椒株高 Plant height 1.0000 2.70810-34.5810-5cm 温室温度 Temperature 1.0000 4.7310-23.9710-4℃ 相对湿度 Relative humidity 1.0000 2.51610-35.1010-5 室内总辐射 Total inner radiation 1.0000 3.44510-37.2010-5 W·m-2本文中时间序列样本经过 EEMD 分解后，得到 5 个独立的本征模函数 IMF1～ IMF5，所以需要建立 5 个子Elman 神经网络，分别预测生长量的各个时间序列值，最后经过 EEMD 重构得到最终的线椒株高生长量预测值。具体预测流程如图 3 所示。 各子网络的预测结果统计性描述如表 2 所示，各个子网络的预测值与目标值之间相关决定系数分别为 0.71、0.88、 0.94、 0.95 和 0.98，平均绝对误差值分别为 0.04、0.09、 0.22、 0.14 和 1.20 cm。由表 2 可知，各尺度预测值与实测值的模型拟合效果均表现出极强的相关性，平均绝对误差均可满足测量精度要求。随着时间尺度的增加，模型相关性决定系数逐渐增强，平均绝对误差逐渐增大， IMF5子网络预测结果的平均绝对误差增大最明显，可能的原因是 IMF5 为趋势项， 对于变量的波动情况表达不明显。 本文对 IMF1～ IMF5 子网络的预测结果进行重构。重构后，线椒株高预测值与实际测量值之间的拟合关系如图 4 所示。结果表明通过 EEMD 多尺度分解，对线椒株高在不同时间尺度上应用 Elman 神经网络预测再合并，预测值与测量值达到极显著相关水平，平均绝对误差值为1.69 cm，可以满足温室中作物生长测量的精度要求。 农业工程学报（ http//www.tcsae.org） 2015 年 172 图 3 线椒株高预测流程图 Fig.3 Flowchart of pepper plant height prediction 表 2 Elman 网络预测线椒株高结果统计分析 Table 2 Summary of Elman neuron network pepper plant height prediction 网络类型 Type of network 决定系数 R2平均绝对误差 Mean absolute error/cm Elman 子网络 _IMF1 Sub-network_Elman_IMF1 0.71 0.04 Elman 子网络 _IMF2 Sub-network_Elman_IMF2 0.88 0.09 Elman 子网络 _IMF3 Sub-network_Elman_IMF3 0.94 0.22 Elman 子网络 _IMF4 Sub-network_Elman_IMF4 0.95 0.14 Elman 子网络 _IMF5 Sub-network_Elman_IMF5 0.98 1.20 注 SE 为标准误， R2为决定系数。 Note SE is standard error, R2 is coefficient of determination. 图 4 线椒株高分解重构预测结果 Fig.4 Results of pepper plant height prediction 分别应用 EEMD 与 BP 神经网络结合方法[30]、 EEMD与 Elman 神经网络结合方法和未经 EEMD 分解的 Elman神经网络方法进行线椒株高生长量的预测。取试验样本验证集数据，对所建模型进行验证。验证结果表明，经过 EEMD 分解的网络模型的预测值和实测值进行比较，其预测值与实测值均可较好的吻合，未经过 EEMD 分解的 Elman 网络预测效果稍差。采用决定系数 R2、标准误SE 和平均绝对误差两个指标分别对 3 个预测模型进行评价，各指标统计结果见表 3。 表 3 3 种预测方法线椒株高结果评价 Table 3 uation of three pepper plant height prediction s 网络类型 Type of network 决定系数 R2标准误 Standard error 平均绝对误差 Mean absolute error/cm EEMD-BP 神经网络 EEMD-BP neural network0.812 7.012 5.40 EEMD-Elman 神经网络 EEMD-Elman neural network0.996 1.104 1.69 Elman 神经网络 Elman neural network 0.908 5.032 8.87 由表 3 可知，经过 EEMD 分解后模型预测结果比未经过分解的模型预测结果准确性高，经过 EEMD 分解后的 Elman 神经网络的预测结果优于 EEMD-BP 神经网络和未经过 EEMD 分解的 Elman 神经网络预测结果。 同时，由图 4 可知， EEMD-Elman 神经网络的预测结果决定系数 R2为 0.996，表明 EEMD-Elman 神经网络输出值与目标值偏差较小，是十分有效的预测方法。 3 结 论 应用集合经验模态分解方法对线椒株高及其环境因子（温室温度、相对湿度、室内总辐射）进行多尺度分解与重构。重构后各个环境因子和线椒株高的时间序列与原始时间序列相关系数在有效数字范围内近似为 1， 标准误接近于 0，平均绝对误差极接近于 0，重构结果十分有效。 表明集合经验模态分解法 （ ensemble empirical mode decomposition, EEMD）分解可以解决多尺度时间序列预测问题。应用 EEMD-Elman 神经网络建立线椒株高生长量预测网络，网络预测值与实测值的平均绝对误差为1.69 cm，相关决定系数为 0.996，标准误为 1.104，表明预测结果与实测值呈极显著性相关。 3 种预测方法的预测结果进行对比，表明与未经过 EEMD 分解的预测网络相比，经过 EEMD 分解的预测网络，能够反映时间序列自身固有的特性，更好地诠释时间序列在不同频率的振荡情况；将经过 EEMD 分解的 EEMD-Elman 神经网络与EEMD-BP 神经网络的预测结果相比， EEMD-Elman 神经网络预测效果优于 EEMD-BP 神经网络预测效果，说明Elman 的双反馈结构可以提高模型的预测精度。 由于本文仅对生长期内的线椒株高进行分析，并未考虑线椒在幼苗期、开花期及坐果期的生长情况，所以本文中预测网络仅适用于定植后线椒株高的预测。试验期间，作物生长未经受任何环境胁迫与营养胁迫，如外界环境发生较大变化时，会对网络预测结果造成一定影响。所以，该预测网络在适应性和稳定性方面存在一定局限性。后续研究应更多关注作物在不同生长阶段，不同生长参数的预测研究以及环境胁迫与营养胁迫对预测结果精度带来的影响。 综上所述，集合经验模态分解法可以根据时间序列自身时间尺度特征，自适应地将时间序列分解到不同时间尺度上，清晰地表达出原始时间序列在不同时间尺度第 18 期 张 瑜等基于集合经验模态分解与 Elman 神经网络的线椒株高预测 173 上的波动情况，有效解决温室中环境因子与作物生长量时间尺度不统一的问题。应用 Elman 神经网络双层反馈结构可以提高时间序列的预测精度。 EEMD 与 Elman 神经网络结合可以准确预测作物实际生长量，为温室环境控制系统控制目标优化提供有效参数，实现有效提高温室环境控制系统的控制效果的目的。 [参 考 文 献] [1] Straten G, van Challa H, Buwalda F. Towards user accepted optimal control of greenhouse climate[J]. Computers and Electronics in Agriculture, 2000, 263 221－ 238. [2] Pucheta J A, Schugurensky C, Fullana R, et al. Optimal greenhouse control of tomato seedling crops[J]. Computers and Electronics in Agriculture, 2006, 501 70－ 82. [3] 倪纪恒，罗卫红，李永秀，等 . 温室番茄干物质分配与产量的模拟分析 [J]. 应用生态学报， 2006， 175 811－ 816. Ni Jiheng, Luo Weihong, Li Yongxiu, et al. Simulation of greenhouse tomato dry matter partitioning and yield prediction[J]. Chinese Journal of Applied Ecology, 2006, 175 811－ 816. in Chinese with English abstract [4] 李永秀，罗卫红，倪纪恒，等 . 温室黄瓜干物质分配与产量预测模拟模型初步研究 [J]. 农业工程学报， 2006， 222116－ 121. Li Yongxiu, Luo Weihong, Ni Jiheng, et al. Simulation of dry matter partitioning and yield prediction in greenhouse cucumber[J]. Transactions of the Chinese Society of Agricultural Engineering Transactions of the CSAE, 2006, 222 116－ 121. in Chinese with English abstract [5] Ma Yuntao, Wen Meiping, Guo Yan, et al. Parameter optimization and field validation of the functional-structural model GREENLAB for maize at different population densities[J]. Annals of Botany, 2008, 1018 1185－ 1194. [6] Yan Hongping, Kang Mengzheng, de Reffye P, et al. A dynamic, architectural plant model simulating resource- dependent growth[J]. Annals of Botany, 2004, 935 591－602. [7] Carvalho S M P, Heuvelink E, Cascals R, et al. Effect of day and night temperature on internode and stem length in chrysanthemum Is everything explained by DIF[J]. Annals of Botany, 2002, 901 111－ 118. [8] Larsen R U, Cecilia H. Predicting leaf unfolding in flower induced shoots of greenhouse grown chrysanthemum[J]. Scientia Horticulturae, 1995, 633/4 225－ 239. [9] Carvalho S M P, Heuvelink E, van Kooten O. Effect of light intensity, plant density and flower bud removal on flower size and number in cut chrysanthemum[J].Aeta Horticulturae, 2002, 5933 33－ 38. [10] NothnagI M, Kosiba A, Larsen R U. Predicting the effect of irradiance and temperature on the flower diameter of greenhouse grown chrysanthemum[J]. Seientia Hortieulturae, 2004, 993/4 319－ 329. [11] 王定成，乔晓军，汪春秀，等 . 基于作物响应的温室环境 SVMR 控制仿真 [J]. 农业工程学报， 2010， 26增刊 2290－ 293. Wang Dingcheng, Qiao Xiaojun, Wang Chunxiu, et al. Simulation of greenhouse environment SVMR control based on plant response[J]. Transactions of the Chinese Society of Agricultural Engineering Transactions of the CSAE, 2010, 26Supp.2 290－ 293. in Chinese with English abstract [12] 李永博，周伟，汪小旵，等 . 基于奇异摄动理论的综合温室控制系统设计 [J]. 农业机械学报， 2012， 435 184－189. Li Yongbo， Zhou Wei， Wang Xiaochan， et al． Greenhouse control system design based on singular perturbation theory[J]. Transactions of the Chinese Society for Agricultural Machinery, 2012, 435 184－ 189. in Chinese with English abstract [13] 张瑜，汪小旵，孙国祥，等 . 基于激光视觉的温室作物茎叶量测方法研究 [J]. 农业机械学报， 2014， 459 254－259. Zhang Yu, Wang Xiaochan, Sun Guoxiang, et al. Leaves and stems measurement of plants based on laser vision in greenhouses[J]. Transactions of the Chinese Society for Agricultural Machinery, 2014, 459 254－ 259. in Chinese with English abstract [14] Field R. Spatial species-richness gradients across scalesA meta-analysis[J]. Journal of biogeography, 2009, 361 132－147. [15] Crawford J. Multi-scale investigations of alpine vascular plant species in the San Juan Mountains of Colorado, USA GLORIA target region. Scientifica Acta, 2008, 22 65－ 69. [16] Foody G M.Spatial non-stationarity and scale-dependency in the relationship between species richness and environmental determinants for the sub-Saharan endemic avifauna[J]. Global Ecology and Biogeography, 2004, 134 315－ 320. [17] Keitt T H, Urban D L. Scale-specific inference using wavelets[J]. Ecology, 2005, 869 2497－ 2504. [18] Grinsted A, Moore J C, Jevrejeva S. Application of cross wavelet trans and wavelet coherence to geophysical time series[J]. Nonlinear Processes in Geophysics, 2004, 11 561－566. [19] Zhou Y, Tao T, Mei X, et al. Feed-axis gearbox condition monitoring using built-in position sensors and EEMD [J]. Robotics and Computer-Integrated Manufacturing, 2011, 274 785－ 793. [20] Ai-Ho L, Che-Chou S, Pai-Chi L. Potential contrast improvement in ultrasound pulse inversion imaging using EMD and EEMD[J]. Ultrasonics, Ferroelectrics and Frequency Control, IEEE Transactions on, 2010, 572 317－ 326. [21] Huang N E, Zheng S, Long S R, et al. The empirical mode decomposition and the Hilbert spectrum for nonlinear and non-stationary time series analysis[J]. Proc R Soc Lond, 1998, 454 903－ 995. [22] Huang Norden E, Attoh-Okine N O.The Hilbert-Huang trans in engineering[M]. England,Taylor 2. Jiangsu Modern Facilities Agricultural Technology and Equipment Engineering Laboratory, Nanjing 210031, China Abstract In order to improve the perance of greenhouse control system and figure out the multi-time scale variable problem between plant growth and environment elements in greenhouses, the ensemble empirical mode decomposition EEMD and Elman neural network were used to predict cayenne pepper height in this paper. Taking the cayenne pepper 8819 and its plant height and environment elements temperature, relative humidity, total inner radiation as research object, plant height and environment elements were decomposed by the of EEMD. The cayenne pepper plants sampled were irrigated with the nutrient solution with electricity conductivity of 2.0 and pH value of 6.5 from Shandong Agricultural University. Five intrinsic mode functions imfs were obtained by the EEMD, named imf1, imf2, imf3, imf4 and imf5. Oscillations of cayenne pepper plant height and environment elements at different frequency were shown by imf1, imf2, imf3 and imf4, whi