• <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期

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

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

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

    中圖分類號:TP391.9 文獻標(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ǔ),以計算機和各種物理效應(yīng)設(shè)備為工具,利用系統(tǒng)模型對實際的或設(shè)想的系統(tǒng)進行試驗研究的一種綜合性技術(shù)[1]。與物理仿真技術(shù)相比,數(shù)字仿真不受被研究系統(tǒng)規(guī)模和結(jié)構(gòu)復(fù)雜性的限制,以計算速度快、使用靈活、成本相對低廉等優(yōu)勢,成為分析、研究電力系統(tǒng)必不可少的工具,是系統(tǒng)仿真的重要研究方向[2]。隨著計算機的發(fā)展,數(shù)字仿真技術(shù)取得突飛猛進的進步[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仿真軟件對無刷直流電機控制系統(tǒng)建模進行仿真[4],實現(xiàn)換熱器穩(wěn)態(tài)傳熱計算[5]。而傳染病模型作為一類發(fā)展十分成熟的數(shù)學(xué)模型,相關(guān)的研究文獻數(shù)不勝數(shù)。文獻[6]劉玉英,肖燕妮等人研究了一類受媒體影響的傳染病模型,文獻[7]肖亞男研究了幾類媒介傳染病模型的穩(wěn)定性。將傳染病模型推廣到實際應(yīng)用中的文獻也并不鮮見,例如趙劍華、萬克文等人基于信息傳播模型-SIR傳染病模型研究了社交網(wǎng)絡(luò)輿情傳播動力學(xué)模型[8]。趙鵬飛,馬民,談依箴等作者曾基于SIR模型研究微信輿情傳播[9]。王長峰,莊文英,于長鉞等人基于改進SIR模型研究了的群體意見競爭演化[10]。雖然他們在研究過程中也有采用相關(guān)數(shù)學(xué)類軟件進行問題分析,但利用數(shù)字仿真技術(shù)對傳染病模型進行系統(tǒng)地仿真還較為鮮見。故本文首先建立傳統(tǒng)的傳染病模型——SIR模型,并建立在SIR模型中加入政府宣傳防控工作效用的傳染病模型,然后利用數(shù)字仿真技術(shù)對傳染病模型進行系統(tǒng)地仿真,得到傳染病傳染過程的仿真結(jié)果,方便人們進一步對傳染病的傳染過程進行研究并推廣到實際應(yīng)用中。

    1 傳染病模型的建立

    1.1 SIR傳染病模型[9]

    假設(shè)傳染病考察區(qū)域的總?cè)藬?shù)為N,由于許多患者在傳染病治愈后會對該病產(chǎn)生很強的免疫力,所以假設(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],對傳染病傳染過程的仿真做進一步優(yōu)化,建立了加入宣傳防控工作的傳染病模型。

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

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

    (2)

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

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

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

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

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

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

    同理,利用Simulink軟件對加入政府宣傳防控工作效用的SIR模型進行仿真,得到滿足公式(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)實生活中的實際情況,達到了對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傳染病模型進行仿真時,設(shè)置了不同的治愈率,結(jié)果表明在治愈率高時病人人數(shù)比例下降速度明顯高于治愈率低時的情況,該仿真結(jié)果符合實際,驗證了仿真系統(tǒng)建立的正確性。在對加入宣傳防控的傳染病模型進行仿真時,本文設(shè)置了相同的治愈率,來對比加入了宣傳防控工作效用與沒有加入宣傳防控效用時的傳染病模型的病人人數(shù)變化趨勢,從圖5來看,加入宣傳防控工作效用后病人人數(shù)明顯下降,與實際結(jié)果吻合。

    4 結(jié)束語

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

    參考文獻:

    [1]黃勇,孟慶鑫.淺談仿真技術(shù)在軍事領(lǐng)域內(nèi)的應(yīng)用[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].機電工程技術(shù),2014,43(12):183-185.

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

    [5]張?zhí)r,安連鎖,劉樹華.基于Simulink仿真實現(xiàn)換熱器穩(wěn)態(tài)傳熱計算[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]王長峰,莊文英,于長鉞.基于改進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.

    bbb黄色大片| 日韩欧美国产一区二区入口| 十八禁网站免费在线| 久久久久久久精品吃奶| 老司机福利观看| 香蕉久久夜色| 久久久国产欧美日韩av| 一个人免费在线观看的高清视频| 很黄的视频免费| 啦啦啦观看免费观看视频高清| 三级男女做爰猛烈吃奶摸视频| 97人妻精品一区二区三区麻豆| 欧美日韩中文字幕国产精品一区二区三区| 日韩有码中文字幕| 亚洲av成人一区二区三| 波多野结衣高清无吗| 哪里可以看免费的av片| 麻豆国产av国片精品| 制服人妻中文乱码| 国产一区二区三区视频了| 1024香蕉在线观看| 欧美乱妇无乱码| 国产aⅴ精品一区二区三区波| 成年版毛片免费区| 他把我摸到了高潮在线观看| 男人的好看免费观看在线视频| 欧洲精品卡2卡3卡4卡5卡区| 精品福利观看| 国产欧美日韩一区二区三| 性色av乱码一区二区三区2| 国产欧美日韩一区二区精品| 又爽又黄无遮挡网站| 麻豆一二三区av精品| 91九色精品人成在线观看| 成人特级黄色片久久久久久久| 黄色女人牲交| 男插女下体视频免费在线播放| 在线播放国产精品三级| 999久久久精品免费观看国产| 国产精品一及| 1024手机看黄色片| 亚洲中文日韩欧美视频| 欧美日本亚洲视频在线播放| 岛国在线免费视频观看| 亚洲国产精品久久男人天堂| 日韩中文字幕欧美一区二区| 亚洲熟妇熟女久久| 国产成人av教育| 国产伦精品一区二区三区四那| 两个人的视频大全免费| 国产91精品成人一区二区三区| 久久天躁狠狠躁夜夜2o2o| 99热这里只有精品一区 | 丰满人妻一区二区三区视频av | 一级作爱视频免费观看| 人人妻人人看人人澡| 好男人电影高清在线观看| 久久香蕉国产精品| 日本黄色片子视频| 黄色日韩在线| 天堂动漫精品| 国产淫片久久久久久久久 | tocl精华| 国产成人精品久久二区二区免费| 在线视频色国产色| www.www免费av| 国产综合懂色| 精品电影一区二区在线| 91av网站免费观看| 特大巨黑吊av在线直播| 国产亚洲av高清不卡| 国产精品一及| 三级毛片av免费| 久久久久久久久中文| 97碰自拍视频| 中文字幕熟女人妻在线| 婷婷丁香在线五月| 亚洲成人免费电影在线观看| 亚洲国产欧美人成| 欧美中文日本在线观看视频| 少妇熟女aⅴ在线视频| 国产一区二区在线观看日韩 | 好男人在线观看高清免费视频| 国产探花在线观看一区二区| 男女午夜视频在线观看| 亚洲国产精品成人综合色| 一卡2卡三卡四卡精品乱码亚洲| 曰老女人黄片| 国产精品免费一区二区三区在线| 两性夫妻黄色片| 人妻丰满熟妇av一区二区三区| 我的老师免费观看完整版| 哪里可以看免费的av片| 国产成人欧美在线观看| av福利片在线观看| 九色成人免费人妻av| 国产成人一区二区三区免费视频网站| 搞女人的毛片| 一进一出好大好爽视频| 午夜福利视频1000在线观看| 国产欧美日韩精品一区二区| 亚洲国产精品999在线| 两个人视频免费观看高清| 亚洲男人的天堂狠狠| 亚洲专区国产一区二区| 九色成人免费人妻av| x7x7x7水蜜桃| 一个人看的www免费观看视频| 91麻豆精品激情在线观看国产| 禁无遮挡网站| 亚洲成av人片在线播放无| 免费看a级黄色片| 国产一区在线观看成人免费| 高清毛片免费观看视频网站| 小蜜桃在线观看免费完整版高清| 美女午夜性视频免费| 免费搜索国产男女视频| 九色国产91popny在线| 一进一出抽搐gif免费好疼| 国产亚洲欧美在线一区二区| 91麻豆av在线| 国产一级毛片七仙女欲春2| 亚洲精品色激情综合| 操出白浆在线播放| 婷婷亚洲欧美| 哪里可以看免费的av片| 欧美成人免费av一区二区三区| 亚洲欧美精品综合一区二区三区| 国产精品九九99| 国产69精品久久久久777片 | 99久久国产精品久久久| 搞女人的毛片| 亚洲国产欧洲综合997久久,| 国产97色在线日韩免费| 九九在线视频观看精品| 极品教师在线免费播放| av黄色大香蕉| 久久久久亚洲av毛片大全| 亚洲欧美一区二区三区黑人| 亚洲欧美日韩高清在线视频| av中文乱码字幕在线| 国产高清有码在线观看视频| 精品国产乱子伦一区二区三区| 亚洲 国产 在线| 少妇裸体淫交视频免费看高清| 亚洲avbb在线观看| 一进一出抽搐gif免费好疼| 国产精品女同一区二区软件 | 日韩三级视频一区二区三区| 婷婷丁香在线五月| 午夜视频精品福利| 欧美zozozo另类| 制服人妻中文乱码| 午夜福利18| 夜夜爽天天搞| 午夜日韩欧美国产| 免费av毛片视频| 国产成人aa在线观看| 日韩国内少妇激情av| 免费电影在线观看免费观看| 国产亚洲欧美98| 婷婷精品国产亚洲av| 丰满人妻一区二区三区视频av | 欧洲精品卡2卡3卡4卡5卡区| 亚洲av成人不卡在线观看播放网| 亚洲av熟女| 国产aⅴ精品一区二区三区波| 国产乱人视频| 午夜福利在线观看吧| 美女被艹到高潮喷水动态| 亚洲乱码一区二区免费版| 国产精品 国内视频| 欧美一级a爱片免费观看看| 久久久久亚洲av毛片大全| 老鸭窝网址在线观看| 精品国产乱子伦一区二区三区| 午夜影院日韩av| 欧美黑人巨大hd| 欧美一级a爱片免费观看看| 国产三级在线视频| 99在线视频只有这里精品首页| 中文资源天堂在线| 亚洲av成人av| 国产精品久久电影中文字幕| 国产午夜福利久久久久久| 97超视频在线观看视频| 19禁男女啪啪无遮挡网站| 18禁黄网站禁片午夜丰满| 中文字幕av在线有码专区| 久久精品人妻少妇| 亚洲av熟女| 日韩大尺度精品在线看网址| 午夜a级毛片| 免费观看的影片在线观看| 好看av亚洲va欧美ⅴa在| 久久99热这里只有精品18| 欧美日韩一级在线毛片| av天堂中文字幕网| 18禁美女被吸乳视频| 九色国产91popny在线| 最好的美女福利视频网| 亚洲欧美日韩高清专用| 国语自产精品视频在线第100页| 亚洲熟妇熟女久久| 午夜久久久久精精品| 亚洲 国产 在线| 夜夜爽天天搞| 国产成+人综合+亚洲专区| 久久久成人免费电影| 三级毛片av免费| 亚洲精品一卡2卡三卡4卡5卡| 视频区欧美日本亚洲| 欧洲精品卡2卡3卡4卡5卡区| 久久精品人妻少妇| 久久午夜亚洲精品久久| 国产精品一区二区精品视频观看| 国内毛片毛片毛片毛片毛片| 久久热在线av| 热99在线观看视频| 亚洲第一欧美日韩一区二区三区| 91字幕亚洲| 老鸭窝网址在线观看| 人人妻人人看人人澡| 色综合亚洲欧美另类图片| 在线看三级毛片| 亚洲中文字幕日韩| 欧美在线一区亚洲| 少妇的逼水好多| 99久久精品热视频| 精品日产1卡2卡| 一进一出抽搐gif免费好疼| 老鸭窝网址在线观看| 在线国产一区二区在线| 十八禁网站免费在线| 熟女少妇亚洲综合色aaa.| 日韩欧美 国产精品| 岛国在线免费视频观看| 国产伦一二天堂av在线观看| 成人特级av手机在线观看| 欧美国产日韩亚洲一区| 成人三级黄色视频| x7x7x7水蜜桃| 精品电影一区二区在线| 国产精华一区二区三区| 网址你懂的国产日韩在线| 90打野战视频偷拍视频| 可以在线观看毛片的网站| 国内揄拍国产精品人妻在线| 亚洲第一欧美日韩一区二区三区| 亚洲精品一区av在线观看| 国产精华一区二区三区| 久久国产精品影院| 夜夜爽天天搞| 欧美激情久久久久久爽电影| 叶爱在线成人免费视频播放| 小蜜桃在线观看免费完整版高清| 长腿黑丝高跟| 国产高清三级在线| 久久精品国产清高在天天线| 99热这里只有精品一区 | 99在线视频只有这里精品首页| 老司机在亚洲福利影院| av视频在线观看入口| 精品久久久久久久人妻蜜臀av| 他把我摸到了高潮在线观看| 亚洲男人的天堂狠狠| 男女午夜视频在线观看| 黄色片一级片一级黄色片| 小说图片视频综合网站| 国产午夜精品论理片| 亚洲国产欧洲综合997久久,| 男女那种视频在线观看| 18禁国产床啪视频网站| 精品久久久久久久久久久久久| 久久国产乱子伦精品免费另类| 国产淫片久久久久久久久 | 亚洲人成网站高清观看| 久久精品国产清高在天天线| 色av中文字幕| 久久国产乱子伦精品免费另类| 精品熟女少妇八av免费久了| 亚洲中文字幕一区二区三区有码在线看 | 亚洲在线观看片| 色噜噜av男人的天堂激情| 国产视频内射| 神马国产精品三级电影在线观看| 成人性生交大片免费视频hd| 黑人巨大精品欧美一区二区mp4| 搡老妇女老女人老熟妇| 麻豆一二三区av精品| 久久久久久久久久黄片| 制服人妻中文乱码| 99久久精品热视频| 日本在线视频免费播放| 男女下面进入的视频免费午夜| e午夜精品久久久久久久| 精品一区二区三区视频在线观看免费| 这个男人来自地球电影免费观看| 三级男女做爰猛烈吃奶摸视频| 亚洲精品中文字幕一二三四区| 亚洲成av人片在线播放无| 19禁男女啪啪无遮挡网站| 真实男女啪啪啪动态图| 1024香蕉在线观看| 久久午夜综合久久蜜桃| 他把我摸到了高潮在线观看| 黄色丝袜av网址大全| 欧美在线黄色| 欧美一级毛片孕妇| 亚洲国产精品sss在线观看| bbb黄色大片| 我的老师免费观看完整版| 久久久久九九精品影院| 熟女电影av网| 免费看十八禁软件| 亚洲专区国产一区二区| 成人无遮挡网站| 久久久久国产精品人妻aⅴ院| 成在线人永久免费视频| 9191精品国产免费久久| 淫秽高清视频在线观看| 午夜日韩欧美国产| 成在线人永久免费视频| 久久精品国产综合久久久| 麻豆国产97在线/欧美| 国产成人av激情在线播放| 久久精品国产清高在天天线| 久久人妻av系列| 黄色女人牲交| 中文字幕最新亚洲高清| 18禁国产床啪视频网站| 久久久久久久久免费视频了| 岛国视频午夜一区免费看| 国产人伦9x9x在线观看| 91麻豆精品激情在线观看国产| 国产久久久一区二区三区| 搡老岳熟女国产| 色视频www国产| av天堂在线播放| 黄色片一级片一级黄色片| 一进一出抽搐gif免费好疼| 真人一进一出gif抽搐免费| 动漫黄色视频在线观看| 欧美乱妇无乱码| 一区二区三区激情视频| 午夜精品在线福利| 亚洲成人久久性| 国产欧美日韩一区二区精品| 老熟妇仑乱视频hdxx| av女优亚洲男人天堂 | 久久久国产成人精品二区| 国产男靠女视频免费网站| 国产成人av教育| 91字幕亚洲| 可以在线观看的亚洲视频| 每晚都被弄得嗷嗷叫到高潮| 久久久国产成人免费| 老汉色∧v一级毛片| 亚洲色图 男人天堂 中文字幕| 亚洲中文字幕日韩| 精品欧美国产一区二区三| 成人特级黄色片久久久久久久| 俺也久久电影网| 亚洲人成网站在线播放欧美日韩| 桃红色精品国产亚洲av| 三级国产精品欧美在线观看 | 日日摸夜夜添夜夜添小说| 久久精品国产99精品国产亚洲性色| 999久久久精品免费观看国产| 禁无遮挡网站| 国产高潮美女av| 1024手机看黄色片| 又粗又爽又猛毛片免费看| 99久久国产精品久久久| 久久人人精品亚洲av| 97超级碰碰碰精品色视频在线观看| 午夜福利在线在线| 性欧美人与动物交配| 97人妻精品一区二区三区麻豆| 久久久久国产一级毛片高清牌| 少妇人妻一区二区三区视频| 91麻豆av在线| 在线永久观看黄色视频| 最近最新中文字幕大全电影3| 1024香蕉在线观看| 一级毛片女人18水好多| 久久精品综合一区二区三区| netflix在线观看网站| or卡值多少钱| 国产免费男女视频| 亚洲国产欧美一区二区综合| 亚洲电影在线观看av| 午夜福利视频1000在线观看| 老汉色∧v一级毛片| 别揉我奶头~嗯~啊~动态视频| 每晚都被弄得嗷嗷叫到高潮| 精品一区二区三区四区五区乱码| 色在线成人网| 最新在线观看一区二区三区| 18禁黄网站禁片免费观看直播| 99久久综合精品五月天人人| 国产免费av片在线观看野外av| 欧洲精品卡2卡3卡4卡5卡区| 岛国视频午夜一区免费看| 国产人伦9x9x在线观看| 色噜噜av男人的天堂激情| 亚洲成a人片在线一区二区| 国产三级在线视频| 美女 人体艺术 gogo| 女同久久另类99精品国产91| 国产精华一区二区三区| 男人舔女人的私密视频| 国产成人欧美在线观看| 久久精品国产亚洲av香蕉五月| 亚洲第一电影网av| 亚洲国产欧洲综合997久久,| 人人妻人人看人人澡| 精品国产乱码久久久久久男人| 成人高潮视频无遮挡免费网站| 啪啪无遮挡十八禁网站| 精品国产超薄肉色丝袜足j| 一级作爱视频免费观看| 国产野战对白在线观看| 欧美另类亚洲清纯唯美| 在线免费观看的www视频| 精品国内亚洲2022精品成人| 天天添夜夜摸| 白带黄色成豆腐渣| 欧美黑人巨大hd| 一区二区三区激情视频| 狂野欧美白嫩少妇大欣赏| 99国产极品粉嫩在线观看| 首页视频小说图片口味搜索| 男女下面进入的视频免费午夜| 级片在线观看| 久久久久久国产a免费观看| 亚洲国产欧美网| 欧美日韩乱码在线| 亚洲精品一卡2卡三卡4卡5卡| 日本一二三区视频观看| 精品电影一区二区在线| 狂野欧美激情性xxxx| 精品国产三级普通话版| 国内精品久久久久精免费| 久久久久久大精品| 色老头精品视频在线观看| 色精品久久人妻99蜜桃| 12—13女人毛片做爰片一| 国产日本99.免费观看| 亚洲国产精品999在线| 欧美国产日韩亚洲一区| 99在线人妻在线中文字幕| 在线观看舔阴道视频| 欧美zozozo另类| 国内久久婷婷六月综合欲色啪| 九九久久精品国产亚洲av麻豆 | 精品久久久久久,| 国产精品久久电影中文字幕| 亚洲熟妇中文字幕五十中出| 香蕉久久夜色| 欧美成人一区二区免费高清观看 | 亚洲成a人片在线一区二区| 色噜噜av男人的天堂激情| xxx96com| 最近视频中文字幕2019在线8| 免费看光身美女| 丝袜人妻中文字幕| 好男人电影高清在线观看| 老司机午夜福利在线观看视频| 国产精品av视频在线免费观看| xxx96com| av在线天堂中文字幕| 欧美日韩精品网址| 精品一区二区三区视频在线 | 日韩欧美免费精品| 一级黄色大片毛片| 热99re8久久精品国产| 国产精品久久久av美女十八| 真实男女啪啪啪动态图| 嫩草影院精品99| 国产乱人伦免费视频| 欧美色视频一区免费| 亚洲乱码一区二区免费版| 国内久久婷婷六月综合欲色啪| 变态另类丝袜制服| 少妇熟女aⅴ在线视频| 黄片大片在线免费观看| 伊人久久大香线蕉亚洲五| 成人三级做爰电影| 深夜精品福利| 91在线观看av| 亚洲精品久久国产高清桃花| 18禁观看日本| 成熟少妇高潮喷水视频| 亚洲av成人av| 亚洲最大成人中文| 久久久色成人| 美女免费视频网站| 国产在线精品亚洲第一网站| 国产精品女同一区二区软件 | 精品一区二区三区四区五区乱码| 国产高清激情床上av| 久久久精品大字幕| 一夜夜www| 精品国产超薄肉色丝袜足j| 国产视频内射| 国产精品,欧美在线| 看黄色毛片网站| 老鸭窝网址在线观看| 老熟妇仑乱视频hdxx| 色综合欧美亚洲国产小说| 黄片小视频在线播放| 亚洲国产日韩欧美精品在线观看 | 黄色丝袜av网址大全| 成人特级av手机在线观看| 免费无遮挡裸体视频| 欧美色视频一区免费| 少妇人妻一区二区三区视频| 精品久久久久久久久久免费视频| 观看美女的网站| 精品一区二区三区视频在线观看免费| 欧美3d第一页| 在线免费观看的www视频| 国产亚洲精品一区二区www| 亚洲精品久久国产高清桃花| 亚洲真实伦在线观看| 国内久久婷婷六月综合欲色啪| 久久久国产精品麻豆| 黄色日韩在线| 免费在线观看成人毛片| 日韩欧美一区二区三区在线观看| 99国产精品99久久久久| 国产精品久久视频播放| 久久久久亚洲av毛片大全| 久久久久精品国产欧美久久久| 亚洲无线在线观看| 亚洲专区国产一区二区| 性色avwww在线观看| 久久久久亚洲av毛片大全| 伦理电影免费视频| 最近视频中文字幕2019在线8| 在线观看免费视频日本深夜| 日本黄色片子视频| 99精品欧美一区二区三区四区| 国产伦一二天堂av在线观看| 亚洲成av人片免费观看| 国产成人系列免费观看| 99热精品在线国产| 亚洲黑人精品在线| 男人舔奶头视频| 精品久久久久久久久久久久久| 舔av片在线| 成人特级黄色片久久久久久久| 性欧美人与动物交配| 波多野结衣高清作品| 亚洲va日本ⅴa欧美va伊人久久| 亚洲色图 男人天堂 中文字幕| 欧美在线黄色| 精品久久久久久久人妻蜜臀av| 久久精品91无色码中文字幕| 精品电影一区二区在线| 日本在线视频免费播放| 极品教师在线免费播放| 久久久久免费精品人妻一区二区| 久久久久九九精品影院| 美女大奶头视频| 国内精品久久久久精免费| av视频在线观看入口| 国产乱人视频| 午夜久久久久精精品| 观看美女的网站| 99久久久亚洲精品蜜臀av| 一卡2卡三卡四卡精品乱码亚洲| 欧美一区二区国产精品久久精品| 在线观看免费视频日本深夜| 国产精品亚洲av一区麻豆| 亚洲国产欧洲综合997久久,| 欧美不卡视频在线免费观看| 精品一区二区三区四区五区乱码| 亚洲,欧美精品.| 成人高潮视频无遮挡免费网站| 亚洲成人久久爱视频| 亚洲国产看品久久| 国产精品野战在线观看| 久久久久国内视频| 亚洲国产精品sss在线观看| 可以在线观看毛片的网站| 久久精品91无色码中文字幕| 亚洲自拍偷在线| 啦啦啦韩国在线观看视频| 亚洲国产欧美一区二区综合| 国产三级在线视频| 免费在线观看日本一区| 国产精品久久视频播放| 亚洲成人精品中文字幕电影| 成人欧美大片| 麻豆av在线久日| 欧美av亚洲av综合av国产av| 日韩大尺度精品在线看网址| 黄频高清免费视频| 日韩人妻高清精品专区| 最新美女视频免费是黄的| 两个人视频免费观看高清| 久久欧美精品欧美久久欧美| 亚洲人成网站高清观看| 国产精品久久久久久人妻精品电影| 精品人妻1区二区| 亚洲欧美日韩高清在线视频| 最近最新免费中文字幕在线| 欧美日本亚洲视频在线播放| 成人一区二区视频在线观看| 欧美xxxx黑人xx丫x性爽| 久久久国产欧美日韩av| 免费一级毛片在线播放高清视频| 久久久久九九精品影院|