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

    基于多特征的APSO-SVR的模擬電路故障預(yù)測

    2021-08-09 02:58:40龔振東
    關(guān)鍵詞:指示器波包特征向量

    王 力,龔振東

    (1.中國民航大學(xué) 職業(yè)技術(shù)學(xué)院,天津 300300;2.中國民航大學(xué) 電子信息與自動化學(xué)院,天津 300300)

    模擬電路是指用來對連續(xù)時(shí)間信號進(jìn)行傳輸、變換、處理和放大等工作的電路,被廣泛用于民航、醫(yī)療和軍事等領(lǐng)域的電子系統(tǒng)中,當(dāng)模擬電路發(fā)生故障時(shí)會導(dǎo)致整個(gè)系統(tǒng)無法正常運(yùn)行.為了提高模擬電路的可靠性,需要對模擬電路進(jìn)行故障預(yù)測與健康管理(Prognostics and Health Management,PHM)的研究[1].模擬電路的PHM 包括故障診斷與故障預(yù)測,在故障診斷領(lǐng)域的研究取得了較好的成果[2].

    模擬電路的故障預(yù)測旨在精準(zhǔn)地估計(jì)其剩余使用壽命(Remaining Useful Life,RUL),故障預(yù)測的主要方法為建立模型的方法與數(shù)據(jù)驅(qū)動的方法.相較而言,數(shù)據(jù)驅(qū)動的方法計(jì)算簡單且能獲得不錯(cuò)的效果[3],其中較典型的有支持向量機(jī)[4]、超限學(xué)習(xí)機(jī)[5]和相關(guān)向量機(jī)[6]等機(jī)器學(xué)習(xí)方法,均被成功地應(yīng)用到模擬電路的故障預(yù)測領(lǐng)域.文獻(xiàn)[7]提取模擬電路的頻率特征進(jìn)行故障預(yù)測.文獻(xiàn)[8]提取被測電路電壓信號并計(jì)算余弦距離進(jìn)行故障預(yù)測.文獻(xiàn)[9]提出了一種多核相關(guān)向量機(jī)方法對測試元件的故障發(fā)生時(shí)間開展預(yù)測.支持向量機(jī)回歸(Support Vector Regression,SVR)是當(dāng)前比較成熟的機(jī)器學(xué)習(xí)方法的理論,憑著實(shí)用性強(qiáng)、訓(xùn)練樣本小、泛化能力好等優(yōu)點(diǎn)[10-11],在故障預(yù)測方面優(yōu)勢突出,已被學(xué)者成功應(yīng)用于多種領(lǐng)域.文獻(xiàn)[12]通過BFOA 算法對SVR 模型進(jìn)行參數(shù)優(yōu)化,最后基于優(yōu)化模型開展溫度預(yù)測.文獻(xiàn)[13]用粒子群(Particle Swarm Optimization,PSO)算法優(yōu)化SVR對多聯(lián)機(jī)功耗進(jìn)行預(yù)測,取得不錯(cuò)的效果.文獻(xiàn)[14]指出SVR 中的懲罰因子C,核函數(shù)參數(shù) γ 對回歸預(yù)測精度有著很大的影響.SVR 通常使用網(wǎng)格搜索法來確定懲罰因子C與核函數(shù)參數(shù) γ,但是網(wǎng)格尋優(yōu)步長選取過大會導(dǎo)致模型預(yù)測精度低,步長過小會使優(yōu)化運(yùn)算時(shí)間過長,實(shí)際的使用效果并不好.

    本文擬將被測電路的輸出的統(tǒng)計(jì)特征和小波包分解能量特征作為特征向量,然后利用歐式距離量化電路退化的狀態(tài)建立故障指示器,最后用自適應(yīng)粒子群算法優(yōu)化的支持向量回歸預(yù)測模型進(jìn)行模擬電路的故障預(yù)測,估算電路的RUL.利用基準(zhǔn)電 路進(jìn)行仿真實(shí)驗(yàn)以驗(yàn)證方法.

    1 模擬電路的故障特征提取

    1.1 特征向量提取及預(yù)處理傳統(tǒng)的特征提取方法設(shè)置模擬電路的輸入為掃頻激勵(lì)源,只考慮統(tǒng)計(jì)特征里的單一響應(yīng)電壓作為輸入特征,無法充分反映電路的實(shí)時(shí)退化狀態(tài)[15].與傳統(tǒng)方法不同,小波包分解能對時(shí)間頻率進(jìn)行局部分析,通過對信號的平移伸縮實(shí)現(xiàn)低頻段的頻率細(xì)分和高頻段的時(shí)間細(xì)分[16],實(shí)現(xiàn)對信號的自適應(yīng)分析.因此,本文用多個(gè)統(tǒng)計(jì)特征和小波包分解后的能量特征融合的方式來構(gòu)造故障指示器.

    首先,選取均值、峰值、峭度、方差4 個(gè)常用統(tǒng)計(jì)特征構(gòu)成模擬電路的統(tǒng)計(jì)特征向量.均值a表示數(shù)據(jù)的平均趨勢;峰值p是電流、電壓或功率的最大瞬間值;峭度ku是歸一化的四階中心距,其大小跟電路的故障程度有關(guān),值越大偏離初始狀態(tài)越遠(yuǎn);方差u用來衡量數(shù)據(jù)和其均值之間偏離的程度.其公式分別定義如下:

    式中,X為模擬電路的輸出響應(yīng),E是取期望值,max 是取最大值.

    其次,利用小波分解對電路輸出數(shù)據(jù)進(jìn)行分解,用4 層小波包將信號分解成如圖1 所示的16 個(gè)子頻帶.電路輸出變化會影響各分量的能量值,這些能量值包含了模擬電路的狀態(tài)信息.因此計(jì)算各子頻帶的能量值用來構(gòu)建故障特征向量,其步驟如下:

    圖1 小波包分解圖Fig.1 Diagram of wavelet packet decomposition

    步驟 1對經(jīng)過4 層小波包分解的子信號能量進(jìn)行如下計(jì)算:

    式中,xjn為各分解信號的幅值.

    步驟 2計(jì)算分解后自信號的總能量

    步驟 3進(jìn)行能量歸一化處理

    步驟 4構(gòu)造模擬電路故障預(yù)測輸出能量的特征向量

    式中,h1到h15分別表示各分量的能量值.最后,將電路的統(tǒng)計(jì)特征和小波包的能量特征分別表示為G={g1,g2,g3,g4}和H={h0,h2,···,h15}.其中g(shù)1,g2,g3,g4分別表示均值、最大值、峭度和方差這4 個(gè)統(tǒng)計(jì)特征向量.將統(tǒng)計(jì)特征和小波包特征進(jìn)行融合,得 到特征向量F={G,H}.

    1.2 故障指示器生成故障指示器(Fault Indicator,F(xiàn)I)是反映電路性能退化的指標(biāo),其數(shù)值會隨著電路性能的衰退而升高.本文采用歐氏距離來衡量模擬電路元件退化前后特征向量間的相似度,構(gòu)建故障指示器以表征電路退化程度.

    2 APSO-SVR 的故障預(yù)測模型

    2.1 SVR 理論SVR 以支持向量機(jī)分類問題中的超平面決策邊界作為回歸模型[17].對于非線性數(shù)據(jù)集{xi,yi}(i=1,2,···,m),xi為樣本輸入,yi為樣本輸出,將樣本從低維空間映射到高維空間的公式如下:

    式中,ω 為法向量,b為位移項(xiàng),?(xi)是xi映射后的特征向量.

    訓(xùn)練支持向量回歸模型的過程就是尋找最優(yōu)ω,b使f(xi)盡量逼近yi,然后轉(zhuǎn)化成一個(gè)凸二次規(guī)劃的問題,其形式為

    最后引入滿足Mercer 條件的核函數(shù)K(x,xi)=?(xi),得到SVR 的回歸模型如下:

    其中,核函數(shù)K(x,xi)對應(yīng)的是x與xi映射到特征空間的內(nèi)積,可以代替特征空間內(nèi)的復(fù)雜運(yùn)算.

    核函數(shù)會對SVR 模型的性能產(chǎn)生較大影響,模擬電路的故障預(yù)測是多因子影響的非線性問題,采用徑向基核函數(shù)的SVR 可以更好地處理模擬電路故障預(yù)測數(shù)據(jù)的非線性問題且擁有良好的泛化能力.因此本文中的SVR 模型采用徑向基核函數(shù),其計(jì)算公式如下:

    式中,γ 是決定映射數(shù)據(jù)在特征空間中位置的核函數(shù)參數(shù),其值與支持向量數(shù)成反比.

    實(shí)際采用SVR 進(jìn)行模擬電路故障預(yù)測時(shí),如使用網(wǎng)格搜索法來確定SVR 的懲罰因子C與核函數(shù)參數(shù) γ,難以確定合理的參數(shù)搜索范圍與步長,對故障預(yù)測的準(zhǔn)確度影響很大,因此本文考慮使用APSO 算法來確定SVR 的懲罰因子與核函數(shù)參數(shù),以提升模型預(yù)測精度.

    2.2 APSO 算法PSO 算法屬于人工智能算法,它有收斂速度快、易實(shí)現(xiàn)、調(diào)整參數(shù)不多、全局搜索能力強(qiáng)等特點(diǎn)[18].但是,PSO 算法也有容易提前收斂、搜索精度不高、后期迭代效率低等不足.因此,本文使用APSO 優(yōu)化算法進(jìn)行參數(shù)尋優(yōu),以得到更高的擬合精度,并避免陷入局部最優(yōu).

    首先,初始化粒子群的種群參數(shù),并計(jì)算每個(gè)粒子的適應(yīng)度作為衡量標(biāo)準(zhǔn).每個(gè)粒子根據(jù)個(gè)體極值和全局極值來更新速度和位置其公式如下:

    式中,k為迭代次數(shù),θ 為用于調(diào)整模型尋優(yōu)性能的慣性因子,id表示第i個(gè)粒子的第d維,gd表示id找到的全局極值,學(xué)習(xí)因子c1、c2分別為調(diào)節(jié)單個(gè)粒子搜尋局部最優(yōu)位置方向和全局最優(yōu)位置方向的步長,r1和r2為服從 [0,1] 均勻分布的隨機(jī)數(shù).

    通過更新粒子當(dāng)前位置得到SVR 中待求參數(shù)C和γ的最優(yōu)解.在這個(gè)過程中,如果粒子的速度一直保持較大值,則在單次迭代中會導(dǎo)致粒子活動的范圍過大,跳出最優(yōu)解的區(qū)域.為使粒子的飛行速度隨著迭代次數(shù)的增加而減小,本文通過改變影響速度的參數(shù)c1,c2和慣性權(quán)重 θ,間接實(shí)現(xiàn)粒子飛行速度的調(diào)整.

    其中,設(shè)定c1從2.5 非線性減小到1.5,c2從2.1 非線性減小到1.1,計(jì)算公式如下:

    式中,n為當(dāng)前迭代次數(shù),nmax為最大迭代次數(shù).

    慣性權(quán)重 θ 是粒子群算法中影響調(diào)整全局與局部尋優(yōu)能力的關(guān)鍵參數(shù),其在較大時(shí)表現(xiàn)出較強(qiáng)的全局尋優(yōu)能力,較小時(shí)表現(xiàn)出較強(qiáng)局部尋優(yōu)能力.為滿足不同進(jìn)化階段的搜索需求,θ 應(yīng)隨著迭代次數(shù)的增加而逐步遞減.因此,本文引入自適應(yīng)慣性權(quán)重調(diào)整法,其公式為:式中,f為每個(gè)粒子的當(dāng)前適應(yīng)度,fmin為所有粒子 的最小適應(yīng)度,favg為所有粒子的平均適應(yīng)度.

    2.3 APSO-SVR 模型使用APSO 算法迭代尋優(yōu)得到SVR 的最優(yōu)參數(shù),得到模擬電路預(yù)測模型,具體步驟如下所示:

    步驟 1將FI 中的訓(xùn)練數(shù)據(jù)作為輸入;

    步驟 2確定初始粒子群的相關(guān)參數(shù),應(yīng)用APSO 算法對SVR 中的參數(shù)C和γ進(jìn)行尋優(yōu),其中參數(shù)C和γ的搜索范圍均為 [2-8,28];

    步驟 3粒子適應(yīng)度的計(jì)算與對比,以SVR的預(yù)測均方差函數(shù)作為目標(biāo)函數(shù),計(jì)算并比較各個(gè)粒子的適應(yīng)度值,找到全局最優(yōu)適應(yīng)度值;

    均方差計(jì)算公式如下:

    式中,η=1,2,···,L,L是用于訓(xùn)練的數(shù)據(jù)量,z(η) 表示真實(shí)數(shù)據(jù),y(η) 為預(yù)測數(shù)據(jù).

    步驟 4根據(jù)式(11)~(14)更新粒子;

    步驟 5滿足條件或達(dá)到迭代上限后輸出SVR最優(yōu)組合參數(shù)C和 γ.

    其模型構(gòu)建流程如圖2 所示.

    圖2 APSO-SVR 預(yù)測模型構(gòu)建流程Fig.2 Construction process of APSO-SVR prognostic model

    3 實(shí)驗(yàn)與分析

    為測試本文構(gòu)建的模擬電路故障預(yù)測模型性能,針對一個(gè)具體模擬電路進(jìn)行故障預(yù)測實(shí)驗(yàn).

    3.1 仿真實(shí)驗(yàn)以Sallen-Key 帶通濾波器電路為故障預(yù)測對象,測試本文構(gòu)建的APSO-SVR 模擬電路故障預(yù)測模型的性能.如圖3 所示,電阻和電容的容差分別設(shè)置為5%和10%,故障閾值設(shè)置為其參數(shù)值偏離標(biāo)稱值的50%.

    圖3 Sallen-key 帶通濾波器電路圖Fig.3 Circuit diagram of sallen-key bandpass filter

    選擇電阻R3 和電容C2 為實(shí)驗(yàn)對象,設(shè)定其在退化過程中的參數(shù)值在每一個(gè)時(shí)間節(jié)點(diǎn)變化一次,電路交流掃頻的起始頻率為1 Hz,終止頻率為100 kHz.電阻R3 從2 kΩ 逐步升高至3.2 kΩ,設(shè)置其以自身正常值的4%遞增,得到150 個(gè)時(shí)間節(jié)點(diǎn).電容C2 從5 nF 逐步升高至8 nF,設(shè)置其以自身正常值的4%遞增,得到150 個(gè)時(shí)間節(jié)點(diǎn).R3、C2 在參數(shù)值分別為3 kΩ、7.5 nF 時(shí)達(dá)到故障閾值,此時(shí)是第125 個(gè)時(shí)間節(jié)點(diǎn),判定電路元件即將發(fā)生故障.

    對樣本數(shù)據(jù)進(jìn)行劃分,將R3 和C2 的前80 個(gè)時(shí)間節(jié)點(diǎn)的樣本劃分為訓(xùn)練樣本,用于訓(xùn)練APSOSVR 故障預(yù)測模型,其余樣本劃分為測試樣本,用于測試模型的故障預(yù)測性能.

    根據(jù)Sallen-key 帶通濾波器電路的輸出數(shù)據(jù),提取統(tǒng)計(jì)特征和小波包能量特征.經(jīng)4 層小波包分解后得到如圖4 所示的子頻帶波形圖,然后計(jì)算其能量值,構(gòu)建小波包能量特征向量,并與統(tǒng)計(jì)特征進(jìn)行融合.

    圖4 子頻帶波形圖Fig.4 Waveform diagram of sub-band

    融合后的特征向量通過歐式距離公式計(jì)算得到故障指示器,其中,電阻R3 在各個(gè)時(shí)間節(jié)點(diǎn)的故障指示器如圖5所示,曲線的變化趨勢顯示使用歐氏距離得到的故障指示器可以較準(zhǔn)確地表征元件的性能退化軌跡.圖5 中電阻R3 在第125 個(gè)時(shí)間節(jié)點(diǎn)時(shí)電路發(fā)生故障,其故障閾值為0.864 45.用相同的方法計(jì)算能得到電容C2 的故障指示器,其故障閾值為0.868 11,當(dāng)故障指示器值大于或等于該閾值時(shí)即判斷模擬電路發(fā)生故障.

    圖5 R3 故障指示器Fig.5 The fault indicator of R3

    設(shè)置電容和電阻2 個(gè)仿真實(shí)驗(yàn):R3-80 和C2-80,即分別以R3 和C2 前80 個(gè)故障指示器數(shù)據(jù)作為訓(xùn)練樣本輸入到SVR 預(yù)測模型中進(jìn)行參數(shù)訓(xùn)練.初始化粒子群參數(shù),將APSO 的種群數(shù)量設(shè)定為20,最大迭代次數(shù)為100,進(jìn)行迭代尋優(yōu)找到S VR 的最優(yōu)組合參數(shù)C和γ,其參數(shù)如表1 所示.

    表1 SVR 組合參數(shù)Tab.1 Combination parameters of SVR

    3.2 預(yù)測結(jié)果與分析評價(jià)指標(biāo)采用剩余使用壽命預(yù)測值Res、絕對誤差EA和均方差EMS,其公式如下:

    式中,Res為預(yù)測剩余使用壽命,R為實(shí)際剩余使用壽命.

    采用APSO-SVR 模型對Sallen-Key 帶通濾波器電路中的電阻R3 和電容C2 進(jìn)行故障預(yù)測,其結(jié)果如圖6 和圖7 所示.其中,均方差分別為6.431×10-5和5.301×10-5,較小的均方差 表明2 個(gè)仿真實(shí)驗(yàn)的預(yù)測數(shù)據(jù)和真實(shí)數(shù)據(jù)之間的偏離度較小,能夠較為精確的預(yù)測出兩元件的退化軌跡.

    圖6 R3-80 的APSO 優(yōu)化后預(yù)測值和真實(shí)值Fig.6 Prediction and real value of R3-80 after APSO optimization

    圖7 C2-80 的APSO 優(yōu)化后預(yù)測值和真實(shí)值Fig.7 Prediction and real value of C2-80 after APSO optimization

    圖8、9 分別為未優(yōu)化的SVR 預(yù)測模型對電阻R3 和電容C2 的預(yù)測結(jié)果,其懲罰因子和核函數(shù)參數(shù)均采用網(wǎng)格搜索法進(jìn)行選擇.圖中預(yù)測值與真實(shí)值的偏差明顯大于APSO-SVR 預(yù)測模型,由此可知APSO-SVR 預(yù)測模型具有更高的預(yù)測精度.

    圖8 SVR 對R3-80 的預(yù)測值和真實(shí)值Fig.8 Prediction and real value of R3-80 by SVR

    對APSO-SVR 預(yù)測模型與未優(yōu)化的SVR 預(yù)測模型進(jìn)行故障預(yù)測得到的RUL 值進(jìn)行比較,如表2所示.APSO-SVR 預(yù)測模型對電阻R3 和電容C2進(jìn)行故障預(yù)測得到的RUL 值都更接近于真實(shí)值,其絕對誤差值分別為0 和1.未優(yōu)化的SVR 預(yù)測模型絕對誤差值分別為8 和9,由減少的絕對誤差值對比真實(shí)RUL 值能知優(yōu)化后的模型預(yù)測精度有較大提升.結(jié)果表明,APSO-SVR 預(yù)測模型能夠更為準(zhǔn)確地預(yù)測模擬電路的剩余使用壽命.

    圖9 SVR 對C2-80 的預(yù)測值和真實(shí)值Fig.9 Predicted and true values of C2-80 by SVR

    表2 模擬電路RUL 預(yù)測實(shí)驗(yàn)結(jié)果Tab.2 Results of analog circuit RUL prediction experimental

    4 結(jié)論

    本文提取模擬電路的輸出響應(yīng)的統(tǒng)計(jì)特征和小波包能量特征,并將兩個(gè)特征融合作為故障預(yù)測的特征向量,然后采用歐式距離構(gòu)建故障指示器來表征模擬電路元件的退化過程.最后利用APSO 算法代替?zhèn)鹘y(tǒng)網(wǎng)格搜索法對SVR 進(jìn)行參數(shù)尋優(yōu),建立APSO-SVR 預(yù)測模型并進(jìn)行模擬電路故障預(yù)測,估計(jì)元件的剩余使用壽命.經(jīng)實(shí)驗(yàn)結(jié)果表明,優(yōu)化后的算法能更加準(zhǔn)確地預(yù)測出模擬電路發(fā)生故障的時(shí)間及其剩余壽命時(shí)間,為模擬電路故障預(yù)測提供了新的方法.

    猜你喜歡
    指示器波包特征向量
    二年制職教本科線性代數(shù)課程的幾何化教學(xué)設(shè)計(jì)——以特征值和特征向量為例
    克羅內(nèi)克積的特征向量
    基于小波包Tsallis熵和RVM的模擬電路故障診斷
    一類特殊矩陣特征向量的求法
    EXCEL表格計(jì)算判斷矩陣近似特征向量在AHP法檢驗(yàn)上的應(yīng)用
    基于小波包變換的電力系統(tǒng)諧波分析
    柴油機(jī)儀表指示器的維修保養(yǎng)
    小波包理論與圖像小波包分解
    接地故障指示器的10kV線路接地故障的判斷與分析研究
    基于非接觸式電流互感器取電的故障指示器設(shè)計(jì)
    亚洲三级黄色毛片| 少妇裸体淫交视频免费看高清| 黄色视频,在线免费观看| 俄罗斯特黄特色一大片| 国产亚洲精品久久久com| 亚洲av电影不卡..在线观看| 欧美日本视频| 小说图片视频综合网站| 99热网站在线观看| 久久久精品94久久精品| 18禁在线无遮挡免费观看视频 | 亚洲国产高清在线一区二区三| 久久精品国产自在天天线| 美女xxoo啪啪120秒动态图| а√天堂www在线а√下载| 一级黄片播放器| 天堂网av新在线| 亚洲精品一卡2卡三卡4卡5卡| 两个人的视频大全免费| 最好的美女福利视频网| a级毛片a级免费在线| 日本a在线网址| 听说在线观看完整版免费高清| 99在线视频只有这里精品首页| 亚洲不卡免费看| 看免费成人av毛片| 国产 一区精品| 噜噜噜噜噜久久久久久91| 插阴视频在线观看视频| 成人美女网站在线观看视频| 国产爱豆传媒在线观看| 亚洲最大成人av| 成人av在线播放网站| 亚洲在线自拍视频| 1000部很黄的大片| 久久久久久久久大av| 校园人妻丝袜中文字幕| 精品日产1卡2卡| 国产精品亚洲美女久久久| 波多野结衣高清无吗| 久久久久久国产a免费观看| 精品一区二区三区视频在线| 国产伦精品一区二区三区四那| 日本黄大片高清| 99热只有精品国产| 麻豆成人午夜福利视频| 麻豆一二三区av精品| 91久久精品电影网| 淫秽高清视频在线观看| 男人狂女人下面高潮的视频| 成年版毛片免费区| 99热只有精品国产| 国产麻豆成人av免费视频| 欧美最黄视频在线播放免费| 97碰自拍视频| 一夜夜www| 欧美中文日本在线观看视频| 成年版毛片免费区| 久久这里只有精品中国| 最近手机中文字幕大全| 国产高清三级在线| 男女啪啪激烈高潮av片| 免费观看人在逋| 美女高潮的动态| 亚洲无线观看免费| 国内揄拍国产精品人妻在线| 18禁在线无遮挡免费观看视频 | 波多野结衣高清作品| 春色校园在线视频观看| 国产精品野战在线观看| 2021天堂中文幕一二区在线观| 噜噜噜噜噜久久久久久91| aaaaa片日本免费| 18+在线观看网站| 97碰自拍视频| 女人十人毛片免费观看3o分钟| 国产v大片淫在线免费观看| 又黄又爽又免费观看的视频| 丰满人妻一区二区三区视频av| 18禁裸乳无遮挡免费网站照片| 国产在视频线在精品| 精品不卡国产一区二区三区| 亚州av有码| 麻豆国产av国片精品| 亚洲丝袜综合中文字幕| 午夜福利18| 午夜亚洲福利在线播放| 欧美高清性xxxxhd video| 变态另类成人亚洲欧美熟女| 久久久成人免费电影| 午夜福利在线观看吧| 亚洲天堂国产精品一区在线| 午夜福利高清视频| 男女之事视频高清在线观看| 国产大屁股一区二区在线视频| 在线天堂最新版资源| 97超视频在线观看视频| 中文字幕精品亚洲无线码一区| 久久久欧美国产精品| 婷婷精品国产亚洲av在线| 免费看美女性在线毛片视频| 成年av动漫网址| 国产高清不卡午夜福利| 麻豆国产97在线/欧美| 午夜激情欧美在线| 不卡视频在线观看欧美| 免费av毛片视频| 男人的好看免费观看在线视频| 精品一区二区三区av网在线观看| 2021天堂中文幕一二区在线观| 亚洲七黄色美女视频| 蜜桃久久精品国产亚洲av| 插逼视频在线观看| 欧美性猛交黑人性爽| 日本黄色片子视频| 插逼视频在线观看| 床上黄色一级片| 午夜激情欧美在线| 波野结衣二区三区在线| 国产一区二区亚洲精品在线观看| 99热全是精品| av天堂中文字幕网| 97在线视频观看| 亚洲一区高清亚洲精品| 国产熟女欧美一区二区| or卡值多少钱| 亚洲第一电影网av| 欧美+亚洲+日韩+国产| 日本免费a在线| 亚洲不卡免费看| 国产极品精品免费视频能看的| 女生性感内裤真人,穿戴方法视频| aaaaa片日本免费| 日韩精品中文字幕看吧| 女生性感内裤真人,穿戴方法视频| 亚洲成a人片在线一区二区| 成人高潮视频无遮挡免费网站| 国产伦精品一区二区三区视频9| 高清毛片免费观看视频网站| 中文资源天堂在线| 国产私拍福利视频在线观看| 99在线视频只有这里精品首页| 国产亚洲精品av在线| 精华霜和精华液先用哪个| 久久久久国内视频| 成人精品一区二区免费| 99久久精品热视频| 麻豆乱淫一区二区| 一个人看的www免费观看视频| 亚洲一级一片aⅴ在线观看| 亚洲国产精品成人久久小说 | av中文乱码字幕在线| 成人午夜高清在线视频| 99精品在免费线老司机午夜| 欧美bdsm另类| 简卡轻食公司| 黄色日韩在线| 久久精品影院6| 99热这里只有是精品在线观看| 看十八女毛片水多多多| 久久久精品94久久精品| 午夜久久久久精精品| 两个人的视频大全免费| 欧美成人a在线观看| 国产亚洲精品久久久久久毛片| 天堂√8在线中文| 久久精品国产亚洲网站| 九九在线视频观看精品| 一个人免费在线观看电影| 又爽又黄a免费视频| 老熟妇仑乱视频hdxx| 综合色丁香网| 黄色欧美视频在线观看| 一区福利在线观看| 日本一二三区视频观看| 亚洲精华国产精华液的使用体验 | 成人午夜高清在线视频| a级一级毛片免费在线观看| 午夜福利18| 性色avwww在线观看| 国产在线男女| 99久久中文字幕三级久久日本| 最近2019中文字幕mv第一页| 精品福利观看| 欧美性猛交╳xxx乱大交人| 国产久久久一区二区三区| 免费在线观看影片大全网站| 91久久精品国产一区二区成人| 中文字幕久久专区| 国产一区二区在线观看日韩| 国产中年淑女户外野战色| 亚洲av二区三区四区| 少妇的逼好多水| 久久精品国产99精品国产亚洲性色| 欧美区成人在线视频| 亚洲精品国产成人久久av| 蜜桃久久精品国产亚洲av| 国产高清视频在线播放一区| 国产精品久久久久久av不卡| 99热精品在线国产| 香蕉av资源在线| 麻豆乱淫一区二区| av在线蜜桃| 观看美女的网站| 亚洲av成人精品一区久久| 国产av在哪里看| 欧美成人一区二区免费高清观看| 亚洲,欧美,日韩| 国产精品一及| 精品99又大又爽又粗少妇毛片| 国内精品美女久久久久久| 少妇裸体淫交视频免费看高清| 丝袜喷水一区| 村上凉子中文字幕在线| 俺也久久电影网| 中文资源天堂在线| 色哟哟·www| 亚洲人成网站在线播| 亚洲成人久久性| 日本爱情动作片www.在线观看 | 欧美性猛交╳xxx乱大交人| 国产在视频线在精品| 最好的美女福利视频网| 国产人妻一区二区三区在| 国产精品国产高清国产av| 少妇的逼好多水| 99热这里只有精品一区| 国产男靠女视频免费网站| 午夜老司机福利剧场| 欧美日本视频| 麻豆国产97在线/欧美| 欧美高清性xxxxhd video| 97碰自拍视频| 精品一区二区三区视频在线| 99视频精品全部免费 在线| 亚洲性夜色夜夜综合| 亚洲天堂国产精品一区在线| 欧美日韩精品成人综合77777| 精品一区二区三区视频在线| 成年版毛片免费区| 国产黄色小视频在线观看| 欧美一级a爱片免费观看看| www日本黄色视频网| 成人性生交大片免费视频hd| 欧美三级亚洲精品| 一边摸一边抽搐一进一小说| 国产毛片a区久久久久| 九九热线精品视视频播放| 天堂影院成人在线观看| 3wmmmm亚洲av在线观看| 久久精品夜色国产| 国产免费一级a男人的天堂| 99国产精品一区二区蜜桃av| 色5月婷婷丁香| 免费高清视频大片| 乱码一卡2卡4卡精品| 亚洲成人精品中文字幕电影| 亚洲国产色片| 免费观看在线日韩| 国模一区二区三区四区视频| 国产毛片a区久久久久| a级毛片a级免费在线| 观看免费一级毛片| 搡老岳熟女国产| 国产黄色视频一区二区在线观看 | 搞女人的毛片| 日韩强制内射视频| 国产欧美日韩精品一区二区| 天堂动漫精品| 级片在线观看| 精品一区二区免费观看| 欧美高清性xxxxhd video| 人人妻人人澡欧美一区二区| 别揉我奶头~嗯~啊~动态视频| 婷婷精品国产亚洲av| 听说在线观看完整版免费高清| 午夜精品一区二区三区免费看| 最新中文字幕久久久久| 97超级碰碰碰精品色视频在线观看| 成人亚洲欧美一区二区av| 我要搜黄色片| 亚洲成av人片在线播放无| 五月玫瑰六月丁香| 中文在线观看免费www的网站| 少妇人妻一区二区三区视频| 一级黄片播放器| 国产女主播在线喷水免费视频网站 | 国产 一区 欧美 日韩| 97人妻精品一区二区三区麻豆| 精品人妻熟女av久视频| 99久久久亚洲精品蜜臀av| 最后的刺客免费高清国语| 国产成人freesex在线 | 在线天堂最新版资源| 欧美在线一区亚洲| 国产精品乱码一区二三区的特点| 免费大片18禁| 亚洲精华国产精华液的使用体验 | 午夜爱爱视频在线播放| 亚洲欧美精品自产自拍| 日本一二三区视频观看| 亚洲中文字幕日韩| 一本久久中文字幕| 亚洲中文字幕日韩| 午夜免费激情av| 人妻制服诱惑在线中文字幕| 午夜亚洲福利在线播放| 久久久久久久久久成人| 观看美女的网站| 亚洲精华国产精华液的使用体验 | 人妻少妇偷人精品九色| 国产精品综合久久久久久久免费| 99久久精品国产国产毛片| 亚洲av五月六月丁香网| 国产成人a∨麻豆精品| 国产探花极品一区二区| 日韩欧美免费精品| 中国美白少妇内射xxxbb| 久久久久久久午夜电影| 亚洲中文日韩欧美视频| 久久精品综合一区二区三区| 国产激情偷乱视频一区二区| 少妇的逼好多水| 亚洲av成人精品一区久久| 丰满的人妻完整版| 久久久久久久久久久丰满| 婷婷色综合大香蕉| 欧美成人免费av一区二区三区| 日产精品乱码卡一卡2卡三| 欧美高清性xxxxhd video| 国产真实伦视频高清在线观看| 狂野欧美白嫩少妇大欣赏| 悠悠久久av| 国内精品宾馆在线| 1000部很黄的大片| 熟女人妻精品中文字幕| 91午夜精品亚洲一区二区三区| 1024手机看黄色片| 日本黄大片高清| 美女cb高潮喷水在线观看| 欧美日韩综合久久久久久| 日韩av在线大香蕉| 老熟妇乱子伦视频在线观看| 99久久无色码亚洲精品果冻| 亚洲专区国产一区二区| 插阴视频在线观看视频| 小说图片视频综合网站| 最近手机中文字幕大全| 欧美最新免费一区二区三区| 18禁裸乳无遮挡免费网站照片| av福利片在线观看| 国产成人影院久久av| 小说图片视频综合网站| 国产精华一区二区三区| 永久网站在线| 国产真实乱freesex| 国产精品一区二区性色av| 国产亚洲精品久久久com| 熟女电影av网| 久久国内精品自在自线图片| 亚洲成人av在线免费| 国产成人a∨麻豆精品| 大香蕉久久网| 18禁在线无遮挡免费观看视频 | 别揉我奶头~嗯~啊~动态视频| 三级国产精品欧美在线观看| 麻豆国产av国片精品| 两个人的视频大全免费| 国产黄色视频一区二区在线观看 | 中文亚洲av片在线观看爽| 18禁在线无遮挡免费观看视频 | 简卡轻食公司| 国产精品久久久久久亚洲av鲁大| 成年免费大片在线观看| eeuss影院久久| 午夜精品国产一区二区电影 | 日韩欧美三级三区| 欧美一级a爱片免费观看看| 亚洲精品色激情综合| 国产精品一区二区免费欧美| 国产高清视频在线观看网站| a级一级毛片免费在线观看| 国产精品日韩av在线免费观看| 可以在线观看毛片的网站| 国产精品三级大全| 亚洲欧美成人精品一区二区| 校园人妻丝袜中文字幕| 天天一区二区日本电影三级| 免费观看人在逋| 91在线观看av| 中文字幕熟女人妻在线| 人人妻人人看人人澡| 少妇的逼好多水| 国产精品野战在线观看| 俄罗斯特黄特色一大片| 又黄又爽又免费观看的视频| 国模一区二区三区四区视频| 寂寞人妻少妇视频99o| 亚洲最大成人手机在线| 人妻夜夜爽99麻豆av| 禁无遮挡网站| 12—13女人毛片做爰片一| 麻豆成人午夜福利视频| 人人妻,人人澡人人爽秒播| 国产精品一二三区在线看| 2021天堂中文幕一二区在线观| 男人狂女人下面高潮的视频| 搡女人真爽免费视频火全软件 | 午夜影院日韩av| 成人鲁丝片一二三区免费| av国产免费在线观看| 亚洲久久久久久中文字幕| 搡女人真爽免费视频火全软件 | www日本黄色视频网| 欧美高清性xxxxhd video| 国产精品久久久久久精品电影| 日韩欧美国产在线观看| 少妇的逼好多水| 亚洲中文日韩欧美视频| 精品免费久久久久久久清纯| .国产精品久久| 午夜福利视频1000在线观看| 一进一出抽搐gif免费好疼| 国产成人aa在线观看| av在线观看视频网站免费| 91精品国产九色| 久久这里只有精品中国| 欧美一级a爱片免费观看看| 色哟哟·www| 在线播放无遮挡| 一个人免费在线观看电影| 欧美不卡视频在线免费观看| 免费观看的影片在线观看| 精品一区二区三区人妻视频| 美女 人体艺术 gogo| 无遮挡黄片免费观看| 极品教师在线视频| 99热网站在线观看| 大又大粗又爽又黄少妇毛片口| a级毛片免费高清观看在线播放| 色5月婷婷丁香| 欧美zozozo另类| 免费av不卡在线播放| 男人的好看免费观看在线视频| 国产高清视频在线播放一区| 成人鲁丝片一二三区免费| 男女下面进入的视频免费午夜| 一本精品99久久精品77| 少妇的逼好多水| 国产精品不卡视频一区二区| 欧美在线一区亚洲| 寂寞人妻少妇视频99o| 午夜激情福利司机影院| 18禁黄网站禁片免费观看直播| 一本一本综合久久| 国产乱人视频| 女人被狂操c到高潮| 草草在线视频免费看| 国产 一区 欧美 日韩| 麻豆国产97在线/欧美| 亚洲无线观看免费| 国产精品久久视频播放| 成人av一区二区三区在线看| 成人二区视频| 91麻豆精品激情在线观看国产| 99精品在免费线老司机午夜| 国产黄色小视频在线观看| 变态另类成人亚洲欧美熟女| 熟女人妻精品中文字幕| 五月伊人婷婷丁香| 亚洲人成网站高清观看| 欧美性猛交╳xxx乱大交人| 看十八女毛片水多多多| 国产片特级美女逼逼视频| 亚洲精品一区av在线观看| 老司机福利观看| 国产精品国产三级国产av玫瑰| av专区在线播放| 国产精品女同一区二区软件| 亚洲成人av在线免费| 能在线免费观看的黄片| eeuss影院久久| 国产高清三级在线| 国产精品精品国产色婷婷| 国产黄片美女视频| 成人性生交大片免费视频hd| 美女免费视频网站| 精品一区二区三区av网在线观看| 美女xxoo啪啪120秒动态图| 国语自产精品视频在线第100页| 国内精品宾馆在线| 午夜福利18| av黄色大香蕉| 一级毛片电影观看 | 老司机福利观看| 99热只有精品国产| 久久久久免费精品人妻一区二区| 精品久久久久久久末码| 日本免费一区二区三区高清不卡| 老司机午夜福利在线观看视频| 丰满的人妻完整版| 亚洲人成网站在线观看播放| av在线播放精品| 国产成人精品久久久久久| 日本撒尿小便嘘嘘汇集6| 欧美bdsm另类| 成人一区二区视频在线观看| 成人午夜高清在线视频| 两个人的视频大全免费| 精品日产1卡2卡| 天堂动漫精品| 国产探花极品一区二区| 日本欧美国产在线视频| 精品人妻视频免费看| 日韩欧美国产在线观看| 舔av片在线| 白带黄色成豆腐渣| 精品人妻熟女av久视频| 精品国产三级普通话版| 免费一级毛片在线播放高清视频| 女人十人毛片免费观看3o分钟| 国内精品美女久久久久久| 成人二区视频| 亚洲熟妇中文字幕五十中出| 国产精品99久久久久久久久| 亚洲最大成人手机在线| 亚洲精品日韩av片在线观看| 黄色欧美视频在线观看| 日日干狠狠操夜夜爽| 精品一区二区三区人妻视频| 欧美成人精品欧美一级黄| 99热只有精品国产| 欧美激情国产日韩精品一区| 亚洲人与动物交配视频| 一区福利在线观看| 亚洲av中文字字幕乱码综合| 亚洲高清免费不卡视频| 91午夜精品亚洲一区二区三区| 九九爱精品视频在线观看| 国产一区二区亚洲精品在线观看| h日本视频在线播放| 亚洲人与动物交配视频| 久久婷婷人人爽人人干人人爱| 波多野结衣巨乳人妻| 一个人免费在线观看电影| 寂寞人妻少妇视频99o| 久久草成人影院| 亚洲av成人av| 成人高潮视频无遮挡免费网站| 国产片特级美女逼逼视频| 日韩中字成人| 欧美潮喷喷水| 亚洲五月天丁香| 91午夜精品亚洲一区二区三区| 色尼玛亚洲综合影院| 91久久精品国产一区二区三区| 人人妻人人澡人人爽人人夜夜 | 免费观看在线日韩| 少妇熟女欧美另类| 亚洲av.av天堂| 亚洲无线观看免费| 一本一本综合久久| 久久精品国产自在天天线| 日日摸夜夜添夜夜添小说| 欧美zozozo另类| 久久精品人妻少妇| 午夜久久久久精精品| 在线观看美女被高潮喷水网站| 国产一级毛片七仙女欲春2| 欧美成人免费av一区二区三区| 麻豆国产97在线/欧美| 黄色日韩在线| 国产精品一二三区在线看| 黑人高潮一二区| 国产成人freesex在线 | 国产一区二区在线观看日韩| 看片在线看免费视频| 亚洲五月天丁香| 免费看a级黄色片| 成人美女网站在线观看视频| av视频在线观看入口| 久久精品人妻少妇| 国产综合懂色| or卡值多少钱| 日本-黄色视频高清免费观看| 中出人妻视频一区二区| 亚洲国产欧美人成| 三级男女做爰猛烈吃奶摸视频| 少妇熟女欧美另类| 深夜精品福利| 亚洲av五月六月丁香网| 日本撒尿小便嘘嘘汇集6| 五月玫瑰六月丁香| 免费无遮挡裸体视频| 麻豆精品久久久久久蜜桃| 在线免费观看的www视频| 国产成人freesex在线 | 日韩精品有码人妻一区| 午夜日韩欧美国产| 久久九九热精品免费| 一区福利在线观看| 日韩成人伦理影院| 老司机影院成人| 一区福利在线观看| 尾随美女入室| 99在线视频只有这里精品首页| 久久婷婷人人爽人人干人人爱| 国产伦在线观看视频一区| 国产片特级美女逼逼视频| www.色视频.com| 国产伦一二天堂av在线观看| 简卡轻食公司| 99久久中文字幕三级久久日本| 看免费成人av毛片| 免费在线观看成人毛片|