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

    基于MATLAB的De Moivre-Laplace中心極限定理的隨機(jī)模擬

    2016-09-09 01:32:28任麗李順初
    關(guān)鍵詞:二項(xiàng)分布正態(tài)分布學(xué)報(bào)

    任麗,李順初

    (西華大學(xué)理學(xué)院應(yīng)用數(shù)學(xué)研究所,成都610039)

    ?

    基于MATLAB的De Moivre-Laplace中心極限定理的隨機(jī)模擬

    任麗,李順初

    (西華大學(xué)理學(xué)院應(yīng)用數(shù)學(xué)研究所,成都610039)

    針對(duì)De Moivre-Laplace中心極限定理,首先利用特征函數(shù)和獨(dú)立同分布下的中心極限定理得到理論上證明,然后利用MATLAB軟件隨機(jī)模擬分析二項(xiàng)分布與正態(tài)分布之間的關(guān)系,并對(duì)圖像進(jìn)行誤差分析,最后證明數(shù)值結(jié)果與實(shí)驗(yàn)結(jié)果的一致性。通過這兩種不同角度的論證,給出De Moivre-Laplace中心極限定理直觀的解釋和說明,這種將數(shù)學(xué)實(shí)驗(yàn)與數(shù)學(xué)原理結(jié)合的方法,不僅使數(shù)學(xué)原理化抽象為具體,便于理解和實(shí)踐,而且擴(kuò)展了數(shù)學(xué)原理和計(jì)算機(jī)軟件的應(yīng)用。

    De Moivre-Laplace中心極限定理;正態(tài)分布;二項(xiàng)分布;隨機(jī)模擬;誤差分析

    引言

    中心極限定理描述了離散型隨機(jī)變量與連續(xù)型隨機(jī)變量的內(nèi)在聯(lián)系,即離散型隨機(jī)變量的極限分布是正態(tài)分布,因此中心極限定理的結(jié)論使正態(tài)分布在概率論和數(shù)理統(tǒng)計(jì)中具有重要作用[1],它的研究促進(jìn)了現(xiàn)代概率論的發(fā)展,并且在統(tǒng)計(jì)分析和近似計(jì)算具有廣泛的應(yīng)用。近年來,國內(nèi)眾多學(xué)者從兩方面對(duì)中心極限定理進(jìn)行了研究:一方面,如文獻(xiàn)[2-9]研究了中心極限定理的實(shí)際應(yīng)用;另一方面,如文獻(xiàn)[10-14]研究了中心極限定理含義和本質(zhì),其研究?jī)?nèi)容是單一的,只用了數(shù)值方法證明中心極限定理,并沒有詳細(xì)闡明常見分布(幾種離散分布和連續(xù)分布)與正態(tài)分布的關(guān)系,也沒有用數(shù)學(xué)軟件模擬這種關(guān)系。

    本文以De Moivre-Laplace中心極限定理為例,在蒙特卡羅思想[15]下,研究中心極限定理中二項(xiàng)分布與正態(tài)分布的關(guān)系,并通過MATLAB軟件[16-20]對(duì)其進(jìn)行模擬,直觀展現(xiàn)De Moivre-Laplace中心極限定理的形成過程,并驗(yàn)證數(shù)值結(jié)果的準(zhǔn)確性,同時(shí)對(duì)實(shí)驗(yàn)結(jié)果進(jìn)行誤差分析,進(jìn)一步得到影響二項(xiàng)分布近似正態(tài)分布的因子,從而對(duì)De Moivre-Laplace中心極限定理作出更合理地詮釋,將為中心極限定理在其他分布中的應(yīng)用提供了一種新的分析手段。這種將計(jì)算機(jī)和數(shù)學(xué)原理結(jié)合的方法,不僅將數(shù)學(xué)原理具體化,而且還實(shí)現(xiàn)了數(shù)學(xué)實(shí)驗(yàn)中計(jì)算機(jī)軟件的應(yīng)用與擴(kuò)展。

    1 預(yù)備知識(shí)

    引理1[21]設(shè){xn}是獨(dú)立同分布的隨機(jī)變量序列,且E(Xi)=μ,Var(Xi)=σ2>0,若記

    則對(duì)任意實(shí)數(shù)y,有

    (1)

    因?yàn)?/p>

    E(Xi)=0,Var(Xi)=σ2

    所以有

    Var(Xn-μ)=σ2=-φ″(0)+φ″(0)2

    由特征函數(shù)的性質(zhì)可得:

    φ′(0)=0,φ″(0)=-σ2

    又因?yàn)樘卣骱瘮?shù)有展開式

    故(1)式得證,引理1證明完畢。

    只假設(shè){Xn}獨(dú)立同分布、方差存在,不管原來的分布是什么,只要當(dāng)n充分大時(shí),Yn才近似服從標(biāo)準(zhǔn)正態(tài)分布N(0,1),當(dāng)n較小時(shí),此種近似不能保證。所以這個(gè)定理有著廣泛的應(yīng)用[21]。

    2 De Moivre-Laplace中心極限定理證明及模擬算法

    2.1De Moivre-Laplace中心極限定理

    設(shè)n重貝努里試驗(yàn)中,事件A在每次實(shí)驗(yàn)中出現(xiàn)的概率為p(0

    Xi~b(1,p)E(Xi)=0,Var(Xi)=σ2,i=1,2,…,n

    且記

    則對(duì)任意實(shí)數(shù)y,有

    (2)

    證明存在n個(gè)相互獨(dú)立,服從同一0-1分布的諸隨機(jī)變量X1,X2,X3,…,Xn,其中Xi(i=1,2,…,n)的分布概率為P(Xi=i)=pi(1-p)1-i,i=0,1,因?yàn)?/p>

    E(Xi)=p,Var(Xi)=p(1-p)

    i=1,2,…,n

    由林德伯格-萊維中心極限定理可得:

    2.2De Moivre-Laplace中心極限定理的隨機(jī)模擬算法設(shè)計(jì)

    (1)從計(jì)算機(jī)中產(chǎn)生100個(gè)b(100,0.9)二項(xiàng)分布的隨機(jī)數(shù),記為x1x2x3,…,x100。

    (2)計(jì)算

    由De Moivre-Laplace中心極限定理知,可將y近似看成來自標(biāo)準(zhǔn)正態(tài)分布N(0,1)的一個(gè)隨機(jī)數(shù)。

    (3)計(jì)算z=μ+σy,則可將z看成來自正態(tài)分布N(μ,σ2)的一個(gè)隨機(jī)數(shù)。

    (4)重復(fù)(1)~(3)k次,就可以得N(μ,σ2)分布的k個(gè)隨機(jī)數(shù)。

    2.3De Moivre-Laplace中心極限定理的隨機(jī)模擬流程圖以及程序

    De Moivre-Laplace中心極限定理隨機(jī)模擬流程圖如圖1所示。

    圖1 De Moivre-Laplace中心極限定理隨機(jī)模擬流圖

    根據(jù)流程圖1,編寫程序如下:

    clear all;

    m=20000000; % m;運(yùn)行次數(shù)

    mu=2;

    sigma=5;

    z=zeros(1,n);

    n=100;

    p=0.9;

    k=100;%隨機(jī)數(shù)的個(gè)數(shù)

    for i=1:n

    x=binornd(n,p,1,k);

    y=(sum(x(1,:))-k*n*p)

    /sqrt(k*n*p*(1-p));

    z(i)=mu+sigma*y;

    end

    zmin=min(z);

    zmax=max(z);

    t=linspace(zmin,zmax,100);

    zz=hist(z,t);

    zz=zz/length(z)/((zmax-zmin)/100);

    plot(t,zz,′r′)

    hold on;

    xxx=-30:0.1:30;

    yyy=normpdf(xxx,mu,sigma);

    plot(xxx,yyy);

    h=legend(′二項(xiàng)分布′,′正態(tài)分布′,1);

    for i=1:k

    e(i)=abs(yyy(i)-zz(i));% 絕對(duì)值誤差

    end

    E1=max(e)% 最大模誤差

    E2=norm(e,2)%平方誤差

    在參量n=100,p=0.9下m=200 000,m=2 000 000,m=2 500 000的運(yùn)行結(jié)果如圖2所示。

    圖2 不同運(yùn)行次數(shù)下的運(yùn)行結(jié)果

    2.4De Moivre-Laplace中心極限定理隨機(jī)模擬圖象誤差分析

    要實(shí)現(xiàn)De Moivre-Laplace中心極限定理的模擬,關(guān)鍵在于盡可能多地進(jìn)行抽樣試驗(yàn),為了解決這一問題,本文在MATLAB R2013a軟件環(huán)境下產(chǎn)生服從二項(xiàng)分布的隨機(jī)數(shù)代替現(xiàn)實(shí)中的抽樣試驗(yàn),由De Moivre-Laplace中心極限定理可知:對(duì)其隨機(jī)模擬的運(yùn)行結(jié)果產(chǎn)生影響的變量只有三個(gè),即運(yùn)行次數(shù)m,二項(xiàng)分布中n和p。對(duì)此,從變量m、n、p的取值角度出發(fā),討論對(duì)De Moivre-Laplace中心極限定理的隨機(jī)模擬圖象誤差分析,當(dāng)n=100,p=0.9時(shí)改變模擬次數(shù)m,隨機(jī)模擬De Moivre-Laplace中心極限定理數(shù)值分析比較見表1。

    表1 不同m下的誤差

    由圖2可知,隨著運(yùn)行次數(shù)的增大,De Moivre-Laplace中心極限定理隨機(jī)模擬圖像效果愈加顯著,二項(xiàng)分布和標(biāo)準(zhǔn)正態(tài)分布重合愈緊密;由表1的二項(xiàng)分布與標(biāo)準(zhǔn)正態(tài)分布圖形之間誤差分析也說明了圖2結(jié)果。

    隨機(jī)模擬運(yùn)行次數(shù)m=2×106,在n取值一定(或p取值一定),改變p的取值(或改變n的取值)隨機(jī)模擬De Moivre-Laplace中心極限定理的數(shù)值分析見表2。

    表2 隨n和p變化的誤差分析

    由表2給出的計(jì)算誤差可知,在隨機(jī)模擬次數(shù)m和二項(xiàng)分布某個(gè)事件發(fā)生概率p一定的情況下,增加進(jìn)行二項(xiàng)分布樣本點(diǎn)的總個(gè)數(shù)n,可減少De Moivre-Laplace中心極限定理隨機(jī)模擬圖像的誤差;同理,在隨機(jī)模擬次數(shù)m和進(jìn)行二項(xiàng)分布樣本點(diǎn)的總個(gè)數(shù)n一定的情況下,增加二項(xiàng)分布某個(gè)事件發(fā)生概率p,可減De Moivre-Laplace中心極限定理隨機(jī)模擬圖像的誤差。但是,也出現(xiàn)某些點(diǎn)的誤差隨著二項(xiàng)分布樣本點(diǎn)的總個(gè)數(shù)n和二項(xiàng)分布某個(gè)事件發(fā)生概率p的增大反而略有變大的情況,這主要是由于蒙特卡羅方法在本質(zhì)上屬于一種非確定性的方法,并且De Moivre-Laplace中心極限定理進(jìn)行的是隨機(jī)模擬,從而會(huì)顯示出一定的隨機(jī)性。

    3 結(jié)論

    (1)在分析De Moivre-Laplace中心極限定理的隨機(jī)模擬圖像后發(fā)現(xiàn):隨機(jī)模擬的運(yùn)行次數(shù)m,二項(xiàng)分布中n和p這三個(gè)因素對(duì)于二項(xiàng)分布的正態(tài)近似都有一定程度的影響。分別增加運(yùn)行次數(shù)m,二項(xiàng)分布中n和p,都可以減少二項(xiàng)分布近似正態(tài)分布的誤差,但由于蒙特卡羅方法本質(zhì)上屬于一種非確定性方法,因而在分析隨機(jī)模擬圖像某些點(diǎn)的誤差時(shí),分別增加運(yùn)行次數(shù)m、二項(xiàng)分布中n和p誤差反而增大。

    (2)基于蒙特卡羅方法的基本思想,通過產(chǎn)生的正態(tài)隨機(jī)數(shù),在MATLAB環(huán)境下隨機(jī)模擬中心極限定理,把傳統(tǒng)教學(xué)中無法實(shí)現(xiàn)的大量實(shí)驗(yàn)變成現(xiàn)實(shí),從而把原本抽象難懂的知識(shí)變得直觀具體,易于理解。在計(jì)算機(jī)的幫助下學(xué)習(xí)數(shù)學(xué)知識(shí)、分析理論、總結(jié)規(guī)律,也實(shí)現(xiàn)數(shù)學(xué)實(shí)驗(yàn)中計(jì)算機(jī)軟件的應(yīng)用與擴(kuò)展,它體現(xiàn)了將數(shù)學(xué)視為一門“試驗(yàn)科學(xué)”的思想。

    [1] 唐莉,李雁如.大數(shù)定律與中心極限定理的實(shí)際應(yīng)用[J].廣東技術(shù)師范學(xué)院學(xué)報(bào),2005(6):75-76.

    [2] 孫大利,王久珂,劉曉陽,等.獨(dú)立同分布中心極限定理在雨量站網(wǎng)規(guī)劃中的應(yīng)用[J].北京大學(xué)學(xué)報(bào):自然科學(xué)版,2015,51(1):35-42.

    [3] 丁健,李紅菊.中心極限定理在統(tǒng)計(jì)推斷中的應(yīng)用[J].長春師范大學(xué)學(xué)報(bào),2015,34(2):12-14.

    [4] 陳學(xué)慧,趙魯濤,張志剛.案例式中心極限定理教學(xué)研究[J].大學(xué)數(shù)學(xué),2015,31(2):115-118.

    [5] 孫蓓.中心極限定理及其在若干實(shí)際問題中的應(yīng)用[J].科教導(dǎo)刊,2012(6):65-67.

    [6] 荊江雁.概率方法在計(jì)算與證明中的應(yīng)用[J].常州工學(xué)院學(xué)報(bào),2011,24(2):30-31.

    [7] 徐姿奕,汪四水.獨(dú)立同分布中心極限定理的應(yīng)用[J].山西大同大學(xué)學(xué)報(bào):自然學(xué)科版,2007,23(1):13-14.

    [8] 杜偉娟,于文娟.中心極限定理及其初步應(yīng)用[J].內(nèi)蒙古電大學(xué)刊,2007(7):107-108.

    [9] 岳金健.利用大數(shù)定律和中心極限定理求解極限[J].龍巖學(xué)院學(xué)報(bào),2007,25(3):95-97.

    [10] 石峰利.幾個(gè)常見統(tǒng)計(jì)模擬問題及其在excel中的實(shí)現(xiàn)和動(dòng)態(tài)展示[D].云南:云南師范大學(xué),2013.

    [11] 蔣觀敏.Mathcad在概率統(tǒng)計(jì)模擬實(shí)驗(yàn)教學(xué)中的應(yīng)用[D].重慶:西南大學(xué),2011.

    [12] 農(nóng)吉夫,黃文寧.隨機(jī)模擬實(shí)驗(yàn)在概率極限理論教學(xué)中的應(yīng)用[J].廣西民族大學(xué)學(xué)報(bào):自然科學(xué)版,2008,14(1):101-104.

    [13] 羅弟亞.大數(shù)定理和中心極限定理可視化[J].四川工業(yè)學(xué)院學(xué)報(bào),2005(6):70-75.

    [14] 林小蘋,吳文杰.用MATLAB模擬大數(shù)定律和中心極限定理[J].運(yùn)城高專學(xué)報(bào),1990(3):12-18.

    [15] 文軍,嚴(yán)忠權(quán).基于matlab的離散事件的隨機(jī)模擬[J].黔南民族師范學(xué)院學(xué)報(bào),2011(6):17-23.

    [16] 劉澤顯,劉志偉,王勤龍.基于Matlab的概率統(tǒng)計(jì)實(shí)驗(yàn)教學(xué)的探索與實(shí)踐[J].賀州學(xué)院學(xué)報(bào),2014,30(4):122-126.

    [17] 莊光明,夏建偉,彭作祥,等.基于Matlab的Poisson分布隨機(jī)數(shù)的Monte carlo模擬[J].數(shù)學(xué)的實(shí)踐與認(rèn)識(shí),2012,42(5):87-92.

    [18] 農(nóng)吉夫,黃文寧.隨機(jī)模擬實(shí)驗(yàn)在概率極限理論教學(xué)中的應(yīng)用[J].廣西民族大學(xué)學(xué)報(bào),2008,14(1):101-104.

    [19] 張杰.MATLAB在計(jì)算物理課程教學(xué)中的應(yīng)用[J].計(jì)算機(jī)應(yīng)用與軟件,2005,22(6):131-132.

    [20] 馬鋒,封盛.在工程風(fēng)險(xiǎn)分析中應(yīng)用MATLAB進(jìn)行蒙特卡洛模擬[J].經(jīng)濟(jì)管理,2005,25(2):156-157.

    [21] 茆詩松,程依明,濮曉龍.概率論與數(shù)理統(tǒng)計(jì)教程[M].2版.北京:高等教育出版社,2011.

    Random Simulation of De Moivre-Laplace Central Limit Theorem

    REN Li,LI Shunchu

    (Institute of Applied Mathematics, College of Science, Xihua University, Chengdu 610039, China)

    Based on De Moivre-Laplace Central Limit Theorem, first, this paper gives a theoretical proof on the theorem by characteristic function and Lindeberg- Levy Theorem. Second, this paper explains and simulates the relationship between normal distribution and binomial distribution by using MATLAB software, and then analyzes the error for the image that is produced by the MATLAB software. Last, prove that the numerical result and the analysis result are consistent. By these two different angles of argument, this paper gives the intuitive explanation and illustration of Moivre-Laplace De central limit theorem. This method combined mathematical experiment and mathematical principles, not only makes mathematical principles specific, easy to understand and practice, but also extends the application of the mathematical principles and computer software.

    De Moivre-Laplace Central Limit Theorem; normal distribution; binomial distribution; random simulation; error analysis

    2016-01-03

    四川省教育廳自然科學(xué)重點(diǎn)項(xiàng)目(12ZA164)

    任麗(1991-),女,四川南充人,碩士生,主要從事微分方程及其應(yīng)用方面的研究,(E-mail)15208416930@163.com;

    李順初(1963-),男,湖北黃岡人,教授,主要從事微分方程及其應(yīng)用、油氣藏工程、滲流力學(xué)方面的研究,(E-mail)lishunchu@163.com

    1673-1549(2016)02-0080-05

    10.11863/j.suse.2016.02.16

    O021;O024

    A

    猜你喜歡
    二項(xiàng)分布正態(tài)分布學(xué)報(bào)
    二項(xiàng)分布與超幾何分布的區(qū)別與聯(lián)系
    深度剖析超幾何分布和二項(xiàng)分布
    概率與統(tǒng)計(jì)(1)——二項(xiàng)分布與超幾何分布
    致敬學(xué)報(bào)40年
    二項(xiàng)分布參數(shù)的E-Bayes估計(jì)及其應(yīng)用
    基于對(duì)數(shù)正態(tài)分布的出行時(shí)長可靠性計(jì)算
    正態(tài)分布及其應(yīng)用
    正態(tài)分布題型剖析
    χ2分布、t 分布、F 分布與正態(tài)分布間的關(guān)系
    學(xué)報(bào)簡(jiǎn)介
    x7x7x7水蜜桃| 日本欧美国产在线视频| 日日干狠狠操夜夜爽| 婷婷精品国产亚洲av| 日本黄大片高清| 国产三级在线视频| 香蕉av资源在线| 国产黄a三级三级三级人| 国产国拍精品亚洲av在线观看| 偷拍熟女少妇极品色| 成人无遮挡网站| 国产真实乱freesex| 亚洲av不卡在线观看| 国产色爽女视频免费观看| 搡老熟女国产l中国老女人| av福利片在线观看| 97碰自拍视频| 99久久精品热视频| 美女高潮喷水抽搐中文字幕| 亚洲精品色激情综合| 少妇的逼好多水| 综合色av麻豆| 丝袜美腿在线中文| 日本一二三区视频观看| 久久99热6这里只有精品| 日韩欧美在线二视频| 真实男女啪啪啪动态图| 老女人水多毛片| 18禁黄网站禁片免费观看直播| 国产熟女欧美一区二区| 内射极品少妇av片p| 亚洲性久久影院| 国产毛片a区久久久久| .国产精品久久| 国产亚洲精品久久久久久毛片| 午夜福利高清视频| 一级a爱片免费观看的视频| 日本色播在线视频| 国产av在哪里看| 最新中文字幕久久久久| 国产精品久久久久久亚洲av鲁大| 大又大粗又爽又黄少妇毛片口| 尤物成人国产欧美一区二区三区| 国产精品电影一区二区三区| 亚洲国产日韩欧美精品在线观看| 在线a可以看的网站| 国产一区二区三区视频了| 丰满人妻一区二区三区视频av| 中国美女看黄片| 亚洲美女搞黄在线观看 | 国产精品人妻久久久久久| 91午夜精品亚洲一区二区三区 | 内射极品少妇av片p| 午夜影院日韩av| 欧美一区二区亚洲| 亚洲不卡免费看| 免费av不卡在线播放| 草草在线视频免费看| 欧美成人a在线观看| 欧美高清性xxxxhd video| 婷婷亚洲欧美| 美女黄网站色视频| 亚洲欧美日韩高清专用| 人人妻人人看人人澡| 国产三级中文精品| 一夜夜www| 国产精品久久久久久久久免| 小说图片视频综合网站| 男女啪啪激烈高潮av片| www日本黄色视频网| 国产爱豆传媒在线观看| 伦理电影大哥的女人| 成人国产综合亚洲| 搞女人的毛片| 国产精品亚洲美女久久久| 中国美白少妇内射xxxbb| 国产乱人伦免费视频| 中文字幕av成人在线电影| 国产精品av视频在线免费观看| 成年女人永久免费观看视频| 3wmmmm亚洲av在线观看| av在线亚洲专区| 国语自产精品视频在线第100页| 成人毛片a级毛片在线播放| 日韩精品有码人妻一区| 亚洲狠狠婷婷综合久久图片| 欧美一区二区精品小视频在线| 免费观看在线日韩| 国产三级在线视频| 国产主播在线观看一区二区| 日韩欧美三级三区| 最近中文字幕高清免费大全6 | 国产v大片淫在线免费观看| 中文字幕久久专区| 两个人的视频大全免费| 午夜福利欧美成人| 亚洲熟妇熟女久久| 国产一级毛片七仙女欲春2| 国产精品一区二区性色av| 99国产极品粉嫩在线观看| 最新中文字幕久久久久| 日韩强制内射视频| 禁无遮挡网站| 国产精品嫩草影院av在线观看 | 国产精品亚洲一级av第二区| 噜噜噜噜噜久久久久久91| 久久久国产成人精品二区| 日韩精品中文字幕看吧| 深夜精品福利| 欧美丝袜亚洲另类 | 日本 欧美在线| 老女人水多毛片| 12—13女人毛片做爰片一| 丰满乱子伦码专区| 亚洲avbb在线观看| 动漫黄色视频在线观看| 欧美一区二区国产精品久久精品| 欧美极品一区二区三区四区| 1000部很黄的大片| 一个人免费在线观看电影| 亚洲精品456在线播放app | 色噜噜av男人的天堂激情| 国产精品不卡视频一区二区| 国产女主播在线喷水免费视频网站 | 国产成年人精品一区二区| www.www免费av| 久久国产精品人妻蜜桃| 日本黄大片高清| 十八禁网站免费在线| 日本欧美国产在线视频| 淫妇啪啪啪对白视频| 成年女人看的毛片在线观看| 99久久精品热视频| a在线观看视频网站| 别揉我奶头 嗯啊视频| 国产乱人伦免费视频| 亚洲精品成人久久久久久| 天美传媒精品一区二区| 午夜激情福利司机影院| 欧美最黄视频在线播放免费| 精品不卡国产一区二区三区| 热99在线观看视频| 久久久久国内视频| 国产高清三级在线| 国产老妇女一区| 精品人妻1区二区| 日本 av在线| 亚洲一级一片aⅴ在线观看| 狠狠狠狠99中文字幕| 国产日本99.免费观看| 99久久精品一区二区三区| 在线a可以看的网站| 香蕉av资源在线| 午夜爱爱视频在线播放| 国产毛片a区久久久久| 亚洲av免费在线观看| 99热6这里只有精品| 国内精品一区二区在线观看| 久久人人爽人人爽人人片va| 欧美另类亚洲清纯唯美| 国产精品久久久久久亚洲av鲁大| 成人无遮挡网站| 窝窝影院91人妻| 精品人妻视频免费看| 一边摸一边抽搐一进一小说| 成人精品一区二区免费| 亚洲美女视频黄频| 亚洲人成伊人成综合网2020| 亚洲人成网站在线播放欧美日韩| 少妇裸体淫交视频免费看高清| 亚洲第一电影网av| 亚洲av中文av极速乱 | 九九在线视频观看精品| 噜噜噜噜噜久久久久久91| 成熟少妇高潮喷水视频| 成人特级黄色片久久久久久久| 露出奶头的视频| 99久国产av精品| 国产老妇女一区| 亚洲内射少妇av| 一级av片app| 国产老妇女一区| 久久人人爽人人爽人人片va| 久久精品国产亚洲av天美| 色哟哟哟哟哟哟| 在线观看美女被高潮喷水网站| 在线播放国产精品三级| 婷婷精品国产亚洲av在线| 欧美zozozo另类| 成人美女网站在线观看视频| 搡女人真爽免费视频火全软件 | 国产亚洲欧美98| 成年版毛片免费区| 久久久久久大精品| 男女那种视频在线观看| 一级a爱片免费观看的视频| 亚洲最大成人手机在线| 日日摸夜夜添夜夜添av毛片 | 国产探花在线观看一区二区| 国产精品福利在线免费观看| 啪啪无遮挡十八禁网站| 国产色婷婷99| 精华霜和精华液先用哪个| 最后的刺客免费高清国语| 免费一级毛片在线播放高清视频| 精品人妻视频免费看| 免费观看的影片在线观看| 人人妻人人看人人澡| 搡女人真爽免费视频火全软件 | 最新在线观看一区二区三区| 国产黄片美女视频| 有码 亚洲区| 一卡2卡三卡四卡精品乱码亚洲| 全区人妻精品视频| 有码 亚洲区| 中出人妻视频一区二区| 国产老妇女一区| 久久久久久久午夜电影| 亚洲av电影不卡..在线观看| 欧美一区二区精品小视频在线| 观看美女的网站| 久久久久精品国产欧美久久久| 日本-黄色视频高清免费观看| or卡值多少钱| 中文字幕免费在线视频6| 亚洲第一区二区三区不卡| 亚洲精品国产成人久久av| 国产精品永久免费网站| 一区二区三区高清视频在线| 久久精品国产清高在天天线| 成人特级av手机在线观看| 两个人的视频大全免费| 久久婷婷人人爽人人干人人爱| 精品免费久久久久久久清纯| 在线观看午夜福利视频| 51国产日韩欧美| 99热这里只有是精品50| 美女被艹到高潮喷水动态| 不卡一级毛片| 日本黄大片高清| 国产伦精品一区二区三区视频9| 精品久久久久久久末码| 国产在视频线在精品| 少妇熟女aⅴ在线视频| 亚洲自拍偷在线| 国内毛片毛片毛片毛片毛片| 欧美另类亚洲清纯唯美| 欧美日韩黄片免| 精品免费久久久久久久清纯| 日本一本二区三区精品| 91av网一区二区| 亚洲四区av| 免费看光身美女| АⅤ资源中文在线天堂| 91精品国产九色| 人妻夜夜爽99麻豆av| 亚洲av熟女| 国产av麻豆久久久久久久| 男人舔女人下体高潮全视频| 一进一出抽搐动态| 一个人看的www免费观看视频| 色综合站精品国产| 校园春色视频在线观看| 国产乱人视频| 男女视频在线观看网站免费| 看黄色毛片网站| 国产精品一及| 国产单亲对白刺激| 日本黄大片高清| 欧美在线一区亚洲| 免费人成在线观看视频色| av国产免费在线观看| 久久午夜亚洲精品久久| 国产精品人妻久久久久久| 免费av毛片视频| xxxwww97欧美| 亚洲中文字幕日韩| 免费大片18禁| 精品久久久噜噜| 久久精品国产自在天天线| 亚洲最大成人手机在线| 日韩,欧美,国产一区二区三区 | 国产v大片淫在线免费观看| 成人国产综合亚洲| 搡老岳熟女国产| 久久久久久久亚洲中文字幕| 伊人久久精品亚洲午夜| 成人三级黄色视频| 亚洲aⅴ乱码一区二区在线播放| 最近最新中文字幕大全电影3| av福利片在线观看| 国产亚洲91精品色在线| 久久午夜福利片| 九色国产91popny在线| 午夜福利高清视频| 最近视频中文字幕2019在线8| 国产v大片淫在线免费观看| 女同久久另类99精品国产91| 色吧在线观看| 成人性生交大片免费视频hd| 波多野结衣巨乳人妻| 亚洲欧美日韩无卡精品| 一区二区三区四区激情视频 | 村上凉子中文字幕在线| 丝袜美腿在线中文| 欧美成人免费av一区二区三区| 婷婷精品国产亚洲av| 床上黄色一级片| 九九爱精品视频在线观看| 男人舔奶头视频| 久久久色成人| 99久久九九国产精品国产免费| or卡值多少钱| 人妻久久中文字幕网| 欧美+日韩+精品| 老女人水多毛片| 免费av观看视频| 午夜福利高清视频| 久久久精品欧美日韩精品| x7x7x7水蜜桃| 午夜免费男女啪啪视频观看 | 内射极品少妇av片p| 日韩精品有码人妻一区| 日韩欧美精品免费久久| 日韩大尺度精品在线看网址| 久久久午夜欧美精品| 欧美日韩乱码在线| 亚洲avbb在线观看| 成人特级黄色片久久久久久久| 变态另类丝袜制服| 成人二区视频| 亚洲专区国产一区二区| 88av欧美| 精品人妻熟女av久视频| 内射极品少妇av片p| 精品午夜福利在线看| 国产v大片淫在线免费观看| 又粗又爽又猛毛片免费看| 美女cb高潮喷水在线观看| 久久久久国产精品人妻aⅴ院| 久久久久久大精品| 成人一区二区视频在线观看| 久久中文看片网| 亚洲第一电影网av| 18禁黄网站禁片免费观看直播| 亚洲成a人片在线一区二区| 麻豆久久精品国产亚洲av| 欧美xxxx黑人xx丫x性爽| 欧美精品啪啪一区二区三区| 美女xxoo啪啪120秒动态图| 天堂av国产一区二区熟女人妻| 欧美黑人欧美精品刺激| 尤物成人国产欧美一区二区三区| 日韩中字成人| 日本撒尿小便嘘嘘汇集6| 日韩一本色道免费dvd| 麻豆国产av国片精品| 国产一区二区三区视频了| 毛片女人毛片| 国产一级毛片七仙女欲春2| 国产亚洲av嫩草精品影院| 一级黄色大片毛片| 最近最新免费中文字幕在线| 欧美人与善性xxx| 琪琪午夜伦伦电影理论片6080| 精品久久久久久久久av| 亚洲精品成人久久久久久| 一边摸一边抽搐一进一小说| 天堂网av新在线| 亚洲精品亚洲一区二区| 国产精品av视频在线免费观看| 久9热在线精品视频| 欧美日韩国产亚洲二区| 国内揄拍国产精品人妻在线| 51国产日韩欧美| 日本成人三级电影网站| 免费大片18禁| 免费在线观看日本一区| www日本黄色视频网| 精品一区二区三区视频在线观看免费| 午夜福利在线观看吧| 国产在线男女| 国产成年人精品一区二区| 黄色丝袜av网址大全| 亚洲美女视频黄频| 99久久成人亚洲精品观看| 日日撸夜夜添| 欧美国产日韩亚洲一区| 久久精品国产自在天天线| 亚洲av电影不卡..在线观看| 日韩中字成人| 日韩欧美在线二视频| 精品国产三级普通话版| 国产精品久久电影中文字幕| 美女大奶头视频| 亚洲精品456在线播放app | 春色校园在线视频观看| 成人特级黄色片久久久久久久| 18禁裸乳无遮挡免费网站照片| 国产精品久久久久久久电影| 男女视频在线观看网站免费| 日韩欧美在线乱码| 国产精品久久久久久精品电影| 国内精品久久久久久久电影| 日本在线视频免费播放| 免费人成在线观看视频色| av天堂在线播放| 久久天躁狠狠躁夜夜2o2o| 欧美日本视频| 国产免费av片在线观看野外av| 国模一区二区三区四区视频| 18+在线观看网站| 日韩中字成人| 日韩 亚洲 欧美在线| 国产高清不卡午夜福利| 男人舔女人下体高潮全视频| 国产伦一二天堂av在线观看| 婷婷六月久久综合丁香| 美女cb高潮喷水在线观看| 内射极品少妇av片p| 久久国产乱子免费精品| 欧美成人一区二区免费高清观看| 两性午夜刺激爽爽歪歪视频在线观看| 黄色女人牲交| 日本五十路高清| 久久人人爽人人爽人人片va| 国产v大片淫在线免费观看| 亚洲第一区二区三区不卡| 俄罗斯特黄特色一大片| av黄色大香蕉| 国产精品无大码| 在线免费十八禁| 久久香蕉精品热| 观看美女的网站| 99精品在免费线老司机午夜| 久久婷婷人人爽人人干人人爱| 桃色一区二区三区在线观看| 精品一区二区三区人妻视频| 啪啪无遮挡十八禁网站| 中文资源天堂在线| 十八禁网站免费在线| 在线播放国产精品三级| 亚洲美女黄片视频| 久久久色成人| 亚洲色图av天堂| 国产精品精品国产色婷婷| 特大巨黑吊av在线直播| 久久久午夜欧美精品| 国产69精品久久久久777片| 嫩草影视91久久| 春色校园在线视频观看| 欧美日韩综合久久久久久 | 美女被艹到高潮喷水动态| 国产精品爽爽va在线观看网站| 欧美高清成人免费视频www| 五月伊人婷婷丁香| 又爽又黄无遮挡网站| 男人舔女人下体高潮全视频| 观看免费一级毛片| 床上黄色一级片| 国产亚洲精品av在线| 一个人免费在线观看电影| 国内精品一区二区在线观看| 天天躁日日操中文字幕| 人妻夜夜爽99麻豆av| 在线免费观看不下载黄p国产 | 国产综合懂色| 麻豆国产97在线/欧美| 欧美日本视频| 欧美zozozo另类| 99视频精品全部免费 在线| 午夜精品久久久久久毛片777| 亚洲精品在线观看二区| 欧美日韩亚洲国产一区二区在线观看| 日韩欧美三级三区| 精品久久久久久久人妻蜜臀av| 制服丝袜大香蕉在线| av天堂中文字幕网| 久久久久久久精品吃奶| 国产午夜福利久久久久久| 麻豆国产av国片精品| 国产男靠女视频免费网站| 亚洲人成网站在线播| 91麻豆精品激情在线观看国产| 啦啦啦观看免费观看视频高清| 91久久精品国产一区二区三区| 国产精品99久久久久久久久| 在线观看午夜福利视频| 久久热精品热| 99久国产av精品| 丰满人妻一区二区三区视频av| 免费一级毛片在线播放高清视频| 久久6这里有精品| 中文亚洲av片在线观看爽| 精品一区二区三区视频在线| 亚洲乱码一区二区免费版| 欧美高清性xxxxhd video| 五月玫瑰六月丁香| 熟妇人妻久久中文字幕3abv| 中亚洲国语对白在线视频| 一进一出抽搐gif免费好疼| 成人二区视频| 国产私拍福利视频在线观看| 在线国产一区二区在线| 女人十人毛片免费观看3o分钟| 国产极品精品免费视频能看的| 久久99热6这里只有精品| 国产精品av视频在线免费观看| 成人一区二区视频在线观看| 精品免费久久久久久久清纯| 不卡一级毛片| 国产精品久久久久久久久免| 深夜精品福利| 欧美绝顶高潮抽搐喷水| 97热精品久久久久久| 人妻夜夜爽99麻豆av| 99热6这里只有精品| 草草在线视频免费看| 亚洲avbb在线观看| 国产成年人精品一区二区| 欧美日本视频| 久久久久国内视频| 此物有八面人人有两片| 亚洲精品乱码久久久v下载方式| 欧美黑人欧美精品刺激| 老司机深夜福利视频在线观看| 国产亚洲av嫩草精品影院| 国产男靠女视频免费网站| 中文在线观看免费www的网站| 国产一区二区亚洲精品在线观看| 久久久久国产精品人妻aⅴ院| 午夜精品在线福利| 精华霜和精华液先用哪个| 国产精品一区二区免费欧美| 国产精品精品国产色婷婷| 亚洲人成网站在线播放欧美日韩| 很黄的视频免费| 他把我摸到了高潮在线观看| 午夜日韩欧美国产| 男女视频在线观看网站免费| 精品人妻偷拍中文字幕| 免费黄网站久久成人精品| 亚洲精品粉嫩美女一区| 亚洲av二区三区四区| 欧美xxxx性猛交bbbb| 男女下面进入的视频免费午夜| 啦啦啦韩国在线观看视频| 天美传媒精品一区二区| 亚洲中文字幕一区二区三区有码在线看| 日本黄色片子视频| 99热网站在线观看| 亚洲在线观看片| 岛国在线免费视频观看| 免费在线观看影片大全网站| 精品久久久久久成人av| 露出奶头的视频| 在线a可以看的网站| 夜夜夜夜夜久久久久| 91麻豆精品激情在线观看国产| 亚洲成av人片在线播放无| 亚洲欧美清纯卡通| 日本与韩国留学比较| 国产高清三级在线| 制服丝袜大香蕉在线| av在线蜜桃| 国产精品久久久久久久久免| 精品久久国产蜜桃| 一个人看的www免费观看视频| 非洲黑人性xxxx精品又粗又长| 国产免费av片在线观看野外av| 国产精品人妻久久久久久| 一边摸一边抽搐一进一小说| av在线观看视频网站免费| 亚洲av成人精品一区久久| 国产黄色小视频在线观看| 久久精品91蜜桃| 啦啦啦啦在线视频资源| 国产欧美日韩精品亚洲av| 久久久久国产精品人妻aⅴ院| 99精品久久久久人妻精品| 久久婷婷人人爽人人干人人爱| 人人妻人人看人人澡| 男插女下体视频免费在线播放| 精华霜和精华液先用哪个| 成人一区二区视频在线观看| 久久午夜福利片| 亚洲真实伦在线观看| 国产又黄又爽又无遮挡在线| 免费看美女性在线毛片视频| 悠悠久久av| 日韩欧美在线二视频| 久久久精品欧美日韩精品| 国产午夜福利久久久久久| 91麻豆av在线| 91麻豆av在线| 日韩欧美一区二区三区在线观看| 国产精品一区二区免费欧美| 国产精品国产高清国产av| 噜噜噜噜噜久久久久久91| 久久久国产成人精品二区| 草草在线视频免费看| 久久久久久久久久久丰满 | 九色成人免费人妻av| 午夜免费男女啪啪视频观看 | 性色avwww在线观看| 成人av一区二区三区在线看| bbb黄色大片| 91av网一区二区| 69av精品久久久久久| 久久精品国产亚洲av涩爱 | 国产白丝娇喘喷水9色精品| 国产成人影院久久av| 久久久午夜欧美精品| 特大巨黑吊av在线直播| 亚洲av免费在线观看|