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

    基于灰色預(yù)測(cè)模型的我國(guó)心腦血管疾病死亡率預(yù)測(cè)

    2015-07-28 12:45:05李論
    現(xiàn)代電子技術(shù) 2015年11期
    關(guān)鍵詞:灰色預(yù)測(cè)心腦血管病死亡率

    李論

    摘 要: 運(yùn)用灰色預(yù)測(cè)模型GM(1,1)對(duì)我國(guó)心腦血管病的死亡率進(jìn)行建模擬合和預(yù)測(cè)分析,為心腦血管病防控工作重點(diǎn)的有效確定及防控措施的合理制定提供科學(xué)依據(jù)。以全國(guó)2003—2012年10年間的心腦血管病死亡率為研究對(duì)象,采用灰色預(yù)測(cè)模型GM(1,1)對(duì)2013—2020年的全國(guó)心腦血管病死亡率進(jìn)行預(yù)測(cè),利用后驗(yàn)差檢驗(yàn),通過(guò)小誤差概率p值,方差比c值,誤差q值,綜合評(píng)價(jià)模型的預(yù)測(cè)精度。建模分析的結(jié)果表明,全國(guó)心腦血管病死亡率呈逐年上升的趨勢(shì),其中,心臟病死亡率上升速度較快,而腦血管病的死亡率上升速度相對(duì)平緩。灰色系統(tǒng)理論可以很好地運(yùn)用到死亡率預(yù)測(cè)分析方面,未來(lái)幾年我國(guó)心腦血管病死亡率仍居高不下,疾病預(yù)防與監(jiān)控部門(mén)應(yīng)做好防治工作,進(jìn)一步提升我國(guó)的醫(yī)療衛(wèi)生水準(zhǔn)。

    關(guān)鍵詞: 灰色預(yù)測(cè); 心腦血管??; 死亡率; Matlab

    中圖分類號(hào): TN802?34 文獻(xiàn)標(biāo)識(shí)碼: A 文章編號(hào): 1004?373X(2015)11?0107?05

    Mortality prediction based on grey forecasting model for cardio?cerebrovascular disease

    LI Lun

    (Business College, China University of Political Science and Law, Beijing 102249, China)

    Abstract: The modeling fitting and prediction analysis on mortality of cardio?cerebrovascular disease in China is proceeded with grey forecasting model GM(1,1), which provides scientific basis to effectively determine the emphasis of prediction and control work, and reasonably formulate the prevention and control measures for cardio?cerebrovascular disease. The mortality of cardio?cerebrovascular disease in China from the year 2003 to 2012 is selected as the research object, and the mortality of cardio?cerebrovascular disease in China from the year 2013 to 2020 is predicted by grey forecasting model GM(1,1). The prediction accuracy of the forecasting model is checked by posterior error, and evaluated by small error probability p, variance ratio c, error q comprehensively. Results of modeling analysis show that the mortality of cardio?cerebrovascular disease is increased year by year, and the mortality of heart disease is increased more obviously while the mortality of the cerebrovascular disease is increased smoothly. The grey system theory is well applied in mortality forecast and analysis. The mortality of cardio?cerebrovascular disease in China should be remain high in the next few years. Disease prevention and monitoring departments should do better in prevention and control work to further promote the medical and health level in China.

    Keywords: grey forecasting; cardio?cerebrovascular disease; mortality; Matlab

    0 引 言

    心腦血管疾病是心血管疾病和腦血管疾病的統(tǒng)稱,泛指由于高脂血癥、血液黏稠、動(dòng)脈粥樣硬化、高血壓等所導(dǎo)致的心臟、大腦及全身組織發(fā)生缺血性或出血性疾病的通稱。在過(guò)去的十幾年間,我國(guó)心腦血管疾病死亡率為全球最高,一直居高不下。由圖1可知,從2003—2012年,我國(guó)心腦血管病死亡率呈明顯上升的趨勢(shì),心臟病死亡率較腦血管病,其增長(zhǎng)速度較快。有效控制心腦血管病死亡率增長(zhǎng)趨勢(shì),成為涉及國(guó)計(jì)民生的重大問(wèn)題。為探索和建立適合中國(guó)人群的心腦血管疾病患者全程管理體系,遏制中國(guó)心腦血管疾病死亡率逐年上升的趨勢(shì),需要對(duì)我國(guó)未來(lái)心腦血管疾病死亡率進(jìn)行科學(xué)準(zhǔn)確的預(yù)測(cè)。本文擬應(yīng)用GM(1,1)模型,根據(jù)國(guó)家衛(wèi)生和計(jì)劃生育委員會(huì)及統(tǒng)計(jì)信息中心發(fā)布的數(shù)據(jù),對(duì)心血管病的死亡率進(jìn)行建模擬合,并預(yù)測(cè)出未來(lái)幾年的心腦血管死亡率。

    1 資料與方法

    1.1 資料

    本文的原始數(shù)據(jù)來(lái)源于《中國(guó)衛(wèi)生和計(jì)劃生育統(tǒng)計(jì)年鑒》,其收集的數(shù)據(jù)質(zhì)量可靠、穩(wěn)定,較為真實(shí)地反映了實(shí)際水平。

    1.2 方法

    利用灰色預(yù)測(cè)模型GM(1,1)對(duì)我國(guó)心腦血管疾病的死亡率進(jìn)行預(yù)測(cè)。GM(1,1)灰色預(yù)測(cè)法由華中科技大學(xué)控制科學(xué)與工程系教授鄧聚龍于20世紀(jì)80年代最先提出的一種新型模型計(jì)算方法,該方法摒棄了傳統(tǒng)的直接在歷史數(shù)據(jù)中尋找統(tǒng)計(jì)規(guī)律及概率分布的觀念,而是將無(wú)規(guī)律的原始數(shù)據(jù)通過(guò)一次或多次累加的數(shù)據(jù)處理方式,使其轉(zhuǎn)化為較有規(guī)律性的一組時(shí)間序列,并利用離散數(shù)列建立微分方程的動(dòng)態(tài)模型,該模型揭示的是歷史數(shù)據(jù)與環(huán)境變量之間存在的動(dòng)態(tài)關(guān)系,即以過(guò)去和現(xiàn)有歷史數(shù)據(jù)為背景,建立起來(lái)的一個(gè)從過(guò)去延伸到未來(lái)的預(yù)測(cè)模型。

    GM(1,1)灰色預(yù)測(cè)法適用于部分信息已知、部分信息未知的灰色系統(tǒng),該預(yù)測(cè)方法認(rèn)為對(duì)既含有已知信息同時(shí)又含有未知信息(或部分信息具有不確定性)的系統(tǒng)進(jìn)行預(yù)測(cè),即對(duì)在一定的方位內(nèi)變化著的、同時(shí)又與時(shí)間有關(guān)的灰色過(guò)程的預(yù)測(cè)。盡管在灰色預(yù)測(cè)過(guò)程中所顯示的現(xiàn)象往往表現(xiàn)是隨機(jī)的、雜亂無(wú)章的,但灰色系統(tǒng)整體卻表現(xiàn)為有序性、有界性。為此灰色系統(tǒng)中的數(shù)據(jù)集合存在著潛在的規(guī)律性,而GM(1,1)灰色預(yù)測(cè)法即是利用這種潛在的規(guī)律性建立起來(lái)的。

    灰色預(yù)測(cè)模型GM(1,1)的基本原理:

    GM(1,1)反映了一個(gè)變量對(duì)時(shí)間的一階微分函數(shù),其相應(yīng)的微分方程為:

    [dx(1)dt+ax(1)=u]

    式中:[x(1)]為經(jīng)過(guò)一次累加生成的數(shù)列;[t]為時(shí)間;[a,u] 為待估參數(shù),分別稱為發(fā)展灰數(shù)和內(nèi)生控制灰數(shù)。

    (1) 建立一次累加生成數(shù)列。設(shè)原始數(shù)列為:

    [x(0)={x(0)(1),x(0)(2),…,x(0)(n)},i=1,2,…,n]

    按下述方法做一次累加,得到生成數(shù)列(為樣本空間):

    [x(1)(i)=m=1ix(0)(m), i=1,2,…,n]

    (2) 利用最小二乘法求參數(shù)[a,u。]設(shè):

    [B=-12[x(1)(1)+x(1)(2)]1-12[x(1)(2)+x(1)(3)]1??-12[x(1)(n-1)+x(1)(n)]1]

    [yn=x(0)(2),x(0)(3),…,x(0)(n)T]

    參數(shù)辨識(shí)[a,u:]

    [a=au=(BTB)-1BTyn]

    (3) 求出GM(1,1)的模型:

    [x(1)(i+1)=(x(0)(1)-ua)e-ai+ua]

    [x(0)(1)=x(1)(1)x(0)(i)=x(1)(i)-x(1)(i-1),i=2,3,…,n]

    (4) 對(duì)模型精度的檢驗(yàn)。檢驗(yàn)的方法有殘差檢驗(yàn)、關(guān)聯(lián)度檢驗(yàn)和后驗(yàn)差檢驗(yàn),本文中采取的是后驗(yàn)差檢驗(yàn)。

    首先計(jì)算原始數(shù)列[x(0)(i)]的均方差[S0。]其定義為:

    [S0=S20n-1,S20=i=1nx(0)(i)-x(0)2,x(0)=1ni=1nx(0)(i)]

    然后計(jì)算殘差數(shù)列[ε(0)(i)=x(0)(i)-x(0)(i)]的均方差[S1]。其定義為:

    [S1=S21n-1,S21=i=1nε(0)(i)-ε(0)2,ε(0)=1ni=1nε(0)(i)]

    由此計(jì)算方差比:[c=S1S0]和小誤差概率:[p=ε(0)(i)-ε(0)<0.674 5?S0。]

    最后根據(jù)預(yù)測(cè)精度等級(jí)劃分表,見(jiàn)表1,檢驗(yàn)得出模型的預(yù)測(cè)精度。

    (5) 如果檢驗(yàn)合格,則可以用模型進(jìn)行預(yù)測(cè)。即用:

    [x(0)(n+1)=x(1)(n+1)-x(1)(n),x(0)(n+2)=x(1)(n+2)-x(1)(n+1),…]

    作為[x(0)(n+1),x(0)(n+2),…] 的預(yù)測(cè)值。

    灰色預(yù)測(cè)模型GM(1,1)的程序框圖如圖2所示。

    2 結(jié)果與分析

    2.1 模型預(yù)測(cè)的準(zhǔn)確性檢驗(yàn)

    利用2003—2010年的數(shù)據(jù)建立灰色預(yù)測(cè)模型,根據(jù)模型得出2011年,2012年的死亡率預(yù)測(cè)值,并將預(yù)測(cè)值與真實(shí)值進(jìn)行比較,對(duì)模型預(yù)測(cè)的準(zhǔn)確性進(jìn)行檢驗(yàn)。

    根據(jù)灰色預(yù)測(cè)的基本原理,計(jì)算出參數(shù)辨識(shí)[a,u,]進(jìn)而求出我國(guó)心腦血管病死亡率的預(yù)測(cè)模型,對(duì)2011年,2012年的死亡率進(jìn)行預(yù)測(cè)。

    以心腦血管病為例,利用Matlab求得參數(shù)[a=-0.067 9,u=310.884 8,]預(yù)測(cè)模型為[x(1)(i+1)=] [4 966.76e-ai-4 578.57。]同理,可得出心臟病死亡率的預(yù)測(cè)模型為[x(1)(i+1)=1 719.44e-ai-1 544.63,]腦血管病的死亡率預(yù)測(cè)模型為[x(1)(i+1)=3 519.63e-ai-3 306.25。]在預(yù)測(cè)模型的基礎(chǔ)上對(duì)歷年心腦血管病的死亡率進(jìn)行預(yù)測(cè),預(yù)測(cè)結(jié)果見(jiàn)表2。

    在表2中,各年度心腦血管病死亡率實(shí)際值等于各年度心臟病和腦血管病的死亡率實(shí)際值之和。

    對(duì)我國(guó)心腦血管病建立的預(yù)測(cè)模型進(jìn)行擬合優(yōu)度檢驗(yàn),結(jié)果見(jiàn)表3。

    對(duì)我國(guó)心腦血管病2011年,2012年死亡率的實(shí)際值與預(yù)測(cè)值進(jìn)行比較,同時(shí)對(duì)2011年,2012年死亡率的實(shí)際值與預(yù)測(cè)值進(jìn)行卡方檢驗(yàn),進(jìn)一步檢驗(yàn)?zāi)P偷臄M合優(yōu)度,檢驗(yàn)結(jié)果見(jiàn)表4。

    從表4中可以看出,根據(jù)2003—2010年的我國(guó)心腦血管病死亡率預(yù)測(cè)的2011年,2012年死亡率和實(shí)際值并沒(méi)有顯著性差異,實(shí)際值與預(yù)測(cè)值通過(guò)了卡方檢驗(yàn)。另一方面,通過(guò)對(duì)總的心腦血管病的死亡率預(yù)測(cè),將其與預(yù)測(cè)所得的心臟病、腦血管病死亡率之和進(jìn)行比較,在一定程度上也可以評(píng)價(jià)模型的擬合效果。經(jīng)查可知,預(yù)測(cè)所得2011年心腦血管病死亡率為561.17,預(yù)測(cè)所得心臟病和腦血管病的死亡率之和為561.73,二者幾乎相等;預(yù)測(cè)所得2012年心腦血管病死亡率為600.58,預(yù)測(cè)所得心臟病和腦血管病的死亡率之和為601.69,二者無(wú)顯著性差異。故可看出,模型的擬合效果較好。

    2.2 預(yù)測(cè)模型的建立及評(píng)價(jià)

    利用2003—2010年的死亡率數(shù)據(jù)建立模型,以心腦血管病為例,求得參數(shù):[a=-0.049 7,][u=337.697 4,]預(yù)測(cè)模型為[x(1)(i+1)=7 182.91e-ai-][6 794.72。]同理,可建立心臟病死亡率的預(yù)測(cè)模型為[x(1)(i+1)=2 237.54e-ai-][2 062.73,]腦血管病的死亡率預(yù)測(cè)模型為[x(1)(i+1)=][5 970.85e-ai-5 757.47。]在預(yù)測(cè)模型的基礎(chǔ)上對(duì)歷年心腦血管病死亡率進(jìn)行預(yù)測(cè),預(yù)測(cè)結(jié)果見(jiàn)表5。

    對(duì)我國(guó)心腦血管病建立的預(yù)測(cè)模型進(jìn)行擬合優(yōu)度檢驗(yàn),結(jié)果見(jiàn)表6。

    從表6可以得知,心腦血管病、心臟病以及腦血管病的擬合優(yōu)度都較好,預(yù)測(cè)精度等級(jí)高。

    2.3 外推預(yù)測(cè)

    對(duì)未來(lái)幾年我國(guó)心腦血管病的死亡率進(jìn)行外推預(yù)測(cè),結(jié)果見(jiàn)表7。

    2.4 對(duì)結(jié)果進(jìn)行檢驗(yàn)和評(píng)價(jià)

    為進(jìn)一步對(duì)我國(guó)心腦血管病死亡率變化分析,利用表6、表7的數(shù)據(jù),繪制2003—2020年我國(guó)心腦血管病死亡率變化趨勢(shì)圖,如圖3所示。由圖3可知,我國(guó)心腦血管病死亡率在未來(lái)幾年仍呈現(xiàn)出增長(zhǎng)趨勢(shì),預(yù)計(jì)2015年心腦血管病、心臟病、腦血管病的死亡率分別為632.54,328.47,307.70;2020年心腦血管病、心臟病、腦血管病的死亡率分別為811.11,460.42,365.88。短短5年時(shí)間,心臟病死亡率的增長(zhǎng)幅度達(dá)到40.17%,而腦血管病的死亡率只增長(zhǎng)了18.91%,由此可以進(jìn)一步得出結(jié)論:心臟病將成為我國(guó)未來(lái)人口死亡的重要誘因,相關(guān)部門(mén)應(yīng)加強(qiáng)心臟病的預(yù)防與控制工作。

    3 結(jié) 論

    GM(1,1)模型是單變量一階線性模型,通過(guò)時(shí)間序列的研究去尋找和探索事物發(fā)展變化的連續(xù)或離散的未來(lái)時(shí)間序列,分析事物發(fā)展變化的內(nèi)部規(guī)律。該理論是灰色預(yù)測(cè)模型中最基本的模型,預(yù)測(cè)效果好,有較高的適用性,應(yīng)用價(jià)值很大。與傳統(tǒng)的數(shù)學(xué)模型相比,在疾病預(yù)測(cè)方面具有所需樣本量小、無(wú)需典型概率分布、計(jì)算簡(jiǎn)單、預(yù)測(cè)效果好等眾多優(yōu)點(diǎn),可以作為疾病監(jiān)測(cè)的有效工具。

    隨著當(dāng)今經(jīng)濟(jì)的飛速發(fā)展,人們的生活方式、飲食結(jié)構(gòu)發(fā)生了很大變化;同時(shí),伴隨著老齡化程度的加速,人均壽命的延長(zhǎng),我國(guó)心腦血管病的死亡率呈現(xiàn)出增長(zhǎng)趨勢(shì)。本文利用《中國(guó)衛(wèi)生和計(jì)劃生育統(tǒng)計(jì)年鑒》的資料,用2003—2010及2003—2012年數(shù)據(jù)分別建立GM(1,1)模型,從中可以看出,全國(guó)心腦血管疾病的死亡率逐年上升。其中,心臟病死亡率增長(zhǎng)速度較快,而腦血管病的死亡率增長(zhǎng)速度相對(duì)平緩。因此,有關(guān)部門(mén)應(yīng)重視心腦血管病的防治工作,加強(qiáng)心腦血管病的防治研究,投入必要的經(jīng)費(fèi)和人員,心腦血管病的防治工作應(yīng)成為現(xiàn)階段疾病控制與預(yù)防的長(zhǎng)期任務(wù)之一;同時(shí),研究適合我國(guó)國(guó)情的心腦血管病社區(qū)綜合防治策略,對(duì)不同地區(qū)、不同人群提出切實(shí)有效的防治方案,以期在不遠(yuǎn)的將來(lái)遏制心腦血管病死亡率上升的趨勢(shì)。

    參考文獻(xiàn)

    [1] HAMZACEBI C, HUSEYIN A E. Forecasting the annual electricity consumption of Turkey using an optimized grey model [J]. Energy, 2014, 70(2): 165?171.

    [2] PAUL H K H. Forecasting construction manpower demand by gray model [J]. Journal of Construction Engineering & Management, 2010, 136(12): 1299?1305.

    [3] WU W Y, HSIAO S W, TSAI C H. Forecasting and evaluating the tourist hotel industry performance in Taiwan based on Grey theory [J]. Tourism & Hospitality Research, 2008, 8(2): 137?152.

    [4] WEN Y H. Shipment forecasting for supply chain collaborative transportation management using grey models with grey numbers [J]. Transportation Planning & Technology, 2011, 34(6): 605?624.

    [5] 陳金宏,何耀,秦銀河.心血管病預(yù)測(cè)模型研究進(jìn)展[J].中華流行病學(xué)雜志,2010(11):1311?1314.

    [6] 方建衛(wèi),王文娟,楚霹.基于最小二乘法的GM(1,1)模型在人口預(yù)測(cè)中的應(yīng)用[J].貴州大學(xué)學(xué)報(bào):自然科學(xué)版,2007,24(4):345?349.

    [7] 高月霞,倪卓健,李小杉,等.1989—2009年海門(mén)市糖尿病死亡分析及預(yù)測(cè)[J].中國(guó)衛(wèi)生統(tǒng)計(jì),2013(1):74?76.

    [8] 黃健元,劉洋.流動(dòng)人口預(yù)測(cè)模型構(gòu)建及其應(yīng)用[J].統(tǒng)計(jì)與決策,2008(23):18?20.

    [9] 李濱,張衛(wèi)民.GM(1,1)殘差修正模型在乙肝發(fā)病預(yù)測(cè)中的應(yīng)用[J].中國(guó)衛(wèi)生統(tǒng)計(jì),2014(2):240?242.

    [10] 李恩臨.灰色預(yù)測(cè)方法的應(yīng)用[J].統(tǒng)計(jì)與決策,2010(14):161?162.

    [11] 劉磊,褚金花,尹愛(ài)田.灰色模型方法在區(qū)域衛(wèi)生人力預(yù)測(cè)中的應(yīng)用[J].中國(guó)衛(wèi)生統(tǒng)計(jì),2011(4):453?454.

    [12] 邵珠艷,王春梅,魏曼莎.灰色GM(1,1)預(yù)測(cè)模型在疾病預(yù)測(cè)中的應(yīng)用[J].中國(guó)醫(yī)院統(tǒng)計(jì),2003,10(3):146?148.

    [13] 汪叢峰,左登敏,洪光烈.灰色系統(tǒng)理論在傳染病優(yōu)勢(shì)病種確定及疫情預(yù)測(cè)的應(yīng)用[J].安徽預(yù)防醫(yī)學(xué)雜志,2011,17(3):180?183.

    [14] 王曉軍,蔡正高.死亡率預(yù)測(cè)模型的新進(jìn)展[J].統(tǒng)計(jì)研究,2008,25(9):80?84.

    [15] 王瑞娜,唐德善.基于改進(jìn)的灰色GM(1,1)模型的人口預(yù)測(cè)[J].統(tǒng)計(jì)與決策,2007(20):93?95.

    [16] 吳菊珍,徐曄,龔新橋.灰色預(yù)測(cè)在股票價(jià)格預(yù)測(cè)中的應(yīng)用[J].統(tǒng)計(jì)與決策,2007(1):139?140.

    [17] 袁鶯楹,董建成.基于數(shù)學(xué)模型的疾病預(yù)測(cè)方法比較研究[J].軟件導(dǎo)刊,2009(5):108?110.

    [18] 周脈耕,楊功煥.灰色數(shù)列GM(1,1)模型在心腦血管疾病死亡預(yù)測(cè)中的應(yīng)用[J].中國(guó)公共衛(wèi)生,2002,18(5):625?627.

    猜你喜歡
    灰色預(yù)測(cè)心腦血管病死亡率
    走路可以降低死亡率
    中老年保健(2021年4期)2021-08-22 07:07:02
    春季養(yǎng)雞這樣降低死亡率
    絡(luò)病理論在心腦血管病治未病中的應(yīng)用初探
    新冠肺炎的死亡率為何難確定?
    急性爛鰓、套腸、敗血癥…一旦治療不及時(shí),死亡率或高達(dá)90%,叉尾鮰真的值得養(yǎng)嗎?
    沉默的殺手——心腦血管病,5秒要人命,該如何預(yù)防?
    《中西醫(yī)結(jié)合心腦血管病雜志》第二屆編委會(huì)名單
    2016年世界園藝博覽會(huì)對(duì)唐山經(jīng)濟(jì)的影響
    卷宗(2016年10期)2017-01-21 02:23:50
    收益還原法在房地產(chǎn)估價(jià)工作中的應(yīng)用與改進(jìn)
    全面放開(kāi)二胎政策對(duì)蚌埠市人口結(jié)構(gòu)的影響的探索
    亚洲综合色网址| 国产成人精品久久二区二区91| 他把我摸到了高潮在线观看| 日本黄色视频三级网站网址 | 精品人妻1区二区| 欧美不卡视频在线免费观看 | 久久精品国产99精品国产亚洲性色 | 热99国产精品久久久久久7| svipshipincom国产片| 大码成人一级视频| 热99re8久久精品国产| 日本vs欧美在线观看视频| 丝袜人妻中文字幕| 精品福利永久在线观看| 性色av乱码一区二区三区2| 90打野战视频偷拍视频| 久久中文看片网| 国产不卡av网站在线观看| 国产免费av片在线观看野外av| 搡老岳熟女国产| 成人国产一区最新在线观看| 久久久久久免费高清国产稀缺| 91精品国产国语对白视频| 制服人妻中文乱码| 午夜久久久在线观看| 国产精品久久电影中文字幕 | 自拍欧美九色日韩亚洲蝌蚪91| 午夜老司机福利片| 亚洲人成伊人成综合网2020| 亚洲人成电影免费在线| 免费看十八禁软件| 最近最新免费中文字幕在线| 18在线观看网站| 久久精品成人免费网站| 麻豆av在线久日| 中国美女看黄片| 成人手机av| 如日韩欧美国产精品一区二区三区| 欧美亚洲 丝袜 人妻 在线| 国产黄色免费在线视频| 一边摸一边抽搐一进一小说 | a级毛片黄视频| www日本在线高清视频| 精品高清国产在线一区| 精品久久久久久久久久免费视频 | 18在线观看网站| 夜夜爽天天搞| 日韩大码丰满熟妇| 欧美亚洲 丝袜 人妻 在线| 交换朋友夫妻互换小说| 交换朋友夫妻互换小说| 国产精品久久视频播放| 精品福利永久在线观看| 嫁个100分男人电影在线观看| 男女下面插进去视频免费观看| 精品免费久久久久久久清纯 | 黄片播放在线免费| 国产成人av激情在线播放| www.熟女人妻精品国产| 亚洲精品国产色婷婷电影| 久久草成人影院| 自拍欧美九色日韩亚洲蝌蚪91| 欧美中文综合在线视频| 欧美丝袜亚洲另类 | 亚洲av成人一区二区三| 天天躁日日躁夜夜躁夜夜| 欧美乱码精品一区二区三区| 国产1区2区3区精品| 9191精品国产免费久久| 又黄又粗又硬又大视频| 丝袜美足系列| 咕卡用的链子| 又黄又粗又硬又大视频| 91九色精品人成在线观看| 美女福利国产在线| 欧美精品亚洲一区二区| 又紧又爽又黄一区二区| 少妇 在线观看| 999久久久精品免费观看国产| ponron亚洲| av线在线观看网站| 韩国精品一区二区三区| 欧美av亚洲av综合av国产av| tube8黄色片| 欧美老熟妇乱子伦牲交| 天堂动漫精品| 国产亚洲欧美在线一区二区| 99在线人妻在线中文字幕 | 欧美日韩亚洲国产一区二区在线观看 | 国产精品成人在线| 久久精品91无色码中文字幕| 欧美精品人与动牲交sv欧美| 中亚洲国语对白在线视频| 欧美成狂野欧美在线观看| 视频区图区小说| 欧美黄色淫秽网站| 丰满人妻熟妇乱又伦精品不卡| 9色porny在线观看| 夜夜躁狠狠躁天天躁| 真人做人爱边吃奶动态| 亚洲国产精品合色在线| 女人被躁到高潮嗷嗷叫费观| 欧美精品高潮呻吟av久久| 他把我摸到了高潮在线观看| 99热只有精品国产| 亚洲九九香蕉| 久久久久久久精品吃奶| 桃红色精品国产亚洲av| 国产精品秋霞免费鲁丝片| 免费不卡黄色视频| 亚洲欧洲精品一区二区精品久久久| 精品国产乱子伦一区二区三区| 一级毛片精品| 免费日韩欧美在线观看| 亚洲成人免费电影在线观看| a级毛片黄视频| 国产乱人伦免费视频| 美女视频免费永久观看网站| 免费在线观看亚洲国产| 69精品国产乱码久久久| 久久这里只有精品19| av片东京热男人的天堂| 免费av中文字幕在线| 久久久久久久久免费视频了| 中文字幕精品免费在线观看视频| 国产精品亚洲一级av第二区| 激情视频va一区二区三区| 自线自在国产av| 国产精品香港三级国产av潘金莲| 欧美老熟妇乱子伦牲交| 五月开心婷婷网| 久久天躁狠狠躁夜夜2o2o| aaaaa片日本免费| 亚洲国产中文字幕在线视频| 精品国产乱子伦一区二区三区| 老司机深夜福利视频在线观看| 黑人巨大精品欧美一区二区mp4| 国产精品亚洲一级av第二区| 99国产极品粉嫩在线观看| 国产又爽黄色视频| 成年人免费黄色播放视频| 国产伦人伦偷精品视频| 人成视频在线观看免费观看| 1024视频免费在线观看| 超色免费av| 久久草成人影院| 精品久久久久久久久久免费视频 | 亚洲成av片中文字幕在线观看| 亚洲精品国产精品久久久不卡| 久久中文字幕人妻熟女| 超色免费av| 色综合欧美亚洲国产小说| 精品国产乱子伦一区二区三区| 一二三四在线观看免费中文在| 女人高潮潮喷娇喘18禁视频| 国产91精品成人一区二区三区| 午夜福利免费观看在线| 久久精品成人免费网站| 亚洲成a人片在线一区二区| 欧美大码av| 国产在视频线精品| 身体一侧抽搐| 一区福利在线观看| 亚洲九九香蕉| 国产av又大| 久久久国产一区二区| 久久青草综合色| 欧美精品亚洲一区二区| 国产野战对白在线观看| 国产欧美日韩一区二区三| 叶爱在线成人免费视频播放| 男人的好看免费观看在线视频 | 国产野战对白在线观看| 欧美在线黄色| 美女高潮到喷水免费观看| 久久精品91无色码中文字幕| 老司机午夜十八禁免费视频| 久久影院123| 欧美 亚洲 国产 日韩一| 女同久久另类99精品国产91| 久久精品国产a三级三级三级| 在线观看免费视频网站a站| 国产成人精品在线电影| 免费看十八禁软件| 日韩一卡2卡3卡4卡2021年| 最近最新中文字幕大全电影3 | 国产精品 欧美亚洲| 黄色毛片三级朝国网站| 老熟妇乱子伦视频在线观看| 亚洲精品在线观看二区| 99精品欧美一区二区三区四区| 一级,二级,三级黄色视频| 99久久综合精品五月天人人| 91老司机精品| 在线国产一区二区在线| 国产视频一区二区在线看| 天天躁狠狠躁夜夜躁狠狠躁| 女人久久www免费人成看片| 免费看十八禁软件| 狂野欧美激情性xxxx| 看黄色毛片网站| 国产一区二区三区视频了| 精品视频人人做人人爽| 99国产综合亚洲精品| 久久久精品国产亚洲av高清涩受| 欧美另类亚洲清纯唯美| 美国免费a级毛片| 亚洲国产欧美日韩在线播放| 欧美黑人精品巨大| 亚洲欧美一区二区三区黑人| 欧美成狂野欧美在线观看| 久久久久精品国产欧美久久久| 女人爽到高潮嗷嗷叫在线视频| 丰满饥渴人妻一区二区三| 91成人精品电影| 国产午夜精品久久久久久| 中亚洲国语对白在线视频| 成年女人毛片免费观看观看9 | 三上悠亚av全集在线观看| 久久精品国产清高在天天线| 亚洲av日韩精品久久久久久密| 男人操女人黄网站| 日韩欧美在线二视频 | 日日爽夜夜爽网站| 一级毛片高清免费大全| 久久婷婷成人综合色麻豆| 在线免费观看的www视频| 69av精品久久久久久| 18禁裸乳无遮挡动漫免费视频| 看片在线看免费视频| ponron亚洲| 久久久久精品人妻al黑| 久久久久国内视频| 久久久国产一区二区| 国产亚洲精品一区二区www | tocl精华| 最近最新中文字幕大全电影3 | 国产精品 欧美亚洲| 午夜91福利影院| 三上悠亚av全集在线观看| 亚洲男人天堂网一区| 国产精品久久久久久人妻精品电影| 99久久精品国产亚洲精品| 国精品久久久久久国模美| 黄网站色视频无遮挡免费观看| 欧美黑人精品巨大| 在线看a的网站| 日本vs欧美在线观看视频| 真人做人爱边吃奶动态| 一进一出抽搐gif免费好疼 | 19禁男女啪啪无遮挡网站| 欧美成狂野欧美在线观看| 国产片内射在线| 757午夜福利合集在线观看| 久久久国产欧美日韩av| 国产有黄有色有爽视频| 精品久久久久久电影网| 啦啦啦视频在线资源免费观看| 热re99久久精品国产66热6| 人妻一区二区av| 超色免费av| 久久热在线av| 亚洲欧美激情在线| 久久精品国产清高在天天线| 在线观看日韩欧美| 美女福利国产在线| 国产精品久久久av美女十八| 精品国产乱码久久久久久男人| 狂野欧美激情性xxxx| 国产乱人伦免费视频| 亚洲视频免费观看视频| 这个男人来自地球电影免费观看| 亚洲免费av在线视频| 又黄又粗又硬又大视频| 午夜福利欧美成人| 亚洲av日韩在线播放| 自拍欧美九色日韩亚洲蝌蚪91| 欧美国产精品va在线观看不卡| 亚洲国产毛片av蜜桃av| 少妇裸体淫交视频免费看高清 | 丰满的人妻完整版| 免费人成视频x8x8入口观看| 多毛熟女@视频| 国产精品亚洲一级av第二区| 人人妻人人爽人人添夜夜欢视频| 久久午夜亚洲精品久久| 人人妻人人添人人爽欧美一区卜| 久久久精品区二区三区| 深夜精品福利| 麻豆国产av国片精品| 亚洲情色 制服丝袜| 老熟妇仑乱视频hdxx| 久久精品人人爽人人爽视色| 欧美丝袜亚洲另类 | 在线国产一区二区在线| 午夜精品国产一区二区电影| 久久久国产精品麻豆| 两人在一起打扑克的视频| 老司机亚洲免费影院| 日韩欧美一区二区三区在线观看 | √禁漫天堂资源中文www| 麻豆成人av在线观看| 亚洲精品国产精品久久久不卡| 操美女的视频在线观看| 又黄又爽又免费观看的视频| 人人妻人人爽人人添夜夜欢视频| 欧美成狂野欧美在线观看| 少妇猛男粗大的猛烈进出视频| 国产男靠女视频免费网站| а√天堂www在线а√下载 | 老熟妇仑乱视频hdxx| 日韩三级视频一区二区三区| 午夜精品久久久久久毛片777| 免费在线观看日本一区| 亚洲专区中文字幕在线| 亚洲av电影在线进入| 亚洲成国产人片在线观看| 美女国产高潮福利片在线看| av视频免费观看在线观看| 欧美日韩福利视频一区二区| 国产免费男女视频| 黄色女人牲交| 国产熟女午夜一区二区三区| 日日摸夜夜添夜夜添小说| 色在线成人网| 欧美亚洲 丝袜 人妻 在线| 亚洲在线自拍视频| 91字幕亚洲| 在线看a的网站| 国产人伦9x9x在线观看| 欧美日韩黄片免| 国产精品一区二区在线观看99| 精品欧美一区二区三区在线| 亚洲精品成人av观看孕妇| 国产精品二区激情视频| 国产麻豆69| 一级a爱视频在线免费观看| 欧美黄色片欧美黄色片| 国产精品美女特级片免费视频播放器 | 亚洲av熟女| 国产精品 欧美亚洲| 激情视频va一区二区三区| 天天添夜夜摸| 精品午夜福利视频在线观看一区| 一二三四在线观看免费中文在| 亚洲久久久国产精品| 国产黄色免费在线视频| 久久精品91无色码中文字幕| 国产在线一区二区三区精| 国产亚洲欧美98| 天堂中文最新版在线下载| 国产野战对白在线观看| 日韩中文字幕欧美一区二区| 欧美在线黄色| 大香蕉久久成人网| 女人爽到高潮嗷嗷叫在线视频| 欧美精品高潮呻吟av久久| 色精品久久人妻99蜜桃| 黄色a级毛片大全视频| 午夜免费成人在线视频| 久久久久久久久久久久大奶| 两个人看的免费小视频| 日韩欧美国产一区二区入口| 99国产精品一区二区蜜桃av | 久久精品aⅴ一区二区三区四区| 91av网站免费观看| 两个人看的免费小视频| 男人舔女人的私密视频| 天天添夜夜摸| 精品午夜福利视频在线观看一区| 超碰97精品在线观看| 极品教师在线免费播放| 国产高清国产精品国产三级| 不卡一级毛片| а√天堂www在线а√下载 | 一区福利在线观看| 欧美日韩视频精品一区| 香蕉久久夜色| 欧美黑人精品巨大| 在线观看免费午夜福利视频| 久久精品国产综合久久久| av中文乱码字幕在线| 亚洲精华国产精华精| 久久久久久久久免费视频了| 亚洲美女黄片视频| 午夜福利在线观看吧| 丝袜美足系列| 啦啦啦免费观看视频1| 免费在线观看黄色视频的| 在线观看舔阴道视频| 一a级毛片在线观看| 亚洲一卡2卡3卡4卡5卡精品中文| 亚洲国产精品合色在线| 亚洲男人天堂网一区| 国产精品亚洲av一区麻豆| 下体分泌物呈黄色| 19禁男女啪啪无遮挡网站| 大码成人一级视频| 婷婷成人精品国产| 狠狠狠狠99中文字幕| 精品国产乱子伦一区二区三区| videos熟女内射| 亚洲国产欧美一区二区综合| 黄色毛片三级朝国网站| 欧美在线黄色| 免费在线观看黄色视频的| 激情在线观看视频在线高清 | 夜夜躁狠狠躁天天躁| 欧美精品高潮呻吟av久久| 中文欧美无线码| 免费一级毛片在线播放高清视频 | 最新美女视频免费是黄的| 老汉色∧v一级毛片| 一夜夜www| 啦啦啦在线免费观看视频4| 午夜福利乱码中文字幕| 18禁黄网站禁片午夜丰满| 亚洲精品成人av观看孕妇| 国产一区二区三区在线臀色熟女 | 我的亚洲天堂| 村上凉子中文字幕在线| 中文欧美无线码| 亚洲午夜精品一区,二区,三区| 国产三级黄色录像| 欧美精品av麻豆av| 天天影视国产精品| av福利片在线| 激情在线观看视频在线高清 | 久久人人爽av亚洲精品天堂| 丁香欧美五月| 视频区图区小说| 在线观看66精品国产| 国内毛片毛片毛片毛片毛片| 人妻丰满熟妇av一区二区三区 | 国产精品av久久久久免费| 精品久久久久久久久久免费视频 | 亚洲五月色婷婷综合| 国产一区二区激情短视频| 精品欧美一区二区三区在线| 久久久久久人人人人人| 色综合欧美亚洲国产小说| 韩国av一区二区三区四区| 亚洲第一青青草原| 一级片'在线观看视频| 日韩制服丝袜自拍偷拍| 久久九九热精品免费| 欧美日韩瑟瑟在线播放| 欧美日韩中文字幕国产精品一区二区三区 | 丰满迷人的少妇在线观看| 国产不卡一卡二| 精品无人区乱码1区二区| 999久久久精品免费观看国产| 人成视频在线观看免费观看| 亚洲中文日韩欧美视频| 一进一出好大好爽视频| 在线av久久热| 18禁观看日本| 嫁个100分男人电影在线观看| 色老头精品视频在线观看| 日韩视频一区二区在线观看| 一级片'在线观看视频| 韩国精品一区二区三区| 久久九九热精品免费| 美女国产高潮福利片在线看| 老司机影院毛片| 看免费av毛片| 国产一卡二卡三卡精品| 777米奇影视久久| 亚洲少妇的诱惑av| 久久久国产成人免费| 亚洲欧美激情在线| 交换朋友夫妻互换小说| 香蕉久久夜色| 久久香蕉国产精品| 黑人欧美特级aaaaaa片| 国产蜜桃级精品一区二区三区 | 国产日韩欧美亚洲二区| 久久午夜亚洲精品久久| 午夜激情av网站| 午夜亚洲福利在线播放| 亚洲国产欧美网| 亚洲精品国产精品久久久不卡| 成人18禁高潮啪啪吃奶动态图| 啦啦啦免费观看视频1| x7x7x7水蜜桃| 亚洲国产欧美日韩在线播放| 十八禁人妻一区二区| 色婷婷久久久亚洲欧美| 国产av精品麻豆| 女人爽到高潮嗷嗷叫在线视频| 日本黄色日本黄色录像| 人人澡人人妻人| 久久精品国产清高在天天线| 午夜成年电影在线免费观看| 国产真人三级小视频在线观看| 国产亚洲欧美98| 乱人伦中国视频| 18禁黄网站禁片午夜丰满| 亚洲精品在线美女| 国产精品 国内视频| 成年人午夜在线观看视频| 亚洲av第一区精品v没综合| 热99国产精品久久久久久7| 国产有黄有色有爽视频| 亚洲av成人一区二区三| 色94色欧美一区二区| 99精品在免费线老司机午夜| 老汉色av国产亚洲站长工具| 我的亚洲天堂| 高清视频免费观看一区二区| 精品一区二区三区av网在线观看| 亚洲精品粉嫩美女一区| 啦啦啦视频在线资源免费观看| 亚洲精品国产精品久久久不卡| 美女福利国产在线| 露出奶头的视频| 国产在视频线精品| 亚洲av日韩在线播放| 黄色毛片三级朝国网站| 亚洲精品美女久久久久99蜜臀| 两个人看的免费小视频| 亚洲,欧美精品.| 欧美大码av| 欧美激情 高清一区二区三区| 成人18禁高潮啪啪吃奶动态图| 亚洲久久久国产精品| 精品第一国产精品| 宅男免费午夜| 黑人猛操日本美女一级片| 巨乳人妻的诱惑在线观看| 国产单亲对白刺激| 国产又爽黄色视频| 国产精品亚洲av一区麻豆| 午夜影院日韩av| av片东京热男人的天堂| 午夜影院日韩av| 亚洲成国产人片在线观看| 日韩成人在线观看一区二区三区| av片东京热男人的天堂| videos熟女内射| 性少妇av在线| 女人高潮潮喷娇喘18禁视频| 女人精品久久久久毛片| 精品亚洲成国产av| 国产精品亚洲av一区麻豆| 亚洲国产精品sss在线观看 | 51午夜福利影视在线观看| 搡老熟女国产l中国老女人| 欧美日韩亚洲综合一区二区三区_| 天天躁日日躁夜夜躁夜夜| 国产亚洲精品久久久久5区| 欧美人与性动交α欧美软件| 51午夜福利影视在线观看| 欧美不卡视频在线免费观看 | 免费黄频网站在线观看国产| 首页视频小说图片口味搜索| 久久久国产欧美日韩av| 欧美另类亚洲清纯唯美| 18禁国产床啪视频网站| 脱女人内裤的视频| 黄色视频,在线免费观看| 超碰成人久久| 亚洲黑人精品在线| 亚洲五月色婷婷综合| 国产黄色免费在线视频| 日本黄色视频三级网站网址 | 建设人人有责人人尽责人人享有的| 又大又爽又粗| 亚洲精品国产区一区二| 午夜精品在线福利| 国产欧美日韩一区二区三区在线| 午夜老司机福利片| 精品亚洲成a人片在线观看| 国产成人欧美在线观看 | 校园春色视频在线观看| 精品人妻在线不人妻| 中文字幕最新亚洲高清| 一进一出好大好爽视频| 国产在线精品亚洲第一网站| 亚洲成av片中文字幕在线观看| 丁香六月欧美| 香蕉丝袜av| 欧美日韩av久久| 超色免费av| 欧美日韩亚洲高清精品| 国产成人av教育| 久久国产乱子伦精品免费另类| 777米奇影视久久| 另类亚洲欧美激情| 午夜免费鲁丝| 色94色欧美一区二区| 国产精品久久久久成人av| 交换朋友夫妻互换小说| 久热爱精品视频在线9| 国产av又大| 免费看a级黄色片| √禁漫天堂资源中文www| 国产一卡二卡三卡精品| 黑人欧美特级aaaaaa片| 久久久国产精品麻豆| 国产免费男女视频| 中国美女看黄片| 欧美日韩黄片免| 国产精品秋霞免费鲁丝片| 美女午夜性视频免费| av天堂久久9| 国产一区二区激情短视频| 日韩免费高清中文字幕av| av超薄肉色丝袜交足视频| 欧美亚洲日本最大视频资源| 亚洲人成伊人成综合网2020| 国产不卡一卡二| 人人澡人人妻人| 久久精品人人爽人人爽视色| 精品国产国语对白av|