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

    基于Copula函數(shù)的列車定位慣性單元剩余壽命預(yù)測方法

    2022-05-10 11:45:10上官偉巴明明孟月月
    鐵道學(xué)報(bào) 2022年4期
    關(guān)鍵詞:性能指標(biāo)加速度計(jì)慣性

    上官偉,巴明明,孟月月

    (1.北京交通大學(xué) 電子信息工程學(xué)院, 北京 100044;2.北京交通大學(xué) 軌道交通控制與安全國家重點(diǎn)實(shí)驗(yàn)室,北京 100044;3.北京市電磁兼容與衛(wèi)星導(dǎo)航工程技術(shù)研究中心, 北京 100044)

    高速鐵路的快速發(fā)展,對列車定位技術(shù)提出了更高的要求。基于全球衛(wèi)星導(dǎo)航系統(tǒng)(Global Navigation Satellite System, GNSS)和慣性測量單元(Inertial Measurement Unit, IMU)的列車組合定位方式,可以實(shí)現(xiàn)列車動態(tài)的精準(zhǔn)定位,保障列車定位的連續(xù)性、自主性,減少列車定位對軌道電路和應(yīng)答器等地面設(shè)備的依賴,降低鐵路建設(shè)投入以及運(yùn)營維護(hù)成本,是列車定位技術(shù)未來的重要發(fā)展方向[1]。作為其中的關(guān)鍵設(shè)備,確保組合定位系統(tǒng)中IMU在服役期間的工作性能,對保障系統(tǒng)的定位精度和運(yùn)行可靠性十分重要。列車定位慣性單元(以下簡稱“慣性單元”)是高集成度、高可靠性的捷聯(lián)式慣性導(dǎo)航系統(tǒng),主要由慣性傳感器和電子部件組成。捷聯(lián)式慣性導(dǎo)航系統(tǒng)直接固定在運(yùn)載體上,其中精密部件慣性傳感器直接承受高速列車的振動、沖擊以及角運(yùn)動。隨著運(yùn)行時(shí)間的積累,慣性單元的性能會發(fā)生退化,當(dāng)退化達(dá)到一定程度時(shí),測量系統(tǒng)不再滿足列車定位精度的要求,意味著慣性單元的有效使用壽命也隨之終結(jié)。

    針對上述問題,國內(nèi)外開展了大量關(guān)于慣性導(dǎo)航系統(tǒng)的壽命預(yù)測研究,主要通過設(shè)備的性能退化數(shù)據(jù),評估、預(yù)測設(shè)備服役期間的健康狀態(tài)或剩余使用壽命(Remaining Useful Life, RUL)。設(shè)備性能評價(jià)指標(biāo)是實(shí)現(xiàn)RUL精確預(yù)測的重要基礎(chǔ)?,F(xiàn)有的數(shù)據(jù)驅(qū)動RUL預(yù)測方法[2],從采用不同元數(shù)性能參量的角度可分為單性能參量和多元性能參量兩大類。

    單性能指標(biāo)的RUL預(yù)測方法極其豐富,根據(jù)目標(biāo)設(shè)備自身性能表征特點(diǎn)、性能數(shù)據(jù)采集的條件以及工程需求不同,進(jìn)行RUL預(yù)測的思路也會不同。Huang等[3]利用狀態(tài)空間刻畫一個帶自適應(yīng)漂移的非線性Wiener過程模型,推導(dǎo)了退化漂移狀態(tài)條件下RUL預(yù)測的顯示表達(dá)式。王瑞峰等[4]建立基于Gamma退化過程的狀態(tài)空間模型,采用EM-PF算法估計(jì)模型參數(shù),并據(jù)此求得RUL的概率密度和分布函數(shù)。Lei等[5]提出一種基于Wiener過程且考慮四源不確定性的隨機(jī)性能模型,通過離線數(shù)據(jù)與在線數(shù)據(jù)的聯(lián)合估計(jì),實(shí)時(shí)獲得軸承的剩余壽命。Wang等[6]進(jìn)一步利用一個基于多隱狀態(tài)帶非線性漂移的擴(kuò)散過程,進(jìn)行可充電電池RUL的預(yù)測。

    目前針對具有多性能評價(jià)指標(biāo)的復(fù)雜設(shè)備,進(jìn)行健康狀態(tài)評估、RUL預(yù)測的方法主要有兩種途徑: ①采用數(shù)據(jù)融合的技術(shù),先將多元性能參數(shù)信息融合后進(jìn)行單特征的提取,再采用單性能參量的RUL預(yù)測方法;②直接對現(xiàn)有性能退化模型或者RUL預(yù)測模型進(jìn)行擴(kuò)展,采用相關(guān)性描述技術(shù),準(zhǔn)確、有效建立多元性能參量間的關(guān)系,確保RUL預(yù)測的精度。金曉航等[7]針對軸承具有二元性能指標(biāo),提出一種基于二元維納過程的RUL預(yù)測方法,并利用AIC準(zhǔn)則選擇Copula函數(shù)進(jìn)行描述兩個性能指標(biāo)間的相關(guān)性。Hao等[8]針對LED燈的二元性能參量采用二元非線性擴(kuò)散過程構(gòu)建退化模型,并基于Copula函數(shù)描述二元性能參量的依賴關(guān)系,利用蒙特卡洛(MCMC)算法進(jìn)行參數(shù)估計(jì),實(shí)現(xiàn)退化設(shè)備的可靠性評估。Peng等[9]提出一種基于逆Gaussian的二元退化模型,采用Copula函數(shù)建立聯(lián)合分布,并利用兩階段貝葉斯方法進(jìn)行退化過程和Copula函數(shù)參數(shù)估計(jì)。

    針對慣性單元的RUL預(yù)測是對其進(jìn)行健康管理、維修維護(hù)時(shí),所需的重要健康狀態(tài)評估信息??v觀當(dāng)前PHM領(lǐng)域技術(shù)的研究重點(diǎn)可知,RUL預(yù)測的關(guān)鍵是如何考慮預(yù)測的不確定性因素,建立合理的性能退化模型,獲得RUL的概率密度分布并提高RUL預(yù)測精度。周源等[10]采用加速度計(jì)的零位電壓百分比增量作為加速度計(jì)的性能退化參數(shù),基于逆Gaussian過程建立性能退化模型,實(shí)現(xiàn)加速度計(jì)的剩余壽命預(yù)測。Miao等[11]針對慣性器件陀螺儀提出一種小波分析和SVM結(jié)合的方法,采用小波分析進(jìn)行陀螺儀性能數(shù)據(jù)的預(yù)處理,利用SVM對處理后的數(shù)據(jù)進(jìn)行建模,實(shí)現(xiàn)陀螺儀的剩余壽命預(yù)測。滕飛等[12]以加速度計(jì)加速退化數(shù)據(jù)為對象,提出了一種較為客觀、實(shí)用的加速退化數(shù)據(jù)統(tǒng)計(jì)分析方法,其特點(diǎn)是基于加速因子不變原則進(jìn)行加速退化數(shù)據(jù)建模,結(jié)合 Matlab軟件中的相關(guān)函數(shù)實(shí)現(xiàn)參數(shù)估計(jì)及可靠壽命預(yù)測。

    根據(jù)上述文獻(xiàn)研究現(xiàn)狀,提出一種列車定位慣性單元的RUL自適應(yīng)預(yù)測研究方法。首先,從IMU結(jié)構(gòu)組成和功能性質(zhì)角度出發(fā),選取慣性傳感器陀螺儀和加速度計(jì)的隨機(jī)誤差指標(biāo)評價(jià)設(shè)備性能狀態(tài),構(gòu)建IMU的性能評價(jià)指標(biāo)體系。其次,采用具有普適性的非線性漂移維納過程建立性能指標(biāo)退化模型?;谑走_(dá)時(shí)間概念獲取RUL的概率密度函數(shù)(Probability Density Function,PDF),在此基礎(chǔ)上,利用Copula理論描述二元性能指標(biāo)間的動態(tài)相關(guān)性,建立二元性能退化模型及RUL的聯(lián)合分布。再次,根據(jù)參數(shù)分步估計(jì)的思想,先對退化模型參數(shù)利用期望最大化(Expectation Maximization, EM)和Rauch-Tung-Striebel(RTS)后向平滑聯(lián)合算法,實(shí)現(xiàn)退化模型參數(shù)的實(shí)時(shí)自適應(yīng)估計(jì);再利用極大似然估計(jì)(Maximum Likelihood Estimation, MLE)算法更新Copula函數(shù)的相關(guān)性參數(shù)。最后,通過對比預(yù)測結(jié)果,證明本文提出的RUL預(yù)測方法可以有效減小預(yù)測的不確定性、提高預(yù)測精度。

    1 列車定位慣性單元的性能指標(biāo)選取

    作為列控系統(tǒng)的重要裝備,列車定位單元普遍使用基于GNSS和IMU的組合定位方式,其結(jié)構(gòu)見圖1。IMU作為其中的重要組成部分,對其開展基于性能狀態(tài)的壽命預(yù)測具有重要意義。慣性單元性能的主要評價(jià)標(biāo)準(zhǔn)是測量信息的精準(zhǔn)度。近年來,中國學(xué)者結(jié)合實(shí)際應(yīng)用,普遍采用隨機(jī)漂移中的隨機(jī)游走系數(shù)作為慣性器件的性能評價(jià)指標(biāo)[13]。

    圖1 基于GNSS和IMU組合的列車定位系統(tǒng)結(jié)構(gòu)

    本文的研究對象為實(shí)驗(yàn)室自主研發(fā)的列車定位系統(tǒng)中所采用的微機(jī)電型慣性測量單元(MEMS-IMU)。該型號慣性單元的產(chǎn)品標(biāo)準(zhǔn)規(guī)格說明中,明確標(biāo)明隨機(jī)游走系數(shù)是該型號慣性單元內(nèi)MEMS慣性器件的主要性能指標(biāo)。因此,本文針對列車定位慣性單元的性能指標(biāo)體系為加速度計(jì)的速率隨機(jī)游走系數(shù)和陀螺儀的角度隨機(jī)游走系數(shù)。

    本文用于研究驗(yàn)證的數(shù)據(jù)是對目標(biāo)慣性單元失效前近2年服役過程中周期性監(jiān)測得到的測量數(shù)據(jù)。根據(jù)實(shí)驗(yàn)需求和實(shí)際應(yīng)用安排,針對目標(biāo)慣性單元每7 d作為一個采樣間隔進(jìn)行性能監(jiān)測。在相同的采集條件下,性能監(jiān)測采集周期共計(jì)85周,總計(jì)85組性能退化數(shù)據(jù)。

    高速列車具有沿軌道水平運(yùn)動的特點(diǎn),根據(jù)列車定位系統(tǒng)車內(nèi)實(shí)驗(yàn)的安裝方式,可知慣性單元的主要敏感軸為x軸與y軸。

    針對目標(biāo)慣性單元退化周期內(nèi)的監(jiān)測數(shù)據(jù),采用Allan方差法[14]對敏感方向上慣性器件的測量數(shù)據(jù)進(jìn)行性能指標(biāo)辨識,所得每個監(jiān)測點(diǎn)的性能退化量見圖2。

    圖2 IMU關(guān)鍵部件性能監(jiān)測數(shù)據(jù)

    由圖2可知,沿敏感方向y軸的慣性傳感器,監(jiān)測獲得的性能指標(biāo)參數(shù)具有明顯的退化趨勢。因此,為有效評估設(shè)備的健康狀態(tài),列車定位慣性單元的性能評價(jià)指標(biāo)參數(shù)為y軸加速度計(jì)和陀螺儀的隨機(jī)游走系數(shù)。

    根據(jù)列車定位系統(tǒng)的精度指標(biāo)要求,結(jié)合實(shí)驗(yàn)室列車組合定位系統(tǒng)中目標(biāo)慣性設(shè)備的性能指標(biāo)評價(jià)體系,慣性單元性能的失效閾值設(shè)定情況如表1所示。

    表1 慣性單元性能指標(biāo)的失效閾值

    2 基于Copula函數(shù)的慣性單元剩余壽命建模方法

    2.1 單性能參量的剩余壽命建模

    2.1.1 單性能參量的性能退化模型

    由第1節(jié)可知,列車定位慣性單元的性能狀態(tài)可由陀螺儀和加速度計(jì)的性能指標(biāo)共同表征,且兩種慣性器件的退化監(jiān)測數(shù)據(jù)均呈現(xiàn)非單調(diào)、連續(xù)波動的特點(diǎn),表明單部件的性能退化過程具有明顯的隨機(jī)性和時(shí)變不確定性。對此本節(jié)將基于具有一般性帶有非線性漂移Wiener過程進(jìn)行單部件的性能退化建模。

    假設(shè)單部件的性能退化過程為{X(t),t≥0},基于帶非線性漂移Wiener過程建立的性能退化模型為

    ( 1 )

    性能退化軌跡的個體差異性見圖3。由圖3可以看出,實(shí)際退化過程中,設(shè)備性能狀態(tài)退化估計(jì),不僅和當(dāng)前時(shí)刻的性能退化量相關(guān),而且依賴于設(shè)備歷史性能退化軌跡。

    圖3 性能退化軌跡的個體差異性

    為了充分利用設(shè)備性能退化軌跡信息,假設(shè)當(dāng)前監(jiān)測時(shí)刻為tk,單部件的性能退化模型可以重構(gòu)為以下狀態(tài)空間模型

    ( 2 )

    依據(jù)式(2)和STF濾波技術(shù)的高斯特性可知,基于X0:k的漂移狀態(tài)λx的概率密度分布為高斯分布,即

    ( 3 )

    2.1.2 單性能參量的RUL分布

    基于性能退化過程首達(dá)時(shí)間的概念,認(rèn)為當(dāng)性能退化過程{X(t),t≥0}初次到達(dá)預(yù)先設(shè)定的性能指標(biāo)失效閾值ω時(shí),判定部件發(fā)生了失效。因此,性能退化模型{X(t),t≥0}對應(yīng)壽命的數(shù)學(xué)定義為

    T=inf{t:X(t)≥ω|X(0)<ω}

    ( 4 )

    根據(jù)式(4),文獻(xiàn)[16]給出退化過程基于首達(dá)時(shí)間的壽命概率密度函數(shù)如下

    ( 5 )

    Lk=inf{lk:X(lk+tk)≥ω|X0:k}

    ( 6 )

    由文獻(xiàn)[2]的相關(guān)推論可知,基于性能退化過程的狀態(tài)空間模型式( 2 )和部件運(yùn)行到當(dāng)前時(shí)刻的全部監(jiān)測數(shù)據(jù)X0:k,監(jiān)測時(shí)刻tk處RUL的概率密度分布fLk|X0:k(lk|X0:k)為

    ( 7 )

    綜上所述,每當(dāng)監(jiān)測獲取新的性能退化數(shù)據(jù)xk,可通過漂移狀態(tài)λk的更新機(jī)制,依賴全部歷史監(jiān)測數(shù)據(jù)X0:k,并利用式(7)自適應(yīng)的更新RUL的概率密度分布,實(shí)現(xiàn)RUL的實(shí)時(shí)預(yù)測。

    2.2 二元性能參量的相關(guān)性分析方法

    2.2.1 二元Copula函數(shù)原理

    Copula函數(shù)是把多個隨機(jī)向量X1,X2,…,XN的聯(lián)合分布函數(shù)F(x1,x2,…,xN)與各自對應(yīng)的邊緣分布函數(shù)u1,u2,…,uN相連接的函數(shù),表達(dá)式為

    F(x1,x2,…,xN)=C(u1,u2,…,uN;θ)

    ( 8 )

    式中:C(u1,u2,…,uN;θ)為Copula函數(shù)。其中,θ為相關(guān)性參數(shù);ui為第i(i=1,2,…,N)個隨機(jī)變量的邊緣分布。

    Sklar定理是用來證明Copula函數(shù)存在性的定理。其中,針對二元變量聯(lián)合分布的Sklar定理為:

    假設(shè)H(x,y)是邊緣分布分別為F(x)、G(y)的二元變量聯(lián)合分布函數(shù),則必定存在一個Copula函數(shù)C(u1,u2)滿足H(x,y)=C[F(x),G(y)]。若F(x),G(y)是連續(xù)函數(shù),那么C(u1,u2)唯一確定;相反,若F(x)、G(y)是一元分布函數(shù),C(u1,u2)是Copula函數(shù),則H(x,y)是具有邊緣分布為F(x)、G(y)的二元聯(lián)合分布函數(shù)。

    2.2.2 二元Copula函數(shù)選擇方法

    根據(jù)Sklar定理,在已知Copula函數(shù)具體表達(dá)形式的前提下,可以更加靈活地建立多元變量間的聯(lián)合分布。

    選擇對多元變量間相關(guān)性描述適配度較高的Copula函數(shù),有助于提高基于二元性能退化過程RUL預(yù)測的精度。本文基于性能指標(biāo)的邊緣分布和Coupla相關(guān)基礎(chǔ)理論,通過邊緣分布的二元直方分布特性及其與經(jīng)驗(yàn)Copula函數(shù)的歐氏距離,綜合評判選取最佳適配Copula函數(shù)來描述兩個性能變量間的動態(tài)相關(guān)性。

    定義經(jīng)驗(yàn)Copula:設(shè)(xi,yi)(i=1,2,…,n)為取自二維總體(X,Y)的樣本。將X,Y的經(jīng)驗(yàn)分布函數(shù)分別記為Fn(x)和Gn(y),得到樣本的經(jīng)驗(yàn)Copula函數(shù)定義如下

    (9)

    2.3 二元性能參量的剩余壽命建模

    2.3.1 二元性能參量退化過程

    假設(shè)基于二元性能指標(biāo)設(shè)備的性能退化過程Z(t)={X1(t),X2(t),t≥0},其中,Xi(t)(i=1,2)為2.1.1節(jié)中描述的單部件性能退化過程,并且滿足以下假設(shè):

    假設(shè)1:設(shè)備兩個性能指標(biāo)的退化過程均可利用帶線性或非線性漂移的Wiener過程進(jìn)行建模,如式(1)。

    假設(shè)2:針對同一設(shè)備,兩個性能指標(biāo)的監(jiān)測時(shí)間和間隔保持同步。

    假設(shè)3:兩個性能指標(biāo)的退化存在相關(guān)性,且相關(guān)性可以利用某種二元Copula函數(shù)進(jìn)行合理描述。

    假設(shè)4:當(dāng)設(shè)備的任意一個性能退化過程的退化量超過其對應(yīng)失效閾值時(shí),視作設(shè)備失效,設(shè)備的失效閾值向量記為W=[ωi],i=1,2。

    2.3.2 RUL聯(lián)合分布模型

    假設(shè)設(shè)備在當(dāng)前監(jiān)測時(shí)刻tk未發(fā)生失效,根據(jù)式(6)和二元性能退化過程失效定義(假設(shè)4),基于二元性能指標(biāo)設(shè)備的RUL定義如下

    (10)

    (11)

    (12)

    3 預(yù)測模型參數(shù)分步自適應(yīng)更新方法

    3.1 性能退化模型參數(shù)自適應(yīng)估計(jì)算法

    假設(shè)單性能指標(biāo)退化模型中“隱含”狀態(tài)λ的歷史數(shù)據(jù)為Υk={λ0,λ1,…,λk},EM算法的基本原理是通過最大化聯(lián)合似然函數(shù)p(X0:k,Υk|ψ)來逼近參數(shù)的極大似然估計(jì)。

    圖4為EM算法過程,L(θ)為完全數(shù)據(jù)的極大似然函數(shù),B(θ,θi)為極大似然函數(shù)的下界函數(shù),當(dāng)且僅當(dāng)θ=θi時(shí),兩函數(shù)值才相等。在E步,已有一個估計(jì)的參數(shù)值(初始化或上次迭代計(jì)算出的參數(shù)),但是此時(shí)這個參數(shù)值對應(yīng)B(θ,θi)的值不等于L(θ)的值。因此在E步以當(dāng)前的估計(jì)參數(shù),極大化B(θ,θi),使得該參數(shù)對應(yīng)的B(θ,θi)與L(θ)的數(shù)值相等。結(jié)果見圖4,θ=θi處極大似然函數(shù)與下界函數(shù)的函數(shù)值相等。在M步,求出使得B(θ,θi)取得極大的參數(shù)θi+1,然后將其作為下一次迭代的當(dāng)前參數(shù)值。

    圖4 EM算法過程圖

    EM算法主要步驟包括:

    E步-對似然函數(shù)求期望得

    (13)

    M步-求極值

    (14)

    最后,通過EM算法循環(huán)迭代可得到ψ的最終估計(jì)值。

    其中,根據(jù)單性能指標(biāo)退化過程的狀態(tài)空間模型式(2),E步-計(jì)算期望可以表示為

    (15)

    Step2后向迭代最優(yōu)平滑估計(jì)。

    form=k:1

    Step3協(xié)方差初始化。

    Step4平滑協(xié)方差的后向迭代計(jì)算。

    Zm|k=Pm|mJm-1+Jm(Zm+1|k-Pm|m)Jm-1

    (16)

    (17)

    (18)

    3.2 Copula函數(shù)未知參數(shù)估計(jì)

    (19)

    于是得到對數(shù)似然函數(shù)為

    (20)

    求解對數(shù)似然函數(shù)的最大值點(diǎn),即可得Copula函數(shù)的未知參數(shù)的估計(jì)為

    (21)

    未知參數(shù)向量Ψ自適應(yīng)的更新估計(jì)流程如下。

    Step1將模型參數(shù)Ψ分成兩部分:ψ退化模型參數(shù)和Copula函數(shù)的參數(shù)θ。

    Step3基于EM算法性能退化模型參數(shù)的自適應(yīng)預(yù)測:

    4 仿真測試與驗(yàn)證

    本節(jié)針對第2、3節(jié)提出的列車定位慣性單元的RUL自適應(yīng)預(yù)測方法,利用第1節(jié)中描述的慣性單元性能退化的實(shí)際監(jiān)測數(shù)據(jù),進(jìn)行實(shí)驗(yàn)分析和模型、算法的有效性驗(yàn)證。

    4.1 慣性單元性能監(jiān)測數(shù)據(jù)分析

    根據(jù)鐵路應(yīng)用RAMS標(biāo)準(zhǔn),列車定位慣性單元中慣性器件的失效閾值如表1所示,陀螺儀的性能退化軌跡于第64個監(jiān)測點(diǎn)達(dá)到失效閾值,加速度計(jì)的性能退化軌跡于第85個監(jiān)測點(diǎn)首次到達(dá)失效閾值。由設(shè)備失效定義推知,該實(shí)驗(yàn)設(shè)備在初始監(jiān)測點(diǎn)的真實(shí)RUL為64周。

    4.2 基于單性能指標(biāo)的剩余壽命預(yù)測

    陀螺儀和加速度計(jì)在每個監(jiān)測點(diǎn)處的實(shí)際退化量與一步預(yù)測結(jié)果的均方誤差(MSE)見圖5。預(yù)測的退化軌跡在前期由于數(shù)據(jù)不夠充足與實(shí)際退化軌跡略有出入。此外,陀螺儀和加速度計(jì)性能退化量預(yù)測MSE的期望分別為0.71×10-2與0.64×10-4。特別的在設(shè)備判定的失效監(jiān)測點(diǎn)處,陀螺儀性能退化量預(yù)測的MSE為0.21×10-2,加速度計(jì)性能退化量預(yù)測的MSE為0.47×10-4。表明本文采用的非線性漂移Wiener過程可以對不同類型的慣性器件的退化過程進(jìn)行有效建模。

    圖5 慣性器件性能退化量一步預(yù)測的MSE

    由式(17)和式(18)計(jì)算可以得到陀螺儀、加速度計(jì)性能退化模型參數(shù)隨著性能監(jiān)測數(shù)據(jù)的累積,在各監(jiān)測點(diǎn)自適應(yīng)估計(jì)、更新的軌跡見圖6、圖7。由圖6、圖7可見,模型參數(shù)可以敏感性能退化過程的狀態(tài)變化保持更新。與參考文獻(xiàn)[10]中使用EM算法對共軛先驗(yàn)分布函數(shù)超參數(shù)值進(jìn)行估計(jì)方法相比,本文提出的參數(shù)分步自適應(yīng)估計(jì)的方法不受先驗(yàn)初始化信息的影響,具有良好的魯棒性,為實(shí)現(xiàn)RUL的自適應(yīng)預(yù)測奠定基礎(chǔ)。

    圖6 陀螺儀性能退化模型參數(shù)更新曲線

    圖7 加速度計(jì)性能退化模型參數(shù)更新曲線

    已知各慣性器件的失效點(diǎn),由此可推出在每個監(jiān)測點(diǎn)處慣性器件的實(shí)際RUL。為進(jìn)一步定量評價(jià)RUL預(yù)測結(jié)果,定義實(shí)際RUL與RUL預(yù)測的期望均方誤差(EMSE)為

    (22)

    圖8 單性能指標(biāo)RUL自適應(yīng)預(yù)測結(jié)果

    圖8(a)、圖8(b)為慣性器件失效前多個關(guān)鍵監(jiān)測點(diǎn)處剩余壽命的PDF曲線,可看出預(yù)測的PDF曲線包含實(shí)際剩余壽命點(diǎn),并隨著監(jiān)測數(shù)據(jù)的累積曲線逐漸變窄,意味著RUL預(yù)測的不確定逐漸減小。

    圖9(a)、圖9(b)為慣性器件在各個監(jiān)測點(diǎn)剩余壽命預(yù)測的EMSE,在慣性器件各自的失效點(diǎn)處,陀螺儀RUL預(yù)測的EMSE為0.72,加速度計(jì)的EMSE為0.33。通過觀察EMSE曲線可以得出隨著數(shù)據(jù)累計(jì)RUL的預(yù)測精度將不斷提高。

    圖9 單性能指標(biāo)RUL的EMSE

    4.3 基于Copula函數(shù)的剩余壽命預(yù)測

    通過將設(shè)備失效前多個關(guān)鍵監(jiān)測點(diǎn)處的二元分布直方圖,與不同形式Copula函數(shù)的分布密度圖進(jìn)行對比,可以初步選定的Copula函數(shù)類型有Gaussian、Frank和Gumbel函數(shù)。接著,根據(jù)極大似然估計(jì)的方法對相應(yīng)Copula函數(shù)的未知參數(shù)進(jìn)行估計(jì),并據(jù)此求解歐氏距離,依據(jù)距離越小相關(guān)性描述越精準(zhǔn)的原則,選定Gumbel函數(shù)作為描述二元變量相關(guān)性的Copula函數(shù)。

    完成Copula函數(shù)選擇和基于式(21)的參數(shù)極大似然估計(jì)后,根據(jù)式(12)計(jì)算設(shè)備的實(shí)時(shí)剩余壽命概率密度函數(shù),得到設(shè)備RUL概率密度曲線見圖10。

    圖10 RUL概率密度曲線

    藍(lán)色實(shí)線是陀螺儀的PDF預(yù)測曲線,綠色點(diǎn)劃線是加速度計(jì)的PDF預(yù)測曲線,紅色虛線為基于陀螺儀和加速度計(jì)二元性能參量的PDF預(yù)測曲線。相較于另外兩條PDF曲線,二元性能參量的預(yù)測曲線方差較小、收斂更加集中、尖銳,具有較高的峰值。

    在失效前的最后6個監(jiān)測時(shí)刻,分別采用陀螺儀性能指標(biāo)、加速度計(jì)性能指標(biāo)和二元性能參量3種不同性能退化指標(biāo),計(jì)算RUL預(yù)測期望值,并計(jì)算預(yù)測的RUL期望值與實(shí)際RUL的均方相對誤差(MRE)為

    (23)

    表5 基于3種性能指標(biāo)方案的RUL預(yù)測的期望和MRE

    由表5可知,在失效前最后6個監(jiān)測時(shí)刻,考慮二者相關(guān)性RUL預(yù)測期望值的平均MRE為9.87%,基于陀螺儀RUL預(yù)測的平均MRE為13.14%??擅黠@看出本文提出的基于Copula函數(shù)的二元性能參量剩余壽命預(yù)測方法可以充分利用設(shè)備組成部件的性能退化信息,提高RUL預(yù)測精度。

    5 結(jié)論與展望

    本文主要針對列車定位慣性單元RUL預(yù)測方法進(jìn)行研究,基于列車定位慣性單元的性能評價(jià)指標(biāo)體系,構(gòu)建了基于Copula函數(shù)的二元性能參量剩余壽命聯(lián)合分布模型,并采用未知參數(shù)自適應(yīng)分步更新的方法,實(shí)現(xiàn)RUL的自適應(yīng)更新預(yù)測。結(jié)果表明,本文提出的基于二元性能參量的RUL預(yù)測方法可以提高預(yù)測精度,對存在多個退化變量且它們關(guān)系未知的設(shè)備的剩余壽命預(yù)測有著重要的參考價(jià)值。

    后續(xù)工作中在列車定位慣性單元具有較多退化、失效數(shù)據(jù)時(shí),可以結(jié)合機(jī)器學(xué)習(xí)與隨機(jī)過程方法預(yù)測設(shè)備RUL,進(jìn)一步豐富RUL預(yù)測結(jié)果,提高RUL的預(yù)測性能。

    猜你喜歡
    性能指標(biāo)加速度計(jì)慣性
    你真的了解慣性嗎
    沖破『慣性』 看慣性
    基于加速度計(jì)的起重機(jī)制動下滑量測量
    瀝青膠結(jié)料基本高溫性能指標(biāo)相關(guān)性研究
    石油瀝青(2021年1期)2021-04-13 01:31:08
    無處不在的慣性
    基于遺傳算法的加速度計(jì)免轉(zhuǎn)臺標(biāo)定方法
    常溫下硅微諧振加速度計(jì)零偏穩(wěn)定性的提高
    普遍存在的慣性
    儲熱水箱分層性能指標(biāo)的研究進(jìn)展
    WebGIS關(guān)鍵性能指標(biāo)測試技術(shù)研究
    综合色av麻豆| 色噜噜av男人的天堂激情| 日韩欧美在线二视频| 宅男免费午夜| 激情在线观看视频在线高清| 最近中文字幕高清免费大全6 | 网址你懂的国产日韩在线| 18禁在线播放成人免费| 三级国产精品欧美在线观看| 亚洲国产色片| 黄片小视频在线播放| av天堂中文字幕网| av黄色大香蕉| 久久久久久大精品| 久久久久久九九精品二区国产| 三级毛片av免费| 成年女人永久免费观看视频| 日韩欧美 国产精品| 日本一本二区三区精品| 成人性生交大片免费视频hd| or卡值多少钱| 一级a爱片免费观看的视频| 日本熟妇午夜| 女人被狂操c到高潮| 亚洲av熟女| 又粗又爽又猛毛片免费看| 日日摸夜夜添夜夜添av毛片 | 偷拍熟女少妇极品色| 久久国产乱子伦精品免费另类| 国产精品美女特级片免费视频播放器| 在线免费观看的www视频| 久久热精品热| 韩国av一区二区三区四区| 亚洲av电影不卡..在线观看| 我要搜黄色片| 老熟妇仑乱视频hdxx| 国产成人啪精品午夜网站| 我要搜黄色片| 亚洲片人在线观看| 国产精品一区二区三区四区免费观看 | 日韩欧美免费精品| 丁香欧美五月| 天天一区二区日本电影三级| 欧美3d第一页| 免费看美女性在线毛片视频| 久久精品人妻少妇| 欧美最新免费一区二区三区 | 99久久精品国产亚洲精品| 亚洲欧美清纯卡通| 狠狠狠狠99中文字幕| 国产欧美日韩一区二区三| 久久这里只有精品中国| 亚洲国产精品999在线| 中文字幕av在线有码专区| 男人和女人高潮做爰伦理| 美女高潮喷水抽搐中文字幕| 一区二区三区激情视频| 欧美潮喷喷水| 可以在线观看的亚洲视频| 欧美zozozo另类| 亚洲内射少妇av| 亚洲av美国av| 色在线成人网| 国产单亲对白刺激| 99久久精品国产亚洲精品| 三级毛片av免费| 精品熟女少妇八av免费久了| 国产精品久久视频播放| 69人妻影院| 两个人视频免费观看高清| 91麻豆精品激情在线观看国产| 人人妻人人澡欧美一区二区| 亚洲成人精品中文字幕电影| 欧美精品啪啪一区二区三区| 天堂av国产一区二区熟女人妻| 亚洲一区高清亚洲精品| 亚洲av熟女| 99在线人妻在线中文字幕| 成人国产一区最新在线观看| 在线播放无遮挡| 免费在线观看日本一区| 国产成人av教育| 51午夜福利影视在线观看| 老鸭窝网址在线观看| 天堂√8在线中文| 精品不卡国产一区二区三区| 十八禁国产超污无遮挡网站| 啪啪无遮挡十八禁网站| 亚洲一区二区三区色噜噜| 看片在线看免费视频| 中国美女看黄片| 久久中文看片网| 在线a可以看的网站| 人妻丰满熟妇av一区二区三区| 男人和女人高潮做爰伦理| 亚洲av免费在线观看| 久久香蕉精品热| 日日摸夜夜添夜夜添小说| 97超级碰碰碰精品色视频在线观看| 有码 亚洲区| 中文字幕av成人在线电影| 日韩欧美在线乱码| 午夜激情欧美在线| 美女高潮喷水抽搐中文字幕| 亚洲中文字幕一区二区三区有码在线看| 精品99又大又爽又粗少妇毛片 | 在线看三级毛片| 亚洲人成电影免费在线| 99国产精品一区二区蜜桃av| 最近视频中文字幕2019在线8| 深爱激情五月婷婷| 淫秽高清视频在线观看| 国产欧美日韩精品亚洲av| 久久精品影院6| 欧美日韩中文字幕国产精品一区二区三区| 露出奶头的视频| 九色成人免费人妻av| 美女黄网站色视频| 国产av在哪里看| 午夜福利在线在线| 国产白丝娇喘喷水9色精品| 国产精品人妻久久久久久| 成人性生交大片免费视频hd| 俄罗斯特黄特色一大片| 在线观看av片永久免费下载| 熟妇人妻久久中文字幕3abv| 乱码一卡2卡4卡精品| 一级黄片播放器| 一区二区三区四区激情视频 | av女优亚洲男人天堂| 天堂√8在线中文| 亚洲精品成人久久久久久| 两个人的视频大全免费| 中文资源天堂在线| 18+在线观看网站| 精品国内亚洲2022精品成人| 欧美色欧美亚洲另类二区| а√天堂www在线а√下载| 国产精品免费一区二区三区在线| 天堂网av新在线| 亚洲成人久久爱视频| 国产一区二区三区在线臀色熟女| 欧美bdsm另类| 亚洲 国产 在线| 91午夜精品亚洲一区二区三区 | 91av网一区二区| 99国产精品一区二区三区| 久久国产精品人妻蜜桃| 97超级碰碰碰精品色视频在线观看| 精品久久久久久久末码| 亚洲内射少妇av| 日本与韩国留学比较| 成年女人看的毛片在线观看| 成人永久免费在线观看视频| 高清在线国产一区| 麻豆成人午夜福利视频| 欧美zozozo另类| 成人特级av手机在线观看| 亚洲片人在线观看| 无遮挡黄片免费观看| 国产精品一区二区免费欧美| 99在线视频只有这里精品首页| 波多野结衣高清作品| 啪啪无遮挡十八禁网站| 香蕉av资源在线| 日本与韩国留学比较| 熟女电影av网| 一个人看视频在线观看www免费| 国产综合懂色| 国产成人欧美在线观看| 给我免费播放毛片高清在线观看| 国产一区二区亚洲精品在线观看| 亚洲熟妇中文字幕五十中出| 男人舔女人下体高潮全视频| av在线观看视频网站免费| 全区人妻精品视频| 91狼人影院| 天堂影院成人在线观看| 国产美女午夜福利| 极品教师在线免费播放| 国产免费av片在线观看野外av| 久久热精品热| 在线观看免费视频日本深夜| 国产熟女xx| 国产久久久一区二区三区| 亚洲精华国产精华精| 国产欧美日韩精品亚洲av| 精华霜和精华液先用哪个| 精品不卡国产一区二区三区| 看黄色毛片网站| 怎么达到女性高潮| av黄色大香蕉| 亚洲avbb在线观看| 国产69精品久久久久777片| 男女之事视频高清在线观看| 波多野结衣高清作品| 亚洲乱码一区二区免费版| 成人毛片a级毛片在线播放| 亚洲中文字幕日韩| 白带黄色成豆腐渣| 亚洲欧美日韩东京热| 可以在线观看毛片的网站| 精品久久久久久久久久久久久| 亚洲成av人片免费观看| 久久午夜福利片| 99riav亚洲国产免费| 久久久国产成人免费| 欧美日韩福利视频一区二区| 天天躁日日操中文字幕| 久久午夜福利片| 久99久视频精品免费| 性欧美人与动物交配| 日韩欧美在线二视频| 99热精品在线国产| 久久久久免费精品人妻一区二区| 日本成人三级电影网站| 欧美一区二区精品小视频在线| 无遮挡黄片免费观看| 黄片小视频在线播放| 九九热线精品视视频播放| 日韩av在线大香蕉| 欧美黑人巨大hd| 麻豆成人午夜福利视频| 国产精品亚洲av一区麻豆| 午夜久久久久精精品| 人人妻人人看人人澡| 黄色日韩在线| 亚洲中文字幕一区二区三区有码在线看| 国产三级中文精品| 国产精品av视频在线免费观看| 色吧在线观看| 丝袜美腿在线中文| 中文字幕精品亚洲无线码一区| 久久草成人影院| 又黄又爽又刺激的免费视频.| 婷婷精品国产亚洲av| 久久精品国产自在天天线| 男女床上黄色一级片免费看| 在线观看美女被高潮喷水网站 | 好男人在线观看高清免费视频| 特大巨黑吊av在线直播| 99热精品在线国产| 动漫黄色视频在线观看| 国产精品一及| 国产精品自产拍在线观看55亚洲| 国产高清有码在线观看视频| 亚洲国产色片| 99热这里只有精品一区| 免费在线观看日本一区| 日韩国内少妇激情av| 女同久久另类99精品国产91| 亚洲成人久久性| 男人舔女人下体高潮全视频| 国产高清视频在线观看网站| 女人十人毛片免费观看3o分钟| 成人av在线播放网站| 成人亚洲精品av一区二区| 99热精品在线国产| 两性午夜刺激爽爽歪歪视频在线观看| 国产精品1区2区在线观看.| 在线观看av片永久免费下载| 男人舔女人下体高潮全视频| 国产一区二区亚洲精品在线观看| 亚洲成人免费电影在线观看| 久久久久亚洲av毛片大全| 宅男免费午夜| 俄罗斯特黄特色一大片| 男插女下体视频免费在线播放| 欧美3d第一页| 欧美一级a爱片免费观看看| 国产一区二区三区视频了| 宅男免费午夜| 亚洲精华国产精华精| 久久草成人影院| 国产亚洲精品综合一区在线观看| 中文亚洲av片在线观看爽| 18禁黄网站禁片免费观看直播| 成年女人永久免费观看视频| 国产成人影院久久av| 蜜桃久久精品国产亚洲av| 18+在线观看网站| 国内精品美女久久久久久| 亚洲成人中文字幕在线播放| 我要看日韩黄色一级片| 真人一进一出gif抽搐免费| 国产成人欧美在线观看| 日日摸夜夜添夜夜添小说| 亚洲专区国产一区二区| av中文乱码字幕在线| 色哟哟·www| 麻豆成人午夜福利视频| 老司机午夜福利在线观看视频| 高潮久久久久久久久久久不卡| 999久久久精品免费观看国产| 久久人人爽人人爽人人片va | 成年人黄色毛片网站| 直男gayav资源| 欧美精品啪啪一区二区三区| 亚洲av成人不卡在线观看播放网| 香蕉av资源在线| 岛国在线免费视频观看| 中文字幕精品亚洲无线码一区| 天堂网av新在线| 欧美最新免费一区二区三区 | 国产一区二区亚洲精品在线观看| 日韩欧美精品免费久久 | 久久香蕉精品热| 亚洲国产精品sss在线观看| 欧美一区二区国产精品久久精品| 嫁个100分男人电影在线观看| 老司机福利观看| 亚洲中文日韩欧美视频| 国产极品精品免费视频能看的| 久久亚洲精品不卡| 国产亚洲欧美在线一区二区| 91在线精品国自产拍蜜月| 精品无人区乱码1区二区| 啦啦啦观看免费观看视频高清| 国产免费男女视频| 国产乱人视频| 十八禁网站免费在线| 18美女黄网站色大片免费观看| 国产日本99.免费观看| 18+在线观看网站| 嫩草影视91久久| 淫秽高清视频在线观看| 很黄的视频免费| 丁香六月欧美| 色哟哟哟哟哟哟| 国产精品免费一区二区三区在线| 97热精品久久久久久| 精品99又大又爽又粗少妇毛片 | 免费观看精品视频网站| 亚洲国产高清在线一区二区三| 久久人人爽人人爽人人片va | 亚洲人成网站高清观看| 亚洲成人久久性| 女生性感内裤真人,穿戴方法视频| 性色avwww在线观看| 国内久久婷婷六月综合欲色啪| 伦理电影大哥的女人| 亚洲国产精品sss在线观看| 免费大片18禁| 国产69精品久久久久777片| 色综合婷婷激情| 91久久精品电影网| 精品久久国产蜜桃| 精品熟女少妇八av免费久了| 日日夜夜操网爽| 亚洲av成人不卡在线观看播放网| 日本a在线网址| 亚洲不卡免费看| 免费无遮挡裸体视频| 亚洲激情在线av| 最后的刺客免费高清国语| 老女人水多毛片| 身体一侧抽搐| 嫩草影院新地址| 中文字幕免费在线视频6| 757午夜福利合集在线观看| 日韩欧美精品v在线| 超碰av人人做人人爽久久| 一级黄片播放器| 亚洲三级黄色毛片| 午夜福利高清视频| 久久久久久国产a免费观看| 可以在线观看的亚洲视频| 看片在线看免费视频| 日本黄大片高清| 国产国拍精品亚洲av在线观看| 最近在线观看免费完整版| 99久久成人亚洲精品观看| 亚洲av二区三区四区| 一a级毛片在线观看| 青草久久国产| 波多野结衣巨乳人妻| 国产精品影院久久| 免费高清视频大片| av女优亚洲男人天堂| 天堂av国产一区二区熟女人妻| 亚洲精品在线美女| 99热只有精品国产| 国产欧美日韩一区二区精品| 成人精品一区二区免费| av在线观看视频网站免费| 午夜激情福利司机影院| 亚洲国产高清在线一区二区三| 欧美黄色片欧美黄色片| 我要搜黄色片| 色综合欧美亚洲国产小说| 亚洲中文字幕一区二区三区有码在线看| 麻豆国产97在线/欧美| 精品熟女少妇八av免费久了| 啦啦啦观看免费观看视频高清| 亚洲精品日韩av片在线观看| 69人妻影院| 免费看日本二区| 我的老师免费观看完整版| 日本a在线网址| 脱女人内裤的视频| av天堂在线播放| 国产高清有码在线观看视频| 日本在线视频免费播放| 亚洲 欧美 日韩 在线 免费| 99在线人妻在线中文字幕| 国产成人av教育| 国产午夜精品论理片| 3wmmmm亚洲av在线观看| 毛片一级片免费看久久久久 | 国产国拍精品亚洲av在线观看| 日韩有码中文字幕| 国产精品人妻久久久久久| 亚洲成人中文字幕在线播放| 夜夜爽天天搞| 日韩欧美国产在线观看| 国内精品美女久久久久久| 成年女人永久免费观看视频| 老司机午夜十八禁免费视频| 99国产精品一区二区三区| 在线观看一区二区三区| 午夜福利在线观看吧| 国产大屁股一区二区在线视频| 身体一侧抽搐| 露出奶头的视频| 男人狂女人下面高潮的视频| .国产精品久久| 久久人人爽人人爽人人片va | 欧美一级a爱片免费观看看| 久久欧美精品欧美久久欧美| 岛国在线免费视频观看| 国内精品久久久久精免费| 成人一区二区视频在线观看| 欧美区成人在线视频| 亚洲国产欧洲综合997久久,| 简卡轻食公司| 欧美精品啪啪一区二区三区| 久久精品国产清高在天天线| 麻豆成人午夜福利视频| 亚洲无线观看免费| 国产高清激情床上av| 麻豆国产av国片精品| 好看av亚洲va欧美ⅴa在| av欧美777| 亚洲av免费高清在线观看| 亚洲国产日韩欧美精品在线观看| 久久亚洲精品不卡| 51国产日韩欧美| 国语自产精品视频在线第100页| 色噜噜av男人的天堂激情| 18禁裸乳无遮挡免费网站照片| 看片在线看免费视频| 欧美绝顶高潮抽搐喷水| 国产免费av片在线观看野外av| 99在线视频只有这里精品首页| 欧美日韩国产亚洲二区| 91久久精品电影网| 亚洲av第一区精品v没综合| 老女人水多毛片| 久久欧美精品欧美久久欧美| 日韩av在线大香蕉| 少妇的逼水好多| 亚洲美女视频黄频| 18禁黄网站禁片午夜丰满| 激情在线观看视频在线高清| 欧美性猛交╳xxx乱大交人| 日日摸夜夜添夜夜添小说| 亚洲精品成人久久久久久| aaaaa片日本免费| 亚洲成人免费电影在线观看| 我要看日韩黄色一级片| 亚洲天堂国产精品一区在线| 国产精品影院久久| 在线观看午夜福利视频| 成人永久免费在线观看视频| 欧美黑人欧美精品刺激| 亚洲欧美精品综合久久99| 成人毛片a级毛片在线播放| 性插视频无遮挡在线免费观看| 国产亚洲精品综合一区在线观看| 亚洲av.av天堂| 一进一出抽搐动态| 波野结衣二区三区在线| 欧美性感艳星| 男女那种视频在线观看| 欧美在线黄色| 午夜激情欧美在线| 国产精品久久久久久久电影| 69人妻影院| 淫妇啪啪啪对白视频| 91av网一区二区| 精品熟女少妇八av免费久了| 变态另类成人亚洲欧美熟女| 午夜亚洲福利在线播放| 脱女人内裤的视频| 亚洲综合色惰| 一级av片app| ponron亚洲| 老鸭窝网址在线观看| 91在线观看av| 美女高潮喷水抽搐中文字幕| 亚洲第一欧美日韩一区二区三区| 12—13女人毛片做爰片一| 女人十人毛片免费观看3o分钟| 久久国产精品人妻蜜桃| 午夜福利18| 国产精品免费一区二区三区在线| 如何舔出高潮| 男人舔奶头视频| 99在线视频只有这里精品首页| 国语自产精品视频在线第100页| 久久中文看片网| 极品教师在线视频| 亚洲精品在线观看二区| 亚洲五月婷婷丁香| 午夜免费男女啪啪视频观看 | 欧美色欧美亚洲另类二区| 日本黄大片高清| 中国美女看黄片| 内射极品少妇av片p| 久久久久久久午夜电影| 亚洲在线自拍视频| 一区福利在线观看| 亚洲人成伊人成综合网2020| 久久精品91蜜桃| 狠狠狠狠99中文字幕| 久久久久久久久中文| 美女高潮喷水抽搐中文字幕| 国产成人啪精品午夜网站| 午夜影院日韩av| 欧美日韩福利视频一区二区| 国产视频一区二区在线看| 欧美在线一区亚洲| 欧美高清成人免费视频www| 他把我摸到了高潮在线观看| 精品午夜福利视频在线观看一区| 色在线成人网| 国产成人av教育| 精品一区二区三区人妻视频| 久久天躁狠狠躁夜夜2o2o| 日本熟妇午夜| 超碰av人人做人人爽久久| 看黄色毛片网站| 哪里可以看免费的av片| 久久伊人香网站| 在线观看一区二区三区| 亚洲欧美精品综合久久99| 日韩成人在线观看一区二区三区| 日韩亚洲欧美综合| 亚洲国产高清在线一区二区三| 一区二区三区免费毛片| 国产精品国产高清国产av| 观看美女的网站| 国产亚洲精品av在线| 欧美日韩中文字幕国产精品一区二区三区| 亚洲在线自拍视频| 99久久精品热视频| 亚洲av不卡在线观看| 久久国产乱子免费精品| 桃色一区二区三区在线观看| 国产av在哪里看| 两人在一起打扑克的视频| 国产午夜精品久久久久久一区二区三区 | 女人十人毛片免费观看3o分钟| 91狼人影院| 免费看美女性在线毛片视频| 看黄色毛片网站| 嫩草影院精品99| 国产亚洲精品av在线| 国产精品自产拍在线观看55亚洲| 看片在线看免费视频| 90打野战视频偷拍视频| 成人精品一区二区免费| 色综合站精品国产| a级一级毛片免费在线观看| 久久久久久大精品| 免费av观看视频| 一级毛片久久久久久久久女| 91狼人影院| 最近中文字幕高清免费大全6 | 色视频www国产| 成人三级黄色视频| 91麻豆av在线| 日韩中字成人| 国产精品久久电影中文字幕| 热99在线观看视频| 欧美国产日韩亚洲一区| 男女做爰动态图高潮gif福利片| 蜜桃久久精品国产亚洲av| 日本 av在线| www.熟女人妻精品国产| 久久久久久久久大av| 精品久久久久久久久久久久久| 男人舔女人下体高潮全视频| 亚洲美女搞黄在线观看 | 国产亚洲欧美98| 国产午夜精品久久久久久一区二区三区 | 美女cb高潮喷水在线观看| 欧美日韩黄片免| 亚洲国产色片| 中文字幕av成人在线电影| 午夜激情欧美在线| av天堂在线播放| 最新在线观看一区二区三区| 亚洲精品日韩av片在线观看| 极品教师在线视频| 精品久久久久久久久亚洲 | 一级av片app| 精华霜和精华液先用哪个| 成熟少妇高潮喷水视频| 草草在线视频免费看| 尤物成人国产欧美一区二区三区| 免费电影在线观看免费观看| 亚洲av成人精品一区久久| 一二三四社区在线视频社区8| 精品一区二区三区av网在线观看|