• <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ù)分布的貝葉斯推斷
    久久国产精品大桥未久av| 黑人巨大精品欧美一区二区蜜桃 | 欧美国产精品va在线观看不卡| 两个人看的免费小视频| 五月天丁香电影| 久久精品国产亚洲av天美| 国产xxxxx性猛交| 久久国产精品男人的天堂亚洲 | 亚洲性久久影院| 日韩成人伦理影院| 久久久久久久久久久免费av| 国产精品久久久久久久电影| 午夜福利在线观看免费完整高清在| 精品少妇黑人巨大在线播放| 黄色 视频免费看| 久久精品夜色国产| 国内精品宾馆在线| 人妻人人澡人人爽人人| 欧美人与性动交α欧美软件| 天天躁日日躁夜夜躁夜夜| 757午夜福利合集在线观看| 妹子高潮喷水视频| 在线观看一区二区三区激情| 无遮挡黄片免费观看| 国产精品久久久久久精品古装| 最新美女视频免费是黄的| 视频区图区小说| 久久天躁狠狠躁夜夜2o2o| 天天操日日干夜夜撸| 最近最新免费中文字幕在线| av超薄肉色丝袜交足视频| 一区二区三区精品91| 热99re8久久精品国产| 国产成人一区二区三区免费视频网站| 亚洲av欧美aⅴ国产| 老鸭窝网址在线观看| 国产成人精品在线电影| 久久精品人人爽人人爽视色| 亚洲久久久国产精品| 久久精品aⅴ一区二区三区四区| 在线永久观看黄色视频| 日韩人妻精品一区2区三区| 国产亚洲欧美精品永久| 9191精品国产免费久久| 国产精品1区2区在线观看. | 性少妇av在线| 国产av又大| 欧美成狂野欧美在线观看| 在线观看日韩欧美| 久久九九热精品免费| 久久性视频一级片| 亚洲av欧美aⅴ国产| 久久亚洲真实| 一进一出抽搐动态| 亚洲成人手机| 涩涩av久久男人的天堂| 69av精品久久久久久| 人妻久久中文字幕网| 国产精品乱码一区二三区的特点 | av一本久久久久| 亚洲成国产人片在线观看| av超薄肉色丝袜交足视频| av中文乱码字幕在线| 身体一侧抽搐| 欧美日韩亚洲高清精品| 黄色成人免费大全| 精品国产超薄肉色丝袜足j| 亚洲美女黄片视频| 交换朋友夫妻互换小说| 制服人妻中文乱码| 国产成人欧美| 国产欧美亚洲国产| 中文字幕人妻熟女乱码| av电影中文网址| 天天躁狠狠躁夜夜躁狠狠躁| 久久久久国内视频| 丰满迷人的少妇在线观看| 精品一区二区三区视频在线观看免费 | 丁香六月欧美| 欧美精品高潮呻吟av久久| 国产男女超爽视频在线观看| 欧美亚洲 丝袜 人妻 在线| 亚洲精华国产精华精| 啪啪无遮挡十八禁网站| 久久天躁狠狠躁夜夜2o2o| 欧美在线一区亚洲| 一区二区三区精品91| 久久香蕉激情| 三级毛片av免费| 国产单亲对白刺激| 成人影院久久| 国产成人精品久久二区二区免费| 欧美性长视频在线观看| 欧美日韩黄片免| 老熟女久久久| 久久人人97超碰香蕉20202| 亚洲国产欧美日韩在线播放| 午夜激情av网站| 国产免费男女视频| www日本在线高清视频| 欧美日韩av久久| 国产亚洲一区二区精品| 午夜影院日韩av| www.999成人在线观看| 男女之事视频高清在线观看| 一级片免费观看大全| 18禁裸乳无遮挡免费网站照片 | 在线观看免费视频日本深夜| 午夜两性在线视频| 午夜精品国产一区二区电影| av天堂在线播放| 两个人免费观看高清视频| 亚洲熟妇中文字幕五十中出 | 久久亚洲真实| 久久久国产成人免费| 一边摸一边抽搐一进一小说 | 精品高清国产在线一区| 91在线观看av| 曰老女人黄片| 少妇粗大呻吟视频| 亚洲av电影在线进入| 女警被强在线播放| 9色porny在线观看| 欧美最黄视频在线播放免费 | 久久中文字幕人妻熟女| 久久人妻福利社区极品人妻图片| 亚洲国产精品一区二区三区在线| 国产成人影院久久av| 国产精品亚洲av一区麻豆| 精品久久久精品久久久| 好男人电影高清在线观看| 国产主播在线观看一区二区| 亚洲美女黄片视频| a级片在线免费高清观看视频| 日本黄色日本黄色录像| 国产亚洲精品久久久久5区| 女警被强在线播放| 国产精品一区二区免费欧美| 正在播放国产对白刺激| 在线国产一区二区在线| 一边摸一边做爽爽视频免费| 国产精品免费一区二区三区在线 | 9191精品国产免费久久| 国产不卡av网站在线观看| 久久午夜亚洲精品久久| 亚洲国产欧美网| 精品国产一区二区三区四区第35| 一级作爱视频免费观看| 亚洲一区二区三区不卡视频| 在线永久观看黄色视频| 女人精品久久久久毛片| a级毛片黄视频| 亚洲男人天堂网一区| 99久久国产精品久久久| 又黄又爽又免费观看的视频| 亚洲精品美女久久久久99蜜臀| 国产成人av教育| 欧美性长视频在线观看| 黄色视频不卡| 国产精品永久免费网站| 亚洲少妇的诱惑av| 国产精品成人在线| 色老头精品视频在线观看| 一进一出好大好爽视频| 亚洲精品乱久久久久久| 久久久久久人人人人人| 交换朋友夫妻互换小说| 亚洲国产欧美网| 午夜福利乱码中文字幕| 99久久精品国产亚洲精品| 欧美日韩av久久| 久久天堂一区二区三区四区| 一级片免费观看大全| 99热网站在线观看| 人成视频在线观看免费观看| 色94色欧美一区二区| 性色av乱码一区二区三区2| 夜夜爽天天搞| 亚洲人成伊人成综合网2020| av一本久久久久| 一区二区三区国产精品乱码| 国产亚洲一区二区精品| 亚洲精品久久成人aⅴ小说| 国产欧美亚洲国产| 久久久久久久午夜电影 | 精品国产亚洲在线| 亚洲午夜精品一区,二区,三区| 好男人电影高清在线观看| videosex国产| 国产麻豆69| 欧美精品av麻豆av| 成人精品一区二区免费| 建设人人有责人人尽责人人享有的| ponron亚洲| 一区二区日韩欧美中文字幕| 久久国产精品大桥未久av| 成人手机av| 91老司机精品| 欧美国产精品va在线观看不卡| 18禁黄网站禁片午夜丰满| 久久国产精品人妻蜜桃| 欧美成人午夜精品| 18禁国产床啪视频网站| 成年人黄色毛片网站| 我的亚洲天堂| 成年版毛片免费区| 乱人伦中国视频| 99在线人妻在线中文字幕 | aaaaa片日本免费| 亚洲一区二区三区不卡视频| 日韩成人在线观看一区二区三区| 一区在线观看完整版| 在线十欧美十亚洲十日本专区| 99re6热这里在线精品视频| 亚洲精品乱久久久久久| www.自偷自拍.com| 欧美乱妇无乱码| 亚洲专区字幕在线| 国产精品二区激情视频| 欧洲精品卡2卡3卡4卡5卡区| 精品亚洲成国产av| 欧美日韩福利视频一区二区| 国产精品亚洲一级av第二区| 久久人妻熟女aⅴ| 人妻丰满熟妇av一区二区三区 | 久久久精品国产亚洲av高清涩受| 下体分泌物呈黄色| 中文字幕另类日韩欧美亚洲嫩草| 免费一级毛片在线播放高清视频 | 国产精品久久久av美女十八| 视频区图区小说| 在线观看66精品国产| 国产男女超爽视频在线观看| 久久人妻av系列| 欧美av亚洲av综合av国产av| 精品国产一区二区久久| 999精品在线视频| 亚洲av第一区精品v没综合| 高潮久久久久久久久久久不卡| 精品一区二区三卡| 黄色片一级片一级黄色片| 亚洲专区国产一区二区| 99国产精品一区二区蜜桃av | 国产精品久久久久久人妻精品电影| 国产1区2区3区精品| 丰满的人妻完整版| 亚洲avbb在线观看| 免费黄频网站在线观看国产| 69精品国产乱码久久久| 一区二区三区激情视频| 男女下面插进去视频免费观看| 久久 成人 亚洲| 另类亚洲欧美激情| 国产精品免费大片| 婷婷精品国产亚洲av在线 | 午夜福利视频在线观看免费| 午夜精品在线福利| 黄色成人免费大全| 亚洲黑人精品在线| 久久久久国产精品人妻aⅴ院 | 激情在线观看视频在线高清 | 九色亚洲精品在线播放| 亚洲av第一区精品v没综合| 一级毛片精品| www.999成人在线观看| 国产激情久久老熟女| 新久久久久国产一级毛片| 亚洲熟女精品中文字幕| 丰满饥渴人妻一区二区三| 亚洲av第一区精品v没综合| 美女福利国产在线| 一进一出好大好爽视频| 亚洲精品乱久久久久久| 精品亚洲成a人片在线观看| 男人操女人黄网站| 亚洲成人免费电影在线观看| 久久精品91无色码中文字幕| 岛国在线观看网站| 久久人妻熟女aⅴ| 99国产精品99久久久久| 国产一区二区激情短视频| 欧美大码av| 午夜福利欧美成人| 精品久久蜜臀av无| 国产一区有黄有色的免费视频| bbb黄色大片| 777久久人妻少妇嫩草av网站| 国产有黄有色有爽视频| 久久精品国产综合久久久| 久久久久精品人妻al黑| 亚洲欧美日韩高清在线视频| 日韩欧美一区二区三区在线观看 | 亚洲全国av大片| 青草久久国产| 动漫黄色视频在线观看| 久久人妻福利社区极品人妻图片| 欧美精品av麻豆av| 一进一出好大好爽视频| 少妇的丰满在线观看| 九色亚洲精品在线播放| 高清黄色对白视频在线免费看| 亚洲国产欧美网| 在线天堂中文资源库| 久久影院123| 国产精品香港三级国产av潘金莲| 777米奇影视久久| 日韩制服丝袜自拍偷拍| 久9热在线精品视频| 亚洲欧美日韩另类电影网站| 久久亚洲真实| 亚洲欧美精品综合一区二区三区| 亚洲精华国产精华精| 亚洲国产精品合色在线| 久久久久久久精品吃奶| 成年人午夜在线观看视频| 少妇裸体淫交视频免费看高清 | 一边摸一边抽搐一进一出视频| 99国产综合亚洲精品| 一级毛片高清免费大全| 国产精品影院久久| 热99re8久久精品国产| 久久婷婷成人综合色麻豆| 一级片免费观看大全| 热99国产精品久久久久久7| 极品教师在线免费播放| 美国免费a级毛片| 午夜福利在线免费观看网站| 国产成人免费无遮挡视频| 天堂动漫精品| 亚洲精品久久午夜乱码| 欧美激情极品国产一区二区三区| 啦啦啦视频在线资源免费观看| 91精品三级在线观看| 国产成人啪精品午夜网站| 波多野结衣av一区二区av| 亚洲成a人片在线一区二区| 欧美+亚洲+日韩+国产| 在线观看日韩欧美| 好男人电影高清在线观看| 欧美日韩福利视频一区二区| 老熟妇仑乱视频hdxx| 久久久国产成人免费| 亚洲精品中文字幕在线视频| 一级作爱视频免费观看| 伊人久久大香线蕉亚洲五| 日韩有码中文字幕| 老熟女久久久| 在线av久久热| 国产精品久久久av美女十八| 大片电影免费在线观看免费| 日韩欧美一区视频在线观看| 亚洲欧洲精品一区二区精品久久久| 久久亚洲精品不卡| 国产片内射在线| 人妻 亚洲 视频| 大片电影免费在线观看免费| 人妻 亚洲 视频| 国产精品永久免费网站| 欧美老熟妇乱子伦牲交| 高潮久久久久久久久久久不卡| 亚洲av美国av| 国产激情欧美一区二区| 91成年电影在线观看| 午夜日韩欧美国产| 少妇裸体淫交视频免费看高清 | 建设人人有责人人尽责人人享有的| 一边摸一边抽搐一进一小说 | 久久这里只有精品19| 自线自在国产av| 欧美日韩福利视频一区二区| 亚洲伊人色综图| 色婷婷av一区二区三区视频| 国产在线观看jvid| 国产精品免费大片| 宅男免费午夜| 国产99久久九九免费精品| 黄色视频不卡| 欧美日韩国产mv在线观看视频| 亚洲性夜色夜夜综合| 香蕉久久夜色| 高清黄色对白视频在线免费看| 国产精品国产高清国产av | 国产欧美日韩一区二区三区在线| x7x7x7水蜜桃| 亚洲欧美色中文字幕在线| 婷婷丁香在线五月| 国产精品国产av在线观看| 国产高清国产精品国产三级| 国产精品二区激情视频| 真人做人爱边吃奶动态| 日韩欧美免费精品| 欧美日韩av久久| 最近最新免费中文字幕在线| 精品国产亚洲在线| 日本一区二区免费在线视频| 精品亚洲成国产av| 国产精品久久电影中文字幕 | 狂野欧美激情性xxxx| 啦啦啦免费观看视频1| 亚洲精品国产一区二区精华液| 精品人妻在线不人妻| 成人影院久久| 国产欧美亚洲国产| 在线观看舔阴道视频| 中文字幕另类日韩欧美亚洲嫩草| 成年版毛片免费区| 久久精品亚洲av国产电影网| 亚洲专区中文字幕在线| 国产单亲对白刺激| 国产av一区二区精品久久| 国产高清国产精品国产三级| 久久久水蜜桃国产精品网| 啦啦啦视频在线资源免费观看| 精品久久久久久久毛片微露脸| 好看av亚洲va欧美ⅴa在| 久久精品国产亚洲av高清一级| 国产乱人伦免费视频| 这个男人来自地球电影免费观看| 欧美亚洲日本最大视频资源| 精品无人区乱码1区二区| 国产野战对白在线观看| 欧美中文综合在线视频| 18禁国产床啪视频网站| 免费一级毛片在线播放高清视频 | 国产精品九九99| 精品国产国语对白av| 午夜亚洲福利在线播放| 欧美日韩av久久| 高清毛片免费观看视频网站 | 欧美av亚洲av综合av国产av| 免费av中文字幕在线| 亚洲,欧美精品.| 精品国产乱子伦一区二区三区| 成人免费观看视频高清| 国产精品乱码一区二三区的特点 | а√天堂www在线а√下载 | 免费在线观看亚洲国产| 一a级毛片在线观看| 正在播放国产对白刺激| 欧美激情极品国产一区二区三区| 久久精品亚洲精品国产色婷小说| 一级,二级,三级黄色视频| 国产99白浆流出| 一进一出抽搐动态| 巨乳人妻的诱惑在线观看| av一本久久久久| 亚洲av成人av| 国产在视频线精品| 久久天躁狠狠躁夜夜2o2o| 黄色毛片三级朝国网站| 亚洲av日韩精品久久久久久密| 好看av亚洲va欧美ⅴa在| cao死你这个sao货| 国产一区二区三区视频了| 国产精品秋霞免费鲁丝片| 亚洲国产精品sss在线观看 | 91精品国产国语对白视频| 成在线人永久免费视频| 国产精品99久久99久久久不卡| 欧美精品一区二区免费开放| 99国产极品粉嫩在线观看| 一级片'在线观看视频| 国产不卡av网站在线观看| 日日夜夜操网爽| 香蕉丝袜av| 国产精品99久久99久久久不卡| 中文欧美无线码| 999久久久国产精品视频| 久久久国产成人精品二区 | 亚洲人成电影观看| 国产深夜福利视频在线观看| 亚洲精品国产一区二区精华液| 亚洲av美国av| 欧美乱码精品一区二区三区| 午夜福利在线观看吧| 久久久久久人人人人人| 777米奇影视久久| 久久中文看片网| 看黄色毛片网站| 大型av网站在线播放| 精品国产超薄肉色丝袜足j| 国产成人一区二区三区免费视频网站| 国产精品香港三级国产av潘金莲| 精品国产乱子伦一区二区三区| 黄色视频,在线免费观看| 色老头精品视频在线观看| 久久性视频一级片| 午夜影院日韩av| av一本久久久久| 五月开心婷婷网| 久久精品亚洲熟妇少妇任你| 久久九九热精品免费| av中文乱码字幕在线| 欧美国产精品一级二级三级| 国产精品九九99| www日本在线高清视频| 国产精华一区二区三区| 女人高潮潮喷娇喘18禁视频| 伦理电影免费视频| 久久国产精品人妻蜜桃| 亚洲人成77777在线视频| 亚洲午夜精品一区,二区,三区| tocl精华| 日本欧美视频一区| 国产精品一区二区免费欧美| videosex国产| 手机成人av网站| 久久香蕉激情| 国产不卡一卡二| 免费观看精品视频网站| а√天堂www在线а√下载 | 国产成人欧美在线观看 | 欧美日韩国产mv在线观看视频| a级毛片在线看网站| 精品电影一区二区在线| 亚洲国产中文字幕在线视频| 国产三级黄色录像| 久久精品亚洲精品国产色婷小说| 一区二区三区精品91| 国产日韩欧美亚洲二区| 91大片在线观看| 99精国产麻豆久久婷婷| 国产亚洲av高清不卡| 老司机午夜十八禁免费视频| 亚洲视频免费观看视频| 天天躁日日躁夜夜躁夜夜| 国产黄色免费在线视频| 人人妻,人人澡人人爽秒播| www.熟女人妻精品国产| 国产男女超爽视频在线观看| 在线观看免费视频日本深夜| 两人在一起打扑克的视频| 国产91精品成人一区二区三区| 亚洲av成人一区二区三| 国产成人影院久久av| 男女午夜视频在线观看| 久久热在线av| 性色av乱码一区二区三区2| 超碰97精品在线观看| 久久久久视频综合| 老汉色av国产亚洲站长工具| a级毛片黄视频| 精品熟女少妇八av免费久了| 精品视频人人做人人爽| 少妇猛男粗大的猛烈进出视频| 亚洲国产欧美一区二区综合| 老熟女久久久| 成人18禁高潮啪啪吃奶动态图| 俄罗斯特黄特色一大片| 动漫黄色视频在线观看| 久久久久国产一级毛片高清牌| 精品人妻熟女毛片av久久网站| 好看av亚洲va欧美ⅴa在| 久久性视频一级片| 在线观看一区二区三区激情| 精品第一国产精品| 51午夜福利影视在线观看| 悠悠久久av| 熟女少妇亚洲综合色aaa.| 女性生殖器流出的白浆| 又紧又爽又黄一区二区| 十分钟在线观看高清视频www| 99精品欧美一区二区三区四区| 国产精品乱码一区二三区的特点 | 欧美日韩成人在线一区二区| 精品久久久久久,| 亚洲精品久久午夜乱码| 欧美精品人与动牲交sv欧美| 亚洲伊人色综图| 久久热在线av| 满18在线观看网站| 精品电影一区二区在线| 亚洲精品国产色婷婷电影| 欧美日韩亚洲综合一区二区三区_| 国产99白浆流出| 高清av免费在线| 人人妻人人澡人人爽人人夜夜| 中文亚洲av片在线观看爽 | 亚洲第一青青草原| 欧美性长视频在线观看| 一进一出抽搐动态| 国产91精品成人一区二区三区| 又黄又粗又硬又大视频| 午夜精品久久久久久毛片777| 精品人妻在线不人妻| 国产亚洲精品一区二区www | 日本欧美视频一区| 在线观看免费高清a一片| 欧美亚洲 丝袜 人妻 在线| 亚洲第一av免费看| 啦啦啦在线免费观看视频4| 麻豆av在线久日| 午夜福利视频在线观看免费| av国产精品久久久久影院| 色精品久久人妻99蜜桃| 久久国产乱子伦精品免费另类| 国产精品电影一区二区三区 | 日韩欧美国产一区二区入口| 亚洲情色 制服丝袜| 交换朋友夫妻互换小说| 国产男靠女视频免费网站| 精品国产乱子伦一区二区三区| 老鸭窝网址在线观看| 91成人精品电影| 窝窝影院91人妻| a在线观看视频网站| 欧美在线黄色| 亚洲午夜理论影院| 欧美色视频一区免费| 另类亚洲欧美激情| 99国产精品一区二区三区| 国产国语露脸激情在线看| 人妻一区二区av| 欧美精品亚洲一区二区| 18禁裸乳无遮挡动漫免费视频|