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

    EM算法在基于區(qū)間數(shù)據(jù)的加速壽命試驗(yàn)統(tǒng)計(jì)分析中的應(yīng)用

    2017-09-03 10:57:06邱雷顰
    關(guān)鍵詞:指數(shù)分布廣義區(qū)間

    邱雷顰

    (閩江師范高等??茖W(xué)校計(jì)算機(jī)系,福州 350002)

    EM算法在基于區(qū)間數(shù)據(jù)的加速壽命試驗(yàn)統(tǒng)計(jì)分析中的應(yīng)用

    邱雷顰

    (閩江師范高等專科學(xué)校計(jì)算機(jī)系,福州 350002)

    在步進(jìn)加速壽命試驗(yàn)中,當(dāng)壽命服從廣義指數(shù)分布且獲得的數(shù)據(jù)是區(qū)間數(shù)據(jù)時(shí),給出試驗(yàn)安排,并通過“各應(yīng)力下產(chǎn)品失效機(jī)理保持不變”等幾個(gè)基本模型假定,得出各應(yīng)力水平下形狀參數(shù)相等的結(jié)論以及不同應(yīng)力水平下的試驗(yàn)時(shí)間t1-t0,t2-t1,…,ti-1-ti-2折算到某一應(yīng)力水平下的時(shí)間τi-1的折算公式,進(jìn)而得出求相關(guān)參數(shù)的極大似然估計(jì)的隱性表達(dá)式。進(jìn)一步討論利用EM算法求解,先由E步求出期望,再由M步求出使得期望極大化的點(diǎn),給出具體的迭代過程。最后采用Monte Carlo數(shù)據(jù)模擬方法分別在大樣本和小樣本場(chǎng)合下給出了參數(shù)真值的估計(jì)值,結(jié)果表明,該方法在大樣本場(chǎng)合下更具有效性。

    區(qū)間數(shù)據(jù);廣義指數(shù)分布;EM算法

    引 言

    加速壽命試驗(yàn),是指將樣品置于超過正常應(yīng)力水平下,通過觀察樣品的失效時(shí)間(即壽命),從而利用統(tǒng)計(jì)方法推斷在正常應(yīng)力條件下產(chǎn)品的各項(xiàng)可靠性指標(biāo)的一種壽命試驗(yàn)。目前常見的類型有恒定應(yīng)力加速壽命試驗(yàn),步進(jìn)應(yīng)力加速壽命試驗(yàn),序進(jìn)應(yīng)力加速壽命試驗(yàn)。文獻(xiàn)[1]中討論了加速壽命試驗(yàn)的類型和理論基礎(chǔ)。本文討論步進(jìn)應(yīng)力加速壽命試驗(yàn),簡(jiǎn)稱步加試驗(yàn),是先選定一組加速應(yīng)力水平S1<…

    然而,針對(duì)Gupta和Kundu提出的廣義指數(shù)分布[11],利用EM算法對(duì)區(qū)間數(shù)據(jù)的統(tǒng)計(jì)分析的研究并不多見。文獻(xiàn)[12]指出,在分析很多壽命數(shù)據(jù)時(shí),廣義指數(shù)分布往往比其它分布更能有效利用。近年來,雙參數(shù)廣義指數(shù)分布已經(jīng)很廣泛的運(yùn)用到分析壽命數(shù)據(jù)。文獻(xiàn)[13]研究了廣義指數(shù)分布的尺度參數(shù)的極大似然估計(jì)基于隨機(jī)截尾模型。文獻(xiàn)[14]介紹了形狀參數(shù)的增加,使得這一分布族產(chǎn)生了更多的靈活性,使它可以被用來分析刪失數(shù)據(jù),如區(qū)間數(shù)據(jù)。

    本文將主要討論當(dāng)獲得的壽命數(shù)據(jù)為區(qū)間數(shù)據(jù),而壽命分布服從廣義指數(shù)分布時(shí)相關(guān)參數(shù)的統(tǒng)計(jì)分析。

    1 試驗(yàn)安排及基本模型假定

    1.1 試驗(yàn)安排

    1.2 基本假定

    假定1 在應(yīng)力水平Si(i=1,2,…,m)下,產(chǎn)品的壽命分布為廣義指數(shù)分布,分布函數(shù)為Fi(t)=(1-e-λit)αi,密度函數(shù)為fi(t)=αiλi(1-e-λit)αi-1·e-λit,其中λi和αi分別是尺度參數(shù)的倒數(shù)和形狀參數(shù),和應(yīng)力水平Si有關(guān)。

    假定2 在應(yīng)力水平Si(i=1,2,…,m)下,產(chǎn)品的失效機(jī)理保持不變[1]。

    假定3 產(chǎn)品在應(yīng)力水平S下,滿足線性加速模型lnλi=a+bφ(Si),其中φ(Si)是與應(yīng)力水平Si有關(guān)的已知函數(shù),下記φi,i=1,2,…,m。

    假定4 產(chǎn)品的剩余壽命僅依賴于當(dāng)時(shí)已累積失效的部分和當(dāng)時(shí)的應(yīng)力水平,而與累積的方式無關(guān)[15]。

    2 參數(shù)估計(jì)

    引理1 在應(yīng)力水平Si(i=1,2,…,m)下,產(chǎn)品的壽命分布為Fi(t)=(1-e-λit)αi,則αi≡α(i=1,2,…,m)。

    由基本假定2,在應(yīng)力水平Si(i=1,2,…,m)下,產(chǎn)品的失效機(jī)理保持不變,即加速系數(shù)是與R無關(guān)的常數(shù),因此要求αi≡α(i=1,2,…,m)。

    引理2 若記τi-1為產(chǎn)品經(jīng)歷了水平S1,S2,…,Si-1的步加試驗(yàn)后,各段試驗(yàn)延續(xù)時(shí)間t1-t0,t2-t1,…,ti-1-ti-2折算到Si應(yīng)力水平下總的折算時(shí)間,則

    證明 由基本假定4知τi-1滿足

    又由假定1和假定2得:

    (1-e-λiτi-1)α=(1-e-λi-1(ti-1-ti-2+τi-2))α

    得到,

    e-λiτi-1=e-λi-1(ti-1-ti-2+τi-2)

    由此得到遞推公式

    又由τ0=0得

    其中,i=2,3,…,m。由基本假定知,

    所以,τi-1為b的函數(shù),下記τi-1(b),

    i=1,2,…,m。

    估計(jì)a,b,α三個(gè)參數(shù)。記n個(gè)樣品的失效時(shí)間為Xk,k=1,2,…,n,它們是獨(dú)立同分布的隨機(jī)變量,密度函數(shù)為f(x)=αλ(1-e-λx)α-1·e-λx,

    x>0,λ>0,α>0。只能觀測(cè)到落在區(qū)間[ti,ti+1)中的樣本數(shù)ni,其中i=0,1,2,…,m,0=t0

    x落在區(qū)間[ti,ti+1)中的概率為

    P(x∈[ti,ti+1))=Fi+1(ti+1)-Fi(ti)=

    [1-e-λi+1(ti+1-ti+τi(b))]α-[1-e-λi(ti-ti-1+τi-1(b))]α

    其中記t-1=τ-1(b)=0,λ0=λm+1=0。

    定理1 參數(shù)a,b,α的極大似然估計(jì)可以由隱性表達(dá)式求解。

    (1)

    (2)

    (3)

    其中,

    證明 容易得到似然函數(shù)為

    對(duì)數(shù)似然函數(shù)為

    3 EM算法求極大似然估計(jì)

    式(1)、式(2)和式(3)無法得出a,b,α的顯性表達(dá)式,于是考慮用EM算法實(shí)現(xiàn)。

    為了便于表述,n個(gè)樣品的失效時(shí)間Xk,k=1,2,…,n的全體記為X,觀測(cè)結(jié)果即落在區(qū)間[ti,ti+1)中的樣本數(shù)目ni,

    i=0,1,2,…,m,

    0=t0

    于是

    [1-exp[-ea(k)+b(k)φi(t-ti-1+τi-1(b(k)))]]α(k)-1·

    exp[-ea(k)+b(k)φi(t-ti-1+τi-1(b(k)))]

    其中,

    [1-exp[-ea(k)+b(k)φi(ti-ti-1+τi-1(b(k)))]]α(k)-

    [1-exp[-ea(k)+b(k)φiτi-1(b(k))]]α(k)

    所以

    其中,

    AiAi(a,b,t)1-exp[-ea+bφi(t-ti-1+

    τi-1(b))]

    (4)

    a=

    (5)

    (6)

    迭代過程:

    (1)當(dāng)k=0時(shí),給定參數(shù)初始值a(0),b(0),α(0),及任意小的正數(shù)ε。

    (2)當(dāng)k≥1時(shí),現(xiàn)有的參數(shù)估計(jì)為a(k),b(k),α(k)。

    ①取式(4)右邊的,b=b(k),即可找到α(k+1),

    ②取式(5)右邊的b=b(k),a=a(k),α=α(k),即可找到a(k+1),

    ③令式(6)為零,即可解出b的顯示表達(dá)式,令表達(dá)式中a=a(k),b=b(k),α=α(k)找到b(k+1),

    其中,

    (t-ti-1+τi-1(b(k)))]

    滿足迭代終止的條件后繼續(xù)進(jìn)行20次迭代,取所得到的20個(gè)參數(shù)迭代值的平均值作為對(duì)真實(shí)參數(shù)的估計(jì)

    4 數(shù)據(jù)模擬

    考慮步進(jìn)加速壽命試驗(yàn)過程中的應(yīng)力為溫度,其中T1=240,T2=350,T3=470,φi=φ(Ti)=1/Ti,t0=0,t1=1.20×104,t2=1.55×104,t3=1.66×104,參數(shù)真值為a=-20,b=104,α=2,取迭代初始值a(0)=-18,b(0)=12000,α(0)=2.5。

    在小樣本場(chǎng)合下,取n0=n1=n2=5,ε=10-3,在大樣本的場(chǎng)合下,取n0=n1=n2=100,ε=10-3運(yùn)用以上迭代過程得到參數(shù)真值的估計(jì),見表1。

    表1 模擬結(jié)果

    由數(shù)據(jù)模擬結(jié)果表1可知,在大樣本場(chǎng)合下,估計(jì)值與參數(shù)真值的相對(duì)誤差有所減小。因此,在產(chǎn)品壽命服從廣義指數(shù)分布的情況下,基于獲得的數(shù)據(jù)為區(qū)間數(shù)據(jù)的步進(jìn)加速壽命試驗(yàn),在大樣本場(chǎng)合下利用EM算法給出參數(shù)的最大似然估計(jì)顯得更為有效。

    [1] YANG Y H,ZHOU Y Q.Theoretical foundation of accelerated life testing[J].Journal of Propulsion and Technology,2001,22(5):353-356.

    [2] ZHAO W,ELSAYED E A.A general accelerated life model for step-stress testing[J].IIE Trans,2005,37:1059-1069.

    [3] MOHAMED T.Step-Stress Accelerated Life Tests[C]//Miodrag Lovric.International Encyclopedia of Statistical Science,Springer,America:2011:1547-1517.

    [4] DEMPSTER A P,LAIRD N M,RUBIN D B.Maximum likelihood from incomplete data via the EM algorithm[J].Math SciNet Google Scholar.JR Stat Soc Ser B,1977,39:1-22.

    [5] WU C F.On the convergence properties of the EM algorithm[J].The Annals of Statistics,1983,11(1):95-103.

    [6] AMBROISE C.The EM Algorithm and Extensions[J].Journal of classification,1988,15(1):154-156.

    [7] CHARLES B,PAUL P B.Handbook of Mathematical Methods in Imaging[M].New York:Springer-Verlag,2015.

    [8] HUANG W H,CHEN Y G.The multiset EM algorithm[J].Statistics & Probability Letters,2017,126:41-48.

    [9] SUVRA P,BALAKRISHNAN N.Destructive negative binomial cure rate model and EM-based likelihood inference under Weibull lifetime[J].Statistics & Probability Letters,2016,116:9-20.

    [10] JEAN-PATRICK B,GILLES C.EM for mixtures[J].Statistics and Computing,2015,25(4):713-726.

    [11] GUPTA R D,KUNDU D.Generalized exponential distributions[J].Australian & New Zealand Journal of Statistics,1999,41:173-188.

    [12] CHIANG-SHENG L,HSINE-JEN T.A note on the generalized linear exponential distribution[J].Statistics & Probability Letters,2017,124(5):49-54.

    [13] NAMHYUN K.Approximate MLE for the scale parameter of the generalized exponential distribution under random censoring[J].Journal of the Korean Statistical Society,2014,43(1):119-131.

    [14] DEBASIS K,RAMESHWAR D G.An extension of the generalized exponential distribution[J].Statistical Methodology,2011,8(6):485-496.

    [15] NELSON W B.Accelerated life testing-step-stress models and data analysis[J].IEEE Transaction on Reliability,1980,29(2):103-108.

    Statistical Analysis Based on EM Algorithm for Accelerated Life Testing Under Interval Censored Samples

    QIULeipin

    (Department of Computer Science, Minjiang Teachers College, Fuzhou 350002, China)

    An accelerated life test carried out by step stress is considered when the the life time follows a generalized exponential distribution and the failure data obtained is interval censored. The test procedure is presented and the exchange formula of time is given by some assumptions. With the help of the expectation-maximization (EM) algorithm which is widely used when the observations can be viewed as incomplete data, the maximum likelihood estimates are computed. Moreover, an example by Mote Carlo data simulation is given to illustrate the procedure and show that this method is available, especially in large sample case.

    interval censored; generalized exponential distribution; EM algorithm

    2017-06-06

    福建省教育廳中青年課題(JAT160827)

    邱雷顰(1980-),女,福建石獅人,講師,碩士,主要從事概率與數(shù)理統(tǒng)計(jì)方面的研究,(E-mail)qiuleipin@163.com

    1673-1549(2017)04-0081-06

    10.11863/j.suse.2017.04.15

    0212.7

    A

    猜你喜歡
    指數(shù)分布廣義區(qū)間
    解兩類含參數(shù)的復(fù)合不等式有解與恒成立問題
    你學(xué)會(huì)“區(qū)間測(cè)速”了嗎
    Rn中的廣義逆Bonnesen型不等式
    從廣義心腎不交論治慢性心力衰竭
    指數(shù)分布抽樣基本定理及在指數(shù)分布參數(shù)統(tǒng)計(jì)推斷中的應(yīng)用
    有限群的廣義交換度
    區(qū)間對(duì)象族的可鎮(zhèn)定性分析
    二元Weinman型指數(shù)分布隨機(jī)變量之和、差、積、商及比率的分布
    指數(shù)分布與其它分布的關(guān)系
    定數(shù)截尾數(shù)據(jù)缺失場(chǎng)合下雙參數(shù)指數(shù)分布的貝葉斯推斷
    国产精品一及| 精品人妻偷拍中文字幕| 在线观看免费视频日本深夜| 精品福利观看| 精品福利观看| 欧美三级亚洲精品| 免费在线观看影片大全网站| 国产一级毛片七仙女欲春2| 国产不卡一卡二| 日韩强制内射视频| 热99re8久久精品国产| 亚洲国产精品久久男人天堂| 欧美不卡视频在线免费观看| 深夜a级毛片| 国产一区二区激情短视频| 欧美xxxx黑人xx丫x性爽| 免费看av在线观看网站| 91狼人影院| 国产熟女欧美一区二区| 免费av毛片视频| 我要搜黄色片| 麻豆av噜噜一区二区三区| 99精品在免费线老司机午夜| 午夜福利在线观看免费完整高清在 | a级毛片免费高清观看在线播放| 国产成人影院久久av| 成熟少妇高潮喷水视频| 亚洲国产精品合色在线| 精品久久久久久,| 色哟哟·www| 禁无遮挡网站| 久久久久久久亚洲中文字幕| 在线观看66精品国产| 成人三级黄色视频| 一区二区三区四区激情视频 | 色播亚洲综合网| 欧美xxxx性猛交bbbb| 日本一本二区三区精品| 成人特级黄色片久久久久久久| 少妇被粗大猛烈的视频| 午夜亚洲福利在线播放| 极品教师在线视频| 国产精品亚洲美女久久久| 精品人妻偷拍中文字幕| 欧美性猛交黑人性爽| 99在线视频只有这里精品首页| 中文在线观看免费www的网站| 欧美日韩瑟瑟在线播放| 成人高潮视频无遮挡免费网站| 国产乱人视频| 国产91精品成人一区二区三区| а√天堂www在线а√下载| xxxwww97欧美| 熟妇人妻久久中文字幕3abv| 日本色播在线视频| 超碰av人人做人人爽久久| 热99re8久久精品国产| 欧美潮喷喷水| 又粗又爽又猛毛片免费看| 亚洲av成人av| 国产人妻一区二区三区在| 亚洲av成人精品一区久久| 熟女电影av网| 99在线人妻在线中文字幕| 男人的好看免费观看在线视频| 亚洲美女搞黄在线观看 | 亚洲国产精品sss在线观看| 国产精品自产拍在线观看55亚洲| 夜夜夜夜夜久久久久| 久久久久久久久中文| 12—13女人毛片做爰片一| 久久久久久久久久黄片| 午夜精品久久久久久毛片777| 日本一本二区三区精品| 国产爱豆传媒在线观看| 亚洲最大成人av| 一进一出抽搐动态| 中文资源天堂在线| 日本色播在线视频| 久久午夜福利片| 久久久久久久久久久丰满 | 俄罗斯特黄特色一大片| 欧美中文日本在线观看视频| 日韩,欧美,国产一区二区三区 | 女的被弄到高潮叫床怎么办 | 午夜老司机福利剧场| 欧美高清性xxxxhd video| 少妇高潮的动态图| 久久久久久久久大av| 中亚洲国语对白在线视频| 中文字幕免费在线视频6| 看免费成人av毛片| 99久久中文字幕三级久久日本| 中国美女看黄片| 国产av一区在线观看免费| 亚洲色图av天堂| 搡女人真爽免费视频火全软件 | 亚洲va在线va天堂va国产| 在线播放国产精品三级| 精品一区二区三区人妻视频| 欧美性猛交黑人性爽| 三级毛片av免费| 免费一级毛片在线播放高清视频| 日本五十路高清| 露出奶头的视频| 中文字幕精品亚洲无线码一区| 国产激情偷乱视频一区二区| 国产免费一级a男人的天堂| 国产高清视频在线观看网站| 成人性生交大片免费视频hd| av在线亚洲专区| 无遮挡黄片免费观看| 欧美一区二区亚洲| 一卡2卡三卡四卡精品乱码亚洲| 亚洲国产精品合色在线| 免费不卡的大黄色大毛片视频在线观看 | 最近最新免费中文字幕在线| 亚洲久久久久久中文字幕| 日韩在线高清观看一区二区三区 | 日韩欧美免费精品| 亚洲国产精品久久男人天堂| 一区二区三区高清视频在线| 亚洲精品国产成人久久av| 午夜激情福利司机影院| 国产精品99久久久久久久久| 美女免费视频网站| 午夜免费激情av| 精品国内亚洲2022精品成人| 91av网一区二区| 久久亚洲真实| 国产伦精品一区二区三区四那| 日韩一本色道免费dvd| 国产精品精品国产色婷婷| 欧美中文日本在线观看视频| 全区人妻精品视频| 我要看日韩黄色一级片| 亚洲美女搞黄在线观看 | 免费电影在线观看免费观看| 好男人在线观看高清免费视频| av专区在线播放| av中文乱码字幕在线| 日韩欧美免费精品| 精品久久久久久成人av| 日韩欧美 国产精品| 长腿黑丝高跟| 嫁个100分男人电影在线观看| 国产精品一区www在线观看 | 在线播放国产精品三级| 欧美极品一区二区三区四区| 男女那种视频在线观看| 久久久久免费精品人妻一区二区| 女同久久另类99精品国产91| 一区二区三区高清视频在线| 99久久中文字幕三级久久日本| av在线天堂中文字幕| 欧美一区二区亚洲| 天天躁日日操中文字幕| 亚洲图色成人| 国产精品电影一区二区三区| 亚洲午夜理论影院| 精品一区二区三区人妻视频| 亚洲av中文av极速乱 | 久久精品国产亚洲网站| 中文字幕免费在线视频6| 国产精品野战在线观看| 全区人妻精品视频| 欧美bdsm另类| 丰满人妻一区二区三区视频av| 久久久成人免费电影| 日韩欧美免费精品| 非洲黑人性xxxx精品又粗又长| 国产大屁股一区二区在线视频| 搡老熟女国产l中国老女人| 老熟妇仑乱视频hdxx| 深夜精品福利| 一区二区三区四区激情视频 | 99九九线精品视频在线观看视频| 精品久久久久久久久久久久久| 人妻夜夜爽99麻豆av| a级毛片a级免费在线| 一卡2卡三卡四卡精品乱码亚洲| 夜夜看夜夜爽夜夜摸| 欧美性猛交黑人性爽| 精品日产1卡2卡| 熟女电影av网| а√天堂www在线а√下载| 国产综合懂色| 色在线成人网| 午夜精品在线福利| 国产探花极品一区二区| 精品久久国产蜜桃| 尤物成人国产欧美一区二区三区| 亚洲真实伦在线观看| 最近最新免费中文字幕在线| 在线播放无遮挡| 搡女人真爽免费视频火全软件 | 美女高潮的动态| 欧美激情久久久久久爽电影| 久久久久久久午夜电影| 国产在视频线在精品| 大型黄色视频在线免费观看| 性色avwww在线观看| 成人精品一区二区免费| 婷婷色综合大香蕉| 可以在线观看的亚洲视频| 欧美绝顶高潮抽搐喷水| 国产男人的电影天堂91| 日韩欧美在线二视频| 一区二区三区高清视频在线| 久久天躁狠狠躁夜夜2o2o| 美女被艹到高潮喷水动态| 婷婷色综合大香蕉| 午夜激情福利司机影院| 精品一区二区三区视频在线观看免费| 精品免费久久久久久久清纯| 看片在线看免费视频| 一a级毛片在线观看| 精品欧美国产一区二区三| 国产探花极品一区二区| 国产男人的电影天堂91| 别揉我奶头 嗯啊视频| 日本黄大片高清| 在线观看美女被高潮喷水网站| 老女人水多毛片| 男人舔奶头视频| 毛片一级片免费看久久久久 | 久久精品国产亚洲av香蕉五月| 精品欧美国产一区二区三| 免费看光身美女| 国产亚洲91精品色在线| 校园春色视频在线观看| 人人妻人人澡欧美一区二区| 国内精品宾馆在线| 亚洲va日本ⅴa欧美va伊人久久| 国产精品,欧美在线| 欧美色欧美亚洲另类二区| 国产白丝娇喘喷水9色精品| 精品久久久久久久久av| 精品人妻一区二区三区麻豆 | 人人妻,人人澡人人爽秒播| 中国美女看黄片| netflix在线观看网站| 日本 av在线| 亚洲欧美日韩卡通动漫| 精品久久久噜噜| 日本免费一区二区三区高清不卡| 国产一区二区亚洲精品在线观看| 特大巨黑吊av在线直播| 日韩欧美国产一区二区入口| 最后的刺客免费高清国语| 在线观看舔阴道视频| 国产精品免费一区二区三区在线| 亚洲欧美精品综合久久99| 三级国产精品欧美在线观看| 久久久久国内视频| 少妇高潮的动态图| 观看美女的网站| 日韩大尺度精品在线看网址| 国产私拍福利视频在线观看| 亚洲av电影不卡..在线观看| 精品人妻1区二区| a在线观看视频网站| 亚洲国产日韩欧美精品在线观看| 亚洲精品粉嫩美女一区| 亚洲在线自拍视频| 久久久成人免费电影| 国产亚洲欧美98| 国产午夜福利久久久久久| 亚洲欧美日韩卡通动漫| 特级一级黄色大片| 亚洲性夜色夜夜综合| 日韩欧美国产在线观看| 成人国产麻豆网| 欧美色欧美亚洲另类二区| 夜夜看夜夜爽夜夜摸| 两个人视频免费观看高清| 成人精品一区二区免费| 亚洲内射少妇av| 在线播放国产精品三级| 亚洲精品亚洲一区二区| 国产高清有码在线观看视频| 日本一本二区三区精品| 日韩精品青青久久久久久| 最近最新免费中文字幕在线| 日韩一区二区视频免费看| 亚洲成人中文字幕在线播放| 国产精品久久久久久av不卡| 久久精品国产鲁丝片午夜精品 | 欧美成人a在线观看| 精品人妻一区二区三区麻豆 | 天堂网av新在线| 国产精品综合久久久久久久免费| 亚洲av成人av| 午夜激情欧美在线| 91狼人影院| 热99在线观看视频| 亚洲四区av| av专区在线播放| 最好的美女福利视频网| 国产高清有码在线观看视频| 国产aⅴ精品一区二区三区波| 亚洲中文字幕一区二区三区有码在线看| 亚洲av.av天堂| 少妇猛男粗大的猛烈进出视频 | 白带黄色成豆腐渣| 97超视频在线观看视频| 国产麻豆成人av免费视频| 欧美人与善性xxx| 国产高清激情床上av| 99国产精品一区二区蜜桃av| 男女做爰动态图高潮gif福利片| 亚洲欧美清纯卡通| 亚洲成av人片在线播放无| 久久精品夜夜夜夜夜久久蜜豆| 中亚洲国语对白在线视频| 久久人妻av系列| 日韩中文字幕欧美一区二区| 亚洲av免费在线观看| 国产亚洲av嫩草精品影院| 国产欧美日韩精品一区二区| 嫩草影视91久久| 久久久久久久精品吃奶| 又爽又黄a免费视频| 成人av一区二区三区在线看| 狂野欧美激情性xxxx在线观看| 成年女人毛片免费观看观看9| 成人高潮视频无遮挡免费网站| 99热网站在线观看| 人妻少妇偷人精品九色| 久久久久免费精品人妻一区二区| 欧美日韩精品成人综合77777| 国产精品久久久久久av不卡| 国产精品一区二区免费欧美| 久久欧美精品欧美久久欧美| 在线免费十八禁| 国产亚洲精品久久久久久毛片| 少妇猛男粗大的猛烈进出视频 | 在线播放国产精品三级| 精品一区二区免费观看| 日本三级黄在线观看| 精华霜和精华液先用哪个| 国产亚洲欧美98| 97碰自拍视频| 韩国av在线不卡| 日韩人妻高清精品专区| 人妻少妇偷人精品九色| 久久热精品热| 午夜激情福利司机影院| 99国产极品粉嫩在线观看| 熟女人妻精品中文字幕| 嫩草影视91久久| 人妻少妇偷人精品九色| 男女啪啪激烈高潮av片| 国产成人av教育| 亚洲成人久久性| 亚洲av免费高清在线观看| 夜夜看夜夜爽夜夜摸| videossex国产| 精品福利观看| 在线观看午夜福利视频| 夜夜看夜夜爽夜夜摸| www.www免费av| 变态另类成人亚洲欧美熟女| 最后的刺客免费高清国语| 欧美最新免费一区二区三区| 久久天躁狠狠躁夜夜2o2o| 看十八女毛片水多多多| 欧美xxxx黑人xx丫x性爽| 亚洲av二区三区四区| 亚洲国产精品成人综合色| 又黄又爽又免费观看的视频| 男女视频在线观看网站免费| 婷婷精品国产亚洲av| 欧美激情在线99| 国内精品一区二区在线观看| 免费搜索国产男女视频| 在线a可以看的网站| 久久午夜福利片| 国语自产精品视频在线第100页| 丝袜美腿在线中文| av天堂中文字幕网| 午夜老司机福利剧场| 欧美日韩综合久久久久久 | 国产伦在线观看视频一区| 极品教师在线免费播放| 内射极品少妇av片p| 亚洲精品国产成人久久av| 少妇的逼水好多| 亚洲精品国产成人久久av| 午夜日韩欧美国产| 制服丝袜大香蕉在线| 欧美绝顶高潮抽搐喷水| 欧美在线一区亚洲| 在线a可以看的网站| 亚洲精品久久国产高清桃花| 亚洲第一区二区三区不卡| 给我免费播放毛片高清在线观看| 1000部很黄的大片| 美女 人体艺术 gogo| 神马国产精品三级电影在线观看| 国内精品久久久久久久电影| 自拍偷自拍亚洲精品老妇| 亚洲一区高清亚洲精品| 久久精品国产自在天天线| 精品99又大又爽又粗少妇毛片 | 亚洲在线观看片| 国产男人的电影天堂91| 波多野结衣巨乳人妻| 国产一区二区三区在线臀色熟女| 91狼人影院| 嫩草影院新地址| 成人高潮视频无遮挡免费网站| 欧美激情在线99| 看十八女毛片水多多多| 观看美女的网站| 免费人成在线观看视频色| 久久人人爽人人爽人人片va| 精品久久久久久久人妻蜜臀av| 床上黄色一级片| 精品久久久久久久人妻蜜臀av| 88av欧美| 久久精品国产鲁丝片午夜精品 | 国产在线精品亚洲第一网站| 日韩精品有码人妻一区| 国产69精品久久久久777片| 国产男靠女视频免费网站| 干丝袜人妻中文字幕| 精华霜和精华液先用哪个| 欧美极品一区二区三区四区| 亚洲欧美日韩无卡精品| 婷婷丁香在线五月| 国产精品久久久久久精品电影| 99久久成人亚洲精品观看| 亚洲一级一片aⅴ在线观看| 欧美一级a爱片免费观看看| 99在线人妻在线中文字幕| 久久香蕉精品热| 天美传媒精品一区二区| 伦理电影大哥的女人| 成人无遮挡网站| 国产精品久久久久久久电影| 久99久视频精品免费| av天堂在线播放| 成年女人毛片免费观看观看9| 欧美bdsm另类| 老女人水多毛片| АⅤ资源中文在线天堂| netflix在线观看网站| 欧美日韩乱码在线| h日本视频在线播放| 春色校园在线视频观看| 日韩一本色道免费dvd| 国产亚洲欧美98| 少妇被粗大猛烈的视频| 久久人妻av系列| 国产一区二区在线av高清观看| 免费av观看视频| 嫩草影院新地址| 国产精品1区2区在线观看.| 国产亚洲精品久久久com| ponron亚洲| 成年免费大片在线观看| 日韩欧美精品v在线| 欧洲精品卡2卡3卡4卡5卡区| 亚洲国产欧洲综合997久久,| 亚洲国产日韩欧美精品在线观看| 丝袜美腿在线中文| 999久久久精品免费观看国产| 国产成人aa在线观看| 精品久久久久久,| 日本三级黄在线观看| 亚洲自偷自拍三级| 成人精品一区二区免费| 亚洲七黄色美女视频| 在线观看免费视频日本深夜| 能在线免费观看的黄片| 无遮挡黄片免费观看| 亚洲国产精品久久男人天堂| 久久人人爽人人爽人人片va| 亚洲性久久影院| 国产精品久久久久久av不卡| 国产免费一级a男人的天堂| 国内精品久久久久久久电影| 久久精品夜夜夜夜夜久久蜜豆| 国产高清三级在线| 久久久久久大精品| 亚洲av五月六月丁香网| 在现免费观看毛片| 亚洲乱码一区二区免费版| 中文字幕人妻熟人妻熟丝袜美| 国产av不卡久久| 国产成人一区二区在线| 中国美女看黄片| 91久久精品电影网| 国产v大片淫在线免费观看| 色综合亚洲欧美另类图片| 国产精品电影一区二区三区| 男女之事视频高清在线观看| 九九在线视频观看精品| 国产美女午夜福利| 少妇丰满av| 99热这里只有是精品50| 美女xxoo啪啪120秒动态图| 国产精品98久久久久久宅男小说| 亚洲欧美日韩高清专用| 国产精品福利在线免费观看| 国产久久久一区二区三区| 99在线人妻在线中文字幕| 久久久久久久久久成人| 嫁个100分男人电影在线观看| 97碰自拍视频| 亚洲精品乱码久久久v下载方式| 波多野结衣高清作品| 91狼人影院| 日韩欧美精品免费久久| 久久欧美精品欧美久久欧美| 亚洲不卡免费看| 精品久久国产蜜桃| 成人国产综合亚洲| 99精品久久久久人妻精品| 国产真实伦视频高清在线观看 | 九九热线精品视视频播放| 国产精品一区二区三区四区免费观看 | 中国美女看黄片| 中文字幕av成人在线电影| 亚洲av电影不卡..在线观看| 99在线人妻在线中文字幕| 日韩国内少妇激情av| 日韩高清综合在线| 婷婷精品国产亚洲av在线| 国产精品永久免费网站| 九九热线精品视视频播放| 网址你懂的国产日韩在线| 简卡轻食公司| 最新中文字幕久久久久| 日韩大尺度精品在线看网址| 日本a在线网址| 不卡一级毛片| 99久国产av精品| 国产蜜桃级精品一区二区三区| 麻豆一二三区av精品| 久久久久久久精品吃奶| 欧美一区二区精品小视频在线| 精品日产1卡2卡| 久久人人爽人人爽人人片va| 日韩,欧美,国产一区二区三区 | 国产三级中文精品| 国内久久婷婷六月综合欲色啪| 久久久久久久午夜电影| 国产黄片美女视频| 午夜激情福利司机影院| 丰满的人妻完整版| 韩国av一区二区三区四区| 国产成人影院久久av| 亚洲精品一区av在线观看| 国产aⅴ精品一区二区三区波| 极品教师在线免费播放| 亚洲国产高清在线一区二区三| 最新在线观看一区二区三区| 美女高潮喷水抽搐中文字幕| 国产成人一区二区在线| 亚洲最大成人手机在线| 美女免费视频网站| 可以在线观看毛片的网站| 亚洲黑人精品在线| 午夜福利高清视频| 日本三级黄在线观看| 国产精品爽爽va在线观看网站| 麻豆成人av在线观看| 亚洲成a人片在线一区二区| 国产精品国产三级国产av玫瑰| 欧美+日韩+精品| 国产精品嫩草影院av在线观看 | 国产麻豆成人av免费视频| 国产伦精品一区二区三区视频9| 日日摸夜夜添夜夜添小说| 亚洲精华国产精华精| 99国产极品粉嫩在线观看| 久久国产精品人妻蜜桃| 精品免费久久久久久久清纯| 国产高清视频在线播放一区| 亚洲不卡免费看| 久久久精品欧美日韩精品| 国产精品久久久久久久电影| 一级av片app| 床上黄色一级片| 黄色日韩在线| 久久99热6这里只有精品| 亚洲自偷自拍三级| 国产毛片a区久久久久| 51国产日韩欧美| 精品一区二区免费观看| 欧美日韩国产亚洲二区| 麻豆国产av国片精品| 一区福利在线观看| 亚洲中文日韩欧美视频| 97热精品久久久久久| 色在线成人网| 联通29元200g的流量卡| 亚洲精华国产精华精| 欧美日本亚洲视频在线播放| 亚洲人成网站高清观看| 亚洲av熟女| 日日啪夜夜撸| 淫妇啪啪啪对白视频| 91久久精品国产一区二区成人| 波多野结衣高清无吗| 国产在线精品亚洲第一网站| 亚洲aⅴ乱码一区二区在线播放| 亚洲av中文av极速乱 | 乱人视频在线观看| 欧美日韩黄片免| 蜜桃久久精品国产亚洲av| av在线天堂中文字幕| 人妻夜夜爽99麻豆av|