• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    采用多源信息融合的妊娠豬舍環(huán)境質(zhì)量評價方法

    2022-02-07 01:07:32郭艷嬌鄭永軍
    農(nóng)業(yè)工程學(xué)報 2022年18期
    關(guān)鍵詞:融合環(huán)境評價

    遲 宇,郭艷嬌,馮 涵,李 寒,鄭永軍,3

    采用多源信息融合的妊娠豬舍環(huán)境質(zhì)量評價方法

    遲 宇1,郭艷嬌1,馮 涵1,李 寒2,鄭永軍1,3※

    (1. 中國農(nóng)業(yè)大學(xué)工學(xué)院,北京 100083;2. 中國農(nóng)業(yè)大學(xué)信息與電氣工程學(xué)院,北京 100083;3. 現(xiàn)代農(nóng)業(yè)裝備與設(shè)施教育部工程研究中心,北京 100083)

    妊娠豬舍作為養(yǎng)殖場豬只繁育的基礎(chǔ)條件,其環(huán)境質(zhì)量對母豬的生產(chǎn)性能有顯著影響。為合理評價妊娠豬舍環(huán)境質(zhì)量,該研究提出一種基于模擬退火的粒子群算法(Simulated Annealing-Particle Swarm Optimization,SA-PSO)、套索算法(Least Absolute Shrinkage and Selection Operator,LASSO)和反向傳播(Back Propagation,BP)神經(jīng)網(wǎng)絡(luò)的環(huán)境質(zhì)量評價模型。利用卡爾曼濾波和分批估計自適應(yīng)加權(quán)融合算法,實現(xiàn)多節(jié)點環(huán)境數(shù)據(jù)的時間與空間序列融合;構(gòu)建豬舍環(huán)境質(zhì)量非線性評價模型,采用LASSO算法,篩選得出與環(huán)境質(zhì)量強(qiáng)相關(guān)的特征參數(shù),實現(xiàn)輸入降維;融合SA-PSO算法實現(xiàn)網(wǎng)絡(luò)初始權(quán)值和閾值的優(yōu)化,形成SA-PSO-LASSO-BP神經(jīng)網(wǎng)絡(luò)評價模型。通過對數(shù)據(jù)采集系統(tǒng)獲取的實際妊娠豬舍環(huán)境數(shù)據(jù)進(jìn)行驗證,結(jié)果表明:提出的環(huán)境質(zhì)量評價模型決定系數(shù)為0.918、總準(zhǔn)確率為95.85%,相比單純使用BP神經(jīng)網(wǎng)絡(luò),加入LASSO和SA-PSO算法后決定系數(shù)與總準(zhǔn)確率分別提高了37.43%、11.09個百分點,具有更高的評價精度和性能,可更好地擬合復(fù)雜環(huán)境參數(shù)與環(huán)境質(zhì)量間的非線性關(guān)系,為妊娠豬舍環(huán)境質(zhì)量評價提供參考。

    模型;環(huán)境;妊娠豬舍;環(huán)境質(zhì)量;BP神經(jīng)網(wǎng)絡(luò);LASSO算法;SA-PSO算法

    0 引 言

    在生豬養(yǎng)殖集約化、規(guī)模化發(fā)展的過程中,母豬的繁殖性能和健康作為影響?zhàn)B殖企業(yè)經(jīng)濟(jì)效益的關(guān)鍵因素,受到妊娠豬舍內(nèi)環(huán)境質(zhì)量(適宜性或舒適度)的直接影響[1-3]。因此建立可監(jiān)測特征環(huán)境參數(shù)與舒適度之間的關(guān)系、及時準(zhǔn)確地評價妊娠豬舍環(huán)境質(zhì)量,是改善環(huán)境調(diào)控措施、實現(xiàn)豬舍環(huán)境舒適度精準(zhǔn)調(diào)控的重要前提,對于減輕母豬生長與繁育過程中受到的環(huán)境脅迫、提高養(yǎng)殖效益有實際意義。

    目前,對于豬舍環(huán)境質(zhì)量評價主要分為2類:模糊綜合評價模型與機(jī)器學(xué)習(xí)模型。由于豬舍環(huán)境是一個由多環(huán)境參數(shù)共同作用所形成的復(fù)雜非線性時變系統(tǒng),且環(huán)境參數(shù)間相互影響,因此舍內(nèi)環(huán)境質(zhì)量與環(huán)境參數(shù)間的關(guān)系難以使用準(zhǔn)確的數(shù)學(xué)模型進(jìn)行映射。模糊綜合評價法[4-6]作為豬舍環(huán)境質(zhì)量常用的評價方法,雖然能夠有效解決多環(huán)境參數(shù)的不確定性和模糊邊界問題,但未能挖掘環(huán)境參數(shù)間的耦合關(guān)系和環(huán)境參數(shù)與環(huán)境質(zhì)量間的非線性關(guān)系,存在欠學(xué)習(xí)的問題。由此,程捷等[7]提出一種基于改進(jìn)D-S(Dempster-Shafer)證據(jù)理論的數(shù)據(jù)融合算法,能夠兼顧各環(huán)境參數(shù)的數(shù)據(jù)特征和內(nèi)在聯(lián)系,有效地識別豬舍環(huán)境狀態(tài)??紤]到環(huán)境參數(shù)與環(huán)境適宜性之間的非線性關(guān)系,孫聰[8]利用改進(jìn)CS(Cuckoo Search)算法優(yōu)化的BP(Back Propagation)神經(jīng)網(wǎng)絡(luò)構(gòu)建豬舍環(huán)境適宜度評價模型;朱保釧[9]通過改進(jìn)的C4.5算法實現(xiàn)對育肥豬舍環(huán)境適宜性的預(yù)測;陳沖等[10]提出了一種基于MSCCS-LSSVR(Mutative Scale Chaos Cuckoo Search-Least Square Support Vector Regression)模型的哺乳母豬舍環(huán)境舒適度評價預(yù)測方法。由于豬舍內(nèi)環(huán)境受到舍外環(huán)境、地理位置等因素影響,不同豬舍內(nèi)環(huán)境參數(shù)的變化范圍和變化趨勢存在實際差異。以上方法雖能擬合豬舍環(huán)境參數(shù)與適宜性之間的非線性關(guān)系,但未對參數(shù)選取原則與其對環(huán)境質(zhì)量的影響程度作深入分析,模型存在應(yīng)用普適性和泛化性差的問題。

    BP神經(jīng)網(wǎng)絡(luò)是一種多層前饋型網(wǎng)絡(luò)模型,具有高度的自學(xué)習(xí)和非線性映射能力,已在畜禽舍內(nèi)環(huán)境參數(shù)評估預(yù)測方面取得良好的效果[11-12]。因此本文采用BP神經(jīng)網(wǎng)絡(luò)構(gòu)建非線性評價模型,但其性能易受初始權(quán)值和閾值隨機(jī)性的影響,仍需采用優(yōu)化算法對其初始參數(shù)進(jìn)行優(yōu)化[12-13],進(jìn)一步提升網(wǎng)絡(luò)性能和預(yù)測精度。針對妊娠豬舍內(nèi)環(huán)境質(zhì)量的特征參數(shù)選擇與BP神經(jīng)網(wǎng)絡(luò)參數(shù)優(yōu)化問題,采用LASSO(Least Absolute Shrinkage and Selection Operator)算法篩選特征參數(shù),簡化多環(huán)境參數(shù)結(jié)構(gòu);利用SA-PSO(Simulated Annealing-Particle Swarm Optimization)算法優(yōu)化BP神經(jīng)網(wǎng)絡(luò)的初始權(quán)值和閾值,構(gòu)建SA-PSO-LASSO-BP神經(jīng)網(wǎng)絡(luò)模型,實現(xiàn)妊娠豬舍環(huán)境質(zhì)量的可靠評價,為后續(xù)環(huán)境質(zhì)量精準(zhǔn)調(diào)控提供決策支持。

    1 數(shù)據(jù)獲取

    1.1 采集設(shè)備與采集方案

    采集現(xiàn)場為湖南省株洲某生豬養(yǎng)殖基地。該基地是一個典型的集生豬培育、種豬繁殖于一體的集約化生豬養(yǎng)殖基地,本研究所選豬舍為密閉環(huán)境,長40 m、寬14 m、高3 m,內(nèi)部為4列式豬欄結(jié)構(gòu),如圖1a所示。該豬舍內(nèi)限位欄長2.2 m、寬0.65 m、高1 m,共養(yǎng)殖228頭長白豬。依據(jù)國家標(biāo)準(zhǔn)[14],豬舍環(huán)境評價主要包括豬舍空氣、豬舍通風(fēng)、豬舍采光、豬舍噪聲等4個方面,且每方面均包含多個環(huán)境參數(shù)。其中溫度、相對濕度、NH3質(zhì)量濃度、CO2質(zhì)量濃度、H2S質(zhì)量濃度、PM2.5質(zhì)量濃度、PM10質(zhì)量濃度、光照強(qiáng)度、噪聲等9個指標(biāo)是豬舍環(huán)境質(zhì)量的主要影響因素[15-17],均可通過傳感器直接獲得。經(jīng)預(yù)檢測,將傳感器安裝于距離地面2 m處,數(shù)值波動較小且可避免豬舍噴水、清潔等管理工作對傳感器造成的損壞。試驗采用的多種傳感器組成1個采集節(jié)點,布置見圖1b,實現(xiàn)數(shù)據(jù)的可持續(xù)監(jiān)測。各傳感器及參數(shù)指標(biāo)見表1。

    圖1 豬舍結(jié)構(gòu)及采集節(jié)點布置

    表1 傳感器參數(shù)

    利用Arduino Mega 2560作為控制核心建立數(shù)據(jù)采集系統(tǒng),總體結(jié)構(gòu)如圖2所示。該系統(tǒng)的10個采集節(jié)點每間隔10 min對環(huán)境數(shù)據(jù)進(jìn)行1次采集,并通過WIFI模塊將各傳感器采集的環(huán)境數(shù)據(jù)傳輸至物聯(lián)網(wǎng)監(jiān)測平臺數(shù)據(jù)中心,實現(xiàn)數(shù)據(jù)的存儲與遠(yuǎn)程監(jiān)控。采集周期為2021年10月1日至2021年12月31日(92 d)。

    圖2 數(shù)據(jù)采集系統(tǒng)

    1.2 數(shù)據(jù)預(yù)處理

    1.2.1 單一傳感器時間序列數(shù)據(jù)融合

    由于環(huán)境參數(shù)傳感器在工作過程中易受到妊娠豬舍日常生產(chǎn)與管理工作的干擾而產(chǎn)生噪聲。采用卡爾曼濾波算法[17]對原始數(shù)據(jù)進(jìn)行濾波處理,能夠有效地抑制系統(tǒng)與環(huán)境噪聲,減小隨機(jī)誤差的影響,較好地實現(xiàn)時間序列數(shù)據(jù)融合。以下以溫度傳感器數(shù)據(jù)為例,說明數(shù)據(jù)預(yù)處理過程。

    本文對2021年10月1日至2日豬舍內(nèi)一個溫度傳感器連續(xù)采集的288個數(shù)據(jù)進(jìn)行預(yù)處理。為填補(bǔ)無線傳輸丟包造成的缺失數(shù)據(jù),將同一時刻其他同質(zhì)傳感器采集數(shù)據(jù)的平均值作為本傳感器缺失的時點數(shù)據(jù)。采集到的原始溫度時序數(shù)據(jù)經(jīng)卡爾曼濾波平滑后消除了峰值突變(圖3),實現(xiàn)異常數(shù)據(jù)修正與環(huán)境噪聲實時動態(tài)計算和補(bǔ)償。濾波前后溫度數(shù)據(jù)的方差分別為1.49和1.43,濾波后數(shù)據(jù)的方差較濾波前下降了4.03%,數(shù)據(jù)波動性降低,確保了數(shù)據(jù)的準(zhǔn)確性和有效性。

    圖3 卡爾曼濾波結(jié)果

    1.2.2 多源同質(zhì)傳感器空間序列數(shù)據(jù)融合

    單一傳感器時間序列數(shù)據(jù)融合后,需要對10個采集節(jié)點的同質(zhì)傳感器數(shù)據(jù)進(jìn)行空間序列融合。本研究以分批估計自適應(yīng)加權(quán)融合算法[18]為核心,首先利用等分法將單一傳感器連續(xù)采集的10個數(shù)據(jù)劃分為2組,根據(jù)分批估計理論分別求取2組數(shù)據(jù)的局部方差;再采用自適應(yīng)加權(quán)融合算法計算多源同質(zhì)傳感器數(shù)據(jù)的整體方差,獲得各傳感器數(shù)據(jù)的加權(quán)因子,從而實現(xiàn)多源同質(zhì)傳感器空間序列融合。

    為驗證分批估計自適應(yīng)加權(quán)融合算法的可行性,首先分別對妊娠豬舍內(nèi)布置的10個采集節(jié)點獲得的原始溫度數(shù)據(jù)進(jìn)行卡爾曼濾波,再對同一時刻不同溫度傳感器數(shù)據(jù)進(jìn)行空間序列融合,融合數(shù)據(jù)的均方誤差(Mean Square Error,MSE)與算術(shù)平均法和自適應(yīng)加權(quán)融合算法對比,結(jié)果如圖4所示。

    注:本文算法為分批估計自適應(yīng)加權(quán)融合算法。

    由圖4可知,利用算術(shù)平均法對傳感器數(shù)據(jù)進(jìn)行空間融合,MSE分布于0.455~2.358,使用自適用應(yīng)權(quán)融合算法融合數(shù)據(jù)的MSE分布于0.006~0.141,而分批估計自適應(yīng)加權(quán)融合算法融合數(shù)據(jù)的MSE最大值僅為0.055。本文采用的分批估計自適應(yīng)加權(quán)融合算法效果最優(yōu),融合數(shù)據(jù)的MSE較算術(shù)平均法下降了92.14%以上,有效提高了融合精度,減少了數(shù)據(jù)冗余。

    按以上缺失時點數(shù)據(jù)同質(zhì)均值補(bǔ)償、卡爾曼濾波、分批估計自適應(yīng)加權(quán)等流程對所有傳感器的數(shù)據(jù)融合預(yù)處理后,進(jìn)行歸一化處理,統(tǒng)一于[0,1]區(qū)間,消除各項參數(shù)的量綱與數(shù)值差異。

    2 環(huán)境質(zhì)量評價模型

    由于豬舍環(huán)境是一個由多種環(huán)境參數(shù)共同作用所形成的時變系統(tǒng),為挖掘環(huán)境參數(shù)間的耦合關(guān)系、擬合環(huán)境參數(shù)與環(huán)境質(zhì)量之間的非線性關(guān)系,本研究采用BP神經(jīng)網(wǎng)絡(luò)對環(huán)境質(zhì)量進(jìn)行評價。通過LASSO算法對特征環(huán)境參數(shù)進(jìn)行篩選,重置BP神經(jīng)網(wǎng)絡(luò)的結(jié)構(gòu)參數(shù),確定BP神經(jīng)網(wǎng)絡(luò)的拓?fù)浣Y(jié)構(gòu),進(jìn)一步提高環(huán)境質(zhì)量評價精度;利用SA-PSO算法優(yōu)化BP神經(jīng)網(wǎng)絡(luò)的初始權(quán)值和閾值,提高BP神經(jīng)網(wǎng)絡(luò)的收斂速度與收斂精度,最終形成的SA-PSO-LASSO-BP神經(jīng)網(wǎng)絡(luò)算法流程如圖 5所示。

    圖5 SA-PSO-LASSO-BP神經(jīng)網(wǎng)絡(luò)算法流程圖

    2.1 評價標(biāo)準(zhǔn)建立

    參考國家標(biāo)準(zhǔn)GB/T 17824.3—2008《規(guī)模豬場環(huán)境參數(shù)及環(huán)境管理》,將環(huán)境質(zhì)量識別框架分為3級:1級(適宜)、2級(較適宜)、3級(不適宜),并確定9種環(huán)境因素在3種環(huán)境質(zhì)量識別框架下的特征值范圍,構(gòu)建妊娠豬舍環(huán)境質(zhì)量評價標(biāo)準(zhǔn)(表2)。

    表2 環(huán)境質(zhì)量評價標(biāo)準(zhǔn)

    注:其他情況為9種環(huán)境參數(shù)處于不同環(huán)境質(zhì)量評價等級特征值范圍。

    Note: Other conditions indicate that nine environmental parameters are in the range of characteristic values of different environmental quality evaluation grades.

    當(dāng)9種環(huán)境參數(shù)均滿足同一等級的特征值要求時,可直接確定其對應(yīng)的等級。對同一時刻環(huán)境參數(shù)處于不同等級特征值范圍的情況,當(dāng)其中某一項環(huán)境參數(shù)嚴(yán)重超標(biāo)時,采用“一項否決”策略評價為3級;否則,通過層次分析法和梯形隸屬度函數(shù)計算確定豬舍環(huán)境質(zhì)量等級[4]。

    基于數(shù)據(jù)采集系統(tǒng)獲得的92 d數(shù)據(jù),預(yù)處理后選取3種環(huán)境質(zhì)量識別框架下的傳感器數(shù)據(jù)共5 110組。采用上述評價標(biāo)準(zhǔn)進(jìn)行環(huán)境質(zhì)量等級評定,最終獲得1等級樣本數(shù)據(jù)1 566組,2等級樣本數(shù)據(jù)1 785組,3等級樣本數(shù)據(jù)1 759組,可作為特征選擇、環(huán)境質(zhì)量評價模型構(gòu)建優(yōu)化的數(shù)據(jù)訓(xùn)練集與測試集。

    2.2 特征選擇

    在采用BP神經(jīng)網(wǎng)絡(luò)進(jìn)行豬舍環(huán)境質(zhì)量評價時,若將9種環(huán)境參數(shù)均作為輸入量會增加模型結(jié)構(gòu)的復(fù)雜性,降低模型的性能和泛化能力;若輸入特征參數(shù)不足,則會降低模型的評價精度。因此,本文采用LASSO算法[19]對9種環(huán)境參數(shù)進(jìn)行特征選擇,優(yōu)化輸入?yún)?shù)量。該方法在回歸系數(shù)縮減的過程中,將特征性較弱的環(huán)境參數(shù)的回歸系數(shù)縮減為0,從而達(dá)到“變量選擇”的目的,其目標(biāo)函數(shù)為

    式中為樣本數(shù)量;為環(huán)境因素個數(shù);y為第個樣本對應(yīng)的環(huán)境質(zhì)量等級;β為第個環(huán)境因素對應(yīng)的回歸系數(shù);x為第個樣本對應(yīng)的第個影響因素;為懲罰因子。

    隨機(jī)從預(yù)處理后的各等級環(huán)境數(shù)據(jù)中抽取700組數(shù)據(jù),共獲得2 100組樣本數(shù)據(jù)。為提高模型的泛化能力和特征選擇的準(zhǔn)確性,采用十折交叉驗證,按照9∶1的比例從2 100組樣本數(shù)據(jù)中每種等級隨機(jī)抽取數(shù)據(jù),共獲得1 890組訓(xùn)練集和210組測試集,求解得到不同值下對應(yīng)的均方根誤差(Root Mean Squared Error,RMSE)如圖6所示。由圖6可知,當(dāng)處于0.1~0.2時,均方根誤差存在最小值。考慮交叉驗證是隨機(jī)分組,為保證模型的泛化性和魯棒性,選取10次交叉驗證結(jié)果的平均值作為LASSO算法最佳懲罰因子,即0.147。

    圖6 不同λ取值下的均方根誤差

    圖7為LASSO算法分析后各環(huán)境參數(shù)所對應(yīng)的回歸系數(shù)軌跡圖,由圖7可知,隨著值的增加,環(huán)境參數(shù)的標(biāo)準(zhǔn)化回歸系數(shù)逐漸收斂。當(dāng)= 0.147時,溫度、相對濕度、NH3質(zhì)量濃度、CO2質(zhì)量濃度的標(biāo)準(zhǔn)化回歸系數(shù)分別為0.57、-0.097、0.138、0.185,H2S質(zhì)量濃度、PM2.5質(zhì)量濃度、PM10質(zhì)量濃度、光照強(qiáng)度、噪聲的標(biāo)準(zhǔn)化回歸系數(shù)均為0,在特征選擇的過程中被淘汰。因此,本文選取溫度、相對濕度、NH3質(zhì)量濃度、CO2質(zhì)量濃度4種環(huán)境參數(shù)作為BP神經(jīng)網(wǎng)絡(luò)模型的輸入量。

    圖7 9種環(huán)境參數(shù)的標(biāo)準(zhǔn)化回歸系數(shù)

    2.3 特征分析

    妊娠豬舍2021年11月24日—2021年11月30日特征環(huán)境參數(shù)的分布如圖8所示。由圖可知妊娠豬舍溫度于14:00—15:00達(dá)到峰值,然后隨時間逐漸下降,在08:00—09:00達(dá)到最低,然后逐漸升高;相反,相對濕度、CO2質(zhì)量濃度和NH3質(zhì)量濃度于08:00—09:00達(dá)到峰值而后逐漸下降,于14:00—15:00降至最低值后波動上升。妊娠豬舍內(nèi)溫度、相對濕度、NH3質(zhì)量濃度、CO2質(zhì)量濃度呈現(xiàn)出明顯的非線性變化趨勢,且具有一定的周期性和相關(guān)性。

    對特征環(huán)境參數(shù)進(jìn)行Pearson相關(guān)性分析,分析結(jié)果如表3所示。由表3可知,溫度與相對濕度、NH3質(zhì)量濃度、CO2質(zhì)量濃度之間的相關(guān)系數(shù)分別為-0.637、-0.871、-0.737,說明溫度與其他特征環(huán)境參數(shù)間存在顯著負(fù)相關(guān)關(guān)系(<0.01);此外,相對濕度與NH3質(zhì)量濃度、CO2質(zhì)量濃度之間的相關(guān)系數(shù)為0.330、0.723,表明相對濕度與NH3質(zhì)量濃度、CO2質(zhì)量濃度間存在顯著正相關(guān)關(guān)系(<0.01),與圖8中各環(huán)境參數(shù)的變化規(guī)律一致。這是因為溫度的上升加速了舍內(nèi)水分的蒸發(fā),從而導(dǎo)致相對濕度下降;隨著溫度的持續(xù)升高,豬在行為上表現(xiàn)出精神萎靡、呼吸深度變淺、采食量減退等現(xiàn)象[20],機(jī)體的排泄量也隨之減少[21],舍內(nèi)CO2質(zhì)量濃度和NH3質(zhì)量濃度降低。因此,溫度、相對濕度、NH3質(zhì)量濃度、CO2質(zhì)量濃度呈現(xiàn)出此種變化趨勢和相關(guān)關(guān)系。

    由于豬舍特征環(huán)境參數(shù)呈非線性變化趨勢,且各參數(shù)間存在相關(guān)關(guān)系,因此難以使用準(zhǔn)確的數(shù)學(xué)模型描述環(huán)境參數(shù)與環(huán)境質(zhì)量之間的關(guān)系。BP神經(jīng)網(wǎng)絡(luò)作為一種多層前饋型網(wǎng)絡(luò)模型,具有高度的自學(xué)習(xí)、自適應(yīng)和泛化能力,較強(qiáng)的非線性映射能力[22],能夠充分挖掘環(huán)境參數(shù)間的耦合關(guān)系,擬合復(fù)雜環(huán)境參數(shù)與環(huán)境質(zhì)量之間的非線性關(guān)系,從而達(dá)到良好的識別與分類目的。因此,本研究使用BP神經(jīng)網(wǎng)絡(luò)進(jìn)行建模,實現(xiàn)妊娠豬舍環(huán)境質(zhì)量的可靠性評價。

    注:7d,2021-11-24至2021-11-30;1d,2021-11-24。

    表3 特征環(huán)境參數(shù)間的相關(guān)系數(shù)

    注(Note): *,<0.01.

    2.4 BP神經(jīng)網(wǎng)絡(luò)模型結(jié)構(gòu)的確定

    BP神經(jīng)網(wǎng)絡(luò)由輸入層、隱含層和輸出層組成,其中輸入層為特征選擇后的4種環(huán)境參數(shù)(4個節(jié)點),輸出層為對應(yīng)的環(huán)境質(zhì)量等級(1個節(jié)點),隱含層的節(jié)點數(shù)需進(jìn)行求解。目前,沒有固定的方法求解隱含層節(jié)點數(shù),但可依據(jù)Kolmogorov定理進(jìn)行隱含層節(jié)點數(shù)的經(jīng)驗估計,其經(jīng)驗公式[12, 23]為

    式中為輸入層節(jié)點數(shù);為輸出層節(jié)點數(shù);為常數(shù),取值范圍為1~10。

    在本研究中,= 4,= 1,最終計算得到隱含層節(jié)點數(shù)的取值范圍為4~12。

    為確定隱含層節(jié)點具體數(shù)量,在I5-7300 HQ、2.5 GHz、8 GB內(nèi)存、Window 10系統(tǒng)集成開發(fā)環(huán)境下,使用MATLAB編寫相關(guān)程序。經(jīng)反復(fù)試驗確定輸入層到隱含層的傳遞函數(shù)為“tansig”,隱含層到輸出層的傳遞函數(shù)為“purelin”“trainlm”作為常用的訓(xùn)練函數(shù)之一,具有收斂速度快、誤差小、訓(xùn)練效果優(yōu)等特點[23],因此本文選取“trainlm”為訓(xùn)練函數(shù);同時,設(shè)置BP神經(jīng)網(wǎng)絡(luò)的迭代次數(shù)為1 000,學(xué)習(xí)率為0.01。在此基礎(chǔ)上,將2100組環(huán)境數(shù)據(jù)按照9∶1的比例劃分為訓(xùn)練集和測試集后,分別訓(xùn)練具有4~12個隱含層節(jié)點個數(shù)的BP神經(jīng)網(wǎng)絡(luò)。

    為了對不同隱含層節(jié)點數(shù)量下模型的預(yù)測結(jié)果進(jìn)行定量分析,采用平均絕對誤差(Mean Absolute Error,MAE)、RMSE、決定系數(shù)(Coefficient of determination,2)、準(zhǔn)確率(Accuracy,ACC)作為評價指標(biāo),不同隱含層節(jié)點數(shù)量下BP神經(jīng)網(wǎng)絡(luò)的性能如圖9所示。由圖9可知,當(dāng)隱含層節(jié)點數(shù)為11時,BP神經(jīng)網(wǎng)絡(luò)的MAE、RMSE、2、ACC分別達(dá)到0.077、0.277、0.878、92.33%,均為最優(yōu)值。相較于其他模型,隱含層節(jié)點數(shù)為11的BP神經(jīng)網(wǎng)絡(luò)2與ACC至少提高了5.73%、2.50%,MAE與RMSE至少下降了3.75%、2.12%,具有更高的性能和評價精度。因此,本研究確定最佳隱含層節(jié)點數(shù)=11。

    注:R2,決定系數(shù);ACC,準(zhǔn)確率,%。

    2.5 BP神經(jīng)網(wǎng)絡(luò)模型參數(shù)的優(yōu)化

    基于單一BP神經(jīng)網(wǎng)絡(luò)對豬舍環(huán)境質(zhì)量進(jìn)行評價時,其性能受到網(wǎng)絡(luò)初始權(quán)值和閾值等參數(shù)的影響,存在收斂速度慢、易陷入局部極小值等問題,從而無法尋得全局最優(yōu)值。為提高妊娠豬舍環(huán)境質(zhì)量評價精度,本研究采用粒子群算法(Particle Swarm Optimization,PSO)對BP神經(jīng)網(wǎng)絡(luò)的初始權(quán)值和閾值進(jìn)行優(yōu)化,既保留神經(jīng)網(wǎng)絡(luò)的非線性擬合能力,又能夠發(fā)揮粒子群算法的全局尋優(yōu)能力,綜合提升BP神經(jīng)網(wǎng)絡(luò)的收斂速度和尋優(yōu)精度[24]。

    PSO算法雖然收斂速度快,但存在易過早收斂、易陷入局部最優(yōu)等缺點[25],因此本文將模擬SA算法嵌入PSO算法中。參考文獻(xiàn)[25],在SA-PSO算法中,采用帶壓縮因子的方法選擇合適的參數(shù)提高PSO算法的收斂性;同時對PSO算法的最優(yōu)解進(jìn)行優(yōu)化,讓粒子有較大的概率接受非最優(yōu)解,從而跳出局部最優(yōu)解。假設(shè)在維空間內(nèi),粒子的種群規(guī)模為,第個粒子在第維空間的位置、速度分別為xv,在算法迭代過程中,每個粒子根據(jù)個體極值P和全局極值G來更新自己的速度和位置

    式中v+1為粒子第+1次的速度;x+1為粒子第+1次的位置;v為粒子第次的速度;x為粒子第次的位置; P為粒子第次的個體極值;G’為粒子第次的全局極值;為慣性權(quán)重;1、2為學(xué)習(xí)因子;1、2為隨機(jī)數(shù),1、2∈[0,1];為當(dāng)前迭代次數(shù);為壓縮因子,max為最大慣性權(quán)值,min為最小慣性權(quán)值,為最大迭代次數(shù)。

    以上述特征選擇和環(huán)境質(zhì)量評價方法為基礎(chǔ),共選取3 010組環(huán)境數(shù)據(jù),1級、2級、3級3種質(zhì)量等級下的傳感器數(shù)據(jù)分別為866、1 085、1 059組,按照9∶1的比例劃分成2 709組訓(xùn)練集和301組測試集。經(jīng)過多次試驗測試,SA-PSO算法的初始參數(shù)設(shè)置如下:粒子的種群規(guī)模= 20,學(xué)習(xí)因子1=2= 2.05,最大慣性權(quán)值max= 0.9,最小慣性權(quán)值min= 0.4,退溫系數(shù)= 0.99,最大迭代次數(shù)= 1 000。

    為驗證SA-PSO算法的優(yōu)化性能,選取遺傳算法(Genetic Algorithm,GA)、麻雀搜索算法(Sparrow Search Algorithm,SSA)、粒子群算法3種優(yōu)化算法進(jìn)行對比。試驗以RMSE為適應(yīng)度函數(shù),迭代終止條件為達(dá)到最大迭代次數(shù),4種算法的適應(yīng)度函數(shù)曲線如圖10所示。由圖10可知,隨著迭代次數(shù)的增加,適應(yīng)度函數(shù)曲線收斂。相較于GA、SSA算法,PSO算法的適應(yīng)度值更低,為0.194,尋優(yōu)效果更佳;加入SA算法后,SA-PSO算法的適應(yīng)度值進(jìn)一步降低,達(dá)到0.187,較PSO算法減小了3.61%。因此,通過SA算法對PSO算法進(jìn)行優(yōu)化,能夠更好地發(fā)揮PSO算法的全局尋優(yōu)能力,提高其尋優(yōu)精度,適合于BP神經(jīng)網(wǎng)絡(luò)初始權(quán)值和閾值的優(yōu)化。

    注:GA,遺傳算法;SSA,麻雀搜索算法;PSO,粒子群算法;SA-PSO,模擬退火-粒子群算法。

    3 結(jié)果與分析

    為驗證SA-PSO-LASSO-BP神經(jīng)網(wǎng)絡(luò)的性能,選取BP神經(jīng)網(wǎng)絡(luò)、LASSO-BP神經(jīng)網(wǎng)絡(luò)、GA-LASSO-BP神經(jīng)網(wǎng)絡(luò)、SSA-LASSO-BP神經(jīng)網(wǎng)絡(luò)、PSO-LASSO-BP神經(jīng)網(wǎng)絡(luò)等5種模型進(jìn)行對比。試驗采用MSE作為網(wǎng)絡(luò)誤差函數(shù),迭代終止條件網(wǎng)絡(luò)誤差目標(biāo)值設(shè)置為0.01,6種模型的3 010組訓(xùn)練集性能測試結(jié)果如表4所示。

    表4 不同神經(jīng)網(wǎng)絡(luò)模型性能

    由表4可知,對于收斂速度,BP神經(jīng)網(wǎng)絡(luò)迭代766次時完成訓(xùn)練,而LASSO-BP神經(jīng)網(wǎng)絡(luò)模型迭代537次完成訓(xùn)練,迭代次數(shù)下降了29.90%,說明LASSO算法能夠降低BP神經(jīng)網(wǎng)絡(luò)結(jié)構(gòu)的復(fù)雜性和訓(xùn)練壓力,從而提高網(wǎng)絡(luò)的收斂速度;通過4種優(yōu)化算法改進(jìn)的LASSO-BP神經(jīng)網(wǎng)絡(luò)的迭代次數(shù)分別為374、192、79、36次,運(yùn)行時間分別為2.65、1.93、1.45、1.20 s,其中建立的SA-PSO-LASSO-BP神經(jīng)網(wǎng)絡(luò)的迭代次數(shù)最少、運(yùn)行時間最短,相較于PSO-LASSO-BP神經(jīng)網(wǎng)絡(luò)迭代次數(shù)減少了54.43%、運(yùn)行時間加快了17.24%,具有更快的收斂速度。

    對于收斂精度,BP神經(jīng)網(wǎng)絡(luò)在誤差為0.090時停止訓(xùn)練,LASSO-BP神經(jīng)網(wǎng)絡(luò)在誤差為0.074時停止訓(xùn)練,網(wǎng)絡(luò)誤差下降了17.78%,驗證了LASSO算法能夠提高BP神經(jīng)網(wǎng)絡(luò)輸入?yún)?shù)的特征性,提高網(wǎng)絡(luò)的收斂精度;通過4種優(yōu)化算法改進(jìn)的LASSO-BP神經(jīng)網(wǎng)絡(luò)的誤差值分別為0.040、0.037、0.036、0.033,本文所提出模型的誤差值最低,較PSO-LASSO-BP神經(jīng)網(wǎng)絡(luò)誤差下降了8.33%,具有更高的評價精度和性能。

    模型訓(xùn)練后,以301組樣本測試6種神經(jīng)網(wǎng)絡(luò)模型的評價效果,結(jié)果如表5所示。由表5可知,SA-PSO-LASSO-BP神經(jīng)網(wǎng)絡(luò)的MAE、RMSE、2、總準(zhǔn)確率分別為0.037、0.176、0.918、95.85%,與其他5種模型相比,效果最優(yōu)。相比單純使用BP神經(jīng)網(wǎng)絡(luò),加入LASSO和SA-PSO算法2與總準(zhǔn)確率分別提高了37.43%、11.09個百分點,MAE與RMSE分別下降了79.10%、44.83%。因此,SA-PSO-LASSO-BP神經(jīng)網(wǎng)絡(luò)通過LASSO算法對9種環(huán)境參數(shù)進(jìn)行特征選擇,有效地提高了BP神經(jīng)網(wǎng)絡(luò)的評價性能;其次,采用SA-PSO算法優(yōu)化網(wǎng)絡(luò)的初始權(quán)值和閾值,更充分地挖掘了環(huán)境參數(shù)間的耦合關(guān)系,更好地擬合了妊娠豬舍環(huán)境參數(shù)與環(huán)境質(zhì)量間的非線性關(guān)系,大幅度增強(qiáng)網(wǎng)絡(luò)的評價效率,提高環(huán)境質(zhì)量評價精度。

    表5 不同神經(jīng)網(wǎng)絡(luò)模型評價結(jié)果

    4 討 論

    通過實際數(shù)據(jù)訓(xùn)練結(jié)果可知,與現(xiàn)有文獻(xiàn)中常用的模糊綜合評價法[4-6]相比,采用的BP神經(jīng)網(wǎng)絡(luò)充分挖掘了環(huán)境參數(shù)間的耦合關(guān)系,擬合復(fù)雜環(huán)境參數(shù)與環(huán)境質(zhì)量間的非線性關(guān)系,可以有效解決模糊綜合評價存在的欠學(xué)習(xí)問題?,F(xiàn)有文獻(xiàn)中的機(jī)器學(xué)習(xí)方法[7-10]進(jìn)行豬舍環(huán)境質(zhì)量評價時,未對環(huán)境參數(shù)選取原則與其對環(huán)境質(zhì)量的影響程度作深入分析,未能考慮不同豬舍內(nèi)環(huán)境參數(shù)受外環(huán)境、地理位置等因素影響所產(chǎn)生的變化范圍和趨勢的實際差異。運(yùn)用LASSO算法選擇特征環(huán)境參數(shù),解決機(jī)器學(xué)習(xí)存在的普適性和泛化性差等問題,同時實現(xiàn)了BP神經(jīng)網(wǎng)絡(luò)模型輸入的降維,有效降低了網(wǎng)絡(luò)結(jié)構(gòu)的復(fù)雜程度。由于BP神經(jīng)網(wǎng)絡(luò)本身容易受到初始權(quán)值和閾值等參數(shù)的影響,存在收斂速度慢、易陷入局部極小值等不足,使用SA-PSO算法對BP神經(jīng)網(wǎng)絡(luò)的初始權(quán)值和閾值進(jìn)行優(yōu)化,有效地提高了網(wǎng)絡(luò)的收斂速度與收斂精度。

    本研究提出的SA-PSO-LASSO-BP神經(jīng)網(wǎng)絡(luò)模型,通過LASSO和SA-PSO算法對BP神經(jīng)網(wǎng)絡(luò)的結(jié)構(gòu)與參數(shù)進(jìn)行優(yōu)化,達(dá)到了提高BP神經(jīng)網(wǎng)絡(luò)預(yù)測精度與性能的目的,實現(xiàn)妊娠豬舍環(huán)境質(zhì)量的可靠評價。采用的特征選擇和BP神經(jīng)網(wǎng)絡(luò)模型,是在現(xiàn)有實際數(shù)據(jù)基礎(chǔ)上構(gòu)建的具有較佳效果的模型,可以進(jìn)一步嘗試采用DEFS(Differential Evolution Feature Selection)算法[13]、LightGBM(Light gradient boosting machine)算法[26]進(jìn)行特征參數(shù)篩選,采用SVM(Support Vector Machine,SVM)[27]、ELM(Extreme Learning Machine,ELM)[28]等方法構(gòu)建評價模型。

    此外,計算流體力學(xué)(Computational Fluid Dynamics ,CFD)近年來被廣泛應(yīng)用于豬舍送風(fēng)降溫[29]、排風(fēng)換氣[30]系統(tǒng)的設(shè)計,以及環(huán)境調(diào)控[31]、通風(fēng)效果[32]的評價。因此,在后續(xù)的研究中可利用CFD對妊娠豬舍通風(fēng)降溫系統(tǒng)進(jìn)行優(yōu)化設(shè)計,結(jié)合本文所提出的豬舍環(huán)境質(zhì)量評價方法,分析不同通風(fēng)風(fēng)速、通風(fēng)溫度等條件下舍內(nèi)的通風(fēng)效果,制定通風(fēng)、溫控等設(shè)備的調(diào)控策略,以滿足豬舍適宜性控制的應(yīng)用要求。

    5 結(jié) 論

    為合理地評價妊娠豬舍環(huán)境質(zhì)量,本研究提出一種基于SA-PSO-LASSO-BP(Simulated Annealing-Particle Swarm Optimization-Least Absolute Shrinkage and Selection Operator-Back Propagation)的環(huán)境質(zhì)量評價模型。利用卡爾曼濾波和分批估計自適應(yīng)加權(quán)融合算法對環(huán)境數(shù)據(jù)進(jìn)行預(yù)處理,采用LASSO算法篩選特征環(huán)境參數(shù),使用SA-PSO算法優(yōu)化BP神經(jīng)網(wǎng)絡(luò)的初始權(quán)值和閾值,得到如下結(jié)論:

    1)利用卡爾曼濾波和分批估計自適應(yīng)加權(quán)融合算法對多采集節(jié)點的原始環(huán)境數(shù)據(jù)進(jìn)行預(yù)處理,單一傳感器時間序列數(shù)據(jù)融合結(jié)果的方差較濾波前下降了4.03%;多源同質(zhì)傳感器時間序列數(shù)據(jù)融合結(jié)果的方差較算術(shù)平均法至少下降了92.14%,能夠有效地抑制噪聲干擾、降低環(huán)境數(shù)據(jù)的冗余性。

    2)采用LASSO算法對預(yù)處理后的環(huán)境數(shù)據(jù)進(jìn)行特征選擇,較全環(huán)境參數(shù)模型,LASSO-BP神經(jīng)網(wǎng)絡(luò)模型迭代次數(shù)下降了29.90%,網(wǎng)絡(luò)誤差降低了17.78%。SA-PSO-LASSO-BP神經(jīng)網(wǎng)絡(luò)模型的迭代次數(shù)為36次,運(yùn)行時間為1.20 s,網(wǎng)絡(luò)誤差為0.033,相較于GA-LASSO- BP神經(jīng)網(wǎng)絡(luò)模型迭代次數(shù)減少了54.43%,運(yùn)行時間加快了17.24%,網(wǎng)絡(luò)誤差下降了8.33%,網(wǎng)絡(luò)收斂精度和收斂速度顯著提高。

    3)采用實際妊娠豬舍環(huán)境數(shù)據(jù)作為測試樣本對SA-PSO-LASSO-BP神經(jīng)網(wǎng)絡(luò)模型及其他5模型的評價效果進(jìn)行驗證,結(jié)果表明:SA-PSO-LASSO-BP神經(jīng)網(wǎng)絡(luò)模型的平均絕對誤差為0.037、均方根誤差為0.176、決定系數(shù)為0.918、總準(zhǔn)確率為95.85%,相比單純使用BP神經(jīng)網(wǎng)絡(luò),加入LASSO和SA-PSO算法之后決定系數(shù)與總準(zhǔn)確率分別提高了37.43%、11.09個百分點,平均絕對誤差與均方根誤差分別下降了79.10%、44.83%,具有更高的評價精度和性能。

    本文提出基于SA-PSO-LASSO-BP神經(jīng)網(wǎng)絡(luò)模型的豬舍環(huán)境質(zhì)量評價方法,可充分挖掘環(huán)境參數(shù)間的耦合關(guān)系,更好地擬合復(fù)雜環(huán)境因素與環(huán)境質(zhì)量間的非線性關(guān)系,為豬舍環(huán)境提供有效的質(zhì)量評價手段。

    [1] 李修松,葉章穎,李國銘,等. 規(guī)?;i場妊娠母豬舍改進(jìn)濕簾降溫系統(tǒng)的環(huán)境特性[J]. 農(nóng)業(yè)工程學(xué)報,2020,36(20):238-245.

    Li Xiusong, Ye Zhangying, Li Guoming, et al. Environmental characteristics of a swine gestation barn with an innovative cooling pad system[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2020, 36(20): 238-245. (in Chinese with English abstract)

    [2] 韓戰(zhàn)強(qiáng),銀嶺,馬亞姣,等. 飼養(yǎng)模式、新型飼料添加劑、豬舍溫?zé)岘h(huán)境對母豬繁殖性能的影響[J]. 養(yǎng)豬,2021(6):35-38.

    Han Zhanqiang, Yin Ling, Ma Yajiao, et al. Research progress on the effects of feeding pattern, feed additives and warm and hot environment on sow reproductive performance[J]. Swine Production, 2021(6): 35-38. (in Chinese with English abstract)

    [3] Li J, Han Q. Effects of environment and breed on growth performance and meat quality of fattening pigs[J]. Animal Welfare, 2020, 29(2): 177-184.

    [4] 謝秋菊,蘇中濱,劉超吉,等. 豬舍環(huán)境適宜性模糊綜合評價[J]. 農(nóng)業(yè)工程學(xué)報,2016,32(16):198-205.

    Xie Qiuju, Su Zhongbin, Liu Chaoji, et al. Fuzzy synthetic assessment of swine house environmental adaptability[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2016, 32(16): 198-205. (in Chinese with English abstract)

    [5] 沈翠鳳. 基于物聯(lián)網(wǎng)的畜禽舍環(huán)境遠(yuǎn)程監(jiān)控和評判系統(tǒng)研究[J]. 中國農(nóng)機(jī)化學(xué)報,2020,41(12):112-118.

    Shen Cuifeng. Research on remote monitoring and evaluation system of livestock and poultry house environment based on Internet of Things[J]. Journal of Chinese Agricultural Mechanization, 2020, 41(12): 112-118. (in Chinese with English abstract)

    [6] 梁帆,劉振濤. 豬舍環(huán)境控制評價指標(biāo)體系構(gòu)建與分析:以河北省規(guī)模豬場為例[J]. 豬業(yè)科學(xué),2021,38(8):86-91.

    Liang Fan, Liu Zhentao. Construction and analysis of evaluation index system of piggery environmental control: A case study of large-scale pig farms in Hebei Province[J]. Swine Industry Science, 2021, 38(8): 86-91. (in Chinese with English abstract)

    [7] 程捷,馮天玉,黃世明,等. 基于D-S證據(jù)理論的豬舍環(huán)境狀態(tài)識別研究[J]. 中國農(nóng)機(jī)化學(xué)報,2021,42(6):50-54.

    Cheng Jie, Feng Tianyu, Huang Shiming, et al. Research on recognition of piggery environmental state based on D-S evidence theory[J]. Journal of Chinese Agricultural Mechanization, 2021, 42(6): 50-54. (in Chinese with English abstract)

    [8] 孫聰. 基于物聯(lián)網(wǎng)豬舍多參數(shù)測控系統(tǒng)設(shè)計及環(huán)境適宜度評價研究[D]. 鎮(zhèn)江:江蘇大學(xué),2019.

    Sun Cong. Design and Environmental Suitability Evaluation of Multi-parameter Measurement and Control System in Pig House Based on Internet of Things[D]. Zhenjiang: Jiangsu University, 2019. (in Chinese with English abstract)

    [9] 朱保釧. 基于物聯(lián)網(wǎng)豬舍參數(shù)測控及環(huán)境適宜度預(yù)測模型研究[D]. 鎮(zhèn)江:江蘇大學(xué),2020.

    Zhu Baochuan. Research on the Prediction Model of Parameters Measurement and Control and Environmental Suitability of Piggery Based on Internet of Things[D]. Zhenjiang: Jiangsu University, 2020. (in Chinese with English abstract)

    [10] 陳沖,劉星橋,劉超吉,等. 哺乳母豬舍環(huán)境舒適度評價預(yù)測模型優(yōu)化[J]. 農(nóng)業(yè)機(jī)械學(xué)報,2020,51(8):311-319.

    Chen Chong, Liu Xingqiao, Liu Chaoji, et al. Optimization of evaluation and prediction model of environmental comfort in lactating sow house[J]. Transactions of the Chinese Society for Agricultural Machinery, 2020, 51(8): 311-319. (in Chinese with English abstract)

    [11] 劉春紅,楊亮,鄧河,等. 基于ARMA和BP神經(jīng)網(wǎng)絡(luò)的豬舍氨氣濃度預(yù)測[J]. 中國環(huán)境科學(xué),2019,39(6):2320-2327.

    Liu Chunhong, Yang Liang, DengHe, et al. Prediction of ammonia concentration in piggery based on ARIMA and BP neural network[J]. China Environmental Science, 2019, 39(6): 2320-2327. (in Chinese with English abstract)

    [12] 楊斷利,李今,陳輝,等. 基于IPSO優(yōu)化BP神經(jīng)網(wǎng)絡(luò)的蛋雞舍有害氣體監(jiān)測系統(tǒng)[J]. 農(nóng)業(yè)機(jī)械學(xué)報,2021,52(4):327-335.

    Yang Duanli, Li Jin, Chen Hui, et al. Monitoring system of harmful gas in layer house based on improved Particle Swarm Optimization BP neural network[J]. Transactions of the Chinese Society for Agricultural Machinery, 2021, 52(4): 327-335. (in Chinese with English abstract)

    [13] 趙建輝,張晨陽,閔林,等. 基于特征選擇和GA-BP神經(jīng)網(wǎng)絡(luò)的多源遙感農(nóng)田土壤水分反演[J]. 農(nóng)業(yè)工程學(xué)報,2021,37(11):112-120.

    Zhao Jianhui, Zhang Chenyang, Min Lin, et al. Retrieval for soil moisture in farmland using multi-source remote sensing data and feature selection with GA-BP neural network[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2021, 37(11): 112-120. (in Chinese with English abstract)

    [14] 中華人民共和國國家質(zhì)量監(jiān)督檢驗檢疫總局,中國國家標(biāo)準(zhǔn)化管理委員會. GB/T 17824.3—2008 《規(guī)模豬場環(huán)境參數(shù)及環(huán)境管理》[S]. 北京:中國標(biāo)準(zhǔn)出版社,2008.

    [15] 申仲健,吳中紅,張錦瑞,等. 新風(fēng)過濾和排風(fēng)除臭裝置條件下樓房豬舍夏季環(huán)境監(jiān)測[J]. 農(nóng)業(yè)工程學(xué)報,2022,38(11):215-222.

    Shen Zhongjian, Wu Zhonghong, Zhang Jinrui, et al. Environmental monitoring of a multi-floor pig building under the condition equipped with fresh air filters and exhaust deodorization devices in summer[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2022, 38(11): 215-222. (in Chinese with English abstract)

    [16] Chen C, Liu X Q. An intelligent monitoring system for a pig breeding environment based on a wireless sensor network[J]. International Journal of Sensor Networks, 2019, 29(4): 275-283.

    [17] 杜鵬,包曉安,胡逸飛,等. 基于卡爾曼濾波的無線傳感網(wǎng)時空數(shù)據(jù)融合算法[J]. 電子科技,2022,35(6):21-27.

    Du Peng, Bao Xiao'an, Hu Yifei, et al. Research on spatio-temporal data fusion algorithm of wireless sensor network based on Kalman filter[J]. Electronic Science and Technology, 2022, 35(6): 21-27. (in Chinese with English abstract)

    [18] Wang H D, Wang D Y. An improved multiple wireless sensor data batch estimation adaptive weighted fusion algorithm[J]. Chinese Journal of Sensors and Actuators, 2015, 28(8): 1239-1243.

    [19] Robert T. Regression Shrinkage and Selection Via the Lasso[J]. Journal of the Royal Statistical Society: Series B (Methodological), 2018, 58(1): 267-288.

    [20] 于瀟瀅,韓蕊,秦貴信. 環(huán)境溫度對豬生產(chǎn)性能、養(yǎng)分消化和產(chǎn)熱量的影響[J]. 黑龍江畜牧獸醫(yī),2017(17):136-138.

    Yu Xiaoying, Han Rui, Qin Guixin. The effect of environmental temperature on growth performance, nutrient digestibility and heat production in pigs[J]. Heilongjiang Animal Science and Veterinary Medicine, 2017(17): 136-138. (in Chinese with English abstract)

    [21] 田宗祥,郭志明,張玲清. 環(huán)境溫度對豬采食量的影響研究進(jìn)展[J]. 國外畜牧學(xué)(豬與禽),2017,37(9):38-39.

    Tian Zongxiang, Guo Zhiming, Zhang Lingqing. Research progress of the effect of ambient temperature on feed intake of pigs[J]. Animal Science Abroad(Pigs and Poultry), 2017, 37(9): 38-39. (in Chinese with English abstract)

    [22]王浩云,宋進(jìn),潘磊慶,等. 優(yōu)化BP神經(jīng)網(wǎng)絡(luò)提高高光譜檢測調(diào)理雞肉菌落總數(shù)精度[J]. 農(nóng)業(yè)工程學(xué)報,2020,36(5):302-309.

    Wang Haoyun, Song Jin, Pan Leiqing, et al. Improving hyperspectral detection accuracy of total bacteria in prepared chicken using optimized BP neural network[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2020, 36(5): 302-309. (in Chinese with English abstract)

    [23]陳嘯,王紅英,孔丹丹,等. 基于粒子群參數(shù)優(yōu)化和BP神經(jīng)網(wǎng)絡(luò)的顆粒飼料質(zhì)量預(yù)測模型[J]. 農(nóng)業(yè)工程學(xué)報,2016,32(14):306-314.

    Chen Xiao, Wang Hongying, Kong Dandan, et al. Quality prediction model of pellet feed basing on BP network using PSO parameters optimization method [J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2016, 32(14): 306-314. (in Chinese with English abstract)

    [24] 劉雪,李亞妹,劉嬌,等. 基于BP神經(jīng)網(wǎng)絡(luò)的鮮雞蛋貨架期預(yù)測模型[J]. 農(nóng)業(yè)機(jī)械學(xué)報,2015,46(10):328-334.

    Liu Xue, Li Yamei, Liu Jiao, et al. BP neural network based prediction model for fresh egg's shelf life[J]. Transactions of the Chinese Society for Agricultural Machinery, 2015, 46(10): 328-334. (in Chinese with English abstract)

    [25] 李平. 基于SA-PSO的Plug-In混合動力汽車模糊控制策略優(yōu)化研究[D]. 濟(jì)南:山東大學(xué),2017.

    Li Ping. Optimization Research on Fuzzy Control Strategy for a Plug-In Hybrid Vehicle Based on SA-PSO Algorithm[D]. Jinan: Shandong University, 2017. (in Chinese with English abstract)

    [26] 尹航,呂佳威,陳耀聰,等. 基于LightGBM-SSA-ELM的新疆羊舍CO2濃度預(yù)測[J]. 農(nóng)業(yè)機(jī)械學(xué)報,2022,53(1):261-270.

    Yin Hang, Lu Jiawei, Chen Yaocong, et al. Prediction of CO2concentration in Xinjiang breeding environment of mutton sheep based on LightGBM-SSA-ELM[J]. Transactions of the Chinese Society for Agricultural Machinery, 2022, 53(1): 261-270. (in Chinese with English abstract)

    [27] 張?zhí)炝?,張東興,崔濤,等. 基于支持向量機(jī)和ReliefF算法的玉米品種抗倒伏預(yù)測[J]. 農(nóng)業(yè)工程學(xué)報,2021,37(20):226-233.

    Zhang Tianliang, Zhang Dongxing, Cui Tao, et al. Lodging resistance prediction of maize varieties based on support vector machine and ReliefF algorithm[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2021, 37(20): 226-233. (in Chinese with English abstract)

    [28] 施珮,匡亮,袁永明,等. 基于改進(jìn)極限學(xué)習(xí)機(jī)的水體溶解氧預(yù)測方法[J]. 農(nóng)業(yè)工程學(xué)報,2020,36(19):225-232.

    Shi Pei, Kuang Liang, Yuan Yongming, et al. Dissolved oxygen prediction for water quality of aquaculture using improved ELM network[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2020, 36(19): 225-232. (in Chinese with English abstract)

    [29] 張校軍,王華,龍沈飛,等. 定點送風(fēng)降溫對母豬周圍溫度和對流換熱的影響[J]. 農(nóng)業(yè)工程學(xué)報,2022,38(5):225-232.

    Zhang Xiaojun, Wang Hua, Long Shenfei, et al. Effects of fixed-point air supply cooling system on airflow and temperature around sows and convective heat transfer of sows[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2022, 38(5): 225-232. (in Chinese with English abstract)

    [30] 吳志東,房俊龍,巴文革,等. 冬季采暖保育豬舍送排風(fēng)管道組合換氣系統(tǒng)設(shè)計與評價[J]. 農(nóng)業(yè)工程學(xué)報,2021,37(10):152-158.

    Wu Zhidong, Fang Junlong, Ba Wenge, et al. Design and evaluation of combined ventilation system with air supply and exhaust ducts for a nursery heating piggery in winter[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2021, 37(10): 152-158. (in Chinese with English abstract)

    [31] 劉營芳,黃子碩. 基于輻射供熱和個性化送風(fēng)的分娩豬舍環(huán)境保障系統(tǒng)[J]. 農(nóng)業(yè)工程學(xué)報,2022,38(2):260-268.

    Liu Yingfang, Huang Zishuo. Thermal environment regulation and energy consumption of delivery pigsty based on radiation and personalized ventilation[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2022, 38(2): 260-268. (in Chinese with English abstract)

    [32] Marco B, Enrica S, Alberto B, et al. Assessment of geometrical and seasonal effects on the natural ventilation of a pig barn using CFD simulations[J]. Computers and Electronics in Agriculture, 2022, 193: 106652.

    Environmental quality evaluation method for swine gestation barns based on multi-source information fusion

    Chi Yu1, Guo Yanjiao1, Feng Han1, Li Han2, Zheng Yongjun1,3※

    (1.,,100083,; 2.,,100083,;3.(),100083,)

    Environmental quality of swine gestation barns can bring a significant impact on the fertility of breeding sows. Therefore, it is crucial to accurately evaluate the environmental quality and then timely trim the conditions, particularly for high breeding efficiency under less environmental stress. In this study, an environmental quality evaluation model of swine gestation barns was proposed using the Simulated Annealing-Particle Swarm Optimization-Least Absolute Shrinkage and Selection Operator-Back Propagation (SA-PSO-LASSO-BP) Neural Network (NN). Firstly, nine parameters were identified using the Chinese National Criteria. A data collection system was then established to collect the environmental data. Secondly, a Kalman filter and a batch estimation adaptive weighted fusion algorithm were introduced to fuse the multi-node environmental data, in order to remove the errors and redundant data from the data collection. Thirdly, a Least Absolute Shrinkage and Selection Operator (LASSO) regression model was selected for the feature selection. There were four feature factors that were closely related to environmental quality, including temperature, relative humidity, NH3concentration, and CO2concentration. Meanwhile, the structural parameters were optimized in the BP-NN , where the number of hidden layer nodes was determined to be 11. Finally, the initial weights and threshold values of the BP NN were optimized by the SA-PSO for the ultimate evaluation model. A comparison was made on the several NNs to verify the evaluation performance of the SA-PSO-LASSO-BP NN, including the BP, LASSO-BP, Genetic Algorithm-LASSO-BP (GA-LASSO-BP), Sparrow Search Algorithm-LASSO-BP (SSA-LASSO-BP), and the PSO-LASSO-BP NN. The training results proved that the convergence accuracy and rate of the SA-PSO-LASSO-BP network were significantly improved by the feature selection with the LASSO regression model. The number of iterations and the network errors of the LASSO-BP NN decreased by 29.9% and 17.78%, respectively, compared with the BP NN. In terms of feature selection, the SA-PSO algorithm implemented by the SA-PSO-LASSO-BP NN was utilized to optimize the initial weights and thresholds of the network, in order to further improve the convergence accuracy and rate of the model. Compared with the PSO-LASSO-BP NN, the number of iterations, running time, and network error were reduced by 54.43%, 17.24%, and 8.33%, respectively. The validation test indicated that the best performance of the model was achieved, with the coefficient of determination (2) of 0.918, an overall accuracy of 95.85%, the Mean Absolute Error (MAE) of 0.037, and the Root Mean Squared Error (RMSE) of 0.176. Consequently, the SA-PSO-LASSO-BP NN model can better fit the nonlinear relationship between complex environmental factors and environmental quality. The finding can serve as a strong reference for the environmental quality evaluation of swine gestation barns.

    models; environment; swine gestation barns; environmental quality; BP neural network; Least Absolute Shrinkage and Selection Operator (LASSO) algorithm; Simulate Anneal-Particle Swarm Optimization (SA-PSO) algorithm

    10.11975/j.issn.1002-6819.2022.18.023

    TU264+. 3

    A

    1002-6819(2022)-18-0212-10

    遲宇,郭艷嬌,馮涵,等. 采用多源信息融合的妊娠豬舍環(huán)境質(zhì)量評價方法[J]. 農(nóng)業(yè)工程學(xué)報,2022,38(18):212-221.doi:10.11975/j.issn.1002-6819.2022.18.023 http://www.tcsae.org

    Chi Yu, Guo Yanjiao, Feng Han, et al. Environmental quality evaluation method for swine gestation barns based on multi-source information fusion[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2022, 38(18): 212-221. (in Chinese with English abstract) doi:10.11975/j.issn.1002-6819.2022.18.023 http://www.tcsae.org

    2022-08-05

    2022-09-07

    國家重點研發(fā)計劃項目(2016YFD0700204)

    遲宇,研究方向為圖像與信息處理技術(shù)。Email:selivia0328@163.com

    鄭永軍,博士,教授,博士生導(dǎo)師,研究方向為農(nóng)業(yè)智能技術(shù)與裝備。Email:zyj@cau.edu.cn

    猜你喜歡
    融合環(huán)境評價
    村企黨建聯(lián)建融合共贏
    SBR改性瀝青的穩(wěn)定性評價
    石油瀝青(2021年4期)2021-10-14 08:50:44
    長期鍛煉創(chuàng)造體內(nèi)抑癌環(huán)境
    融合菜
    一種用于自主學(xué)習(xí)的虛擬仿真環(huán)境
    從創(chuàng)新出發(fā),與高考數(shù)列相遇、融合
    《融合》
    孕期遠(yuǎn)離容易致畸的環(huán)境
    環(huán)境
    基于Moodle的學(xué)習(xí)評價
    天天一区二区日本电影三级| 亚洲精品久久成人aⅴ小说| 国产97色在线日韩免费| 国产精品av久久久久免费| www日本黄色视频网| 亚洲熟女毛片儿| 国产av麻豆久久久久久久| 精品久久久久久久人妻蜜臀av| 岛国在线观看网站| 亚洲av片天天在线观看| 免费观看精品视频网站| 日韩高清综合在线| 51午夜福利影视在线观看| 色综合站精品国产| 欧美黄色片欧美黄色片| 精品久久久久久,| 999精品在线视频| 9191精品国产免费久久| 国产真实乱freesex| 欧美在线一区亚洲| 长腿黑丝高跟| 国产真人三级小视频在线观看| 成人特级黄色片久久久久久久| 一本大道久久a久久精品| 高清在线国产一区| 国内精品久久久久久久电影| 精品久久久久久久末码| 国产午夜精品久久久久久| 亚洲成人精品中文字幕电影| 亚洲精品在线美女| 国产免费男女视频| 免费在线观看视频国产中文字幕亚洲| 久久久久久久久免费视频了| 国产视频内射| 中文字幕精品亚洲无线码一区| 一级片免费观看大全| 69av精品久久久久久| 久久中文字幕一级| АⅤ资源中文在线天堂| 十八禁网站免费在线| 99re在线观看精品视频| 一区二区三区高清视频在线| 超碰成人久久| 精品一区二区三区av网在线观看| 国产69精品久久久久777片 | 免费在线观看完整版高清| 日本黄色视频三级网站网址| 激情在线观看视频在线高清| 看黄色毛片网站| 中文在线观看免费www的网站 | 熟女电影av网| 亚洲自偷自拍图片 自拍| 后天国语完整版免费观看| 淫秽高清视频在线观看| 日韩欧美三级三区| 在线免费观看的www视频| 99热只有精品国产| 一二三四社区在线视频社区8| 久久久久久久久久黄片| 在线观看一区二区三区| 无遮挡黄片免费观看| 久久亚洲精品不卡| 午夜免费成人在线视频| 亚洲av五月六月丁香网| 两性午夜刺激爽爽歪歪视频在线观看 | 可以在线观看的亚洲视频| 国产亚洲av嫩草精品影院| 国产真实乱freesex| 国产日本99.免费观看| 人妻久久中文字幕网| 国产又黄又爽又无遮挡在线| 精品久久久久久久人妻蜜臀av| 亚洲专区字幕在线| 曰老女人黄片| 日本成人三级电影网站| 日本黄大片高清| 久久精品91蜜桃| av免费在线观看网站| 精品少妇一区二区三区视频日本电影| 超碰成人久久| 久久精品亚洲精品国产色婷小说| 香蕉久久夜色| 亚洲午夜理论影院| 色综合婷婷激情| 亚洲精品av麻豆狂野| 国产1区2区3区精品| 国产亚洲精品一区二区www| 欧美精品亚洲一区二区| 999久久久精品免费观看国产| ponron亚洲| 欧美av亚洲av综合av国产av| 国产精品久久久久久精品电影| 国产精品香港三级国产av潘金莲| 久久人人精品亚洲av| 久久久久亚洲av毛片大全| 2021天堂中文幕一二区在线观| 在线观看www视频免费| 麻豆成人av在线观看| 久久精品人妻少妇| 欧美av亚洲av综合av国产av| 在线永久观看黄色视频| 午夜影院日韩av| 日韩成人在线观看一区二区三区| 日韩成人在线观看一区二区三区| 搡老岳熟女国产| 久久亚洲精品不卡| 深夜精品福利| 舔av片在线| 精品国产乱码久久久久久男人| 最近最新中文字幕大全免费视频| 欧美zozozo另类| 欧美成人一区二区免费高清观看 | 欧美黑人欧美精品刺激| 国内揄拍国产精品人妻在线| 国产在线观看jvid| 免费电影在线观看免费观看| 日韩欧美三级三区| 久久精品aⅴ一区二区三区四区| 亚洲片人在线观看| 18禁黄网站禁片免费观看直播| 国产又色又爽无遮挡免费看| 18禁黄网站禁片免费观看直播| 亚洲成人中文字幕在线播放| 日韩精品青青久久久久久| 亚洲成人久久性| 亚洲自拍偷在线| 少妇粗大呻吟视频| 免费在线观看影片大全网站| 黄色 视频免费看| 999久久久国产精品视频| 999精品在线视频| 黄色 视频免费看| 12—13女人毛片做爰片一| 国产一级毛片七仙女欲春2| 日韩 欧美 亚洲 中文字幕| 午夜福利欧美成人| 国产不卡一卡二| 国产三级在线视频| 日韩精品青青久久久久久| 成人国产一区最新在线观看| 亚洲专区中文字幕在线| 少妇粗大呻吟视频| 12—13女人毛片做爰片一| 成人av一区二区三区在线看| 人妻夜夜爽99麻豆av| 久久婷婷成人综合色麻豆| 又粗又爽又猛毛片免费看| av天堂在线播放| 丁香欧美五月| 一区福利在线观看| 亚洲国产欧洲综合997久久,| 午夜福利在线观看吧| 国产真实乱freesex| 欧美一级毛片孕妇| 精品午夜福利视频在线观看一区| 欧美国产日韩亚洲一区| 欧美高清成人免费视频www| 国产精品99久久99久久久不卡| 搡老熟女国产l中国老女人| 久久精品亚洲精品国产色婷小说| 国产精品一区二区三区四区免费观看 | 一级片免费观看大全| 亚洲五月天丁香| 色播亚洲综合网| 好男人在线观看高清免费视频| 成人av在线播放网站| 久久精品国产99精品国产亚洲性色| 很黄的视频免费| 亚洲第一欧美日韩一区二区三区| av福利片在线| 久久草成人影院| 精品国产乱子伦一区二区三区| 中文字幕高清在线视频| 亚洲va日本ⅴa欧美va伊人久久| 777久久人妻少妇嫩草av网站| 亚洲 欧美 日韩 在线 免费| 欧美三级亚洲精品| 国产亚洲精品第一综合不卡| 欧美黄色片欧美黄色片| 久久人妻福利社区极品人妻图片| 最近最新中文字幕大全免费视频| 午夜福利在线观看吧| 国产成人aa在线观看| 亚洲人成伊人成综合网2020| 国产片内射在线| 中文字幕高清在线视频| 色老头精品视频在线观看| 欧美不卡视频在线免费观看 | 丝袜美腿诱惑在线| 中国美女看黄片| 欧美色视频一区免费| 亚洲人与动物交配视频| 特大巨黑吊av在线直播| 国内久久婷婷六月综合欲色啪| 日日摸夜夜添夜夜添小说| 超碰成人久久| 我的老师免费观看完整版| 在线永久观看黄色视频| 制服丝袜大香蕉在线| 精品久久久久久久末码| 免费在线观看成人毛片| 欧美大码av| 成人高潮视频无遮挡免费网站| 欧美成人午夜精品| 婷婷精品国产亚洲av在线| 日本 av在线| 亚洲国产欧美人成| 免费搜索国产男女视频| 最近在线观看免费完整版| 国产高清videossex| 在线国产一区二区在线| 亚洲av五月六月丁香网| 又粗又爽又猛毛片免费看| 国产蜜桃级精品一区二区三区| 天堂√8在线中文| 一夜夜www| 一本一本综合久久| 久久久国产成人精品二区| 精品一区二区三区av网在线观看| 国产午夜精品论理片| 国产精品亚洲一级av第二区| 亚洲精品一区av在线观看| 国产免费av片在线观看野外av| 制服丝袜大香蕉在线| 久久久精品国产亚洲av高清涩受| 最新在线观看一区二区三区| 一级毛片精品| 手机成人av网站| 日韩免费av在线播放| 欧美黑人欧美精品刺激| 午夜福利在线观看吧| 欧美日韩亚洲国产一区二区在线观看| 看黄色毛片网站| 一区二区三区高清视频在线| 亚洲精品粉嫩美女一区| 成年人黄色毛片网站| 日本在线视频免费播放| netflix在线观看网站| 久久香蕉精品热| 97碰自拍视频| 色综合亚洲欧美另类图片| 亚洲无线在线观看| 男人的好看免费观看在线视频 | 每晚都被弄得嗷嗷叫到高潮| 午夜福利欧美成人| 久久人妻福利社区极品人妻图片| 桃红色精品国产亚洲av| 老司机福利观看| 欧美日本视频| 999精品在线视频| 欧美黑人精品巨大| 国产伦人伦偷精品视频| e午夜精品久久久久久久| 国产激情欧美一区二区| 日本黄大片高清| 丰满的人妻完整版| 亚洲av熟女| 色精品久久人妻99蜜桃| 久久性视频一级片| 国产精品亚洲av一区麻豆| 国产又黄又爽又无遮挡在线| 亚洲va日本ⅴa欧美va伊人久久| 国内少妇人妻偷人精品xxx网站 | 亚洲乱码一区二区免费版| 欧美黄色片欧美黄色片| 精品免费久久久久久久清纯| av福利片在线| 久久午夜综合久久蜜桃| 免费在线观看视频国产中文字幕亚洲| 亚洲电影在线观看av| 国产一区二区在线观看日韩 | e午夜精品久久久久久久| 国产av麻豆久久久久久久| 一区二区三区激情视频| 99久久精品热视频| 免费人成视频x8x8入口观看| 婷婷丁香在线五月| 欧美性长视频在线观看| 色在线成人网| 黄频高清免费视频| 欧美成人免费av一区二区三区| 9191精品国产免费久久| 欧美成人午夜精品| 国内揄拍国产精品人妻在线| 久久天躁狠狠躁夜夜2o2o| 国产激情久久老熟女| 亚洲午夜精品一区,二区,三区| 男人舔奶头视频| 淫秽高清视频在线观看| 国产成人av教育| 亚洲国产欧美网| av中文乱码字幕在线| 听说在线观看完整版免费高清| av福利片在线| 18禁裸乳无遮挡免费网站照片| 久久久国产成人免费| 无限看片的www在线观看| 亚洲成人精品中文字幕电影| 日韩av在线大香蕉| 久久久久精品国产欧美久久久| 成人高潮视频无遮挡免费网站| 国产三级在线视频| 国产高清视频在线观看网站| 日本一二三区视频观看| 成人18禁高潮啪啪吃奶动态图| 国产视频一区二区在线看| 1024香蕉在线观看| 日本 欧美在线| 欧洲精品卡2卡3卡4卡5卡区| www日本在线高清视频| 人人妻,人人澡人人爽秒播| 久久香蕉国产精品| 国产麻豆成人av免费视频| 国产片内射在线| 亚洲精品一卡2卡三卡4卡5卡| 国产真人三级小视频在线观看| 久久香蕉国产精品| 午夜精品久久久久久毛片777| 丁香六月欧美| 波多野结衣高清无吗| 丝袜美腿诱惑在线| 夜夜夜夜夜久久久久| 欧美性长视频在线观看| 美女免费视频网站| 亚洲一码二码三码区别大吗| 法律面前人人平等表现在哪些方面| 久热爱精品视频在线9| 99久久精品热视频| 日韩av在线大香蕉| 国产视频内射| 欧美乱码精品一区二区三区| 老司机靠b影院| 国产精品,欧美在线| 在线观看免费日韩欧美大片| 观看免费一级毛片| 一区二区三区激情视频| 免费av毛片视频| 欧美黑人巨大hd| √禁漫天堂资源中文www| 亚洲色图av天堂| 美女免费视频网站| av在线天堂中文字幕| 亚洲人与动物交配视频| 男男h啪啪无遮挡| 久久精品人妻少妇| 色综合站精品国产| 亚洲欧美激情综合另类| 中文资源天堂在线| 女人高潮潮喷娇喘18禁视频| 视频区欧美日本亚洲| 免费在线观看亚洲国产| 日本熟妇午夜| 国产蜜桃级精品一区二区三区| 久久性视频一级片| 欧美久久黑人一区二区| 中文亚洲av片在线观看爽| 国产精品,欧美在线| 99在线视频只有这里精品首页| 精品国产乱子伦一区二区三区| 亚洲中文字幕一区二区三区有码在线看 | 午夜激情福利司机影院| 亚洲av中文字字幕乱码综合| 色哟哟哟哟哟哟| 亚洲熟妇中文字幕五十中出| 可以在线观看的亚洲视频| 国产亚洲精品一区二区www| 亚洲自拍偷在线| 一个人观看的视频www高清免费观看 | 人妻夜夜爽99麻豆av| 国产亚洲精品久久久久5区| √禁漫天堂资源中文www| 欧美精品亚洲一区二区| 五月伊人婷婷丁香| 国产伦人伦偷精品视频| 亚洲成a人片在线一区二区| 亚洲色图av天堂| 国产精品久久久久久精品电影| 一本久久中文字幕| 免费一级毛片在线播放高清视频| 高清在线国产一区| 动漫黄色视频在线观看| 免费高清视频大片| 在线看三级毛片| 久久久久国产精品人妻aⅴ院| 99精品欧美一区二区三区四区| 在线观看日韩欧美| 精品日产1卡2卡| 国产精品国产高清国产av| 免费av毛片视频| 欧美成狂野欧美在线观看| 久久 成人 亚洲| 亚洲欧美日韩无卡精品| 男女视频在线观看网站免费 | 午夜福利高清视频| 舔av片在线| 欧美日韩亚洲国产一区二区在线观看| 成人特级黄色片久久久久久久| 床上黄色一级片| 老司机午夜十八禁免费视频| 白带黄色成豆腐渣| 久久中文看片网| 国产成人一区二区三区免费视频网站| 女人爽到高潮嗷嗷叫在线视频| 色综合亚洲欧美另类图片| 后天国语完整版免费观看| 老司机深夜福利视频在线观看| 香蕉av资源在线| 国产欧美日韩一区二区精品| 中文字幕av在线有码专区| 午夜福利18| 午夜免费激情av| 在线观看舔阴道视频| 国产一区二区在线观看日韩 | 美女扒开内裤让男人捅视频| 亚洲中文字幕一区二区三区有码在线看 | 日韩精品中文字幕看吧| 欧美午夜高清在线| 午夜老司机福利片| 欧美精品亚洲一区二区| 九九热线精品视视频播放| 日本撒尿小便嘘嘘汇集6| 日韩欧美在线乱码| 波多野结衣高清无吗| 亚洲人成77777在线视频| 99久久综合精品五月天人人| 久久久久久免费高清国产稀缺| av有码第一页| 成人三级做爰电影| 国产精品久久久久久久电影 | 国产激情欧美一区二区| 国内少妇人妻偷人精品xxx网站 | 俄罗斯特黄特色一大片| 国产片内射在线| 人人妻,人人澡人人爽秒播| 丝袜人妻中文字幕| netflix在线观看网站| 黄色女人牲交| 久久久精品国产亚洲av高清涩受| 日韩精品中文字幕看吧| 国产91精品成人一区二区三区| 国产精品美女特级片免费视频播放器 | 国模一区二区三区四区视频 | 老熟妇乱子伦视频在线观看| 一本大道久久a久久精品| 国内揄拍国产精品人妻在线| 欧美黑人巨大hd| 国产99久久九九免费精品| 欧美日本亚洲视频在线播放| 精品久久久久久久末码| 人人妻人人看人人澡| 国产蜜桃级精品一区二区三区| 亚洲人与动物交配视频| 午夜激情av网站| 欧美日韩中文字幕国产精品一区二区三区| 99久久综合精品五月天人人| 亚洲成人久久性| 大型av网站在线播放| 窝窝影院91人妻| 久久亚洲精品不卡| 国产成人av教育| 黑人巨大精品欧美一区二区mp4| 九九热线精品视视频播放| 国产av一区在线观看免费| 午夜免费成人在线视频| www.www免费av| 91九色精品人成在线观看| 国产精品,欧美在线| 少妇的丰满在线观看| 身体一侧抽搐| 欧美在线一区亚洲| 人人妻人人看人人澡| 宅男免费午夜| 久久草成人影院| 精品电影一区二区在线| 精品久久久久久久末码| 天天一区二区日本电影三级| 最近视频中文字幕2019在线8| 99热这里只有精品一区 | 国产激情久久老熟女| 97碰自拍视频| 十八禁人妻一区二区| 在线国产一区二区在线| 岛国视频午夜一区免费看| 成人三级做爰电影| 成人国产一区最新在线观看| 日本一二三区视频观看| 成人高潮视频无遮挡免费网站| 国产一区二区激情短视频| 国产av一区在线观看免费| 亚洲乱码一区二区免费版| 亚洲成av人片免费观看| 欧美激情久久久久久爽电影| 久久久久久大精品| 久久热在线av| 亚洲成人中文字幕在线播放| 三级毛片av免费| 动漫黄色视频在线观看| cao死你这个sao货| 欧美乱码精品一区二区三区| 欧美黄色片欧美黄色片| 亚洲精品国产一区二区精华液| 欧美日韩中文字幕国产精品一区二区三区| 国产精品 欧美亚洲| 亚洲成人精品中文字幕电影| av福利片在线| 全区人妻精品视频| 欧美成人性av电影在线观看| 一级毛片高清免费大全| 97超视频在线观看视频| 蜜桃久久精品国产亚洲av| 欧美日本亚洲视频在线播放| 精品一区二区三区视频在线| 在线观看66精品国产| 久久6这里有精品| 尾随美女入室| 一级毛片电影观看 | 99riav亚洲国产免费| 黄色配什么色好看| 深爱激情五月婷婷| 午夜视频国产福利| 久久精品综合一区二区三区| 中文字幕久久专区| 又粗又硬又长又爽又黄的视频 | 国产黄片美女视频| 日本爱情动作片www.在线观看| 亚洲人成网站在线观看播放| 麻豆国产97在线/欧美| 国内精品久久久久精免费| 亚洲一级一片aⅴ在线观看| 日本在线视频免费播放| 91aial.com中文字幕在线观看| 在线免费观看不下载黄p国产| 12—13女人毛片做爰片一| 禁无遮挡网站| 亚洲人成网站在线播放欧美日韩| 天天躁日日操中文字幕| 亚洲色图av天堂| 国产成人a∨麻豆精品| 国产免费男女视频| 成人特级av手机在线观看| 国产精品野战在线观看| 麻豆国产97在线/欧美| 国产日韩欧美在线精品| 99热网站在线观看| 久久久久久久久久久丰满| 亚洲美女视频黄频| 亚洲国产精品久久男人天堂| 3wmmmm亚洲av在线观看| 蜜桃亚洲精品一区二区三区| 精品人妻视频免费看| 一级av片app| 国产一区二区三区av在线 | 男的添女的下面高潮视频| 好男人在线观看高清免费视频| 日本-黄色视频高清免费观看| 一本久久精品| 夜夜看夜夜爽夜夜摸| 日韩高清综合在线| 欧美区成人在线视频| 欧美xxxx黑人xx丫x性爽| 国语自产精品视频在线第100页| 久久韩国三级中文字幕| 国产白丝娇喘喷水9色精品| 欧美在线一区亚洲| 少妇高潮的动态图| 中国美女看黄片| 三级男女做爰猛烈吃奶摸视频| 午夜福利在线观看吧| 国产精品1区2区在线观看.| 久久午夜福利片| 一级黄色大片毛片| 最后的刺客免费高清国语| av.在线天堂| 国产老妇伦熟女老妇高清| 1000部很黄的大片| 亚洲成人久久爱视频| 久久精品人妻少妇| 国产美女午夜福利| 波多野结衣高清无吗| 成人鲁丝片一二三区免费| 亚洲乱码一区二区免费版| 性欧美人与动物交配| av国产免费在线观看| 欧美成人免费av一区二区三区| 精品国产三级普通话版| 成年av动漫网址| 99久久精品热视频| 国产在视频线在精品| 国产成年人精品一区二区| 成人欧美大片| 亚洲最大成人av| 两性午夜刺激爽爽歪歪视频在线观看| 亚洲精品影视一区二区三区av| 三级经典国产精品| 欧美xxxx黑人xx丫x性爽| 全区人妻精品视频| 给我免费播放毛片高清在线观看| 成年女人看的毛片在线观看| 欧美色视频一区免费| 国产精品综合久久久久久久免费| 最近的中文字幕免费完整| 日本与韩国留学比较| 成熟少妇高潮喷水视频| 菩萨蛮人人尽说江南好唐韦庄 | 国产精华一区二区三区| 亚洲成人久久性| 国产日本99.免费观看| 精品久久久久久久久亚洲| 又爽又黄a免费视频| АⅤ资源中文在线天堂| 尾随美女入室| 99久久无色码亚洲精品果冻| 久久精品91蜜桃| 国产成人aa在线观看| 亚洲欧美精品专区久久|