999精品在线视频,手机成人午夜在线视频,久久不卡国产精品无码,中日无码在线观看,成人av手机在线观看,日韩精品亚洲一区中文字幕,亚洲av无码人妻,四虎国产在线观看 ?

粒子群優化的細菌覓食算法在變風量空調機組中的應用

2015-01-08 02:13:53孟顯嬌朱棟華殷文雪李思
中國科技縱橫 2014年12期

孟顯嬌++朱棟華++殷文雪++李思

【摘 要】 本文主要研究對象是變風量空調機組中的送風管道靜壓控制回路,將具有良好全局尋優能力的細菌覓食算法應用到該控制回路中,同時為了改善細菌覓食算法收斂速度較慢的缺點,將粒子群算法引入到細菌覓食算法中,對細菌覓食算法中的細菌位置更新進行優化,并利用優化后的細菌覓食算法對PID控制器的三個參數進行整定,將整定后的參數應用到控制回路中。通過matlab仿真,并將細菌覓食算法、粒子群算法,粒子群優化的細菌覓食算法進行比較,結果表明經過粒子群優化的細菌覓食算法收斂速度明顯加快。

【關鍵詞】 細菌覓食算法 粒子群算法 PID 送風管道靜壓

變風量空調憑借優越的性能在近年來得到了廣泛的關注,空調機組作為變風量空調系統中的重要設備,其運行直接影響到變風量空調的品質和節能效果。本文針對變風量空調機組送風管道靜壓控制回路進行研究,控制方法主要有定靜壓控制和變靜壓控制,其中定靜壓控制方法運行穩定,控制簡單,易于實現,是應用比較廣泛的方法。這種方法主要是通過調節風機轉速來實現風道靜壓恒定。PID控制器廣泛應用于控制領域,但是由于控制對象越來越復雜,人們對PID的控制要求也在不斷升高,本文利用粒子群優化的細菌覓食算法來整定PID參數,將優化后的參數應用到變風量空調機組送風管道靜壓控制回路中,取得滿意效果。

1 變風量空調機組概述

變風量空氣處理機組是變風量空調系統的關鍵設備,其主要功能是為空調系統提供足量的新風,維持送風溫度和送風管道靜壓在其設定值等。本文將針對送風管道靜壓控制回路進行研究,期望的目標是維持送風管道靜壓在設定值,這一過程是通過調節空氣處理機組中的風機轉速來實現,當靜壓實際值小于設定值,通過控制器調節變頻器,使風機轉速加大,反之,則通過控制器調節變頻器使風機轉速減小。送風管道靜壓控制回路框圖如圖1所示。

2 細菌覓食算法

細菌覓食算法由Kevin M. Passino于2002年提出,該算法模仿大腸桿菌在人體腸道內覓食行為,屬于仿生類優化算法,基本原理是將待優化的問題進行編碼,并定義待優化問題的解對應于搜索空間中的細菌狀態。針對具體問題的求解過程為:產生初始解群體、計算評價函數的值、利用群體的相互影響和作用機制進行優化[1]。細菌覓食算法主要通過趨向性操作、復制操作和驅散操作這三種方式的迭代計算來求解問題。

趨向性操作:趨向性操作是大腸桿菌向著食物豐富的地區前進的過程,這個過程的實現是通過大腸桿菌鞭毛的游動,當鞭毛逆時針時,大腸桿菌就會朝一個方向向前游動,這個過程稱為前進;當鞭毛順時針時,它就會減速直至停止然后重新選擇新的前進方向,這個過程稱為翻轉。大腸桿菌正是通過不斷的前進和翻轉來實現位置的更新,最終找到最優位置。細菌的趨向性操作可表示為如下形式:

(1)

(2)

其中,表示個體i的位置,其中j表示第j代趨向性循環,k表示第k代復制循環,l表示第l代驅散循環,表示選定的前進步長,為生成的隨即向量,表示細菌選定的前進方向。

復制操作:生物進化的規律是優勝劣汰,細菌在進行一定次數的趨向性操作之后會進行健康值排序,這里的健康值就是細菌在之前經歷的不同位置的適度值的累積值,根據排序結果,淘汰掉尋優效果不好的一半數量的細菌,再將剩下的細菌進行分裂,子細菌將繼承母細菌的步長及方向,通過這種方式既保證了新產生的細菌具有良好的尋優能力,同時也保證了菌群規模不變。

驅散操作:細菌的趨向性操作過程有時會導致算法陷入局部最優,為了解決這一問題在細菌覓食優化算法中引入了驅散操作。驅散操作就是將細菌以一定的概率驅散到搜索空間任意一個新的位置,將細菌在新位置的適度值與原位置適度值比較,若優于原位置,則說明經過驅散操作后細菌的全局尋優能力加強。

3 利用粒子群算法優化細菌覓食算法

細菌覓食算法可以在搜索最優解的過程中任意改變搜索方向,這樣大大的提高了算法局部搜索能力及搜索精度,但是細菌在趨向性操作中隨機翻轉,缺少細菌之間的相互學習,因此收斂速度較慢,針對這一不足,本文將粒子群中的粒子更新引入到細菌覓食算法中的細菌位置更新,由于在粒子群算法中粒子是根據本身的最佳位置和整個群體最佳位置來進行位置更新,這樣就可以利用粒子群算法中粒子群的記憶功能提高細菌覓食算法的搜索效率。在粒子群算法中,粒子群更新方式如下所示:

(3)

(4)

其中表示粒子本身最優解,表示整個粒子群的最優解,、表示加速因子,、為[0,1]之間的隨機數。

4 粒子群優化細菌覓食算法的具體步驟

(1)初始化細菌覓食算法以及粒子群算法相關參數,計算細菌的初始適度值。(2)趨向性操作:細菌按照公式(1)前進,如果前進后細菌所在位置的適度值優于前一個位置,則保持在該方向上繼續前進,直至達到設定游動次數,否則,細菌將通過翻轉改變原來前進方向,在每次的趨向性操作中記錄細菌本身經歷過的位置最優值以及全局位置最優值,并根據公式(3)來確定細菌的今后前進方向。(3)復制操作:將細菌在趨向性環節得到的適度值進行累加,累加結果記為,其中,對進行排序,淘汰掉一半細菌,將剩下的細菌進行分裂。(4)驅散操作:將符合驅散操作條件的細菌進行驅散。(5)完成以上流程,將最終的最優適度值及其所對應的細菌所在位置作為結果輸出。

5 系統的方案設計

本文利用粒子群優化的細菌覓食算法來整定PID參數,將優化后的參數應用到變風量空調機組送風管道靜壓控制回路中,具體的控制框圖如圖2所示。

6 系統仿真

粒子群優化的細菌覓食算法各個參數設置如下:細菌數量S設為20,參數維數P為3維,驅散次數為2次,繁殖次數為4次,趨化次數為40次,驅散概率為0.25,游動步數為4次,粒子群算法中的取2,取2。本次仿真選用的尋優函數為絕對偏差積分:

(5)

變風量空調機組的送風管道靜壓控制回路采用一階傳遞函數模型,具體表達式如下:

(6)

本次仿真的采樣周期為1s,采樣點數為100個,根據采樣時間以及傳遞函數表達式可以知道系統傳輸過程中有延遲,延遲約五個采樣周期,用Matlab編寫控制方法的m程序,得到相關PID參數進行仿真,其仿真結果如圖3所示。

粒子群算法,細菌覓食算法,粒子群優化的細菌覓食算法對應的PID值以及在本次仿真過程中的調節時間如表1所示。

根據以上圖表,在調節過程中粒子群算法出現超調量,細菌覓食算法及粒子群優化的細菌覓食算法沒有出現超調量。從調節時間角度看,粒子群算法和細菌覓食優化算法調節時間基本一致,經過粒子群優化的細菌覓食算法的調節時間較粒子群算法及細菌覓食算法有明顯縮短。

7 結語

通過以上分析,粒子群優化的細菌覓食算法與細菌覓食算法以及粒子群算法比較,經過優化后的細菌覓食算法在收斂速度上有明顯提升,調節時間縮短。

參考文獻:

[1]黃偉峰,林衛星,范懷科等.細菌覓食優化的智能PID控制[J].計算機工程與應用,2011,47(21):82-85.

[2]K. M. Passino.Biomimicry of Bacterial Foraging for Distributed Optimization and Control.IEEE.Control System Magazine,June,2002,pp:52-67.

[3]楊世忠,任慶昌.變風量空調靜壓的魯棒PID控制[J].控制工程,2013,20(6):1176-1180.

[4]劉小龍.細菌覓食優化算法的改進及應用[D].廣東:華南理工大學,2011.endprint

【摘 要】 本文主要研究對象是變風量空調機組中的送風管道靜壓控制回路,將具有良好全局尋優能力的細菌覓食算法應用到該控制回路中,同時為了改善細菌覓食算法收斂速度較慢的缺點,將粒子群算法引入到細菌覓食算法中,對細菌覓食算法中的細菌位置更新進行優化,并利用優化后的細菌覓食算法對PID控制器的三個參數進行整定,將整定后的參數應用到控制回路中。通過matlab仿真,并將細菌覓食算法、粒子群算法,粒子群優化的細菌覓食算法進行比較,結果表明經過粒子群優化的細菌覓食算法收斂速度明顯加快。

【關鍵詞】 細菌覓食算法 粒子群算法 PID 送風管道靜壓

變風量空調憑借優越的性能在近年來得到了廣泛的關注,空調機組作為變風量空調系統中的重要設備,其運行直接影響到變風量空調的品質和節能效果。本文針對變風量空調機組送風管道靜壓控制回路進行研究,控制方法主要有定靜壓控制和變靜壓控制,其中定靜壓控制方法運行穩定,控制簡單,易于實現,是應用比較廣泛的方法。這種方法主要是通過調節風機轉速來實現風道靜壓恒定。PID控制器廣泛應用于控制領域,但是由于控制對象越來越復雜,人們對PID的控制要求也在不斷升高,本文利用粒子群優化的細菌覓食算法來整定PID參數,將優化后的參數應用到變風量空調機組送風管道靜壓控制回路中,取得滿意效果。

1 變風量空調機組概述

變風量空氣處理機組是變風量空調系統的關鍵設備,其主要功能是為空調系統提供足量的新風,維持送風溫度和送風管道靜壓在其設定值等。本文將針對送風管道靜壓控制回路進行研究,期望的目標是維持送風管道靜壓在設定值,這一過程是通過調節空氣處理機組中的風機轉速來實現,當靜壓實際值小于設定值,通過控制器調節變頻器,使風機轉速加大,反之,則通過控制器調節變頻器使風機轉速減小。送風管道靜壓控制回路框圖如圖1所示。

2 細菌覓食算法

細菌覓食算法由Kevin M. Passino于2002年提出,該算法模仿大腸桿菌在人體腸道內覓食行為,屬于仿生類優化算法,基本原理是將待優化的問題進行編碼,并定義待優化問題的解對應于搜索空間中的細菌狀態。針對具體問題的求解過程為:產生初始解群體、計算評價函數的值、利用群體的相互影響和作用機制進行優化[1]。細菌覓食算法主要通過趨向性操作、復制操作和驅散操作這三種方式的迭代計算來求解問題。

趨向性操作:趨向性操作是大腸桿菌向著食物豐富的地區前進的過程,這個過程的實現是通過大腸桿菌鞭毛的游動,當鞭毛逆時針時,大腸桿菌就會朝一個方向向前游動,這個過程稱為前進;當鞭毛順時針時,它就會減速直至停止然后重新選擇新的前進方向,這個過程稱為翻轉。大腸桿菌正是通過不斷的前進和翻轉來實現位置的更新,最終找到最優位置。細菌的趨向性操作可表示為如下形式:

(1)

(2)

其中,表示個體i的位置,其中j表示第j代趨向性循環,k表示第k代復制循環,l表示第l代驅散循環,表示選定的前進步長,為生成的隨即向量,表示細菌選定的前進方向。

復制操作:生物進化的規律是優勝劣汰,細菌在進行一定次數的趨向性操作之后會進行健康值排序,這里的健康值就是細菌在之前經歷的不同位置的適度值的累積值,根據排序結果,淘汰掉尋優效果不好的一半數量的細菌,再將剩下的細菌進行分裂,子細菌將繼承母細菌的步長及方向,通過這種方式既保證了新產生的細菌具有良好的尋優能力,同時也保證了菌群規模不變。

驅散操作:細菌的趨向性操作過程有時會導致算法陷入局部最優,為了解決這一問題在細菌覓食優化算法中引入了驅散操作。驅散操作就是將細菌以一定的概率驅散到搜索空間任意一個新的位置,將細菌在新位置的適度值與原位置適度值比較,若優于原位置,則說明經過驅散操作后細菌的全局尋優能力加強。

3 利用粒子群算法優化細菌覓食算法

細菌覓食算法可以在搜索最優解的過程中任意改變搜索方向,這樣大大的提高了算法局部搜索能力及搜索精度,但是細菌在趨向性操作中隨機翻轉,缺少細菌之間的相互學習,因此收斂速度較慢,針對這一不足,本文將粒子群中的粒子更新引入到細菌覓食算法中的細菌位置更新,由于在粒子群算法中粒子是根據本身的最佳位置和整個群體最佳位置來進行位置更新,這樣就可以利用粒子群算法中粒子群的記憶功能提高細菌覓食算法的搜索效率。在粒子群算法中,粒子群更新方式如下所示:

(3)

(4)

其中表示粒子本身最優解,表示整個粒子群的最優解,、表示加速因子,、為[0,1]之間的隨機數。

4 粒子群優化細菌覓食算法的具體步驟

(1)初始化細菌覓食算法以及粒子群算法相關參數,計算細菌的初始適度值。(2)趨向性操作:細菌按照公式(1)前進,如果前進后細菌所在位置的適度值優于前一個位置,則保持在該方向上繼續前進,直至達到設定游動次數,否則,細菌將通過翻轉改變原來前進方向,在每次的趨向性操作中記錄細菌本身經歷過的位置最優值以及全局位置最優值,并根據公式(3)來確定細菌的今后前進方向。(3)復制操作:將細菌在趨向性環節得到的適度值進行累加,累加結果記為,其中,對進行排序,淘汰掉一半細菌,將剩下的細菌進行分裂。(4)驅散操作:將符合驅散操作條件的細菌進行驅散。(5)完成以上流程,將最終的最優適度值及其所對應的細菌所在位置作為結果輸出。

5 系統的方案設計

本文利用粒子群優化的細菌覓食算法來整定PID參數,將優化后的參數應用到變風量空調機組送風管道靜壓控制回路中,具體的控制框圖如圖2所示。

6 系統仿真

粒子群優化的細菌覓食算法各個參數設置如下:細菌數量S設為20,參數維數P為3維,驅散次數為2次,繁殖次數為4次,趨化次數為40次,驅散概率為0.25,游動步數為4次,粒子群算法中的取2,取2。本次仿真選用的尋優函數為絕對偏差積分:

(5)

變風量空調機組的送風管道靜壓控制回路采用一階傳遞函數模型,具體表達式如下:

(6)

本次仿真的采樣周期為1s,采樣點數為100個,根據采樣時間以及傳遞函數表達式可以知道系統傳輸過程中有延遲,延遲約五個采樣周期,用Matlab編寫控制方法的m程序,得到相關PID參數進行仿真,其仿真結果如圖3所示。

粒子群算法,細菌覓食算法,粒子群優化的細菌覓食算法對應的PID值以及在本次仿真過程中的調節時間如表1所示。

根據以上圖表,在調節過程中粒子群算法出現超調量,細菌覓食算法及粒子群優化的細菌覓食算法沒有出現超調量。從調節時間角度看,粒子群算法和細菌覓食優化算法調節時間基本一致,經過粒子群優化的細菌覓食算法的調節時間較粒子群算法及細菌覓食算法有明顯縮短。

7 結語

通過以上分析,粒子群優化的細菌覓食算法與細菌覓食算法以及粒子群算法比較,經過優化后的細菌覓食算法在收斂速度上有明顯提升,調節時間縮短。

參考文獻:

[1]黃偉峰,林衛星,范懷科等.細菌覓食優化的智能PID控制[J].計算機工程與應用,2011,47(21):82-85.

[2]K. M. Passino.Biomimicry of Bacterial Foraging for Distributed Optimization and Control.IEEE.Control System Magazine,June,2002,pp:52-67.

[3]楊世忠,任慶昌.變風量空調靜壓的魯棒PID控制[J].控制工程,2013,20(6):1176-1180.

[4]劉小龍.細菌覓食優化算法的改進及應用[D].廣東:華南理工大學,2011.endprint

【摘 要】 本文主要研究對象是變風量空調機組中的送風管道靜壓控制回路,將具有良好全局尋優能力的細菌覓食算法應用到該控制回路中,同時為了改善細菌覓食算法收斂速度較慢的缺點,將粒子群算法引入到細菌覓食算法中,對細菌覓食算法中的細菌位置更新進行優化,并利用優化后的細菌覓食算法對PID控制器的三個參數進行整定,將整定后的參數應用到控制回路中。通過matlab仿真,并將細菌覓食算法、粒子群算法,粒子群優化的細菌覓食算法進行比較,結果表明經過粒子群優化的細菌覓食算法收斂速度明顯加快。

【關鍵詞】 細菌覓食算法 粒子群算法 PID 送風管道靜壓

變風量空調憑借優越的性能在近年來得到了廣泛的關注,空調機組作為變風量空調系統中的重要設備,其運行直接影響到變風量空調的品質和節能效果。本文針對變風量空調機組送風管道靜壓控制回路進行研究,控制方法主要有定靜壓控制和變靜壓控制,其中定靜壓控制方法運行穩定,控制簡單,易于實現,是應用比較廣泛的方法。這種方法主要是通過調節風機轉速來實現風道靜壓恒定。PID控制器廣泛應用于控制領域,但是由于控制對象越來越復雜,人們對PID的控制要求也在不斷升高,本文利用粒子群優化的細菌覓食算法來整定PID參數,將優化后的參數應用到變風量空調機組送風管道靜壓控制回路中,取得滿意效果。

1 變風量空調機組概述

變風量空氣處理機組是變風量空調系統的關鍵設備,其主要功能是為空調系統提供足量的新風,維持送風溫度和送風管道靜壓在其設定值等。本文將針對送風管道靜壓控制回路進行研究,期望的目標是維持送風管道靜壓在設定值,這一過程是通過調節空氣處理機組中的風機轉速來實現,當靜壓實際值小于設定值,通過控制器調節變頻器,使風機轉速加大,反之,則通過控制器調節變頻器使風機轉速減小。送風管道靜壓控制回路框圖如圖1所示。

2 細菌覓食算法

細菌覓食算法由Kevin M. Passino于2002年提出,該算法模仿大腸桿菌在人體腸道內覓食行為,屬于仿生類優化算法,基本原理是將待優化的問題進行編碼,并定義待優化問題的解對應于搜索空間中的細菌狀態。針對具體問題的求解過程為:產生初始解群體、計算評價函數的值、利用群體的相互影響和作用機制進行優化[1]。細菌覓食算法主要通過趨向性操作、復制操作和驅散操作這三種方式的迭代計算來求解問題。

趨向性操作:趨向性操作是大腸桿菌向著食物豐富的地區前進的過程,這個過程的實現是通過大腸桿菌鞭毛的游動,當鞭毛逆時針時,大腸桿菌就會朝一個方向向前游動,這個過程稱為前進;當鞭毛順時針時,它就會減速直至停止然后重新選擇新的前進方向,這個過程稱為翻轉。大腸桿菌正是通過不斷的前進和翻轉來實現位置的更新,最終找到最優位置。細菌的趨向性操作可表示為如下形式:

(1)

(2)

其中,表示個體i的位置,其中j表示第j代趨向性循環,k表示第k代復制循環,l表示第l代驅散循環,表示選定的前進步長,為生成的隨即向量,表示細菌選定的前進方向。

復制操作:生物進化的規律是優勝劣汰,細菌在進行一定次數的趨向性操作之后會進行健康值排序,這里的健康值就是細菌在之前經歷的不同位置的適度值的累積值,根據排序結果,淘汰掉尋優效果不好的一半數量的細菌,再將剩下的細菌進行分裂,子細菌將繼承母細菌的步長及方向,通過這種方式既保證了新產生的細菌具有良好的尋優能力,同時也保證了菌群規模不變。

驅散操作:細菌的趨向性操作過程有時會導致算法陷入局部最優,為了解決這一問題在細菌覓食優化算法中引入了驅散操作。驅散操作就是將細菌以一定的概率驅散到搜索空間任意一個新的位置,將細菌在新位置的適度值與原位置適度值比較,若優于原位置,則說明經過驅散操作后細菌的全局尋優能力加強。

3 利用粒子群算法優化細菌覓食算法

細菌覓食算法可以在搜索最優解的過程中任意改變搜索方向,這樣大大的提高了算法局部搜索能力及搜索精度,但是細菌在趨向性操作中隨機翻轉,缺少細菌之間的相互學習,因此收斂速度較慢,針對這一不足,本文將粒子群中的粒子更新引入到細菌覓食算法中的細菌位置更新,由于在粒子群算法中粒子是根據本身的最佳位置和整個群體最佳位置來進行位置更新,這樣就可以利用粒子群算法中粒子群的記憶功能提高細菌覓食算法的搜索效率。在粒子群算法中,粒子群更新方式如下所示:

(3)

(4)

其中表示粒子本身最優解,表示整個粒子群的最優解,、表示加速因子,、為[0,1]之間的隨機數。

4 粒子群優化細菌覓食算法的具體步驟

(1)初始化細菌覓食算法以及粒子群算法相關參數,計算細菌的初始適度值。(2)趨向性操作:細菌按照公式(1)前進,如果前進后細菌所在位置的適度值優于前一個位置,則保持在該方向上繼續前進,直至達到設定游動次數,否則,細菌將通過翻轉改變原來前進方向,在每次的趨向性操作中記錄細菌本身經歷過的位置最優值以及全局位置最優值,并根據公式(3)來確定細菌的今后前進方向。(3)復制操作:將細菌在趨向性環節得到的適度值進行累加,累加結果記為,其中,對進行排序,淘汰掉一半細菌,將剩下的細菌進行分裂。(4)驅散操作:將符合驅散操作條件的細菌進行驅散。(5)完成以上流程,將最終的最優適度值及其所對應的細菌所在位置作為結果輸出。

5 系統的方案設計

本文利用粒子群優化的細菌覓食算法來整定PID參數,將優化后的參數應用到變風量空調機組送風管道靜壓控制回路中,具體的控制框圖如圖2所示。

6 系統仿真

粒子群優化的細菌覓食算法各個參數設置如下:細菌數量S設為20,參數維數P為3維,驅散次數為2次,繁殖次數為4次,趨化次數為40次,驅散概率為0.25,游動步數為4次,粒子群算法中的取2,取2。本次仿真選用的尋優函數為絕對偏差積分:

(5)

變風量空調機組的送風管道靜壓控制回路采用一階傳遞函數模型,具體表達式如下:

(6)

本次仿真的采樣周期為1s,采樣點數為100個,根據采樣時間以及傳遞函數表達式可以知道系統傳輸過程中有延遲,延遲約五個采樣周期,用Matlab編寫控制方法的m程序,得到相關PID參數進行仿真,其仿真結果如圖3所示。

粒子群算法,細菌覓食算法,粒子群優化的細菌覓食算法對應的PID值以及在本次仿真過程中的調節時間如表1所示。

根據以上圖表,在調節過程中粒子群算法出現超調量,細菌覓食算法及粒子群優化的細菌覓食算法沒有出現超調量。從調節時間角度看,粒子群算法和細菌覓食優化算法調節時間基本一致,經過粒子群優化的細菌覓食算法的調節時間較粒子群算法及細菌覓食算法有明顯縮短。

7 結語

通過以上分析,粒子群優化的細菌覓食算法與細菌覓食算法以及粒子群算法比較,經過優化后的細菌覓食算法在收斂速度上有明顯提升,調節時間縮短。

參考文獻:

[1]黃偉峰,林衛星,范懷科等.細菌覓食優化的智能PID控制[J].計算機工程與應用,2011,47(21):82-85.

[2]K. M. Passino.Biomimicry of Bacterial Foraging for Distributed Optimization and Control.IEEE.Control System Magazine,June,2002,pp:52-67.

[3]楊世忠,任慶昌.變風量空調靜壓的魯棒PID控制[J].控制工程,2013,20(6):1176-1180.

[4]劉小龍.細菌覓食優化算法的改進及應用[D].廣東:華南理工大學,2011.endprint

主站蜘蛛池模板: 亚洲成aⅴ人在线观看| 欧美中文字幕在线视频| 欧洲极品无码一区二区三区| 亚洲中文字幕国产av| 国产好痛疼轻点好爽的视频| 国产白浆在线| 国产成人综合亚洲网址| 欧美黄色a| 免费a级毛片18以上观看精品| 亚洲AV一二三区无码AV蜜桃| 99久久精品无码专区免费| 欧美午夜理伦三级在线观看 | 免费A级毛片无码免费视频| 日韩精品中文字幕一区三区| 国产欧美在线| 国产成人一区| 19国产精品麻豆免费观看| 一级一级特黄女人精品毛片| 亚洲综合一区国产精品| 老色鬼久久亚洲AV综合| 少妇人妻无码首页| 欧美一级色视频| 青青青视频免费一区二区| 亚洲女同一区二区| 婷婷中文在线| 中文成人在线| 欧美yw精品日本国产精品| 日韩高清欧美| 超碰精品无码一区二区| 麻豆精品在线视频| 中文字幕欧美成人免费| 无码国产伊人| 国产老女人精品免费视频| 欧美亚洲一区二区三区在线| 激情亚洲天堂| 日韩 欧美 小说 综合网 另类| 国产精品第5页| 午夜在线不卡| 免费一级大毛片a一观看不卡| 欧美国产精品不卡在线观看| 天天摸夜夜操| 精品国产免费人成在线观看| 国产精品色婷婷在线观看| 国产流白浆视频| 热久久国产| 亚洲成人网在线播放| 四虎成人精品| 伊人中文网| 午夜福利网址| 国产午夜在线观看视频| 欧美国产日韩在线播放| a级毛片毛片免费观看久潮| 99久久精品视香蕉蕉| 19国产精品麻豆免费观看| 99视频有精品视频免费观看| 国产精品成人AⅤ在线一二三四 | 国语少妇高潮| 成人福利视频网| 久久香蕉国产线| 一级成人a毛片免费播放| 精品日韩亚洲欧美高清a| 在线无码av一区二区三区| 精品国产成人a在线观看| 亚洲有无码中文网| 狼友视频国产精品首页| 亚洲国产精品不卡在线| 久久综合亚洲鲁鲁九月天| 97人妻精品专区久久久久| 久久综合AV免费观看| 五月天久久综合| 99在线视频精品| 老司机精品99在线播放| 香蕉久久国产超碰青草| 国产00高中生在线播放| 国产精品美乳| jizz国产视频| 中字无码av在线电影| 久久久91人妻无码精品蜜桃HD| 国产欧美自拍视频| www.youjizz.com久久| 亚洲天堂精品视频| 日日拍夜夜嗷嗷叫国产|