蘭興利，鄒志榮, 宛 金

(西北農(nóng)林科技大學(xué) 園藝學(xué)院，陜西 楊凌 712100)

0 引言

在溫室下培育作物的方式在我國(guó)已經(jīng)得到了大面積的推廣，這樣的方式最大的優(yōu)勢(shì)就是讓人們?cè)谌魏渭竟?jié)都能吃到新鮮的反季節(jié)蔬菜甚至是水果。然而，溫室條件下作物生長(zhǎng)環(huán)境的調(diào)控需要嚴(yán)格地技術(shù)支持，這也導(dǎo)致了溫室蔬果的產(chǎn)能遠(yuǎn)遠(yuǎn)滿足不了人們的實(shí)際需求。

正交試驗(yàn)法非常適合多因素組合試驗(yàn)尋求最優(yōu)方案的情形。正交試驗(yàn)方法問(wèn)世以來(lái)備受科研工作者青睞是因?yàn)槠湓跍p少試驗(yàn)次數(shù)的同時(shí)又能保證試驗(yàn)結(jié)果的均勻分布和良好的可靠性。然而，正交試驗(yàn)法也有其弊端,比如得出的最優(yōu)方案只能是各變量因子在事先所規(guī)定的水平內(nèi)的線性組合，其量化精度不足。這樣的精度在溫室環(huán)境管控中是不能接受的，因?yàn)闇厥覂?nèi)環(huán)境控制設(shè)備的調(diào)節(jié)是連續(xù)進(jìn)行的，因此各個(gè)因子的組合方案理論上是有無(wú)數(shù)種的，即使是進(jìn)行全面試驗(yàn)也不能對(duì)所有方案進(jìn)行試驗(yàn)驗(yàn)證。筆者采取的方法是在二次正交試驗(yàn)的基礎(chǔ)上，根據(jù)最終的二次回歸方程對(duì)參量因子重新量化編碼后采用改進(jìn)后的遺傳算法進(jìn)行參數(shù)尋優(yōu)，得到精確度更高的多環(huán)境因子配置方案

1 材料和試驗(yàn)方法介紹

1.1 材料和試驗(yàn)地點(diǎn)

試驗(yàn)地點(diǎn)為溫室，選培的是 “粉提一號(hào)”番茄。該品種是西安番茄研究所培育的無(wú)限生長(zhǎng)型櫻桃番茄。

1.2 試驗(yàn)設(shè)計(jì)

采用四因素三水平正交試驗(yàn)，分別對(duì)番茄的苗期、開花座果期和結(jié)果期的環(huán)境調(diào)控配置方案進(jìn)行了試驗(yàn)研究，環(huán)境調(diào)控因子選擇的是溫度、光照、濕度和CO2濃度這四個(gè)因素。

試驗(yàn)采用分批育苗，前期在人工氣候箱中育苗。挑選大小均勻、飽滿的番茄種子經(jīng)浸種催芽后播于 72 孔的塑料穴盤中進(jìn)行育苗。對(duì)苗期、開花坐果期和結(jié)果期整個(gè)流程的環(huán)境調(diào)控因子數(shù)據(jù)和番茄的生長(zhǎng)狀態(tài)數(shù)據(jù)進(jìn)行記錄和分析。

1.3 結(jié)果期試驗(yàn)結(jié)果展示

首先，設(shè)計(jì)試驗(yàn)因素水平和編碼，如表1所示：

表1 試驗(yàn)因素水平及編碼

接下來(lái)，給出環(huán)境因子組合試驗(yàn)方案，如表2所示：

表2 環(huán)境因子組合試驗(yàn)方案

1.4 測(cè)定指標(biāo)

由于番茄培育中，生長(zhǎng)苗期、開花座果期和結(jié)果期的環(huán)境調(diào)控配置方案具有差異性。比如，生長(zhǎng)苗期、開花座果期對(duì)于光合速率、蒸騰速率等指標(biāo)的要求比較高，測(cè)定的時(shí)候需要測(cè)定番茄植株的地上部鮮重和地下部鮮重，然后根據(jù)壯苗指數(shù)公式計(jì)算出的結(jié)果來(lái)評(píng)估環(huán)境調(diào)節(jié)配置方案的優(yōu)良。而結(jié)果期，僅僅需要測(cè)定果實(shí)的產(chǎn)量，就能得到選定的環(huán)境調(diào)控配置因子與最終的果實(shí)產(chǎn)量之間關(guān)系，從而為結(jié)果期的番茄培育提供指導(dǎo)。

2 試驗(yàn)結(jié)果展示與建模

各個(gè)環(huán)境調(diào)控配置方案的產(chǎn)量結(jié)果和各水平下的優(yōu)先順序計(jì)算分別如表3和表4所示所示：

表3 方案產(chǎn)量結(jié)果展示

表4 主次順序計(jì)算

按照四因素三水平這樣的試驗(yàn)方案，最終得到結(jié)果期番茄的最優(yōu)培育方案是A2B3C1D2，結(jié)果期產(chǎn)量影響的因子主次順序?yàn)镈、C、A、B。四因素三水平試驗(yàn)的回歸模型以溫度(X1)、濕度(X2)、光照(X3)和CO2濃度(X4)為自變量，以最終的番茄產(chǎn)量Y為因變量。茄產(chǎn)量Y與環(huán)境調(diào)控因子的擬合關(guān)系為：

(1)

3 模型優(yōu)化

四因素三水平正交試驗(yàn)得出的最優(yōu)的環(huán)境因子組合為A2B3C1D2，即通過(guò)溫室環(huán)境調(diào)控設(shè)備的管控，使得結(jié)果期的番茄所處環(huán)境的溫度控制在26.4℃、空氣濕度控制在95.6%、光照控制在200(PAR)、CO2濃度控制在1 400 mol·m-2·s-1就能確保該實(shí)驗(yàn)條件下的理論產(chǎn)量最高值，最終的擬合回歸方程也能反映各因子組合與產(chǎn)量的關(guān)系。然而，該多因子輸入輸出模型的最大弊端在于溫室內(nèi)環(huán)境因子的管控精度不高，最終只有數(shù)量極其有限的幾個(gè)因子組合方案的有效性和優(yōu)異性有理論的支撐。實(shí)際中，環(huán)境因子的調(diào)節(jié)可以視為連續(xù)的，即因子的方案組合理論上有無(wú)數(shù)個(gè)，因此正交試驗(yàn)得出的最優(yōu)因子組合方案的精度是不夠的。采用遺傳算法來(lái)對(duì)正交試驗(yàn)得出的多因子輸入輸出模型求解，得到精度更好的高產(chǎn)量環(huán)境因子組合方案。

3.1 編碼

編碼采用二進(jìn)制編碼，見表5。

表5 編碼表

表中n1為每個(gè)環(huán)境控制變量因子的編碼位數(shù)?？筛鶕?jù)設(shè)備傳感器的精度和實(shí)際環(huán)境管控要求的精度來(lái)確定因子編碼的位數(shù)。比如文章中溫度因子的水平設(shè)計(jì)和編碼的跳躍間距為8.40，在筆者的試驗(yàn)中溫度因子的變化范圍是18～34.8，為了得到更為精切的溫度因子配置方案，這里對(duì)溫度因子進(jìn)行二進(jìn)制編碼。假設(shè)溫室內(nèi)，環(huán)境管控設(shè)備溫度傳感器的最小探測(cè)精度為0.1℃，此時(shí)設(shè)置編碼位數(shù)為8，則編碼的精度為16.8/255=0.066，這是滿足精度要求的編碼。

3.2 適應(yīng)度函數(shù)的選擇

四因素三水平正交試驗(yàn)得出的擬合回歸方程可以表述各個(gè)環(huán)境因子間的內(nèi)在關(guān)系，因此筆者選擇正交試驗(yàn)得出的擬合方程作為遺傳算法尋優(yōu)的適應(yīng)度方程。

3.3 選擇、交叉、變異機(jī)制的設(shè)計(jì)

(1)選擇操作。種群進(jìn)行每一次迭代時(shí)，選擇操作采用經(jīng)典“輪盤賭”機(jī)制，適應(yīng)度值高的染色體具有更高的被選擇概率。

(2)交叉操作。進(jìn)行交叉操作時(shí)是以種群中的每個(gè)個(gè)體為單位，采用單點(diǎn)交叉的方式。先產(chǎn)生一個(gè)隨機(jī)數(shù)，假如則隨機(jī)產(chǎn)生一個(gè)交叉點(diǎn)，兩兩配對(duì)的個(gè)體交叉點(diǎn)以后的所有元素進(jìn)行互換。

(3)變異操作。進(jìn)行變異操作時(shí)是以解空間里的每個(gè)元素為單位，先產(chǎn)生一個(gè)隨機(jī)數(shù)，假如則進(jìn)行變異操作，將原先元素進(jìn)行“取反”操作即可。

(4)模擬退火操作。每次變異操作之后，需要執(zhí)行模擬退火操作，以此來(lái)提高算法的全局搜索能力。在到達(dá)冷卻溫度之前，以概率接受新的最優(yōu)解，其中表示新個(gè)體適應(yīng)度值與當(dāng)前迭代次數(shù)內(nèi)的最優(yōu)解的適應(yīng)度值之差。

4 仿真試驗(yàn)

(1)這里設(shè)定溫度、濕度和光照輻射的編碼位數(shù)為8，CO2濃度因子的編碼位數(shù)為10，這樣保證每個(gè)因子的量化精度都比較合理。下面給出每一步適應(yīng)度值和迭代次數(shù)的關(guān)系圖，如圖1所示：

結(jié)果分析:(1)未引入模擬退火時(shí)算法收斂性很差，引入模擬退火算法后整體算法的收斂性和搜索效率有了很大幅度的提升。

圖1 搜索效率對(duì)比

(2)本次運(yùn)行的可行最優(yōu)解為溫度因子，濕度因子，光照輻射，CO2濃度因子。此時(shí)對(duì)應(yīng)的理論產(chǎn)量是169899.55。

(3)尋優(yōu)的結(jié)果相對(duì)比較接近正交試驗(yàn)得出的最優(yōu)的環(huán)境因子組合為，一定程度上也驗(yàn)證了正交試驗(yàn)的正確性。

5 結(jié)論

筆者針對(duì)溫室條件下環(huán)境因子管控的問(wèn)題進(jìn)行了一些創(chuàng)新性的研究。首先利用正交試驗(yàn)得出初步的環(huán)境因子配置方案和擬合方程；隨后利用遺傳算法和相應(yīng)的改進(jìn)算法對(duì)擬合方程進(jìn)行參數(shù)尋優(yōu)，針對(duì)性設(shè)計(jì)了編碼機(jī)制和選擇機(jī)制；最后，從尋優(yōu)的結(jié)果可以看出采用智能算法得出的解和正交試驗(yàn)得出的解是不矛盾的，這也驗(yàn)證了正交試驗(yàn)的正確性和科學(xué)性。另外，智能算法的尋優(yōu)結(jié)果反映出結(jié)果期的番茄產(chǎn)量保障對(duì)CO2濃度的需求程度比正交試驗(yàn)得出的需求量略高。如何將這一套環(huán)境因子數(shù)據(jù)分析處理機(jī)制和算法融入到實(shí)際的溫室設(shè)備管控系統(tǒng)中是下一階段的研究?jī)?nèi)容。