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

    基于QR分解的Duffing系統(tǒng)Lyapunov指數(shù)求解方法*

    2014-09-05 03:30:08劉耀峰
    關(guān)鍵詞:吉首振子軌跡

    劉耀峰,楊 喜,舒 婷

    (吉首大學(xué)物理與機(jī)電工程學(xué)院,湖南 吉首 416000)

    基于QR分解的Duffing系統(tǒng)Lyapunov指數(shù)求解方法*

    劉耀峰,楊 喜,舒 婷

    (吉首大學(xué)物理與機(jī)電工程學(xué)院,湖南 吉首 416000)

    Lyapunov指數(shù)衡量了非線性軌跡的穩(wěn)定性和非線性系統(tǒng)的動(dòng)力學(xué)特性,常作為判斷Duffing系統(tǒng)混沌態(tài)和大尺度周期態(tài)的依據(jù).根據(jù)Lyapunov指數(shù)的特征,Duffing系統(tǒng)的策動(dòng)項(xiàng)采用正弦函數(shù)替代余弦函數(shù)方式,提出了一種基于QR分解的Duffing系統(tǒng)Lyapunov指數(shù)的求解方法.Matlab仿真結(jié)果表明了該算法的正確性、可靠性和有效性.

    Duffing系統(tǒng);Lyapunov指數(shù);QR分解

    近年來,隨著弱信號(hào)檢測(cè)技術(shù)和混沌理論的快速發(fā)展,一些研究者已經(jīng)成功地將混沌理論用于信號(hào)檢測(cè)領(lǐng)域,并取得了一些成果.Lyapunov指數(shù)作為混沌系統(tǒng)的一個(gè)重要參量,它既是混沌系統(tǒng)的一個(gè)重要判據(jù),也是信號(hào)檢測(cè)模型中的重要參考量.目前,大多數(shù)關(guān)于Lyapunov指數(shù)的求解采用QR分解方法.Benettin等[1]最早利用GS標(biāo)準(zhǔn)正交化過程對(duì)系統(tǒng)進(jìn)行QR 分解,求解了系統(tǒng)的Lyapunov指數(shù).Lai等[2]提出了用Jacobian法求解系統(tǒng)的Lyapunov指數(shù),這種方法適應(yīng)于有噪聲的環(huán)境,與p-范數(shù)法相比,它的計(jì)算量比較小.然而重復(fù)正交化會(huì)帶來龐大地計(jì)算量,影響計(jì)算效率,Rangaraian等提出了QR法的改進(jìn)算法RHR算法,有效地避免了重復(fù)正交化.Udwadia等[3-4]針對(duì)RHR算法,提出了RHR改進(jìn)算法,提高了計(jì)算效率.張賓[5]對(duì)低維系統(tǒng)Lyapunov指數(shù)的求解做了深入研究,并把Lyapunov指數(shù)引入到微弱信號(hào)檢測(cè)領(lǐng)域.宋春云等[6]把最大Lyapunov指數(shù)作為混沌判據(jù)引入到Duffing系統(tǒng)信號(hào)檢測(cè)模型中,金天等[7]分析了Lyapunov指數(shù)存在統(tǒng)計(jì)特性,確定了系統(tǒng)檢測(cè)概率與誤警概率等計(jì)算方法,李琳等[8]利用Lyapunov指數(shù)來確定系統(tǒng)的檢測(cè)閾值.

    筆者根據(jù)QR分解的基本思想,提出了基于QR分解的Duffing系統(tǒng)Lyapunov指數(shù)方法,詳細(xì)分析了該方法設(shè)計(jì)思路,并根據(jù)該算法編寫程序并仿真,仿真結(jié)果驗(yàn)證了該算法的正確性、可靠性和有效性,為下一步研究弱信號(hào)檢測(cè)提供了平臺(tái).

    1 Lyapunov指數(shù)

    混沌系統(tǒng)的動(dòng)力學(xué)特征可以通過系統(tǒng)的Lyapunov指數(shù)[5]和相軌跡來描述.當(dāng)相空間中相鄰軌跡f(t,x0)和f(t,x0+Δx)隨著時(shí)間推移時(shí),其軌跡將按f(t,x0)=f(t,x0+Δx)eλt規(guī)律相互吸引或離開,把這種相互吸引或離開的平均變化率稱為Lyapunov指數(shù).該指數(shù)衡量了非線性軌跡的穩(wěn)定性,并從統(tǒng)計(jì)特性上反映非線性系統(tǒng)的動(dòng)力學(xué)特性.同時(shí),Lyapunov指數(shù)也量度了混沌系統(tǒng)對(duì)初始值敏感這一特性,即一個(gè)微小擾動(dòng),都將使混沌系統(tǒng)的Lyapunov指數(shù)發(fā)生改變,因此可以利用Lyapunov指數(shù)作為判斷系統(tǒng)混沌態(tài)和大尺度周期態(tài)的依據(jù).

    對(duì)于N維相空間中的連續(xù)動(dòng)力學(xué)系統(tǒng),考慮相圖中第i個(gè)維度方向上相鄰的2個(gè)點(diǎn).在時(shí)間0處和時(shí)間t處,設(shè)該2個(gè)點(diǎn)的間距分別為Pi(0)和Pi(t),那么該系統(tǒng)在相圖的第i個(gè)維度方向上的無量綱值Lyapunov指數(shù)為

    其中:λi為Lyapunov指數(shù);t為演變時(shí)間.

    Lyapunov指數(shù)與相圖中隨時(shí)間演化得到的軌線之間收縮和擴(kuò)張是息息相關(guān)地,在Lyapunov指數(shù)為負(fù)值方向上的軌道是收縮的,運(yùn)動(dòng)相對(duì)來說是穩(wěn)定地且有規(guī)則的,對(duì)初始條件不敏感;相反,在Lyapunov指數(shù)為正值方向上的軌道是分離的,對(duì)初始條件敏感.N維系統(tǒng)擁有N個(gè)Lyapunov指數(shù),系統(tǒng)的所有Lyapunov指數(shù)組成一個(gè)集合,稱為Lyapunov指數(shù)譜.對(duì)于Duffing系統(tǒng)是否存在動(dòng)力學(xué)混沌,可以先求解Duffing系統(tǒng)的Lyapunov指數(shù)譜,再借助于指數(shù)譜中最大Lyapunov指數(shù)是否大于零來直接判斷[9],當(dāng)系統(tǒng)最大Lyapunov指數(shù)大于零時(shí),系統(tǒng)隨時(shí)間演化得到的相軌跡上的相鄰點(diǎn)有排斥離開趨勢(shì),最終處于不穩(wěn)定狀態(tài),即混沌狀態(tài);當(dāng)系統(tǒng)最大Lyapunov指數(shù)小于零時(shí),系統(tǒng)隨時(shí)間演化得到的相軌跡上的相鄰點(diǎn)有吸引收縮趨勢(shì),最終系統(tǒng)會(huì)處于穩(wěn)定狀態(tài),即周期狀態(tài);當(dāng)系統(tǒng)的最大Lyapunov指數(shù)等于零時(shí),系統(tǒng)將處于混沌和周期狀態(tài)之間的臨界狀態(tài)[10].

    2 基于QR分解的Duffing系統(tǒng)Lyapunov指數(shù)求解方法

    QR分解法即Jacobian方法,其基本思想是將已知?jiǎng)恿ο到y(tǒng)的基本解矩陣,分解成正交矩陣Q和對(duì)角線元素都為正的上三角矩陣R的乘積,其中矩陣R是正交矩陣Q和系統(tǒng)Jacobian矩陣的函數(shù),最終根據(jù)上三角矩陣R的對(duì)角線元素可以求得系統(tǒng)的Lyapunov指數(shù).

    考慮混沌Duffing振子的系統(tǒng)方程為

    (1)

    其中:k為阻尼比;-ax(t)+bx3(t)為非線性恢復(fù)力;γsin(ωt)為周期策動(dòng)力;γ為周期策動(dòng)力振幅.

    (2)

    其中令y1(t),y2(t),y3(t)的初始值分別為y1(0)=y10,y2(0)=y20,y3(0)=0.則(2)式的線性變分方程為

    (3)

    其中:I3是3階的單位陣;J(t)是三維自治系統(tǒng)的Jacobian矩陣;Y(t)是(3)式的基本解矩陣,且為

    (4)

    對(duì)矩陣Y(t)進(jìn)行QR分解,記作:Y(t)=Q(t)R(t).經(jīng)過推導(dǎo)可以得到文獻(xiàn)[6,10,12]所示的滿足(3)式的Lyapunov指數(shù)表達(dá)式為

    (5)

    其中Rii(t)是上三角矩陣R(t)的主對(duì)角線元素.利用文獻(xiàn)[12]中的連續(xù)系統(tǒng)中QR分解算法可得

    (6)

    對(duì)ln(Rii(t))求導(dǎo)可得

    (7)

    (7)式兩邊同時(shí)做積分運(yùn)算后,結(jié)合(6)式可得

    因此(5)式可進(jìn)一步化為

    (8)

    其中(QT(t)J(t)Q(t))ii就是矩陣QT(t)J(t)Q(t)的主對(duì)角線元素.

    從上述推導(dǎo)過程可知,利用QR分解算法求解系統(tǒng)的Lyapunov指數(shù),則按照下列步驟求解:

    (1) 將(1)式轉(zhuǎn)化成三維自治系統(tǒng)(2),求出Duffing系統(tǒng)的Jacobian矩陣;

    (2) 求解正交矩陣Q(t);

    (3) 對(duì)(8)式進(jìn)行數(shù)值積分,從而求得Lyapunov指數(shù).

    3 仿真實(shí)驗(yàn)與討論

    混沌Duffing振子方程(1)在k=0.5,a=b=1,ω=1時(shí),系統(tǒng)狀態(tài)將隨著策動(dòng)力振幅γ的變化而變化,同時(shí)Lyapunov指數(shù)也會(huì)發(fā)生相應(yīng)變化.從文獻(xiàn)[6-8,13]可知,Duffing系統(tǒng)存在一個(gè)閾值γd,該閾值是系統(tǒng)從混沌狀態(tài)變化到周期狀態(tài)的臨界值.文中混沌Duffing振子的系統(tǒng)閾值為γd=0.826 0,當(dāng)系統(tǒng)的策動(dòng)力幅值γ<γd時(shí),可以通過QR分解算法求得系統(tǒng)最大Lyapunov指數(shù)大于零,此時(shí)可以判斷系統(tǒng)處于混沌狀態(tài),從相應(yīng)的相軌跡圖中可得到驗(yàn)證;當(dāng)γ>γd時(shí),同樣可以用QR分解算法求得系統(tǒng)最大Lyapunov指數(shù)小于零,此時(shí)可以判斷系統(tǒng)處于周期狀態(tài),同樣也可以從相應(yīng)的相軌跡圖中得到驗(yàn)證.

    設(shè)仿真的時(shí)間步長為0.01,每次迭代步數(shù)為10,總的循環(huán)次數(shù)為100 000次.利用庫塔算法得到Duffing系統(tǒng)的Lyapunov指數(shù)譜.為了減少誤差,提高系統(tǒng)的檢測(cè)概率,應(yīng)刪除不穩(wěn)定迭代,因?yàn)榉抡娴玫降闹笖?shù)譜中均存在一定的過渡帶[14-15],在這些過渡帶中的Lyapunov指數(shù)的數(shù)值,呈現(xiàn)正負(fù)交差狀態(tài).為了準(zhǔn)確的判定系統(tǒng)的狀態(tài),要選取穩(wěn)定的Lyapunov指數(shù)值.文中選取i=10 000點(diǎn)處的穩(wěn)定的Lyapunov指數(shù)值.仿真結(jié)果如圖1~6所示.當(dāng)γ=0.825 9時(shí),相軌跡仿真結(jié)果為圖1的混沌狀態(tài).從圖2的仿真結(jié)果可知,穩(wěn)定時(shí)的最大Lyapunov指數(shù)值約為0.078 5;當(dāng)γ=0.826 0時(shí),相軌跡仿真結(jié)果為圖3的間歇混沌狀態(tài).從圖4的仿真結(jié)果可知,穩(wěn)定時(shí)的最大Lyapunov指數(shù)值約為0.003 1;當(dāng)γ=0.826 1時(shí),相軌跡仿真結(jié)果為圖5的周期狀態(tài).從圖6的仿真結(jié)果可知,穩(wěn)定時(shí)最大Lyapunov指數(shù)值約為-0.012 2,與文獻(xiàn)[9]中的結(jié)論一致.

    圖2 γ=0.825 9時(shí)的混沌狀態(tài)的相軌跡曲線

    圖3 γ=0.826 0時(shí)的Lyapunov指數(shù)演化曲線

    圖4 γ=0.8260時(shí)的臨界狀態(tài)的相軌跡曲線

    圖5 γ=0.826 1時(shí)的Lyapunov指數(shù)演化曲線

    圖6 γ=0.826 1時(shí)的周期狀態(tài)的相軌跡曲線

    4 結(jié)語

    本文分析了用QR分解求解Duffing系統(tǒng)Lyapunov指數(shù)的方法,仿真結(jié)果表明,當(dāng)Duffing系統(tǒng)處于臨界狀態(tài)時(shí),系統(tǒng)再加入一個(gè)微弱的信號(hào)后,系統(tǒng)狀態(tài)將由混沌態(tài)躍遷到周期狀態(tài),且圖形變化明顯.

    [1] BENETTIN G,GALGANI L,GIORGILLI A,et al.Lyapunov Characteristic Exponents for Smooth Dynamical Systems and for Hamiltonian Systems;A Method For Computing All of Them.Part 1:Theory[J].Meccanica,1980,15(1):9-20.

    [2] LAI D,CHEN G.Statistical Analysis of Lyapunov Exponents from Time Series:A Jacobian Approach[J].Mathematical and Computer Modelling,1998,27(7):1-9.

    [3] UDWADIA F E,VON BREMEN H F.An Efficient and Stable Approach for Computation of Lyapunov Characteristic Exponents of Continuous Dynamical Systems[J].Applied Mathematics and Computation,2001,121(2):219-259.

    [4] UDWADIA F E,VON BREMEN H F.Computation of Lyapunov Characteristic Exponents for Continuous Dynamical Systems[J].Zeitschrift FüR Angewandte Mathematik Und Physik ZAMP,2002,53(1):123-146.

    [5] 張 賓.Lyapunov特性指數(shù)的算法研究及其在弱信號(hào)混沌檢測(cè)中的應(yīng)用[D].吉林:吉林大學(xué),2002.

    [6] 宋春云.最大Lyapunov特性指數(shù)在微弱信號(hào)檢測(cè)中的應(yīng)用[J].聲學(xué)技術(shù),2007,26(1):126-129.

    [7] 金 天,張 驊.基于統(tǒng)計(jì)方法的混沌Duffing振子弱信號(hào)檢測(cè)與估計(jì)[J].中國科學(xué):信息科學(xué),2011,41(10):1 184-1 199.

    [8] 李 琳,劉春剛,石 碩,等.基于混沌振子和Lyapunov指數(shù)的微弱信號(hào)檢測(cè)方法[J].黑龍江大學(xué)學(xué)報(bào):自然科學(xué)版,2012,29(4):556-560.

    [9] 李 月,楊寶俊.混沌振子系統(tǒng)(L-Y)與檢測(cè)[M].北京:科學(xué)出版社,2005:83-101.

    [10] 王曉亮.基于多小波變換與Duffing振子的微弱信號(hào)檢測(cè)與估計(jì)[D].哈爾濱:哈爾濱工程大學(xué),2012.

    [11] DIECI L,RUSSELL R D,VAN VLECK E S.On the Compuation of Lyapunov Exponents for Continuous Dynamical Systems[J].SIAM Journal on Numerical Analysis,1997,34(1):402-423.

    [12] MCDONALD E J,HIGHAM D J.Error Analysis of QR Algorithms for Computing Lyapunov Exponents[J].Electronic Transactions on Numerical Analysis,2001(12):234-251.

    [13] 楊 淼,安建平,陳 寧,等.基于Duffing混沌系統(tǒng)的頻譜感知算法[J].北京理工大學(xué)學(xué)報(bào),2011,31(3):329-332.

    [14] 楊紅英,葉 昊,王桂增,等.Duffing振子的Lyapunov指數(shù)與Floquet指數(shù)研究[J].儀器儀表學(xué)報(bào),2008,29(5):927-931.

    [15] 劉海波,吳德偉,金 偉,等.Duffing振子微弱信號(hào)檢測(cè)方法研究[J].物理學(xué)報(bào),2013,62(5):42-47.

    (責(zé)任編輯 陳炳權(quán))

    QRDecompositionBasedMethodfortheComputationofLyapunovExponentinDuffingSystems

    LIU Yaofeng,YANG Xi,SHU Ting

    (College of Physics and Mechanical & Electrical Eengineering,Jishou University,Jishou 416000,Hunan China)

    In order to compute the Lyapunov exponent of Duffing system,this paper adopts a kind of classical method based on QR decomposition to calculate the Lyapunov exponent of Duffing system.According to the basic characteristics of the Lyapunov exponent,we presents a kind of Lyapunov exponent algorithm based on QR decomposition for the known Duffing system.By using software simulation,we know that the Duffing system Lyapunov exponent algorithm based on QR decomposition is correctness.

    Duffing system;Lyapunov exponent;QR decomposition

    1007-2985(2014)01-0058-05

    2013-10-15

    劉耀峰(1985-),男,湖南婁底人,吉首大學(xué)物理與機(jī)電工程學(xué)院無線電物理碩士研究生,主要從事信號(hào)檢測(cè)研究

    楊 喜(1978-),男,湖南湘陰人,吉首大學(xué)信息科學(xué)與工程學(xué)院副教授,博士,碩士生導(dǎo)師,主要從事認(rèn)知無線電研究.

    TP391.41

    A

    10.3969/j.issn.1007-2985.2014.01.014

    猜你喜歡
    吉首振子軌跡
    吉首大學(xué)美術(shù)學(xué)院作品精選
    聲屏世界(2022年15期)2022-11-08 10:58:04
    彈簧振子問題的分析與求解
    湘粵專家學(xué)者相聚吉首研討聲樂套曲《四季如歌》
    軌跡
    軌跡
    吉首美術(shù)館
    軌跡
    非線性Duffing擾動(dòng)振子共振機(jī)制的研究
    進(jìn)化的軌跡(一)——進(jìn)化,無盡的適應(yīng)
    中國三峽(2017年2期)2017-06-09 08:15:29
    基于近似熵和混沌振子的電力諧波檢測(cè)與估計(jì)
    亚洲国产精品合色在线| 十八禁网站免费在线| 日韩精品有码人妻一区| 亚洲 国产 在线| 国产爱豆传媒在线观看| 中文字幕久久专区| 蜜桃亚洲精品一区二区三区| 国产aⅴ精品一区二区三区波| 最新在线观看一区二区三区| 一级av片app| 国产亚洲av嫩草精品影院| 国产真实乱freesex| 少妇裸体淫交视频免费看高清| 亚洲av成人av| 少妇裸体淫交视频免费看高清| 久久久久国内视频| 最新在线观看一区二区三区| 一本精品99久久精品77| 天堂影院成人在线观看| 日韩欧美在线乱码| 身体一侧抽搐| 91久久精品电影网| 久久久久久久久大av| 在线观看美女被高潮喷水网站| 97超视频在线观看视频| 国产91精品成人一区二区三区| 国产一区二区在线av高清观看| 日韩欧美免费精品| 亚洲va在线va天堂va国产| 午夜精品久久久久久毛片777| 成人av在线播放网站| 成人性生交大片免费视频hd| 国内久久婷婷六月综合欲色啪| 又爽又黄无遮挡网站| 男插女下体视频免费在线播放| 长腿黑丝高跟| 欧美日本亚洲视频在线播放| 少妇人妻一区二区三区视频| 欧美+亚洲+日韩+国产| 欧美三级亚洲精品| 亚洲精华国产精华液的使用体验 | 一级a爱片免费观看的视频| 变态另类丝袜制服| 国产高清有码在线观看视频| 熟女人妻精品中文字幕| av黄色大香蕉| 91在线观看av| 亚洲电影在线观看av| 3wmmmm亚洲av在线观看| 国产精品女同一区二区软件 | 婷婷亚洲欧美| 亚洲第一区二区三区不卡| 日本色播在线视频| 久久人人爽人人爽人人片va| 极品教师在线免费播放| 99九九线精品视频在线观看视频| 长腿黑丝高跟| 日韩欧美免费精品| 亚洲五月天丁香| 午夜爱爱视频在线播放| 成年免费大片在线观看| 女的被弄到高潮叫床怎么办 | 日韩欧美精品v在线| 日韩欧美一区二区三区在线观看| 床上黄色一级片| 日日干狠狠操夜夜爽| 真人一进一出gif抽搐免费| 在线看三级毛片| 国产久久久一区二区三区| 亚洲avbb在线观看| 免费人成在线观看视频色| 不卡视频在线观看欧美| 中文亚洲av片在线观看爽| 免费av不卡在线播放| 国语自产精品视频在线第100页| 男女那种视频在线观看| 国产精品不卡视频一区二区| 麻豆久久精品国产亚洲av| 欧美又色又爽又黄视频| 少妇人妻一区二区三区视频| 97超视频在线观看视频| 女人十人毛片免费观看3o分钟| 一区二区三区四区激情视频 | 国产在线男女| 少妇熟女aⅴ在线视频| 午夜福利在线观看免费完整高清在 | 国产高清视频在线播放一区| 国产91精品成人一区二区三区| 日韩欧美免费精品| 欧美色视频一区免费| 一级av片app| 啦啦啦观看免费观看视频高清| 看黄色毛片网站| 成年版毛片免费区| 久久午夜福利片| 亚洲欧美日韩高清专用| 一级av片app| 亚洲av中文字字幕乱码综合| 亚洲国产精品合色在线| 国产一区二区三区视频了| 男女啪啪激烈高潮av片| 精华霜和精华液先用哪个| 国产精品一区二区三区四区久久| 长腿黑丝高跟| 1024手机看黄色片| 真实男女啪啪啪动态图| 伦理电影大哥的女人| 欧美+亚洲+日韩+国产| 我的老师免费观看完整版| 久久热精品热| 此物有八面人人有两片| 午夜福利在线在线| 一级a爱片免费观看的视频| 国产男人的电影天堂91| 欧美性猛交黑人性爽| 欧美激情在线99| 亚洲av免费高清在线观看| 在线免费观看不下载黄p国产 | 日韩一本色道免费dvd| 午夜影院日韩av| 成年女人永久免费观看视频| 国产成人影院久久av| 久久草成人影院| 日韩在线高清观看一区二区三区 | 在线观看av片永久免费下载| 国产私拍福利视频在线观看| 搡老妇女老女人老熟妇| 久久久久性生活片| 午夜视频国产福利| 悠悠久久av| 搡女人真爽免费视频火全软件 | eeuss影院久久| 狂野欧美激情性xxxx在线观看| 国产黄色小视频在线观看| 在现免费观看毛片| 蜜桃久久精品国产亚洲av| 欧美日韩国产亚洲二区| 国内少妇人妻偷人精品xxx网站| 国内揄拍国产精品人妻在线| av在线天堂中文字幕| 在线免费十八禁| 尤物成人国产欧美一区二区三区| 成人精品一区二区免费| 久久香蕉精品热| 国内揄拍国产精品人妻在线| 日韩欧美免费精品| 国产中年淑女户外野战色| 国产精品野战在线观看| 久久亚洲精品不卡| 麻豆国产av国片精品| 在线播放无遮挡| 18禁黄网站禁片免费观看直播| 欧美在线一区亚洲| 国产激情偷乱视频一区二区| 亚洲美女视频黄频| 在线看三级毛片| 网址你懂的国产日韩在线| 国产精品免费一区二区三区在线| 亚洲人成伊人成综合网2020| av天堂中文字幕网| 乱人视频在线观看| av黄色大香蕉| 人人妻,人人澡人人爽秒播| 欧美成人性av电影在线观看| 99热这里只有是精品50| 在线观看av片永久免费下载| 国产国拍精品亚洲av在线观看| 亚洲av日韩精品久久久久久密| 国产免费一级a男人的天堂| 久久久久国产精品人妻aⅴ院| 日本与韩国留学比较| 美女被艹到高潮喷水动态| 成人国产一区最新在线观看| 久久精品影院6| 久久久久国内视频| 精品日产1卡2卡| 全区人妻精品视频| 丰满乱子伦码专区| 午夜久久久久精精品| 亚洲 国产 在线| 丝袜美腿在线中文| 国产精品不卡视频一区二区| 亚洲av免费在线观看| 免费在线观看成人毛片| 久久精品91蜜桃| 18禁在线播放成人免费| 99热精品在线国产| 男女做爰动态图高潮gif福利片| 国产精华一区二区三区| 日本熟妇午夜| 亚洲国产欧美人成| 天堂影院成人在线观看| 婷婷色综合大香蕉| x7x7x7水蜜桃| 精品不卡国产一区二区三区| 国内毛片毛片毛片毛片毛片| 人妻丰满熟妇av一区二区三区| 午夜福利在线在线| 1024手机看黄色片| 国产高清不卡午夜福利| 日本免费a在线| 国产综合懂色| 黄片wwwwww| 久久久国产成人免费| 内射极品少妇av片p| 99热只有精品国产| 校园春色视频在线观看| 精品人妻1区二区| 成人无遮挡网站| 中文字幕av在线有码专区| 搡女人真爽免费视频火全软件 | 丝袜美腿在线中文| 日日摸夜夜添夜夜添av毛片 | 久久久精品大字幕| 草草在线视频免费看| 国产成人av教育| 久久久久九九精品影院| 欧美+日韩+精品| 88av欧美| 日本在线视频免费播放| 国产三级在线视频| 色播亚洲综合网| 成人精品一区二区免费| 天堂动漫精品| 99热这里只有精品一区| 午夜免费成人在线视频| 人妻少妇偷人精品九色| 三级男女做爰猛烈吃奶摸视频| 精品久久国产蜜桃| 国产高清三级在线| 国产不卡一卡二| 国产真实乱freesex| 亚洲色图av天堂| a级毛片a级免费在线| 精品久久久久久久久久久久久| 久久热精品热| 精品欧美国产一区二区三| 又黄又爽又免费观看的视频| or卡值多少钱| 精品人妻一区二区三区麻豆 | 人人妻,人人澡人人爽秒播| 午夜福利在线观看吧| 亚洲成av人片在线播放无| 在线免费观看的www视频| 亚洲欧美激情综合另类| 91久久精品电影网| 精品久久久久久久末码| 国产探花极品一区二区| av在线天堂中文字幕| 午夜精品一区二区三区免费看| 国产精品久久久久久亚洲av鲁大| 午夜老司机福利剧场| 欧美激情国产日韩精品一区| 日韩强制内射视频| 成人欧美大片| 五月玫瑰六月丁香| 中文在线观看免费www的网站| 人人妻,人人澡人人爽秒播| 男人的好看免费观看在线视频| 在线播放无遮挡| 精品午夜福利在线看| 亚洲久久久久久中文字幕| 人人妻,人人澡人人爽秒播| 免费观看的影片在线观看| 久久久久国产精品人妻aⅴ院| 国产国拍精品亚洲av在线观看| 亚洲人成伊人成综合网2020| 国产高潮美女av| 69av精品久久久久久| 午夜免费男女啪啪视频观看 | 亚洲精华国产精华精| 日韩欧美在线乱码| 久久精品综合一区二区三区| 99精品在免费线老司机午夜| 亚洲精品色激情综合| 一卡2卡三卡四卡精品乱码亚洲| 成人毛片a级毛片在线播放| 天天躁日日操中文字幕| 亚洲国产精品合色在线| 在线观看美女被高潮喷水网站| 国产 一区 欧美 日韩| 婷婷色综合大香蕉| 色视频www国产| 久久热精品热| 人妻少妇偷人精品九色| 两个人的视频大全免费| 日本一二三区视频观看| 女同久久另类99精品国产91| 男女视频在线观看网站免费| 欧美激情久久久久久爽电影| 一夜夜www| 小蜜桃在线观看免费完整版高清| 精品乱码久久久久久99久播| 亚洲成a人片在线一区二区| 1000部很黄的大片| 1024手机看黄色片| 欧美区成人在线视频| 久久6这里有精品| 中文字幕高清在线视频| 国产一区二区激情短视频| 尤物成人国产欧美一区二区三区| 精品久久久久久久久亚洲 | 久久人人精品亚洲av| 久久久国产成人免费| 欧美区成人在线视频| 69av精品久久久久久| 日日摸夜夜添夜夜添av毛片 | 在线看三级毛片| 亚洲第一区二区三区不卡| 国产精品日韩av在线免费观看| 简卡轻食公司| 国产高潮美女av| 国产亚洲精品久久久com| 黄色配什么色好看| 两个人视频免费观看高清| 日日摸夜夜添夜夜添av毛片 | 久久久久久久久久成人| 极品教师在线免费播放| 在线观看舔阴道视频| 久久精品综合一区二区三区| 中文在线观看免费www的网站| 久久精品91蜜桃| 观看美女的网站| 在线播放国产精品三级| 午夜福利高清视频| 一个人免费在线观看电影| 中出人妻视频一区二区| 国产视频一区二区在线看| 天堂√8在线中文| 桃色一区二区三区在线观看| 亚洲色图av天堂| 国产男人的电影天堂91| xxxwww97欧美| 亚洲欧美日韩东京热| 国产精品,欧美在线| 亚洲四区av| 亚洲午夜理论影院| 精品久久久久久久久亚洲 | 色综合色国产| 天堂√8在线中文| av在线老鸭窝| 丰满的人妻完整版| 久久人妻av系列| www.色视频.com| 99热这里只有精品一区| 国产伦精品一区二区三区视频9| 久久精品国产99精品国产亚洲性色| 色哟哟哟哟哟哟| 亚洲av免费高清在线观看| 黄色配什么色好看| 露出奶头的视频| 国产aⅴ精品一区二区三区波| 22中文网久久字幕| 亚洲黑人精品在线| 小蜜桃在线观看免费完整版高清| 一a级毛片在线观看| 人妻少妇偷人精品九色| 婷婷亚洲欧美| 久久人人爽人人爽人人片va| 成年免费大片在线观看| 精品福利观看| 亚洲国产欧美人成| 校园春色视频在线观看| 国产高清视频在线播放一区| 日本爱情动作片www.在线观看 | 国产色爽女视频免费观看| 村上凉子中文字幕在线| 免费在线观看成人毛片| 免费av不卡在线播放| 亚洲无线在线观看| 欧美成人免费av一区二区三区| 成人特级黄色片久久久久久久| 欧美一区二区国产精品久久精品| 午夜精品一区二区三区免费看| 91麻豆精品激情在线观看国产| 九九爱精品视频在线观看| 国产高清激情床上av| 亚洲avbb在线观看| 免费av毛片视频| 黄色丝袜av网址大全| 91麻豆精品激情在线观看国产| 色av中文字幕| 国产精品av视频在线免费观看| 91精品国产九色| 黄片wwwwww| 日韩欧美在线二视频| 国产免费av片在线观看野外av| 欧美日韩亚洲国产一区二区在线观看| 日韩中字成人| 亚洲图色成人| 亚洲av美国av| 窝窝影院91人妻| 国产一区二区亚洲精品在线观看| 亚洲av中文字字幕乱码综合| 偷拍熟女少妇极品色| 亚洲,欧美,日韩| 国产伦人伦偷精品视频| 成人国产综合亚洲| 日韩欧美精品免费久久| 悠悠久久av| bbb黄色大片| 男人舔奶头视频| 18禁黄网站禁片午夜丰满| 又紧又爽又黄一区二区| 欧美bdsm另类| 精品久久国产蜜桃| 午夜久久久久精精品| 亚洲欧美日韩无卡精品| 别揉我奶头~嗯~啊~动态视频| 91狼人影院| 日本精品一区二区三区蜜桃| 夜夜夜夜夜久久久久| 三级男女做爰猛烈吃奶摸视频| 99久久无色码亚洲精品果冻| 女同久久另类99精品国产91| 国产淫片久久久久久久久| 大型黄色视频在线免费观看| 又爽又黄无遮挡网站| 少妇丰满av| 国产在视频线在精品| 国产精品人妻久久久影院| 日日摸夜夜添夜夜添小说| 国产精品爽爽va在线观看网站| 久久久久久久久中文| 国内毛片毛片毛片毛片毛片| 成人国产综合亚洲| 免费观看的影片在线观看| 91久久精品电影网| 免费看av在线观看网站| 久久亚洲精品不卡| 成人综合一区亚洲| 极品教师在线免费播放| 精品福利观看| 好男人在线观看高清免费视频| 国内精品美女久久久久久| 国产麻豆成人av免费视频| av在线观看视频网站免费| 婷婷丁香在线五月| 91麻豆精品激情在线观看国产| 国产 一区 欧美 日韩| 少妇的逼好多水| 欧美高清成人免费视频www| 亚洲av一区综合| h日本视频在线播放| 久久国内精品自在自线图片| 国产三级中文精品| 如何舔出高潮| 亚洲av美国av| 免费av观看视频| 嫩草影院新地址| 三级男女做爰猛烈吃奶摸视频| 俺也久久电影网| 成人毛片a级毛片在线播放| 日本撒尿小便嘘嘘汇集6| www日本黄色视频网| 直男gayav资源| 国产男人的电影天堂91| 色5月婷婷丁香| 91久久精品国产一区二区成人| 国产高清不卡午夜福利| 久久久国产成人免费| 国产男靠女视频免费网站| 老司机福利观看| 国产精品免费一区二区三区在线| 国产探花在线观看一区二区| 在线观看美女被高潮喷水网站| 成人高潮视频无遮挡免费网站| 日韩av在线大香蕉| 我的老师免费观看完整版| 国产视频内射| 亚洲美女视频黄频| 国产av麻豆久久久久久久| 午夜亚洲福利在线播放| 99视频精品全部免费 在线| 国产黄色小视频在线观看| 国产精品一区二区三区四区免费观看 | 99热网站在线观看| 波多野结衣巨乳人妻| 免费观看人在逋| 麻豆成人av在线观看| 99精品久久久久人妻精品| 欧美日本亚洲视频在线播放| 精品久久久久久久久久久久久| 88av欧美| 亚洲性久久影院| 欧美xxxx黑人xx丫x性爽| 成人毛片a级毛片在线播放| 国产单亲对白刺激| 国产精品,欧美在线| 人人妻人人看人人澡| 亚洲aⅴ乱码一区二区在线播放| 精品不卡国产一区二区三区| 午夜激情福利司机影院| 国产高清激情床上av| 国内精品久久久久精免费| 国内毛片毛片毛片毛片毛片| 国产熟女欧美一区二区| 国产精品一区www在线观看 | 欧美日韩国产亚洲二区| 国产成人一区二区在线| 一区福利在线观看| 简卡轻食公司| 亚洲国产精品合色在线| 欧美日本亚洲视频在线播放| 嫩草影院入口| 国产精品人妻久久久久久| 国内揄拍国产精品人妻在线| 在线免费观看不下载黄p国产 | 欧美日韩精品成人综合77777| 日韩欧美免费精品| 国产单亲对白刺激| 一个人看的www免费观看视频| 国产亚洲91精品色在线| 老司机午夜福利在线观看视频| 国产亚洲精品久久久久久毛片| 免费高清视频大片| 在线播放国产精品三级| 亚洲精品成人久久久久久| 国产精品一区www在线观看 | 日韩强制内射视频| 日本一本二区三区精品| 亚洲综合色惰| 啦啦啦啦在线视频资源| 欧美日韩综合久久久久久 | 国产精华一区二区三区| 成人高潮视频无遮挡免费网站| 天堂√8在线中文| 欧美国产日韩亚洲一区| 国产精品无大码| 国产成人av教育| 免费在线观看影片大全网站| 免费搜索国产男女视频| 真人做人爱边吃奶动态| 国产私拍福利视频在线观看| 欧美日韩亚洲国产一区二区在线观看| 国产老妇女一区| 亚洲精华国产精华精| 亚洲欧美日韩东京热| 亚洲狠狠婷婷综合久久图片| 最近最新中文字幕大全电影3| 国产 一区精品| 97超级碰碰碰精品色视频在线观看| 国产一区二区三区视频了| 啪啪无遮挡十八禁网站| 国产精品99久久久久久久久| 国国产精品蜜臀av免费| 久久久久久久亚洲中文字幕| 草草在线视频免费看| 蜜桃亚洲精品一区二区三区| 国产高清视频在线播放一区| 欧美xxxx黑人xx丫x性爽| 国产精品1区2区在线观看.| 欧美在线一区亚洲| 亚洲va日本ⅴa欧美va伊人久久| 丰满人妻一区二区三区视频av| 久久久久性生活片| 成人国产一区最新在线观看| 国产精品久久久久久av不卡| 日本熟妇午夜| 极品教师在线免费播放| 国产伦精品一区二区三区四那| 精品一区二区免费观看| 免费看日本二区| 精品人妻视频免费看| 国内精品久久久久久久电影| 夜夜看夜夜爽夜夜摸| 天天躁日日操中文字幕| 特级一级黄色大片| 99九九线精品视频在线观看视频| 男女那种视频在线观看| 日韩精品有码人妻一区| 在线看三级毛片| 村上凉子中文字幕在线| 在线看三级毛片| 人人妻人人澡欧美一区二区| 麻豆av噜噜一区二区三区| 亚洲真实伦在线观看| 麻豆av噜噜一区二区三区| 国产真实伦视频高清在线观看 | 亚洲国产精品合色在线| 成人午夜高清在线视频| 特级一级黄色大片| 欧美成人性av电影在线观看| 亚洲最大成人中文| 免费观看人在逋| 日韩欧美在线乱码| 噜噜噜噜噜久久久久久91| av在线亚洲专区| 校园人妻丝袜中文字幕| 人妻制服诱惑在线中文字幕| 最近中文字幕高清免费大全6 | 人妻夜夜爽99麻豆av| 熟女人妻精品中文字幕| 国产精品久久久久久亚洲av鲁大| bbb黄色大片| 又粗又爽又猛毛片免费看| 亚洲欧美激情综合另类| 日本与韩国留学比较| 桃色一区二区三区在线观看| 国产午夜精品久久久久久一区二区三区 | 99热精品在线国产| 麻豆国产97在线/欧美| 午夜免费男女啪啪视频观看 | 少妇的逼好多水| 在线国产一区二区在线| 久久久久久国产a免费观看| 日韩在线高清观看一区二区三区 | 国产欧美日韩一区二区精品| 三级国产精品欧美在线观看| 人妻少妇偷人精品九色| 欧美+亚洲+日韩+国产| 国产亚洲欧美98| 日本色播在线视频| 精品人妻一区二区三区麻豆 | 夜夜看夜夜爽夜夜摸|