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

    基于Simulink的傳染病模型仿真

    2018-02-23 15:21:29石云鳳,楊洪明,胡宇,張洪銀,梁攀
    科技創(chuàng)新與應(yīng)用 2018年6期

    石云鳳,楊洪明,胡宇,張洪銀,梁攀

    摘 要:隨著計算機(jī)的發(fā)展,數(shù)字仿真技術(shù)發(fā)展迅速且應(yīng)用領(lǐng)域十分廣泛。而傳染病模型作為一類發(fā)展已十分成熟的數(shù)學(xué)模型,對其進(jìn)行研究的相關(guān)文獻(xiàn)數(shù)不勝數(shù),傳染病模型的實際應(yīng)用也十分廣泛,但很少有人將數(shù)字仿真技術(shù)應(yīng)用于傳染病模型。文章利用Simulink仿真軟件對傳統(tǒng)的SIR模型和在SIR模型中加入政府宣傳防控效用的傳染病模型進(jìn)行仿真。首先建立傳統(tǒng)SIR模型與在SIR模型中加入政府宣傳防控工作效用的數(shù)學(xué)模型,然后針對這兩類數(shù)學(xué)模型建立仿真系統(tǒng),最后通過設(shè)置合適的參數(shù)值,得出仿真結(jié)果并與實際情況相對比,結(jié)果表明該仿真的結(jié)果與實際情況相吻合。

    關(guān)鍵詞:數(shù)字仿真;Simulink;SIR模型;宣傳防控;過程仿真

    中圖分類號:TP391.9 文獻(xiàn)標(biāo)志碼:A 文章編號:2095-2945(2018)06-0039-03

    Abstract: With the development of computer, digital simulation technology develops rapidly and is widely used in many fields. As a very mature mathematical model, infectious disease model has been studied in numerous literatures, and the practical application of infectious disease model is also very extensive. But few people apply digital simulation technology to infectious disease model. In this paper, the Simulink simulation software is used to simulate the traditional SIR model and the infectious disease model which is added to the SIR model to promote the effectiveness of prevention and control. Firstly, the traditional SIR model and the mathematical model of the effectiveness of government propaganda and control are established in the SIR model, and then the simulation system is established for these two kinds of mathematical models. Finally, the appropriate parameters are set up. The simulation results are compared with the actual situation, and the results show that the simulation results are consistent with the actual situation.

    Keywords: digital simulation; Simulink; SIR model; publicity and control; process simulation

    引言

    現(xiàn)代數(shù)字仿真技術(shù)是以相似原理、信息技術(shù)、系統(tǒng)技術(shù)及其應(yīng)用領(lǐng)域有關(guān)的專業(yè)技術(shù)為基礎(chǔ),以計算機(jī)和各種物理效應(yīng)設(shè)備為工具,利用系統(tǒng)模型對實際的或設(shè)想的系統(tǒng)進(jìn)行試驗研究的一種綜合性技術(shù)[1]。與物理仿真技術(shù)相比,數(shù)字仿真不受被研究系統(tǒng)規(guī)模和結(jié)構(gòu)復(fù)雜性的限制,以計算速度快、使用靈活、成本相對低廉等優(yōu)勢,成為分析、研究電力系統(tǒng)必不可少的工具,是系統(tǒng)仿真的重要研究方向[2]。隨著計算機(jī)的發(fā)展,數(shù)字仿真技術(shù)取得突飛猛進(jìn)的進(jìn)步[3]。雖然數(shù)字仿真技術(shù)當(dāng)前的應(yīng)用領(lǐng)域很廣,但是在該技術(shù)領(lǐng)域依然存在一些尚未攻克的技術(shù)難點,作為下一世紀(jì)重要技術(shù)競爭點,數(shù)字仿真技術(shù)需要人們給予更多的重視。

    MATLAB中的一種可視化仿真工具Simulink,是一種基于MATLAB的框圖設(shè)計環(huán)境,實現(xiàn)動態(tài)系統(tǒng)建模、仿真和分析的一個軟件包,被廣泛應(yīng)用于線性系統(tǒng)、非線性系統(tǒng)、數(shù)字控制及數(shù)字信號處理的建模和仿真中。Simulink軟件在科學(xué)研究中應(yīng)用十分廣泛,它提供一個動態(tài)系統(tǒng)建模、仿真和綜合分析的集成環(huán)境。在該環(huán)境中,無需大量書寫程序,而只需要通過簡單直觀的鼠標(biāo)操作,就可構(gòu)造出復(fù)雜的系統(tǒng)。曾有人利用MATLAB中的Simulink仿真軟件對無刷直流電機(jī)控制系統(tǒng)建模進(jìn)行仿真[4],實現(xiàn)換熱器穩(wěn)態(tài)傳熱計算[5]。而傳染病模型作為一類發(fā)展十分成熟的數(shù)學(xué)模型,相關(guān)的研究文獻(xiàn)數(shù)不勝數(shù)。文獻(xiàn)[6]劉玉英,肖燕妮等人研究了一類受媒體影響的傳染病模型,文獻(xiàn)[7]肖亞男研究了幾類媒介傳染病模型的穩(wěn)定性。將傳染病模型推廣到實際應(yīng)用中的文獻(xiàn)也并不鮮見,例如趙劍華、萬克文等人基于信息傳播模型-SIR傳染病模型研究了社交網(wǎng)絡(luò)輿情傳播動力學(xué)模型[8]。趙鵬飛,馬民,談依箴等作者曾基于SIR模型研究微信輿情傳播[9]。王長峰,莊文英,于長鉞等人基于改進(jìn)SIR模型研究了的群體意見競爭演化[10]。雖然他們在研究過程中也有采用相關(guān)數(shù)學(xué)類軟件進(jìn)行問題分析,但利用數(shù)字仿真技術(shù)對傳染病模型進(jìn)行系統(tǒng)地仿真還較為鮮見。故本文首先建立傳統(tǒng)的傳染病模型——SIR模型,并建立在SIR模型中加入政府宣傳防控工作效用的傳染病模型,然后利用數(shù)字仿真技術(shù)對傳染病模型進(jìn)行系統(tǒng)地仿真,得到傳染病傳染過程的仿真結(jié)果,方便人們進(jìn)一步對傳染病的傳染過程進(jìn)行研究并推廣到實際應(yīng)用中。

    1 傳染病模型的建立

    1.1 SIR傳染病模型[9]

    假設(shè)傳染病考察區(qū)域的總?cè)藬?shù)為N,由于許多患者在傳染病治愈后會對該病產(chǎn)生很強(qiáng)的免疫力,所以假設(shè)治愈者已經(jīng)退出感染系統(tǒng),故將人群分為健康者、病人和治愈者,則t時刻的這三類人群所占的比例分別為s(t)、i(t)、r(t),病人傳染給健康者的傳染率為?姿,日治愈率為?滋。設(shè)t=0時健康者、病人和治愈者所占比例分別為s0、i0、r0=0,于是便有微分方程的初值問題:

    (1)

    其中s(t)+i(t)+r(t)=1。

    此傳染病模型已是當(dāng)前發(fā)展較為成熟的數(shù)學(xué)模型,也是比較傳統(tǒng)典型的傳染病微分方程模型,對傳染病的傳染過程能起到抽象描述作用。但在實際生活中,若突發(fā)傳染病,政府為減少健康者染病的幾率,會采取相應(yīng)的宣傳與防控措施。本文通過查閱相關(guān)資料[11],對傳染病傳染過程的仿真做進(jìn)一步優(yōu)化,建立了加入宣傳防控工作的傳染病模型。

    1.2 加入宣傳防控工作的傳染病模型

    在SIR傳染病模型的基礎(chǔ)上,加入政府的宣傳防控工作效用,假設(shè)宣傳防控工作的開展將使病人人數(shù)減少,減少的速度與傳染病考察區(qū)域的總?cè)藬?shù)N成正比,其比例常數(shù)a取決于宣傳防控的強(qiáng)度,故a也被稱為宣傳防控強(qiáng)度系數(shù)。設(shè)從t=t0>0開始做持續(xù)的宣傳防控工作,于是便有微分方程的初值問題:

    (2)

    其中 稱為Heaviside函數(shù)[12]。

    此傳染病模型在SIR模型的基礎(chǔ)上加入宣傳防控工作的效用,既考慮病人被治愈的情況,又考慮到政府的宣傳防控工作,從而建立起更加符合實際情況的傳染病模型,因而有助于本文對傳染病的傳染過程進(jìn)行進(jìn)一步的仿真,使得仿真結(jié)果更加真實可信。

    2 仿真系統(tǒng)的建立

    2.1 基于SIR模型的傳染過程仿真系統(tǒng)

    利用MATLAB軟件中的Simulink仿真軟件對SIR傳染病模型進(jìn)行仿真。此仿真系統(tǒng)通過人為設(shè)置日傳染率a_1(即?姿)和日治愈率b(即?滋)的值、病人和健康者所占比例的初值,然后進(jìn)行傳染病傳染過程的仿真,得到滿足公式(1)所建立的SIR傳染病模型,最終通過運行仿真系統(tǒng)得到基于SIR模型的傳染病傳染過程的仿真結(jié)果??偡抡嫦到y(tǒng)圖如圖1。

    2.2 加入宣傳防控工作效用的SIR傳染病模型仿真系統(tǒng)

    同理,利用Simulink軟件對加入政府宣傳防控工作效用的SIR模型進(jìn)行仿真,得到滿足公式(2)所建立模型的仿真系統(tǒng),通過設(shè)置日傳染率a_1(即?姿)和日治愈率b(即?滋)的值、病人和健康者所占比例的合適初值,計算t時刻病人所占的比例,再將計算所得的t時刻病人比率轉(zhuǎn)換為病人人數(shù),最終得到加入防控宣傳效用的SIR傳染病模型的仿真結(jié)果??偡抡嫦到y(tǒng)圖如圖2。

    3 仿真結(jié)果

    3.1 SIR傳染病模型仿真結(jié)果

    設(shè)置第一組參數(shù):?姿=0.5 ?滋=0.05 i0=0.02 S0=0.98,三類人口所占比例的仿真結(jié)果如下圖3所示。設(shè)置第二組參數(shù):?姿=0.5 ?滋=0.2 i0=0.02 S0=0.98,三類人口所占比例的仿真結(jié)果如下圖4所示。

    在此傳染病模型中,設(shè)置的兩組參數(shù)中僅有治愈率這一個參數(shù)不同,其余各參數(shù)值均相同,因此圖3、圖4起相互對照作用。故由上面的仿真結(jié)果可知,當(dāng)治愈率越高,健康者所占比例下降得越慢,病人在染病后便會更快被完全治愈,成為治愈者。因此治愈者所占比例增加速率就越快,且當(dāng)t→∞時,病人人數(shù)所占比例趨于0,符合現(xiàn)實生活中的實際情況,達(dá)到了對SIR傳染病模型傳染過程仿真的效果。

    3.2 加入宣傳防控的傳染病模型仿真結(jié)果

    設(shè)置參數(shù):?姿=0.5 S0=0.98 ?滋=0.05 i0=0.02 a=0.2 t=20 N=10000,仿真結(jié)果如下圖5所示。

    在Heaviside函數(shù)中,本文設(shè)置在t=20時加入宣傳和防控工作。按照Heaviside函數(shù)的定義[12],當(dāng)020時,由于已經(jīng)實行宣傳與防控工作,產(chǎn)生了宣傳與防控工作的效用,即Heaviside函數(shù)值為1,即T=1。此模型設(shè)置的參數(shù)與SIR傳染病模型仿真的第一組參數(shù)是相同的,并且在t=20時,加入宣傳防控系數(shù)a=0.2作用于SIR傳染病模型,病人人數(shù)快速下降,以此與SIR傳染病模型形成對比,更加符合生活中的實際情況,并更有利于實際應(yīng)用。

    在對傳統(tǒng)的SIR傳染病模型進(jìn)行仿真時,設(shè)置了不同的治愈率,結(jié)果表明在治愈率高時病人人數(shù)比例下降速度明顯高于治愈率低時的情況,該仿真結(jié)果符合實際,驗證了仿真系統(tǒng)建立的正確性。在對加入宣傳防控的傳染病模型進(jìn)行仿真時,本文設(shè)置了相同的治愈率,來對比加入了宣傳防控工作效用與沒有加入宣傳防控效用時的傳染病模型的病人人數(shù)變化趨勢,從圖5來看,加入宣傳防控工作效用后病人人數(shù)明顯下降,與實際結(jié)果吻合。

    4 結(jié)束語

    利用數(shù)字仿真技術(shù)直觀、形象地反映傳染病傳染機(jī)制,從而對更好地理解傳染病模型起到幫助作用?,F(xiàn)代仿真技術(shù)不僅在傳統(tǒng)的工程領(lǐng)域上應(yīng)用廣泛,而且正日益廣泛地應(yīng)用于如交通管理、環(huán)境治理、城市規(guī)劃等領(lǐng)域,尤其對于制造業(yè)利用數(shù)字仿真技術(shù)具有更為重要的意義[3]。通過利用數(shù)字仿真技術(shù)對傳染病模型進(jìn)行仿真,掌握傳染病模型的數(shù)學(xué)規(guī)律,可更加有利于對傳染病的預(yù)防控制、經(jīng)濟(jì)行為以及其他社會實際問題提出實際建議。

    參考文獻(xiàn):

    [1]黃勇,孟慶鑫.淺談仿真技術(shù)在軍事領(lǐng)域內(nèi)的應(yīng)用[J].計算機(jī)與網(wǎng)絡(luò),2015,41(06):63-67.

    [2]彭沖.電力系統(tǒng)數(shù)字仿真技術(shù)現(xiàn)狀及展望[J].河北電力技術(shù),2013,32(04):23-25+54.

    [3]丁明瑞,施松新,周華民,等.基于數(shù)字仿真技術(shù)的家具產(chǎn)品的虛擬展示[J].機(jī)電工程技術(shù),2014,43(12):183-185.

    [4]紀(jì)志成,沈艷霞,姜建國.基于MATLAB無刷直流電機(jī)系統(tǒng)仿真建模的新方法[J].系統(tǒng)仿真學(xué)報,2003(12):1745-1749+1758.

    [5]張?zhí)r,安連鎖,劉樹華.基于Simulink仿真實現(xiàn)換熱器穩(wěn)態(tài)傳熱計算[J].水利電力機(jī)械,2006(08):30-31+46.

    [6]劉玉英,肖燕妮.一類受媒體影響的傳染病模型的研究[J].應(yīng)用數(shù)學(xué)和力學(xué),2013,34(04):399-407.

    [7]肖亞男.幾類媒介傳染病模型的穩(wěn)定性分析[D].中北大學(xué),2015.

    [8]趙劍華,萬克文.基于信息傳播模型-SIR傳染病模型的社交網(wǎng)絡(luò)輿情傳播動力學(xué)模型研究[J].情報科學(xué),2017,35(12):34-38.

    [9]趙鵬飛,馬民,談依箴.基于SIR模型的微信輿情傳播研究[J].情報探索,2017(10):1-6.

    [10]王長峰,莊文英,于長鉞.基于改進(jìn)SIR模型的群體意見競爭演化研究[J].情報雜志,2017,36(10):97-103+136.

    [11]李漢龍,艾瑛,韓婷,等.Mathematica基礎(chǔ)及其在數(shù)學(xué)建模中的應(yīng)用[M].北京:國防工業(yè)出版社,2016:221-225.

    [12]王海林,王吉民.Heaviside函數(shù)的Laplace變換建立超靜定梁撓曲線方程[J].浙江科技學(xué)院學(xué)報,2012,24(06):466-469+487.

    国产成人aa在线观看| 一级黄片播放器| 中文字幕免费在线视频6| 亚洲丝袜综合中文字幕| 国产免费福利视频在线观看| 久久久久久久久大av| 天美传媒精品一区二区| 久久精品夜色国产| 久热这里只有精品99| 人人妻人人澡人人爽人人夜夜| av国产免费在线观看| 18禁在线播放成人免费| 婷婷色av中文字幕| 国产欧美亚洲国产| 一区二区av电影网| 最后的刺客免费高清国语| 熟妇人妻不卡中文字幕| 国产精品一区二区在线观看99| 国产精品无大码| 精品一区在线观看国产| 国产亚洲一区二区精品| 精品视频人人做人人爽| 91午夜精品亚洲一区二区三区| 自拍欧美九色日韩亚洲蝌蚪91 | 在线观看av片永久免费下载| 成人毛片a级毛片在线播放| 22中文网久久字幕| 天天躁夜夜躁狠狠久久av| 日韩人妻高清精品专区| 天天躁夜夜躁狠狠久久av| 黄片wwwwww| 99热网站在线观看| 久久久精品94久久精品| av国产精品久久久久影院| 禁无遮挡网站| 麻豆乱淫一区二区| 熟妇人妻不卡中文字幕| 国产免费又黄又爽又色| 国产精品99久久久久久久久| av福利片在线观看| 亚洲精品久久久久久婷婷小说| 亚洲国产av新网站| 777米奇影视久久| 少妇猛男粗大的猛烈进出视频 | 亚洲国产精品专区欧美| 午夜福利视频1000在线观看| 国产探花在线观看一区二区| 久久久亚洲精品成人影院| 精品久久久久久久久av| 亚洲欧美一区二区三区国产| 日韩在线高清观看一区二区三区| 精品人妻视频免费看| 国产精品女同一区二区软件| 亚洲av不卡在线观看| www.av在线官网国产| 亚洲天堂av无毛| 久久综合国产亚洲精品| 简卡轻食公司| 少妇的逼水好多| 国产永久视频网站| 成人漫画全彩无遮挡| 韩国av在线不卡| 人人妻人人看人人澡| 亚洲精品日韩在线中文字幕| 精品少妇久久久久久888优播| 国产一级毛片在线| av天堂中文字幕网| 大片电影免费在线观看免费| 最近的中文字幕免费完整| 久久久精品免费免费高清| 涩涩av久久男人的天堂| 免费av不卡在线播放| 精品国产三级普通话版| 婷婷色综合www| 日本一本二区三区精品| 男女边摸边吃奶| 精品一区在线观看国产| 中文字幕免费在线视频6| 亚洲成色77777| 精品久久久久久电影网| 国产老妇伦熟女老妇高清| 欧美97在线视频| 久久热精品热| 日韩成人av中文字幕在线观看| 亚洲av不卡在线观看| 看十八女毛片水多多多| 午夜激情久久久久久久| 在线免费十八禁| 午夜视频国产福利| 免费高清在线观看视频在线观看| av在线播放精品| 内地一区二区视频在线| 三级男女做爰猛烈吃奶摸视频| 亚洲av.av天堂| 午夜福利视频1000在线观看| 免费看不卡的av| 亚洲精品,欧美精品| 色视频在线一区二区三区| 欧美 日韩 精品 国产| 久久鲁丝午夜福利片| 日韩伦理黄色片| 国产免费一区二区三区四区乱码| 男女无遮挡免费网站观看| 久久久亚洲精品成人影院| 国产视频首页在线观看| 成年av动漫网址| 插逼视频在线观看| 激情五月婷婷亚洲| 免费黄网站久久成人精品| 免费av不卡在线播放| 各种免费的搞黄视频| 精品一区二区三区视频在线| 国产欧美日韩精品一区二区| 亚洲精品成人久久久久久| 国产成人a∨麻豆精品| 直男gayav资源| av.在线天堂| 高清av免费在线| 久久精品夜色国产| 波野结衣二区三区在线| 久久久久久久午夜电影| 国产精品久久久久久久久免| 国内精品美女久久久久久| 成人特级av手机在线观看| 干丝袜人妻中文字幕| 色综合色国产| 国产成人福利小说| 美女高潮的动态| 五月开心婷婷网| 久久久久久伊人网av| 老司机影院成人| 国产精品人妻久久久久久| 国产真实伦视频高清在线观看| 久久6这里有精品| 伊人久久国产一区二区| 神马国产精品三级电影在线观看| 亚洲婷婷狠狠爱综合网| 亚洲人成网站在线观看播放| 99热国产这里只有精品6| 成人无遮挡网站| 大又大粗又爽又黄少妇毛片口| 神马国产精品三级电影在线观看| 午夜老司机福利剧场| 国产淫语在线视频| 在线观看一区二区三区| 水蜜桃什么品种好| 国产男女超爽视频在线观看| 久久ye,这里只有精品| 亚洲国产成人一精品久久久| 国产高清不卡午夜福利| 激情五月婷婷亚洲| 久久精品国产a三级三级三级| 1000部很黄的大片| 卡戴珊不雅视频在线播放| 最近最新中文字幕大全电影3| 国产午夜福利久久久久久| 日韩中字成人| 男女边吃奶边做爰视频| 男男h啪啪无遮挡| 成人无遮挡网站| 干丝袜人妻中文字幕| 成人欧美大片| 九九在线视频观看精品| 久久精品久久久久久噜噜老黄| 少妇的逼好多水| 男插女下体视频免费在线播放| 狠狠精品人妻久久久久久综合| 18+在线观看网站| av网站免费在线观看视频| 日韩大片免费观看网站| 91精品伊人久久大香线蕉| 欧美3d第一页| 看免费成人av毛片| 亚洲欧洲日产国产| 国产在线一区二区三区精| 国产爱豆传媒在线观看| 免费观看性生交大片5| 国产 一区 欧美 日韩| 国内少妇人妻偷人精品xxx网站| 免费大片18禁| 男人添女人高潮全过程视频| 国产真实伦视频高清在线观看| 精品久久久精品久久久| 久久久欧美国产精品| 26uuu在线亚洲综合色| 色婷婷久久久亚洲欧美| 国产探花极品一区二区| 一边亲一边摸免费视频| 丝瓜视频免费看黄片| 麻豆久久精品国产亚洲av| 少妇裸体淫交视频免费看高清| 亚洲国产精品国产精品| 国产精品无大码| 大话2 男鬼变身卡| 亚洲精品影视一区二区三区av| 国产探花在线观看一区二区| 国内揄拍国产精品人妻在线| 国产精品国产三级国产专区5o| 亚洲成色77777| 久久精品综合一区二区三区| 99热这里只有精品一区| 男人舔奶头视频| av卡一久久| 秋霞伦理黄片| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 大片免费播放器 马上看| 成人鲁丝片一二三区免费| 老司机影院毛片| 麻豆国产97在线/欧美| 国产精品久久久久久精品古装| 欧美三级亚洲精品| 亚洲aⅴ乱码一区二区在线播放| 联通29元200g的流量卡| 一级毛片 在线播放| 三级国产精品片| 插逼视频在线观看| 欧美一区二区亚洲| 国产探花在线观看一区二区| 99热网站在线观看| 中文字幕亚洲精品专区| 午夜福利在线在线| 简卡轻食公司| 成年女人看的毛片在线观看| 精品久久久久久久末码| 亚洲,欧美,日韩| 亚州av有码| 婷婷色av中文字幕| 日韩一区二区视频免费看| 亚洲自拍偷在线| 午夜福利视频1000在线观看| 一个人看的www免费观看视频| 国产亚洲最大av| 高清在线视频一区二区三区| 欧美bdsm另类| 色视频www国产| 亚洲电影在线观看av| 青青草视频在线视频观看| 国产成人精品一,二区| 亚洲欧美精品专区久久| 亚洲成人久久爱视频| 国产精品一区二区在线观看99| 综合色av麻豆| 免费人成在线观看视频色| 国内揄拍国产精品人妻在线| 亚洲精品456在线播放app| 1000部很黄的大片| 免费观看的影片在线观看| 亚洲性久久影院| 亚洲精品乱久久久久久| 熟妇人妻不卡中文字幕| 亚洲av中文av极速乱| 日本一二三区视频观看| 成人亚洲欧美一区二区av| 国产精品精品国产色婷婷| 观看免费一级毛片| 国产伦精品一区二区三区视频9| 水蜜桃什么品种好| 国产熟女欧美一区二区| 久久午夜福利片| 亚洲无线观看免费| 国产男人的电影天堂91| 蜜桃久久精品国产亚洲av| 免费高清在线观看视频在线观看| 最近2019中文字幕mv第一页| 亚洲精品aⅴ在线观看| 91精品伊人久久大香线蕉| 91精品一卡2卡3卡4卡| 亚洲自偷自拍三级| 国产色爽女视频免费观看| 亚洲熟女精品中文字幕| 久久精品久久久久久久性| 亚洲欧美一区二区三区黑人 | 91aial.com中文字幕在线观看| 在线 av 中文字幕| 成年女人看的毛片在线观看| videossex国产| 在线精品无人区一区二区三 | 一级黄片播放器| 国产免费福利视频在线观看| 午夜福利视频1000在线观看| 精品国产露脸久久av麻豆| 久久精品久久精品一区二区三区| 高清视频免费观看一区二区| 一级二级三级毛片免费看| 一级毛片aaaaaa免费看小| 美女国产视频在线观看| 日本三级黄在线观看| 国产精品无大码| 最近的中文字幕免费完整| 天堂中文最新版在线下载 | 内射极品少妇av片p| 国产一区二区亚洲精品在线观看| 精品一区二区三卡| 国产亚洲91精品色在线| 精品人妻一区二区三区麻豆| 精品久久久久久久久亚洲| 青春草国产在线视频| 国产永久视频网站| 99久国产av精品国产电影| av天堂中文字幕网| 欧美xxxx黑人xx丫x性爽| 亚洲在久久综合| 久久ye,这里只有精品| 老女人水多毛片| 国产乱人视频| 精品少妇久久久久久888优播| 性色av一级| 久久亚洲国产成人精品v| 少妇人妻一区二区三区视频| 欧美zozozo另类| 中文在线观看免费www的网站| 国产老妇女一区| 国产黄a三级三级三级人| 在线免费观看不下载黄p国产| 日韩国内少妇激情av| 国产成人免费观看mmmm| 国产老妇伦熟女老妇高清| 欧美日本视频| 国产伦精品一区二区三区视频9| 在现免费观看毛片| 国产又色又爽无遮挡免| 成人国产麻豆网| 国产片特级美女逼逼视频| 老司机影院毛片| 一区二区三区四区激情视频| 伊人久久精品亚洲午夜| 黄色怎么调成土黄色| 97超碰精品成人国产| 精品久久久久久久久av| 人妻 亚洲 视频| 亚洲国产欧美在线一区| 国产成人精品一,二区| 一级毛片 在线播放| 日本三级黄在线观看| 26uuu在线亚洲综合色| 久热久热在线精品观看| 日韩中字成人| 欧美高清性xxxxhd video| 久久国产乱子免费精品| 晚上一个人看的免费电影| 街头女战士在线观看网站| 少妇人妻精品综合一区二区| 日日啪夜夜撸| 天堂俺去俺来也www色官网| 久久久亚洲精品成人影院| 国产精品久久久久久精品电影小说 | 国产探花在线观看一区二区| 国产一区二区亚洲精品在线观看| 久久鲁丝午夜福利片| 综合色丁香网| 一个人看视频在线观看www免费| 在线观看国产h片| 精品久久久久久久人妻蜜臀av| 色哟哟·www| 午夜精品一区二区三区免费看| 国产精品国产三级国产av玫瑰| 97在线人人人人妻| 一本色道久久久久久精品综合| 国产日韩欧美亚洲二区| 亚洲精品乱久久久久久| 男人添女人高潮全过程视频| 欧美国产精品一级二级三级 | 亚洲美女搞黄在线观看| 一区二区av电影网| 亚洲国产精品国产精品| 国产精品国产三级国产av玫瑰| 亚洲av不卡在线观看| 少妇被粗大猛烈的视频| 亚洲精品乱码久久久v下载方式| 又大又黄又爽视频免费| av在线app专区| 最新中文字幕久久久久| 伊人久久国产一区二区| 亚洲自偷自拍三级| 日韩欧美一区视频在线观看 | 日韩欧美 国产精品| 美女cb高潮喷水在线观看| 在线天堂最新版资源| 午夜福利视频1000在线观看| 亚洲精品aⅴ在线观看| 国产精品一区www在线观看| 亚洲av欧美aⅴ国产| 熟女人妻精品中文字幕| 日韩av在线免费看完整版不卡| 在线免费观看不下载黄p国产| 午夜福利在线在线| 午夜激情福利司机影院| 男人舔奶头视频| 热99国产精品久久久久久7| 插逼视频在线观看| 亚洲国产色片| 久久99精品国语久久久| 干丝袜人妻中文字幕| 亚洲天堂av无毛| 亚洲av不卡在线观看| 精品国产乱码久久久久久小说| 国产91av在线免费观看| www.色视频.com| 精品人妻一区二区三区麻豆| 亚洲精品色激情综合| 夜夜爽夜夜爽视频| 狂野欧美激情性xxxx在线观看| 另类亚洲欧美激情| 免费大片18禁| 日本熟妇午夜| 天天躁日日操中文字幕| 九九在线视频观看精品| 日韩亚洲欧美综合| 91午夜精品亚洲一区二区三区| 乱系列少妇在线播放| 婷婷色麻豆天堂久久| 一级黄片播放器| 国产午夜精品久久久久久一区二区三区| 韩国高清视频一区二区三区| 亚洲欧美日韩无卡精品| 蜜桃久久精品国产亚洲av| 九九久久精品国产亚洲av麻豆| 国产大屁股一区二区在线视频| 国产精品国产三级国产av玫瑰| 国产精品一区二区三区四区免费观看| 亚洲美女视频黄频| 国国产精品蜜臀av免费| 美女内射精品一级片tv| 久久久久久久久久成人| 777米奇影视久久| 亚洲精品456在线播放app| 久久女婷五月综合色啪小说 | 日产精品乱码卡一卡2卡三| 听说在线观看完整版免费高清| 一级爰片在线观看| 看非洲黑人一级黄片| 狂野欧美激情性bbbbbb| 全区人妻精品视频| 视频中文字幕在线观看| 可以在线观看毛片的网站| 99久久精品热视频| 美女视频免费永久观看网站| 久久久久久九九精品二区国产| 一区二区三区免费毛片| 日本熟妇午夜| 精品一区在线观看国产| 国产探花在线观看一区二区| 亚洲av成人精品一二三区| 亚洲精品成人av观看孕妇| 制服丝袜香蕉在线| 一区二区三区四区激情视频| 久久久欧美国产精品| 国产成人精品婷婷| 国产日韩欧美在线精品| 久久久成人免费电影| 中文字幕人妻熟人妻熟丝袜美| 国产成年人精品一区二区| 久久这里有精品视频免费| 国产成人精品一,二区| 一级片'在线观看视频| 欧美变态另类bdsm刘玥| 观看美女的网站| 三级男女做爰猛烈吃奶摸视频| 国产亚洲av嫩草精品影院| 国产爱豆传媒在线观看| 婷婷色av中文字幕| 精品一区二区免费观看| 日韩一区二区视频免费看| 听说在线观看完整版免费高清| 黄色配什么色好看| 日韩免费高清中文字幕av| 欧美精品人与动牲交sv欧美| 亚洲最大成人手机在线| 老司机影院毛片| 久久久久久久午夜电影| 在线观看国产h片| 亚洲人成网站高清观看| 国产一区二区三区av在线| 97精品久久久久久久久久精品| 日本黄大片高清| 毛片一级片免费看久久久久| 欧美激情国产日韩精品一区| 日韩强制内射视频| 国产色婷婷99| 成人国产av品久久久| 欧美日韩在线观看h| 少妇熟女欧美另类| 国产av码专区亚洲av| 亚洲高清免费不卡视频| 九草在线视频观看| 久久精品久久精品一区二区三区| 欧美老熟妇乱子伦牲交| 久久久久久久午夜电影| 男人添女人高潮全过程视频| 国产成年人精品一区二区| 亚洲三级黄色毛片| 日韩av免费高清视频| 欧美高清性xxxxhd video| 中国三级夫妇交换| 日韩成人伦理影院| 免费观看无遮挡的男女| 永久免费av网站大全| 干丝袜人妻中文字幕| 我的女老师完整版在线观看| 国产精品人妻久久久久久| av卡一久久| 精品人妻偷拍中文字幕| 18禁动态无遮挡网站| xxx大片免费视频| 人妻一区二区av| 能在线免费看毛片的网站| 亚洲精品自拍成人| 亚洲国产精品国产精品| 国产亚洲最大av| 亚洲成色77777| 国产精品99久久99久久久不卡 | 大陆偷拍与自拍| 亚洲精品日本国产第一区| 亚洲熟女精品中文字幕| 丝袜喷水一区| av在线天堂中文字幕| 国产精品一区二区在线观看99| 国产精品久久久久久精品电影小说 | 男女下面进入的视频免费午夜| 久久久久久久亚洲中文字幕| 色视频在线一区二区三区| av.在线天堂| 麻豆乱淫一区二区| av免费在线看不卡| kizo精华| 亚洲综合精品二区| 日本与韩国留学比较| 亚洲精品一区蜜桃| 综合色av麻豆| 欧美最新免费一区二区三区| 免费黄色在线免费观看| 一区二区av电影网| 80岁老熟妇乱子伦牲交| 国产免费一区二区三区四区乱码| 国产极品天堂在线| 久久久久久久久久久丰满| 精品一区在线观看国产| .国产精品久久| 韩国av在线不卡| 我的女老师完整版在线观看| 在线播放无遮挡| 亚洲精品国产av成人精品| 黄片无遮挡物在线观看| 中文字幕人妻熟人妻熟丝袜美| 国产精品精品国产色婷婷| 最近中文字幕2019免费版| 国产精品一二三区在线看| 国产成人免费观看mmmm| 日日摸夜夜添夜夜添av毛片| 久久国内精品自在自线图片| av免费观看日本| 国产色婷婷99| 色播亚洲综合网| 国产爱豆传媒在线观看| 九九爱精品视频在线观看| 国产视频首页在线观看| 夜夜看夜夜爽夜夜摸| 日韩一区二区三区影片| 各种免费的搞黄视频| 欧美+日韩+精品| 国产精品嫩草影院av在线观看| 黄片wwwwww| 可以在线观看毛片的网站| 大片免费播放器 马上看| 蜜臀久久99精品久久宅男| 日日啪夜夜爽| 视频中文字幕在线观看| 欧美精品一区二区大全| 成人亚洲精品av一区二区| 欧美另类一区| 九九在线视频观看精品| 欧美丝袜亚洲另类| xxx大片免费视频| 精品酒店卫生间| 亚洲人与动物交配视频| 一级毛片 在线播放| 亚洲无线观看免费| 国产老妇伦熟女老妇高清| 亚洲欧洲国产日韩| 另类亚洲欧美激情| 男插女下体视频免费在线播放| 97人妻精品一区二区三区麻豆| 国产成人精品一,二区| 国产成人免费无遮挡视频| freevideosex欧美| 丝袜喷水一区| 欧美zozozo另类| 2022亚洲国产成人精品| 免费黄频网站在线观看国产| 免费观看性生交大片5| 男女啪啪激烈高潮av片| 97在线人人人人妻| 久久99热这里只频精品6学生| 天堂中文最新版在线下载 | 日日啪夜夜撸| 免费少妇av软件| 六月丁香七月| 国产午夜精品久久久久久一区二区三区| 亚洲在线观看片| 免费观看无遮挡的男女| 国产极品天堂在线| 亚洲人成网站高清观看| 国产成人a∨麻豆精品| 久久99蜜桃精品久久| 成人一区二区视频在线观看| 色视频在线一区二区三区| 天天躁夜夜躁狠狠久久av| 国产 一区 欧美 日韩| 青青草视频在线视频观看| 午夜福利高清视频| 老女人水多毛片| 啦啦啦中文免费视频观看日本| 男人舔奶头视频| 黄色欧美视频在线观看| 国产成人a∨麻豆精品|