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

    定數(shù)截尾場合三參數(shù)pareto分布參數(shù)的最優(yōu)置信區(qū)間

    2022-04-27 02:22:32李云飛
    內(nèi)江師范學院學報 2022年4期
    關鍵詞:樞軸定數(shù)置信區(qū)間

    劉 璐, 李云飛

    (西華師范大學 數(shù)學與信息學院, 四川 南充 637009)

    0 引言

    Pareto分布是研究收入模型的一種分布,經(jīng)過不斷發(fā)展,由兩參數(shù)pareto分布推廣到三參數(shù)pareto分布和廣義pareto分布.Pareto 分布在社會經(jīng)濟、軍事以及可靠性統(tǒng)計分析等很多領域中得到廣泛的應用[1].因此,研究 pareto 分布,尤其是研究其參數(shù)估計問題,具有重要的理論意義和實用價值.

    關于pareto分布參數(shù)的最優(yōu)置信區(qū)間,許多專家和學者對其進行了研究,已經(jīng)形成了豐富的研究成果.Neyman[2]的理論,在給定置信水平,以保證有一定的可靠度下,精度往往用區(qū)間長度來衡量,區(qū)間長度越短說明對參數(shù)估計得越準確.Ali等[3]提出了一種推導廣義pareto分布形狀參數(shù)統(tǒng)計檢驗的一般方法,并證明了這種方法優(yōu)于經(jīng)典的極大似然比檢驗.Singh等[4]利用最大熵原理,提出了三參數(shù)廣義pareto分布參數(shù)估計一種新的方法.佟毅[5]研究并給出了定數(shù)截尾下指數(shù)分布平均壽命參數(shù)和平均失效率參數(shù)的最優(yōu)置信區(qū)間,通過算例分析并得出了小樣本場合下用最優(yōu)置信區(qū)間作為未知參數(shù)的置信區(qū)間更好.王娟等[6]給出Pareto分布中尺度參數(shù)的幾種區(qū)間估計方法,重點研究之前不常見的極大似然估計的漸近正態(tài)性法和輪廓似然函數(shù)法,并說明這幾種方法的適應范圍及優(yōu)缺點.徐美萍等[7]通過對比weibull分布尺度參數(shù)的最短置信區(qū)間和常用置信區(qū)間,發(fā)現(xiàn)研究小樣本情形時參數(shù)最短置信區(qū)間的重要性和必要性.王芳等[8]將兩參數(shù)的pareto分布情形推廣到了三參數(shù),并研究了尺度參數(shù)和形狀參數(shù)估計的漸進性質.劉榮玄等[9]在逐步增加首失效截尾樣本下,研究了三參數(shù)pareto分布族形狀參數(shù)的一致最小方差無偏估計.秦祖啟等[10]討論了樞軸量服從F分布,而F分布的密度函數(shù)是極不對稱的,基于此兩正態(tài)分布方差比優(yōu)化后的置信區(qū)間.蔡潔等[11]用最優(yōu)假設檢驗的統(tǒng)計量構造出兩正態(tài)總體方差比的樞軸量,并從最短置信區(qū)間的本質意義出發(fā),構造出求解最短置信區(qū)間的條件和證明出其解的存在唯一性,得出用最短置信區(qū)間作為方差比的區(qū)間估計,精度將會得到顯著提高.劉璐等[12]通過構造置信區(qū)間的一個樞軸量的方法,求定數(shù)截尾場合Pareto分布形狀參數(shù)的最短置信區(qū)間,并將其與傳統(tǒng)方法求得的置信區(qū)間進行對比分析得出最優(yōu)置信區(qū)間.

    本文將在定數(shù)截尾壽命試驗場合下,針對三參數(shù)pareto分布單總體參數(shù)的最優(yōu)置信區(qū)間和兩總體形狀參數(shù)比的最優(yōu)置信區(qū)間進行研究.

    1 單總體尺度參數(shù)σ的置信區(qū)間

    若隨機變量X服從三參數(shù)Pareto分布,其密度函數(shù)和分布函數(shù)分別為

    f(x;σ,θ,μ)=θσθ(x-μ)-(1+θ),x>σ+μ,

    (1)

    F(x;σ,θ,μ)=1-σθ(x-μ)-θ,

    (2)

    記為X~Pareto(θ,σ,μ),其中θ>0為形狀參數(shù),σ>0為尺度參數(shù),μ∈R為位置參數(shù),并假設為已知常數(shù).

    抽取來自總體Pareto(θ,σ,μ)的容量為n的一個樣本X1,X2,…,Xn,當有r個樣品失效時停止壽命試驗,失效數(shù)據(jù)為X(1)≤X(2)≤…≤X(r).

    1.1 用樞軸量法求σ的最短置信區(qū)間

    定理1設T1,T2,…,Tn是來自雙參數(shù)指數(shù)分布總體E(δ,θ)的容量為n的樣本,定數(shù)截尾數(shù)是r,T(1)≤T(2)≤…≤T(r)是前r個次序統(tǒng)計量,則:

    證明令H1=rT(1),Hi=(r-i+1)[T(i)-T(i-1)],i=2,3,…,r,則H1,H2,…,Hr相互獨立,且H1~E(r,δ,θ),H2,H3,…,Hr同分布于E(θ),那么H1-rδ,H2,…,Hr獨立同分布于E(θ),則:

    為了求δ的最短置信區(qū)間,假設存在a,b(0

    此時δ的置信區(qū)間為:

    這個區(qū)間的平均長度為:

    定理2當α0=0時,δ的置信度為1-α的最短置信區(qū)間為:

    設h(α0)=b-a,則

    (3)

    其中

    1.2 用輪廓似然函數(shù)法求σ的置信區(qū)間

    因為總體Pareto(θ,σ,μ)含有三個參數(shù),而現(xiàn)在只對尺度參數(shù)σ進行估計,所以可以通過輪廓似然函數(shù)法求其置信區(qū)間.

    令η=(θ,σ,μ)Τ,根據(jù)(1)和(2)得η在定數(shù)截尾壽命試驗下的似然函數(shù)為:

    從而η的對數(shù)似然函數(shù)為:

    (4)

    (5)

    要使lnL(xi,η)取得最大值,則σ的極大似然估計為:

    (6)

    θ的極大似然估計為:

    (7)

    由(4)式可以得到參數(shù)σ的輪廓似然函數(shù)為:

    將(5)式代入lnpL(xi,η)中得到:

    又因為參數(shù)σ的偏差度函數(shù)為:

    將(6)和(7)代入上式得到:

    所以,對于給定的α∈(0,1),構造σ的置信區(qū)間:

    (8)

    2 兩總體形狀參數(shù)比θ2/θ1的置信區(qū)間

    定理3設X1,X2,…,Xn1是來自總體Pareto(θ1,σ1,μ1)的樣本,Y1,Y2,…,Yn2來自總體Pareto(θ2,σ2,μ2)的樣本,兩個總體相互獨立.對它們做定數(shù)截尾壽命試驗后,得到失效數(shù)據(jù)分別為:X(1)≤X(2)≤…≤X(r),Y(1)≤Y(2)≤…≤Y(r),其中θ1,σ1,θ2,σ2均未知.

    分別令T=ln(X-μ1),Z=ln(Y-μ2),則:

    證明令T=ln(X-μ1),Z=ln(Y-μ2),則fT(t;δ1,θ1)=θ1e-θ1(t-δ1),其中δ1=lnσ1;fY(y;δ2,θ2)=θ2e-θ2(t-δ2),其中δ2=lnσ2.根據(jù)定理1的證明,同理可得:

    ~F(2(r1-1),2(r2-1)).

    因為Q2的分布F(2(r1-1),2(r2-1))是已知的,且此分布與θ1,θ2無關.故將Q2作為樞軸量來構造形狀參數(shù)比θ2/θ1的置信區(qū)間.因為樞軸量服從F分布,F(xiàn)分布的概率密度函數(shù)非對稱分布,利用分位數(shù)

    構造的傳統(tǒng)置信區(qū)間

    就并非θ2/θ1的最短置信區(qū)間.其中

    為了求θ2/θ1的最短置信區(qū)間,假設存在φ1,φ2滿足0<φ1<φ2,使得:

    此時θ2/θ1的置信區(qū)間為:

    這個區(qū)間的平均長度為:

    當r1=2時,F(xiàn)(2(r1-1),2(r2-1))分布的第一自由度≤2,其密度函數(shù)為單調遞減,θ2/θ1的置信度為1-α的最短置信區(qū)間為:

    成立,其中g2(x)為F(2(r1-1),2(r2-1))分布的概率密度函數(shù),那么

    (9)

    為θ2/θ1的最短置信區(qū)間.

    3 算例分析

    3.1 隨機模擬

    利用Matlab軟件可以用隨機模擬的方法產(chǎn)生一個總體為X~Pareto(θ,σ,μ)的截尾樣本,具體步驟如下:

    1)產(chǎn)生一個容量為n=30且服從均勻分布U(0,1)的獨立同分布樣本T1,T2,…,Tn;

    2)當σ=b時,給定參數(shù)θ=a和μ=c,其中a,b,c均為不等于零的常數(shù).

    3.2 單總體定數(shù)截尾下對比樞軸量法和輪廓似然函數(shù)法所求得σ的最優(yōu)置信區(qū)間

    假設產(chǎn)品的壽命服從形狀參數(shù)為θ1=2.22、尺度參數(shù)為σ1=68和位置參數(shù)為μ1=10的pareto分布,現(xiàn)在從這批產(chǎn)品中隨機抽取n1個樣品進行定數(shù)截尾壽命試驗.當n1個樣品中出現(xiàn)第r1個失效時,所得到的壽命數(shù)據(jù)(小時)按照從小到大的順序排列如下:

    當r1=1時,X1=79.1231;

    當r1=2時,X1=79.1231,X2=81.2175;

    當r1=3時,X1=79.1231,X2=81.2175,X3=82.2896;

    當r1=30時,X1=79.1231,X2=81.2175,X3=82.2896,X4=82.8524,X5=83.4619,X6=84.0014,X7=88.7709,X8=95.0987,X9=97.0298,X10=101.7209,X11=107.1334,X12=116.7245,X13=119.8929,X14=119.9307,X15=123.4077,X16=135.4312,X17=138.7833,X18=148.0011,X19=150.4866,X20=155.3181,X21=169.4068,X22=207.0749,X23=214.668,X24=217.2301,X25=241.3268,X26=291.0766,X27=292.0871,X28=298.2338,X29=317.4058,X30=342.9619.

    給定置信水平為1-α=0.95,以下通過算例分析給出尺度參數(shù)的最優(yōu)置信區(qū)間.具體算例結果如表1所示,其中置信區(qū)間1表示運用樞軸量法求得的尺度參數(shù)σ1的最短置信區(qū)間;置信區(qū)間2表示運用輪廓似然函數(shù)法求得的尺度參數(shù)σ1的置信區(qū)間;相對縮短比率ε%表示置信區(qū)間2的長度比置信區(qū)間1的長度相對縮短百分比.

    表1 尺度參數(shù)的置信區(qū)間的精度分析

    從表1中可以看出:對于定數(shù)截尾三參數(shù)pareto分布,當2≤r≤7和r≥14時,輪廓似然函數(shù)法所求得σ1的置信區(qū)間比樞軸量法求得的最短置信區(qū)間更優(yōu);當8≤r≤13時,樞軸量法所求得σ1的最短置信區(qū)間比輪廓似然函數(shù)法求得的置信區(qū)間更優(yōu).

    3.3 兩總體定數(shù)截尾下形狀參數(shù)比θ2/θ1的最優(yōu)置信區(qū)間

    在3.1的基礎上再進行壽命試驗,假設產(chǎn)品的壽命服從形狀參數(shù)為θ2=4.20,尺度參數(shù)為σ2=128和位置參數(shù)為μ2=25的pareto分布,現(xiàn)在從這批產(chǎn)品中隨機抽取n2=30個產(chǎn)品進行壽命試驗.當n2個樣品中出現(xiàn)第r2個失效時所得到的壽命數(shù)據(jù)(小時)按照從小到大的順序排列如下:

    當r2=1時,X1=153.9897;

    當r2=2時,X1=153.9897,X2=154.0728;

    當r2=3時,X1=153.9897,X2=154.0728,X3=154.4488;

    當r2=30時,X1=153.9897,X2=154.0728,X3=154.4488,X4=156.1522,X5=156.9200,X6=158.5244,X7=159.4624,X8=160.9591,X9=163.2322,X10=163.2730,X11=165.1679,X12=166.3306,X13=168.5161,X14=171.8702,X15=172.2997,X16=175.2410,X17=175.8502,X18=182.8402,X19=188.9387,X20=189.9235,X21=192.8554,X22=194.8000,X23=196.3208,X24=196.7846,X25=203.8605,X26=205.7937,X27=211.7149,X28=218.4335,X29=286.4805,X30=300.0382.

    給定置信水平為1-α=0.95,以下通過算例分析給出兩總體形狀參數(shù)比的最優(yōu)置信區(qū)間.具體算例結果如下表所示:

    因為,對于F(2(r1-1),2(r2-1))分布,當?shù)谝蛔杂啥萺1≥3時,其密度函數(shù)才為單峰非對稱.由于篇幅的問題,這里只考慮r1=3,4,5,6,7,11,r2=4,5,6,7,11時,兩總體形狀參數(shù)比的置信區(qū)間.

    表2 形狀參數(shù)比θ2/θ1置信區(qū)間的精度分析

    從表2可以看出固定r1的取值時,隨著r2的取值不斷增大,最短置信區(qū)間下邊界值不斷增大(φ1的取值不斷增大),最短置信區(qū)間上邊界值不斷減小(φ2的取值不斷減小),形狀參數(shù)比的置信區(qū)間的精度越來越高;同時還可以看出,兩種方法算得的置信區(qū)間的長度隨著r2的增大差異越來越小,縮小的速度也越來越慢.因此,在樣本容量較小時,形狀參數(shù)比的最短置信區(qū)間即為最優(yōu)置信區(qū)間.

    4 結論

    利用“樞軸量”和“輪廓似然函數(shù)”兩種方法,求出了定數(shù)截尾場合下三參數(shù)pareto分布單總體尺度參數(shù)σ的置信區(qū)間,通過算例分析的結果,得到當2≤r≤7和r≥14時,輪廓似然函數(shù)法所求得σ的置信區(qū)間比樞軸量法求得的置信區(qū)間更短,所以此時用輪廓似然函數(shù)法求得的置信區(qū)間即為最優(yōu)置信區(qū)間;當8≤r≤13時,樞軸量法求得的置信區(qū)間長度更短,這時用樞軸量法求得的置信區(qū)間即為最優(yōu)置信區(qū)間.然后,根據(jù)樞軸量法得到了兩總體形狀參數(shù)比θ2/θ1的最優(yōu)置信區(qū)間,由于容易構造形狀參數(shù)比θ2/θ1置信區(qū)間所需的樞軸量,且樞軸量所服從分布的密度函數(shù)曲線是單峰非對稱,在樣本容量較小時,用最短置信區(qū)間作為目標函數(shù)得到的置信區(qū)間比按傳統(tǒng)概率對稱方法求得的置信區(qū)間更優(yōu),參數(shù)估計精度也得到顯著提高.

    猜你喜歡
    樞軸定數(shù)置信區(qū)間
    WK-35 電鏟中央樞軸液氮冷裝工藝研究
    面向神經(jīng)機器翻譯的樞軸方法研究綜述
    探討參數(shù)區(qū)間估計中樞軸量的選取——以單個正態(tài)總體均值為例
    雙定數(shù)混合截尾下兩參數(shù)Pareto分布的統(tǒng)計分析
    Maxwell分布參數(shù)的最短置信區(qū)間研究
    p-范分布中參數(shù)的置信區(qū)間
    多個偏正態(tài)總體共同位置參數(shù)的Bootstrap置信區(qū)間
    食有定數(shù)
    意林(2019年9期)2019-06-02 16:59:40
    這邊日落
    世界文化(2019年2期)2019-03-27 16:30:18
    列車定位中置信區(qū)間的確定方法
    国产私拍福利视频在线观看| 久久精品国产清高在天天线| 91在线精品国自产拍蜜月| 春色校园在线视频观看| 婷婷色av中文字幕| 国产精品人妻久久久久久| 两个人的视频大全免费| 亚洲av成人精品一区久久| 黄色配什么色好看| 亚洲国产精品sss在线观看| 日韩欧美三级三区| 欧美色视频一区免费| 久久人妻av系列| 亚洲国产精品合色在线| 国产高清激情床上av| 校园人妻丝袜中文字幕| 精品日产1卡2卡| 国产激情偷乱视频一区二区| 男女啪啪激烈高潮av片| 日韩欧美一区二区三区在线观看| 亚洲精品久久国产高清桃花| 亚洲图色成人| 国产精品久久久久久久电影| 亚洲18禁久久av| 能在线免费看毛片的网站| 中文字幕免费在线视频6| 伦精品一区二区三区| 此物有八面人人有两片| 欧美性感艳星| 99久国产av精品| 国产伦一二天堂av在线观看| av在线观看视频网站免费| 久久久久久国产a免费观看| 在现免费观看毛片| 黄色欧美视频在线观看| 边亲边吃奶的免费视频| 淫秽高清视频在线观看| 成年免费大片在线观看| 亚洲欧美日韩东京热| 欧美+亚洲+日韩+国产| 神马国产精品三级电影在线观看| 日本欧美国产在线视频| 97热精品久久久久久| 国产成人影院久久av| 亚洲五月天丁香| 在线观看午夜福利视频| 精品欧美国产一区二区三| 人妻少妇偷人精品九色| 简卡轻食公司| 久久久久性生活片| 嘟嘟电影网在线观看| 色尼玛亚洲综合影院| 中文在线观看免费www的网站| 久久这里有精品视频免费| 人妻制服诱惑在线中文字幕| 12—13女人毛片做爰片一| 97超视频在线观看视频| 国产精品一区www在线观看| 免费人成在线观看视频色| 亚洲成人久久性| 国产一区二区在线av高清观看| 99在线视频只有这里精品首页| 成人综合一区亚洲| 国产av在哪里看| 欧美性感艳星| 国内久久婷婷六月综合欲色啪| 一个人看的www免费观看视频| 男女做爰动态图高潮gif福利片| 色播亚洲综合网| 一卡2卡三卡四卡精品乱码亚洲| 久久鲁丝午夜福利片| 少妇丰满av| 天天躁日日操中文字幕| 久久人妻av系列| 午夜福利在线观看免费完整高清在 | 高清在线视频一区二区三区 | 国产午夜福利久久久久久| 久久久久久久久久成人| 亚洲最大成人中文| 夫妻性生交免费视频一级片| 高清午夜精品一区二区三区 | 91久久精品国产一区二区三区| 最近中文字幕高清免费大全6| 小说图片视频综合网站| 国产老妇女一区| 99久久久亚洲精品蜜臀av| 只有这里有精品99| 久久精品国产亚洲av香蕉五月| 国产高清视频在线观看网站| 嫩草影院新地址| 精品欧美国产一区二区三| 免费看av在线观看网站| av天堂中文字幕网| 天堂√8在线中文| 久99久视频精品免费| 中文字幕av在线有码专区| 男人舔奶头视频| 欧美性猛交╳xxx乱大交人| 男人狂女人下面高潮的视频| 99视频精品全部免费 在线| 97超视频在线观看视频| 能在线免费观看的黄片| 国产精品乱码一区二三区的特点| 国产精品,欧美在线| 99久国产av精品国产电影| 中文字幕久久专区| 欧美又色又爽又黄视频| 夜夜爽天天搞| 亚洲丝袜综合中文字幕| 成人三级黄色视频| 午夜福利视频1000在线观看| 亚洲欧美成人精品一区二区| 日韩中字成人| 99久国产av精品国产电影| a级一级毛片免费在线观看| 久久久午夜欧美精品| 免费观看人在逋| 国产精品久久电影中文字幕| 国产高潮美女av| 最近视频中文字幕2019在线8| 女人被狂操c到高潮| 国产精品国产高清国产av| 免费av毛片视频| 国内精品宾馆在线| 免费一级毛片在线播放高清视频| av国产免费在线观看| 色综合亚洲欧美另类图片| 一边摸一边抽搐一进一小说| 小蜜桃在线观看免费完整版高清| 男人舔女人下体高潮全视频| 国产淫片久久久久久久久| 看黄色毛片网站| 日本色播在线视频| 精品久久久久久久末码| 青春草视频在线免费观看| 一进一出抽搐gif免费好疼| 亚洲综合色惰| 亚洲av第一区精品v没综合| av视频在线观看入口| 国产精品嫩草影院av在线观看| 久久精品国产99精品国产亚洲性色| 亚洲av电影不卡..在线观看| 一本精品99久久精品77| 在线国产一区二区在线| 国产熟女欧美一区二区| 久久99精品国语久久久| 日韩三级伦理在线观看| 在线观看美女被高潮喷水网站| 嫩草影院入口| 久久久久性生活片| 天堂影院成人在线观看| 天天一区二区日本电影三级| 丰满的人妻完整版| 99热这里只有是精品在线观看| 国产毛片a区久久久久| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 黄色配什么色好看| 国产国拍精品亚洲av在线观看| 一边摸一边抽搐一进一小说| 欧美性感艳星| 少妇人妻一区二区三区视频| 色吧在线观看| 一夜夜www| 国产片特级美女逼逼视频| 成熟少妇高潮喷水视频| 国内精品久久久久精免费| 国产毛片a区久久久久| 欧美又色又爽又黄视频| 国产私拍福利视频在线观看| 亚洲av熟女| 一级黄片播放器| 中文字幕久久专区| 精品一区二区三区人妻视频| 中文精品一卡2卡3卡4更新| 欧美xxxx性猛交bbbb| 亚洲最大成人中文| 久久久久久久午夜电影| 天堂影院成人在线观看| 成人三级黄色视频| 最近2019中文字幕mv第一页| 日本成人三级电影网站| 久久综合国产亚洲精品| 啦啦啦韩国在线观看视频| 深夜a级毛片| 亚洲欧美清纯卡通| 国产午夜精品一二区理论片| 丰满人妻一区二区三区视频av| 99热全是精品| 国产极品精品免费视频能看的| 美女xxoo啪啪120秒动态图| 亚洲国产精品国产精品| 91aial.com中文字幕在线观看| 成人无遮挡网站| 国产在线男女| 国产成年人精品一区二区| 日本撒尿小便嘘嘘汇集6| .国产精品久久| 国产精品三级大全| 国产色爽女视频免费观看| 只有这里有精品99| 亚洲欧美中文字幕日韩二区| 日本成人三级电影网站| 噜噜噜噜噜久久久久久91| 男人狂女人下面高潮的视频| 日本三级黄在线观看| 卡戴珊不雅视频在线播放| 晚上一个人看的免费电影| 三级经典国产精品| av免费在线看不卡| 最后的刺客免费高清国语| 国产久久久一区二区三区| 尾随美女入室| 三级国产精品欧美在线观看| 少妇丰满av| 亚洲精品国产成人久久av| 久久久久网色| 久久人妻av系列| 三级经典国产精品| 色噜噜av男人的天堂激情| 激情 狠狠 欧美| 可以在线观看的亚洲视频| 观看免费一级毛片| 一个人看的www免费观看视频| 自拍偷自拍亚洲精品老妇| 亚洲最大成人av| 国产中年淑女户外野战色| 国产成人午夜福利电影在线观看| 日本五十路高清| videossex国产| 日韩成人伦理影院| 日韩欧美三级三区| 久久精品国产亚洲av天美| 久久精品人妻少妇| 男人和女人高潮做爰伦理| 岛国毛片在线播放| 在线观看美女被高潮喷水网站| 亚洲在线观看片| 3wmmmm亚洲av在线观看| 久久99热6这里只有精品| 亚洲国产高清在线一区二区三| 国产黄色小视频在线观看| 大又大粗又爽又黄少妇毛片口| 精品久久久久久久久亚洲| 简卡轻食公司| 国产成人精品久久久久久| 91在线精品国自产拍蜜月| 久久这里有精品视频免费| 99久久精品热视频| 久久99热6这里只有精品| 寂寞人妻少妇视频99o| av视频在线观看入口| 精品久久久久久久人妻蜜臀av| 日韩欧美一区二区三区在线观看| 国产单亲对白刺激| 国产av在哪里看| 精品欧美国产一区二区三| 白带黄色成豆腐渣| 99九九线精品视频在线观看视频| 日韩一区二区视频免费看| 九九久久精品国产亚洲av麻豆| 桃色一区二区三区在线观看| 22中文网久久字幕| 最近手机中文字幕大全| 黄色一级大片看看| 一级毛片久久久久久久久女| 在线播放无遮挡| av又黄又爽大尺度在线免费看 | 成年免费大片在线观看| 大型黄色视频在线免费观看| 1024手机看黄色片| 国产午夜福利久久久久久| 99九九线精品视频在线观看视频| 久久精品国产亚洲网站| 青春草视频在线免费观看| 国产熟女欧美一区二区| 午夜a级毛片| 老司机影院成人| h日本视频在线播放| 亚洲最大成人av| 成人美女网站在线观看视频| 国产私拍福利视频在线观看| 国产乱人偷精品视频| 国产 一区精品| 三级毛片av免费| 亚洲中文字幕日韩| 成人无遮挡网站| 看免费成人av毛片| 欧美日韩乱码在线| 久久久久免费精品人妻一区二区| 激情 狠狠 欧美| 精品国内亚洲2022精品成人| 亚洲18禁久久av| 尤物成人国产欧美一区二区三区| 看非洲黑人一级黄片| 国产精品一二三区在线看| 国产精品综合久久久久久久免费| 1024手机看黄色片| 国产精品一区二区性色av| 欧美成人一区二区免费高清观看| 大香蕉久久网| 亚洲av熟女| 色5月婷婷丁香| 麻豆国产97在线/欧美| 老师上课跳d突然被开到最大视频| 日本色播在线视频| 99久国产av精品| 波多野结衣高清无吗| 精品一区二区三区视频在线| 色哟哟哟哟哟哟| 夜夜爽天天搞| 熟妇人妻久久中文字幕3abv| 热99在线观看视频| 国产淫片久久久久久久久| 女人十人毛片免费观看3o分钟| 麻豆成人午夜福利视频| 亚洲三级黄色毛片| 男的添女的下面高潮视频| 狂野欧美激情性xxxx在线观看| 亚洲一区高清亚洲精品| 人人妻人人澡人人爽人人夜夜 | 精品99又大又爽又粗少妇毛片| 寂寞人妻少妇视频99o| 超碰av人人做人人爽久久| 国产精品久久电影中文字幕| 亚洲国产精品成人久久小说 | 一夜夜www| 国产色爽女视频免费观看| 不卡一级毛片| 听说在线观看完整版免费高清| 日韩大尺度精品在线看网址| 少妇的逼好多水| 国产日韩欧美在线精品| 亚洲av成人av| 中文字幕久久专区| 有码 亚洲区| 中文字幕久久专区| 色吧在线观看| 在线免费十八禁| 日本av手机在线免费观看| 免费在线观看成人毛片| 欧美潮喷喷水| 此物有八面人人有两片| 欧美zozozo另类| 亚洲美女搞黄在线观看| 国产精品综合久久久久久久免费| 狠狠狠狠99中文字幕| 久久久精品大字幕| 亚洲自拍偷在线| 日本欧美国产在线视频| 国产伦在线观看视频一区| 亚洲综合色惰| 久久久久性生活片| 国产精品久久久久久久电影| 久久久精品欧美日韩精品| 色尼玛亚洲综合影院| 波多野结衣高清无吗| 亚洲va在线va天堂va国产| 国产成人一区二区在线| 亚洲国产欧洲综合997久久,| 午夜福利成人在线免费观看| 一区二区三区高清视频在线| 日韩成人av中文字幕在线观看| 国产亚洲av片在线观看秒播厂 | 免费不卡的大黄色大毛片视频在线观看 | 久久久久久久久大av| 97热精品久久久久久| 18+在线观看网站| 色综合亚洲欧美另类图片| 国产精品美女特级片免费视频播放器| 久久久久国产网址| 久久人人爽人人爽人人片va| 联通29元200g的流量卡| 亚洲欧洲日产国产| 中国美白少妇内射xxxbb| 成人漫画全彩无遮挡| 欧美日本视频| 日日干狠狠操夜夜爽| 欧美激情久久久久久爽电影| 国模一区二区三区四区视频| 日本-黄色视频高清免费观看| 国产亚洲av嫩草精品影院| 99视频精品全部免费 在线| 最近手机中文字幕大全| 99久久久亚洲精品蜜臀av| 美女cb高潮喷水在线观看| 一个人看视频在线观看www免费| 三级国产精品欧美在线观看| 三级国产精品欧美在线观看| 欧美一区二区亚洲| 校园春色视频在线观看| 欧美又色又爽又黄视频| 亚洲精品久久国产高清桃花| 国产激情偷乱视频一区二区| 久久精品久久久久久噜噜老黄 | 国产精品久久视频播放| 国产老妇女一区| 亚洲欧美清纯卡通| 国产女主播在线喷水免费视频网站 | 99热网站在线观看| 午夜精品一区二区三区免费看| 国产亚洲91精品色在线| kizo精华| av视频在线观看入口| 网址你懂的国产日韩在线| 久久午夜福利片| 欧美日韩一区二区视频在线观看视频在线 | 婷婷六月久久综合丁香| 亚洲欧美日韩高清专用| 亚洲国产日韩欧美精品在线观看| 色吧在线观看| 成年av动漫网址| 欧美日韩综合久久久久久| 日本爱情动作片www.在线观看| 波多野结衣巨乳人妻| 久久久久久久久大av| 亚洲久久久久久中文字幕| 国产精品无大码| 久久99热6这里只有精品| 日韩一本色道免费dvd| 能在线免费看毛片的网站| 亚洲国产精品合色在线| 中文字幕熟女人妻在线| 国产成人精品一,二区 | 婷婷色av中文字幕| 久久人人精品亚洲av| 国产黄色视频一区二区在线观看 | 亚洲精品乱码久久久久久按摩| 性插视频无遮挡在线免费观看| 97超碰精品成人国产| 天堂√8在线中文| 亚洲天堂国产精品一区在线| 26uuu在线亚洲综合色| 精品欧美国产一区二区三| 亚洲成人av在线免费| 综合色av麻豆| 日韩欧美国产在线观看| 九九久久精品国产亚洲av麻豆| 久久久国产成人精品二区| 只有这里有精品99| 精品无人区乱码1区二区| 悠悠久久av| 在现免费观看毛片| 永久网站在线| 亚洲国产精品成人综合色| 不卡视频在线观看欧美| 少妇裸体淫交视频免费看高清| 国产成人午夜福利电影在线观看| 国产成人精品婷婷| 国产综合懂色| 99久久中文字幕三级久久日本| 久久久久久久久久久丰满| 深夜a级毛片| 久久久久性生活片| 亚洲一区二区三区色噜噜| 国产午夜精品一二区理论片| 在线免费十八禁| 不卡一级毛片| 精品一区二区三区视频在线| 亚洲第一区二区三区不卡| 国产视频首页在线观看| 亚洲婷婷狠狠爱综合网| 国产成人freesex在线| av在线天堂中文字幕| 精品国内亚洲2022精品成人| 国产精品,欧美在线| 午夜精品国产一区二区电影 | 久久久久免费精品人妻一区二区| 精品久久久久久久久av| 精品久久久噜噜| 国产精品国产三级国产av玫瑰| 级片在线观看| 亚洲一级一片aⅴ在线观看| 啦啦啦韩国在线观看视频| 美女xxoo啪啪120秒动态图| 网址你懂的国产日韩在线| 在现免费观看毛片| 亚洲自拍偷在线| 精品久久久久久久久久久久久| АⅤ资源中文在线天堂| 搞女人的毛片| 一级毛片我不卡| 22中文网久久字幕| 欧美另类亚洲清纯唯美| 深夜精品福利| 欧美人与善性xxx| 日韩高清综合在线| 亚洲精品久久久久久婷婷小说 | 大又大粗又爽又黄少妇毛片口| 一级二级三级毛片免费看| 国产精品一区二区在线观看99 | 国产视频首页在线观看| h日本视频在线播放| 麻豆成人午夜福利视频| 亚洲精品成人久久久久久| 国产 一区 欧美 日韩| 成年女人永久免费观看视频| 黄色一级大片看看| av在线播放精品| 色综合色国产| 亚洲国产高清在线一区二区三| 天堂√8在线中文| 午夜久久久久精精品| 精品久久久久久久久亚洲| 女人被狂操c到高潮| 国产真实伦视频高清在线观看| 女人被狂操c到高潮| 成人三级黄色视频| 国产精品伦人一区二区| 免费观看a级毛片全部| 桃色一区二区三区在线观看| 精品国产三级普通话版| 99热精品在线国产| 亚洲av熟女| 成人毛片60女人毛片免费| 日本撒尿小便嘘嘘汇集6| 免费av毛片视频| 国产成人午夜福利电影在线观看| 久久午夜福利片| 国产三级中文精品| 欧洲精品卡2卡3卡4卡5卡区| 色综合站精品国产| 黄色日韩在线| 精品无人区乱码1区二区| 菩萨蛮人人尽说江南好唐韦庄 | 黄色欧美视频在线观看| 搡老妇女老女人老熟妇| 我的老师免费观看完整版| 婷婷色av中文字幕| 波多野结衣高清作品| 国产精品乱码一区二三区的特点| 国产精品嫩草影院av在线观看| 亚洲精品日韩av片在线观看| 中文字幕久久专区| 亚洲aⅴ乱码一区二区在线播放| av在线蜜桃| 成熟少妇高潮喷水视频| 亚洲国产欧美人成| 夫妻性生交免费视频一级片| 久久精品国产亚洲网站| 免费电影在线观看免费观看| 成人永久免费在线观看视频| 午夜福利视频1000在线观看| 亚洲国产精品sss在线观看| 午夜精品一区二区三区免费看| 69人妻影院| 边亲边吃奶的免费视频| 欧美高清性xxxxhd video| 久久久精品94久久精品| 女同久久另类99精品国产91| 国产亚洲欧美98| 爱豆传媒免费全集在线观看| 国产探花极品一区二区| av女优亚洲男人天堂| 男人舔女人下体高潮全视频| 小蜜桃在线观看免费完整版高清| 中国美女看黄片| 2021天堂中文幕一二区在线观| 国产精品麻豆人妻色哟哟久久 | 亚洲精品国产成人久久av| 一级黄色大片毛片| 欧美3d第一页| 久久精品影院6| 久久久久久久久久久免费av| 欧美xxxx性猛交bbbb| 亚洲av男天堂| 久久精品国产亚洲av涩爱 | 18禁黄网站禁片免费观看直播| 99精品在免费线老司机午夜| 亚洲五月天丁香| av在线观看视频网站免费| 99热这里只有是精品50| 日日干狠狠操夜夜爽| 日韩人妻高清精品专区| 高清在线视频一区二区三区 | 日本爱情动作片www.在线观看| 一个人看的www免费观看视频| 深夜精品福利| 中文字幕精品亚洲无线码一区| 国产午夜精品久久久久久一区二区三区| 男人狂女人下面高潮的视频| 久久精品国产亚洲av香蕉五月| 老女人水多毛片| 亚州av有码| 一级毛片aaaaaa免费看小| 中国国产av一级| 伦理电影大哥的女人| 一本久久中文字幕| 精品人妻偷拍中文字幕| 成人二区视频| 亚洲综合色惰| 久久这里有精品视频免费| 午夜亚洲福利在线播放| 国产精品一及| 特级一级黄色大片| 熟女人妻精品中文字幕| 日日撸夜夜添| 久久鲁丝午夜福利片| 最近的中文字幕免费完整| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 日韩一本色道免费dvd| 中文字幕人妻熟人妻熟丝袜美| 日韩成人av中文字幕在线观看| 国内少妇人妻偷人精品xxx网站| 天堂影院成人在线观看| 91在线精品国自产拍蜜月| 男人舔奶头视频| 亚洲在线观看片| 亚洲高清免费不卡视频| 国产真实伦视频高清在线观看| 国产激情偷乱视频一区二区| 亚洲欧美清纯卡通| 男女下面进入的视频免费午夜| 可以在线观看毛片的网站| 国产精品久久视频播放| 成人三级黄色视频|