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

    具有標準發(fā)生率和脈沖干擾的SIRS模型穩(wěn)定性分析

    2017-07-01 15:47:10馬艷麗李海霞褚正清聶東明安徽新華學院通識教育部安徽合肥230088
    長春師范大學學報 2017年6期
    關鍵詞:染病傳染病脈沖

    馬艷麗,李海霞,褚正清,聶東明(安徽新華學院通識教育部,安徽合肥 230088)

    具有標準發(fā)生率和脈沖干擾的SIRS模型穩(wěn)定性分析

    馬艷麗,李海霞,褚正清,聶東明
    (安徽新華學院通識教育部,安徽合肥 230088)

    本文同時考慮脈沖接種和脈沖剔除策略,建立一個具有標準發(fā)生率的SIRS傳染病模型,從理論分析和數(shù)值模擬方面研究了SIRS傳染病模型的動力學性質。首先,得到模型無病T周期解的存在性和疾病流行與否的閾值-基本再生數(shù)R0;其次,應用Floquet定理證明了無病T周期解是局部漸近穩(wěn)定的;然后,利用脈沖微分不等式證明了無病T周期解是全局漸近穩(wěn)定的;最后,進行計算機數(shù)值模擬來進一步驗證理論結果的正確性。

    SIRS模型;標準發(fā)生率;基本再生數(shù);全局漸近穩(wěn)定性;脈沖微分方程

    傳染病是人類的健康大敵,世界衛(wèi)生組織2014年發(fā)布《世界衛(wèi)生統(tǒng)計報告》顯示,現(xiàn)在每年新增艾滋病毒感染者約250萬,這與20年前的每年300萬相比有所下降,但新增感染者仍比死于艾滋病的人數(shù)多8萬,因此艾滋病人的絕對人數(shù)仍在增加.同時,艾滋病對高感染率國家成年人的死亡率產(chǎn)生了重要影響,比如在南非,該國的男性和女性平均預期壽命在1990年均為63歲,但到2011年下降為58歲,另一個非洲國家津巴布韋的人均壽命下降幅度更大,為6歲.近年來,隨著環(huán)境的污染,如肺結核、性病、禽流感等疾病再次抬頭蔓延,而一些新出現(xiàn)的傳染病,如艾滋病(AIDS)、傳染性非典型肺炎(SARS)等也來勢兇猛,全球范圍內(nèi)的傳染病防控形勢依然嚴峻.

    為了對傳染病進行預防和控制,人們常采用免疫接種策略,免疫接種有兩種方式:連續(xù)接種和脈沖接種.麻疹、病毒性肝炎、脊髓灰質炎等人類長期面對的傳染病,常采用連續(xù)接種策略.對于突發(fā)傳染病,如重癥急性呼吸綜合征(SARS)、流行性感冒等,則適合采用脈沖免疫接種策略.關于脈沖免疫接種模型已有大量文獻研究[1-15].剔除也是預防和控制傳染病的重要措施,例如人畜共患的傳染病禽流感、結核病、破傷風和輪狀病毒感染等都可以采取剔除的方法加以控制,關于剔除的傳染病模型在文獻[16-17]中被研究過.但是將這些預防策略共同考慮的模型鮮見報道,受上述文獻[1-17]的啟發(fā),本文同時考慮脈沖接種策略和脈沖剔除策略,建立了一個具有標準發(fā)生率的SIRS傳染病模型,討論了無病T周期解的局部穩(wěn)定性和全局穩(wěn)定性,并利用計算機軟件進行了數(shù)值模擬.

    1 基本模型

    模型假設.將種群分為三類:易感者類,記為S(t),表示t時刻時易感者的數(shù)量;染病者類,記為I(t),表示t時刻時染病者的數(shù)量;恢復者類,記為R(t),表示t時刻時恢復者的數(shù)量.b表示出生率;β表示感染率;μ表示自然死亡率;α表示染病者因病死亡率;γ表示染病者到恢復者的轉移率;δ表示失去免疫率;p表示比例接種系數(shù);q表示染病者的剔除系數(shù);設b,β,μ,α,γ,δ的值均大于零,T為脈沖周期,根據(jù)流行病動力學倉室建模思想建立如下的脈沖微分方程:

    (1)

    (2)

    2 模型分析

    令N(t)=S(t)+I(t)+R(t),由方程組(1)和(2)得到關于N(t)的方程

    N′(t)=(b-μ)N-αI.

    (3)

    (4)

    (5)

    N′=(b-μ-εi)N.

    因為s+i+r=1,所以僅須考慮下面的方程組

    (6)

    (7)

    2.1 無病T周期解的存在性

    研究無病T周期解的存在性就是尋找當i=0時,滿足方程組(6)和(7)的T周期解.注意到,當i=0時,方程組(6)和(7)變?yōu)?/p>

    其在區(qū)間tn≤t≤tn+1上的解為

    rn+1=p+(1-p)rnexp{-(b+δ)}.

    設F:rn→rn+1是一個映射,滿足

    rn+1=F(rn)=p+(1-p)rnexp{-(b+δ)}.

    該映射有唯一的不動點

    2.2 無病T周期解的穩(wěn)定性

    設Φ(t)是該線性系統(tǒng)的基解矩陣,且滿足Φ(0)=I,I為單位陣,基解矩陣為

    其中,

    φ22(t)=exp{-(b+δ)t},

    當t=tn時,得

    由Floquet定理得,無病T周期解穩(wěn)定的充分必要條件是矩陣M的特征值的模都小于1,即只須φ11(T)<1,也就是

    定義

    則R0是基本再生數(shù),并由上面的討論可得下面的定理.

    定理1 當R0<1時,系統(tǒng)(1)和(2)的無病T周期解是局部漸近穩(wěn)定的.

    定理2 當R0<1時,系統(tǒng)(1)和(2)的無病T周期解是全局漸近穩(wěn)定的.

    為簡化定理2的證明,引入如下引理.

    引理1 設f(t),g(t)∈C1[0,+),且則

    由方程(1)和(2)的第一個方程得

    (8)

    對(8)應用微分脈沖不等式可得

    從而有

    由方程(1)和(2)的第二方程得

    (9)

    (10)

    對(10)應用脈沖微分方程不等式得

    因為

    這里

    且有

    (11)

    將式(7)和(8)代入到式(6)可得

    (12)

    則可以得到

    其中

    D1(t)是正的且有上界,當R0<1時,有I(t)→0(t→+∞).

    事實上,令

    則當t≠tn時,有

    當t=tk時,

    應用微分脈沖等式,由(11)和(12)可得

    (13)

    對于式(13)的第二項,有

    綜上所述,當R0<1時,無病T周期解是全局漸近穩(wěn)定的.

    3 數(shù)值模擬

    通過對具有脈沖接種與脈沖剔除的SIRS流行病模型的動力學性態(tài)研究,得到了無病平衡T周期解的穩(wěn)定性.以下用數(shù)值模擬當基本再生數(shù)R0<1時,系統(tǒng)(1)和(2)的解最終穩(wěn)定于無病T周期解時的情況,從而驗證理論結果的正確性,計算機的模擬結果如圖1、圖2和圖3所示.

    圖1 R0<1時,易感者人數(shù)S隨時間t的變化曲線

    圖2 R0<1時,染病者人數(shù)I隨時間t的變化曲線

    圖3 R0<1時,恢復且免疫者人數(shù)R隨時間t的變化曲線

    4 結論與建議

    本文同時考慮脈沖接種策略和脈沖剔除策略,建立一個具有標準發(fā)生率的SIRS傳染病模型,從理論分析和數(shù)值模擬方面研究了SIRS傳染病模型的動力學性質.通過對模型的穩(wěn)定性分析,得到了控制疾病再生和消亡的重要閾值——基本再生數(shù)R0,證明了無病T周期解的存在性,并證明了當R0<1時,無病T周期解是全局漸近穩(wěn)定的,疾病最終消亡.通過對具有脈沖接種與脈沖剔除的SIRS流行病模型的穩(wěn)定性研究,為該類傳染病的防治決策提供了理論基礎和數(shù)量依據(jù).在傳染病模型中,考慮疾病發(fā)生率為標準發(fā)生率,對其他形式的發(fā)生率如一般形式接觸率或非線性等問題有待進一步研究.

    [1]D’Onofrio A. Stability properties of pulse vaccination strategy in SEIR epidemic model[J]. Mathematics Biosciences,2002,179(1):57-72.

    [2]Stone L.Theoretical examination of the pulse vaccination policy in the SIR epidemic model[J].Math Modeling,2000,31(4-5):207-215.

    [3]章培軍,李維德,朱凌峰.SIRS傳染病模型的連續(xù)接種和脈沖接種的比較[J].蘭州大學學報:自然科學版,2011,47(1):82-86.

    [4]朱璣,李維德,朱凌峰.基于SIR傳染病模型的不同控制策略比較[J].北華大學學報:自然科學版,2011,12(3):15-21.

    [5]Jianquan Li,Yanli Yang.SIR-SVS epidemic models with continuous and impulsive vaccination stratgies [J].Journal of Theoretical Biology,2011,280(1):108-116.

    [6]Lijie Hao, Guirong Jiang, Suyu Liu,et al.Global dynamics of an SIRS epidemic model with saturation incidence[J].Biosystems,2013,114(1):56-63.

    [7]G.P.Samanta.Analysis of a delayed epidemic model with pulse vaccination[J].Chaos,Solitions & Fractals,2014,66(1):74-85.

    [8]黃燦云,楊小光,汪訓洋.一類具有時滯和雙線性發(fā)生率的脈沖接種SVIR傳染病模型[J].蘭州理工大學學報,2012,38(4):131-135.

    [9]楊金根,李學志.帶脈沖接種和時滯的乙肝模型的穩(wěn)定性分析[J].數(shù)學的實踐與認識,2013,12(1):141-149.

    [10]Peilin Shi,Lingzhen Dong.Dynamical models for infectious diseases with varying population size and vaccination[J].Journal of Applied Mathematics,2012(1):253-273.

    [11]Wenjie Qin,Sanyi Tang,Robert A.Nonlinear pulse vaccination in an SIR epidemic model with resource limitation[J].Abstravt and Applied Analysis,2013(1):24-37.

    [12]Pei Yong-zhen,Li Shu-ping,Li Chang-guo,et al.The effect of constant and pulse vaccination on an SIR epidemic model with infectious period[J].Applied Mathematical Modeling,2011,35(8):3866-3878.

    [13]Zeng Guang-chao,Chen Lan-sun,Sun Li-hua.Complexity of an SIR epidemic dynamic model with impulsive vaccination control[J].Chaos Solitions and Fractals,2005,26(2):495-505.

    [14]Lu Zhong-hua,Cui Xue-bin,Chen Lan-sun.The effect of constant and pulse vaccination on an SIR epidemic model with horizontal and vertical transmission[J].Mathematical and Computer Modeling,2002, 36(9/10):1039-1057.

    [15]Helong Liu,Jinyuan Yu,Guangtian Zhu.Global stability of an age-structured SIR epidemic model with impulsive vaccination strategy[J].Journal of Systems Science Complexity,2012,25(3):417-429.

    [16]馬艷麗,徐文雄,張仲華.具有一般形式接觸率的SEIR模型的穩(wěn)定性分析[J].中國科學技術大學學報,2015,45(1):737-744.

    [17]Zhang Zhen,Jin Zhen.Stability properties in an SIR epidemic model with asynchronous pulse vaccination and pulse elimination[J].Journal of Taiyuan Normal University(Natural Science Edition), 2006,54(4):8-10.

    Stability Analysis of an SIRS Model With Standard Incidence and Impulsive Perturbations

    MA Yan-li,LI Hai-xia,CHU Zheng-qing,NIE Dong-ming

    (General Education Department,Anhui Xinhua University,Hefei Anhui 230088,China)

    In this paper, impulsive vaccination and impulsive elimination are considered in an SIRS model. The dynamical behavior of an SIRS epidemic model with standard incidence is discussed by means of both theoretical and numerical ways. Firstly, the disease-freeTperiodic solution and the threshold, basic reproductive numberR0which determines whether a disease is extinct or not , are obtained.Secondly, the disease-freeTperiodic solution is locally asymptotically stable by Floquet theorem. Thirdly, the disease-freeTperiodic solution is globally asymptotically stable by impulsive differention in equation. Finally, numerical simulation is given to illustrate the theoretical analysis.

    SIRS model; standard incidence; basic reproductive number; global stability; impulsive differention equation

    2016-12-30

    安徽省高校優(yōu)秀青年人才支持計劃項目“幾類具有潛伏期和混合控制策略傳染病模型的動力學性態(tài)研究”(gxyq2017125);安徽省高校自然科學重點研究項目“關于常曲率空間中基本凸體的幾何不等式理論研究”(KJ2016A310);安徽新華學院校級自然科學重點研究項目“若干分子圖的拓撲指標及其逆問題研究”(2016zr003);安徽省教學研究項目“以數(shù)學建模競賽為契機——高等數(shù)學分層模塊化教學在民辦應用型本科院校中的探索與實踐”(2016jyxm0481);安徽省精品資源共享課程“高等數(shù)學”(2016gxk061);安徽新華學院校級教學團隊“高等數(shù)學”(2016jxtdx03)。

    馬艷麗(1983- ),女,講師,碩士,從事生物數(shù)學研究。

    O175

    A

    2095-7602(2017)06-0001-08

    猜你喜歡
    染病傳染病脈沖
    他們使阿秒光脈沖成為可能
    《傳染病信息》簡介
    傳染病信息(2022年3期)2022-07-15 08:25:08
    偶感
    傳染病的預防
    肝博士(2022年3期)2022-06-30 02:48:50
    脈沖離散Ginzburg-Landau方程組的統(tǒng)計解及其極限行為
    3種傳染病出沒 春天要格外提防
    呼吸道傳染病為何冬春多發(fā)
    均勻網(wǎng)絡上SIR模型三種不同逼近方法比較
    黃芩苷脈沖片的制備
    中成藥(2017年12期)2018-01-19 02:06:54
    愛 情
    詩選刊(2016年9期)2016-11-26 13:47:43
    成人亚洲精品一区在线观看| 亚洲在线自拍视频| 久久香蕉国产精品| 99re在线观看精品视频| 在线观看www视频免费| 久久精品亚洲精品国产色婷小说| 亚洲精品国产一区二区精华液| 在线观看免费午夜福利视频| 在线天堂中文资源库| 一级黄色大片毛片| 男人操女人黄网站| 99国产极品粉嫩在线观看| 中文字幕精品免费在线观看视频| 黄色a级毛片大全视频| 欧美日韩亚洲国产一区二区在线观看| 视频在线观看一区二区三区| 香蕉丝袜av| www.999成人在线观看| 欧美在线黄色| 欧美+亚洲+日韩+国产| 精品不卡国产一区二区三区| 久久久久免费精品人妻一区二区 | 欧美一级毛片孕妇| 九色国产91popny在线| 99国产精品99久久久久| 黄色成人免费大全| 亚洲国产精品sss在线观看| 亚洲激情在线av| a级毛片a级免费在线| 久久亚洲真实| 亚洲天堂国产精品一区在线| av免费在线观看网站| 日本五十路高清| 久久99热这里只有精品18| 成人亚洲精品一区在线观看| 在线观看午夜福利视频| 国产在线精品亚洲第一网站| 亚洲色图av天堂| 成人av一区二区三区在线看| 色综合亚洲欧美另类图片| 国产精品日韩av在线免费观看| 国产成人精品久久二区二区91| 后天国语完整版免费观看| 久久久国产成人精品二区| 亚洲第一电影网av| 国产成人精品久久二区二区免费| 久久久国产欧美日韩av| 亚洲精品美女久久久久99蜜臀| 成在线人永久免费视频| 一边摸一边抽搐一进一小说| 无人区码免费观看不卡| 成人18禁在线播放| 韩国av一区二区三区四区| av免费在线观看网站| 久久欧美精品欧美久久欧美| 好男人电影高清在线观看| 国产成人系列免费观看| 波多野结衣高清作品| 日本免费a在线| 麻豆av在线久日| 看黄色毛片网站| 国产乱人伦免费视频| 亚洲专区中文字幕在线| 男女床上黄色一级片免费看| 黄频高清免费视频| 99精品久久久久人妻精品| 97超级碰碰碰精品色视频在线观看| 黄频高清免费视频| 脱女人内裤的视频| 午夜激情福利司机影院| 啪啪无遮挡十八禁网站| 日本免费a在线| 性色av乱码一区二区三区2| 国产91精品成人一区二区三区| 中文字幕人成人乱码亚洲影| 国产亚洲精品第一综合不卡| 999久久久精品免费观看国产| 激情在线观看视频在线高清| 欧美日本视频| 国产成人精品久久二区二区免费| 高潮久久久久久久久久久不卡| 久久久久免费精品人妻一区二区 | 日韩免费av在线播放| 国产不卡一卡二| 女生性感内裤真人,穿戴方法视频| 久久亚洲真实| 亚洲欧美精品综合久久99| 黑人欧美特级aaaaaa片| 两性午夜刺激爽爽歪歪视频在线观看 | 国产高清videossex| 欧美精品亚洲一区二区| 村上凉子中文字幕在线| 午夜影院日韩av| 日韩欧美三级三区| 18美女黄网站色大片免费观看| 欧美日韩一级在线毛片| 成人18禁高潮啪啪吃奶动态图| 看黄色毛片网站| 欧美成人午夜精品| 午夜影院日韩av| 午夜福利欧美成人| av视频在线观看入口| 久9热在线精品视频| 日日爽夜夜爽网站| 亚洲 欧美 日韩 在线 免费| 精品久久久久久,| 欧美黑人欧美精品刺激| 搡老妇女老女人老熟妇| 操出白浆在线播放| 日本一本二区三区精品| 一本综合久久免费| 999精品在线视频| 国产高清激情床上av| 免费看美女性在线毛片视频| 最近在线观看免费完整版| 人妻久久中文字幕网| 午夜日韩欧美国产| 一进一出抽搐gif免费好疼| 在线永久观看黄色视频| 欧美黑人巨大hd| 午夜福利在线观看吧| 久久久国产成人精品二区| 波多野结衣巨乳人妻| 少妇 在线观看| 人人妻人人澡欧美一区二区| 变态另类丝袜制服| 成人手机av| 成人三级做爰电影| 亚洲精品av麻豆狂野| 亚洲 欧美 日韩 在线 免费| 国产精品一区二区三区四区久久 | 亚洲激情在线av| 人成视频在线观看免费观看| 国产v大片淫在线免费观看| 国产99久久九九免费精品| 欧美日韩精品网址| 亚洲国产精品sss在线观看| 俄罗斯特黄特色一大片| 夜夜夜夜夜久久久久| 欧美绝顶高潮抽搐喷水| 久久狼人影院| 国产成人精品无人区| 免费在线观看亚洲国产| 精品国产国语对白av| 不卡一级毛片| 亚洲av美国av| 给我免费播放毛片高清在线观看| 亚洲一卡2卡3卡4卡5卡精品中文| 亚洲欧美日韩高清在线视频| 不卡av一区二区三区| 神马国产精品三级电影在线观看 | 午夜免费观看网址| 久久久国产成人免费| 亚洲天堂国产精品一区在线| 亚洲电影在线观看av| 免费av毛片视频| 亚洲一码二码三码区别大吗| 亚洲国产中文字幕在线视频| 一本大道久久a久久精品| 人人妻,人人澡人人爽秒播| 久久久水蜜桃国产精品网| 国产成人影院久久av| 熟女少妇亚洲综合色aaa.| 超碰成人久久| 国产不卡一卡二| 亚洲精品中文字幕一二三四区| aaaaa片日本免费| 中文字幕av电影在线播放| 欧美三级亚洲精品| 亚洲精品美女久久久久99蜜臀| 大香蕉久久成人网| 成人18禁高潮啪啪吃奶动态图| 久久国产乱子伦精品免费另类| 波多野结衣高清作品| 91av网站免费观看| av片东京热男人的天堂| 我的亚洲天堂| 精品第一国产精品| 国产一级毛片七仙女欲春2 | 99在线人妻在线中文字幕| 欧美乱妇无乱码| 最新美女视频免费是黄的| 少妇熟女aⅴ在线视频| 伊人久久大香线蕉亚洲五| 久久中文看片网| 亚洲七黄色美女视频| 老司机午夜十八禁免费视频| 亚洲午夜理论影院| 男女午夜视频在线观看| 满18在线观看网站| 国产亚洲欧美精品永久| 日日摸夜夜添夜夜添小说| 91在线观看av| 999久久久精品免费观看国产| 久久久久久久久久黄片| 久9热在线精品视频| 欧美国产精品va在线观看不卡| 午夜免费成人在线视频| 国内毛片毛片毛片毛片毛片| 在线av久久热| 久久精品亚洲精品国产色婷小说| 亚洲国产精品sss在线观看| 日本a在线网址| 精品一区二区三区四区五区乱码| 国产成+人综合+亚洲专区| 午夜免费成人在线视频| 免费女性裸体啪啪无遮挡网站| 久久精品国产99精品国产亚洲性色| 国产熟女xx| 国产男靠女视频免费网站| 免费高清在线观看日韩| 十分钟在线观看高清视频www| 欧美乱色亚洲激情| 国产黄a三级三级三级人| 久久久久国产一级毛片高清牌| 亚洲电影在线观看av| 欧美激情高清一区二区三区| 国产又黄又爽又无遮挡在线| 欧美在线黄色| 桃色一区二区三区在线观看| 精品一区二区三区视频在线观看免费| 一个人观看的视频www高清免费观看 | 哪里可以看免费的av片| 免费看日本二区| 久久精品91蜜桃| 丝袜美腿诱惑在线| 窝窝影院91人妻| av在线播放免费不卡| 在线观看日韩欧美| 久久精品亚洲精品国产色婷小说| 精品国产一区二区三区四区第35| 给我免费播放毛片高清在线观看| 国产1区2区3区精品| www.自偷自拍.com| a级毛片a级免费在线| 亚洲一区高清亚洲精品| 岛国在线观看网站| 男女做爰动态图高潮gif福利片| 亚洲国产精品成人综合色| 久久草成人影院| 岛国在线观看网站| 色综合亚洲欧美另类图片| 国产三级黄色录像| 国产成人影院久久av| 国产成人av激情在线播放| 一进一出抽搐gif免费好疼| 视频区欧美日本亚洲| 亚洲男人的天堂狠狠| 国产黄片美女视频| 日韩欧美一区视频在线观看| 久久精品成人免费网站| 一本精品99久久精品77| 国产亚洲av高清不卡| 精品乱码久久久久久99久播| 国产极品粉嫩免费观看在线| 国产亚洲精品一区二区www| 天天躁狠狠躁夜夜躁狠狠躁| 最近最新免费中文字幕在线| 黄片播放在线免费| 成人国产一区最新在线观看| 国产黄色小视频在线观看| 法律面前人人平等表现在哪些方面| 亚洲专区国产一区二区| 国产精品久久久av美女十八| 亚洲精品久久国产高清桃花| 性欧美人与动物交配| 看黄色毛片网站| 亚洲国产高清在线一区二区三 | 99riav亚洲国产免费| √禁漫天堂资源中文www| av免费在线观看网站| 亚洲成av片中文字幕在线观看| 无人区码免费观看不卡| 欧美日韩瑟瑟在线播放| 免费在线观看完整版高清| 免费在线观看黄色视频的| 国产成+人综合+亚洲专区| 一区二区日韩欧美中文字幕| 俺也久久电影网| 夜夜躁狠狠躁天天躁| 亚洲成人久久性| 国产精品一区二区三区四区久久 | 久久精品91无色码中文字幕| 黑人欧美特级aaaaaa片| av在线播放免费不卡| 久久午夜亚洲精品久久| 欧美绝顶高潮抽搐喷水| 18禁国产床啪视频网站| 亚洲成人久久爱视频| 亚洲一区二区三区不卡视频| 国产精品香港三级国产av潘金莲| 99riav亚洲国产免费| 黄色丝袜av网址大全| 国产99久久九九免费精品| 国产成人啪精品午夜网站| 99久久精品国产亚洲精品| 级片在线观看| 九色国产91popny在线| 亚洲五月色婷婷综合| 亚洲欧美精品综合一区二区三区| 亚洲片人在线观看| 久久精品成人免费网站| 老汉色av国产亚洲站长工具| 国产成人影院久久av| 一本精品99久久精品77| 又黄又爽又免费观看的视频| 久久精品91蜜桃| 免费电影在线观看免费观看| 日韩高清综合在线| 桃色一区二区三区在线观看| 色播在线永久视频| 久久婷婷成人综合色麻豆| 宅男免费午夜| 亚洲精品粉嫩美女一区| 精品久久久久久久毛片微露脸| 老鸭窝网址在线观看| 午夜福利视频1000在线观看| 亚洲第一电影网av| av电影中文网址| 国内少妇人妻偷人精品xxx网站 | 国产精品一区二区精品视频观看| 久热这里只有精品99| 国产精品自产拍在线观看55亚洲| 后天国语完整版免费观看| 国产在线精品亚洲第一网站| 黄色成人免费大全| 欧美在线黄色| 极品教师在线免费播放| 在线看三级毛片| 亚洲午夜理论影院| 老熟妇仑乱视频hdxx| 亚洲av成人av| 欧美另类亚洲清纯唯美| 亚洲成人久久性| 久久人妻福利社区极品人妻图片| 亚洲天堂国产精品一区在线| 欧美精品啪啪一区二区三区| 女性生殖器流出的白浆| 亚洲第一电影网av| 国产亚洲精品第一综合不卡| 国产精品久久视频播放| 在线观看舔阴道视频| 久久99热这里只有精品18| 黄色片一级片一级黄色片| 白带黄色成豆腐渣| 成人免费观看视频高清| 午夜成年电影在线免费观看| 91字幕亚洲| 最近最新免费中文字幕在线| 一区二区三区国产精品乱码| 久久99热这里只有精品18| 久久这里只有精品19| 久久国产精品影院| 久久久国产精品麻豆| 国产精品永久免费网站| 精品人妻1区二区| 日韩高清综合在线| 亚洲va日本ⅴa欧美va伊人久久| 亚洲国产欧洲综合997久久, | 熟女少妇亚洲综合色aaa.| 色老头精品视频在线观看| 欧美黑人欧美精品刺激| 久99久视频精品免费| 国产精品乱码一区二三区的特点| 老汉色av国产亚洲站长工具| 性欧美人与动物交配| 免费无遮挡裸体视频| 国产精品爽爽va在线观看网站 | 成人三级做爰电影| 一本大道久久a久久精品| 欧美黄色片欧美黄色片| 亚洲第一欧美日韩一区二区三区| 成人av一区二区三区在线看| 1024视频免费在线观看| 久久人妻av系列| 女人高潮潮喷娇喘18禁视频| cao死你这个sao货| 精品国产乱子伦一区二区三区| 亚洲国产毛片av蜜桃av| 国产麻豆成人av免费视频| 曰老女人黄片| 女生性感内裤真人,穿戴方法视频| 午夜久久久在线观看| 成人永久免费在线观看视频| 美女午夜性视频免费| 国产欧美日韩一区二区精品| 欧美日韩瑟瑟在线播放| 欧美性猛交黑人性爽| 日本在线视频免费播放| 国产av在哪里看| 亚洲精华国产精华精| 国产精品98久久久久久宅男小说| a级毛片在线看网站| 日韩大码丰满熟妇| 人妻久久中文字幕网| 国内揄拍国产精品人妻在线 | 免费在线观看影片大全网站| 女警被强在线播放| 国产精品久久久久久亚洲av鲁大| av福利片在线| 国产精品影院久久| 中出人妻视频一区二区| 在线观看一区二区三区| 最新美女视频免费是黄的| 中文字幕人妻丝袜一区二区| 高潮久久久久久久久久久不卡| 亚洲熟妇熟女久久| 亚洲国产日韩欧美精品在线观看 | 在线观看www视频免费| 亚洲色图 男人天堂 中文字幕| 国产高清视频在线播放一区| 精品一区二区三区四区五区乱码| 亚洲国产精品sss在线观看| 亚洲av电影不卡..在线观看| 日韩欧美一区二区三区在线观看| 两个人看的免费小视频| 变态另类丝袜制服| 搡老熟女国产l中国老女人| 丝袜美腿诱惑在线| 免费无遮挡裸体视频| 一边摸一边抽搐一进一小说| 亚洲成人久久性| 精品无人区乱码1区二区| 人人妻人人澡欧美一区二区| 亚洲男人的天堂狠狠| 国产欧美日韩一区二区三| 黄片播放在线免费| 国产三级黄色录像| 又黄又粗又硬又大视频| 大型av网站在线播放| 在线观看66精品国产| 成人手机av| 天堂动漫精品| 亚洲精品中文字幕一二三四区| 国产高清videossex| 欧美+亚洲+日韩+国产| av视频在线观看入口| 99久久99久久久精品蜜桃| 亚洲精品色激情综合| 一卡2卡三卡四卡精品乱码亚洲| 久久久久亚洲av毛片大全| 久久九九热精品免费| 婷婷精品国产亚洲av| 少妇裸体淫交视频免费看高清 | 久久久国产成人免费| 国产亚洲精品久久久久久毛片| 亚洲aⅴ乱码一区二区在线播放 | 免费无遮挡裸体视频| 中文字幕另类日韩欧美亚洲嫩草| 免费在线观看完整版高清| 美女大奶头视频| 操出白浆在线播放| 动漫黄色视频在线观看| 色精品久久人妻99蜜桃| 欧美性长视频在线观看| 亚洲狠狠婷婷综合久久图片| 热re99久久国产66热| 变态另类丝袜制服| 淫妇啪啪啪对白视频| 美女免费视频网站| 非洲黑人性xxxx精品又粗又长| 亚洲成人免费电影在线观看| 欧美在线一区亚洲| 欧美性长视频在线观看| 精品久久久久久久毛片微露脸| 国产99久久九九免费精品| 免费电影在线观看免费观看| 高清在线国产一区| 国产不卡一卡二| 欧美色欧美亚洲另类二区| 999久久久精品免费观看国产| 夜夜爽天天搞| 欧美绝顶高潮抽搐喷水| 久久久久久大精品| 夜夜躁狠狠躁天天躁| 久久中文看片网| 国产激情久久老熟女| 亚洲自拍偷在线| 美女 人体艺术 gogo| 身体一侧抽搐| 18美女黄网站色大片免费观看| 欧美不卡视频在线免费观看 | 精品国产一区二区三区四区第35| 成人亚洲精品av一区二区| 女人高潮潮喷娇喘18禁视频| 麻豆国产av国片精品| 亚洲专区国产一区二区| 国产成人精品无人区| 亚洲va日本ⅴa欧美va伊人久久| 天天添夜夜摸| 99国产精品99久久久久| 欧美在线黄色| 天天添夜夜摸| 成人一区二区视频在线观看| 看黄色毛片网站| 91成人精品电影| 一级作爱视频免费观看| 一二三四在线观看免费中文在| 51午夜福利影视在线观看| 国产在线精品亚洲第一网站| 男人舔女人下体高潮全视频| 久久中文字幕一级| 在线看三级毛片| 无遮挡黄片免费观看| 99久久综合精品五月天人人| 日韩精品青青久久久久久| 久久久国产成人免费| 国产精品综合久久久久久久免费| 欧美成狂野欧美在线观看| 日韩中文字幕欧美一区二区| 两个人视频免费观看高清| 在线观看舔阴道视频| 亚洲无线在线观看| 看黄色毛片网站| 国产精品精品国产色婷婷| 久久久久免费精品人妻一区二区 | 久久久久九九精品影院| 国产精品1区2区在线观看.| 啪啪无遮挡十八禁网站| 欧美精品亚洲一区二区| 1024视频免费在线观看| 久久香蕉精品热| 久久性视频一级片| 丝袜在线中文字幕| 欧美成狂野欧美在线观看| 老司机靠b影院| 日本一区二区免费在线视频| 97碰自拍视频| 国产亚洲欧美98| www日本黄色视频网| 午夜成年电影在线免费观看| 丰满的人妻完整版| 日本一本二区三区精品| 亚洲精品美女久久av网站| 国产1区2区3区精品| 满18在线观看网站| ponron亚洲| 国产真实乱freesex| 久久狼人影院| 精品久久久久久久人妻蜜臀av| cao死你这个sao货| 亚洲全国av大片| 日本在线视频免费播放| 嫩草影院精品99| 一级a爱视频在线免费观看| 国产精品美女特级片免费视频播放器 | www国产在线视频色| 亚洲,欧美精品.| 搞女人的毛片| 久久香蕉激情| av天堂在线播放| 丰满人妻熟妇乱又伦精品不卡| 露出奶头的视频| 国产视频一区二区在线看| 国产97色在线日韩免费| 在线免费观看的www视频| 老司机深夜福利视频在线观看| 真人一进一出gif抽搐免费| 一级a爱片免费观看的视频| 我的亚洲天堂| 欧美在线一区亚洲| 久久精品人妻少妇| 日本三级黄在线观看| videosex国产| 中文亚洲av片在线观看爽| 久久久国产成人精品二区| 老汉色av国产亚洲站长工具| 成年版毛片免费区| 老司机午夜十八禁免费视频| 亚洲精品色激情综合| 日韩精品中文字幕看吧| 午夜老司机福利片| 丰满的人妻完整版| 欧美乱色亚洲激情| 国产真实乱freesex| 国产av一区二区精品久久| 韩国av一区二区三区四区| 欧美乱妇无乱码| av在线播放免费不卡| 中文字幕高清在线视频| 嫁个100分男人电影在线观看| 国产在线观看jvid| 国产1区2区3区精品| www.www免费av| 搡老岳熟女国产| 麻豆一二三区av精品| 十八禁人妻一区二区| 精华霜和精华液先用哪个| 麻豆一二三区av精品| 高清毛片免费观看视频网站| www日本黄色视频网| 中文字幕高清在线视频| 午夜福利视频1000在线观看| 国产成人系列免费观看| 免费一级毛片在线播放高清视频| 亚洲精品中文字幕一二三四区| 天天躁夜夜躁狠狠躁躁| 国产片内射在线| 丝袜美腿诱惑在线| 少妇粗大呻吟视频| 国产高清videossex| 18禁美女被吸乳视频| 久久精品人妻少妇| 黑人欧美特级aaaaaa片| 色播亚洲综合网| 国产精品影院久久| 亚洲欧美日韩高清在线视频| 国产精品 国内视频| 女人高潮潮喷娇喘18禁视频| bbb黄色大片| 亚洲第一青青草原| 久久精品亚洲精品国产色婷小说| 亚洲人成电影免费在线| 久久中文看片网| 欧美不卡视频在线免费观看 |