亚洲免费av电影一区二区三区,日韩爱爱视频,51精品视频一区二区三区,91视频爱爱,日韩欧美在线播放视频,中文字幕少妇AV,亚洲电影中文字幕,久久久久亚洲av成人网址,久久综合视频网站,国产在线不卡免费播放

        ?

        基于支持向量機的城市小區(qū)燃氣管網(wǎng)日負荷預測模型

        2014-11-10 12:35:57張建秋
        科技創(chuàng)新導報 2014年12期
        關(guān)鍵詞:支持向量機

        張建秋

        摘 要:該文使用支持向量機中的兩種核函數(shù),采用grid-search算法、遺傳算法、粒子群算法優(yōu)化參數(shù),建立對吉林市某小區(qū)燃氣管網(wǎng)日負荷預測的支持向量機模型。將日最高溫度、日最低溫度、日平均溫度、小區(qū)人員最高年齡、小區(qū)人員最低年齡、小區(qū)人員平均年齡作為燃氣管網(wǎng)日負荷變化密切相關(guān)的主要影響因素,分別作為支持向量機的輸入量,將小區(qū)人員臨時出差、小區(qū)臨時增加暫住人口等隨機因素作為燃氣管網(wǎng)日負荷變化密切相關(guān)的次要影響因素,將隨機因素統(tǒng)一歸為支持向量機的一個輸入量。采用[0,1]歸一化方法,對作為影響因素的輸入量數(shù)據(jù)與日負荷預測輸出量數(shù)據(jù)進行歸一化處理。對節(jié)假日和工作日的燃氣管網(wǎng)日負荷預測采用獨立處理方法,避免了相互之間的干擾影響。試驗結(jié)果表明,采用徑向基核函數(shù)的支持向量機預測模型對燃氣管網(wǎng)日負荷預測擬合程度達到90%以上。

        關(guān)鍵詞:燃氣管網(wǎng) 支持向量機 日負荷預測

        中圖分類號:TP391.4 文獻標識碼:A 文章編號:1674-098X(2014)04(c)-0099-04

        城市小區(qū)燃氣管網(wǎng)的日負荷指標是城市燃氣的基礎(chǔ)工作,對分析燃氣管網(wǎng)的年負荷、月負荷具有重要指導意義,并且燃氣經(jīng)營企業(yè)對燃氣管網(wǎng)的日負荷及預測情況日益重視[1]。燃氣管網(wǎng)日負荷數(shù)據(jù)的變化規(guī)律十分復雜,主要同天氣、氣溫、終端用戶人口數(shù)量等因素密切相關(guān),并且國家法定節(jié)假日的日負荷變化與工作日亦不相同。傳統(tǒng)的燃氣管網(wǎng)日負荷預測方法主要有回歸綜合位移平均模型、多元主要有線性回歸法、三角函數(shù)模型、BP神經(jīng)網(wǎng)絡法等[2],試驗結(jié)果表明這些預測方法具有一定的精度[3-5]。但上述方法所需樣本數(shù)據(jù)較多,并且樣本數(shù)據(jù)越多,上述方法得到的預測模型泛化能力越強,可信度越高。

        支持向量機(SVM)是機器學習研究重大成果,其具有泛化能力強、全局尋優(yōu)的特點,尤其具有所需樣本數(shù)據(jù)少的特點[6]。SVM利用松弛變量和核函數(shù)[7],針對樣本數(shù)據(jù)線性不可分的情況,在高維空間尋找其最優(yōu)分類面。

        該文將吉林市某小區(qū)燃氣管網(wǎng)的日負荷作為研究對象,將燃氣管網(wǎng)日負荷變化密切相關(guān)的主要影響因素——日最高溫度、日最低溫度、日平均溫度、小區(qū)人員最高年齡、小區(qū)人員最低年齡、小區(qū)人員平均年齡,作為SVM的6個輸入量,將燃氣管網(wǎng)日負荷變化密切相關(guān)的次要影響因素——小區(qū)人員臨時出差、小區(qū)臨時增加暫住人口等隨機因素,統(tǒng)一作為SVM的1個輸入量,利用不同核函數(shù)、不同參數(shù)的支持向量機建立吉林市某小區(qū)燃氣管網(wǎng)日負荷預測模型,通過該預測模型得到該小區(qū)的燃氣管網(wǎng)日負荷預測值。

        1 支持向量機

        1.1 支持向量機原理

        設(shè)表示線性函數(shù),其表達式為: (1.1)

        其中:為偏置,為可調(diào)權(quán)值向量,決策規(guī)則:,將稱為超平面(Hyper Plane)。由式定義超平面L,L將輸入空間X劃分成兩部分,如圖1所示。

        顯然,有許多分類超平面可以將圖1中的“圓圈”和“五角星”兩類點正確區(qū)分開。

        取訓練樣本{},其中:為第i個樣本,。對于線性可分訓練集存在超平面L:,即式1.2和式1.3。

        , (1.2)

        , (1.3)

        定義超平面L1:,,設(shè)為超平面L1上的一點;超平面L2:,,為超平面L2上的一點。和滿足式1.4和式1.5。

        (1.4)

        (1.5)

        超平面L1與超平面L2之間間隔為式1.6。

        (1.6)

        尋找超平面L1和L2間隔最大化的平面——最優(yōu)超平面,即為二次規(guī)劃問題,如式1.7和式1.8所示。

        (1.7)

        s.t (1.8)

        利用lagrange乘子法解決上述二次規(guī)劃問題,建立lagrange函數(shù)如式1.9所示。

        (1.9)

        其中為lagrange乘子。

        將lagrange函數(shù)對,求其最小值,對求其最大值,解、和在函數(shù)的鞍點上滿足式1.10和式1.11。

        (1.10)

        (1.11)

        將式1.10代入lagrange函數(shù)式1.9,利用式1.11可得原優(yōu)化問題的對偶問題,可構(gòu)造出最優(yōu)超平面,如式1.12、式1.13和式1.14所示。

        (1.12)

        s.t. (1.13)

        (1.14)

        基于最優(yōu)超平面的分類規(guī)則如式1.15所示。

        (1.15)

        1.2 核函數(shù)

        對于二維空間的線性不可分數(shù)據(jù)群,為構(gòu)造最優(yōu)超平面,其函數(shù)表達式如式1.16所示。

        (1.16)

        構(gòu)造向量y和a,如式1.17和式1.18所示。

        (1.17)

        (1.18)

        可轉(zhuǎn)化為,即式1.19。

        (1.19)

        可見,映射到四維空間后,原來在二維空間中一個線性不可分的問題變成了線性可分,構(gòu)造函數(shù),其滿足式1.20和式1.21。

        (1.20)

        (1.21)

        其中,為低維空間多維常量,為低維空間多維變量,為由映射到高維空間的對應常量,為由映射到高維空間的對應向量。

        若只需將低維空間的輸入代入式1.20,計算即可得到高維空間的線性分類函數(shù),則函數(shù)稱為核函數(shù),能使得上述式1.20和式1.21中的兩個函數(shù)計算結(jié)果完全一致。

        2 燃氣管網(wǎng)日負荷的影響因素及數(shù)據(jù)處理

        本文以吉林市某小區(qū)燃氣管網(wǎng)日負荷為研究對象,時間為2012年9月份非節(jié)假日和節(jié)假日,非節(jié)假日的時間跨度為2012年9月份工作日,即9月3日至9月7日、9月10日至9月14日、9月17日至9月21日、9月24日至9月29日,節(jié)假日的時間跨度為2012年9月份中的部分法定節(jié)假日,取9月8日、9日、15日、16日、22日、23日、30日。將日最高溫度、日最低溫度、日平均溫度、小區(qū)人員最高年齡、小區(qū)人員最低年齡、小區(qū)人員平均年齡作為影響小區(qū)燃氣管網(wǎng)日負荷的主要因素,分別作為SVM的6個輸入量。將小區(qū)人員臨時出差、小區(qū)臨時增加暫住人口等隨機因素作為影響小區(qū)燃氣管網(wǎng)日負荷的次要因素,采用隨機變量法統(tǒng)一作為SVM的1個輸入量。將小區(qū)燃氣管網(wǎng)日負荷實際消耗量作為SVM的輸出量。由于溫度、人員數(shù)量、人員年齡等數(shù)據(jù)單位不同,會引起數(shù)據(jù)數(shù)量級差異,本文對日最高溫度、日最低溫度、日平均溫度、小區(qū)人員最高年齡、小區(qū)人員最低年齡、小區(qū)人員平均年齡、隨機因素、日負荷實際消耗量的數(shù)據(jù)進行[0,1]歸一化,使所測數(shù)據(jù)歸為一個數(shù)量級。工作日和節(jié)假日中小區(qū)燃氣管網(wǎng)日負荷影響因素的歸一化數(shù)據(jù),分別如圖2和圖3所示。endprint

        3 建立小區(qū)燃氣管網(wǎng)日負荷的預測模型

        3.1 支持向量機建模

        本文采用多項式核函數(shù)(polynomial核函數(shù)),徑向基核函數(shù)(radial basis function核函數(shù)),懲罰參數(shù)c和gamma值采用grid-search算法、遺傳算法(GA)、粒子群算法(PSO)進行參數(shù)尋優(yōu)。對以上組合所建立的不同支持向量模型進行交叉試驗,尋求最優(yōu)預測模型。利用該模型對工作日和節(jié)假日的小區(qū)燃氣管網(wǎng)日負荷進行預測。

        3.2 試驗結(jié)果分析

        尋求工作日小區(qū)燃氣管網(wǎng)日負荷最優(yōu)預測模型階段,訓練集時間跨度為9月3日至9月7日、9月10日至9月14日、9月17日至9月21日,預測集時間跨度為9月24日至9月29日;尋求節(jié)假日小區(qū)燃氣管網(wǎng)日負荷最優(yōu)預測模型階段,訓練集時間跨度為9月8日、9日、15日、16日、22日,預測集時間跨度為9月23日、30日。訓練集和預測集的數(shù)據(jù)為日最高溫度、日最低溫度、日平均溫度、小區(qū)人員最高年齡、小區(qū)人員最低年齡、小區(qū)人員平均年齡、隨機因素、日負荷實際值。利用支持向量機中polynomial核函數(shù)和radial basis function核函數(shù)(RBF)、優(yōu)化方法采用grid-search算法、遺傳算法(GA)、粒子群算法(PSO)。對工作日、節(jié)假日小區(qū)燃氣管網(wǎng)日負荷產(chǎn)生的預測效果分別見表1和表2。

        由表1可見,采用radial basis function核函數(shù)和PSO優(yōu)化算法的模型,可以實現(xiàn)工作日小區(qū)燃氣管網(wǎng)日負荷預測集的擬合程度最高,達到了90.781 %。預測結(jié)果如圖4、圖5所示。

        由表2可見,采用polynomial核函數(shù)和PSO優(yōu)化算法的模型,可以實現(xiàn)節(jié)假日小區(qū)燃氣管網(wǎng)日負荷預測集的擬合程度最高,達到了92.538 %。預測結(jié)果如圖6、圖7所示。

        可見,支持向量機采用radial basis function核函數(shù)和PSO優(yōu)化算法時,對工作日小區(qū)燃氣管網(wǎng)日負荷預測集的擬合程度最高;采用polynomial核函數(shù)和PSO優(yōu)化算法時,對節(jié)假日小區(qū)燃氣管網(wǎng)日負荷預測集的擬合程度最高,均高于90%。

        4 結(jié)語

        該文采用多項式核函數(shù)、徑向基核函數(shù),利用支持向量機對吉林市某小區(qū)工作日、節(jié)假日的燃氣管網(wǎng)日負荷建立預測模型,模型參數(shù)采用grid-search算法、遺傳算法、粒子群算法進行優(yōu)化,將日最高溫度、日最低溫度、日平均溫度、小區(qū)人員最高年齡、小區(qū)人員最低年齡、小區(qū)人員平均年齡作為影響小區(qū)燃氣管網(wǎng)日負荷的主要因素,將小區(qū)人員臨時出差、小區(qū)臨時增加暫住人口等隨機因素作為影響小區(qū)燃氣管網(wǎng)日負荷的次要因素,主要因素與次要因素統(tǒng)一作為SVM的輸入量,日負荷預測值作為SVM的輸出量。試驗結(jié)果表明,采用支持向量機建立的城市小區(qū)燃氣管網(wǎng)日負荷預測模型,其預測精度達到90%以上,為城市小區(qū)燃氣管網(wǎng)年負荷、月負荷的進一步分析提供了理論依據(jù)、技術(shù)支持和試驗方法。

        參考文獻

        [1] 譚羽非.城市燃氣管網(wǎng)日負荷預測的灰色神經(jīng)網(wǎng)絡模型[J].哈爾濱工業(yè)大學學報,2003,35(6):679-682.

        [2] 朱剛.城市燃氣管網(wǎng)負荷預測的研究[D].天津大學,2009.

        [3] Potocnik, Primoz,Govekar. Short-term natural gas consumption forecasting[J].Proceeding of the International Conference on Applied Simulation and Modelling,2007(5):353-357.

        [4] 杜元順.煤氣日負荷用的回歸分析方法[J].煤氣與熱力,1982(4):26-28.

        [5] 田一梅,趙元,趙新華.城市煤氣負荷預測[J].煤氣與熱力,1998,18(14):20-23.

        [6] 武海巍,于海業(yè),張蕾.光合有效輻射預測模型的核函數(shù)組合優(yōu)化[J].農(nóng)業(yè)機械學報,2011,42(6):167-173.

        [7] 武海巍.核函數(shù)與仿生智能算法在林下參光環(huán)境評價系統(tǒng)中的研究[D].吉林大學,2012.endprint

        3 建立小區(qū)燃氣管網(wǎng)日負荷的預測模型

        3.1 支持向量機建模

        本文采用多項式核函數(shù)(polynomial核函數(shù)),徑向基核函數(shù)(radial basis function核函數(shù)),懲罰參數(shù)c和gamma值采用grid-search算法、遺傳算法(GA)、粒子群算法(PSO)進行參數(shù)尋優(yōu)。對以上組合所建立的不同支持向量模型進行交叉試驗,尋求最優(yōu)預測模型。利用該模型對工作日和節(jié)假日的小區(qū)燃氣管網(wǎng)日負荷進行預測。

        3.2 試驗結(jié)果分析

        尋求工作日小區(qū)燃氣管網(wǎng)日負荷最優(yōu)預測模型階段,訓練集時間跨度為9月3日至9月7日、9月10日至9月14日、9月17日至9月21日,預測集時間跨度為9月24日至9月29日;尋求節(jié)假日小區(qū)燃氣管網(wǎng)日負荷最優(yōu)預測模型階段,訓練集時間跨度為9月8日、9日、15日、16日、22日,預測集時間跨度為9月23日、30日。訓練集和預測集的數(shù)據(jù)為日最高溫度、日最低溫度、日平均溫度、小區(qū)人員最高年齡、小區(qū)人員最低年齡、小區(qū)人員平均年齡、隨機因素、日負荷實際值。利用支持向量機中polynomial核函數(shù)和radial basis function核函數(shù)(RBF)、優(yōu)化方法采用grid-search算法、遺傳算法(GA)、粒子群算法(PSO)。對工作日、節(jié)假日小區(qū)燃氣管網(wǎng)日負荷產(chǎn)生的預測效果分別見表1和表2。

        由表1可見,采用radial basis function核函數(shù)和PSO優(yōu)化算法的模型,可以實現(xiàn)工作日小區(qū)燃氣管網(wǎng)日負荷預測集的擬合程度最高,達到了90.781 %。預測結(jié)果如圖4、圖5所示。

        由表2可見,采用polynomial核函數(shù)和PSO優(yōu)化算法的模型,可以實現(xiàn)節(jié)假日小區(qū)燃氣管網(wǎng)日負荷預測集的擬合程度最高,達到了92.538 %。預測結(jié)果如圖6、圖7所示。

        可見,支持向量機采用radial basis function核函數(shù)和PSO優(yōu)化算法時,對工作日小區(qū)燃氣管網(wǎng)日負荷預測集的擬合程度最高;采用polynomial核函數(shù)和PSO優(yōu)化算法時,對節(jié)假日小區(qū)燃氣管網(wǎng)日負荷預測集的擬合程度最高,均高于90%。

        4 結(jié)語

        該文采用多項式核函數(shù)、徑向基核函數(shù),利用支持向量機對吉林市某小區(qū)工作日、節(jié)假日的燃氣管網(wǎng)日負荷建立預測模型,模型參數(shù)采用grid-search算法、遺傳算法、粒子群算法進行優(yōu)化,將日最高溫度、日最低溫度、日平均溫度、小區(qū)人員最高年齡、小區(qū)人員最低年齡、小區(qū)人員平均年齡作為影響小區(qū)燃氣管網(wǎng)日負荷的主要因素,將小區(qū)人員臨時出差、小區(qū)臨時增加暫住人口等隨機因素作為影響小區(qū)燃氣管網(wǎng)日負荷的次要因素,主要因素與次要因素統(tǒng)一作為SVM的輸入量,日負荷預測值作為SVM的輸出量。試驗結(jié)果表明,采用支持向量機建立的城市小區(qū)燃氣管網(wǎng)日負荷預測模型,其預測精度達到90%以上,為城市小區(qū)燃氣管網(wǎng)年負荷、月負荷的進一步分析提供了理論依據(jù)、技術(shù)支持和試驗方法。

        參考文獻

        [1] 譚羽非.城市燃氣管網(wǎng)日負荷預測的灰色神經(jīng)網(wǎng)絡模型[J].哈爾濱工業(yè)大學學報,2003,35(6):679-682.

        [2] 朱剛.城市燃氣管網(wǎng)負荷預測的研究[D].天津大學,2009.

        [3] Potocnik, Primoz,Govekar. Short-term natural gas consumption forecasting[J].Proceeding of the International Conference on Applied Simulation and Modelling,2007(5):353-357.

        [4] 杜元順.煤氣日負荷用的回歸分析方法[J].煤氣與熱力,1982(4):26-28.

        [5] 田一梅,趙元,趙新華.城市煤氣負荷預測[J].煤氣與熱力,1998,18(14):20-23.

        [6] 武海巍,于海業(yè),張蕾.光合有效輻射預測模型的核函數(shù)組合優(yōu)化[J].農(nóng)業(yè)機械學報,2011,42(6):167-173.

        [7] 武海巍.核函數(shù)與仿生智能算法在林下參光環(huán)境評價系統(tǒng)中的研究[D].吉林大學,2012.endprint

        3 建立小區(qū)燃氣管網(wǎng)日負荷的預測模型

        3.1 支持向量機建模

        本文采用多項式核函數(shù)(polynomial核函數(shù)),徑向基核函數(shù)(radial basis function核函數(shù)),懲罰參數(shù)c和gamma值采用grid-search算法、遺傳算法(GA)、粒子群算法(PSO)進行參數(shù)尋優(yōu)。對以上組合所建立的不同支持向量模型進行交叉試驗,尋求最優(yōu)預測模型。利用該模型對工作日和節(jié)假日的小區(qū)燃氣管網(wǎng)日負荷進行預測。

        3.2 試驗結(jié)果分析

        尋求工作日小區(qū)燃氣管網(wǎng)日負荷最優(yōu)預測模型階段,訓練集時間跨度為9月3日至9月7日、9月10日至9月14日、9月17日至9月21日,預測集時間跨度為9月24日至9月29日;尋求節(jié)假日小區(qū)燃氣管網(wǎng)日負荷最優(yōu)預測模型階段,訓練集時間跨度為9月8日、9日、15日、16日、22日,預測集時間跨度為9月23日、30日。訓練集和預測集的數(shù)據(jù)為日最高溫度、日最低溫度、日平均溫度、小區(qū)人員最高年齡、小區(qū)人員最低年齡、小區(qū)人員平均年齡、隨機因素、日負荷實際值。利用支持向量機中polynomial核函數(shù)和radial basis function核函數(shù)(RBF)、優(yōu)化方法采用grid-search算法、遺傳算法(GA)、粒子群算法(PSO)。對工作日、節(jié)假日小區(qū)燃氣管網(wǎng)日負荷產(chǎn)生的預測效果分別見表1和表2。

        由表1可見,采用radial basis function核函數(shù)和PSO優(yōu)化算法的模型,可以實現(xiàn)工作日小區(qū)燃氣管網(wǎng)日負荷預測集的擬合程度最高,達到了90.781 %。預測結(jié)果如圖4、圖5所示。

        由表2可見,采用polynomial核函數(shù)和PSO優(yōu)化算法的模型,可以實現(xiàn)節(jié)假日小區(qū)燃氣管網(wǎng)日負荷預測集的擬合程度最高,達到了92.538 %。預測結(jié)果如圖6、圖7所示。

        可見,支持向量機采用radial basis function核函數(shù)和PSO優(yōu)化算法時,對工作日小區(qū)燃氣管網(wǎng)日負荷預測集的擬合程度最高;采用polynomial核函數(shù)和PSO優(yōu)化算法時,對節(jié)假日小區(qū)燃氣管網(wǎng)日負荷預測集的擬合程度最高,均高于90%。

        4 結(jié)語

        該文采用多項式核函數(shù)、徑向基核函數(shù),利用支持向量機對吉林市某小區(qū)工作日、節(jié)假日的燃氣管網(wǎng)日負荷建立預測模型,模型參數(shù)采用grid-search算法、遺傳算法、粒子群算法進行優(yōu)化,將日最高溫度、日最低溫度、日平均溫度、小區(qū)人員最高年齡、小區(qū)人員最低年齡、小區(qū)人員平均年齡作為影響小區(qū)燃氣管網(wǎng)日負荷的主要因素,將小區(qū)人員臨時出差、小區(qū)臨時增加暫住人口等隨機因素作為影響小區(qū)燃氣管網(wǎng)日負荷的次要因素,主要因素與次要因素統(tǒng)一作為SVM的輸入量,日負荷預測值作為SVM的輸出量。試驗結(jié)果表明,采用支持向量機建立的城市小區(qū)燃氣管網(wǎng)日負荷預測模型,其預測精度達到90%以上,為城市小區(qū)燃氣管網(wǎng)年負荷、月負荷的進一步分析提供了理論依據(jù)、技術(shù)支持和試驗方法。

        參考文獻

        [1] 譚羽非.城市燃氣管網(wǎng)日負荷預測的灰色神經(jīng)網(wǎng)絡模型[J].哈爾濱工業(yè)大學學報,2003,35(6):679-682.

        [2] 朱剛.城市燃氣管網(wǎng)負荷預測的研究[D].天津大學,2009.

        [3] Potocnik, Primoz,Govekar. Short-term natural gas consumption forecasting[J].Proceeding of the International Conference on Applied Simulation and Modelling,2007(5):353-357.

        [4] 杜元順.煤氣日負荷用的回歸分析方法[J].煤氣與熱力,1982(4):26-28.

        [5] 田一梅,趙元,趙新華.城市煤氣負荷預測[J].煤氣與熱力,1998,18(14):20-23.

        [6] 武海巍,于海業(yè),張蕾.光合有效輻射預測模型的核函數(shù)組合優(yōu)化[J].農(nóng)業(yè)機械學報,2011,42(6):167-173.

        [7] 武海巍.核函數(shù)與仿生智能算法在林下參光環(huán)境評價系統(tǒng)中的研究[D].吉林大學,2012.endprint

        猜你喜歡
        支持向量機
        基于支持向量回歸機的電能質(zhì)量評估
        基于智能優(yōu)化算法選擇特征的網(wǎng)絡入侵檢測
        數(shù)據(jù)挖掘技術(shù)在電廠經(jīng)濟性分析系統(tǒng)中的應用Q
        基于改進支持向量機的船舶縱搖預報模型
        中國水運(2016年11期)2017-01-04 12:26:47
        基于SVM的煙草銷售量預測
        軟件導刊(2016年11期)2016-12-22 21:52:38
        動態(tài)場景中的視覺目標識別方法分析
        論提高裝備故障預測準確度的方法途徑
        價值工程(2016年32期)2016-12-20 20:36:43
        基于熵技術(shù)的公共事業(yè)費最優(yōu)組合預測
        價值工程(2016年29期)2016-11-14 00:13:35
        基于支持向量機的金融數(shù)據(jù)分析研究
        管理類研究生支持向量機預測決策實驗教學研究
        考試周刊(2016年53期)2016-07-15 09:08:21
        国语对白做受xxxxx在线| 一本之道日本熟妇人妻| 一区二区三区在线观看高清视频| 国产91色在线|亚洲| 亚洲av成本人无码网站| 国产无码swag专区| 国产一级三级三级在线视| 久久精品国产6699国产精| 456亚洲人成在线播放网站| 亚洲香蕉毛片久久网站老妇人| 亚洲av永久无码精品成人| 网友自拍人妻一区二区三区三州| 一区二区三区国产精品| 亚洲国产精品美女久久久| 亚洲成av人片在久久性色av| 亚洲精品中文字幕一二三| 97精品熟女少妇一区二区三区 | 日本av一区二区在线| 人妻少妇av中文字幕乱码| 久久久久久人妻无码| 亚洲av无码一区东京热| 国产日产精品_国产精品毛片| 国产成人精品无码一区二区三区| 亚洲处破女av日韩精品| 欧美亚洲日韩国产人成在线播放| 国产精品一区二区久久| 亚洲电影久久久久久久9999| 中文字幕亚洲视频三区| 久久老熟女一区二区三区福利| 人人妻人人澡人人爽人人精品av | 午夜无码一区二区三区在线| 激情 一区二区| 国产人成在线免费视频| 午夜男女靠比视频免费| 国产av自拍视频在线观看| 亚洲av色欲色欲www| 欧美性xxxx极品高清| 在线精品国产一区二区三区| 处破痛哭a√18成年片免费| 国内久久婷婷精品人双人| 成人免费毛片在线播放|