• <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)定性與Hopf分岔分析*

    2014-09-17 06:00:50陳方方洪靈
    動力學與控制學報 2014年1期
    關(guān)鍵詞:特征方程平衡點時滯

    陳方方 洪靈

    (西安交通大學機械結(jié)構(gòu)強度與振動國家重點實驗室,西安 710049)

    一類具有時滯和非線性發(fā)生率的SIRS傳染病模型穩(wěn)定性與Hopf分岔分析*

    陳方方 洪靈?

    (西安交通大學機械結(jié)構(gòu)強度與振動國家重點實驗室,西安 710049)

    研究了一類具有時滯及非線性特性發(fā)生率的SIRS傳染病模型,首先利用特征值理論分析了無病平衡點和地方病平衡點的局部穩(wěn)定性;并以時滯τ作為分岔參數(shù),分析了模型的Hopf分岔行為,運用中心流形定理和規(guī)范型理論給出了分岔方向及分岔周期解穩(wěn)定性的計算公式;最后,數(shù)值模擬驗證了理論分析結(jié)果.

    穩(wěn)定性, 時滯, 非線性發(fā)生率, 階段結(jié)構(gòu), Hopf分岔

    引言

    K-M在1927年建立了所謂的“倉室”模型[1]后,“倉室”模型的基本思想方法一直被廣泛的使用并不斷地發(fā)展著,截至目前,學者已根據(jù)不同疾病的發(fā)病機理建立了不同類型的傳染病模型[2],并得出了一些結(jié)論.

    在以前討論的大多數(shù)模型中,我們通常將人群分成易感染者、染病者和移出者等若干類,并假定類型中的個體無差異,即他們在被感染、隔離、接種等方面完全相同,針對這類模型的研究已有很多研究結(jié)果.但是這類模型假設(shè)的合理性與實際疾病的傳播仍存在一定的差距,特別是對某些傳染病或接觸性疾病,在不同的年齡階段,其傳播概率有很大不同,某些類型的疾病在成人中的傳播概率很大,如淋病,AIDS等,而有些疾病,如麻疹,水痘等則在兒童中的傳播概率較大.因此,為了更清楚的描述疾病的傳播機理,就需在相應(yīng)的傳染病模型中考慮階段結(jié)構(gòu).于是,在刻畫傳染病模型時,具有階段結(jié)構(gòu)[3-4]的傳染病模型能更好地反映生物個體的生理特征和疾病的傳播機理,引起了很多學者的關(guān)注,且具有重要的生物學意義.

    依靠媒介傳染病毒的傳染病也非常多,它通過病毒、原蟲或細菌等生物載體傳播的,如瘧疾通過蚊子來傳播.Cooke在文[5]中給出了一種通過媒介傳播疾病的數(shù)學模型,該建模思想在以后的建模中被廣泛應(yīng)用.

    為方便起見,我們假定由于疾病的影響,染病的成熟個體沒有繁殖能力;易感染人群僅僅被媒介染病者所感染,媒介和人群充分混合.在該假設(shè)下,我們基于文[6],建立了具有階段結(jié)構(gòu)和非線性發(fā)生率的SIRS時滯傳染病模型,研究模型的平衡點的局部穩(wěn)定性,分析其復雜的動力學行為.

    1 模型的建立

    本章所建立的SIRS時滯傳染病模型:

    其中,X(t)為t時刻處于不成熟階段的人口密度,S(t),I(t),R(t)為t時刻成年人中易染者、感染者和移除者的密度,ω為不成熟個體到成熟個體的轉(zhuǎn)換率,μ0為不成熟個體的自然死亡率,μ為成熟個體的自然死亡率,σ為感染率,α為抑制率,ε為感染者的因病死亡率,γ為治愈率,δ為已免疫的個體重新獲得易感染能力的概率,τ為疾病的潛伏期.

    2 平衡點的存在性

    由于模型(1)的形式較復雜,為方便討論,我們對模型進行簡化.對(1)進行無量綱化和變量替換:

    得簡化后的模型為:

    模型(2)的初始條件為:

    其中:?4+={x1,x2,x3,x4∈?4:xi≥0,i=1,2,3,4}

    定義基本再生數(shù):

    經(jīng)直接計算得,系統(tǒng)(2)存在無病平衡點:

    對地方病平衡點的存在性,我們有如下性質(zhì):

    性質(zhì) 模型(2)有唯一正平衡點E*=(u*,s*,i*,r*)當且僅當以下條件成立:

    3 平衡點的穩(wěn)定性

    3.1 無病平衡點的局部穩(wěn)定性

    由平衡點的生物學意義得,若存在無病平衡點,則無病平衡點應(yīng)為正,即b1>a.因此,有如下定理:

    定理1 假設(shè)R0<1,對任意的τ≥0,E0局部漸近穩(wěn)定的;R0>1成立,對任意的τ≥0,E0不穩(wěn)定的.證明: 系統(tǒng)(2)在無病平衡點E0附近線性化系統(tǒng)的特征方程為:

    易證,當R0<1時,特征根均具有負實部;當R0>1時,至少存在一個具有正實部的特征根,因此得證.

    3.2 地方病平衡點的局部穩(wěn)定性

    定理2 假設(shè)R0>1成立,且τ=0時,系統(tǒng)(2)的地方病平衡點E*局部漸近穩(wěn)定.

    證明: 當τ=0時,系統(tǒng)(2)在地方病平衡點E*附近線性化系統(tǒng)的特征方程為:

    根據(jù) Routh-Hurwitz判據(jù),易證:當R0>1 時,a1>0,且有H1,H2,H3>0.故 τ=0 時,特征根均具有負實部,地方病平衡點E*局部漸近穩(wěn)定.

    定理3 當 τ>0時,若aβ+b1b2<ab2,則E*條件穩(wěn)定.

    證明: 當τ>0時,系統(tǒng)(2)在地方病平衡點E*附近線性化系統(tǒng)的特征方程為:

    由aβ+b1b2<ab2得R0>1,已有結(jié)論,當 τ=0時系統(tǒng)特征根均具有負實部.以下證明特征方程(3)有唯一一對純虛根 ±iω0(ω0>0),假定對 τ >0,iω0(ω0>0)是(3)的根,代入(3)得:

    分離實部虛部得:

    將上式兩邊分別平方相加得:

    當R0>1時,我們可證得:p>0,q>0.由條件aβ+b1b2<ab2,易得s<0,由引理 2.1[6]得:(6)式存在唯一的正根ω0,即特征方程(3)存在唯一一對純虛根,由此可得到τk>0,使得有一對純虛根.

    由引理2.3[7],于是我們完成了定理3的證明.以下我們將證明:

    這意味著,在τ>τk時,至少存在一個具有正實部的特征根.將特征方程(3)對τ微分可得:

    根據(jù) Rouché's定理[8]得,時滯量 τ 由小于 τk的值增加到大于τk的值時,特征方程(3)的特征根,從虛軸左側(cè)穿過虛軸到達虛軸右側(cè),橫截條件成立,因此,滿足 Hopf分岔條件[9],在 τ=τk處發(fā)生了Hopf分岔.

    4 Hopf分岔方向和分岔穩(wěn)定性

    對模型(2),運用中心流形定理和規(guī)范型理論[10]給出系統(tǒng)(2)的Hopf分岔方向及分岔周期解的穩(wěn)定性和周期計算公式.令:

    為方便起見,去掉“—”,則系統(tǒng)(2)可以寫成C=C([-1,0],R4)上的泛函微分方程:

    f:R×C→R4,Lμ:C→R4分別表示為:由Riesz表示引理得:對 θ∈[0,1],存在一個有界變差函數(shù) η(θ,μ),使得:

    實際上,可選取:

    其中δ(.)表示Dirac Delta函數(shù).

    對于 φ =(φ1,φ2,φ3,φ4)T∈C[-1,0],R4),定義:

    因此,系統(tǒng)(7)可化為:

    和雙線性內(nèi)積:

    顯然,A和A*為共軛算子.±iω0τk是A(0)的特征值,也是A*的特征值.下面計算A和A*關(guān)于iω0τk和 -iω0τk的特征向量.q(θ)=(1,q1,q2,q3)Teiω0τkθ是A(0)特征向量,有A(0)q(θ)=iω0τkq(θ),根據(jù)A(0)定義及式(9)(10)(11)可解得:

    A*對應(yīng)特征向量為q*(s)=D(1,)

    確?!磓*(s),q(θ)〉=1.接下來計算在 μ=0決定中心流形C0的坐標,令μ=0時(12)的解為ut,定義:

    在中心流形C0上,我們有

    z和表示q*和上中心流形C0的局部坐標.

    由式(13)和(14)得:

    連同式(8)得:

    下面需要計算W20(θ)和W11(θ),由式(12)(13)得:

    將式(17)代入式(16),并比較系數(shù)得:

    比較式(15)和式(17)的系數(shù)得:

    由A(0)的定義及式(18)和(19)得:

    類似地,由式(18)和式(19)得:

    接下來,計算E1和E2的值.由A(0)及式(18),有:

    由式(16)得,當θ=0時,

    此時,將式(20)和式(23)代入式(22),可求得E1.類似將式(21)和(24)代入式(22),可求得E2.因此,可得如下計算公式:

    根據(jù)計算所得結(jié)果,給出如下結(jié)論:

    定理4 (1)μ2決定了Hopf分岔的方向:若μ2<0(>0),則系統(tǒng)(2)產(chǎn)生次臨界(超臨界)的Hopf分岔,且在τ<τ0(τ>τ0)時存在分岔周期解;(2)β2決定了Hopf分岔周期解的穩(wěn)定性,若β2>0(<0),則周期解是不穩(wěn)定(穩(wěn)定)的;(3)T2決定了分岔周期解的周期,若T2<0(>0),則分岔周期解的周期隨τ的增加而減少(增加)的.

    5 數(shù)值模擬

    在系統(tǒng)(2)中,取a=2.5,b0=0.1,b1=0.5,b2=2.5,β =1.2,c=0.5,d=0.6 時,通過式(25)給出的計算公式,結(jié)合推導過程,可得:

    根據(jù)定理4,模型(2)在τ0處產(chǎn)生超臨界Hopf分岔,且分岔周期解的周期是穩(wěn)定的,分岔周期解的周期將隨τ的增加而增加.進一步的數(shù)值模擬結(jié)果為:

    圖1 地方病平衡點的時間歷程圖(τ=17.5)Fig.1 time history of the endemic equilibrium when τ=17.5

    圖2 地方病平衡點的相圖(τ=17.5)Fig.2 phase diagram of the endemic equilibrium when τ=17.5

    圖3 地方病平衡點的時間歷程圖(τ=19)Fig.3 time history of the endemic equilibrium when τ=19

    圖4 地方病平衡點的相圖(τ=19)Fig.4 phase diagram of the endemic equilibrium τ=19

    圖5 不同非線性影響因素下染病個體的時間歷程圖Fig.5 time history of infected individuals under different nonlinear factors

    6 結(jié)論

    本文建立了具有非線性發(fā)生率和階段結(jié)構(gòu)的時滯SIRS模型,分析了平衡點的局部穩(wěn)定性.得出:當R0<1時,對任意的τ,無病平衡點局部漸近穩(wěn)定;當R0>1時,對任意的τ,無病平衡點是不穩(wěn)定.對地方病平衡點穩(wěn)定性分析時,在R0>1的前提下,當τ=0時,得到了地方病平衡點局部穩(wěn)定的;當τ>0時,地方病平衡點條件穩(wěn)定,且在τk處發(fā)生了Hopf分岔,即疾病會出現(xiàn)周期性爆發(fā).

    對Hopf分岔方向和分岔周期解的穩(wěn)定性及分岔周期的討論,我們運用了中心流形定理和規(guī)范型理論進行研究分析,得出系統(tǒng)(2)在參數(shù)滿足一定的條件下,發(fā)生了超臨界的Hopf分岔,分岔周期解是穩(wěn)定的,且分岔周期解的周期隨著時滯量τ的增加而增大.進一步的數(shù)值模擬結(jié)果,也驗證了理論分析的正確性.在此基礎(chǔ)上,為了驗證非線性發(fā)生率對疾病的影響效果,通過一定的數(shù)值模擬得出,適當?shù)脑黾油饨缫种茝姸?,能使染病個體的數(shù)量明顯的下降,更利于治療疾病,因此非線性發(fā)生率的引入具有非常重要的生物學意義.

    1 Kermack W O,McKendrick A G.Contribution to the mathematical theory of epidemics.Proceedings of the Royal Society of London Series A,Containing Papers of a Mathematical and Physical Character,1927,115(772):700~721.

    2 邵英英,劉孟,付新楚.無標度網(wǎng)絡(luò)上分片線性傳染力與免疫作用下的流行病動力學.動力學與控制學報,2009,7(3):264~269(Shao Y Y,Liu M,F(xiàn)u X C.Epidemic dynamics of SIR model with piecewise linear infectivity and immunizations on scale-free networks.Journal of Dynamics and Control,2009,7(3):264 ~ 269(in Chinese))

    3 Jia J W,Li Q Y.Qualitative analysis of an SIR epidemic model with stage structure.Applied Mathematics and Computation,2007,193(1):106~115

    4 Lu Z H,Gao S J.Analysis of an SI epidemic model with nonlinear transmission and stage structure.Acta Mathematica Scientia,2003,23(4):440~446

    5 Cooke K L.Stability analysis for a vector disease model.Journal of Mathematics,1979,9(1):31~42

    6 Zhang T L,Liu J L,Teng Z D.Stability of Hopf bifurcation of a delayed SIRS epidemic model with stage structure.Nonlinear Analysis:Real World Applications,2010,11(1):293~306

    7 Li X L,Wei J J.On the zeros of a fourth degree exponential polynomial with applications to a neural network model with delays.Chaos,Solitons&Fractals,2005,26(2):519~526

    8 Dieudonné J A,Dieudonné J,Mathematician F.Foundations of modern analysis.New York:Academic Press ,1960

    9 Hale J K,Lunel SMV.Introduction to functional differential equations.Springer,1993

    10 Hassard B D,Kazarinoff N D,Wan Y H.Theory of applications of Hopf bifurcation.CUP Archive,1981

    *The project supported by the National Natural Science Foundation of China(11332008)

    ? Corresponding author E-mail:hongling@mail.xjtu.edu.cn

    STABILITY AND HOPF BIFURCATION ANALYSIS OF A DELAYED SIRS EPIDEMIC MODEL WITH NONLINEAR SATURATION INCIDENCE*

    Chen Fangfang Hong Ling?
    (State Key Laboratory for Strength and Vibration of Mechanical Structures,Xi'an Jiaotong University,Xi'an710049,China)

    An SIRS epidemic model with nonlinear saturation incidence rate and time delay was investigated.By analyzing the corresponding characteristic equations,the local stability of disease-free equilibrium and endemic equilibrium was discussed.The bifurcation property was obtained as the time delay passed through a critical value.Applying the center manifold argument and normal form theory,some local bifurcation results were obtained and the formulas for determining the bifurcation direction and stability of the bifurcated periodic solution were derived.Numerical simulations were presented to illustrate the theoretical analysis.

    stability, time delay, nonlinear incidence rate, stage structure, Hopf bifurcation

    31 May 2013,

    8 June 2013.

    10.6052/1672-6553-2013-070

    2013-05-31 收到第 1 稿,2013-06-08 收到修改稿.

    *國家自然科學基金資助項目(11332008)

    E-mail:hongling@mail.xjtu.edu.cn

    猜你喜歡
    特征方程平衡點時滯
    相鄰三項線性遞推關(guān)系數(shù)列通項的簡便求法
    一些常系數(shù)非齊次線性微分方程的復數(shù)解法
    帶有時滯項的復Ginzburg-Landau方程的拉回吸引子
    探尋中國蘋果產(chǎn)業(yè)的產(chǎn)銷平衡點
    煙臺果樹(2019年1期)2019-01-28 09:34:58
    電視庭審報道,如何找到媒體監(jiān)督與司法公正的平衡點
    傳媒評論(2018年7期)2018-09-18 03:45:52
    在給專車服務(wù)正名之前最好找到Uber和出租車的平衡點
    IT時代周刊(2015年7期)2015-11-11 05:49:56
    一類n階非齊次線性微分方程特解的證明及應(yīng)用*
    一階非線性時滯微分方程正周期解的存在性
    高階齊次線性遞歸數(shù)列特征方程的由來
    考試周刊(2014年69期)2014-10-13 05:44:44
    一類時滯Duffing微分方程同宿解的存在性
    日本vs欧美在线观看视频| 亚洲成人国产一区在线观看| 亚洲av日韩精品久久久久久密| 国产av精品麻豆| www.999成人在线观看| 亚洲性夜色夜夜综合| 国产亚洲欧美在线一区二区| 国产熟女午夜一区二区三区| 国产片内射在线| 9191精品国产免费久久| ponron亚洲| 人妻一区二区av| 精品免费久久久久久久清纯 | 日韩免费高清中文字幕av| 亚洲精品一卡2卡三卡4卡5卡| 国产成人免费无遮挡视频| 午夜福利在线免费观看网站| 欧美午夜高清在线| 免费看a级黄色片| 18禁观看日本| 嫩草影视91久久| 美女视频免费永久观看网站| 麻豆乱淫一区二区| 丰满饥渴人妻一区二区三| 99国产精品一区二区三区| 9热在线视频观看99| 国产成人精品久久二区二区免费| 亚洲色图综合在线观看| 高潮久久久久久久久久久不卡| 黄色视频不卡| 日韩中文字幕欧美一区二区| www日本在线高清视频| 日日爽夜夜爽网站| 国产激情欧美一区二区| 国产精品一区二区在线不卡| 色综合欧美亚洲国产小说| 夫妻午夜视频| 亚洲情色 制服丝袜| 窝窝影院91人妻| 国产亚洲精品一区二区www | 日韩视频一区二区在线观看| 久久久国产成人精品二区 | 久久精品国产a三级三级三级| 欧美日本中文国产一区发布| tocl精华| 热re99久久国产66热| 天天躁夜夜躁狠狠躁躁| 啦啦啦 在线观看视频| 欧美色视频一区免费| 精品一区二区三区av网在线观看| 精品国产亚洲在线| 五月开心婷婷网| 波多野结衣一区麻豆| 我的亚洲天堂| 久久久久久人人人人人| 国产av精品麻豆| 久久久久久亚洲精品国产蜜桃av| 欧美不卡视频在线免费观看 | 欧美精品人与动牲交sv欧美| 9色porny在线观看| 女人被躁到高潮嗷嗷叫费观| 久久久水蜜桃国产精品网| 久9热在线精品视频| 人人妻人人添人人爽欧美一区卜| 好男人电影高清在线观看| 亚洲中文av在线| bbb黄色大片| 少妇猛男粗大的猛烈进出视频| 18禁观看日本| 国产伦人伦偷精品视频| av一本久久久久| 亚洲精品成人av观看孕妇| 国产精品国产av在线观看| 十八禁人妻一区二区| 欧美久久黑人一区二区| 999久久久精品免费观看国产| 亚洲av美国av| 精品国产一区二区久久| 成年人午夜在线观看视频| 久久狼人影院| 亚洲一区二区三区欧美精品| 日日夜夜操网爽| 久久人人97超碰香蕉20202| 成人免费观看视频高清| 亚洲人成伊人成综合网2020| 欧美日韩瑟瑟在线播放| 一区福利在线观看| 国产熟女午夜一区二区三区| 午夜日韩欧美国产| 久久狼人影院| 制服诱惑二区| 香蕉久久夜色| av视频免费观看在线观看| 99riav亚洲国产免费| 久久影院123| 不卡一级毛片| 中文字幕人妻丝袜制服| 中出人妻视频一区二区| 午夜两性在线视频| 国产一区二区三区在线臀色熟女 | 香蕉久久夜色| 热99久久久久精品小说推荐| 国产极品粉嫩免费观看在线| 又黄又爽又免费观看的视频| 精品人妻在线不人妻| 十八禁人妻一区二区| xxxhd国产人妻xxx| 国产深夜福利视频在线观看| 一级片免费观看大全| 亚洲精品国产色婷婷电影| 激情视频va一区二区三区| 国产精品一区二区免费欧美| 天堂俺去俺来也www色官网| 国产精品二区激情视频| e午夜精品久久久久久久| 亚洲自偷自拍图片 自拍| 国产又色又爽无遮挡免费看| 免费看十八禁软件| 色播在线永久视频| 日韩精品免费视频一区二区三区| 啪啪无遮挡十八禁网站| 黄色a级毛片大全视频| 男人舔女人的私密视频| av网站免费在线观看视频| 精品第一国产精品| 午夜精品国产一区二区电影| 99国产综合亚洲精品| 人妻久久中文字幕网| 精品少妇一区二区三区视频日本电影| 黄色视频不卡| 777米奇影视久久| 国产精品98久久久久久宅男小说| 久久久水蜜桃国产精品网| aaaaa片日本免费| 岛国在线观看网站| 久久久水蜜桃国产精品网| 亚洲午夜精品一区,二区,三区| tube8黄色片| 一级a爱视频在线免费观看| 侵犯人妻中文字幕一二三四区| 午夜福利乱码中文字幕| 成在线人永久免费视频| 欧美日韩福利视频一区二区| av天堂久久9| 欧美激情久久久久久爽电影 | e午夜精品久久久久久久| 香蕉久久夜色| 天堂俺去俺来也www色官网| 在线永久观看黄色视频| 欧美日本中文国产一区发布| 中亚洲国语对白在线视频| 在线观看免费午夜福利视频| 日韩欧美在线二视频 | 在线永久观看黄色视频| 精品国产乱码久久久久久男人| 久久久久国内视频| 国产成人av教育| 91麻豆av在线| 免费一级毛片在线播放高清视频 | 香蕉国产在线看| 国产成人av教育| 50天的宝宝边吃奶边哭怎么回事| 丝袜美足系列| 在线观看www视频免费| 热99久久久久精品小说推荐| 大陆偷拍与自拍| 中文字幕人妻熟女乱码| 亚洲第一欧美日韩一区二区三区| 亚洲一码二码三码区别大吗| 国产欧美日韩一区二区精品| 精品少妇久久久久久888优播| 亚洲欧美一区二区三区久久| 韩国精品一区二区三区| 一个人免费在线观看的高清视频| 18禁黄网站禁片午夜丰满| 成人三级做爰电影| 国产精品久久久久成人av| 精品久久久久久,| 老司机影院毛片| 亚洲精品中文字幕一二三四区| 免费不卡黄色视频| 法律面前人人平等表现在哪些方面| 成年人午夜在线观看视频| 国产淫语在线视频| 最近最新免费中文字幕在线| 久久人妻av系列| 在线观看www视频免费| 国产在线观看jvid| 亚洲黑人精品在线| 午夜精品在线福利| 熟女少妇亚洲综合色aaa.| 久久精品亚洲精品国产色婷小说| 国产xxxxx性猛交| 岛国毛片在线播放| 无人区码免费观看不卡| 久久人妻av系列| 不卡一级毛片| 久久中文看片网| 亚洲av日韩在线播放| 精品国产一区二区久久| 中文字幕最新亚洲高清| 日韩免费高清中文字幕av| 国产欧美日韩一区二区三区在线| 91麻豆精品激情在线观看国产 | 99热只有精品国产| 老汉色av国产亚洲站长工具| 亚洲精华国产精华精| 国产精品亚洲av一区麻豆| 国产精品一区二区免费欧美| 亚洲午夜精品一区,二区,三区| 婷婷精品国产亚洲av在线 | 黄色丝袜av网址大全| 大型黄色视频在线免费观看| 超色免费av| 女人爽到高潮嗷嗷叫在线视频| 欧美老熟妇乱子伦牲交| 欧美丝袜亚洲另类 | 午夜福利影视在线免费观看| 国产一区二区三区综合在线观看| 国产精品免费大片| 村上凉子中文字幕在线| 操出白浆在线播放| 黄网站色视频无遮挡免费观看| 亚洲九九香蕉| 少妇猛男粗大的猛烈进出视频| av视频免费观看在线观看| 国产麻豆69| 又紧又爽又黄一区二区| 亚洲七黄色美女视频| 人妻丰满熟妇av一区二区三区 | 黑人巨大精品欧美一区二区蜜桃| 午夜精品在线福利| 水蜜桃什么品种好| 90打野战视频偷拍视频| 欧美精品av麻豆av| 美女午夜性视频免费| 欧美色视频一区免费| 手机成人av网站| 欧美日韩中文字幕国产精品一区二区三区 | 手机成人av网站| 女人久久www免费人成看片| 欧美人与性动交α欧美精品济南到| 国产激情欧美一区二区| 热99久久久久精品小说推荐| tocl精华| 成年人午夜在线观看视频| 女性生殖器流出的白浆| 国产成人av激情在线播放| 日本五十路高清| 91国产中文字幕| 两性午夜刺激爽爽歪歪视频在线观看 | 黄色视频不卡| 成人18禁在线播放| 亚洲第一av免费看| 在线av久久热| 老熟女久久久| 嫁个100分男人电影在线观看| 精品一区二区三区av网在线观看| 国精品久久久久久国模美| 国产亚洲精品久久久久5区| 18禁黄网站禁片午夜丰满| 精品久久久久久久久久免费视频 | 99国产精品一区二区蜜桃av | 精品一区二区三卡| 黄色视频,在线免费观看| x7x7x7水蜜桃| 成人18禁高潮啪啪吃奶动态图| 欧美 亚洲 国产 日韩一| 精品亚洲成a人片在线观看| 啦啦啦 在线观看视频| 亚洲精华国产精华精| 精品第一国产精品| 国产精品综合久久久久久久免费 | 欧美精品一区二区免费开放| 精品国产一区二区三区四区第35| 99re6热这里在线精品视频| 天堂俺去俺来也www色官网| www.999成人在线观看| av福利片在线| 丝袜美足系列| 婷婷精品国产亚洲av在线 | 欧洲精品卡2卡3卡4卡5卡区| 亚洲专区国产一区二区| 亚洲精品久久成人aⅴ小说| 中文字幕人妻丝袜制服| 日本撒尿小便嘘嘘汇集6| 国产成+人综合+亚洲专区| 涩涩av久久男人的天堂| 老司机影院毛片| 中文字幕色久视频| 窝窝影院91人妻| 国产一区二区三区视频了| 亚洲熟妇熟女久久| x7x7x7水蜜桃| 岛国在线观看网站| 欧美黑人精品巨大| 欧美丝袜亚洲另类 | 国产精品99久久99久久久不卡| 亚洲九九香蕉| 美国免费a级毛片| 女人久久www免费人成看片| 亚洲在线自拍视频| 人妻一区二区av| 丰满人妻熟妇乱又伦精品不卡| 搡老岳熟女国产| 99久久精品国产亚洲精品| 黄色 视频免费看| 视频在线观看一区二区三区| 在线观看66精品国产| 大型av网站在线播放| 咕卡用的链子| 国产男女超爽视频在线观看| netflix在线观看网站| 啦啦啦 在线观看视频| 国产在线精品亚洲第一网站| 人人妻人人爽人人添夜夜欢视频| 久久热在线av| 日韩欧美三级三区| 淫妇啪啪啪对白视频| 亚洲精品一二三| 法律面前人人平等表现在哪些方面| 国产精品久久久av美女十八| 久久久久视频综合| 亚洲成人免费av在线播放| 无人区码免费观看不卡| 99久久人妻综合| 久99久视频精品免费| 国产精品久久久久久精品古装| 1024视频免费在线观看| 一边摸一边做爽爽视频免费| 又大又爽又粗| 在线观看免费视频日本深夜| 18在线观看网站| 亚洲午夜精品一区,二区,三区| 男女之事视频高清在线观看| 女性被躁到高潮视频| 黄频高清免费视频| 女性被躁到高潮视频| 性色av乱码一区二区三区2| 天天躁日日躁夜夜躁夜夜| 大香蕉久久网| 色播在线永久视频| 久久中文字幕人妻熟女| 亚洲专区字幕在线| 午夜亚洲福利在线播放| 在线十欧美十亚洲十日本专区| 中文字幕人妻丝袜制服| 免费日韩欧美在线观看| av网站免费在线观看视频| 岛国在线观看网站| av超薄肉色丝袜交足视频| 色老头精品视频在线观看| а√天堂www在线а√下载 | 久9热在线精品视频| 国产无遮挡羞羞视频在线观看| 18在线观看网站| 黑人巨大精品欧美一区二区mp4| 亚洲精华国产精华精| 777久久人妻少妇嫩草av网站| 亚洲欧美一区二区三区黑人| 欧美精品啪啪一区二区三区| 电影成人av| 水蜜桃什么品种好| 日韩大码丰满熟妇| 亚洲精品国产区一区二| 久久香蕉激情| 日韩制服丝袜自拍偷拍| 两人在一起打扑克的视频| 乱人伦中国视频| 精品久久久久久电影网| 免费黄频网站在线观看国产| 成熟少妇高潮喷水视频| 国产av一区二区精品久久| 国产精品 欧美亚洲| 久久青草综合色| 18禁观看日本| 久久亚洲精品不卡| 精品一品国产午夜福利视频| 久久天躁狠狠躁夜夜2o2o| 色播在线永久视频| 亚洲精品一二三| 中出人妻视频一区二区| 好男人电影高清在线观看| av网站免费在线观看视频| 日日爽夜夜爽网站| 亚洲成av片中文字幕在线观看| 中出人妻视频一区二区| 亚洲全国av大片| a级毛片在线看网站| 亚洲成人手机| 在线永久观看黄色视频| 十分钟在线观看高清视频www| 女人被躁到高潮嗷嗷叫费观| 黄色 视频免费看| 日韩成人在线观看一区二区三区| 亚洲国产欧美一区二区综合| 又黄又粗又硬又大视频| 国产精品九九99| 久久影院123| 久久久久久免费高清国产稀缺| 久久久久久久精品吃奶| 午夜日韩欧美国产| 国产一区二区三区综合在线观看| 国产精品 欧美亚洲| 在线观看一区二区三区激情| 久久久国产成人免费| 国产亚洲精品一区二区www | 少妇的丰满在线观看| 12—13女人毛片做爰片一| 麻豆av在线久日| 色婷婷久久久亚洲欧美| 久久国产精品大桥未久av| xxxhd国产人妻xxx| 一进一出好大好爽视频| 久久久国产一区二区| av中文乱码字幕在线| 欧美日韩一级在线毛片| 国产亚洲一区二区精品| 大型av网站在线播放| 美国免费a级毛片| 久久久久久久午夜电影 | 韩国av一区二区三区四区| 十八禁网站免费在线| av天堂久久9| 欧美黑人精品巨大| 可以免费在线观看a视频的电影网站| 丝瓜视频免费看黄片| 欧美人与性动交α欧美精品济南到| 亚洲专区国产一区二区| 男女免费视频国产| 亚洲国产欧美网| 大香蕉久久网| 国产精品一区二区免费欧美| 两个人免费观看高清视频| 精品久久蜜臀av无| ponron亚洲| 成人黄色视频免费在线看| 亚洲三区欧美一区| 亚洲国产欧美一区二区综合| 国产成+人综合+亚洲专区| 国产男女超爽视频在线观看| 国产1区2区3区精品| 深夜精品福利| 丰满人妻熟妇乱又伦精品不卡| 日韩 欧美 亚洲 中文字幕| 69精品国产乱码久久久| 久久香蕉精品热| av超薄肉色丝袜交足视频| 免费日韩欧美在线观看| 露出奶头的视频| 在线观看免费高清a一片| 欧美激情极品国产一区二区三区| 国产精品一区二区免费欧美| 国产又色又爽无遮挡免费看| 首页视频小说图片口味搜索| 国产精品综合久久久久久久免费 | 在线观看66精品国产| 久久久久久久精品吃奶| 好看av亚洲va欧美ⅴa在| 亚洲av电影在线进入| 人妻丰满熟妇av一区二区三区 | 欧美丝袜亚洲另类 | 亚洲成人免费av在线播放| 无人区码免费观看不卡| 12—13女人毛片做爰片一| 久久九九热精品免费| 午夜两性在线视频| 麻豆av在线久日| 国产精品免费视频内射| 69av精品久久久久久| 岛国在线观看网站| 狠狠婷婷综合久久久久久88av| 天天躁日日躁夜夜躁夜夜| 国产精品一区二区在线观看99| 一级毛片女人18水好多| 国产91精品成人一区二区三区| 国内久久婷婷六月综合欲色啪| 精品一区二区三卡| 中出人妻视频一区二区| 亚洲欧美色中文字幕在线| 国产成人欧美在线观看 | 久久久久视频综合| 亚洲美女黄片视频| 国产成人免费观看mmmm| 免费看a级黄色片| 女人精品久久久久毛片| 国产野战对白在线观看| 极品人妻少妇av视频| 精品人妻1区二区| 99久久国产精品久久久| 自线自在国产av| 欧美日韩一级在线毛片| 婷婷成人精品国产| 国产av又大| 热re99久久精品国产66热6| 男人的好看免费观看在线视频 | 精品卡一卡二卡四卡免费| av有码第一页| 91精品国产国语对白视频| 美女高潮到喷水免费观看| 在线观看免费高清a一片| 亚洲av日韩精品久久久久久密| 久久国产亚洲av麻豆专区| 日韩三级视频一区二区三区| 久久人人爽av亚洲精品天堂| 免费黄频网站在线观看国产| 精品久久久久久电影网| 91字幕亚洲| 国产成人av教育| 亚洲精品在线观看二区| 精品久久久久久电影网| videos熟女内射| 可以免费在线观看a视频的电影网站| 在线视频色国产色| 亚洲第一欧美日韩一区二区三区| 丁香欧美五月| 自拍欧美九色日韩亚洲蝌蚪91| 国产精品久久久av美女十八| 成年版毛片免费区| 人成视频在线观看免费观看| 精品福利观看| 一级,二级,三级黄色视频| 精品少妇一区二区三区视频日本电影| a在线观看视频网站| 老汉色∧v一级毛片| 一级作爱视频免费观看| 亚洲全国av大片| av片东京热男人的天堂| 精品久久久久久,| 老司机福利观看| 国产蜜桃级精品一区二区三区 | 18禁国产床啪视频网站| 91老司机精品| 亚洲av成人av| 在线国产一区二区在线| 91av网站免费观看| netflix在线观看网站| 99re6热这里在线精品视频| 成人影院久久| 午夜福利,免费看| 亚洲国产中文字幕在线视频| 精品亚洲成a人片在线观看| 国产成人欧美在线观看 | 欧美国产精品一级二级三级| 国产免费av片在线观看野外av| 欧美久久黑人一区二区| 亚洲熟妇熟女久久| 黑人操中国人逼视频| 亚洲精品中文字幕一二三四区| 老司机在亚洲福利影院| 最新的欧美精品一区二区| 亚洲精品国产区一区二| 视频区图区小说| 久久久国产欧美日韩av| 成人三级做爰电影| 在线永久观看黄色视频| 91老司机精品| 久久久国产精品麻豆| 好看av亚洲va欧美ⅴa在| 高清黄色对白视频在线免费看| 老司机深夜福利视频在线观看| 天天操日日干夜夜撸| 精品国产超薄肉色丝袜足j| 国产精品偷伦视频观看了| 黄网站色视频无遮挡免费观看| 欧美人与性动交α欧美精品济南到| 韩国精品一区二区三区| 久久久久视频综合| 久久精品aⅴ一区二区三区四区| www.精华液| 亚洲精品一二三| 黄色怎么调成土黄色| 国产成人精品在线电影| 亚洲专区字幕在线| 国产精品久久久久久人妻精品电影| 国产精品亚洲一级av第二区| 午夜精品久久久久久毛片777| 国产在线一区二区三区精| 欧美成人免费av一区二区三区 | 视频在线观看一区二区三区| 身体一侧抽搐| 日韩一卡2卡3卡4卡2021年| 99re在线观看精品视频| 成人18禁在线播放| 99国产精品免费福利视频| 国产精品电影一区二区三区 | 91九色精品人成在线观看| 亚洲av日韩在线播放| 黑人巨大精品欧美一区二区mp4| 中文欧美无线码| 午夜精品国产一区二区电影| 亚洲成国产人片在线观看| 80岁老熟妇乱子伦牲交| 人人妻人人添人人爽欧美一区卜| 亚洲av电影在线进入| 成年人黄色毛片网站| 久久久久久人人人人人| 国产亚洲欧美98| 久久久国产欧美日韩av| svipshipincom国产片| 亚洲男人天堂网一区| 在线国产一区二区在线| 亚洲一码二码三码区别大吗| 在线观看免费视频日本深夜| 国产欧美日韩综合在线一区二区| 久久国产乱子伦精品免费另类| 国产精品免费大片| 国产精品二区激情视频| 在线观看舔阴道视频| 高清黄色对白视频在线免费看| 丰满饥渴人妻一区二区三| 日韩一卡2卡3卡4卡2021年| 日韩欧美一区视频在线观看| 看免费av毛片| 亚洲人成77777在线视频|