於沈剛,馬明舟,岳雪峰,萬(wàn) 衡,王運(yùn)圣
(1.上海應(yīng)用技術(shù)大學(xué),上海200235;2.大連理工大學(xué)城市學(xué)院,大連116600;3.上海市農(nóng)業(yè)科學(xué)院農(nóng)業(yè)科技信息研究所,上海201403;4.上海數(shù)字農(nóng)業(yè)工程技術(shù)研究中心,上海 201403)
20世紀(jì)90年代以來(lái)我國(guó)節(jié)水農(nóng)業(yè)發(fā)展進(jìn)入了新的時(shí)期,如今的節(jié)水灌溉技術(shù)已發(fā)展成節(jié)水灌溉工程技術(shù)和節(jié)水灌溉農(nóng)業(yè)技術(shù)為主并相互結(jié)合的技術(shù)體系[2]。由于人工智能技術(shù)的發(fā)展,一些控制理論如專家系統(tǒng)、模糊邏輯控制、神經(jīng)網(wǎng)絡(luò)控制、協(xié)同控制等被廣泛地應(yīng)用到工程問(wèn)題的預(yù)測(cè)和建模。這些理論基于自動(dòng)化傳感器設(shè)備,監(jiān)測(cè)土壤水分, 利用氣象數(shù)據(jù),估算植物蒸騰量,然后經(jīng)過(guò)控制理論的分析和處理,進(jìn)行邏輯判斷,進(jìn)一步將信號(hào)送到執(zhí)行機(jī)構(gòu),控制一些電磁閥的關(guān)閉,更好、更及時(shí)、更準(zhǔn)確地預(yù)測(cè)環(huán)境參數(shù),及時(shí)調(diào)整輸出響應(yīng)以適應(yīng)作物不同生長(zhǎng)時(shí)期對(duì)水分的需求。
溫室是隨著農(nóng)村經(jīng)濟(jì)結(jié)構(gòu)和作物種植結(jié)構(gòu)調(diào)整發(fā)展起來(lái)的新型設(shè)施農(nóng)業(yè),相對(duì)于大田農(nóng)業(yè)而言,其具有占地面積小,種植靈活,易于控制作物生長(zhǎng)環(huán)境等優(yōu)點(diǎn)。為了解決我國(guó)北方冬天蔬菜供應(yīng)短缺問(wèn)題,北方開(kāi)始大規(guī)模推廣和發(fā)展溫室農(nóng)業(yè)[3]。溫室是一個(gè)復(fù)雜的大慣性、純延時(shí)和非線性系統(tǒng),難以對(duì)其建立精確的數(shù)學(xué)模型,傳統(tǒng)的控制方法在實(shí)時(shí)性和控制精度上已不能滿足其控制要求[4]。模糊控制是一種非線性的控制方式,不需要建立對(duì)象的精確數(shù)學(xué)模型,根據(jù)領(lǐng)域?qū)<抑R(shí)或操作人員經(jīng)驗(yàn)就可制定針對(duì)對(duì)象有效的控制策略,響應(yīng)速度快,但控制誤差較大[5],PID控制器對(duì)于時(shí)變系統(tǒng)具有很好的魯棒性,適應(yīng)性強(qiáng)?;谝陨?種控制方法優(yōu)點(diǎn),針對(duì)溫室灌溉系統(tǒng)具有非線性、嚴(yán)重滯后的特性[6],本論文提出了一種PID控制和模糊控制相結(jié)合的模糊PID控制方法,并利用MATLAB搭建仿真模型證實(shí)其可行性。
通常溫室內(nèi)的水分環(huán)境包括空氣濕度和土壤濕度2個(gè)方面,而對(duì)植物的生長(zhǎng)發(fā)育起到關(guān)鍵作用的是土壤濕度。本文設(shè)計(jì)的控制器通過(guò)對(duì)土壤濕度變化的控制,使得土壤濕度一直穩(wěn)定在植物最佳土壤濕度值附近,從而實(shí)現(xiàn)對(duì)植物實(shí)時(shí)、適量的灌溉。
系統(tǒng)采用的是兩輸入三輸出的二維模糊控制策略,通過(guò)模糊控制調(diào)節(jié)PID控制器的3個(gè)修正參數(shù),提高系統(tǒng)的動(dòng)態(tài)響應(yīng)速度,增強(qiáng)系統(tǒng)對(duì)于外界條件變化的適應(yīng)性,同時(shí)PID控制器能夠有效地降低穩(wěn)態(tài)誤差,提高精度。模糊控制器的輸入量為傳感器測(cè)得的土壤濕度值y測(cè)與設(shè)定的最佳土壤濕度值y目的偏差E和偏差的變化率EC,輸出量為PID控制器的3個(gè)輸入修正參數(shù)ΔKp、ΔKi、ΔKd。系統(tǒng)的結(jié)構(gòu)框圖見(jiàn)圖1。
圖1 控制系統(tǒng)結(jié)構(gòu)框圖
圖1中誤差E值由式E=y目-y測(cè)得到,其中y目為系統(tǒng)設(shè)定的土壤濕度最佳值;y測(cè)為測(cè)得土壤濕度實(shí)時(shí)值。誤差的變化率由EC=E(k)-E(k-1)得到,E(k)為k時(shí)刻的土壤濕度值;E(k-1)為k-1時(shí)刻的土壤濕度值。將E與EC作為模糊控制器的2個(gè)輸入,通過(guò)設(shè)定變量因子進(jìn)行模糊化處理,制定模糊控制規(guī)則表,進(jìn)行模糊決策,輸出為PID控制器輸入的3個(gè)修正參數(shù)ΔKp、ΔKi、ΔKd,通過(guò)比例因子進(jìn)行解模糊化處理并作用于PID控制器輸出 值,最后作用于控制對(duì)象輸出灌溉水量,完成作物灌溉。
模糊PID控制器主要由模糊控制器和參數(shù)可調(diào)節(jié)PID控制器2部分組成,模糊控制器完成對(duì)參數(shù)可調(diào)節(jié)PID控制器3個(gè)參數(shù)的整定,參數(shù)可調(diào)節(jié)控制器的輸出完成對(duì)系統(tǒng)的控制[7]。PID控制器由比例環(huán)節(jié)(Proportional)、積分環(huán)節(jié)(Integral)和微分環(huán)節(jié)(Differential)組成,連續(xù)PID控制器的一般形式為[8]:
(1)
式中:Kp為比例系數(shù);Ki為積分系數(shù);Kd為微分系數(shù);u(t)為系統(tǒng)t時(shí)刻的控制作用;E(t)為誤差,EC(t)為誤差的變化率,這2個(gè)變量同時(shí)作為模糊控制器的2個(gè)輸入。
實(shí)驗(yàn)對(duì)象是甜瓜,考慮節(jié)水灌溉的實(shí)際條件和控制精度,定義偏差E的基本論域是[-5%,5%],偏差變化率EC的變化范圍是[-3%,3%],甜瓜的整個(gè)生長(zhǎng)期內(nèi)的最佳土壤濕度范圍為[13%,23%][9];根據(jù)滴灌實(shí)驗(yàn)可測(cè)的在深度20 cm處的土壤濕度從13%升高到23%用時(shí)約30 min,則設(shè)定輸出變量U的基本論域?yàn)閇0,30]。模糊子集的數(shù)量決定了對(duì)所研究對(duì)象的模糊論域的分割程度,模糊子集太少分割太粗,影響系統(tǒng)控制精確度;模糊子集太多,則模糊規(guī)則制定較多,系統(tǒng)控制較復(fù)雜,但控制精確度明顯提高?;谀:蛹贫ㄒ?guī)則以及控制系統(tǒng)的實(shí)際控制要求,定義E和EC的模糊語(yǔ)言值各有7個(gè):{負(fù)大,負(fù)中,負(fù)小,零,正小,正中,正大},簡(jiǎn)記為{NB,NM,NS,Z,PS,PM,PB},其對(duì)應(yīng)的模糊子集的論域?yàn)閧-3,-2,-1,0,1,2,3}。輸出變量u包括PID控制器的3個(gè)修正參數(shù)ΔKp、ΔKi、ΔKd,定義3個(gè)修正參數(shù)的模糊子集為{負(fù)大,負(fù)中,負(fù)小,零,正小,正中,正大},簡(jiǎn)記為{NB,NM,NS,Z,PS,PM,PB},其對(duì)應(yīng)的模糊子集論域?yàn)閧-10,-6.67,-3.33,0,3.33,6.67,10},考慮實(shí)際溫室灌溉系統(tǒng)還沒(méi)有措施可以快速有效降低土壤濕度,模糊PID控制器輸出U的上限為4,下限為0,避免了因人為或系統(tǒng)的故障導(dǎo)致土壤濕度遠(yuǎn)遠(yuǎn)高于設(shè)定閾值而帶來(lái)的農(nóng)業(yè)經(jīng)濟(jì)損失。為了進(jìn)行模糊化處理,必須利用誤差E的量化因子ke和誤差變化率EC的量化因子kec將輸入量由基本論域轉(zhuǎn)化到模糊論域。誤差E的量化因子ke和誤差變化率EC的量化因子kec以及輸出變化量U的比例因子ku由下式表示[10]:
電噴柴油機(jī)試車臺(tái)監(jiān)控系統(tǒng)按通用化、系列化和組合化的要求設(shè)計(jì),選用SRI-VC2100標(biāo)準(zhǔn)的控制器,主要由集控臺(tái)、監(jiān)測(cè)數(shù)據(jù)采集箱和車鐘復(fù)示箱組成(見(jiàn)圖3)。該系統(tǒng)主要配置的模塊是遙控控制模塊、安保控制模塊和數(shù)據(jù)測(cè)量模塊等(見(jiàn)表1)。
ke=n1/|emax|=3/0.05=60
(2)
kec=n2/|ecmax|=3/0.03=100
(3)
ku=|Umax|/m=30/4=7.5
(4)
式中:n1,n2,m分別為偏差E、偏差變化率EC和輸出變量U的模糊論域最大值;emax、ecmax、umax分別為偏差E、偏差變化率EC和輸出變量U的基本論域最大值。
隸屬度函數(shù)的選擇遵循的規(guī)則是[11]:在誤差為零附近區(qū)域,采用分辨率較高的隸屬度函數(shù);在誤差較大的區(qū)域,為使系統(tǒng)具有更好的魯棒性能,選擇分辨率較低的隸屬度函數(shù)。針對(duì)實(shí)際情況,為了更好地解決系統(tǒng)靈敏度和穩(wěn)定性相矛盾的關(guān)系,對(duì)濕度偏差E語(yǔ)言值{NB,NM,NS,Z,S,PM,PB}分2種制定規(guī)則,隸屬度函數(shù)NB、NM、PM、PB采用高斯型隸屬度函數(shù),NS、Z、PS采用三角隸屬度函數(shù)。濕度變化率采用三角隸屬度函數(shù),輸出3個(gè)修正參數(shù)值采用三角隸屬度函,具體隸屬度函數(shù)曲線見(jiàn)圖2、圖3。
圖2 E隸屬度函數(shù)曲線
圖3 EC、ΔKp、ΔKi、ΔKd隸屬度函數(shù)曲線
(1)當(dāng)誤差|E|較大時(shí),為使系統(tǒng)具有較好的快速收斂性能,不論誤差的變化趨勢(shì)如何都應(yīng)取較大的Kp和較小的Kd,同時(shí)為了避免積分飽和使得系統(tǒng)響應(yīng)出現(xiàn)較大超調(diào),應(yīng)降低積分作用,Ki取值較小,通常取Ki=0。
(2)當(dāng)誤差|E|和|EC|處于中等大小時(shí),為防止超調(diào)過(guò)大,保證系統(tǒng)的收斂速度,Kp的取值應(yīng)小一些,Ki取值中等大小。如果此時(shí)EEC<0,則應(yīng)取較大的Kd防止系統(tǒng)超調(diào);反之,EEC>0,應(yīng)取適中的Kd,阻止誤差繼續(xù)增大。
(3)當(dāng)誤差|E|較小或?yàn)榱銜r(shí),為保證系統(tǒng)的穩(wěn)態(tài)性能,應(yīng)取較大的Kp、Ki值,為了避免輸出相應(yīng)在設(shè)定值附近出現(xiàn)劇烈震蕩以及增加調(diào)節(jié)時(shí)間,Kd值須要遵循原則是:當(dāng)偏差|EC|值較大時(shí),Kd取值較??;當(dāng)偏差|EC|值較小時(shí),Kd取值較大,通常Kd取值中等大小。
模糊控制的核心是總結(jié)工程操作人員的技術(shù)知識(shí)和實(shí)際操作經(jīng)驗(yàn),制定符合工程項(xiàng)目需求的模糊控制規(guī)則表。根據(jù)以上PID參數(shù)整定規(guī)則可制定模糊PID控制邏輯規(guī)則表,見(jiàn)表1~表3。
表1 ΔKp的模糊控制規(guī)則
表2ΔKi的模糊控制規(guī)則
EECNBNMNSZOPSPMPBNBNBZOPBNMPBPBNBNMNBZOPBNSPBZONBNSNBNSPBNSPBPSNBZONBNBPBZOPBPSNBPSNBNSPBZOPBPMNBPMNBZOPBPSPBZONBPBNBZOPBPSPBPBNB
表3 ΔKd的模糊控制規(guī)則
為了驗(yàn)證所設(shè)計(jì)控制器的有效性,應(yīng)用MATLAB環(huán)境下的Simulink仿真軟件對(duì)溫室灌溉系統(tǒng)進(jìn)行建模與仿真,具體模型見(jiàn)圖4。
圖4 simulink搭建溫室灌溉系統(tǒng)仿真模型
模型中考慮到實(shí)際土壤水勢(shì)曲線為非線性,采用目標(biāo)函數(shù)[13]Δy=sin(2 πu/T),由于土壤濕度從13%變化到23%用時(shí)為30 min,計(jì)算周期T=1 882 min??紤]到土壤濕度值受到多種因素影響,在灌溉過(guò)程中土壤中的水分會(huì)通過(guò)基質(zhì)滲透或植物蒸騰作用而消耗,因此設(shè)定系統(tǒng)的損耗函數(shù)Fun2=exp(-u/500)。溫室灌溉系統(tǒng)是純時(shí)延系統(tǒng),添加一個(gè)Transport Delay模塊,模擬延時(shí)環(huán)節(jié)。通過(guò)單獨(dú)PID控制器仿真分析,當(dāng)PID Controller中的3個(gè)參數(shù)Kp、Ki、Kd的初始值設(shè)為Kp0=15、Ki0=0.5、Kd0=0.08且Kp的變化范圍為[10,20],Ki的變化范圍為[0,1],Kd的變化范圍為[0,0.16]時(shí),系統(tǒng)的控制效果較好。PID控制器的3個(gè)參數(shù)Kp、Ki、Kd及其對(duì)應(yīng)的輸出比例因子ap、ai、ad由以下公式表示:
Kp=ΔKpap+Kp0
(5)
Ki=ΔKiai+Ki0
(6)
Kd=ΔKdad+Kd0
(7)
(8)
(9)
(10)
式中:Kpmax、Kimax、Kdmax分別對(duì)應(yīng)Kp、Ki、Kd變化范圍內(nèi)的最大值;n3、n4、n5分別對(duì)應(yīng)Kp、Ki、Kd模糊論域最大值。
仿真結(jié)果曲線見(jiàn)圖5,系統(tǒng)動(dòng)態(tài)穩(wěn)定局部放大曲線見(jiàn)圖6。
圖5 濕度仿真曲線
圖6 濕度仿真曲線局部放大
從仿真曲線可以看出,本文設(shè)計(jì)的PID模糊控制器能夠使目標(biāo)濕度和控制濕度的誤差穩(wěn)定 以內(nèi),與單一的PID控制器相比,誤差精度提高了2%左右;與單一的而模糊控制器相比,仿真曲線收斂速度快了20 s,且超調(diào)量減少了5%。仿真結(jié)果達(dá)到預(yù)期,系統(tǒng)設(shè)計(jì)符合溫室灌溉要求。
通過(guò)比較模糊控制器和PID控制器的優(yōu)缺點(diǎn),綜合設(shè)計(jì)模糊PID控制器,揚(yáng)長(zhǎng)避短,用Simulink搭建溫室灌溉控制系統(tǒng),考慮影響植物生長(zhǎng)的主要因素,控制土壤濕度穩(wěn)定在植物生長(zhǎng)所需最佳土壤濕度誤差范圍內(nèi)。仿真曲線表明控制誤差可以穩(wěn)定在5%范圍內(nèi),超調(diào)量為6.7%,仿真曲線收斂時(shí)長(zhǎng)80 s,能夠達(dá)到溫室節(jié)水灌溉要求。該控制器對(duì)非線性,純時(shí)延系統(tǒng)有很好的魯棒性,穩(wěn)態(tài)精度高,控制效果好,對(duì)實(shí)驗(yàn)平臺(tái)的搭建、灌溉控制器策略的選擇有一定的借鑒意義。
:
[1] 鄧 忠,翟國(guó)亮,呂謀超,等. 我國(guó)農(nóng)業(yè)應(yīng)對(duì)干旱災(zāi)害的技術(shù)研究現(xiàn)狀及展望[J]. 節(jié)水灌溉,2016,(8):162-165.
[2] 徐文靜,王翔翔,施六林,等. 中國(guó)節(jié)水灌溉技術(shù)現(xiàn)狀與發(fā)展趨勢(shì)研究[J]. 中國(guó)農(nóng)學(xué)通報(bào),2016,32(11):184-187.
[3] 杜云明,蓋麗娜,顏兵兵. 基于單片機(jī)的溫室灌溉控制系統(tǒng)設(shè)計(jì)[J]. 農(nóng)機(jī)化究,2012,34 (12) 88-91.
[4] Zhu Zhijian, Zao Remu, Ni Jiaoti.Autocontrol variable voltage frequency type of irrigation Feri- gation device[J].Transaction of the CSAE,2005,21(9):94-97.
[5] 田思慶,曹 宇,魏 強(qiáng),等. 基于模糊控制的智能滴灌控制系統(tǒng)設(shè)計(jì)[J]. 節(jié)水灌溉,2017,(6):101-104.
[6] 謝宋和,甘 勇. 單片機(jī)模糊控制系統(tǒng)設(shè)計(jì)與應(yīng)用實(shí)例[M]. 北京: 電子工業(yè)出版社, 1999.
[7] 張育斌,魏正英,馬勝利,等.灰色預(yù)測(cè)模糊PID灌溉控制技術(shù)開(kāi)發(fā)[J].中國(guó)農(nóng)村水利水電,2016,(2):5-8.
[8] 謝仕宏.MATLAB控制系統(tǒng)動(dòng)態(tài)仿真實(shí)例教程[M]. 北京:化學(xué)工業(yè)出版社, 2008.
[9] 陳年來(lái),黃海霞,高慧娟,等.甜瓜葉片氣體交換特性和幼苗生長(zhǎng)對(duì)土壤水分和大氣濕度的響應(yīng)[J].蘭州大學(xué)學(xué)報(bào)(自然科學(xué)版),2009,45(4):73-77.
[10] 孫增圻,鄧志東,張?jiān)倥d.智能控制理論與技術(shù)[M].北京:清華大學(xué)出版社,2011.
[11] 梁月云,崔天時(shí),何亞非. 溫室節(jié)水灌溉系統(tǒng)模糊控制器設(shè)計(jì)及MATLAB仿真[J]. 農(nóng)機(jī)化研究,2014,36(6):202-205,209.
[12] 李 健,王冬青,王麗美. 模糊PID控制器設(shè)計(jì)及MATLAB仿真[J]. 工業(yè)控制計(jì)算機(jī),2011,24(5):56-57,95.
[13] 徐立鴻. 一種新型溫室滴灌控制算法[C]∥ 中國(guó)農(nóng)業(yè)工程學(xué)會(huì).中國(guó)農(nóng)業(yè)工程學(xué)會(huì)第七次全國(guó)會(huì)員代表大會(huì)及學(xué)術(shù)年會(huì)論文集.中國(guó)農(nóng)業(yè)工程學(xué)會(huì),2004:3.