• <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ì)
    法律面前人人平等表现在哪些方面| 一个人免费在线观看的高清视频| 亚洲精品一卡2卡三卡4卡5卡| 成人一区二区视频在线观看| 高清毛片免费观看视频网站| 欧美激情 高清一区二区三区| 在线观看66精品国产| 精品午夜福利视频在线观看一区| www.999成人在线观看| 不卡一级毛片| 精品久久久久久久久久久久久 | 免费女性裸体啪啪无遮挡网站| 久久人人精品亚洲av| 女性被躁到高潮视频| 精品日产1卡2卡| 人成视频在线观看免费观看| 日韩成人在线观看一区二区三区| 亚洲av成人一区二区三| 老司机午夜福利在线观看视频| 午夜激情福利司机影院| 嫩草影院精品99| 久久国产乱子伦精品免费另类| 在线十欧美十亚洲十日本专区| 亚洲avbb在线观看| 99riav亚洲国产免费| 亚洲,欧美精品.| 亚洲自拍偷在线| 国产亚洲精品av在线| 校园春色视频在线观看| 亚洲aⅴ乱码一区二区在线播放 | 不卡av一区二区三区| 美女扒开内裤让男人捅视频| 国产高清有码在线观看视频 | 精品国产乱子伦一区二区三区| 最近在线观看免费完整版| 亚洲国产欧美一区二区综合| av中文乱码字幕在线| 91九色精品人成在线观看| 熟女少妇亚洲综合色aaa.| 亚洲国产欧美一区二区综合| 亚洲熟妇熟女久久| 国产男靠女视频免费网站| 韩国av一区二区三区四区| 精品欧美一区二区三区在线| 久久天躁狠狠躁夜夜2o2o| 宅男免费午夜| 亚洲av成人不卡在线观看播放网| 一本大道久久a久久精品| 老司机福利观看| 啪啪无遮挡十八禁网站| 久久中文看片网| 嫩草影视91久久| 欧美性长视频在线观看| 亚洲一卡2卡3卡4卡5卡精品中文| 亚洲av电影在线进入| 国产亚洲精品第一综合不卡| 国产精品免费视频内射| 国产99久久九九免费精品| 99久久精品国产亚洲精品| av超薄肉色丝袜交足视频| 国产精品日韩av在线免费观看| 搡老岳熟女国产| 国产成人欧美| 久久久久免费精品人妻一区二区 | 男女那种视频在线观看| e午夜精品久久久久久久| 69av精品久久久久久| 天堂√8在线中文| 人人妻人人澡欧美一区二区| 熟女电影av网| 成人永久免费在线观看视频| www.熟女人妻精品国产| 国产单亲对白刺激| 国产精品一区二区三区四区久久 | 亚洲国产高清在线一区二区三 | 国内精品久久久久久久电影| 男女午夜视频在线观看| 老司机深夜福利视频在线观看| svipshipincom国产片| 69av精品久久久久久| 国产私拍福利视频在线观看| 999久久久精品免费观看国产| 色播在线永久视频| 91麻豆av在线| tocl精华| 天天添夜夜摸| 国产91精品成人一区二区三区| 国产精品1区2区在线观看.| 久久人妻av系列| 在线观看日韩欧美| 亚洲狠狠婷婷综合久久图片| 亚洲精品美女久久av网站| 国产精品一区二区三区四区久久 | 欧美+亚洲+日韩+国产| 国产黄片美女视频| netflix在线观看网站| 国产在线精品亚洲第一网站| 色综合欧美亚洲国产小说| 精品熟女少妇八av免费久了| 满18在线观看网站| 天堂√8在线中文| 国产精品,欧美在线| 精华霜和精华液先用哪个| 自线自在国产av| 亚洲 欧美 日韩 在线 免费| 99久久99久久久精品蜜桃| 黄片播放在线免费| 亚洲成a人片在线一区二区| 青草久久国产| 在线观看午夜福利视频| 91av网站免费观看| 在线天堂中文资源库| 国产私拍福利视频在线观看| 成人三级黄色视频| 久久婷婷成人综合色麻豆| 亚洲熟女毛片儿| 侵犯人妻中文字幕一二三四区| 久久久久国产精品人妻aⅴ院| 久9热在线精品视频| 国产三级黄色录像| 精品乱码久久久久久99久播| 久久狼人影院| 日本黄色视频三级网站网址| 久久久久国产一级毛片高清牌| 成人欧美大片| 国产高清视频在线播放一区| 午夜激情av网站| aaaaa片日本免费| 欧美激情极品国产一区二区三区| 91麻豆精品激情在线观看国产| 久久精品夜夜夜夜夜久久蜜豆 | 日韩欧美国产在线观看| 91麻豆av在线| 最新在线观看一区二区三区| 最近在线观看免费完整版| 黄片大片在线免费观看| 性色av乱码一区二区三区2| 亚洲精品国产精品久久久不卡| 美女高潮喷水抽搐中文字幕| 久久99热这里只有精品18| 97人妻精品一区二区三区麻豆 | 国内揄拍国产精品人妻在线 | 国产av不卡久久| 亚洲第一青青草原| 丰满的人妻完整版| 一区二区日韩欧美中文字幕| 亚洲欧美日韩无卡精品| 91字幕亚洲| 嫩草影视91久久| 嫩草影视91久久| 999久久久国产精品视频| 妹子高潮喷水视频| 欧美三级亚洲精品| 亚洲精品一卡2卡三卡4卡5卡| 老司机深夜福利视频在线观看| 国产亚洲精品一区二区www| 国产成人精品久久二区二区免费| 一级毛片精品| 亚洲在线自拍视频| 国产人伦9x9x在线观看| 午夜久久久在线观看| 精品欧美一区二区三区在线| 首页视频小说图片口味搜索| x7x7x7水蜜桃| 人人妻人人澡人人看| 欧美色欧美亚洲另类二区| 老司机午夜十八禁免费视频| 中文字幕人妻熟女乱码| 成人手机av| 一区二区三区精品91| 日本a在线网址| 精品国产乱码久久久久久男人| 日本精品一区二区三区蜜桃| 久久中文看片网| 91麻豆精品激情在线观看国产| 久热这里只有精品99| 国产一区二区激情短视频| 啦啦啦 在线观看视频| 88av欧美| 国产亚洲精品久久久久5区| 99久久无色码亚洲精品果冻| 给我免费播放毛片高清在线观看| 国产成人精品无人区| 黑丝袜美女国产一区| 亚洲人成77777在线视频| 久久久久久大精品| 哪里可以看免费的av片| 欧美黄色片欧美黄色片| 国产一卡二卡三卡精品| 欧美激情 高清一区二区三区| 宅男免费午夜| 老汉色av国产亚洲站长工具| 亚洲熟妇中文字幕五十中出| 88av欧美| 哪里可以看免费的av片| 久久天堂一区二区三区四区| 国产免费av片在线观看野外av| 日本 av在线| 国产欧美日韩一区二区三| 精品国产美女av久久久久小说| 久久精品aⅴ一区二区三区四区| 久久精品夜夜夜夜夜久久蜜豆 | 国产精品永久免费网站| 成人亚洲精品av一区二区| 成人三级做爰电影| 美女扒开内裤让男人捅视频| 伦理电影免费视频| 国产av一区在线观看免费| 欧美性长视频在线观看| 亚洲国产日韩欧美精品在线观看 | av视频在线观看入口| 久久久久久久精品吃奶| 欧美zozozo另类| 国产精品久久视频播放| 一本一本综合久久| 亚洲免费av在线视频| 欧美另类亚洲清纯唯美| 中文字幕最新亚洲高清| 夜夜躁狠狠躁天天躁| 国产成人精品久久二区二区免费| 91在线观看av| 侵犯人妻中文字幕一二三四区| 久久婷婷成人综合色麻豆| 老熟妇仑乱视频hdxx| 午夜影院日韩av| 男女下面进入的视频免费午夜 | 久久久久精品国产欧美久久久| 在线av久久热| 熟女少妇亚洲综合色aaa.| svipshipincom国产片| 天堂影院成人在线观看| 久久 成人 亚洲| av电影中文网址| 在线观看www视频免费| 人人妻,人人澡人人爽秒播| 亚洲国产精品999在线| 久久中文看片网| 国产精品爽爽va在线观看网站 | 黄色片一级片一级黄色片| 人妻久久中文字幕网| 中文字幕最新亚洲高清| 看片在线看免费视频| 老熟妇乱子伦视频在线观看| 久久久久久人人人人人| 欧美av亚洲av综合av国产av| 国产成年人精品一区二区| 中文字幕另类日韩欧美亚洲嫩草| 麻豆av在线久日| aaaaa片日本免费| 亚洲欧洲精品一区二区精品久久久| 国产成人精品久久二区二区免费| 亚洲av五月六月丁香网| 伊人久久大香线蕉亚洲五| 成人三级黄色视频| 99精品欧美一区二区三区四区| 国产精品av久久久久免费| 国产精品二区激情视频| 久久欧美精品欧美久久欧美| 一二三四在线观看免费中文在| av电影中文网址| 91老司机精品| 婷婷丁香在线五月| 久久这里只有精品19| 久久亚洲真实| 俄罗斯特黄特色一大片| 19禁男女啪啪无遮挡网站| avwww免费| 精品久久久久久,| 欧美一区二区精品小视频在线| 欧美成人午夜精品| 日韩精品中文字幕看吧| 麻豆一二三区av精品| 级片在线观看| 欧美另类亚洲清纯唯美| 成人18禁在线播放| 亚洲人成网站高清观看| 成年女人毛片免费观看观看9| 中文字幕精品亚洲无线码一区 | 国产成人精品久久二区二区免费| 亚洲色图av天堂| 久久欧美精品欧美久久欧美| 国产蜜桃级精品一区二区三区| 国产精品自产拍在线观看55亚洲| 俺也久久电影网| 国产精品亚洲av一区麻豆| 精品一区二区三区av网在线观看| 少妇熟女aⅴ在线视频| 伊人久久大香线蕉亚洲五| 国产91精品成人一区二区三区| 久久久国产精品麻豆| 日韩大码丰满熟妇| 欧美色欧美亚洲另类二区| 国产精品一区二区免费欧美| 国产精品亚洲美女久久久| 亚洲 欧美 日韩 在线 免费| av在线播放免费不卡| 午夜福利免费观看在线| 久久九九热精品免费| 老司机午夜福利在线观看视频| 免费一级毛片在线播放高清视频| 白带黄色成豆腐渣| 亚洲美女黄片视频| 美女高潮到喷水免费观看| av超薄肉色丝袜交足视频| 人人妻人人澡人人看| 午夜免费鲁丝| 侵犯人妻中文字幕一二三四区| 久久久国产欧美日韩av| 成人亚洲精品一区在线观看| 色精品久久人妻99蜜桃| 首页视频小说图片口味搜索| 日日摸夜夜添夜夜添小说| 国产亚洲精品久久久久久毛片| 亚洲国产精品sss在线观看| 黄色视频不卡| 中国美女看黄片| 老司机福利观看| 中文资源天堂在线| 久久久久久九九精品二区国产 | 久久精品人妻少妇| 免费av毛片视频| 成人午夜高清在线视频 | 成人av一区二区三区在线看| 免费电影在线观看免费观看| 色播亚洲综合网| 99久久综合精品五月天人人| 一级毛片高清免费大全| 亚洲成人免费电影在线观看| 国产不卡一卡二| 久久人妻福利社区极品人妻图片| 国产熟女午夜一区二区三区| 国产久久久一区二区三区| 好看av亚洲va欧美ⅴa在| 亚洲国产欧美日韩在线播放| 一进一出好大好爽视频| 午夜亚洲福利在线播放| 人人妻人人澡人人看| 日韩大尺度精品在线看网址| 中文亚洲av片在线观看爽| 亚洲精品久久成人aⅴ小说| 大型av网站在线播放| 久久久久国内视频| 亚洲成国产人片在线观看| 麻豆一二三区av精品| 熟女电影av网| 在线观看一区二区三区| 国产精品av久久久久免费| 国产极品粉嫩免费观看在线| 最近最新中文字幕大全电影3 | 91字幕亚洲| 一二三四社区在线视频社区8| 18禁裸乳无遮挡免费网站照片 | 国产99久久九九免费精品| 在线观看66精品国产| 国产在线观看jvid| 成人午夜高清在线视频 | 变态另类成人亚洲欧美熟女| 精品久久久久久,| 亚洲人成77777在线视频| 99在线视频只有这里精品首页| 精品一区二区三区av网在线观看| 婷婷精品国产亚洲av| 搡老岳熟女国产| 精品国产乱码久久久久久男人| 亚洲色图 男人天堂 中文字幕| 欧美国产日韩亚洲一区| 色婷婷久久久亚洲欧美| 2021天堂中文幕一二区在线观 | 国产精品永久免费网站| 国产精品99久久99久久久不卡| 一a级毛片在线观看| 国产精品av久久久久免费| bbb黄色大片| 在线观看日韩欧美| 黄色女人牲交| e午夜精品久久久久久久| 国产久久久一区二区三区| 国产精品亚洲一级av第二区| 亚洲av电影在线进入| 在线播放国产精品三级| 日本免费一区二区三区高清不卡| 欧美一区二区精品小视频在线| 午夜两性在线视频| 国产爱豆传媒在线观看 | av欧美777| 69av精品久久久久久| www.精华液| 精品福利观看| 日本 av在线| 99riav亚洲国产免费| 国内精品久久久久精免费| 妹子高潮喷水视频| 国产午夜精品久久久久久| 欧美日韩福利视频一区二区| 精品国内亚洲2022精品成人| 国产成人欧美| 国产精品爽爽va在线观看网站 | 亚洲欧美日韩高清在线视频| 悠悠久久av| 19禁男女啪啪无遮挡网站| 非洲黑人性xxxx精品又粗又长| 精品国产国语对白av| 美女午夜性视频免费| 1024视频免费在线观看| 久久人妻av系列| 999久久久精品免费观看国产| 天堂√8在线中文| 午夜激情福利司机影院| 一本大道久久a久久精品| 男人舔女人下体高潮全视频| 色精品久久人妻99蜜桃| 亚洲人成网站高清观看| 天天躁夜夜躁狠狠躁躁| 国产精品自产拍在线观看55亚洲| 成年女人毛片免费观看观看9| 午夜激情av网站| 国产精品免费一区二区三区在线| 欧美日韩亚洲综合一区二区三区_| 欧美成人一区二区免费高清观看 | 中文字幕精品亚洲无线码一区 | 国内精品久久久久久久电影| 国产精品亚洲av一区麻豆| 亚洲天堂国产精品一区在线| 亚洲av成人一区二区三| 中文亚洲av片在线观看爽| 此物有八面人人有两片| 久久精品aⅴ一区二区三区四区| 久久香蕉国产精品| 夜夜爽天天搞| 最新在线观看一区二区三区| 亚洲欧美日韩高清在线视频| 色婷婷久久久亚洲欧美| 黄色 视频免费看| 男女做爰动态图高潮gif福利片| 久久久久久久午夜电影| 色精品久久人妻99蜜桃| 亚洲 欧美一区二区三区| 国产精品久久久av美女十八| 久久久久久久精品吃奶| 久久热在线av| 亚洲真实伦在线观看| 一级片免费观看大全| 午夜久久久在线观看| 午夜激情av网站| 99在线人妻在线中文字幕| 在线av久久热| 欧美成人午夜精品| 91成年电影在线观看| 婷婷丁香在线五月| 亚洲中文字幕一区二区三区有码在线看 | 日韩精品免费视频一区二区三区| 国产亚洲av嫩草精品影院| 亚洲午夜理论影院| 久热这里只有精品99| 一进一出抽搐动态| 欧美绝顶高潮抽搐喷水| 欧美激情 高清一区二区三区| 非洲黑人性xxxx精品又粗又长| 久久精品亚洲精品国产色婷小说| 欧美午夜高清在线| 非洲黑人性xxxx精品又粗又长| 免费无遮挡裸体视频| 长腿黑丝高跟| 少妇的丰满在线观看| 身体一侧抽搐| av视频在线观看入口| 我的亚洲天堂| 亚洲国产毛片av蜜桃av| 久久午夜亚洲精品久久| 亚洲国产日韩欧美精品在线观看 | 999久久久国产精品视频| 18禁美女被吸乳视频| 人人澡人人妻人| 亚洲国产精品sss在线观看| 国产成人精品久久二区二区91| 中文字幕精品亚洲无线码一区 | 一级毛片女人18水好多| 久久精品91蜜桃| 免费av毛片视频| 一本一本综合久久| xxx96com| 亚洲 欧美 日韩 在线 免费| 草草在线视频免费看| 大型av网站在线播放| 性欧美人与动物交配| 9191精品国产免费久久| 午夜视频精品福利| 久久久水蜜桃国产精品网| 午夜精品在线福利| 中文字幕最新亚洲高清| 又紧又爽又黄一区二区| 久久99热这里只有精品18| 精品免费久久久久久久清纯| 亚洲成a人片在线一区二区| 可以在线观看毛片的网站| 丝袜在线中文字幕| 黑丝袜美女国产一区| 成年女人毛片免费观看观看9| 日韩中文字幕欧美一区二区| 久久久久久大精品| 精品无人区乱码1区二区| 国产99久久九九免费精品| 亚洲精品国产区一区二| 精品乱码久久久久久99久播| 亚洲五月天丁香| 午夜久久久在线观看| 欧美亚洲日本最大视频资源| 人成视频在线观看免费观看| 亚洲国产欧美网| 亚洲成人国产一区在线观看| www国产在线视频色| 亚洲精品国产区一区二| 国产精品 国内视频| 人人妻人人看人人澡| 日韩欧美三级三区| 成熟少妇高潮喷水视频| 神马国产精品三级电影在线观看 | 少妇裸体淫交视频免费看高清 | 国产精品久久久久久人妻精品电影| 亚洲av中文字字幕乱码综合 | 成人三级做爰电影| 亚洲av成人不卡在线观看播放网| 高清毛片免费观看视频网站| 国产精品乱码一区二三区的特点| 19禁男女啪啪无遮挡网站| 欧美一级a爱片免费观看看 | 身体一侧抽搐| 免费搜索国产男女视频| 免费在线观看日本一区| 午夜a级毛片| 操出白浆在线播放| 国产精品亚洲一级av第二区| 久久中文字幕一级| 成人一区二区视频在线观看| 午夜福利一区二区在线看| 午夜激情福利司机影院| 身体一侧抽搐| 亚洲精品久久国产高清桃花| 老司机午夜十八禁免费视频| 久久精品91蜜桃| 中亚洲国语对白在线视频| 亚洲中文日韩欧美视频| 午夜免费观看网址| 久久草成人影院| 国产精品影院久久| 欧美日韩一级在线毛片| 动漫黄色视频在线观看| 国产成人精品久久二区二区91| 色播在线永久视频| 免费高清视频大片| 伦理电影免费视频| 中文字幕人成人乱码亚洲影| av片东京热男人的天堂| 日韩欧美 国产精品| 日韩高清综合在线| 精品福利观看| 精品国产乱码久久久久久男人| 亚洲 欧美一区二区三区| 国产精品免费一区二区三区在线| 亚洲avbb在线观看| 91成人精品电影| 99热6这里只有精品| 免费在线观看影片大全网站| 高清毛片免费观看视频网站| 天天一区二区日本电影三级| 国产精品九九99| 男女那种视频在线观看| 久久99热这里只有精品18| 美女扒开内裤让男人捅视频| 亚洲成人精品中文字幕电影| 精品一区二区三区四区五区乱码| 亚洲国产欧美日韩在线播放| 国产色视频综合| 国产成人欧美| 真人一进一出gif抽搐免费| 亚洲精品中文字幕在线视频| 大香蕉久久成人网| 亚洲精品在线美女| 1024手机看黄色片| 丝袜在线中文字幕| 欧美日韩福利视频一区二区| 国产精品综合久久久久久久免费| 久9热在线精品视频| 久久久久久免费高清国产稀缺| 久久九九热精品免费| 妹子高潮喷水视频| 1024手机看黄色片| 一区二区三区高清视频在线| 欧美日韩一级在线毛片| 男人舔女人的私密视频| 可以在线观看的亚洲视频| 亚洲专区中文字幕在线| 热99re8久久精品国产| 欧美+亚洲+日韩+国产| 黄色片一级片一级黄色片| 男人舔奶头视频| 亚洲成a人片在线一区二区| 黄网站色视频无遮挡免费观看| 日本免费a在线| 99精品在免费线老司机午夜| 久久精品夜夜夜夜夜久久蜜豆 | 少妇 在线观看| 天天添夜夜摸| 人妻丰满熟妇av一区二区三区| 搡老岳熟女国产| 18禁黄网站禁片午夜丰满| 黑人巨大精品欧美一区二区mp4| 日韩欧美免费精品| 可以在线观看毛片的网站| 午夜免费观看网址| 99久久综合精品五月天人人| 午夜福利高清视频| 麻豆av在线久日| 国产精品乱码一区二三区的特点| 色综合站精品国产|