• <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ù)分布的貝葉斯推斷
    一边亲一边摸免费视频| 18禁裸乳无遮挡动漫免费视频| 精品欧美一区二区三区在线| 欧美变态另类bdsm刘玥| 免费在线观看视频国产中文字幕亚洲 | 宅男免费午夜| 精品视频人人做人人爽| 97在线人人人人妻| 最近最新中文字幕大全免费视频 | 精品福利观看| 久久99一区二区三区| 各种免费的搞黄视频| 久久久久久久大尺度免费视频| 日韩av不卡免费在线播放| 中文字幕最新亚洲高清| 99精国产麻豆久久婷婷| 十八禁人妻一区二区| 日本一区二区免费在线视频| 国产精品亚洲av一区麻豆| 免费av中文字幕在线| xxx大片免费视频| 国产日韩一区二区三区精品不卡| 黑人猛操日本美女一级片| 久久国产精品影院| 黄片播放在线免费| 涩涩av久久男人的天堂| 最近最新中文字幕大全免费视频 | 一二三四社区在线视频社区8| 9热在线视频观看99| 热99国产精品久久久久久7| 精品高清国产在线一区| 久久青草综合色| 黄色怎么调成土黄色| 50天的宝宝边吃奶边哭怎么回事| 人妻 亚洲 视频| 精品国产一区二区三区四区第35| 国产真人三级小视频在线观看| kizo精华| 国产色视频综合| 男女国产视频网站| 免费看av在线观看网站| 五月开心婷婷网| 日本午夜av视频| 中国国产av一级| 九色亚洲精品在线播放| 久久久欧美国产精品| 男女免费视频国产| 国产精品免费大片| 亚洲五月色婷婷综合| 在线观看免费视频网站a站| 免费观看av网站的网址| 大码成人一级视频| 9热在线视频观看99| 欧美97在线视频| 少妇裸体淫交视频免费看高清 | 十分钟在线观看高清视频www| 大片免费播放器 马上看| 久久久久久久国产电影| 黑丝袜美女国产一区| 国产极品粉嫩免费观看在线| 午夜日韩欧美国产| 国产在线观看jvid| 久久久亚洲精品成人影院| 国产精品三级大全| 欧美成人午夜精品| 亚洲,欧美,日韩| 欧美亚洲 丝袜 人妻 在线| svipshipincom国产片| 97人妻天天添夜夜摸| 成年动漫av网址| 最近最新中文字幕大全免费视频 | 欧美黑人欧美精品刺激| 欧美日韩一级在线毛片| 一级毛片我不卡| 精品人妻一区二区三区麻豆| 日本a在线网址| 国产精品人妻久久久影院| 日韩免费高清中文字幕av| 97在线人人人人妻| 交换朋友夫妻互换小说| 国产精品一区二区在线不卡| 国产成人欧美| 精品久久久久久久毛片微露脸 | 免费久久久久久久精品成人欧美视频| 麻豆乱淫一区二区| 久久鲁丝午夜福利片| 免费观看a级毛片全部| 久久久久精品国产欧美久久久 | 久久人人97超碰香蕉20202| 啦啦啦 在线观看视频| 最黄视频免费看| 国产成人a∨麻豆精品| 成年av动漫网址| 精品久久久久久久毛片微露脸 | 亚洲欧洲精品一区二区精品久久久| 欧美精品亚洲一区二区| 男女床上黄色一级片免费看| 人人妻人人澡人人看| 少妇猛男粗大的猛烈进出视频| 18禁观看日本| 日本a在线网址| 大片电影免费在线观看免费| 国产精品国产三级国产专区5o| 亚洲第一青青草原| 精品国产一区二区久久| 婷婷色综合www| 国产真人三级小视频在线观看| 人成视频在线观看免费观看| 久久久久久久久免费视频了| 精品亚洲乱码少妇综合久久| 亚洲国产欧美一区二区综合| 久久99热这里只频精品6学生| 国产高清视频在线播放一区 | e午夜精品久久久久久久| 丰满少妇做爰视频| 热99久久久久精品小说推荐| 久久精品久久久久久噜噜老黄| 午夜福利视频精品| 亚洲欧美色中文字幕在线| 交换朋友夫妻互换小说| av在线播放精品| 18禁观看日本| 伊人亚洲综合成人网| 欧美另类一区| 91字幕亚洲| 国产av精品麻豆| 中文字幕色久视频| 国产日韩欧美视频二区| 亚洲成色77777| 女人精品久久久久毛片| 亚洲欧美一区二区三区国产| 黄色视频不卡| 2018国产大陆天天弄谢| 日韩中文字幕欧美一区二区 | 欧美 亚洲 国产 日韩一| 欧美在线黄色| 亚洲色图 男人天堂 中文字幕| 午夜福利视频精品| 精品人妻在线不人妻| 日韩大码丰满熟妇| netflix在线观看网站| 在线观看一区二区三区激情| 国产精品国产三级专区第一集| 亚洲情色 制服丝袜| 女人被躁到高潮嗷嗷叫费观| 久久久久久久久久久久大奶| 50天的宝宝边吃奶边哭怎么回事| 久久国产精品男人的天堂亚洲| 十分钟在线观看高清视频www| 亚洲精品乱久久久久久| 久久 成人 亚洲| 韩国高清视频一区二区三区| 国产亚洲精品第一综合不卡| 久久女婷五月综合色啪小说| 一边亲一边摸免费视频| 成人国产一区最新在线观看 | 欧美精品亚洲一区二区| videosex国产| 伦理电影免费视频| 涩涩av久久男人的天堂| 91国产中文字幕| 亚洲欧美日韩高清在线视频 | 老汉色∧v一级毛片| 在线天堂中文资源库| 欧美黄色淫秽网站| 一区二区三区激情视频| 大香蕉久久网| 国产一区亚洲一区在线观看| 悠悠久久av| 成人亚洲精品一区在线观看| 国产片内射在线| 晚上一个人看的免费电影| 国产av精品麻豆| 一级毛片电影观看| 精品亚洲乱码少妇综合久久| 国产精品成人在线| 国产精品久久久久成人av| 欧美变态另类bdsm刘玥| 成人国语在线视频| 晚上一个人看的免费电影| 中文字幕另类日韩欧美亚洲嫩草| 9色porny在线观看| 丁香六月天网| 国产一区二区三区av在线| 丁香六月欧美| 久久久精品免费免费高清| 最黄视频免费看| 中文乱码字字幕精品一区二区三区| 亚洲第一青青草原| 色综合欧美亚洲国产小说| 国产精品久久久久久精品古装| 国产免费现黄频在线看| 免费av中文字幕在线| 18禁观看日本| 国产精品偷伦视频观看了| 精品高清国产在线一区| 人人妻人人爽人人添夜夜欢视频| 成年女人毛片免费观看观看9 | 日本a在线网址| 9色porny在线观看| 精品一品国产午夜福利视频| 亚洲国产av新网站| 99热全是精品| 狠狠婷婷综合久久久久久88av| 国产亚洲午夜精品一区二区久久| 嫩草影视91久久| 一级,二级,三级黄色视频| 99国产精品99久久久久| 少妇被粗大的猛进出69影院| 十分钟在线观看高清视频www| 亚洲熟女毛片儿| 欧美日韩综合久久久久久| 国产男人的电影天堂91| www.自偷自拍.com| 亚洲美女黄色视频免费看| 咕卡用的链子| 黄色视频不卡| www.熟女人妻精品国产| 国产在线观看jvid| 一区二区av电影网| 亚洲精品成人av观看孕妇| 成人亚洲精品一区在线观看| 亚洲欧洲精品一区二区精品久久久| 午夜免费成人在线视频| 香蕉国产在线看| 一级毛片我不卡| 亚洲国产精品999| 国产精品亚洲av一区麻豆| 免费看av在线观看网站| 久久精品成人免费网站| 在线 av 中文字幕| 99精国产麻豆久久婷婷| 国产1区2区3区精品| 赤兔流量卡办理| 国产成人免费观看mmmm| 亚洲 国产 在线| 这个男人来自地球电影免费观看| 在线 av 中文字幕| 久久久国产精品麻豆| 亚洲精品国产色婷婷电影| 伊人久久大香线蕉亚洲五| 精品福利永久在线观看| 热re99久久国产66热| 另类亚洲欧美激情| 啦啦啦在线观看免费高清www| 久久久久久久国产电影| 亚洲熟女精品中文字幕| 国产视频一区二区在线看| 久久女婷五月综合色啪小说| 国精品久久久久久国模美| av视频免费观看在线观看| 少妇人妻 视频| 又紧又爽又黄一区二区| 在线看a的网站| 久久 成人 亚洲| 亚洲精品一二三| 一本色道久久久久久精品综合| 久久久久久免费高清国产稀缺| 男女免费视频国产| 99久久综合免费| 国产亚洲av片在线观看秒播厂| 成人亚洲精品一区在线观看| 色精品久久人妻99蜜桃| 久久ye,这里只有精品| 亚洲欧美中文字幕日韩二区| 国产精品麻豆人妻色哟哟久久| 国产男女超爽视频在线观看| 欧美久久黑人一区二区| 一区福利在线观看| 国产精品秋霞免费鲁丝片| 秋霞在线观看毛片| 天天操日日干夜夜撸| 午夜福利视频精品| 黄色一级大片看看| 精品一区二区三区av网在线观看 | 国产片特级美女逼逼视频| 国产欧美亚洲国产| 菩萨蛮人人尽说江南好唐韦庄| 亚洲精品国产一区二区精华液| 午夜91福利影院| 精品人妻一区二区三区麻豆| 叶爱在线成人免费视频播放| 免费久久久久久久精品成人欧美视频| 黄色视频在线播放观看不卡| a级片在线免费高清观看视频| 香蕉国产在线看| 国产av国产精品国产| 亚洲欧美日韩另类电影网站| 1024香蕉在线观看| 久9热在线精品视频| 亚洲av片天天在线观看| 乱人伦中国视频| 精品国产乱码久久久久久男人| 国产亚洲欧美精品永久| 一二三四在线观看免费中文在| 一本一本久久a久久精品综合妖精| 国产欧美日韩综合在线一区二区| 欧美乱码精品一区二区三区| 激情视频va一区二区三区| 中文字幕人妻丝袜一区二区| 麻豆国产av国片精品| 精品少妇内射三级| 男人添女人高潮全过程视频| 十八禁人妻一区二区| 国产免费视频播放在线视频| 日韩 欧美 亚洲 中文字幕| 日韩人妻精品一区2区三区| 好男人视频免费观看在线| 男女床上黄色一级片免费看| 91老司机精品| 亚洲久久久国产精品| 亚洲欧洲精品一区二区精品久久久| 男人添女人高潮全过程视频| 国产亚洲欧美精品永久| 青春草视频在线免费观看| 桃花免费在线播放| 亚洲伊人色综图| 精品熟女少妇八av免费久了| 少妇粗大呻吟视频| 日韩伦理黄色片| 超碰97精品在线观看| cao死你这个sao货| 久久久久网色| 满18在线观看网站| 国产精品久久久久久精品电影小说| 黄频高清免费视频| av不卡在线播放| 人人妻人人澡人人爽人人夜夜| 人人妻人人添人人爽欧美一区卜| 欧美人与善性xxx| 激情五月婷婷亚洲| 国产免费又黄又爽又色| 亚洲色图 男人天堂 中文字幕| 亚洲国产欧美网| 亚洲 国产 在线| 亚洲七黄色美女视频| 国产日韩一区二区三区精品不卡| 精品亚洲乱码少妇综合久久| 男女免费视频国产| 欧美日韩视频高清一区二区三区二| 亚洲国产精品国产精品| 一级片免费观看大全| 国产成人免费无遮挡视频| 欧美在线黄色| 久9热在线精品视频| 亚洲欧美清纯卡通| 亚洲精品久久成人aⅴ小说| 国产黄频视频在线观看| 丰满少妇做爰视频| kizo精华| 菩萨蛮人人尽说江南好唐韦庄| 亚洲专区国产一区二区| 国产精品九九99| 亚洲精品国产一区二区精华液| 18禁观看日本| 亚洲一区二区三区欧美精品| 国语对白做爰xxxⅹ性视频网站| 色94色欧美一区二区| 午夜福利视频精品| 久久女婷五月综合色啪小说| av电影中文网址| 叶爱在线成人免费视频播放| 不卡av一区二区三区| 国产熟女欧美一区二区| 欧美人与性动交α欧美软件| 国产女主播在线喷水免费视频网站| 欧美人与性动交α欧美软件| 国产女主播在线喷水免费视频网站| 中文乱码字字幕精品一区二区三区| 在线看a的网站| 色婷婷久久久亚洲欧美| 亚洲黑人精品在线| 亚洲欧美精品综合一区二区三区| 亚洲激情五月婷婷啪啪| 在线观看www视频免费| 只有这里有精品99| 欧美 日韩 精品 国产| 最黄视频免费看| 国产爽快片一区二区三区| 国产国语露脸激情在线看| 韩国精品一区二区三区| 最新在线观看一区二区三区 | 久久99精品国语久久久| 每晚都被弄得嗷嗷叫到高潮| 欧美人与性动交α欧美精品济南到| 午夜视频精品福利| 亚洲伊人久久精品综合| 免费黄频网站在线观看国产| 高清av免费在线| 男女无遮挡免费网站观看| 2018国产大陆天天弄谢| 1024香蕉在线观看| 亚洲色图综合在线观看| 国产男女超爽视频在线观看| 欧美人与性动交α欧美软件| 一级片免费观看大全| 97精品久久久久久久久久精品| 国产成人影院久久av| 精品国产超薄肉色丝袜足j| 日本av免费视频播放| 飞空精品影院首页| 日韩一卡2卡3卡4卡2021年| 在线观看人妻少妇| 色综合欧美亚洲国产小说| 久久亚洲精品不卡| 久久午夜综合久久蜜桃| 国产欧美亚洲国产| 咕卡用的链子| 国产精品久久久久久人妻精品电影 | 日韩大片免费观看网站| 美女中出高潮动态图| 一级毛片电影观看| 精品少妇黑人巨大在线播放| 大片免费播放器 马上看| 欧美激情高清一区二区三区| 黄色视频在线播放观看不卡| 久久九九热精品免费| 欧美日韩亚洲高清精品| 久久久精品国产亚洲av高清涩受| 国产在线观看jvid| 欧美精品av麻豆av| 岛国毛片在线播放| 丰满少妇做爰视频| 欧美精品人与动牲交sv欧美| 热re99久久精品国产66热6| 美女主播在线视频| 涩涩av久久男人的天堂| 国产精品三级大全| www.999成人在线观看| 91精品伊人久久大香线蕉| 国产熟女欧美一区二区| 九色亚洲精品在线播放| 日本色播在线视频| 国产亚洲精品第一综合不卡| 多毛熟女@视频| 久久久欧美国产精品| 校园人妻丝袜中文字幕| 最近中文字幕2019免费版| 老司机深夜福利视频在线观看 | 国产片内射在线| 99国产精品免费福利视频| 啦啦啦在线免费观看视频4| 国产在线免费精品| 一本久久精品| 人成视频在线观看免费观看| 亚洲精品在线美女| 777米奇影视久久| 精品人妻1区二区| 777米奇影视久久| 热99久久久久精品小说推荐| 啦啦啦中文免费视频观看日本| 成人免费观看视频高清| 老司机亚洲免费影院| 日韩一本色道免费dvd| 国产无遮挡羞羞视频在线观看| 18禁黄网站禁片午夜丰满| 天天躁狠狠躁夜夜躁狠狠躁| 男女高潮啪啪啪动态图| 岛国毛片在线播放| 国产精品熟女久久久久浪| 久久精品成人免费网站| 中文字幕制服av| 操美女的视频在线观看| 日韩欧美一区视频在线观看| 欧美日韩视频高清一区二区三区二| 90打野战视频偷拍视频| 欧美日韩av久久| 欧美日韩视频精品一区| 久久久国产欧美日韩av| 国产高清国产精品国产三级| videosex国产| 丝袜脚勾引网站| 两性夫妻黄色片| 两性夫妻黄色片| www.自偷自拍.com| 日韩av免费高清视频| 男女床上黄色一级片免费看| xxxhd国产人妻xxx| 99国产精品免费福利视频| 黄片播放在线免费| 99re6热这里在线精品视频| 久久ye,这里只有精品| kizo精华| 激情五月婷婷亚洲| 男人添女人高潮全过程视频| 成年女人毛片免费观看观看9 | 丝袜喷水一区| 国产欧美日韩综合在线一区二区| 国产成人欧美| 91国产中文字幕| 热re99久久精品国产66热6| 免费人妻精品一区二区三区视频| 水蜜桃什么品种好| 欧美在线一区亚洲| 亚洲午夜精品一区,二区,三区| 日本a在线网址| 热99国产精品久久久久久7| 欧美xxⅹ黑人| 老司机亚洲免费影院| 中文字幕人妻丝袜制服| 国产在视频线精品| 交换朋友夫妻互换小说| 亚洲国产成人一精品久久久| 成年女人毛片免费观看观看9 | 夫妻午夜视频| 成人免费观看视频高清| 欧美精品亚洲一区二区| 亚洲欧美成人综合另类久久久| 免费少妇av软件| 天堂8中文在线网| 久久中文字幕一级| 男人爽女人下面视频在线观看| 男的添女的下面高潮视频| 亚洲国产精品成人久久小说| 日韩一本色道免费dvd| 9热在线视频观看99| 欧美精品一区二区大全| 97人妻天天添夜夜摸| 日韩一本色道免费dvd| 久久国产精品影院| 一级,二级,三级黄色视频| 一二三四在线观看免费中文在| 嫁个100分男人电影在线观看 | 欧美亚洲 丝袜 人妻 在线| 久久精品久久久久久久性| 免费观看a级毛片全部| 欧美 日韩 精品 国产| 国产极品粉嫩免费观看在线| 亚洲 国产 在线| 精品久久久久久电影网| 新久久久久国产一级毛片| 国产成人欧美在线观看 | 亚洲九九香蕉| 免费在线观看视频国产中文字幕亚洲 | 国产真人三级小视频在线观看| 美女扒开内裤让男人捅视频| 人妻 亚洲 视频| av一本久久久久| 午夜福利,免费看| 丁香六月欧美| 亚洲av日韩精品久久久久久密 | 国产日韩欧美在线精品| 国产激情久久老熟女| 亚洲成av片中文字幕在线观看| av不卡在线播放| 成人黄色视频免费在线看| 免费少妇av软件| av在线播放精品| 久久久久久久久久久久大奶| 乱人伦中国视频| 成人18禁高潮啪啪吃奶动态图| 久久人妻福利社区极品人妻图片 | 亚洲免费av在线视频| 午夜福利,免费看| 2021少妇久久久久久久久久久| 亚洲天堂av无毛| 男人添女人高潮全过程视频| 国产免费视频播放在线视频| 中国美女看黄片| 女人高潮潮喷娇喘18禁视频| 人人澡人人妻人| 满18在线观看网站| 男人舔女人的私密视频| 欧美精品av麻豆av| 免费观看av网站的网址| 欧美日韩亚洲高清精品| 欧美少妇被猛烈插入视频| 夫妻性生交免费视频一级片| 天堂俺去俺来也www色官网| 精品久久蜜臀av无| 美女高潮到喷水免费观看| 侵犯人妻中文字幕一二三四区| 国产成人免费观看mmmm| www.999成人在线观看| 手机成人av网站| 亚洲国产av影院在线观看| 久久久欧美国产精品| 无限看片的www在线观看| 日本av免费视频播放| avwww免费| 真人做人爱边吃奶动态| 国产成人影院久久av| 在线观看一区二区三区激情| 久久 成人 亚洲| 亚洲精品久久久久久婷婷小说| 在线观看人妻少妇| 国产av国产精品国产| 亚洲一区二区三区欧美精品| 波多野结衣av一区二区av| 日日夜夜操网爽| 丝袜人妻中文字幕| 国产一区有黄有色的免费视频| www.熟女人妻精品国产| 黄色视频不卡| av不卡在线播放| 丝袜人妻中文字幕| 在线亚洲精品国产二区图片欧美| 七月丁香在线播放| 男人添女人高潮全过程视频| 欧美久久黑人一区二区| 国产成人精品在线电影| 亚洲午夜精品一区,二区,三区| 一区二区三区精品91| 99久久人妻综合| 少妇的丰满在线观看| 午夜免费观看性视频| 国产精品一区二区免费欧美 | 久久天躁狠狠躁夜夜2o2o | 欧美大码av| 国产成人精品久久二区二区91| 91字幕亚洲| 青青草视频在线视频观看| 国产精品久久久久久人妻精品电影 | 大片免费播放器 马上看| 欧美成狂野欧美在线观看|