模糊pid控制方法和裝置及系統的制作方法

            文檔序號:10593169閱讀:640來源:國知局
            模糊pid控制方法和裝置及系統的制作方法
            【專利摘要】本發明公開了一種模糊PID控制方法和裝置及系統。該方法包括:通過第一預設算法計算PID控制參數;通過第二預設算法計算模糊算法參數,其中,模糊算法參數用于計算模糊控制參數,其中,第一預設算法和/或第二預設算法為單變量分布估計算法;獲取被控對象的輸出量和目標量,其中,被控對象為模糊PID控制器的控制對象,被控對象的目標量為被控對象預設的目標輸出量;根據被控對象的輸出量和目標量、PID控制參數和模糊算法參數確定被控對象的控制參數;以及向被控對象輸入確定的控制參數。通過本發明,解決了相關技術中的模糊PID控制方法的控制性能較差的問題。
            【專利說明】
            模糊PID控制方法和裝置及系統
            技術領域
            [0001 ]本發明設及控制領域,具體而言,設及一種模糊PID控制方法和裝置及系統。
            【背景技術】
            [0002] 模糊PID控制方法是控制科學中的經典控制方法,模糊PID控制方法通過模糊PID 控制器對被控量進行控制,模糊PID控制器的輸入量為被控量與目標量的差值,模糊PID控 制器通過輸入量的變化對模糊PID控制器的控制參數進行調整,達到使模糊PID控制器的被 控量跟隨目標量的目的。但是,現有的模糊PID控制方法仍存在一定的弊端,對大滯后的被 控量的控制性能較差。
            [0003] 例如,在對大滯后的火力發電廠機組主蒸汽溫度和壓力進行控制時,存在響應時 間較長、超調量過大、穩態誤差較大等缺點。火力發電廠是通過燃燒產生電能的工廠,燃料 在燃燒時加熱水生成蒸汽,蒸汽通過壓力推動汽輪機旋轉,汽輪機帶動發電機旋轉產生電 能,實現燃料的化學能一熱能一機械能一電能的轉換過程。火力發電廠通過火電機組主蒸 汽溫度壓力控制系統對火電機組蒸汽的溫度和壓力進行控制。
            [0004] 鍋爐汽包出口的飽和蒸汽經由過熱器與爐膛的高溫煙氣進行熱交換,最終加熱為 汽輪發電機組所需的具有合格的溫度壓力品質的過熱蒸汽,其中,過熱器是指從汽包出口 到汽輪機入口前的蒸汽加熱管道。大型鍋爐的過熱器工作于高溫、高壓的惡劣環境中,尤其 是鍋爐過熱器出口處的過熱蒸汽在整個汽水行程中溫度最高。對鍋爐蒸汽溫度的控制品質 是考核火電廠機組運行質量的首選評價指標之一,對火力發電機組的安全經濟生產有著重 要的影響。如果主蒸汽溫度過高,可能會使過熱器管壁的金屬強度降低,損傷過熱器管道的 高溫段,甚至可能會燒壞汽輪機高壓缸設備,嚴重影響安全;如果主蒸汽溫度過低,經過汽 輪機高壓缸的過熱蒸汽濕度增加,會腐蝕高壓缸的葉片,降低汽輪機的輸出功率。一般說 來,大型火電機組對蒸汽溫度的暫時偏差要求為設定值±10°c W內,長期偏差要求為±5°C W內,運對火電廠機組主蒸汽溫度控制系統的要求相當高。由于鍋爐爐膛內燃燒工況的擾 動,如風量擾動、燃料量擾動等,造成沿著過熱器整個管道的煙氣傳熱量幾乎同時變化,所 W過熱器出口主蒸汽溫度的動態特性反應較快,其特征時間和特征滯后時間均比其它擾動 小。因此,在工業生產過程中,主蒸汽溫度控制系統需滿足克服各種干擾的能力,使主蒸汽 溫度在允許的范圍內變動,W保證過熱器溫度不超過允許值。但是主蒸汽溫度控制系統的 本質是利用一個動態特性反應較緩慢的減溫水擾動來克服動態特性反應較快的過熱器吸 熱擾動,運是鍋爐蒸汽溫度品質通常難W達到較高品質的根本原因。
            [0005] 鍋爐蒸汽壓力作為表征鍋爐運行狀態的重要參數,不僅直接關系到鍋爐設備的安 全運行,而且其穩定狀況反映了燃燒過程中能量的供求關系,主蒸汽壓力與機組的發電功 率。通常采用協調控制鍋爐跟隨方式(Coordinated Control of Boiler Following Mode, 簡稱CCBF)或協調控制汽輪機跟隨方式(Coordinated Control of Steam Turbine Following Mode,簡稱CCTF)兩種方式對機組主蒸汽壓力進行控制。
            [0006] 由于火電廠中機組復雜度增加,主蒸汽溫度控制系統的不確定性,非線性,大慣 性,大時延的特征更加突出,常規的PID控制已經難W滿足主蒸汽溫度控制系統對控制品質 的要求。
            [0007] 而蒸汽壓力控制系統具有時變、大滯后等特點,常規的PID算法由于每次輸出均與 過去的狀態有關,計算時要對誤差進行累加,會導致計算輸出的控制量大幅度變化而引起 執行機構的大幅度變化,運種情況在實際生產過程中是不允許的。
            [0008] 針對相關技術中的模糊PID控制方法的控制性能較差的問題,目前尚未提出有效 的解決方案。

            【發明內容】

            [0009] 本發明的主要目的在于提供一種模糊PID控制方法和裝置及系統,W解決相關技 術中的模糊PID控制方法的控制性能較差的問題。
            [0010] 為了實現上述目的,根據本發明的一個方面,提供了一種模糊PID控制方法。該方 法包括:通過第一預設算法計算PID控制參數;通過第二預設算法計算模糊算法參數,其中, 模糊算法參數用于計算模糊控制參數,其中,第一預設算法和/或第二預設算法為單變量分 布估計算法;獲取被控對象的輸出量和目標量,其中,被控對象為模糊PID控制器的控制對 象,被控對象的目標量為被控對象預設的目標輸出量;W及根據被控對象的輸出量和目標 量、PID控制參數和模糊算法參數確定被控對象的控制參數;向被控對象輸入確定的控制參 數。
            [0011] 進一步地,第一預設算法為單變量分布估計算法,在通過第一預設算法計算PID控 制參數之前,該方法還包括:確定PID控制參數的初始染色體種群,其中,PID控制參數的初 始染色體種群包括多個染色體,多個染色體用于在通過預設解碼算法進行解碼之后確定 PID控制參數的初始值,通過第一預設算法計算PID控制參數包括:通過適應度函數計算PID 控制參數的初始染色體種群中各個染色體的適應度值;根據適應度值的大小對PID控制參 數的初始染色體種群中各個染色體進行排序;通過預設更新算法更新在PID控制參數的初 始染色體種群中排序在預設序號之后的染色體,得到PID控制參數的優化染色體種群;判斷 是否滿足停止單變量分布估計算法的條件;如果判斷出滿足停止單變量分布估計算法的條 件,通過預設解碼算法對PID控制參數的優化染色體種群進行解碼,得到PID控制參數;W及 如果判斷出不滿足停止單變量分布估計算法的條件,將PID控制參數的優化染色體種群作 為PID控制參數的初始染色體種群,重新執行通過第一預設算法計算PID控制參數。
            [0012] 進一步地,第二預設算法為單變量分布估計算法,在通過第二預設算法計算模糊 算法參數之前,該方法還包括:確定模糊算法參數的初始染色體種群,其中,模糊算法參數 的初始染色體種群包括多個染色體,多個染色體用于在通過預設解碼算法進行解碼之后確 定模糊算法參數的初始值,通過第二預設算法計算模糊算法參數包括:通過適應度函數計 算模糊算法參數的初始染色體種群中各個染色體的適應度值;根據適應度值的大小對模糊 算法參數的初始染色體種群中各個染色體進行排序;通過預設更新算法更新在模糊算法參 數的初始染色體種群中排序在預設序號之后的染色體,得到模糊算法參數的優化染色體種 群;判斷是否滿足停止單變量分布估計算法的條件;如果判斷出滿足停止單變量分布估計 算法的條件,通過預設解碼算法對模糊算法參數的優化染色體種群進行解碼,得到模糊算 法參數;W及如果判斷出不滿足停止單變量分布估計算法的條件,將模糊算法參數的優化 染色體種群作為模糊算法參數的初始染色體種群,重新執行通過第二預設算法計算模糊算 法參數。
            [0013] 進一步地,在通過適應度函數計算模糊算法參數的初始染色體種群中各個染色體 的適應度值之后,該方法還包括:更新模糊算法參數的歷史最優染色體種群中的染色體,其 中,模糊算法參數的歷史最優染色體種群中的各個染色體的值為模糊算法參數的初始染色 體種群中各個染色體的歷史記錄中適應度值最大的值,確定模糊算法參數的初始染色體種 群包括:根據模糊算法參數的歷史最優染色體種群確定模糊算法參數的初始染色體種群。
            [0014] 進一步地,在根據被控對象的輸出量和目標量、PID控制參數和模糊算法參數向被 控對象輸入總控制量之前,該方法還包括:通過模糊規則系統生成模糊規則庫,其中,模糊 規則系統用于生成模糊規則庫,模糊規則庫包括多個模糊規則,模糊規則庫用于根據被控 對象的輸出量和目標量確定對應的模糊規則,對應的模糊規則用于確定模糊控制參數的模 糊量,模糊控制參數的模糊量用于根據模糊算法參數計算模糊控制參數。
            [0015] 進一步地,獲取被控對象的輸出量和目標量包括:獲取主蒸汽參數采樣值,其中, 主蒸汽參數采樣值是對主蒸汽參數進行采樣得到的采樣值,其中,主蒸汽參數是鍋爐的主 蒸汽的溫度參數或壓力參數,主蒸汽參數是被控對象,主蒸汽參數采樣值是被控對象的輸 出量;W及獲取主蒸汽參數設定值,其中,主蒸汽參數設定值為被控對象的目標量。
            [0016] 進一步地,根據被控對象的輸出量和目標量、PID控制參數和模糊算法參數確定被 控對象的控制參數包括:根據被控對象的輸出量和目標量和模糊算法參數計算模糊控制參 數;根據被控對象的輸出量和目標量和模糊控制參數確定被控對象的模糊控制參數;根據 被控對象的輸出量和目標量和PID控制參數確定被控對象的PID控制參數;W及將確定的模 糊控制參數和PID控制參數相加,得到被控對象的控制參數。
            [0017] 進一步地,在通過第一預設算法計算PID控制參數之前,該方法還包括:獲取被控 對象的傳遞函數;確定PID控制參數的初始值,在通過第二預設算法計算模糊算法參數之 前,該方法還包括:確定模糊算法參數的初始值,通過第一預設算法計算PID控制參數包括: 通過第一評價算法計算PID控制參數的初始值對被控對象的傳遞函數的控制指標;通過第 一優化算法對PID控制參數的初始值進行優化,得到PID控制參數,通過第二預設算法計算 模糊算法參數包括:通過第二評價算法計算模糊算法參數的初始值對被控對象的傳遞函數 的控制指標;通過第二優化算法對模糊算法參數的初始值進行優化,得到模糊算法參數。
            [0018] 為了實現上述目的,根據本發明的一個方面,提供了一種模糊PID控制裝置。該裝 置包括:第一計算單元,用于通過第一預設算法計算PID控制參數;第二計算單元,用于通過 第二預設算法計算模糊算法參數,其中,模糊算法參數用于計算模糊控制參數,其中,第一 預設算法和/或第二預設算法為單變量分布估計算法;獲取單元,用于獲取被控對象的輸出 量和目標量,其中,被控對象為模糊PID控制器的控制對象,被控對象的目標量為被控對象 預設的目標輸出量;確定單元,用于根據被控對象的輸出量和目標量、PID控制參數和模糊 算法參數確定被控對象的控制參數;W及輸入單元,用于向被控對象輸入確定的控制參數。
            [0019] 為了實現上述目的,根據本發明的一個方面,提供了一種鍋爐主蒸汽溫度和壓力 的模糊PID控制系統。該系統包括本發明提供的模糊PID控制裝置,用于向被控對象輸入確 定的控制參數,其中,被控對象為鍋爐的主蒸汽溫度參數和壓力參數,W及采樣器,用于采 集被控對象的輸出量。
            [0020] 本發明通過單變量分布估計算法計算PID控制參數和/或模糊算法參數,W對PID 控制器和模糊控制器的控制參數進行自整定優化,解決了相關技術中的模糊PID控制方法 的控制性能較差的問題,進而達到了提高模糊PID控制方法的控制性能的效果。
            【附圖說明】
            [0021] 構成本申請的一部分的附圖用來提供對本發明的進一步理解,本發明的示意性實 施例及其說明用于解釋本發明,并不構成對本發明的不當限定。在附圖中:
            [0022] 圖1是根據本發明實施例的模糊PID控制方法的流程圖;
            [0023] 圖2是根據本發明實施例的單變量分布估計算法的流程圖;
            [0024] 圖3是根據本發明實施例的模糊算法的流程圖;
            [0025] 圖4是根據本發明實施例的PID控制參數和模糊算法參數優化方法的流程圖;
            [0026] 圖5是根據PID控制方法的仿真控制結果示意圖;
            [0027] 圖6是根據常規的模糊PID控制方法的仿真控制結果示意圖;
            [002引圖7是根據本發明實施例的模糊PID控制方法的仿真控制結果示意圖;
            [0029] 圖8是根據本發明實施例的模糊PID控制裝置的示意圖;W及
            [0030] 圖9是根據本發明實施例的模糊PID控制系統的示意圖。
            【具體實施方式】
            [0031] 需要說明的是,在不沖突的情況下,本申請中的實施例及實施例中的特征可W相 互組合。下面將參考附圖并結合實施例來詳細說明本發明。
            [0032] 為了使本技術領域的人員更好地理解本申請方案,下面將結合本申請實施例中的 附圖,對本申請實施例中的技術方案進行清楚、完整地描述,顯然,所描述的實施例僅僅是 本申請一部分的實施例,而不是全部的實施例。基于本申請中的實施例,本領域普通技術人 員在沒有做出創造性勞動前提下所獲得的所有其他實施例,都應當屬于本申請保護的范 圍。
            [0033] 需要說明的是,本申請的說明書和權利要求書及上述附圖中的術語"第一"、"第 二"等是用于區別類似的對象,而不必用于描述特定的順序或先后次序。應該理解運樣使用 的數據在適當情況下可W互換,W便運里描述的本申請的實施例。此外,術語"包括"和"具 有"W及他們的任何變形,意圖在于覆蓋不排他的包含,例如,包含了一系列步驟或單元的 過程、方法、系統、產品或設備不必限于清楚地列出的那些步驟或單元,而是可包括沒有清 楚地列出的或對于運些過程、方法、產品或設備固有的其它步驟或單元。
            [0034] 本發明的實施例提供了一種模糊PID控制方法。
            [0035] 圖1是根據本發明實施例的模糊PID控制方法的流程圖。如圖1所示,該方法包括W 下步驟:
            [0036] 步驟SlOl,通過第一預設算法計算PID控制參數。
            [0037] PID控制參數是PID控制器的系數,PID控制器是一種常用的控制器,根據被控對象 的給定輸出與實際輸出構成的偏差作為PID控制器的輸入,將偏差分別通過比例環節、積分 環節和微分環節的線性組合構成控制器,根據PID控制器的輸入輸出控制參數對被控對象 進行控制。
            [0038] PID控制器的輸出用于輸入被控對象,也即,輸入PID控制器的控制對象,W使被控 對象根據該輸入量相應的輸出一個量,達到對被控對象的控制,其中,PID控制器的輸出為 被控對象的PID控制參數。PID控制器的輸入量是被控對象的輸出量與目標量的偏差值e,其 中,目標量是預設的被控對象的輸出的目標量。
            [0039] PID控制器的連續數學模型如公式1所示,
            [0040]
            (公式 1)
            [0041] 公式1中,PIDnut(t)為PID控制器的輸出,e(t)是被控對象的實際輸出與目標輸出 的偏差,kp是比例系數,Tl是積分系數,Td是微分系數,由被控對象的特性決定,其中,比例系 數、積分系數和微分系數為PID控制器的系數,也即,PID控制參數。
            [0042] 在工程應用中,按照采樣周期T對被控對象進行采樣,則得到PID控制器的離散數 學模型如公式2所示,
            [0044] 公中,UoUj足束k/|、米巧町刻的FiU粒制帯的輸出,e化)是第k個采樣時刻的被 控對象的實際輸出與目標輸出的偏差,積分系I
            [0043] (公式 2)
            [0045] PID控制參數是指PID控制器的數學模型的系數,例如,離散數學模型中的kp、ki和 kd。通過第一預設算法計算PID控制參數即是通過算法確定出PID控制器的數學模型的系 數,W控制被控對象輸出的量。
            [0046] 在通過第一預設算法計算PID控制參數之前,還需要獲取被控對象的傳遞函數。通 過第一預設算法計算PID控制參數可W是通過試驗法,通過人工試驗不同的PID控制參數對 被控對象的控制效果,根據被控對象的輸出與目標輸出的差值來人工調試PID控制參數。
            [0047] 優選地,第一預設算法也可W是在通過第一預設算法計算PID控制參數之前,確定 PID控制參數的初始值,通過第一評價算法計算PID控制參數的初始值對被控對象的傳遞函 數的控制指標;通過第一優化算法對PID控制參數的初始值進行優化,得到PID控制參數。 [004引優選地,上述第一預設算法可W是單變量分布估計算法。通過單變量分布估計算 法確定一個數學模型中待定參數的取值的方法如圖2所示,單變量分布估計算法包括:
            [0049] 步驟SlOll,初始化染色體種群q(0),q(0)是一個包含m個數值的矩陣,每個數值稱 為染色體,每個染色體的取值為0或1。初始化時,優化次數t記為0。
            [0050] 步驟S1012,按照預設解碼方式對染色體種群q(t)解碼。解碼過程可W簡述為:q (t)在解碼之后對應n個參數,其中,第i個參數Wi的取值范圍是WiG [0,1],假設Wi的染色體 是[qi Q2…qn],貝賄Wi ' = qi X 2〇+q2 X 21+... +Qn X ,再將Wi ' 映射到區間[0,1 ],菊
            在步驟S1012之后,得到了待優化的數學模型。
            [0051] 步驟S1013,計算染色體種群q(t)中m個染色體的適應度值。根據待優化的數學模 型計算染色體種群q(t)中各個染色體的適應度值,適應度值是通過適應度函數計算的值, 用于評價各個染色體的適應度。
            [0052] 步驟S1014,選取適應度值最高的前r個染色體保留。在對m個染色體進行排序之 后,保留排序在前r的染色體,其中,r<m。
            [0053] 步驟S1015,用分布估計產生后m-r個染色體。對于排序在r之后的染色體,采用分 布估計算法重新生成m-r個染色體來替換排序在r之后的染色體。
            [0054] 步驟S1016,令優化次數t = t+l,在步驟S1015之后產生的染色體種群為q(t)。
            [0化5] 步驟S1017,判斷t是否等于N,如果t等于N,結束優化過程,如果t不等于N,返回步 驟S1012。
            [0056] 在完成優化過程之后,單變量分布估計算法結束,得到了優化好的參數。
            [0057] 如果通過單變量分布估計算法計算PID控制參數,則在通過第一預設算法計算PID 控制參數之前,確定PID控制參數的初始染色體種群,其中,PID控制參數的初始染色體種群 包括多個染色體,多個染色體用于在通過預設解碼算法進行解碼之后確定PID控制參數的 初始值。
            [0化引例如,PID控制參數包括kp、ki和kdS個參數,對應的染色體種群分別包括化、化和 Nd位二進制碼,染色體種群分別對應kp、ki和kd取值大小。
            [0059] PID控制參數的初始染色體種群可W是通過隨機算法產生的。在確定PID控制參數 的初始染色體種群之后,通過適應度函數計算PID控制參數的初始染色體種群中各個染色 體的適應度值;根據適應度值的大小對PID控制參數的初始染色體種群中各個染色體進行 排序;通過預設更新算法更新在PID控制參數的初始染色體種群中排序在預設序號之后的 染色體,得到PID控制參數的優化染色體種群。
            [0060] 在得到PID控制參數的優化染色體種群之后,判斷是否滿足停止單變量分布估計 算法的條件。判斷條件可W是判斷優化次數是否到達預設次數,也可W判斷優化染色體種 群在通過預設解碼算法解碼之后得到的PID控制參數是否達到預設的控制指標,其中,控制 指標是指響應時間、超調量、穩態誤差等對控制器的控制性能進行評價的指標。
            [0061] 在判斷是否滿足停止單變量分布估計算法的條件之后,如果判斷出滿足停止單變 量分布估計算法的條件,通過預設解碼算法對PID控制參數的優化染色體種群進行解碼,得 至化ID控制參數;如果判斷出不滿足停止單變量分布估計算法的條件,將PID控制參數的優 化染色體種群作為PID控制參數的初始染色體種群,重新執行通過第一預設算法計算PID控 制參數,并將優化次數累加一次。
            [0062] 步驟S102,通過第二預設算法計算模糊算法參數。
            [0063] 模糊算法參數用于計算模糊控制參數,計算模糊控制參數的方式為根據模糊算法 計算模糊控制參數。
            [0064] 模糊控制參數是模糊控制器的參數,是PID控制參數的增量,模糊控制器用于調整 PID控制器的系數。
            [0065] 模糊控制器的輸入量是被控對象的輸出量與目標量的偏差值e,模糊控制器的輸 出用于輸入被控對象,模糊控制器的輸出為被控對象的模糊控制參數,被控對象同時受PID 控制器和模糊控制器的控制,接收模糊控制器輸出的被控對象的控制參數和PID控制器輸 出的被控對象的控制參數的總和。
            [0066] 模糊控制器的數學模型如公式3所示,
            [0067]
            (公式 3)
            [006引公式3中,山化)是第k個采樣時刻的模糊控制器的輸出,A kp、A ki和A kd是kp、ki和 kd的增量。通過第二預設算法計算模糊控制參數即是通過算法確定出模糊控制器的數學模 型的系數,W調整被控對象輸出的量,其中,模糊控制參數可W是A kp、A ki和A kd。
            [0069] 模糊算法是一種常用的算法,包括模糊化、模糊推理和反模糊化=個計算過程。在 確定出模糊算法參數之后,可W通過模糊算法計算出模糊控制參數。
            [0070] 模糊化過程是對模糊控制器的輸入量進行模糊化得到輸入量的模糊量的過程,模 糊化過程需要兩個輸入量,一個是被控對象的輸出量與目標量的偏差值e,另一個是被控對 象的輸出量與目標量的偏差的變化率e。。對于離散系統來說,e。是當前采樣時刻的被控對象 的輸出量與目標量的偏差值e化)與當前時刻的上一時刻的被控對象的輸出量與目標量的 偏差值e化-1)的差值。
            [0071] 首先確定對兩個輸入量模糊化的目標,也即,將兩個輸入量定義到n個模糊子集 中,每個模糊子集用一個模糊語言值來表示。例如,分別將e和e。定義到7個模糊子集{NB, 醒,NS,ZO,PS,PM,PB}中,其中,NB,醒,NS,ZO,PS,PM和PB為模糊語言值,分別表示負大,負 中,負小,零,正小,正中,正大。將兩個輸入量定義到n個模糊子集中是采用隸屬度函數進行 定義的,隸屬度函數可W是均勻隸屬度函數(公式4)、S角形隸屬度函數(公式5)或者高斯 隸屬度函數(公式6)。
            [0072] (:勞式 4):
            [0073] (公式 5)
            [麵] (公式6)
            [0075] 其中,X是輸入變量,a、b、c是不同隸屬度函數的參數,通過調整參數可W改變函數 形狀,W適應不同系統不同輸入變量在不同值域時的控制要求。模糊算法參數可W包括隸 屬度函數標識和/或隸屬度函數參數。隸屬度函數標識為隸屬度函數的標識,可W用數字代 碼來表示,例如,00表示均勻隸屬度函數,Ol表示=角形隸屬度函數,10表示高斯隸屬度函 數。隸屬度函數參數為隸屬度函數的參數。
            [0076] 在將輸入量模糊化之后,根據模糊算法的輸入量的模糊量通過模糊推理過程推出 模糊控制參數的模糊量。模糊算法的輸入量的模糊量為采用模糊語言值表示的模糊算法的 輸入量,模糊控制參數的模糊量為采用模糊語言值表示的模糊控制參數,模糊控制參數的 模糊量可W在反模糊化過程之后求解出模糊控制參數,其中,反模糊化過程需要用到隸屬 度函數和隸屬度函數參數。
            [0077] 模糊推理過程為根據模糊規則庫確定模糊控制參數的模糊量的過程。模糊規則庫 包括多個模糊規則,模糊規則庫用于根據模糊算法的輸入量的模糊量確定對應的模糊規 貝1J,對應的模糊規則用于確定模糊控制參數的模糊量。模糊規則庫可W是通過模糊規則系 統生成的模糊規則庫,模糊規則系統可W是一個專家系統,該專家系統用于生成模糊規則。 [007引舉例而言,分別將e和ec定義到7個模糊子集{NB,匪,NS,ZO,PS,PM,PB}中,其中, NB,醒,NS,ZO,PS,PM和PB分別表示負大,負中,負小,零,正小,正中,正大,模糊規則庫的二 維控制表可W如表1所示。
            [0079] 表1模糊規則庫的A kp的二維控制表
            [0080]
            [0081] 表1為模糊控制參數Akp的二維控制表,根據e和e。的模糊量可W找到對應的模糊 控制參數Akp的模糊量。
            [0082] 表2模糊規則庫的A ki的二維控制表
            [0083]
            [0084] 表2為模糊控制參數A ki的二維控制表,根據e和ec的模糊量可W找到對應的模糊 控制參數Aki的模糊量。
            [0085] 表3模糊規則庫的A kd的二維控制表
            [0086]
            [0087] 表3為模糊控制參數A kd的二維控制表,根據e和ec的模糊量可W找到對應的模糊 控制參數Akd的模糊量。
            [0088] 在確定模糊控制參數的模糊量之后,通過反模糊化過程求解除模糊控制參數。反 模糊化的過程如公式7所示,
            [0089]
            (公式 7)
            [0090] 式中,Ak表示任意一個模糊控制參數,G為該模糊控制參數對應的比例因子,G等 于模糊控制參數的范圍與模糊子集的個數之比,Wi為模糊控制參數的模糊量,化為Wi在各個 模糊子集中對應的隸屬度函數值。
            [0091] 圖3是根據本發明實施例的模糊算法的流程圖,通過模糊算法求解模糊控制參數 的流程如圖3所示。
            [0092] 步驟S1021為通過第二預設算法確定模糊算法參數,步驟S1022在步驟S1021之后 執行,步驟S1022為根據模糊算法參數對e和ec執行模糊化過程,得到e和ec的模糊量,模糊算 法參數可W包括隸屬度函數標識和隸屬度函數參數。步驟S1023為模糊推理過程,具體地, 根據e和e。的模糊量在模糊規則庫中確定模糊控制參數的模糊量,在確定模糊控制參數的 模糊量之后,執行步驟S1024,解模糊化過程,根據模糊算法參數求解模糊控制參數的模糊 量,得到模糊控制參數。
            [0093] 在通過第二預設算法計算模糊算法參數之前,還需要獲取被控對象的傳遞函數。 通過第二預設算法計算模糊算法參數可W是通過試驗法,通過人工試驗不同的模糊算法參 數計算出來的模糊控制參數對被控對象的控制效果,根據被控對象的輸出與目標輸出的差 值來人工調試模糊算法參數。
            [0094] 優選地,在通過第二預設算法計算模糊算法參數之前,確定模糊算法參數的初始 值,通過第二預設算法計算模糊算法參數包括:通過第二評價算法計算模糊算法參數的初 始值對被控對象的傳遞函數的控制指標;通過第二優化算法對模糊算法參數的初始值進行 優化,得到模糊算法參數。
            [00M]在確定PID控制參數和模糊算法參數之后,PID控制參數和模糊算法參數不變,PID 控制器的數學模型的系數不變,模糊控制器的數學模型的系數根據模糊控制器的輸入量的 變化而變化。
            [0096] 優選地,第二預設算法為單變量分布估計算法,在通過第二預設算法計算模糊算 法參數之前,確定模糊算法參數的初始染色體種群,其中,模糊算法參數的初始染色體種群 包括多個染色體,多個染色體用于在通過預設解碼算法進行解碼之后確定模糊算法參數的 初始值。
            [0097] 在確定模糊算法參數的初始染色體種群之后,通過適應度函數計算模糊算法參數 的初始染色體種群中各個染色體的適應度值;根據適應度值的大小對模糊算法參數的初始 染色體種群中各個染色體進行排序;通過預設更新算法更新在模糊算法參數的初始染色體 種群中排序在預設序號之后的染色體,得到模糊算法參數的優化染色體種群;判斷是否滿 足停止單變量分布估計算法的條件;如果判斷出滿足停止單變量分布估計算法的條件,通 過預設解碼算法對模糊算法參數的優化染色體種群進行解碼,得到模糊算法參數;W及如 果判斷出不滿足停止單變量分布估計算法的條件,將模糊算法參數的優化染色體種群作為 模糊算法參數的初始染色體種群,重新執行通過第二預設算法計算模糊算法參數。
            [0098] 優選地,在通過適應度函數計算模糊算法參數的初始染色體種群中各個染色體的 適應度值之后,更新模糊算法參數的歷史最優染色體種群中的染色體,其中,模糊算法參數 的歷史最優染色體種群中的各個染色體的值為模糊算法參數的初始染色體種群中各個染 色體的歷史記錄中適應度值最大的值,確定模糊算法參數的初始染色體種群包括:根據模 糊算法參數的歷史最優染色體種群確定模糊算法參數的初始染色體種群。
            [0099] 步驟S103,獲取被控對象的輸出量和目標量,其中,被控對象為模糊PID控制器的 控制對象,被控對象的目標量為被控對象預設的目標輸出量。
            [0100] 優選地,獲取被控對象的輸出量和目標量可W是獲取主蒸汽參數采樣值,其中,主 蒸汽參數采樣值是對主蒸汽參數進行采樣得到的采樣值,其中,主蒸汽參數是鍋爐的主蒸 汽的溫度參數或壓力參數,主蒸汽參數是被控對象,主蒸汽參數采樣值是被控對象的輸出 量;W及獲取主蒸汽參數設定值,其中,主蒸汽參數設定值為被控對象的目標量。
            [0101] 步驟S104,根據被控對象的輸出量和目標量、PID控制參數和模糊算法參數確定被 控對象的控制參數。
            [0102] 優選地,根據被控對象的輸出量和目標量、PID控制參數和模糊算法參數確定被控 對象的控制參數可W是根據被控對象的輸出量和目標量和模糊算法參數計算模糊控制參 數;根據被控對象的輸出量和目標量和模糊控制參數確定被控對象的模糊控制參數;根據 被控對象的輸出量和目標量和PID控制參數確定被控對象的PID控制參數;W及將確定的模 糊控制參數和PID控制參數相加,得到被控對象的控制參數。
            [0103] 步驟S105,向被控對象輸入確定的控制參數。
            [0104] 向被控對象輸入確定的控制參數,確定的控制參數為確定的模糊控制參數和PID 控制參數的總和,確定的控制參數的值是模糊控制器和PID控制器分別根據被控對象的目 標量與輸出量的偏差所輸出的值的總和。
            [0105] 該實施例提供的模糊PID控制方法,通過第一預設算法計算PID控制參數;通過第 二預設算法計算模糊算法參數,其中,模糊算法參數用于計算模糊控制參數,其中,第一預 設算法和/或第二預設算法為單變量分布估計算法;獲取被控對象的輸出量和目標量,其 中,被控對象為模糊PID控制器的控制對象,被控對象的目標量為被控對象預設的目標輸出 量;W及根據被控對象的輸出量和目標量、PID控制參數和模糊算法參數確定被控對象的控 制參數;向被控對象輸入確定的控制參數,解決了相關技術中的模糊PID控制方法的控制性 能較差的問題,進而達到了提高模糊PID控制方法的控制性能的效果。
            [0106] 圖4是根據本發明實施例的PID控制參數和模糊算法參數優化方法的流程圖。通過 PID控制參數和模糊算法參數優化方法確定出PID控制參數和模糊算法參數之后,可W根據 被控對象的輸出量與目標量的偏差對被控對象執行控制。
            [0107] 如圖4所示,PID控制參數和模糊算法參數優化方法包括W下步驟:
            [0108] 步驟一,產生PID控制參數和模糊算法參數的染色體種群。染色體種群為一個矩 陣,包括多個染色體,每個染色體為一個數值。PID控制參數和模糊算法參數的染色體種群 中包括PID控制參數基因和模糊算法參數基因,其中,基因為染色體種群中的染色體。
            [0109] 步驟二,將PID控制參數基因和模糊算法參數基因分別轉換為基因顯型,也即,通 過預設解碼算法分別對PID控制參數基因和模糊算法參數基因進行解碼,得到PID控制參數 和模糊算法參數。
            [0110] 步驟=,將PID控制參數和模糊算法參數分別賦值到PID控制模塊和模糊控制模 塊,得到參數確定的模糊PID控制系統。
            [0111] 步驟四,在步驟=之前,產生模擬信號。其中,模擬信號分為訓練信號和測試信號。 訓練信號用于作為PID控制器和模糊控制器的輸入信號對PID控制參數和模糊算法參數進 行訓練。
            [0112] 步驟五,在步驟=和步驟四之后,將訓練信號輸入到參數確定的模糊PID控制系 統,根據模糊PID控制系統對訓練新號的響應對PID控制參數基因和模糊算法參數基因進行 適應度評估。
            [0113] 步驟六,判斷是否滿足終止條件。
            [0114] 步驟屯,在步驟六之后,如果判斷出不滿足終止條件,更新概率分布向量,其中, PID控制參數和模糊算法參數的染色體種群根據分布估計算法計算概率分布向量的結果進 行更新,W調整PID控制參數和模糊算法參數。
            [0115] 步驟八,在步驟六之后,如果判斷出滿足終止條件,則停止更新,得到調整好的模 糊PID控制系統。
            [0116] 步驟九,在步驟六之后,如果判斷出滿足終止條件,則更新最優染色體種群,其中, 最優染色體種群中的各個染色體的值是歷史中適應度最高的染色體值。
            [0117] 步驟九,在步驟八之后,將測試信號輸入調整好的模糊PID控制系統,對調整好的 模糊PID控制系統對測試信號的控制性能進行適應度評估。
            [0118] W具有不確定性、非線性、大慣性、大時延特性的火電廠鍋爐的主蒸汽溫度作為被 控對象,被控對象的傳遞函數為
            '對其進行仿真控制,輸入的測試信 號為方波,采用不同的控制方法進行控制,得到結果如下:
            [0119] 圖5是根據PID控制方法的仿真控制結果示意圖。如圖5所示,僅采用PID控制時,雖 然PID控制方法可W快速地進行響應,最終也能夠穩定在目標值上,但是由此帶來了較顯著 的超調,并且震蕩明顯,延長了調整時間。
            [0120] 圖6是根據常規的模糊PID控制方法的仿真控制結果示意圖。常規的模糊PID控制 方法在PID控制的基礎上疊加了模糊控制修正,如圖6所示,常規的模糊PID控制方法不僅快 速地進行了響應,最終也能夠穩定在目標值上,而且基本消除了超調量,與僅采用PID控制 相比明顯縮短了調整時間,但是在模糊控制消除慣性的過程中帶來了一次明顯的震蕩。
            [0121] 圖7是根據本發明實施例的模糊PID控制方法的仿真控制結果示意圖。如圖7所示, 采用單變量分布估計算法對PID控制參數和模糊算法參數進行自整定優化后,不僅可W快 速地進行響應,最終也能夠穩定在目標值上,而且基本消除了超調量,與PID控制方法和常 規的模糊PID控制方法相比明顯縮短了調整時間,并且也消除了模糊PID控制方法在消除慣 性的過程中帶來的一次明顯的震蕩。
            [0122] 需要說明的是,在附圖的流程圖示出的步驟可W在諸如一組計算機可執行指令的 計算機系統中執行,并且,雖然在流程圖中示出了邏輯順序,但是在某些情況下,可W W不 同于此處的順序執行所示出或描述的步驟。
            [0123] 本發明的實施例還提供了一種模糊PID控制裝置。需要說明的是,本發明實施例的 模糊PID控制裝置可W用于執行本發明的模糊PID控制方法。
            [0124] 圖8是根據本發明實施例的模糊PID控制裝置的示意圖。如圖8所示,該裝置包括第 一計算單元10,第二計算單元20,獲取單元30,確定單元40和輸入單元50。
            [0125] 第一計算單元10用于通過第一預設算法計算PID控制參數。第二計算單元20用于 通過第二預設算法計算模糊算法參數,其中,模糊算法參數用于計算模糊控制參數,其中, 第一預設算法和/或第二預設算法為單變量分布估計算法。獲取單元30用于獲取被控對象 的輸出量和目標量,其中,被控對象為模糊PID控制器的控制對象,被控對象的目標量為被 控對象預設的目標輸出量。確定單元40用于根據被控對象的輸出量和目標量、PID控制參數 和模糊算法參數確定被控對象的控制參數。輸入單元50用于向被控對象輸入確定的控制參 數。
            [0126] 該實施例提供的模糊PID控制裝置解決了相關技術中的模糊PID控制方法的控制 性能較差的問題,進而達到了提高模糊PID控制方法的控制性能的效果。
            [0127] 本發明的實施例還提供了一種鍋爐主蒸汽溫度和壓力的模糊PID控制系統。需要 說明的是,本發明實施例的鍋爐主蒸汽溫度和壓力的模糊PID控制系統包括本發明實施例 的模糊PID控制裝置,其中,被控對象為鍋爐的主蒸汽溫度參數和壓力參數,用于向被控對 象輸入確定的控制參數,該系統還包括采樣器,用于采集被控對象的輸出量。
            [0128] 圖9是根據本發明實施例的鍋爐主蒸汽溫度和壓力的模糊PID控制系統的示意圖。 如圖9所示,該系統包括模糊PID控制裝置100和采樣器201。模糊PID控制裝置100包括數字 PID控制模塊101、模糊控制器102、后向差分器1021和單變量分布估計優化模塊103。
            [0129] 數字PID控制模塊101為PID控制器,單變量分布估計優化模塊103對數字PID控制 模塊101的控制參數和計算模糊控制器102的控制參數的模糊算法參數進行自整定優化,確 定出數字PID控制模塊101的控制參數和模糊算法參數。
            [0130] 采樣器201采集被控對象200的輸出量y(t),與被控對象200的目標量r(t)相減,得 到被控對象200的輸出量y(t)與目標量r(t)的偏差e(t)。被控對象200可W是鍋爐的主蒸汽 溫度參數和壓力參數。
            [0131] 數字PID控制模塊101對偏差e(t)進行處理得到被控對象的第一控制參數U0。
            [0132] 后向差分器1021對偏差e(t)進行處理,將當前采樣時刻的偏差e化)和上一采樣時 刻的偏差e化-1)相減得到偏差變化率e。,模糊控制器102根據模糊算法參數、偏差e (t)和偏 差e化-1)得到被控對象的第二控制參數m。
            [0133] 將第一控制參數UO和第二控制參數Ui相加,得到被控對象的控制參數U,將控制參 數U輸入被控對象200 W實現對被控對象200的控制。
            [0134] 顯然,本領域的技術人員應該明白,上述的本發明的各模塊或各步驟可W用通用 的計算裝置來實現,它們可W集中在單個的計算裝置上,或者分布在多個計算裝置所組成 的網絡上,可選地,它們可W用計算裝置可執行的程序代碼來實現,從而,可W將它們存儲 在存儲裝置中由計算裝置來執行,或者將它們分別制作成各個集成電路模塊,或者將它們 中的多個模塊或步驟制作成單個集成電路模塊來實現。運樣,本發明不限制于任何特定的 硬件和軟件結合。
            [0135] W上所述僅為本發明的優選實施例而已,并不用于限制本發明,對于本領域的技 術人員來說,本發明可W有各種更改和變化。凡在本發明的精神和原則之內,所作的任何修 改、等同替換、改進等,均應包含在本發明的保護范圍之內。
            【主權項】
            1. 一種模糊PID控制方法,其特征在于,包括: 通過第一預設算法計算PID控制參數; 通過第二預設算法計算模糊算法參數,其中,所述模糊算法參數用于計算模糊控制參 數,其中,所述第一預設算法和/或所述第二預設算法為單變量分布估計算法; 獲取被控對象的輸出量和目標量,其中,所述被控對象為模糊PID控制器的控制對象, 所述被控對象的目標量為所述被控對象預設的目標輸出量; 根據所述被控對象的輸出量和目標量、所述PID控制參數和所述模糊算法參數確定所 述被控對象的控制參數;以及 向所述被控對象輸入確定的控制參數。2. 根據權利要求1所述的方法,其特征在于,所述第一預設算法為所述單變量分布估計 算法, 在通過第一預設算法計算PID控制參數之前,所述方法還包括:確定所述PID控制參數 的初始染色體種群,其中,所述PID控制參數的初始染色體種群包括多個染色體,所述多個 染色體用于在通過預設解碼算法進行解碼之后確定所述PID控制參數的初始值, 通過第一預設算法計算PID控制參數包括:通過適應度函數計算所述PID控制參數的初 始染色體種群中各個染色體的適應度值;根據所述適應度值的大小對所述PID控制參數的 初始染色體種群中各個染色體進行排序;通過預設更新算法更新在所述PID控制參數的初 始染色體種群中排序在預設序號之后的染色體,得到所述PID控制參數的優化染色體種群; 判斷是否滿足停止所述單變量分布估計算法的條件;如果判斷出滿足停止所述單變量分布 估計算法的條件,通過所述預設解碼算法對所述PID控制參數的優化染色體種群進行解碼, 得到所述PID控制參數;以及如果判斷出不滿足停止所述單變量分布估計算法的條件,將所 述PID控制參數的優化染色體種群作為所述PID控制參數的初始染色體種群,重新執行通過 第一預設算法計算PID控制參數。3. 根據權利要求1所述的方法,其特征在于,所述第二預設算法為所述單變量分布估計 算法, 在通過第二預設算法計算模糊算法參數之前,所述方法還包括:確定所述模糊算法參 數的初始染色體種群,其中,所述模糊算法參數的初始染色體種群包括多個染色體,所述多 個染色體用于在通過預設解碼算法進行解碼之后確定所述模糊算法參數的初始值, 通過第二預設算法計算模糊算法參數包括:通過適應度函數計算所述模糊算法參數的 初始染色體種群中各個染色體的適應度值;根據所述適應度值的大小對所述模糊算法參數 的初始染色體種群中各個染色體進行排序;通過預設更新算法更新在所述模糊算法參數的 初始染色體種群中排序在預設序號之后的染色體,得到所述模糊算法參數的優化染色體種 群;判斷是否滿足停止所述單變量分布估計算法的條件;如果判斷出滿足停止所述單變量 分布估計算法的條件,通過所述預設解碼算法對所述模糊算法參數的優化染色體種群進行 解碼,得到所述模糊算法參數;以及如果判斷出不滿足停止所述單變量分布估計算法的條 件,將所述模糊算法參數的優化染色體種群作為所述模糊算法參數的初始染色體種群,重 新執行通過第二預設算法計算模糊算法參數。4. 根據權利要求3所述的方法,其特征在于, 在通過適應度函數計算所述模糊算法參數的初始染色體種群中各個染色體的適應度 值之后,所述方法還包括:更新所述模糊算法參數的歷史最優染色體種群中的染色體,其 中,所述模糊算法參數的歷史最優染色體種群中的各個染色體的值為所述模糊算法參數的 初始染色體種群中各個染色體的歷史記錄中適應度值最大的值, 確定所述模糊算法參數的初始染色體種群包括:根據所述模糊算法參數的歷史最優染 色體種群確定所述模糊算法參數的初始染色體種群。5. 根據權利要求1所述的方法,其特征在于,在根據所述被控對象的輸出量和目標量、 所述PID控制參數和所述模糊算法參數向所述被控對象輸入總控制量之前,所述方法還包 括: 通過模糊規則系統生成模糊規則庫,其中,所述模糊規則系統用于生成所述模糊規則 庫,所述模糊規則庫包括多個模糊規則,所述模糊規則庫用于根據所述被控對象的輸出量 和目標量確定對應的模糊規則,所述對應的模糊規則用于確定所述模糊控制參數的模糊 量,所述模糊控制參數的模糊量用于根據所述模糊算法參數計算所述模糊控制參數。6. 根據權利要求1所述的方法,其特征在于,獲取所述被控對象的輸出量和目標量包 括: 獲取主蒸汽參數采樣值,其中,所述主蒸汽參數采樣值是對主蒸汽參數進行采樣得到 的采樣值,其中,所述主蒸汽參數是鍋爐的主蒸汽的溫度參數或壓力參數,所述主蒸汽參數 是所述被控對象,所述主蒸汽參數采樣值是所述被控對象的輸出量;以及 獲取主蒸汽參數設定值,其中,所述主蒸汽參數設定值為所述被控對象的目標量。7. 根據權利要求1所述的方法,其特征在于,根據所述被控對象的輸出量和目標量、所 述PID控制參數和所述模糊算法參數確定所述被控對象的控制參數包括: 根據所述被控對象的輸出量和目標量和所述模糊算法參數計算所述模糊控制參數; 根據所述被控對象的輸出量和目標量和所述模糊控制參數確定所述被控對象的模糊 控制參數; 根據所述被控對象的輸出量和目標量和所述PID控制參數確定所述被控對象的PID控 制參數;以及 將確定的模糊控制參數和PID控制參數相加,得到所述被控對象的控制參數。8. 根據權利要求1所述的方法,其特征在于, 在通過第一預設算法計算PID控制參數之前,所述方法還包括:獲取所述被控對象的傳 遞函數;確定所述PID控制參數的初始值, 在通過第二預設算法計算模糊算法參數之前,所述方法還包括:確定所述模糊算法參 數的初始值, 通過第一預設算法計算PID控制參數包括:通過第一評價算法計算所述PID控制參數的 初始值對所述被控對象的傳遞函數的控制指標;通過第一優化算法對所述PID控制參數的 初始值進行優化,得到所述PID控制參數, 通過第二預設算法計算模糊算法參數包括:通過第二評價算法計算所述模糊算法參數 的初始值對所述被控對象的傳遞函數的控制指標;通過第二優化算法對所述模糊算法參數 的初始值進行優化,得到所述模糊算法參數。9. 一種模糊PID控制裝置,其特征在于,包括: 第一計算單元,用于通過第一預設算法計算PID控制參數; 第二計算單元,用于通過第二預設算法計算模糊算法參數,其中,所述模糊算法參數用 于計算模糊控制參數,其中,所述第一預設算法和/或所述第二預設算法為單變量分布估計 算法; 獲取單元,用于獲取被控對象的輸出量和目標量,其中,所述被控對象為模糊PID控制 器的控制對象,所述被控對象的目標量為所述被控對象預設的目標輸出量; 確定單元,用于根據所述被控對象的輸出量和目標量、所述PID控制參數和所述模糊算 法參數確定所述被控對象的控制參數;以及 輸入單元,用于向所述被控對象輸入確定的控制參數。10.-種鍋爐主蒸汽溫度和壓力的模糊PID控制系統,其特征在于,包括權利要求9所述 的模糊PID控制裝置,用于向被控對象輸入確定的控制參數,其中,所述被控對象為鍋爐的 主蒸汽溫度參數和壓力參數,所述系統還包括: 采樣器,用于采集所述被控對象的輸出量。
            【文檔編號】G05B13/04GK105955026SQ201610371439
            【公開日】2016年9月21日
            【申請日】2016年5月30日
            【發明人】王天堃, 黃波波, 彭大天
            【申請人】神華集團有限責任公司, 神華國能集團有限公司, 陜西德源府谷能源有限公司
            網友詢問留言 已有0條留言
            • 還沒有人留言評論。精彩留言會獲得點贊!
            1
            婷婷六月激情在线综合激情,亚洲国产大片,久久中文字幕综合婷婷,精品久久久久久中文字幕,亚洲一区二区三区高清不卡,99国产精品热久久久久久夜夜嗨 ,欧美日韩亚洲综合在线一区二区,99国产精品电影,伊人精品线视天天综合,精品伊人久久久大香线蕉欧美
            亚洲精品1区 国产成人一级 91精品国产欧美一区二区 亚洲精品乱码久久久久久下载 国产精品久久久久久久伊一 九色国产 国产精品九九视频 伊人久久成人爱综合网 欧美日韩亚洲区久久综合 欧美日本一道免费一区三区 夜夜爽一区二区三区精品 欧美日韩高清一区二区三区 国产成人av在线 国产精品对白交换绿帽视频 国产视频亚洲 国产在线欧美精品 国产精品综合网 国产日韩精品欧美一区色 国产日韩精品欧美一区喷 欧美日韩在线观看区一二 国产区精品 欧美视频日韩视频 中文字幕天天躁日日躁狠狠躁97 视频一二三区 欧美高清在线精品一区二区不卡 国产精品揄拍一区二区久久 99久久综合狠狠综合久久aⅴ 亚洲乱码视频在线观看 日韩在线第二页 亚洲精品无码专区在线播放 成人亚洲网站www在线观看 欧美三级一区二区 99久久精品免费看国产高清 91麻豆国产在线观看 最新日韩欧美不卡一二三区 成人在线观看不卡 日韩国产在线 在线亚洲精品 亚洲午夜久久久久中文字幕 国产精品成人久久久久久久 精品国产一区二区在线观看 欧美精品国产一区二区三区 中文在线播放 亚洲第一页在线视频 国产午夜精品福利久久 九色国产 精品国产九九 国产永久视频 久久精品人人做人人综合试看 国产一区二区三区免费观看 亚洲精品国产电影 9999热视频 国产精品资源在线 麻豆久久婷婷国产综合五月 国产精品免费一级在线观看 亚洲国产一区二区三区青草影视 中文在线播放 国产成人综合在线 国产在线观看色 国产亚洲三级 国产片一区二区三区 久久99精品久久久久久牛牛影视 亚洲欧美日韩国产 四虎永久免费网站 国产一毛片 国产精品视频在 九九热在线精品 99精品福利视频 色婷婷色99国产综合精品 97成人精品视频在线播放 精品久久久久久中文字幕 亚洲欧美一区二区三区孕妇 亚洲欧美成人网 日韩高清在线二区 国产尤物在线观看 在线不卡一区二区 91网站在线看 韩国精品福利一区二区 欧美日韩国产成人精品 99热精品久久 国产精品免费视频一区 高清视频一区 精品九九久久 欧美日韩在线观看免费 91欧美激情一区二区三区成人 99福利视频 亚洲国产精品91 久热国产在线 精品久久久久久中文字幕女 国产精品久久久久久久久99热 成人自拍视频网 国产精品视频久久久久久 久久影院国产 国产玖玖在线观看 99精品在线免费 亚洲欧美一区二区三区导航 久久久久久久综合 国产欧美日韩精品高清二区综合区 国产精品视频自拍 亚洲一级片免费 久久久久久九九 国产欧美自拍视频 视频一区二区在线观看 欧美日韩一区二区三区久久 中文在线亚洲 伊人热人久久中文字幕 日韩欧美亚洲国产一区二区三区 欧美亚洲国产成人高清在线 欧美日韩国产码高清综合人成 国产性大片免费播放网站 亚洲午夜综合网 91精品久久一区二区三区 国产无套在线播放 国产精品视频网站 国产成人亚洲精品老王 91在线网站 国产视频97 欧美黑人欧美精品刺激 国产一区二区三区免费在线视频 久久久国产精品免费看 99re6久精品国产首页 久久精品91 国产成人一级 国产成人精品曰本亚洲 日本福利在线观看 伊人成综合网 久久综合一本 国产综合久久久久久 久久精品成人免费看 久久福利 91精品国产91久久久久久麻豆 亚洲精品成人在线 亚洲伊人久久精品 欧美日本二区 国产永久视频 国产一区二 一区二区福利 国产一毛片 亚洲精品1区 毛片一区二区三区 伊人久久大香线蕉综合影 国产欧美在线观看一区 亚洲国产欧洲综合997久久 国产一区二区免费视频 国产91精品对白露脸全集观看 久久亚洲国产伦理 欧美成人伊人久久综合网 亚洲性久久久影院 久久99国产精一区二区三区! 91精品国产欧美一区二区 欧美日韩亚洲区久久综合 日韩精品一二三区 久久久夜色精品国产噜噜 国产在线精品福利91香蕉 久久久久久久亚洲精品 97se色综合一区二区二区 91国语精品自产拍在线观看性色 91久久国产综合精品女同我 日韩中文字幕a 国产成人亚洲日本精品 久久国产精品-国产精品 久久国产经典视频 久久国产精品伦理 亚洲第一页在线视频 国产精品久久久久三级 日韩毛片网 久久免费高清视频 麻豆国产在线观看一区二区 91麻豆国产福利在线观看 国产成人精品男人的天堂538 一区二区三区中文字幕 免费在线视频一区 欧美日韩国产成人精品 国产综合网站 国产资源免费观看 亚洲精品亚洲人成在线播放 精品久久久久久中文字幕专区 亚洲人成人毛片无遮挡 国产一起色一起爱 国产香蕉精品视频在 九九热免费观看 日韩亚洲欧美一区 九九热精品在线观看 精品久久久久久中文字幕专区 亚洲欧美自拍偷拍 国产精品每日更新 久久久久国产一级毛片高清板 久久天天躁狠狠躁夜夜中文字幕 久久精品片 日韩在线毛片 国产成人精品本亚洲 国产成人精品一区二区三区 九九热在线观看 国产r级在线观看 国产欧美日韩精品高清二区综合区 韩国电影一区二区 国产精品毛片va一区二区三区 五月婷婷伊人网 久久一区二区三区免费 一本色道久久综合狠狠躁篇 亚洲综合色站 国产尤物在线观看 亚洲一区亚洲二区 免费在线视频一区 欧洲精品视频在线观看 日韩中文字幕a 中文字幕日本在线mv视频精品 91精品在线免费视频 精品国产免费人成在线观看 精品a级片 中文字幕日本在线mv视频精品 日韩在线精品视频 婷婷丁香色 91精品国产高清久久久久 国产成人精品日本亚洲直接 五月综合视频 欧美日韩在线亚洲国产人 精液呈暗黄色 亚洲乱码一区 久久精品中文字幕不卡一二区 亚洲天堂精品在线 激情婷婷综合 国产免费久久精品久久久 国产精品亚洲二区在线 久久免费播放视频 五月婷婷丁香综合 在线亚洲欧美日韩 久久免费精品高清麻豆 精品久久久久久中文字幕 亚洲一区网站 国产精品福利社 日韩中文字幕免费 亚洲综合丝袜 91精品在线播放 国产精品18 亚洲日日夜夜 伊人久久大香线蕉综合影 亚洲精品中文字幕乱码影院 亚洲一区二区黄色 亚洲第一页在线视频 一区二区在线观看视频 国产成人福利精品视频 亚洲高清二区 国内成人免费视频 精品亚洲性xxx久久久 国产精品合集一区二区三区 97av免费视频 国产一起色一起爱 国产区久久 国产资源免费观看 99精品视频免费 国产成人一级 国产精品九九免费视频 欧美91精品久久久久网免费 99热国产免费 久久精品色 98精品国产综合久久 久久精品播放 中文字幕视频免费 国产欧美日韩一区二区三区在线 精品久久蜜桃 国产小视频精品 一本色道久久综合狠狠躁篇 91在线免费观看 亚洲精品区 伊人成综合网 伊人热人久久中文字幕 伊人黄色片 99国产精品热久久久久久夜夜嗨 久久免费精品视频 亚洲一区二区三区高清不卡 久久久久国产一级毛片高清板 国产片一区二区三区 久久狠狠干 99久久婷婷国产综合精品电影 国产99区 国产精品成人久久久久 久久狠狠干 青青国产在线观看 亚洲高清国产拍精品影院 国产精品一区二区av 九九热在线免费视频 伊人久久国产 国产精品久久久久久久久久一区 在线观看免费视频一区 国产精品自在在线午夜区app 国产精品综合色区在线观看 国产毛片久久久久久国产毛片 97国产免费全部免费观看 国产精品每日更新 国产尤物视频在线 九九视频这里只有精品99 一本一道久久a久久精品综合 久久综合给会久久狠狠狠 国产成人精品男人的天堂538 欧美一区二区高清 毛片一区二区三区 国产欧美日韩在线观看一区二区三区 在线国产二区 欧美不卡网 91在线精品中文字幕 在线国产福利 国内精品91久久久久 91亚洲福利 日韩欧美国产中文字幕 91久久精品国产性色也91久久 亚洲性久久久影院 欧美精品1区 国产热re99久久6国产精品 九九热免费观看 国产精品欧美日韩 久久久久国产一级毛片高清板 久久国产经典视频 日韩欧美亚洲国产一区二区三区 欧美亚洲综合另类在线观看 国产精品自在在线午夜区app 97中文字幕在线观看 视频一二三区 精品国产一区在线观看 国产欧美日韩在线一区二区不卡 欧美一区二三区 伊人成人在线观看 国内精品91久久久久 97在线亚洲 国产在线不卡一区 久久久全免费全集一级全黄片 国产精品v欧美精品∨日韩 亚洲毛片网站 在线不卡一区二区 99re热在线视频 久久激情网 国产毛片一区二区三区精品 久久亚洲综合色 中文字幕视频免费 国产视频亚洲 婷婷伊人久久 国产一区二区免费播放 久久99国产精品成人欧美 99国产在线视频 国产成人免费视频精品一区二区 国产不卡一区二区三区免费视 国产码欧美日韩高清综合一区 久久精品国产主播一区二区 国产一区电影 久久精品国产夜色 国产精品国产三级国产 日韩一区二区三区在线 久久97久久97精品免视看 久久国产免费一区二区三区 伊人久久大香线蕉综合电影网 99re6久精品国产首页 久久激情网 亚洲成人高清在线 国产精品网址 国产成人精品男人的天堂538 香蕉国产综合久久猫咪 国产专区中文字幕 91麻豆精品国产高清在线 久久国产经典视频 国产精品成人va在线观看 国产精品爱啪在线线免费观看 日本精品久久久久久久久免费 亚洲综合一区二区三区 久久五月网 精品国产网红福利在线观看 久久综合亚洲伊人色 亚洲国产精品久久久久久网站 在线日韩国产 99国产精品热久久久久久夜夜嗨 国产综合精品在线 国产区福利 精品亚洲综合久久中文字幕 国产制服丝袜在线 毛片在线播放网站 在线观看免费视频一区 国产精品久久久精品三级 亚洲国产电影在线观看 最新日韩欧美不卡一二三区 狠狠综合久久综合鬼色 日本精品1在线区 国产日韩一区二区三区在线播放 欧美日韩精品在线播放 亚洲欧美日韩国产一区二区三区精品 久久综合久久网 婷婷六月激情在线综合激情 亚洲乱码一区 国产专区91 97av视频在线观看 精品久久久久久中文字幕 久久五月视频 国产成人福利精品视频 国产精品网址 中文字幕视频在线 精品一区二区三区免费视频 伊人手机在线视频 亚洲精品中文字幕乱码 国产在线视频www色 色噜噜国产精品视频一区二区 精品亚洲成a人在线观看 国产香蕉尹人综合在线 成人免费一区二区三区在线观看 国产不卡一区二区三区免费视 欧美精品久久天天躁 国产专区中文字幕 久久精品国产免费中文 久久精品国产免费一区 久久无码精品一区二区三区 国产欧美另类久久久精品免费 欧美精品久久天天躁 亚洲精品在线视频 国产视频91在线 91精品福利一区二区三区野战 日韩中文字幕免费 国产精品99一区二区三区 欧美成人高清性色生活 国产精品系列在线观看 亚洲国产福利精品一区二区 国产成人在线小视频 国产精品久久久久免费 99re热在线视频 久久久久久久综合 一区二区国产在线播放 成人国产在线视频 亚洲精品乱码久久久久 欧美日韩一区二区综合 精品久久久久免费极品大片 中文字幕视频二区 激情粉嫩精品国产尤物 国产成人精品一区二区视频 久久精品中文字幕首页 亚洲高清在线 国产精品亚洲一区二区三区 伊人久久艹 中文在线亚洲 国产精品一区二区在线播放 国产精品九九免费视频 亚洲二区在线播放 亚洲狠狠婷婷综合久久久久网站 亚洲欧美日韩网站 日韩成人精品 亚洲国产一区二区三区青草影视 91精品国产福利在线观看 国产精品久久久久久久久99热 国产一区二区精品尤物 久碰香蕉精品视频在线观看 亚洲日日夜夜 在线不卡一区二区 国产午夜亚洲精品 九九热在线视频观看这里只有精品 伊人手机在线视频 91免费国产精品 日韩欧美中字 91精品国产91久久久久 国产全黄三级播放 视频一区二区三区免费观看 国产开裆丝袜高跟在线观看 国产成人欧美 激情综合丝袜美女一区二区 国产成人亚洲综合无 欧美精品一区二区三区免费观看 欧美亚洲国产日韩 日韩亚州 国产欧美日韩精品高清二区综合区 亚洲午夜国产片在线观看 精品久久久久久中文字幕 欧美精品1区 久久伊人久久亚洲综合 亚洲欧美日韩精品 国产成人精品久久亚洲高清不卡 久久福利影视 国产精品99精品久久免费 久久久久免费精品视频 国产日产亚洲精品 亚洲国产午夜电影在线入口 精品无码一区在线观看 午夜国产精品视频 亚洲一级片免费 伊人久久大香线蕉综合影 国产精品久久影院 久碰香蕉精品视频在线观看 www.欧美精品 在线小视频国产 亚洲国产天堂久久综合图区 欧美一区二区三区不卡 日韩美女福利视频 九九精品免视频国产成人 不卡国产00高中生在线视频 亚洲第一页在线视频 欧美日韩在线播放成人 99re视频这里只有精品 国产精品91在线 精品乱码一区二区三区在线 国产区久久 91麻豆精品国产自产在线观看一区 日韩精品成人在线 九九热在线观看 国产精品久久不卡日韩美女 欧美一区二区三区综合色视频 欧美精品免费一区欧美久久优播 国产精品网址 国产专区中文字幕 国产精品欧美亚洲韩国日本久久 日韩美香港a一级毛片 久久精品123 欧美一区二区三区免费看 99r在线视频 亚洲精品国产字幕久久vr 国产综合激情在线亚洲第一页 91免费国产精品 日韩免费小视频 亚洲国产精品综合一区在线 国产亚洲第一伦理第一区 在线亚洲精品 国产精品一区二区制服丝袜 国产在线成人精品 九九精品免视频国产成人 亚洲国产网 欧美日韩亚洲一区二区三区在线观看 在线亚洲精品 欧美一区二区三区高清视频 国产成人精品男人的天堂538 欧美日韩在线观看区一二 亚洲欧美一区二区久久 久久精品中文字幕首页 日本高清www午夜视频 久久精品国产免费 久久999精品 亚洲国产精品欧美综合 88国产精品视频一区二区三区 91久久偷偷做嫩草影院免费看 国产精品夜色视频一区二区 欧美日韩导航 国产成人啪精品午夜在线播放 一区二区视频在线免费观看 99久久精品国产自免费 精液呈暗黄色 久久99国产精品 日本精品久久久久久久久免费 精品国产97在线观看 99re视频这里只有精品 国产视频91在线 999av视频 亚洲美女视频一区二区三区 久久97久久97精品免视看 亚洲国产成人久久三区 99久久亚洲国产高清观看 日韩毛片在线视频 综合激情在线 91福利一区二区在线观看 一区二区视频在线免费观看 激情粉嫩精品国产尤物 国产成人精品曰本亚洲78 国产成人精品本亚洲 国产精品成人免费视频 国产成人啪精品视频免费软件 久久精品国产亚洲妲己影院 国产精品成人久久久久久久 久久大香线蕉综合爱 欧美一区二区三区高清视频 99热国产免费 在线观看欧美国产 91精品视频在线播放 国产精品福利社 欧美精品一区二区三区免费观看 国产一区二区免费视频 国产午夜精品一区二区 精品视频在线观看97 91精品福利久久久 国产一区福利 国产综合激情在线亚洲第一页 国产精品久久久久久久久久久不卡 九色国产 在线日韩国产 黄网在线观看 亚洲一区小说区中文字幕 中文字幕丝袜 日本二区在线观看 日本国产一区在线观看 欧美日韩一区二区三区久久 欧美精品亚洲精品日韩专 国产日产亚洲精品 久久综合九色综合欧美播 亚洲国产欧美无圣光一区 欧美视频区 亚洲乱码视频在线观看 久久无码精品一区二区三区 九九热精品免费视频 久久99精品久久久久久牛牛影视 国产精品成久久久久三级 国产一区福利 午夜国产精品视频 日本二区在线观看 99久久网站 国产亚洲天堂 精品国产一区二区三区不卡 亚洲国产日韩在线一区 国产成人综合在线观看网站 久久免费高清视频 欧美在线导航 午夜精品久久久久久99热7777 欧美久久综合网 国产小视频精品 国产尤物在线观看 亚洲国产精品综合一区在线 欧美一区二区三区不卡视频 欧美黑人欧美精品刺激 日本福利在线观看 久久国产偷 国产手机精品一区二区 国产热re99久久6国产精品 国产高清啪啪 欧美亚洲国产成人高清在线 国产在线第三页 亚洲综合一区二区三区 99r在线视频 99精品久久久久久久婷婷 国产精品乱码免费一区二区 国产在线精品福利91香蕉 国产尤物视频在线 五月婷婷亚洲 中文字幕久久综合伊人 亚洲精品一级毛片 99国产精品电影 在线视频第一页 久久99国产精品成人欧美 国产白白视频在线观看2 成人精品一区二区www 亚洲成人网在线观看 麻豆91在线视频 色综合合久久天天综合绕视看 久久精品国产免费高清 国产不卡一区二区三区免费视 欧美国产中文 99精品欧美 九九在线精品 国产中文字幕在线免费观看 国产一区中文字幕在线观看 国产成人一级 国产精品一区二区制服丝袜 国产一起色一起爱 亚洲精品成人在线 亚洲欧美精品在线 国产欧美自拍视频 99精品久久久久久久婷婷 久99视频 国产热re99久久6国产精品 视频一区亚洲 国产精品视频分类 国产精品成在线观看 99re6久精品国产首页 亚洲在成人网在线看 亚洲国产日韩在线一区 久久国产三级 日韩国产欧美 欧美在线一区二区三区 国产精品美女一级在线观看 成人午夜免费福利视频 亚洲天堂精品在线 91精品国产手机 欧美日韩视频在线播放 狠狠综合久久综合鬼色 九一色视频 青青视频国产 亚洲欧美自拍一区 中文字幕天天躁日日躁狠狠躁97 日韩免费大片 996热视频 伊人成综合网 亚洲天堂欧美 日韩精品亚洲人成在线观看 久久综合给会久久狠狠狠 日韩精品亚洲人成在线观看 日韩国产欧美 亚洲成aⅴ人片在线影院八 亚洲精品1区 99久久精品免费 国产精品高清在线观看 国产精品久久久免费视频 在线亚洲欧美日韩 91在线看视频 国产精品96久久久久久久 欧美日韩国产成人精品 91在线亚洲 热久久亚洲 国产精品美女免费视频观看 日韩在线毛片 亚洲永久免费视频 九九免费在线视频 亚洲一区网站 日本高清二区视频久二区 精品国产美女福利在线 伊人久久艹 国产精品久久久久三级 欧美成人精品第一区二区三区 99久久精品国产自免费 在线观看日韩一区 国产中文字幕一区 成人免费午夜视频 欧美日韩另类在线 久久99国产精品成人欧美 色婷婷中文网 久久天天躁夜夜躁狠狠躁2020 欧美成人伊人久久综合网 国产精品福利资源在线 国产伦精品一区二区三区高清 国产精品亚洲综合色区韩国 亚洲一区欧美日韩 色综合视频 国语自产精品视频在线区 国产高清a 成人国内精品久久久久影 国产在线精品香蕉综合网一区 国产不卡在线看 国产成人精品精品欧美 国产欧美日韩综合精品一区二区三区 韩国电影一区二区 国产在线视频www色 91中文字幕在线一区 国产人成午夜免视频网站 亚洲综合一区二区三区 色综合视频一区二区观看 久久五月网 九九热精品在线观看 国产一区二区三区国产精品 99久热re在线精品996热视频 亚洲国产网 在线视频亚洲一区 日韩字幕一中文在线综合 国产高清一级毛片在线不卡 精品国产色在线 国产高清视频一区二区 精品日本久久久久久久久久 亚洲国产午夜精品乱码 成人免费国产gav视频在线 日韩欧美一区二区在线观看 欧美曰批人成在线观看 韩国电影一区二区 99re这里只有精品6 日韩精品一区二区三区视频 99re6久精品国产首页 亚洲欧美一区二区三区导航 欧美色图一区二区三区 午夜精品视频在线观看 欧美激情在线观看一区二区三区 亚洲热在线 成人国产精品一区二区网站 亚洲一级毛片在线播放 亚洲一区小说区中文字幕 亚洲午夜久久久久影院 国产自产v一区二区三区c 国产精品视频免费 久久调教视频 国产成人91激情在线播放 国产精品欧美亚洲韩国日本久久 久久亚洲日本不卡一区二区 91中文字幕网 成人国产在线视频 国产视频91在线 欧美成人精品第一区二区三区 国产精品福利在线 久久综合九色综合精品 欧美一区二区三区精品 久久国产综合尤物免费观看 久久99青青久久99久久 日韩精品免费 久久国产精品999 91亚洲视频在线观看 国产精品igao视频 色综合区 在线亚洲欧国产精品专区 国产一区二区三区在线观看视频 亚洲精品成人在线 一区二区国产在线播放 中文在线亚洲 亚洲精品第一国产综合野 国产一区二区精品久久 一区二区三区四区精品视频 99热精品久久 中文字幕视频二区 国产成人精品男人的天堂538 99精品影视 美女福利视频一区二区 久久午夜夜伦伦鲁鲁片 综合久久久久久久综合网 国产精品国产欧美综合一区 国产99视频在线观看 国产亚洲女在线精品 婷婷影院在线综合免费视频 国产亚洲3p一区二区三区 91成人爽a毛片一区二区 亚洲一区二区高清 国产欧美亚洲精品第二区首页 欧美日韩导航 亚洲高清二区 欧美激情观看一区二区久久 日韩毛片在线播放 亚洲欧美日韩高清中文在线 亚洲日本在线播放 国产精品一区二区制服丝袜 精品国产一区二区三区不卡 国产不卡在线看 国产欧美网站 四虎永久在线观看视频精品 国产黄色片在线观看 夜夜综合 一本色道久久综合狠狠躁篇 欧美亚洲综合另类在线观看 国产91在线看 伊人久久国产 欧美一区二区在线观看免费网站 国产精品久久久久三级 久久福利 日韩中文字幕a 亚洲午夜久久久久影院 91在线高清视频 国产亚洲一区二区三区啪 久久人精品 国产精品亚洲午夜一区二区三区 综合久久久久久 久久伊人一区二区三区四区 国产综合久久久久久 日韩一区精品视频在线看 国产精品日韩欧美制服 日本精品1在线区 99re视频 无码av免费一区二区三区试看 国产视频1区 日韩欧美中文字幕一区 日本高清中文字幕一区二区三区a 亚洲国产欧美无圣光一区 国产在线视频一区二区三区 欧美国产第一页 在线亚洲欧美日韩 日韩中文字幕第一页 在线不卡一区二区 伊人久久青青 国产精品一区二区在线播放 www.五月婷婷 麻豆久久婷婷国产综合五月 亚洲精品区 久久国产欧美另类久久久 99在线视频免费 伊人久久中文字幕久久cm 久久精品成人免费看 久久这里只有精品首页 88国产精品视频一区二区三区 中文字幕日本在线mv视频精品 国产在线精品成人一区二区三区 伊人精品线视天天综合 亚洲一区二区黄色 国产尤物视频在线 亚洲精品99久久久久中文字幕 国产一区二区三区免费观看 伊人久久大香线蕉综合电影网 国产成人精品区在线观看 日本精品一区二区三区视频 日韩高清在线二区 久久免费播放视频 一区二区成人国产精品 国产精品免费精品自在线观看 亚洲精品视频二区 麻豆国产精品有码在线观看 精品日本一区二区 亚洲欧洲久久 久久中文字幕综合婷婷 中文字幕视频在线 国产成人精品综合在线观看 91精品国产91久久久久福利 精液呈暗黄色 香蕉国产综合久久猫咪 国产专区精品 亚洲精品无码不卡 国产永久视频 亚洲成a人片在线播放观看国产 一区二区国产在线播放 亚洲一区二区黄色 欧美日韩在线观看视频 亚洲精品另类 久久国产综合尤物免费观看 国产一区二区三区国产精品 高清视频一区 国产精品igao视频 国产精品资源在线 久久综合精品国产一区二区三区 www.五月婷婷 精品色综合 99热国产免费 麻豆福利影院 亚洲伊人久久大香线蕉苏妲己 久久电影院久久国产 久久精品伊人 在线日韩理论午夜中文电影 亚洲国产欧洲综合997久久 伊人国产精品 久草国产精品 欧美一区精品二区三区 亚洲成人高清在线 91免费国产精品 日韩精品福利在线 国产一线在线观看 国产不卡在线看 久久99青青久久99久久 亚洲精品亚洲人成在线播放 99久久免费看国产精品 国产日本在线观看 青草国产在线视频 麻豆久久婷婷国产综合五月 国产中文字幕一区 91久久精品国产性色也91久久 国产一区a 国产欧美日韩成人 国产亚洲女在线精品 一区二区美女 中文字幕在线2021一区 在线小视频国产 久久这里只有精品首页 国产在线第三页 欧美日韩中文字幕 在线亚洲+欧美+日本专区 精品国产一区二区三区不卡 久久这里精品 欧美在线va在线播放 精液呈暗黄色 91精品国产手机 91在线免费播放 欧美视频亚洲色图 欧美国产日韩精品 日韩高清不卡在线 精品视频免费观看 欧美日韩一区二区三区四区 国产欧美亚洲精品第二区首页 亚洲韩精品欧美一区二区三区 国产精品视频免费 在线精品小视频 久久午夜夜伦伦鲁鲁片 国产无套在线播放 久热这里只精品99re8久 欧美久久久久 久久香蕉国产线看观看精品蕉 国产成人精品男人的天堂538 亚洲人成网站色7799在线观看 日韩在线第二页 一本色道久久综合狠狠躁篇 国产一区二区三区不卡在线观看 亚洲乱码在线 在线观看欧美国产 久久福利青草精品资源站免费 国产玖玖在线观看 在线亚洲精品 亚洲成aⅴ人在线观看 精品91在线 欧美一区二三区 日韩中文字幕视频在线 日本成人一区二区 日韩免费专区 国内精品在线观看视频 久久国产综合尤物免费观看 国产精品系列在线观看 一本一道久久a久久精品综合 亚洲免费播放 久久精品国产免费 久久人精品 亚洲毛片网站 亚洲成a人一区二区三区 韩国福利一区二区三区高清视频 亚洲精品天堂在线 一区二区三区中文字幕 亚洲国产色婷婷精品综合在线观看 亚洲国产成人久久笫一页 999国产视频 国产精品香港三级在线电影 欧美日韩一区二区三区四区 日韩国产欧美 国产精品99一区二区三区 午夜国产精品理论片久久影院 亚洲精品中文字幕麻豆 亚洲国产高清视频 久久免费手机视频 日韩a在线观看 五月婷婷亚洲 亚洲精品中文字幕麻豆 中文字幕丝袜 www国产精品 亚洲天堂精品在线 亚洲乱码一区 国产日韩欧美三级 久久999精品 伊人热人久久中文字幕 久热国产在线视频 国产欧美日韩在线观看一区二区三区 国产一二三区在线 日韩国产欧美 91精品国产91久久久久 亚洲一区小说区中文字幕 精品一区二区免费视频 国产精品视频免费 国产精品亚洲综合色区韩国 亚洲国产精品成人午夜在线观看 欧美国产日韩精品 中文字幕精品一区二区精品