• <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
    大陆偷拍与自拍| 日韩在线高清观看一区二区三区| 飞空精品影院首页| 国产亚洲精品久久久com| 国产免费一区二区三区四区乱码| 男女国产视频网站| 欧美日韩视频精品一区| 成年人免费黄色播放视频| 国产在线免费精品| av.在线天堂| 亚洲精品视频女| 欧美日韩视频精品一区| 我的女老师完整版在线观看| 色婷婷av一区二区三区视频| 国产精品免费大片| 久久韩国三级中文字幕| 中文字幕精品免费在线观看视频 | 国产高清国产精品国产三级| 免费观看在线日韩| 天天操日日干夜夜撸| 美女主播在线视频| 大片电影免费在线观看免费| 亚洲av电影在线观看一区二区三区| 一区二区三区四区激情视频| 亚洲综合精品二区| 欧美3d第一页| 在线观看人妻少妇| 欧美日本中文国产一区发布| 边亲边吃奶的免费视频| 91aial.com中文字幕在线观看| 亚洲综合色惰| 亚洲av电影在线观看一区二区三区| 一区二区三区免费毛片| 午夜激情福利司机影院| 青春草国产在线视频| 晚上一个人看的免费电影| 国产av一区二区精品久久| 伦精品一区二区三区| 久久亚洲国产成人精品v| 国产精品一区二区三区四区免费观看| 大陆偷拍与自拍| 欧美日韩在线观看h| 国产成人91sexporn| 国产一区亚洲一区在线观看| 亚洲精品久久午夜乱码| 王馨瑶露胸无遮挡在线观看| 欧美日韩av久久| 久热久热在线精品观看| 伊人久久国产一区二区| 夜夜爽夜夜爽视频| 老女人水多毛片| 国产成人精品一,二区| 久久久a久久爽久久v久久| 精品人妻在线不人妻| 一个人免费看片子| videos熟女内射| 在线亚洲精品国产二区图片欧美 | 久久久久国产网址| 亚洲av二区三区四区| 大陆偷拍与自拍| 久久精品人人爽人人爽视色| 精品熟女少妇av免费看| 成人毛片a级毛片在线播放| 交换朋友夫妻互换小说| 国产精品成人在线| 亚洲精品色激情综合| 自拍欧美九色日韩亚洲蝌蚪91| 大片免费播放器 马上看| 亚洲精品成人av观看孕妇| 18禁观看日本| 人人妻人人澡人人爽人人夜夜| av免费观看日本| 9色porny在线观看| 久久久精品区二区三区| 国产成人午夜福利电影在线观看| 欧美日韩视频精品一区| 丝袜脚勾引网站| 一个人免费看片子| 男女啪啪激烈高潮av片| 精品一区在线观看国产| 亚洲av综合色区一区| 日日摸夜夜添夜夜添av毛片| 美女主播在线视频| 性色avwww在线观看| 亚洲精品一区蜜桃| 日本黄大片高清| 亚州av有码| av电影中文网址| 新久久久久国产一级毛片| 免费高清在线观看日韩| 国产欧美另类精品又又久久亚洲欧美| 亚洲精品国产av成人精品| 中文字幕人妻熟人妻熟丝袜美| .国产精品久久| 亚洲怡红院男人天堂| av国产久精品久网站免费入址| 乱人伦中国视频| 最近中文字幕2019免费版| 亚洲欧洲国产日韩| 亚洲欧美日韩卡通动漫| 国产爽快片一区二区三区| 成人亚洲精品一区在线观看| 久久久久视频综合| 久久国内精品自在自线图片| 久久精品熟女亚洲av麻豆精品| 两个人的视频大全免费| 日日撸夜夜添| 草草在线视频免费看| 丁香六月天网| 久久精品国产亚洲网站| 啦啦啦啦在线视频资源| 国产亚洲精品第一综合不卡 | 18禁动态无遮挡网站| 精品一区二区三卡| 久久久国产精品麻豆| 日韩人妻高清精品专区| www.色视频.com| 久久影院123| 我要看黄色一级片免费的| 晚上一个人看的免费电影| 国产不卡av网站在线观看| 国产精品一区二区在线观看99| 亚洲av日韩在线播放| 2018国产大陆天天弄谢| 亚洲精品中文字幕在线视频| 国产男女超爽视频在线观看| 中文精品一卡2卡3卡4更新| 国产成人aa在线观看| a级片在线免费高清观看视频| 只有这里有精品99| 大片免费播放器 马上看| 涩涩av久久男人的天堂| 18在线观看网站| 黄片播放在线免费| 极品少妇高潮喷水抽搐| 天天躁夜夜躁狠狠久久av| 制服人妻中文乱码| 久久久久久久久久久丰满| 国产精品一区二区在线观看99| 国产精品国产三级国产av玫瑰| 精品国产国语对白av| 亚洲精品中文字幕在线视频| 精品人妻在线不人妻| 国产熟女欧美一区二区| 五月玫瑰六月丁香| 性色av一级| 国产精品.久久久| 亚洲综合色惰| 青春草国产在线视频| 看免费成人av毛片| 又粗又硬又长又爽又黄的视频| 欧美日韩av久久| 亚洲人成网站在线观看播放| 草草在线视频免费看| 2021少妇久久久久久久久久久| 精品视频人人做人人爽| 一级,二级,三级黄色视频| 中文字幕最新亚洲高清| 久久精品夜色国产| 赤兔流量卡办理| 国产精品久久久久久久电影| 男人操女人黄网站| 一级毛片电影观看| 精品一区二区三区视频在线| 午夜免费观看性视频| 交换朋友夫妻互换小说| 精品久久久久久久久av| 能在线免费看毛片的网站| 51国产日韩欧美| 日本欧美视频一区| 亚洲经典国产精华液单| 成人黄色视频免费在线看| 欧美另类一区| 2021少妇久久久久久久久久久| 熟妇人妻不卡中文字幕| 高清视频免费观看一区二区| 母亲3免费完整高清在线观看 | 午夜激情福利司机影院| 日韩一区二区视频免费看| 亚洲精品亚洲一区二区| 国产高清有码在线观看视频| 久久狼人影院| av在线观看视频网站免费| 亚洲色图综合在线观看| 天天影视国产精品| 久久精品久久久久久噜噜老黄| 国产精品久久久久成人av| 国产又色又爽无遮挡免| 亚洲精品视频女| 亚洲国产av新网站| 色94色欧美一区二区| 亚洲av.av天堂| 啦啦啦视频在线资源免费观看| 日本黄色片子视频| 午夜福利,免费看| 男女边摸边吃奶| 国产成人精品婷婷| 在线 av 中文字幕| 精品亚洲成a人片在线观看| 十八禁网站网址无遮挡| 最近中文字幕2019免费版| 亚洲三级黄色毛片| 久久国产精品大桥未久av| 美女内射精品一级片tv| 天美传媒精品一区二区| 精品99又大又爽又粗少妇毛片| 99国产精品免费福利视频| 亚洲色图 男人天堂 中文字幕 | 老司机影院毛片| 五月玫瑰六月丁香| 国产精品嫩草影院av在线观看| 97在线人人人人妻| 午夜福利在线观看免费完整高清在| 久久久国产一区二区| 亚洲精品成人av观看孕妇| 成年人免费黄色播放视频| 69精品国产乱码久久久| 大片免费播放器 马上看| 国产深夜福利视频在线观看| 色网站视频免费| www.色视频.com| 校园人妻丝袜中文字幕| 中文字幕制服av| 亚洲在久久综合| 亚洲国产av影院在线观看| 亚洲综合精品二区| 青春草视频在线免费观看| 黑丝袜美女国产一区| 99热国产这里只有精品6| 少妇丰满av| 9色porny在线观看| 免费黄色在线免费观看| 久久久久久久久久成人| 69精品国产乱码久久久| 嫩草影院入口| av在线老鸭窝| a级毛片在线看网站| 男女国产视频网站| 久久久精品区二区三区| 国产精品.久久久| 婷婷色综合大香蕉| 国产国语露脸激情在线看| 国产成人a∨麻豆精品| 婷婷成人精品国产| 伊人久久精品亚洲午夜| 国产精品一国产av| 777米奇影视久久| 色5月婷婷丁香| 欧美日本中文国产一区发布| 国产精品 国内视频| 亚洲精品色激情综合| 大香蕉97超碰在线| 欧美老熟妇乱子伦牲交| 免费观看a级毛片全部| 啦啦啦在线观看免费高清www| 日本黄大片高清| 成人影院久久| 免费看av在线观看网站| 日韩,欧美,国产一区二区三区| 视频区图区小说| 久久国内精品自在自线图片| 热99久久久久精品小说推荐| av免费观看日本| 91在线精品国自产拍蜜月| 丝袜美足系列| 丝瓜视频免费看黄片| 欧美老熟妇乱子伦牲交| 一区在线观看完整版| 国内精品宾馆在线| 一个人免费看片子| 大话2 男鬼变身卡| 欧美精品国产亚洲| 在线 av 中文字幕| 黄色欧美视频在线观看| 国产高清不卡午夜福利| 色婷婷av一区二区三区视频| 国产精品久久久久成人av| 青春草视频在线免费观看| 母亲3免费完整高清在线观看 | 色婷婷久久久亚洲欧美| 一区二区日韩欧美中文字幕 | 99国产精品免费福利视频| 3wmmmm亚洲av在线观看| 久久午夜福利片| 精品亚洲成国产av| 纵有疾风起免费观看全集完整版| 亚洲内射少妇av| 亚洲精品中文字幕在线视频| 亚洲成人av在线免费| 国产精品免费大片| 欧美亚洲日本最大视频资源| 男男h啪啪无遮挡| 在线观看国产h片| 99久久人妻综合| 这个男人来自地球电影免费观看 | a级毛片黄视频| 80岁老熟妇乱子伦牲交| 精品久久久久久电影网| 免费日韩欧美在线观看| 伊人久久精品亚洲午夜| 国产成人91sexporn| 少妇 在线观看| 久久精品人人爽人人爽视色| 大片电影免费在线观看免费| 人妻系列 视频| 国产精品一区二区三区四区免费观看| 精品久久久久久电影网| 一本久久精品| 国产又色又爽无遮挡免| 少妇被粗大猛烈的视频| 亚洲四区av| 久久久精品免费免费高清| 亚洲精品,欧美精品| 美女xxoo啪啪120秒动态图| 久久精品国产亚洲网站| 免费看不卡的av| 日本av免费视频播放| 国产日韩一区二区三区精品不卡 | 伊人久久国产一区二区| 精品熟女少妇av免费看| 91久久精品国产一区二区三区| 999精品在线视频| 亚洲人成77777在线视频| 日韩中文字幕视频在线看片| 国产片特级美女逼逼视频| 久久av网站| 男的添女的下面高潮视频| 日韩,欧美,国产一区二区三区| 中文字幕人妻丝袜制服| 3wmmmm亚洲av在线观看| 亚洲国产最新在线播放| 久久狼人影院| 欧美日韩视频精品一区| av电影中文网址| 午夜免费男女啪啪视频观看| 国产欧美日韩综合在线一区二区| 精品一区二区三区视频在线| 欧美+日韩+精品| 亚洲国产av影院在线观看| a 毛片基地| 日韩,欧美,国产一区二区三区| 久久精品熟女亚洲av麻豆精品| 在线天堂最新版资源| 青青草视频在线视频观看| 精品一区二区三卡| 精品人妻偷拍中文字幕| 欧美性感艳星| 国产日韩欧美亚洲二区| 亚洲精品日韩av片在线观看| 在线 av 中文字幕| 亚洲国产日韩一区二区| 中国三级夫妇交换| 丰满饥渴人妻一区二区三| 夜夜看夜夜爽夜夜摸| 午夜免费鲁丝| 精品久久久噜噜| 最黄视频免费看| 嫩草影院入口| 亚洲久久久国产精品| 777米奇影视久久| 999精品在线视频| a级毛片在线看网站| 国产在视频线精品| 久久久久久久久久成人| 爱豆传媒免费全集在线观看| 中文字幕亚洲精品专区| 日韩一本色道免费dvd| 国产成人a∨麻豆精品| 91精品一卡2卡3卡4卡| 亚洲av国产av综合av卡| 日日摸夜夜添夜夜爱| 中文字幕久久专区| 午夜老司机福利剧场| 日韩 亚洲 欧美在线| 80岁老熟妇乱子伦牲交| 日本与韩国留学比较| 成人漫画全彩无遮挡| 欧美xxxx性猛交bbbb| 高清视频免费观看一区二区| 午夜日本视频在线| 精品少妇久久久久久888优播| 丁香六月天网| 成人手机av| 春色校园在线视频观看| 中国美白少妇内射xxxbb| 国产伦精品一区二区三区视频9| 国产 精品1| av有码第一页| 精品久久国产蜜桃| 高清毛片免费看| 亚洲激情五月婷婷啪啪| 国产在视频线精品| 亚洲国产日韩一区二区| 国国产精品蜜臀av免费| 久久精品国产a三级三级三级| a级毛片在线看网站| 亚洲婷婷狠狠爱综合网| kizo精华| 亚洲一级一片aⅴ在线观看| 亚洲精品国产av蜜桃| 亚洲精品久久久久久婷婷小说| 日韩免费高清中文字幕av| 日日摸夜夜添夜夜添av毛片| 欧美人与性动交α欧美精品济南到 | 婷婷色av中文字幕| 免费大片黄手机在线观看| 亚洲av福利一区| 久久午夜福利片| 亚洲精品乱码久久久v下载方式| 久久久久久久国产电影| 汤姆久久久久久久影院中文字幕| 大片免费播放器 马上看| 国产毛片在线视频| 国产国拍精品亚洲av在线观看| 久久久a久久爽久久v久久| 一区二区av电影网| 飞空精品影院首页| 国产黄频视频在线观看| 亚洲国产色片| 国产精品一区二区在线观看99| 中国国产av一级| 精品视频人人做人人爽| 日本wwww免费看| 中文字幕制服av| 91午夜精品亚洲一区二区三区| 五月天丁香电影| 欧美日韩综合久久久久久| 91精品国产国语对白视频| 亚洲久久久国产精品| 久久影院123| 秋霞伦理黄片| 女人精品久久久久毛片| 国产爽快片一区二区三区| 免费高清在线观看视频在线观看| av女优亚洲男人天堂| 性高湖久久久久久久久免费观看| 天美传媒精品一区二区| 国产精品不卡视频一区二区| 国产亚洲午夜精品一区二区久久| 国产色爽女视频免费观看| 少妇人妻久久综合中文| 欧美成人精品欧美一级黄| 亚洲精品日韩在线中文字幕| av女优亚洲男人天堂| 亚洲国产精品一区二区三区在线| 午夜91福利影院| 人妻制服诱惑在线中文字幕| 久久久精品94久久精品| 丁香六月天网| 亚洲美女黄色视频免费看| 久久人人爽人人片av| 天天操日日干夜夜撸| 热99国产精品久久久久久7| 久久久国产精品麻豆| xxxhd国产人妻xxx| 欧美老熟妇乱子伦牲交| 激情五月婷婷亚洲| 久久久国产欧美日韩av| 亚洲av中文av极速乱| 国产片特级美女逼逼视频| 亚洲经典国产精华液单| 国产综合精华液| 99久久精品一区二区三区| 国产成人免费观看mmmm| 九九爱精品视频在线观看| 国产精品无大码| 建设人人有责人人尽责人人享有的| 日产精品乱码卡一卡2卡三| 亚洲精品成人av观看孕妇| 亚洲一区二区三区欧美精品| 精品人妻熟女av久视频| 天天操日日干夜夜撸| 五月天丁香电影| 美女脱内裤让男人舔精品视频| 18禁在线无遮挡免费观看视频| 女人久久www免费人成看片| 另类精品久久| freevideosex欧美| 黄色视频在线播放观看不卡| 欧美三级亚洲精品| 国语对白做爰xxxⅹ性视频网站| 午夜免费男女啪啪视频观看| 五月玫瑰六月丁香| 18禁动态无遮挡网站| 黄片无遮挡物在线观看| 国产日韩欧美视频二区| 精品人妻偷拍中文字幕| 美女国产高潮福利片在线看| 久久99一区二区三区| 最近手机中文字幕大全| 大陆偷拍与自拍| 亚洲成人手机| 99国产精品免费福利视频| 成人影院久久| 久久精品熟女亚洲av麻豆精品| 日本vs欧美在线观看视频| 亚洲精品一区蜜桃| 欧美bdsm另类| 婷婷成人精品国产| 男人操女人黄网站| 欧美 日韩 精品 国产| 最黄视频免费看| 中文字幕免费在线视频6| 国产一区二区在线观看日韩| 亚洲国产欧美日韩在线播放| 国产精品成人在线| 秋霞在线观看毛片| 22中文网久久字幕| 国产精品99久久99久久久不卡 | 18+在线观看网站| 国国产精品蜜臀av免费| 久久这里有精品视频免费| 777米奇影视久久| 少妇精品久久久久久久| 久久久久久久久大av| 精品人妻一区二区三区麻豆| 少妇人妻久久综合中文| 97在线人人人人妻| 桃花免费在线播放| 亚洲欧美色中文字幕在线| 亚洲精品色激情综合| 欧美97在线视频| 一本色道久久久久久精品综合| 日日啪夜夜爽| 日韩制服骚丝袜av| 国产爽快片一区二区三区| 3wmmmm亚洲av在线观看| 如何舔出高潮| 在线观看免费日韩欧美大片 | av.在线天堂| 夜夜骑夜夜射夜夜干| 亚洲av二区三区四区| 王馨瑶露胸无遮挡在线观看| 看免费成人av毛片| 大片免费播放器 马上看| 国产精品一区二区在线观看99| 91精品伊人久久大香线蕉| 精品少妇黑人巨大在线播放| 午夜免费鲁丝| 久久精品国产亚洲av涩爱| 亚洲,一卡二卡三卡| 成年人免费黄色播放视频| 亚洲精品av麻豆狂野| 看十八女毛片水多多多| 国产成人a∨麻豆精品| 国产男女超爽视频在线观看| 91在线精品国自产拍蜜月| 麻豆成人av视频| 观看美女的网站| 国产成人freesex在线| 校园人妻丝袜中文字幕| 99re6热这里在线精品视频| 99九九线精品视频在线观看视频| kizo精华| 久久ye,这里只有精品| 日韩电影二区| 精品久久久久久久久av| 亚洲av成人精品一区久久| 五月玫瑰六月丁香| 日韩精品免费视频一区二区三区 | 我的老师免费观看完整版| 人妻少妇偷人精品九色| 久久免费观看电影| 国产日韩欧美视频二区| 观看av在线不卡| 国产日韩欧美亚洲二区| 欧美精品国产亚洲| 99久久精品一区二区三区| 国产精品女同一区二区软件| 午夜视频国产福利| 观看美女的网站| 免费av不卡在线播放| 婷婷色综合www| a级毛片在线看网站| 亚洲成人一二三区av| 亚洲国产av影院在线观看| 精品一区二区三卡| 9色porny在线观看| 欧美日本中文国产一区发布| 国语对白做爰xxxⅹ性视频网站| 人妻夜夜爽99麻豆av| 简卡轻食公司| 视频在线观看一区二区三区| 嘟嘟电影网在线观看| 欧美日韩精品成人综合77777| 嘟嘟电影网在线观看| 99热网站在线观看| 亚洲第一av免费看| 一级毛片电影观看| 99re6热这里在线精品视频| 热re99久久国产66热| 寂寞人妻少妇视频99o| 水蜜桃什么品种好| 久久久午夜欧美精品| 久久精品国产亚洲网站| 欧美人与性动交α欧美精品济南到 | 女性被躁到高潮视频| 国产永久视频网站| 国产精品国产三级国产专区5o| 26uuu在线亚洲综合色| 午夜福利视频在线观看免费| 各种免费的搞黄视频| 国产淫语在线视频| 国产精品免费大片| 精品人妻在线不人妻| 伦理电影免费视频| 日韩精品免费视频一区二区三区 | 亚洲经典国产精华液单| 国产精品久久久久成人av| 久久99一区二区三区| 韩国av在线不卡| 26uuu在线亚洲综合色| 涩涩av久久男人的天堂| 精品午夜福利在线看| 亚洲丝袜综合中文字幕|