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

    一類具有標(biāo)準(zhǔn)發(fā)生率和雙垂直傳播的媒介傳染病模型分析

    2021-07-23 09:40:50劉晨竇霽虹李玉峰趙婷婷
    關(guān)鍵詞:感者染病平衡點(diǎn)

    劉晨, 竇霽虹, 李玉峰, 趙婷婷

    (西北大學(xué)數(shù)學(xué)學(xué)院, 陜西 西安 710127)

    1 引言

    媒介傳染病是通過一種被稱為生物媒介載體傳播的傳染性疾病, 通常由病毒、細(xì)菌、原生動物或立克次體等引起的, 特別是蚊媒疾病, 例如瘧疾、登革熱、西尼羅河熱等[1], 給人們的生活帶來巨大的影響, 時刻威脅著人們的生命健康. 隨著科學(xué)技術(shù)的迅速發(fā)展, 醫(yī)學(xué)界不斷嘗試去解決長久以來困擾人們的疾病問題. 然而, 理論是實(shí)踐的基礎(chǔ), 大量的傳染病動力學(xué)模型推動著對傳染病問題的研究. 人們更多地關(guān)注決定疾病爆發(fā)后是持續(xù)流行還是逐漸消亡的閾值, 即基本再生數(shù)[2], 從而能更好地控制疾病的傳播.

    根據(jù)文獻(xiàn)[3-4] 可知染病媒介的繁殖后代中有部分?jǐn)y帶病毒, 能夠進(jìn)行疾病的傳播;在2014 年關(guān)于新生兒登革熱的報告中, 表明新生兒的病因正是由于母嬰垂直傳播[5].因此, 在實(shí)際的病例中確實(shí)存在人群和媒介均具有垂直傳播這種雙垂直傳播傳染病的傳播機(jī)制, 但現(xiàn)今對雙方同時考慮垂直傳播特征的疾病研究還為數(shù)不多.

    因此, 本文建立了一類具有標(biāo)準(zhǔn)發(fā)生率和雙垂直傳播的媒介傳染病模型, 并對其動力學(xué)形態(tài)進(jìn)行理論分析證明與數(shù)值模擬, 進(jìn)一步完善對媒介傳染病模型的研究.

    2 模型建立

    假設(shè)在傳染病流行期間, 該地區(qū)人群的遷移忽略不計, 同時不考慮因病死亡, 且人群的出生率和死亡率相等用μH表示, 則總?cè)巳壕S持于一個固定的常數(shù)N.SH(t),IH(t) 及RH(t) 分別表示易感者, 帶病毒人群(簡稱患者) 及康復(fù)者在t時刻相應(yīng)的數(shù)量. 對于媒介(如蚊子) 而言, 繁殖周期短, 產(chǎn)量大, 每日的出生量和死亡量都十分巨大,在無人為干擾的情況下, 媒介總數(shù)的變化量相對于該種群總數(shù)可忽略不計, 因此假設(shè)媒介總數(shù)也維持于一個固定常數(shù)M, 媒介的出生率和死亡率均為μM.SM(t),IM(t) 分別表示易感媒介和染病媒介在t時刻相應(yīng)的數(shù)量.

    根據(jù)媒介習(xí)性的不同, 將易感媒介每日接觸人的頻率記為a1, 染病媒介每日接觸人的頻率相對更高記為a2, 同時感染效率也有所不同[6-7], 將易感媒介接觸一次患者獲病的概率記為b1, 易感者接觸一次染病媒介病的概率相對更高記為b2, 疾病的發(fā)生率均為標(biāo)準(zhǔn)發(fā)生率.r表示患者的康復(fù)率, 患者一旦康復(fù)可對該疾病終身免疫, 人群和媒介各自的垂直傳染率分別為qH及qM. 可得此類媒介傳染病傳播的流程圖如圖1.

    圖1 媒介傳染病的傳播流程圖

    根據(jù)建模思想, 本文建立如下動力學(xué)模型:

    根據(jù)Rh(t) = 1?Sh(t)?Ih(t), Sm(t) = 1?Im(t), 對模型(1) 進(jìn)行降維, 因此只需考慮如下系統(tǒng):

    考慮到實(shí)際背景, 模型的解均為非負(fù)數(shù), 本文僅在可行域Ω 中進(jìn)行討論,

    3 基本再生數(shù)的計算

    令模型(2) 右邊各等式等于0, 易得該模型的無病平衡點(diǎn)E0(1,0,0). 本文采用Driessche 和Watmough 的第二代生成矩陣法[8], 給出基本再生數(shù)的表示.

    令X=(Ih,Im,Sm)T,β1=a1b1,β2=a2b2, 則模型(2) 可等價為

    其中,

    φ,ψ在無病平衡點(diǎn)E0處的雅克比矩陣分別如下:

    當(dāng)基本再生數(shù)R0<1 時, 意味著平均每位患者在其患病期內(nèi)導(dǎo)致新增患病人數(shù)不足一人, 疾病逐漸消亡; 反之, 當(dāng)基本再生數(shù)R0>1 時, 意味著疾病持續(xù)蔓延.

    4 平衡點(diǎn)的存在性

    定理4.1若R0<1, 模型(2) 只存在無病平衡點(diǎn)E0(1,0,0); 若R0>1, 模型(2)存在唯一的地方病平衡點(diǎn)E?(S?h,I?h,I?m).

    證明令模型(2) 右邊各等式等于0, 計算得

    當(dāng)R0<1 時,Ih=0, 模型(2) 只存在無病平衡點(diǎn)E0.

    5 平衡點(diǎn)的穩(wěn)定性

    定理5.1當(dāng)R0< 1, 模型(2) 的無病平衡點(diǎn)E0(1,0,0) 是局部漸近穩(wěn)定的;當(dāng)R0>1 時, 則不穩(wěn)定.

    證明模型(2) 在無病平衡點(diǎn)E0(1,0,0) 處的雅克比矩陣為

    對應(yīng)的特征方程為

    解得λ1=?μH<0.另外兩個特征根由下面的方程決定

    整理得

    其中,c1=μM(1?qM)+μH(1?qH)+r, c2=(1?R20)μM(1?qM)[μH(1?qH)+r].

    由c1> 0, 且當(dāng)R0< 1 時,c2> 0, 根據(jù)Hurwitz 判據(jù)可知, 方程(7) 所有特征根均具有負(fù)實(shí)部. 因此, 無病平衡點(diǎn)E0是局部漸近穩(wěn)定的; 當(dāng)R0> 1 時,c2< 0, 則方程(9) 有一正根, 故無病平衡點(diǎn)E0不穩(wěn)定.

    定理5.2當(dāng)R0<1, 模型(2) 的無病平衡點(diǎn)E0(1,0,0) 是全局漸近穩(wěn)定.

    證明構(gòu)造正定的Lyapunov 函數(shù)

    其沿模型(2) 軌線的導(dǎo)數(shù)為

    當(dāng)R0<1 時, ˙V ≤0, 當(dāng)且僅當(dāng)Ih=0 時, ˙V=0. 此時Sh=1,Im=0, 故無病平衡點(diǎn)為唯一的不動點(diǎn),{E0}是{Sh,Ih,Im ∈Ω : ˙V= 0}的最大不變集. 由LaSalle 不變集原理[9]可得, 當(dāng)R0<1 時, 模型(2) 的無病平衡點(diǎn)E0是全局漸近穩(wěn)定的.

    定理5.3在R0>1 的基礎(chǔ)上,當(dāng)滿足條件(H1)R20[μH(1?I?m)?(μH+r)I?h)]<μH時, 模型(2) 的地方病平衡點(diǎn)E?(S?h,I?h,I?m) 是局部漸近穩(wěn)定的.

    證明當(dāng)R0>1 時, 模型(2) 在地方病平衡點(diǎn)E?處的雅克比矩陣經(jīng)矩陣變換為

    其中

    經(jīng)整理

    由條件(H1), 得

    則B<0. 因此, 地方病平衡點(diǎn)E?(S?h,I?h,I?m) 是局部漸近穩(wěn)定的.

    定理5.4在R0> 1 的基礎(chǔ)上, 當(dāng)滿足條件(H2)μH+r ?2qHμH> 0 時, 模型(2) 的地方病平衡點(diǎn)E?(S?h,I?h,I?m) 是全局漸近穩(wěn)定的.

    證明根據(jù)文獻(xiàn)[10] 及定理5.2 知, 當(dāng)R0>1 時,E0不穩(wěn)定, 則在Ω 內(nèi)存在一個緊吸引子集Γ, 且模型(2) 存在唯一的地方病平衡點(diǎn)E?, 故滿足文獻(xiàn)[11] 中定理3.3.7的條件, 下面證明q<0.

    模型(2) 的雅克比矩陣的第二加性復(fù)合矩陣為

    其中,

    其中,

    令(u,v,l)∈R3, 其范數(shù)‖·‖定義為‖(u,v,l)‖= max{|u|,|v|+|l|}, 相應(yīng)于范數(shù)‖·‖的Lozinskiˇi 測度是μ(B),μ(B)≤sup{g1,g2}, 其中

    將B22的每一列非對角元素取絕對值后加到相應(yīng)列的對角元素上, 取其兩個對角元素的最大值, 則

    由條件(H2) 知,η=min{μH,μH+r ?2qHμH}>0. 利用

    則對所有滿足初值的X0=(Sh(0),Ih(0),Im(0))∈Γ, 當(dāng)t>t?時, 有

    又由Ω?R3是單連通, 故根據(jù)文獻(xiàn)[11] 中的定理3.3.7 證畢.

    6 數(shù)值模擬

    為了對以上理論結(jié)果進(jìn)行驗(yàn)證, 本文參考文獻(xiàn)[6] 對參數(shù)進(jìn)行適當(dāng)?shù)剡x取, 借助Matlab 軟件進(jìn)行數(shù)值模擬, 針對不同的基本再生數(shù)計算結(jié)果, 描繪出疾病的發(fā)展趨勢.

    (1) 選取參數(shù)K= 10,a1= 0.12,a2= 0.28,b1= 0.5,b2= 0.8,qH= 0.3,qM= 0.36,μH= 0.2,μM= 0.6,r= 0.25. 此時計算得R0= 0.9473< 1. 圖2 為易感者比例Sh(t)、患者比例I(t) 及染病媒介比例Im(t) 隨時間t的變化情況.

    圖2 無病平衡點(diǎn)的穩(wěn)定性(R0 =0.9473)

    (2)選取參數(shù)K=10,a1=0.12,a2=0.2,b1=0.5,b2=0.8,qH=0.3,qM=0.36,μH=0.2,μM=0.6,r=0.4. 此時計算得R0=0.6211<1. 圖3 為易感者比例Sh(t)、患者比例I(t) 及染病媒介比例Im(t) 隨時間t的變化情況.

    圖3 無病平衡點(diǎn)的穩(wěn)定性(R0 =0.6211)

    通過圖2 和圖3 可以直觀地看出, 當(dāng)R0<1 時, 易感者比例逐漸趨于1, 患者比例和染病媒介比例逐漸趨于0, 無病平衡點(diǎn)在可行域內(nèi)是全局漸近穩(wěn)定的. 另外, 通過比較兩圖, 可以發(fā)現(xiàn)當(dāng)減小人群與媒介之間的接觸率、增大疾病的康復(fù)率, 會引起基本再生數(shù)減小, 從而導(dǎo)致疾病更快的走向消亡.

    (3) 選取參數(shù)K= 10,a1= 0.25,a2= 0.37,b1= 0.5,b2= 0.8,qH= 0.3,qM= 0.36,μH= 0.2,μM= 0.6,r= 0.08. 此時計算得R0= 2.0928> 1. 圖4 為易感者比例Sh(t)、患者比例I(t) 及染病媒介比例Im(t) 隨時間t的變化情況.

    圖4 地方病平衡點(diǎn)的穩(wěn)定性(R0 =2.0928)

    (4) 選取參數(shù)K= 10,a1= 0.12,a2= 0.28,b1= 0.5,b2= 0.8,qH= 0.3,qM=0.36,μH=0.2,μM=0.6,r=0.1. 此時計算得R0=1.2076>1. 圖5 為易感者比例Sh(t)、患者比例I(t) 及染病媒介比例Im(t) 隨時間t的變化情況.

    圖5 地方病平衡點(diǎn)的穩(wěn)定性(R0 =1.2076)

    通過圖4 和圖5 可以直觀地看出, 當(dāng)R0>1 時, 隨著時間的推移, 易感者比例、患者比例和染病媒介的比例均趨于各自的一個穩(wěn)定值, 可見地方病平衡點(diǎn)在可行域內(nèi)是全局漸近穩(wěn)定的. 另外, 通過比較兩圖, 可以發(fā)現(xiàn)若減小人群與媒介之間的接觸率、增大疾病的康復(fù)率, 從而引起基本再生數(shù)減小, 則會導(dǎo)致當(dāng)疾病發(fā)展到地方病平衡點(diǎn)時, 相應(yīng)的患者比例會穩(wěn)定在更小的數(shù)值.

    進(jìn)一步根據(jù)基本再生數(shù)的表達(dá)式

    可知,R0與人群的垂直傳染率qH及媒介的垂直傳染率qM正相關(guān). 因此, 相比宿主和媒介中一方具有垂直傳播, 或雙方均不考慮垂直傳播的情況而言, 當(dāng)加入雙垂直傳播特征時, 會引起基本再生數(shù)增大, 從而R0<1 時, 疾病消亡的時間會相對延長;R0>1 時,疾病到達(dá)地方病平衡點(diǎn)相應(yīng)的患者比例有所提高.

    7 結(jié)論

    本文考慮到在媒介傳染病中存在傳播媒介和宿主均具有母嬰垂直傳播的實(shí)際情況,建立了一類具有標(biāo)準(zhǔn)發(fā)生率和雙垂直傳播的媒介傳染病模型, 并對模型進(jìn)行了定性分析和定量模擬, 從而達(dá)到了對媒介傳染病進(jìn)一步完善研究的目的. 最后, 通過計算機(jī)不僅實(shí)現(xiàn)了對模型的數(shù)值模擬, 驗(yàn)證了理論結(jié)果的有效性, 而且通過比較不同基本再生數(shù)下的疾病發(fā)展趨勢, 直觀地說明了采取減小人群與媒介之間接觸率的相應(yīng)措施, 或者采取有效的治療方法提高康復(fù)率, 能夠使得基本再生數(shù)降低, 從而可以達(dá)到控制疾病蔓延的效果, 此為實(shí)際應(yīng)用提供了理論依據(jù). 另一方面, 說明了相對于不考慮雙垂直傳播因素,增加該因素之后, 會引起基本再生數(shù)的增大, 從而導(dǎo)致疾病消亡的時間延長, 地方病平衡點(diǎn)對應(yīng)的患者比例提高.

    猜你喜歡
    感者染病平衡點(diǎn)
    考慮媒體影響的一類時滯傳染病模型的分岔周期解
    一類具有年齡結(jié)構(gòu)和接種干預(yù)的手足口病模型動力學(xué)分析
    偶感
    分析采取措施對性病傳播動態(tài)的影響
    探尋中國蘋果產(chǎn)業(yè)的產(chǎn)銷平衡點(diǎn)
    煙臺果樹(2019年1期)2019-01-28 09:34:58
    均勻網(wǎng)絡(luò)上SIR模型三種不同逼近方法比較
    電視庭審報道,如何找到媒體監(jiān)督與司法公正的平衡點(diǎn)
    傳媒評論(2018年7期)2018-09-18 03:45:52
    愛 情
    詩選刊(2016年9期)2016-11-26 13:47:43
    在給專車服務(wù)正名之前最好找到Uber和出租車的平衡點(diǎn)
    IT時代周刊(2015年7期)2015-11-11 05:49:56
    食餌染病生態(tài)流行病系統(tǒng)的穩(wěn)定性與最優(yōu)收獲
    97超碰精品成人国产| 国产精品一区二区三区四区免费观看| 美女国产视频在线观看| 免费看不卡的av| 一级黄片播放器| 国产亚洲91精品色在线| 国产白丝娇喘喷水9色精品| 一级av片app| 超碰av人人做人人爽久久| 久久久成人免费电影| 熟女电影av网| 日韩在线高清观看一区二区三区| 美女高潮的动态| 日韩成人av中文字幕在线观看| 2021天堂中文幕一二区在线观| 久久久久久久久久成人| 国产黄色视频一区二区在线观看| av在线天堂中文字幕| 26uuu在线亚洲综合色| 午夜福利在线观看免费完整高清在| 日日啪夜夜爽| 免费在线观看成人毛片| 精品少妇久久久久久888优播| 国产成人免费观看mmmm| 国产一区二区在线观看日韩| 不卡视频在线观看欧美| 久热久热在线精品观看| 国产成人a∨麻豆精品| 美女国产视频在线观看| 久久女婷五月综合色啪小说 | 成人特级av手机在线观看| 国产精品一区二区性色av| 成人一区二区视频在线观看| 秋霞在线观看毛片| 亚洲精品国产成人久久av| 黄片wwwwww| 亚洲性久久影院| 极品少妇高潮喷水抽搐| 18禁在线无遮挡免费观看视频| 亚洲欧洲国产日韩| 成人亚洲精品一区在线观看 | 欧美bdsm另类| 亚洲成色77777| 亚洲国产成人一精品久久久| 天天躁日日操中文字幕| 国产精品秋霞免费鲁丝片| 国产高清不卡午夜福利| 成年女人在线观看亚洲视频 | 18禁在线无遮挡免费观看视频| 一本色道久久久久久精品综合| 纵有疾风起免费观看全集完整版| 国产免费又黄又爽又色| 久久久久精品性色| 色哟哟·www| 国产精品成人在线| 日韩国内少妇激情av| 韩国av在线不卡| 久久午夜福利片| 欧美精品国产亚洲| 国产真实伦视频高清在线观看| 高清午夜精品一区二区三区| 小蜜桃在线观看免费完整版高清| 欧美 日韩 精品 国产| 卡戴珊不雅视频在线播放| 日韩三级伦理在线观看| 91aial.com中文字幕在线观看| 网址你懂的国产日韩在线| 欧美97在线视频| a级毛色黄片| 国产亚洲av嫩草精品影院| 舔av片在线| 免费不卡的大黄色大毛片视频在线观看| 免费高清在线观看视频在线观看| 精品国产一区二区三区久久久樱花 | 女的被弄到高潮叫床怎么办| 下体分泌物呈黄色| 91在线精品国自产拍蜜月| 国产av不卡久久| 亚洲伊人久久精品综合| 亚洲av国产av综合av卡| av福利片在线观看| 国产毛片在线视频| av专区在线播放| 久久久a久久爽久久v久久| 少妇 在线观看| 97热精品久久久久久| 综合色av麻豆| 亚洲精品色激情综合| 蜜臀久久99精品久久宅男| 国产有黄有色有爽视频| www.av在线官网国产| 欧美日本视频| 国产淫语在线视频| av在线app专区| 国产精品一二三区在线看| 51国产日韩欧美| 亚洲高清免费不卡视频| 偷拍熟女少妇极品色| 亚洲av福利一区| 亚洲av男天堂| 中文字幕av成人在线电影| 亚洲国产精品成人综合色| 亚洲自拍偷在线| 日本av手机在线免费观看| 日本-黄色视频高清免费观看| 最近中文字幕2019免费版| 亚洲不卡免费看| 男人爽女人下面视频在线观看| 国产亚洲av嫩草精品影院| 一级毛片我不卡| 在线观看免费高清a一片| 日韩成人av中文字幕在线观看| 国产精品精品国产色婷婷| 舔av片在线| 成人漫画全彩无遮挡| 国产视频首页在线观看| 欧美国产精品一级二级三级 | 久久久久久久大尺度免费视频| 极品教师在线视频| 亚洲精品中文字幕在线视频 | 亚洲三级黄色毛片| 中文字幕久久专区| 成人欧美大片| 成人一区二区视频在线观看| 久久久久久久亚洲中文字幕| 日本三级黄在线观看| 内地一区二区视频在线| 晚上一个人看的免费电影| 美女视频免费永久观看网站| 91精品国产九色| 日本熟妇午夜| 欧美97在线视频| 免费人成在线观看视频色| 久热这里只有精品99| 国产在线男女| 中国美白少妇内射xxxbb| 国产又色又爽无遮挡免| 亚洲精华国产精华液的使用体验| 国产欧美日韩精品一区二区| 精品久久久噜噜| 久久久久久久大尺度免费视频| 青春草国产在线视频| 黄色视频在线播放观看不卡| av女优亚洲男人天堂| 插阴视频在线观看视频| 欧美xxxx黑人xx丫x性爽| 欧美日韩精品成人综合77777| 日韩欧美一区视频在线观看 | 一本久久精品| 精品国产三级普通话版| 伊人久久精品亚洲午夜| 麻豆精品久久久久久蜜桃| 激情五月婷婷亚洲| 精品人妻偷拍中文字幕| av黄色大香蕉| 欧美成人a在线观看| av又黄又爽大尺度在线免费看| 一个人看的www免费观看视频| 亚洲国产成人一精品久久久| videossex国产| 丰满乱子伦码专区| 一级二级三级毛片免费看| 国产亚洲最大av| 国产视频首页在线观看| 搡老乐熟女国产| 午夜免费鲁丝| 国国产精品蜜臀av免费| 97在线视频观看| 超碰av人人做人人爽久久| 亚洲aⅴ乱码一区二区在线播放| 久久6这里有精品| 嫩草影院入口| 春色校园在线视频观看| 国产亚洲av嫩草精品影院| 水蜜桃什么品种好| 色视频www国产| 蜜桃亚洲精品一区二区三区| 熟女人妻精品中文字幕| 久久久久网色| 99热这里只有是精品在线观看| 国产精品国产三级国产专区5o| 色视频www国产| 亚洲av不卡在线观看| 丝瓜视频免费看黄片| 精品人妻一区二区三区麻豆| 又爽又黄a免费视频| 男人添女人高潮全过程视频| 免费少妇av软件| 熟女av电影| 欧美+日韩+精品| www.av在线官网国产| 你懂的网址亚洲精品在线观看| 交换朋友夫妻互换小说| xxx大片免费视频| 亚洲精品国产色婷婷电影| 一边亲一边摸免费视频| 亚洲怡红院男人天堂| 草草在线视频免费看| 亚洲伊人久久精品综合| 国产成年人精品一区二区| 国产一区二区亚洲精品在线观看| 91久久精品国产一区二区三区| 国产精品三级大全| 国产精品爽爽va在线观看网站| 校园人妻丝袜中文字幕| 久久人人爽av亚洲精品天堂 | 日韩伦理黄色片| 各种免费的搞黄视频| 国产高清国产精品国产三级 | 亚洲欧洲日产国产| 综合色av麻豆| 久久人人爽av亚洲精品天堂 | 国产亚洲91精品色在线| 七月丁香在线播放| 国产欧美另类精品又又久久亚洲欧美| 国产黄a三级三级三级人| 丰满乱子伦码专区| 激情 狠狠 欧美| 国产黄片视频在线免费观看| 91精品伊人久久大香线蕉| 男插女下体视频免费在线播放| 国产精品一区www在线观看| 好男人在线观看高清免费视频| 男人和女人高潮做爰伦理| 性色avwww在线观看| 春色校园在线视频观看| 最后的刺客免费高清国语| 五月天丁香电影| 麻豆乱淫一区二区| 极品教师在线视频| 国产日韩欧美亚洲二区| 亚洲成人一二三区av| 极品教师在线视频| 免费黄频网站在线观看国产| 欧美潮喷喷水| 亚洲国产精品999| 一个人看的www免费观看视频| 国产精品国产三级专区第一集| 久久久久久伊人网av| av在线天堂中文字幕| 最近手机中文字幕大全| 国产亚洲一区二区精品| 亚洲av成人精品一区久久| 伦精品一区二区三区| 制服丝袜香蕉在线| 亚洲精品成人久久久久久| 亚洲伊人久久精品综合| 99九九线精品视频在线观看视频| 精品久久久精品久久久| 少妇熟女欧美另类| 国产在线一区二区三区精| 80岁老熟妇乱子伦牲交| 亚洲经典国产精华液单| 亚洲av国产av综合av卡| 国产精品国产三级国产av玫瑰| 成人特级av手机在线观看| 国语对白做爰xxxⅹ性视频网站| 又爽又黄无遮挡网站| 精品一区二区三卡| 亚洲色图av天堂| av在线观看视频网站免费| 亚洲国产欧美在线一区| 亚洲aⅴ乱码一区二区在线播放| 最近的中文字幕免费完整| 亚洲人成网站在线播| 看黄色毛片网站| 国产精品久久久久久久电影| 国产 一区精品| 九色成人免费人妻av| 成人无遮挡网站| 美女高潮的动态| 新久久久久国产一级毛片| 久久久欧美国产精品| 亚洲av日韩在线播放| 久久久久精品久久久久真实原创| 国产男人的电影天堂91| 亚洲最大成人av| 亚洲高清免费不卡视频| 国产成人aa在线观看| 色视频www国产| 久久精品熟女亚洲av麻豆精品| 黄色配什么色好看| 国产av码专区亚洲av| 天天躁夜夜躁狠狠久久av| av免费观看日本| 国产国拍精品亚洲av在线观看| 综合色丁香网| 最近手机中文字幕大全| 啦啦啦啦在线视频资源| 精品一区二区三卡| 一个人看的www免费观看视频| av播播在线观看一区| av又黄又爽大尺度在线免费看| 欧美激情在线99| 日韩一区二区三区影片| 欧美一级a爱片免费观看看| 国产成人精品福利久久| 男人舔奶头视频| 五月玫瑰六月丁香| 久久久久精品性色| 极品少妇高潮喷水抽搐| 熟女av电影| 寂寞人妻少妇视频99o| 国产精品秋霞免费鲁丝片| 国产综合精华液| 亚洲成人久久爱视频| 国产 一区 欧美 日韩| 白带黄色成豆腐渣| 亚洲精品视频女| 国产精品国产三级国产av玫瑰| 一区二区三区四区激情视频| 99热6这里只有精品| 一级爰片在线观看| 国产伦在线观看视频一区| 欧美精品一区二区大全| 国产精品秋霞免费鲁丝片| 最近中文字幕高清免费大全6| 97精品久久久久久久久久精品| 蜜桃亚洲精品一区二区三区| 免费av观看视频| 身体一侧抽搐| 七月丁香在线播放| av卡一久久| 国产亚洲一区二区精品| 18+在线观看网站| 久久久久久久久久人人人人人人| 日韩精品有码人妻一区| 国产黄片视频在线免费观看| 少妇 在线观看| 久久久精品免费免费高清| av.在线天堂| 一级毛片 在线播放| 国精品久久久久久国模美| 成年免费大片在线观看| 国产又色又爽无遮挡免| 久久久亚洲精品成人影院| 国产精品人妻久久久影院| 超碰av人人做人人爽久久| 天堂俺去俺来也www色官网| 伊人久久精品亚洲午夜| 在线天堂最新版资源| 亚洲伊人久久精品综合| 国精品久久久久久国模美| 在线免费观看不下载黄p国产| 免费大片18禁| 国产欧美亚洲国产| 下体分泌物呈黄色| 久久久久久伊人网av| 性色av一级| 啦啦啦啦在线视频资源| 国产亚洲一区二区精品| 午夜精品一区二区三区免费看| 夜夜看夜夜爽夜夜摸| 色吧在线观看| 插阴视频在线观看视频| 亚洲国产精品国产精品| 免费看光身美女| 美女主播在线视频| 免费观看无遮挡的男女| 亚洲不卡免费看| 大香蕉久久网| 亚洲国产av新网站| av国产久精品久网站免费入址| 99热这里只有精品一区| 精品酒店卫生间| 亚洲国产精品专区欧美| av国产久精品久网站免费入址| 丝瓜视频免费看黄片| 女人久久www免费人成看片| 街头女战士在线观看网站| 美女xxoo啪啪120秒动态图| 美女高潮的动态| 国产成人一区二区在线| 自拍偷自拍亚洲精品老妇| 久久综合国产亚洲精品| 欧美3d第一页| 大香蕉久久网| 少妇丰满av| 性色avwww在线观看| 99热这里只有是精品50| 两个人的视频大全免费| 久久精品综合一区二区三区| 日韩一区二区三区影片| 尾随美女入室| 国产 一区精品| 一二三四中文在线观看免费高清| 伊人久久国产一区二区| 插逼视频在线观看| 亚洲一级一片aⅴ在线观看| 久久久久精品性色| 国产一区有黄有色的免费视频| 夜夜看夜夜爽夜夜摸| 国产成年人精品一区二区| 亚洲精品国产av成人精品| 国产色爽女视频免费观看| 永久网站在线| 欧美国产精品一级二级三级 | 麻豆成人av视频| 国产精品国产三级专区第一集| 狂野欧美激情性xxxx在线观看| 少妇人妻精品综合一区二区| 久久亚洲国产成人精品v| 欧美日韩在线观看h| 欧美激情在线99| 亚洲精品乱码久久久v下载方式| 一个人看的www免费观看视频| 久久99热这里只有精品18| 成人无遮挡网站| 青春草视频在线免费观看| 午夜日本视频在线| 偷拍熟女少妇极品色| 欧美 日韩 精品 国产| 国产精品一区二区在线观看99| 一区二区av电影网| 久久午夜福利片| 国产色婷婷99| 亚洲av中文av极速乱| 校园人妻丝袜中文字幕| 美女视频免费永久观看网站| 国产成人精品福利久久| 国产精品av视频在线免费观看| 一二三四中文在线观看免费高清| 欧美日韩亚洲高清精品| 国产一区二区三区综合在线观看 | 99热6这里只有精品| 欧美人与善性xxx| 日韩制服骚丝袜av| 国产精品熟女久久久久浪| 97超碰精品成人国产| 亚洲欧美日韩无卡精品| 久久97久久精品| 涩涩av久久男人的天堂| 成人毛片60女人毛片免费| 久久久精品欧美日韩精品| 国产亚洲一区二区精品| 在线a可以看的网站| 黑人高潮一二区| 久久精品久久久久久久性| 亚洲经典国产精华液单| 国产成人福利小说| 特大巨黑吊av在线直播| 国产精品一区二区在线观看99| 又粗又硬又长又爽又黄的视频| 美女高潮的动态| 欧美日韩综合久久久久久| 国产国拍精品亚洲av在线观看| 久久精品久久久久久久性| 久久女婷五月综合色啪小说 | 大片免费播放器 马上看| 国产日韩欧美在线精品| 亚洲精品亚洲一区二区| 少妇猛男粗大的猛烈进出视频 | 日本wwww免费看| 大香蕉97超碰在线| 国产在线男女| 菩萨蛮人人尽说江南好唐韦庄| 国产美女午夜福利| 大香蕉久久网| 国产黄a三级三级三级人| 欧美日韩亚洲高清精品| 国产成人精品久久久久久| 午夜亚洲福利在线播放| 午夜精品一区二区三区免费看| 只有这里有精品99| 91在线精品国自产拍蜜月| 国模一区二区三区四区视频| 午夜视频国产福利| 亚洲熟女精品中文字幕| 好男人在线观看高清免费视频| 国产伦理片在线播放av一区| 国产成人91sexporn| 欧美xxxx性猛交bbbb| 成人免费观看视频高清| 啦啦啦在线观看免费高清www| 超碰97精品在线观看| 制服丝袜香蕉在线| 久久99热这里只频精品6学生| 一级爰片在线观看| av免费观看日本| 国产精品嫩草影院av在线观看| 一个人看视频在线观看www免费| 黑人高潮一二区| 国产精品伦人一区二区| 国产乱人视频| 97热精品久久久久久| 尾随美女入室| 一本色道久久久久久精品综合| 天堂俺去俺来也www色官网| 欧美xxⅹ黑人| 在线观看美女被高潮喷水网站| 国产精品麻豆人妻色哟哟久久| 亚洲精品乱码久久久久久按摩| 大片免费播放器 马上看| 日韩在线高清观看一区二区三区| 99热国产这里只有精品6| 最近2019中文字幕mv第一页| 久久女婷五月综合色啪小说 | 欧美激情在线99| 赤兔流量卡办理| 亚洲av免费在线观看| 亚洲精品国产av蜜桃| freevideosex欧美| 婷婷色综合www| www.av在线官网国产| 狂野欧美激情性xxxx在线观看| 丝袜喷水一区| eeuss影院久久| 韩国av在线不卡| 制服丝袜香蕉在线| 天天躁日日操中文字幕| 久久久欧美国产精品| 亚洲自偷自拍三级| 亚洲图色成人| 国产精品福利在线免费观看| 一本一本综合久久| 国产精品一区二区在线观看99| 深爱激情五月婷婷| 免费人成在线观看视频色| 久久久久精品久久久久真实原创| 日韩一区二区三区影片| 纵有疾风起免费观看全集完整版| 成人无遮挡网站| 亚洲最大成人中文| 成人特级av手机在线观看| 在线观看一区二区三区激情| 亚洲av成人精品一二三区| 热re99久久精品国产66热6| 纵有疾风起免费观看全集完整版| 一区二区三区免费毛片| 久久鲁丝午夜福利片| 久久综合国产亚洲精品| 大又大粗又爽又黄少妇毛片口| 亚洲熟女精品中文字幕| 亚洲欧美日韩另类电影网站 | 国产久久久一区二区三区| 日本与韩国留学比较| 欧美xxxx性猛交bbbb| 永久网站在线| 能在线免费看毛片的网站| 麻豆成人av视频| 97在线视频观看| 免费少妇av软件| 国产成人免费观看mmmm| 亚洲av不卡在线观看| 久久99精品国语久久久| 精品人妻熟女av久视频| 插逼视频在线观看| 在线观看免费高清a一片| 欧美人与善性xxx| 精品一区在线观看国产| 男人添女人高潮全过程视频| 一级毛片 在线播放| 高清毛片免费看| 丰满人妻一区二区三区视频av| 如何舔出高潮| eeuss影院久久| 美女视频免费永久观看网站| 蜜臀久久99精品久久宅男| 日韩 亚洲 欧美在线| 精品久久久久久久人妻蜜臀av| 全区人妻精品视频| 好男人视频免费观看在线| 日本熟妇午夜| 午夜亚洲福利在线播放| 亚洲欧美成人精品一区二区| 夫妻午夜视频| 午夜免费观看性视频| 自拍欧美九色日韩亚洲蝌蚪91 | 婷婷色麻豆天堂久久| 99久国产av精品国产电影| 亚洲成人av在线免费| 欧美zozozo另类| 亚洲欧美日韩东京热| 色视频在线一区二区三区| 亚洲性久久影院| 我要看日韩黄色一级片| 日韩在线高清观看一区二区三区| a级毛片免费高清观看在线播放| av专区在线播放| 丰满乱子伦码专区| 国产在线男女| 激情 狠狠 欧美| 色5月婷婷丁香| 观看美女的网站| 岛国毛片在线播放| av在线播放精品| 黄色一级大片看看| 大片免费播放器 马上看| 欧美日韩综合久久久久久| 性插视频无遮挡在线免费观看| 人人妻人人澡人人爽人人夜夜| 天天躁日日操中文字幕| 国产高清不卡午夜福利| 精品久久久精品久久久| 最近中文字幕2019免费版| 亚洲综合精品二区| 国产成人精品久久久久久| 在线 av 中文字幕| 日韩av不卡免费在线播放| 少妇熟女欧美另类| 亚洲自拍偷在线| 女的被弄到高潮叫床怎么办| 精品国产乱码久久久久久小说| 中文字幕免费在线视频6| 九九久久精品国产亚洲av麻豆| 精品亚洲乱码少妇综合久久| 免费av毛片视频| 久久久久久久亚洲中文字幕| 亚洲av福利一区| 波多野结衣巨乳人妻| 国产精品久久久久久久电影| videossex国产| 精品熟女少妇av免费看| 亚洲欧美中文字幕日韩二区| 亚洲精品亚洲一区二区|