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

    一個類Lorenz系統(tǒng)的Hopf分岔分析及分岔控制

    2015-12-19 09:15:56張中華付景超鄧冠男
    關(guān)鍵詞:波形圖平衡點(diǎn)線性

    張中華,付景超,鄧冠男

    (東北電力大學(xué)理學(xué)院,吉林 吉林132012)

    0 引言

    1963年,美國著名氣象學(xué)家E.N.Lorenz在刻畫熱對流不穩(wěn)定性時發(fā)現(xiàn)了第一個混沌系統(tǒng),即Lorenz系統(tǒng)[1],其動力學(xué)方程表示為

    其中,x為對流強(qiáng)度,y為上升流和下降流的溫差,z為鉛直方向溫度分布的非線性強(qiáng)度。c為Rayleigh數(shù),為系統(tǒng)的主要控制參數(shù),a是Prandt數(shù),b是外形比。當(dāng)參數(shù)時,上述系統(tǒng)有混沌吸引子[2-3]。

    自Lorenz發(fā)現(xiàn)第一個混沌系統(tǒng)之后,人們開始不斷尋找新的混沌系統(tǒng),并發(fā)現(xiàn)了許多新的混沌系統(tǒng),如Chua系統(tǒng)、Chen系統(tǒng)、Qi系統(tǒng)等,對這方面的研究也很多[4-7]。1994年,Sprott利用計(jì)算搜索的方法提出了19種簡單形式的三階二次混沌系統(tǒng)[8],但對其中的系統(tǒng)進(jìn)行分析的文獻(xiàn)不多。文獻(xiàn)[9]介紹了三維二次多項(xiàng)式自治混沌系統(tǒng)中不滿足Shil′nikov定理?xiàng)l件的一種特殊系統(tǒng),即著名的Sprott C系統(tǒng):并在此基礎(chǔ)上提出了一個新的混沌系統(tǒng)(稱廣義Sprott系統(tǒng)):進(jìn)一步研究了該系統(tǒng)的Hopf分岔及周期解穩(wěn)定性。2008年,楊啟貴和陳關(guān)榮提出了一個新的三維混沌系統(tǒng)[10]文獻(xiàn)[11]在此基礎(chǔ)上構(gòu)造了一個新的四維自治超混沌系統(tǒng),研究了系統(tǒng)的基本動力學(xué)行為,并設(shè)計(jì)自適應(yīng)追蹤控制器,實(shí)現(xiàn)對各種不同參考信號的單變量追蹤控制。文獻(xiàn)[12]研究了新的三維自治連續(xù)時間類Lorenz系統(tǒng)平衡點(diǎn)的存在性與穩(wěn)定性,并給出了系統(tǒng)發(fā)生Hopf分岔和余維二退化Hopf分岔的參數(shù)條件。但上述研究主要集中在混沌的同步和控制方面,對分岔控制的研究相對較少。

    本文針對文獻(xiàn)[12]中的模型,主要研究該模型的Hopf分岔行為及分岔控制問題。

    1 新的類Lorenz系統(tǒng)模型

    考慮模型:

    其中,x,y,z為狀態(tài)變量,a,b,c為系統(tǒng)參數(shù)。

    1.1 平衡點(diǎn)穩(wěn)定性分析

    引理 由文獻(xiàn)[12]中的結(jié)論知,當(dāng)a≠0,b≠0時,若bc<0,系統(tǒng)(1)只有一個平衡點(diǎn)O(0,0,0);當(dāng)bc≥0,系統(tǒng)(1)有3個平衡點(diǎn)且當(dāng)a>0,c<0時,平衡點(diǎn)O漸近穩(wěn)定;b≠0,ac>0時,平衡點(diǎn)O不穩(wěn)定;當(dāng)bc>0,a>0,a+b>0,ab(a+b-2c)>0時,平衡點(diǎn)E1,E2漸近穩(wěn)定;當(dāng)bc>0,ab時,系統(tǒng)在平衡點(diǎn)E1發(fā)生Hopf分岔,且

    取c為分岔參數(shù),由上面的結(jié)論知,當(dāng)a=7,b=2,c=c0=4.5時,系統(tǒng)(1)在平衡點(diǎn)E1(3,3,4.5)和E2(-3,-3,4.5)處發(fā)生Hopf分岔,分岔位置如圖1所示,分岔圖如圖2~圖3所示。另外,從圖2和圖3可看出,系統(tǒng)(1)在E1和E2處的Hopf分岔方向均發(fā)生在臨界點(diǎn)右側(cè)。

    1.2 Hopf分岔類型及周期解

    首先,討論系統(tǒng)在平衡點(diǎn)E1處的分岔類型。經(jīng)計(jì)算得α′(0)=Re{λ′(c0)}=0.147 4>0,ω′(0)=lm{λ′(c0)}=0.354 5。做變換

    將平衡點(diǎn)E1移到原點(diǎn)。整理得

    圖1 平衡點(diǎn)分岔位置圖Fig.1 Bifurcation position diagram near equilibrium

    圖2 系統(tǒng)(1)在E1處的分岔圖Fig.2 Bifurcation diagram of system (1)for E1

    圖3 系統(tǒng)(1)在E2處的分岔圖Fig.3 Bifurcation diagram of system (1)for E2

    則系統(tǒng)(2)的線性化矩陣為

    矩陣P有一對純虛特征根λ1,2=±3.741 7i和一個負(fù)實(shí)根λ3=-9,設(shè)α1,α3分別為特征值λ1=+3.741 7i和λ3=-9對應(yīng)的特征向量。通過計(jì)算得

    定義矩陣

    作變換

    系統(tǒng)(2)變?yōu)?/p>

    其中

    由文獻(xiàn)[13]及上式計(jì)算得

    將上述結(jié)果代入式(5)

    得到

    又因?yàn)棣痢洌?)=0.147 4>0,由Hopf分岔理論知,系統(tǒng)(1)在分岔點(diǎn)E1處的分岔周期解是軌道漸近穩(wěn)定的,產(chǎn)生與初始條件無關(guān)的穩(wěn)定極限環(huán),發(fā)生的Hopf分岔為超臨界分岔,分岔方向?yàn)閏>c0=4.5,這點(diǎn)從圖2也可以看出。當(dāng)c<c0=4.5且時,平衡點(diǎn)漸近穩(wěn)定,如圖4所示;當(dāng)c>c0=4.5且時,平衡點(diǎn)不穩(wěn)定并在其附近產(chǎn)生穩(wěn)定的極限環(huán),如圖5所示。其中,仿真初值為(x0,y0,z0)=(4.5,4.5,4.5)。

    根據(jù)上面的計(jì)算結(jié)果,可以計(jì)算出重要特征向量

    進(jìn)一步得到系統(tǒng)的分岔周期解為

    其中

    周期和特征性指數(shù)分別為

    用同樣的方法可判斷系統(tǒng)(1)在平衡點(diǎn)E2處的分岔類型為超臨界Hopf分岔,系統(tǒng)受擾離開平衡點(diǎn)將產(chǎn)生等幅振蕩,在平衡點(diǎn)E2附近產(chǎn)生一個穩(wěn)定的極限環(huán),如圖6所示。

    圖4 c=4.2時,系統(tǒng)(1)在E1處的相圖Fig.4 Phase chart of system(1)at E1for c=4.2

    圖5 c=4.6時,系統(tǒng)(1)的相圖Fig.5 Phase chart of system(1)at E1for c=4.6

    圖6 c=4.6時,系統(tǒng)(1)在E2處的相圖Fig.6 Phase chart of system (1)at E2for c=4.6

    2 系統(tǒng)Hopf分岔控制

    2.1 線性控制

    不改變原系統(tǒng)的Hopf分岔點(diǎn),采用線性控制器,對系統(tǒng)(1)作如式(8)的控制:

    其中,k1,k2,k3為待定控制參數(shù)。當(dāng)a=7,b=2,c=4.5時,(x0,y0,z0)==(3,3,4.5)。此時系統(tǒng)(8)變?yōu)?/p>

    則系統(tǒng)(9)轉(zhuǎn)化為

    系統(tǒng)(10)在平衡點(diǎn)(0,0,0)的穩(wěn)定性情況相當(dāng)于系統(tǒng)(9)在平衡點(diǎn)(3,3,4.5)的穩(wěn)定性情況。系統(tǒng)(10)在平衡點(diǎn)(0,0,0)的線性化矩陣為

    下面分3種情況討論線性控制對Hopf分岔位置的影響。為對結(jié)果進(jìn)行比較,以下的仿真初值均選為(15,3.5,7.5)。

    2.1.1 k2≠0,k1=0,k3=0時,受控系統(tǒng) Hopf分岔分析

    此時矩陣A對應(yīng)特征方程為

    根據(jù)Routh-Hurwitz原理,方程所有的根具有負(fù)實(shí)部的充要條件是9-k2>0,-14k2-126>0,(9-k2)(14-9k2)>-14k2-126。因此,當(dāng)且僅當(dāng)k2<-9(k1=0,k3=0)時,方程所有根的實(shí)部為負(fù)數(shù),在這個條件下,系統(tǒng)(10)的平衡點(diǎn)(0,0,0)是漸近穩(wěn)定的。如取k2=-10,k1=0,k3=0時,矩陣A對應(yīng)特征根為λ1=-0.138 1,λ2,3=-9.431 0±3.528 8i。所有特征根具有負(fù)實(shí)部,系統(tǒng)(10)在平衡點(diǎn)(0,0,0)漸近穩(wěn)定,進(jìn)而得到系統(tǒng)(9)在平衡點(diǎn)(3,3,4.5)漸近穩(wěn)定,故系統(tǒng)(1)在施加線性控制器后,Hopf分岔行為得到控制(見圖7)。這表明可以采用只在系統(tǒng)(8)的第2個方程中施加線性控制器的方法來抑制Hopf分岔的出現(xiàn)。但有時為了表明不同線性控制器參數(shù)對Hopf分岔位置影響的不同以及它們之間的相互影響,還可以在系統(tǒng)的其他方程中施加線性控制器,討論施加單個線性控制器和施加多個線性控制器對抑制Hopf分岔發(fā)生的不同影響。

    2.1.2 k2=-10,k1=0,k3≠0時,受控系統(tǒng) Hopf分岔分析

    此時矩陣A對應(yīng)特征方程為

    根據(jù)Routh-Hurwitz原理,當(dāng)且僅當(dāng)k3<0.2(k2=-10,k1=0)時,方程所有根具有負(fù)實(shí)部,此時,系統(tǒng)(10)的平衡點(diǎn)(0,0,0)漸近穩(wěn)定。如取k2=-10,k1=0,k3=0.1,矩陣A 對應(yīng)特征根為λ1=-0.069 3,λ2,3=-9.415 3±3.513 7i;如取k2=-10,k1=0,k3=-8,矩陣A 對應(yīng)特征根為λ1=-3.762 0,λ2,3=-11.619 0±4.192 6i。所有特征根具有負(fù)實(shí)部。所以,系統(tǒng)(10)在平衡點(diǎn)(0,0,0)漸近穩(wěn)定,進(jìn)一步說明原系統(tǒng)Hopf分岔行為得到控制(見圖8~9)。

    圖7 k2=-10,k1=0,k3=0時,系統(tǒng)(8)波形圖Fig.7 Waveform chart of system (8)for k2=-10,k1=0,k3=0

    圖8 k2=-10,k1=0,k3=0.1時,系統(tǒng)(8)波形圖Fig.8 Waveform chart of system (8)for k2=-10,k1=0,k3=0.1

    圖9 k2=-10,k1=0,k3=-8時,系統(tǒng)(8)波形圖Fig.9 Waveform chart of system (8)for k2=-10,k1=0,k3=-8

    由圖8和圖9可知,當(dāng)k1=0,k2=-10時,系統(tǒng)Hopf分岔點(diǎn)趨于穩(wěn)定的速度隨著k3的減小而加快。由圖7和圖8知,施加雙線性反饋控制比單線性反饋控制的效果好,系統(tǒng)分岔點(diǎn)趨于穩(wěn)定的速度更快。

    2.1.3 k2=-10,k1≠0,k3=-8時,受控系統(tǒng) Hopf分岔分析

    此時矩陣A對應(yīng)特征方程為

    方程所有的根具有負(fù)實(shí)部的充要條件是k1<5.74(k2=-10,k3=-8),在此條件下,系統(tǒng)(10)的平衡點(diǎn)(0,0,0)漸近穩(wěn)定。如取k2=-10,k1=3,k3=-8,矩陣A對應(yīng)特征根為λ1=-2.020 9,λ2,3=-10.989 5±3.848 6i;如取k2=-10,k1=-9,k3=-8時,矩陣A對應(yīng)特征根為λ1=-6.380 8,λ2,3=-14.809 6±3.417 9i。所有特征根具有負(fù)實(shí)部,系統(tǒng)(10)在平衡點(diǎn)(0,0,0)漸近穩(wěn)定(見圖10和圖11)。從圖10和圖11可看出,當(dāng)k2,k3固定時,系統(tǒng)Hopf分岔點(diǎn)趨于穩(wěn)定的速度隨著k1的減小而加快。

    綜上可知,當(dāng)線性控制參數(shù)滿足一定條件時,可抑制原系統(tǒng)Hopf分岔的發(fā)生,進(jìn)而抑制混沌的出現(xiàn)。

    2.2 極限環(huán)幅值控制

    2.2.1 基于Normal Form方法的立方非線性反饋控制

    不改變原系統(tǒng)的Hopf分岔點(diǎn),對系統(tǒng)(1)施加如下非線性控制器:

    因?yàn)橄到y(tǒng)在平衡點(diǎn)的穩(wěn)定性主要由系統(tǒng)在平衡點(diǎn)的線性化矩陣對應(yīng)特征方程的特征根性質(zhì)來判定,又因?yàn)槲闹杏懻摰钠胶恻c(diǎn)是O(0,0,0),所以,影響系統(tǒng)在平衡點(diǎn)O的穩(wěn)定性因素是系統(tǒng)的線性項(xiàng)部分,非線性項(xiàng)影響的是系統(tǒng)的非線性特性,比如周期解的穩(wěn)定性情況及周期解幅值大小等。下面討論非線性控制器參數(shù)對極限環(huán)幅值的影響。

    則系統(tǒng)(14)化為

    其中,系統(tǒng)(15)的線性和非線性部分分別為

    矩陣B有一對純虛特征根λ1,2=±3.741 7i,相應(yīng)的特征向量為φ,。其中

    根據(jù)文獻(xiàn)[14]中提出的計(jì)算Normal Form的直接方法,引入變換:

    即可求得方程的Hopf分岔Normal Form

    其中,C = 〈ψ,F(xiàn)21〉,ψ 滿足(BT-3.7417iI)ψ =0,〈ψ,φ〉=1,

    將非線性變換(17)代入非線性項(xiàng)中整理成關(guān)于u,uˉ的多項(xiàng)式形式

    其中,F(xiàn)21為u2對應(yīng)的系數(shù)向量。

    上式中,Hjk,m為向量Hjk中的第m 個元素。其中

    由式(20)~ (23)計(jì)算C得

    所以,式(14)的Normal Form為

    2.2.2 仿真研究

    由 Hopf分岔理論知,當(dāng)k<0.243 4時,Re(C)<0,又因?yàn)棣痢洌?)=0.147 4>0,所以,系統(tǒng)(14)的 Hopf分岔類型為超臨界分岔且Hopf分岔方向?yàn)閏>c0=4.5,極限環(huán)穩(wěn)定,極限環(huán)的幅控關(guān)系為

    如圖12所示,極限環(huán)幅值r與控制系數(shù)k是增函數(shù)關(guān)系。另外從圖13和圖14也可看出,振動幅值(極限環(huán)幅值)隨著k的減小而減小。并且,k的范圍恰好包含了k=0的情況,與原系統(tǒng)的分岔類型相對應(yīng)。其中,仿真初值為(x0,y0,z0)=(4.5,4.5,4.5)。

    當(dāng)Re(C)>0,即k>0.243 4時,系統(tǒng)(14)的Hopf分岔類型為亞臨界分岔,極限環(huán)不穩(wěn)定(見圖15)。

    圖10 k2=-10,k1=3,k3=-8時,系統(tǒng)(8)波形圖Fig.10 Waveform chart of system (8)for k2=-10,k1=3,k3=-8

    圖11 k2=-10,k1=-9,k3=-8時,系統(tǒng)(8)波形圖Fig.11 Waveform chart of system (8)for k2=-10,k1=-9,k3=-8

    圖12 受控系統(tǒng)(14)的限極環(huán)幅值曲線(c=4.6)Fig.12 Gain amplitude curves for the controlled system (14)(with c=4.6)

    圖14 c=4.6,k=-5時,系統(tǒng)(14)波形圖Fig.14 Phase chart of system (14)with c=4.6,k=-5

    圖15 c=4.6,k=1.2時,系統(tǒng)(14)相圖Fig.15 Phase chart of system (14)for c=4.6,k=1.2

    由以上分析知,當(dāng)非線性控制參數(shù)滿足一定條件時,可以改變原系統(tǒng)的分岔類型和分岔周期解振幅大小,從而實(shí)現(xiàn)類Lorenz系統(tǒng)的極限環(huán)幅值控制。

    3 結(jié)論

    本文主要研究了一個新的類Lorenz系統(tǒng)的Hopf分岔行為和分岔控制問題。首先,通過規(guī)范形方法計(jì)算出系統(tǒng)的穩(wěn)定性指標(biāo),進(jìn)而判定系統(tǒng)的Hopf分岔類型;然后,對系統(tǒng)分別施加線性和非線性控制器,并詳細(xì)討論了線性控制器和非線性控制器對Hopf分岔位置、分岔類型和分岔周期解幅值(極限環(huán)幅值)的影響,進(jìn)而實(shí)現(xiàn)類Lorenz系統(tǒng)的分岔控制問題。

    [1] Lorenz E N.Deterministic non-periodic flows[J].J Atmos Sci,1963,20:130-141.

    [2] Tucker W.The Lorenz attractor exists[J].C R Acad Sci Pairis,1999,328:1197-1202.

    [3]Stewart I.The Lorenz attractor exists[J].Nature,2002,406:948-949.

    [4] Pang S Q,Liu Y J.A new hyperchaotic system from the Lüsystem and its control[J].Journal of Computational and Applied Mathematics,2011,235(8):2775-2789.

    [5] Mahmoud E E.Dynamics and synchronization of new hyperchaotic complex Lorenz system[J].Mathematical and Computer Modelling,2012,55(7/8):1951-1962.

    [6] Niu Y J,Wang X Y.A new hyperchaotic system and its circuit implementation[J].Commun Nonlinear Sci Numer Simulat,2010,15(11):3518-3524.

    [7] Mkaouar H,Boubaker O.Chaos synchronization for master slave piecewise linear systems:application to Chua′s circuit[J].Commun Nonlinear Sci Numer Simulat,2012,17(3):1292-1302.

    [8]Sprott J C.Some simple chaotic flows[J].Physical Review E,1994,50(2):647-650.

    [9] 魏周超,楊啟貴.基于Lorenz型系統(tǒng)族的三維系統(tǒng)的復(fù)雜動力學(xué)研究[D].廣州:華南理工大學(xué),2011.Wei Zhouchao,Yang Qigui.Research on complicated dynamics of three-dimensional systems based on a family of Lorenz-type systems[D].Guangzhou:South China University of Thecnology,2011.

    [10]Yang Q G,Chen G R.A chaotic system with one saddle and two stable node-foci[J].Int J Bifurc Chaos,2008,18(5):1393-1414.

    [11]李春來,禹思敏.一個新的超混沌系統(tǒng)及其自適應(yīng)追蹤控制[J].物理學(xué)報(bào),2012,61(4):1-7.Li Chunlai,Yu Simin.A new hyperchaotic system and its adaptive tracking control[J].Acta Phys Sin,2012,61(4):1-7.

    [12]李群宏,徐德貴.一個類Lorenz系統(tǒng)的動力學(xué)分析[J].重慶理工大學(xué)學(xué)報(bào),2011,25(2):112-116.Li Qunhong,Xu Degui.Dynamics analysis of a Lorenz-like system[J].Journal of Chongqing University of Technology,2011,25(2):112-116.

    [13]Hassard B D.Theory and Application of Hopf Bifurcation[M].New York:Cambridge University,1981.

    [14]吳志強(qiáng),孫立明.基于washout濾波器的R¨ossler系統(tǒng) Hopf分岔控制[J].物理學(xué)報(bào),2011,60(5):1-5.Wu Zhiqiang,Sun Liming.Hopf bifurcation control of the R¨ossler system based on washout filter controller[J].Acta Physical Sinica,2011,60(5):1-5.

    猜你喜歡
    波形圖平衡點(diǎn)線性
    交流接觸器典型通斷過程分析
    漸近線性Klein-Gordon-Maxwell系統(tǒng)正解的存在性
    由波形圖領(lǐng)悟聲音特性
    淺析CH-γ方程中解的求法
    線性回歸方程的求解與應(yīng)用
    用DIS數(shù)字化信息系統(tǒng)測量金屬棒中聲速的有效方法
    探尋中國蘋果產(chǎn)業(yè)的產(chǎn)銷平衡點(diǎn)
    煙臺果樹(2019年1期)2019-01-28 09:34:58
    二階線性微分方程的解法
    電視庭審報(bào)道,如何找到媒體監(jiān)督與司法公正的平衡點(diǎn)
    傳媒評論(2018年7期)2018-09-18 03:45:52
    在給專車服務(wù)正名之前最好找到Uber和出租車的平衡點(diǎn)
    IT時代周刊(2015年7期)2015-11-11 05:49:56
    丰满的人妻完整版| 亚洲av日韩精品久久久久久密| 香蕉久久夜色| 亚洲一区二区三区欧美精品| 少妇被粗大的猛进出69影院| 久久久久久久国产电影| 9热在线视频观看99| 亚洲色图综合在线观看| avwww免费| 在线国产一区二区在线| 一级作爱视频免费观看| 亚洲欧美一区二区三区黑人| 老司机在亚洲福利影院| 国产成人影院久久av| 岛国毛片在线播放| 国产亚洲精品第一综合不卡| www日本在线高清视频| 99热网站在线观看| 这个男人来自地球电影免费观看| 国产免费av片在线观看野外av| 18禁观看日本| 欧美激情久久久久久爽电影 | 波多野结衣一区麻豆| 村上凉子中文字幕在线| 免费在线观看视频国产中文字幕亚洲| 女人被狂操c到高潮| 女性生殖器流出的白浆| 久久久久久人人人人人| 精品高清国产在线一区| 窝窝影院91人妻| 成人影院久久| 高清av免费在线| 无人区码免费观看不卡| 欧美成狂野欧美在线观看| 国产黄色免费在线视频| 大陆偷拍与自拍| 亚洲第一av免费看| 免费一级毛片在线播放高清视频 | 欧美性长视频在线观看| 欧美久久黑人一区二区| 国产成人精品久久二区二区免费| 欧美激情久久久久久爽电影 | 国产黄色免费在线视频| 午夜精品在线福利| 母亲3免费完整高清在线观看| 69av精品久久久久久| 天天操日日干夜夜撸| 又大又爽又粗| 免费看十八禁软件| 国产av精品麻豆| 少妇裸体淫交视频免费看高清 | 亚洲人成伊人成综合网2020| 精品国产一区二区久久| 精品熟女少妇八av免费久了| 免费观看a级毛片全部| 性少妇av在线| 亚洲第一欧美日韩一区二区三区| 超色免费av| 精品免费久久久久久久清纯 | 欧美日韩亚洲国产一区二区在线观看 | 久久久久久免费高清国产稀缺| 老熟女久久久| 欧美 亚洲 国产 日韩一| 日韩免费高清中文字幕av| 一级片'在线观看视频| 在线观看免费日韩欧美大片| 黄片大片在线免费观看| 99国产精品一区二区蜜桃av | 香蕉国产在线看| 午夜福利免费观看在线| 国产又爽黄色视频| 一区二区三区激情视频| 91九色精品人成在线观看| 最近最新中文字幕大全电影3 | 久久久久久久午夜电影 | 亚洲久久久国产精品| 精品少妇久久久久久888优播| 自拍欧美九色日韩亚洲蝌蚪91| 久久午夜亚洲精品久久| 黑丝袜美女国产一区| 亚洲中文字幕日韩| 欧美久久黑人一区二区| x7x7x7水蜜桃| 久久中文字幕一级| 久久中文字幕人妻熟女| 国产精品二区激情视频| 国产欧美日韩一区二区三| 精品国产一区二区三区四区第35| 国产真人三级小视频在线观看| 日本欧美视频一区| 国产高清激情床上av| 精品人妻1区二区| 午夜日韩欧美国产| 欧美中文综合在线视频| 日韩有码中文字幕| 老司机午夜十八禁免费视频| 91麻豆av在线| 视频区图区小说| 高清视频免费观看一区二区| 搡老岳熟女国产| 黑人巨大精品欧美一区二区mp4| 久久久久国产精品人妻aⅴ院 | 三上悠亚av全集在线观看| 又紧又爽又黄一区二区| 一二三四在线观看免费中文在| √禁漫天堂资源中文www| 国产精品成人在线| 51午夜福利影视在线观看| 中文字幕人妻丝袜一区二区| 国产一区有黄有色的免费视频| 丁香欧美五月| 国产精品香港三级国产av潘金莲| 无限看片的www在线观看| 人人妻人人添人人爽欧美一区卜| 精品欧美一区二区三区在线| 精品国产一区二区久久| 如日韩欧美国产精品一区二区三区| а√天堂www在线а√下载 | www.熟女人妻精品国产| 首页视频小说图片口味搜索| av福利片在线| 高清毛片免费观看视频网站 | 99久久人妻综合| 国产成人欧美| 亚洲精品成人av观看孕妇| 欧美精品啪啪一区二区三区| netflix在线观看网站| 国产亚洲欧美98| 18在线观看网站| 国产成人精品无人区| a在线观看视频网站| 51午夜福利影视在线观看| 国产成人精品在线电影| 久久国产精品大桥未久av| 久久人妻福利社区极品人妻图片| 90打野战视频偷拍视频| 90打野战视频偷拍视频| videos熟女内射| 日本五十路高清| 人妻 亚洲 视频| 1024香蕉在线观看| 在线天堂中文资源库| 侵犯人妻中文字幕一二三四区| 757午夜福利合集在线观看| 国产成人精品久久二区二区免费| 国产1区2区3区精品| 亚洲成人手机| 成人永久免费在线观看视频| 精品福利观看| 日本黄色视频三级网站网址 | 日韩欧美国产一区二区入口| 成人特级黄色片久久久久久久| 亚洲色图av天堂| 中文字幕高清在线视频| 韩国av一区二区三区四区| 超碰成人久久| 国产精品1区2区在线观看. | 欧美激情久久久久久爽电影 | 露出奶头的视频| 成人亚洲精品一区在线观看| 99re6热这里在线精品视频| 黄片小视频在线播放| 欧美老熟妇乱子伦牲交| 国产1区2区3区精品| 成人18禁高潮啪啪吃奶动态图| 中文字幕最新亚洲高清| 色尼玛亚洲综合影院| 老司机福利观看| 一级片免费观看大全| 日韩欧美免费精品| 亚洲熟妇熟女久久| 一进一出抽搐gif免费好疼 | 法律面前人人平等表现在哪些方面| av网站在线播放免费| av网站在线播放免费| 亚洲第一欧美日韩一区二区三区| 女人精品久久久久毛片| 精品久久久久久电影网| 三级毛片av免费| 老司机午夜十八禁免费视频| 99久久精品国产亚洲精品| 99国产综合亚洲精品| 美女高潮到喷水免费观看| 国产精品.久久久| 亚洲一区二区三区欧美精品| 国产精品秋霞免费鲁丝片| 99久久综合精品五月天人人| 人妻丰满熟妇av一区二区三区 | 真人做人爱边吃奶动态| 精品午夜福利视频在线观看一区| 国产亚洲欧美98| 正在播放国产对白刺激| 欧美日韩国产mv在线观看视频| 欧美乱妇无乱码| 欧美在线一区亚洲| 69av精品久久久久久| 亚洲一区二区三区欧美精品| 欧美日韩福利视频一区二区| av超薄肉色丝袜交足视频| 后天国语完整版免费观看| 婷婷丁香在线五月| 高潮久久久久久久久久久不卡| 黄色怎么调成土黄色| 亚洲国产中文字幕在线视频| 91大片在线观看| 两人在一起打扑克的视频| 成熟少妇高潮喷水视频| 最新美女视频免费是黄的| 黄色丝袜av网址大全| 亚洲情色 制服丝袜| 一级毛片女人18水好多| 亚洲午夜理论影院| 欧美性长视频在线观看| 国产精品免费视频内射| 热re99久久精品国产66热6| 欧美 日韩 精品 国产| 欧美黄色片欧美黄色片| 久久亚洲精品不卡| 叶爱在线成人免费视频播放| 如日韩欧美国产精品一区二区三区| 欧美中文综合在线视频| 欧美久久黑人一区二区| 在线观看免费高清a一片| 国产成+人综合+亚洲专区| 精品福利观看| 999久久久精品免费观看国产| 999精品在线视频| 91字幕亚洲| av天堂在线播放| 日韩精品免费视频一区二区三区| 丰满的人妻完整版| 在线天堂中文资源库| 欧美日韩亚洲综合一区二区三区_| 精品免费久久久久久久清纯 | 我的亚洲天堂| 看片在线看免费视频| 天天躁日日躁夜夜躁夜夜| 国产视频一区二区在线看| 热99久久久久精品小说推荐| 亚洲精品中文字幕一二三四区| 在线观看www视频免费| av欧美777| 一边摸一边抽搐一进一小说 | 身体一侧抽搐| 国产精品九九99| 男人的好看免费观看在线视频 | av中文乱码字幕在线| 日韩欧美一区二区三区在线观看 | 亚洲av成人av| 人人妻人人澡人人爽人人夜夜| 黄频高清免费视频| 男人操女人黄网站| av福利片在线| 免费av中文字幕在线| 两性午夜刺激爽爽歪歪视频在线观看 | 建设人人有责人人尽责人人享有的| 超碰97精品在线观看| 三上悠亚av全集在线观看| 久久国产精品男人的天堂亚洲| 他把我摸到了高潮在线观看| videosex国产| 亚洲精华国产精华精| 久久精品亚洲av国产电影网| 日本vs欧美在线观看视频| 在线播放国产精品三级| 黑人操中国人逼视频| 9热在线视频观看99| 婷婷成人精品国产| 黑人巨大精品欧美一区二区蜜桃| 在线视频色国产色| 精品人妻熟女毛片av久久网站| 母亲3免费完整高清在线观看| 欧美黄色淫秽网站| 黑人操中国人逼视频| 国产深夜福利视频在线观看| 国产欧美日韩精品亚洲av| 99国产精品免费福利视频| 大香蕉久久网| 国产有黄有色有爽视频| 免费在线观看亚洲国产| 在线观看免费高清a一片| 欧美日韩亚洲高清精品| 国产精品国产av在线观看| 国产午夜精品久久久久久| 母亲3免费完整高清在线观看| 免费不卡黄色视频| 美女扒开内裤让男人捅视频| 黑人欧美特级aaaaaa片| 欧美不卡视频在线免费观看 | 亚洲aⅴ乱码一区二区在线播放 | 在线观看舔阴道视频| 777久久人妻少妇嫩草av网站| 男女午夜视频在线观看| 国产亚洲欧美在线一区二区| 后天国语完整版免费观看| 精品久久蜜臀av无| 亚洲精品国产色婷婷电影| 91九色精品人成在线观看| 999久久久精品免费观看国产| 欧美激情 高清一区二区三区| 人人妻人人爽人人添夜夜欢视频| 交换朋友夫妻互换小说| 两性夫妻黄色片| 亚洲第一欧美日韩一区二区三区| 久久久久久久久免费视频了| 成人永久免费在线观看视频| 成人av一区二区三区在线看| av视频免费观看在线观看| 精品国产乱码久久久久久男人| 成人黄色视频免费在线看| 一级片'在线观看视频| 亚洲av成人一区二区三| 999久久久精品免费观看国产| 极品少妇高潮喷水抽搐| 欧美激情 高清一区二区三区| 满18在线观看网站| 国产国语露脸激情在线看| 久久人妻av系列| 手机成人av网站| 免费黄频网站在线观看国产| 久久香蕉国产精品| 一区二区三区精品91| 免费久久久久久久精品成人欧美视频| 丁香六月欧美| 天堂俺去俺来也www色官网| 久久精品熟女亚洲av麻豆精品| 黄色视频,在线免费观看| 别揉我奶头~嗯~啊~动态视频| 中文字幕最新亚洲高清| 欧美日韩瑟瑟在线播放| 男女午夜视频在线观看| 国产精品九九99| 日韩有码中文字幕| 国产成人av激情在线播放| 久热爱精品视频在线9| 精品久久久久久,| 欧美午夜高清在线| 亚洲精品自拍成人| 精品第一国产精品| 丝袜在线中文字幕| 亚洲中文av在线| 两性夫妻黄色片| 黑人操中国人逼视频| 午夜福利视频在线观看免费| 亚洲精品久久午夜乱码| 99国产精品一区二区三区| 老汉色∧v一级毛片| 日本精品一区二区三区蜜桃| 成人黄色视频免费在线看| 久久国产亚洲av麻豆专区| 午夜视频精品福利| 国产高清激情床上av| 少妇裸体淫交视频免费看高清 | 欧美 日韩 精品 国产| 黄色女人牲交| 国产成人精品无人区| 国产精品欧美亚洲77777| 久久香蕉精品热| 国产精品国产av在线观看| 这个男人来自地球电影免费观看| 欧美大码av| 亚洲国产欧美日韩在线播放| 亚洲全国av大片| 亚洲av美国av| 久久香蕉激情| 欧美日韩国产mv在线观看视频| 国产又爽黄色视频| 午夜福利在线免费观看网站| 黄片大片在线免费观看| 国产97色在线日韩免费| 国产亚洲精品一区二区www | 又大又爽又粗| 欧美大码av| 少妇粗大呻吟视频| 午夜福利,免费看| 国产精品成人在线| 看片在线看免费视频| 人人妻,人人澡人人爽秒播| 成在线人永久免费视频| 首页视频小说图片口味搜索| 欧美日韩av久久| 高清视频免费观看一区二区| 国产成人av教育| 国产成+人综合+亚洲专区| 两个人免费观看高清视频| 久久久久久久久久久久大奶| 免费av中文字幕在线| 黄色成人免费大全| 亚洲,欧美精品.| 黄片大片在线免费观看| 国产成人免费无遮挡视频| 精品国产亚洲在线| 精品视频人人做人人爽| 亚洲av欧美aⅴ国产| 国产精品一区二区在线不卡| 999久久久精品免费观看国产| 成在线人永久免费视频| 久久久久国产一级毛片高清牌| 国产在视频线精品| 久久久国产欧美日韩av| 国产精品久久视频播放| 久久久久精品国产欧美久久久| 欧美亚洲日本最大视频资源| 亚洲成人国产一区在线观看| 男女下面插进去视频免费观看| 国内久久婷婷六月综合欲色啪| 亚洲三区欧美一区| 亚洲成av片中文字幕在线观看| 国产成人一区二区三区免费视频网站| 热99国产精品久久久久久7| 久久精品国产亚洲av香蕉五月 | 香蕉久久夜色| 丝瓜视频免费看黄片| 日韩熟女老妇一区二区性免费视频| 成年人黄色毛片网站| xxx96com| 每晚都被弄得嗷嗷叫到高潮| 午夜福利免费观看在线| 十分钟在线观看高清视频www| 美女福利国产在线| 99国产极品粉嫩在线观看| 90打野战视频偷拍视频| 在线av久久热| 色老头精品视频在线观看| 久久亚洲精品不卡| 很黄的视频免费| 久久久久国产精品人妻aⅴ院 | bbb黄色大片| 成人永久免费在线观看视频| 亚洲第一青青草原| av片东京热男人的天堂| 侵犯人妻中文字幕一二三四区| 国产男女超爽视频在线观看| 久久人人爽av亚洲精品天堂| 午夜视频精品福利| 欧美日韩中文字幕国产精品一区二区三区 | 国产黄色免费在线视频| 久久 成人 亚洲| 十八禁高潮呻吟视频| 久久亚洲真实| 久久人人爽av亚洲精品天堂| 欧美成人午夜精品| 亚洲自偷自拍图片 自拍| 欧美黑人精品巨大| 欧美成狂野欧美在线观看| 99精品在免费线老司机午夜| 老熟女久久久| 亚洲一码二码三码区别大吗| 欧美成人午夜精品| av有码第一页| 精品久久蜜臀av无| 侵犯人妻中文字幕一二三四区| 99热只有精品国产| 99热国产这里只有精品6| 亚洲精品久久午夜乱码| 韩国精品一区二区三区| 欧美黑人欧美精品刺激| 亚洲黑人精品在线| 伦理电影免费视频| 亚洲欧美精品综合一区二区三区| 免费黄频网站在线观看国产| 视频在线观看一区二区三区| 国产99久久九九免费精品| 亚洲少妇的诱惑av| 亚洲一区二区三区欧美精品| 少妇猛男粗大的猛烈进出视频| 巨乳人妻的诱惑在线观看| 亚洲五月婷婷丁香| 热re99久久国产66热| 精品视频人人做人人爽| 一区在线观看完整版| 超碰成人久久| 视频区欧美日本亚洲| svipshipincom国产片| 欧美 日韩 精品 国产| 亚洲av熟女| 亚洲午夜理论影院| 午夜福利,免费看| 欧美日韩亚洲国产一区二区在线观看 | 看片在线看免费视频| 欧美午夜高清在线| 两人在一起打扑克的视频| 国产精品久久久久久精品古装| 人成视频在线观看免费观看| 欧美日韩国产mv在线观看视频| 久久精品国产亚洲av香蕉五月 | 久久久久国产精品人妻aⅴ院 | 欧美国产精品一级二级三级| 涩涩av久久男人的天堂| tocl精华| 欧美日韩精品网址| 国产三级黄色录像| 国产单亲对白刺激| 一进一出抽搐gif免费好疼 | 久久久久久久午夜电影 | 国产激情久久老熟女| av中文乱码字幕在线| 欧美亚洲 丝袜 人妻 在线| 午夜91福利影院| 亚洲成av片中文字幕在线观看| 高清av免费在线| 国产亚洲精品第一综合不卡| 别揉我奶头~嗯~啊~动态视频| 欧美日韩精品网址| 日韩欧美免费精品| 啦啦啦在线免费观看视频4| 精品人妻在线不人妻| 99热网站在线观看| 久久久久久免费高清国产稀缺| 亚洲欧美一区二区三区黑人| 午夜日韩欧美国产| 99re6热这里在线精品视频| 免费av中文字幕在线| 国产xxxxx性猛交| 操美女的视频在线观看| 欧美 日韩 精品 国产| 免费看十八禁软件| 超碰成人久久| 亚洲精品粉嫩美女一区| 高清在线国产一区| tube8黄色片| 操出白浆在线播放| 久久香蕉精品热| 亚洲av日韩在线播放| 在线观看免费高清a一片| 1024香蕉在线观看| 欧美成人免费av一区二区三区 | 91麻豆精品激情在线观看国产 | 无人区码免费观看不卡| 一边摸一边做爽爽视频免费| 天堂中文最新版在线下载| 国产成人免费无遮挡视频| 亚洲精品一二三| 高清av免费在线| 国产日韩欧美亚洲二区| 别揉我奶头~嗯~啊~动态视频| 变态另类成人亚洲欧美熟女 | 国产亚洲精品一区二区www | 免费在线观看日本一区| 国产欧美日韩一区二区三| 美女扒开内裤让男人捅视频| 18禁观看日本| 久久九九热精品免费| 国产精品1区2区在线观看. | 欧美日韩一级在线毛片| 日韩中文字幕欧美一区二区| 免费女性裸体啪啪无遮挡网站| 久久 成人 亚洲| 亚洲久久久国产精品| 成人18禁在线播放| 99re在线观看精品视频| 国产成人免费观看mmmm| 久久久久久久国产电影| 一区在线观看完整版| 黄色 视频免费看| 最近最新免费中文字幕在线| 超碰成人久久| 国产伦人伦偷精品视频| 国产成人精品久久二区二区免费| 五月开心婷婷网| 国产成人av激情在线播放| 亚洲视频免费观看视频| 精品欧美一区二区三区在线| 国产男靠女视频免费网站| 亚洲人成伊人成综合网2020| 韩国精品一区二区三区| 国产精品美女特级片免费视频播放器 | 男人操女人黄网站| 国产亚洲一区二区精品| 男女免费视频国产| 看黄色毛片网站| 亚洲精品国产精品久久久不卡| 看免费av毛片| 欧美黑人欧美精品刺激| 最新美女视频免费是黄的| 久久久精品免费免费高清| 亚洲avbb在线观看| 午夜福利乱码中文字幕| 99久久综合精品五月天人人| 女人高潮潮喷娇喘18禁视频| 无人区码免费观看不卡| 嫩草影视91久久| 三级毛片av免费| 日韩免费av在线播放| 99精品久久久久人妻精品| 欧美亚洲日本最大视频资源| 麻豆av在线久日| 欧美日韩精品网址| 免费女性裸体啪啪无遮挡网站| 无限看片的www在线观看| 亚洲精品中文字幕一二三四区| 午夜91福利影院| 老司机影院毛片| 日韩免费高清中文字幕av| 国产精品av久久久久免费| 男男h啪啪无遮挡| 一个人免费在线观看的高清视频| 一个人免费在线观看的高清视频| 国产精品一区二区免费欧美| 深夜精品福利| 夜夜躁狠狠躁天天躁| 少妇被粗大的猛进出69影院| www.精华液| 久久天堂一区二区三区四区| svipshipincom国产片| 国产不卡av网站在线观看| 丰满人妻熟妇乱又伦精品不卡| 久久久精品免费免费高清| 黄色丝袜av网址大全| 1024香蕉在线观看| 制服人妻中文乱码| 国产精品二区激情视频| 亚洲在线自拍视频| 亚洲精品美女久久av网站| 一级毛片女人18水好多|