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

    Duffing振子變尺度未知多頻微弱信號(hào)提取方法研究

    2019-06-06 08:28:06黃繼堯陳長(zhǎng)興凌云飛藺向陽
    儀表技術(shù)與傳感器 2019年5期
    關(guān)鍵詞:振子尺度軌跡

    黃繼堯,陳長(zhǎng)興,凌云飛,藺向陽

    (空軍工程大學(xué)基礎(chǔ)部,陜西西安 710043)

    0 引言

    Duffing振子是一種用于信號(hào)處理的非線性混沌振子模型,它在混沌臨界狀態(tài)下對(duì)噪聲具有免疫性,但對(duì)初始條件極為敏感[1]。可以利用Duffing振子的這些特性來檢測(cè)強(qiáng)噪聲背景下的微弱信號(hào)。自20世紀(jì)后期引入混沌振子系統(tǒng)以檢測(cè)微弱信號(hào)[2-4]以來,諸多學(xué)者在此領(lǐng)域進(jìn)行了深入的研究。謝濤等證明了混沌振子檢測(cè)微弱信號(hào)的可靠性,為混沌檢測(cè)微弱信號(hào)提供了理論論據(jù)[5]。

    李月等分析了不同類型噪聲對(duì)混沌系統(tǒng)的影響,證明了混沌系統(tǒng)對(duì)噪聲的免疫性[6];秦衛(wèi)陽等分析了Duffing方程的派生系統(tǒng)在同步過程中能夠?qū)崿F(xiàn)檢測(cè)2種信號(hào)中的微小差別[7];劉海波等提出一種利用單一振子實(shí)現(xiàn)以Duffing振子大周期態(tài)為基礎(chǔ)的測(cè)量頻率的方法[8]。常見的混沌狀態(tài)判定方法有Lyapunov指數(shù)法、龐加萊截面法、分維數(shù)計(jì)算法、Kolmogorov熵法、Melnikov法等[9]。與傳統(tǒng)的線性信號(hào)檢測(cè)方法相比,基于Duffing振子的微弱信號(hào)檢測(cè)方法可以實(shí)現(xiàn)非常低的信噪比(SNR),信噪比門限最低可達(dá)到-100 dB,并可得到傳統(tǒng)檢測(cè)方法無法得到的微弱待測(cè)信號(hào),這為低信噪比的微弱信號(hào)檢測(cè)提供了一種新的方法[10]。

    然而目前大多數(shù)基于Duffing振子的微弱信號(hào)檢測(cè)研究都是針對(duì)單頻諧波信號(hào)的檢測(cè)識(shí)別[11-14],而很少有研究關(guān)注微弱多頻信號(hào)的檢測(cè)。因此對(duì)基于Duffing振子的微弱多頻信號(hào)檢測(cè)進(jìn)行深入研究是十分有必要的。本文提出了一種結(jié)合信號(hào)頻率截取預(yù)處理和Duffing振子檢測(cè)微弱多頻信號(hào)的方法。

    1 Duffing振子微弱信號(hào)檢測(cè)原理

    Holmes型Duffing方程[15]一般形式為

    (1)

    式中:k為阻尼比;F為周期性策動(dòng)力幅值;ω為角頻率;Fcosωt為周期性策動(dòng)力。

    當(dāng)k保持不變時(shí),隨著F的增大,存在一個(gè)臨界策動(dòng)力振幅Fd。當(dāng)F小于Fd時(shí),系統(tǒng)處于混沌狀態(tài);當(dāng)F大于Fd時(shí),系統(tǒng)處于大周期態(tài)。在Duffing振子方程式(1)中,臨界值Fd≈0.825?;煦绾痛笾芷趹B(tài)的變化及策動(dòng)力臨界值如圖1所示。

    (a)混沌態(tài)

    (b)周期態(tài)

    (c)臨界點(diǎn)圖1 系統(tǒng)歷經(jīng)的狀態(tài)變化及分岔臨界圖

    當(dāng)Duffing系統(tǒng)處于從混沌狀態(tài)轉(zhuǎn)變?yōu)榇笾芷趹B(tài)的臨界狀態(tài)時(shí),就統(tǒng)計(jì)特性而言,它不受高斯白噪聲的影響。由于對(duì)參數(shù) 的敏感性和對(duì)噪聲的免疫性,式(1)可以進(jìn)一步表示為

    (2)

    式中:s(t)為待測(cè)信號(hào);n(t)為噪聲;s(t)+n(t)為帶有噪聲的待測(cè)信號(hào)。

    由式(2)可知,當(dāng)式(1)中的F略小于臨界幅值Fd時(shí),s(t)+n(t)輸入Duffing振子系統(tǒng),n(t)作為擾動(dòng)項(xiàng),則可以使用式(2)檢測(cè)待測(cè)信號(hào)。sn(t)=s(t)+n(t)是含噪待測(cè)信號(hào)。當(dāng)s(t)包含角頻率為ω的信號(hào)分量時(shí),系統(tǒng)將從混沌狀態(tài)轉(zhuǎn)換為大周期態(tài),這可以從系統(tǒng)響應(yīng)的相平面軌跡獲得,然后便可以獲得待測(cè)信號(hào)的頻率。

    傅里葉變換是一種分析信號(hào)的方法,可以把具有周期特征的非正弦信號(hào)分解為多個(gè)疊加的正弦信號(hào),然后通過對(duì)時(shí)間連續(xù)信號(hào)的等間隔采樣,得到采樣頻率fs,并把采樣值轉(zhuǎn)換為數(shù)字序列。連續(xù)時(shí)間序列x(t)與其對(duì)應(yīng)的連續(xù)非周期頻譜函數(shù)X(f)之間的轉(zhuǎn)換關(guān)系為:

    (3)

    當(dāng)設(shè)定一個(gè)待測(cè)信號(hào)為hcos(ωt+φ),其采樣頻率為fs。ω為確定的角頻率,φ是初始相位。假設(shè)采樣信號(hào)向坐標(biāo)軸左側(cè)移動(dòng)dt=1/fs。根據(jù)公式T/dt=2π/dθ,就可以得到dθ=dt·2π/T=ω/fs,T為信號(hào)的周期,dθ為對(duì)應(yīng)相位的增量。因此,將采樣信號(hào)向左挪動(dòng)時(shí)間j等同于將待測(cè)信號(hào)由hcos(ωt+φ)變成hcos(ωt+φ+jdθ)。

    在不考慮噪聲,即n(t)=0的情況下,設(shè)定待測(cè)信號(hào)尺度變換后(尺度變換系數(shù)R=ω)與策動(dòng)力有相同的角頻率,則檢測(cè)待測(cè)信號(hào)的Duffing方程為

    (4)

    改變式(6)右側(cè)項(xiàng)的形式

    (5)

    從式(4)中,可以推出只有當(dāng)(φ+jdθ)=2mπ(m是整數(shù))時(shí),與待測(cè)信號(hào)和策動(dòng)力相加的待測(cè)信號(hào)的振幅是最大的。

    文獻(xiàn)[16]提出了一種基于Duffing振子的尺度變換微弱信號(hào)檢測(cè)模型,有效地解決了這些問題。檢測(cè)模型如式(6)所示:

    (6)

    這4個(gè)方程的輸入信號(hào)初始相位的檢測(cè)范圍能夠覆蓋整段頻率(-π,π),避免了單一方程只能檢測(cè)某一段初始相位的限制。

    式(6)解決了基于Duffing振子微弱信號(hào)檢測(cè)的一些關(guān)鍵問題,然而目前僅對(duì)含噪聲的單頻諧波信號(hào)進(jìn)行分析,而不涉及實(shí)際工程應(yīng)用中常見的微弱多頻信號(hào)的檢測(cè)。因此本文將研究基于Duffing振子的微弱多頻信號(hào)檢測(cè)。

    2 基于頻率截獲Duffing振子的微弱多頻信號(hào)檢測(cè)

    2.1 信號(hào)頻率截取方法

    為了將式(6)中的檢測(cè)模型應(yīng)用于微弱多頻信號(hào),提出了一種信號(hào)頻率截取的預(yù)處理方法。其原理如下:根據(jù)Shannon采樣定理,試驗(yàn)信號(hào)sn(t)的采樣頻率為fs,則其截止頻率為fs/2,為了便于后續(xù)分析,首先在頻域中截取待測(cè)信號(hào)。截取頻率的下限為fmin,上限為fmax,其滿足0

    首先,利用快速傅里葉變換(FFT)將試驗(yàn)信號(hào)sn(t)變換為頻域信號(hào)sn(f);其次,信號(hào)sn′(f)是由sn(f)通過頻率截取得到的,即離散譜數(shù)據(jù)的頻率范圍為(fmin,fmax);最后,將信號(hào)sn′(f)通過快速傅里葉逆變換(IFFT)轉(zhuǎn)變?yōu)闀r(shí)域信號(hào)sn′(t)。

    上述頻率截獲過程可以表示為:

    (7)

    (8)

    (9)

    頻率截取過程后,sn′(t)只包含滿足fmin

    2.2 待測(cè)信號(hào)提取模型

    將信號(hào)頻率截獲方法與式(6)檢測(cè)模型相結(jié)合,針對(duì)微弱多頻信號(hào),提出了一種新的待測(cè)信號(hào)提取模型。其過程具體如下:

    (1)對(duì)原始待測(cè)信號(hào)進(jìn)行加噪處理,得到被噪聲污染的待測(cè)信號(hào)sn(t),采樣頻率為fs, 由強(qiáng)度為D的背景噪聲和微弱待測(cè)信號(hào)s1(t),s2(t),…,sn(t)構(gòu)成,微弱待測(cè)信號(hào)的角頻率和振幅分別為ω1,ω2,…,ωn和h1,h2,…,hn。

    (2)確定原始待測(cè)信號(hào)的頻率測(cè)試范圍[ωmin,ωmax]。然后設(shè)定ω=ωmin,接著設(shè)置頻率增量dω,令ωi+1=ωi+dω。

    (3)然后對(duì)待測(cè)信號(hào)sn(t)進(jìn)行傅里葉變換后得到信號(hào)sn(f)再進(jìn)行頻段為[ωmin,ωmax]的頻率截取預(yù)處理。新的信號(hào)在預(yù)處理后為sn′(f) ;

    (4)設(shè)定尺度系數(shù)Ri=ωi(即ω的值每增長(zhǎng)一次,R的值都會(huì)相應(yīng)的增長(zhǎng)并與之相等),然后對(duì)信號(hào)進(jìn)行頻率/時(shí)間尺度變換,之后得到在時(shí)間尺度t′下新的待測(cè)信號(hào)sn′(t)。變換后的采樣頻率為fsr=fs/Ri;

    (5)頻率截取預(yù)處理和尺度變換后,將新獲得的信號(hào)sn′(t)輸入到式(6)中。然后利用四階Runge-Kutta算法求解各方程,繪制各輸出響應(yīng)的相平面軌跡圖。只要任何一個(gè)相平面軌跡處于大周期態(tài),則意味著待測(cè)信號(hào)sn(t)中的頻率分量ωk被檢測(cè)到,然后重復(fù)迭代過程,直到頻率滿足ω≥ωmax為止,過程如圖2所示。

    圖2 信號(hào)檢測(cè)框圖

    從而檢測(cè)到該頻率范圍內(nèi)的所有信號(hào)頻率分量,實(shí)現(xiàn)對(duì)未知多頻率信號(hào)的檢測(cè)。如果dω太小,檢測(cè)未知信號(hào)的計(jì)算量將大幅增加。選擇角頻率為ω=20 rad/s的待測(cè)信號(hào)進(jìn)行測(cè)試,設(shè)定Δω為0.2 rad/s。選擇不同的Δω的倍數(shù)與計(jì)算時(shí)間的關(guān)系如圖3所示,選擇dω為2Δω有較快的運(yùn)算速度,因此選擇dω為2Δω。

    圖3 不同的Δω倍數(shù)與計(jì)算時(shí)間的關(guān)系圖

    所以當(dāng)檢測(cè)過程完成后,可以覆蓋完整的頻段[ωmin,ωmax]。迭代過程如圖4所示:

    圖4 未知多頻信號(hào)檢測(cè)的迭代過程

    3 仿真與分析

    設(shè)定待測(cè)信號(hào)的方程為

    (10)

    式中:ω1,ω2,ω3和h1,h2,h3分別為sn(t)中包含的3個(gè)頻率分量的角頻率和振幅;g(t)為均值為0,方差為1的高斯白噪聲。

    設(shè)定待測(cè)信號(hào)的初始相位為0,待測(cè)信號(hào)的采樣頻率fs=100 Hz,h1=h2=h3=0.1,ω1、ω2、ω3分別為20,30,50 rad/s,噪聲強(qiáng)度D=1,Δω為0.2 rad/s。設(shè)定[ωmin,ωmax]為(0,60)。繪制待測(cè)信號(hào)的時(shí)域圖和頻域圖如圖5所示,波形顯示3 000個(gè)采樣點(diǎn)。

    由于在實(shí)際工程問題中噪聲一般是不可忽略的,因此在時(shí)域圖和頻域圖中,可以輕易發(fā)現(xiàn)由于噪聲的污染,很難找到待測(cè)信號(hào)的頻率的存在。而相圖的檢測(cè)也同樣容易受到噪聲的干擾,如圖6所示。

    將ω=20 rad/s的待測(cè)信號(hào)送入式(6)中的第一個(gè)方程,可以發(fā)現(xiàn)圖6(a)所顯示的相圖在噪聲的影響下,較難分辨是出于混沌狀態(tài)還是大周期態(tài)。而圖6(b)所顯示的則是在經(jīng)過尺度變換后消除了噪聲的影響,相圖是處于大周期態(tài)的。在下文中,都將以經(jīng)過尺度變換后的相圖闡述。

    采用本文提出的優(yōu)化模型進(jìn)行檢測(cè),在ω=10 rad/s,20 rad/s,30 rad/s,40 rad/s,50 rad/s,60 rad/s,6個(gè)點(diǎn)處進(jìn)行取值,如圖7所示。

    (a)時(shí)域圖

    (a)未經(jīng)過尺度變換的相平面軌跡

    (b)經(jīng)過尺度變換的相平面軌跡 圖6 ω=20 rad/s的待測(cè)信號(hào)的相平面軌跡

    通過實(shí)驗(yàn)可以得知,圖7(b)、(c)、(e)所顯示的是分別為ω=20 rad/s、ω=30 rad/s和ω=50 rad/s的頻率分量所對(duì)應(yīng)的相圖,相圖是明顯處于大周期態(tài)的。而圖7(a)、(d)、(f)所顯示的分別是在ω=10 rad/s、ω=40 rad/s和ω=60 rad/s的頻點(diǎn)所對(duì)應(yīng)的相圖,相圖是處于混沌狀態(tài)的,符合實(shí)驗(yàn)設(shè)定的預(yù)期條件。

    實(shí)驗(yàn)發(fā)現(xiàn):

    (1)未使用尺度變換,就很難將噪聲污染的相圖分辨出來,存在對(duì)相圖的誤判可能,尺度變換可以有效地消除噪聲對(duì)相圖的干擾,尺度系數(shù)R的存在有利于待測(cè)信號(hào)的提取。

    (2)使用本方法成功檢測(cè)到了ω=20 rad/s,ω=30 rad/s,ω=50 rad/s的頻率分量。而在ω=60 rad/s這兩個(gè)頻點(diǎn)則相圖未發(fā)生狀態(tài)轉(zhuǎn)移,證明不存在該頻率的分量,具有普遍性和實(shí)用性。

    (a)ω=10 rad/s的相平面軌跡

    (b)ω=20 rad/s的相平面軌跡

    (c)ω=30 rad/s的相平面軌跡

    (d)ω=40 rad/s的相平面軌跡

    (e)ω=50 rad/s的相平面軌跡

    (f)ω=60 rad/s的相平面軌跡

    (3)本方法實(shí)現(xiàn)的微弱信號(hào)檢測(cè)最低信噪比門限可以達(dá)到-24.2 dB,低于文獻(xiàn)[16]中提出的原檢測(cè)模型最低信噪比門限的-13.01 dB。

    (4)Δω的選取如果過大,則會(huì)出現(xiàn)漏檢測(cè)現(xiàn)象;Δω的選取如果過小,則檢測(cè)效率就較低,因此需要找一個(gè)平衡點(diǎn)。一般采用比待測(cè)信號(hào)頻率低3個(gè)數(shù)量級(jí)的Δω可以實(shí)現(xiàn)限定頻率范圍內(nèi)全部頻率的檢測(cè)。本文所提出的方法在實(shí)際實(shí)驗(yàn)中采用的Δω是0.2 rad/s,檢測(cè)到了3個(gè)待測(cè)信號(hào)的頻率。

    (5)本文所提出的方法與Duffing振子并行處理方式的區(qū)別在于,本文所提出的方法是對(duì)頻域加帶通濾波器的基礎(chǔ)上進(jìn)行而且僅用一個(gè)Duffing振子便可檢測(cè)到限定頻域內(nèi)全部的待測(cè)信號(hào),而Duffing振子并行處理方式一般采用陣列的方式進(jìn)行檢測(cè),即需要較多振子。本文提出的方法節(jié)省了振子的數(shù)量,對(duì)運(yùn)算條件需求較低,在實(shí)際實(shí)驗(yàn)處理中也可以達(dá)到較為理想的效果。

    4 結(jié)束語

    本文研究了一種未知多頻微弱信號(hào)檢測(cè)方法,該方法針對(duì)原檢測(cè)模型只能同時(shí)檢測(cè)某一頻率的待測(cè)信號(hào),克服了噪聲誘導(dǎo)對(duì)檢測(cè)結(jié)果的不利影響,提出了一種可同時(shí)檢測(cè)多個(gè)未知頻率微弱信號(hào)的Duffing振子優(yōu)化模型,可以實(shí)現(xiàn)提取未知多頻微弱信號(hào)各信號(hào)分量的頻率參數(shù),具有實(shí)用性和普遍性。仿真結(jié)果表明:該方法為微弱多頻信號(hào)處理提供了一種解決方案,拓寬了Duffing振子在微弱信號(hào)檢測(cè)領(lǐng)域的應(yīng)用。

    猜你喜歡
    振子尺度軌跡
    彈簧振子問題的分析與求解
    財(cái)產(chǎn)的五大尺度和五重應(yīng)對(duì)
    軌跡
    軌跡
    軌跡
    非線性Duffing擾動(dòng)振子共振機(jī)制的研究
    進(jìn)化的軌跡(一)——進(jìn)化,無盡的適應(yīng)
    宇宙的尺度
    太空探索(2016年5期)2016-07-12 15:17:55
    基于近似熵和混沌振子的電力諧波檢測(cè)與估計(jì)
    電磁彈簧振子實(shí)驗(yàn)裝置的改進(jìn)
    一个人看的www免费观看视频| 亚洲av免费高清在线观看| a级毛片免费高清观看在线播放| 亚洲精品影视一区二区三区av| 禁无遮挡网站| 搡老乐熟女国产| 在线a可以看的网站| 国产人妻一区二区三区在| 男女国产视频网站| 观看免费一级毛片| 日韩三级伦理在线观看| 纵有疾风起免费观看全集完整版| 亚洲色图av天堂| 看免费成人av毛片| 久久久久久久国产电影| 成人漫画全彩无遮挡| 激情五月婷婷亚洲| 日韩不卡一区二区三区视频在线| 免费观看a级毛片全部| 午夜免费男女啪啪视频观看| 如何舔出高潮| 青春草国产在线视频| 亚洲欧美日韩卡通动漫| 国产成人a∨麻豆精品| 一级毛片我不卡| 丝袜脚勾引网站| 日本黄色片子视频| 亚洲av在线观看美女高潮| 亚洲av男天堂| 国产av码专区亚洲av| 国产白丝娇喘喷水9色精品| 久久久久久久久久久免费av| 亚洲电影在线观看av| 亚洲av国产av综合av卡| 成人国产麻豆网| 啦啦啦在线观看免费高清www| 搡老乐熟女国产| 亚洲国产最新在线播放| 秋霞在线观看毛片| 熟妇人妻不卡中文字幕| 欧美97在线视频| 国产成人一区二区在线| 欧美区成人在线视频| 寂寞人妻少妇视频99o| 精华霜和精华液先用哪个| 大码成人一级视频| 啦啦啦中文免费视频观看日本| 国产淫片久久久久久久久| 欧美 日韩 精品 国产| 中文在线观看免费www的网站| 亚洲av二区三区四区| 免费播放大片免费观看视频在线观看| 好男人视频免费观看在线| 国产白丝娇喘喷水9色精品| 26uuu在线亚洲综合色| 亚洲,一卡二卡三卡| 国产v大片淫在线免费观看| 伊人久久精品亚洲午夜| 亚洲怡红院男人天堂| 一二三四中文在线观看免费高清| 欧美成人午夜免费资源| 日本免费在线观看一区| 91久久精品国产一区二区三区| 免费黄色在线免费观看| 99re6热这里在线精品视频| 哪个播放器可以免费观看大片| 超碰av人人做人人爽久久| 天堂俺去俺来也www色官网| 欧美老熟妇乱子伦牲交| 免费播放大片免费观看视频在线观看| 久久6这里有精品| 亚洲成人久久爱视频| 热re99久久精品国产66热6| 国产乱来视频区| 免费黄网站久久成人精品| 久热这里只有精品99| 日韩国内少妇激情av| 欧美区成人在线视频| av播播在线观看一区| 国产成年人精品一区二区| 色视频在线一区二区三区| 性插视频无遮挡在线免费观看| 日本午夜av视频| 欧美3d第一页| 亚洲综合精品二区| 久久久久性生活片| 精品99又大又爽又粗少妇毛片| 亚洲av在线观看美女高潮| 下体分泌物呈黄色| 免费看av在线观看网站| 欧美变态另类bdsm刘玥| 三级国产精品欧美在线观看| 欧美日韩综合久久久久久| 久久99蜜桃精品久久| xxx大片免费视频| 一区二区三区乱码不卡18| 欧美成人精品欧美一级黄| 久久精品熟女亚洲av麻豆精品| 99视频精品全部免费 在线| 天天躁日日操中文字幕| 欧美激情在线99| 男人狂女人下面高潮的视频| 久久久久久久久久成人| 只有这里有精品99| 91久久精品国产一区二区三区| 国产淫片久久久久久久久| 国产一区有黄有色的免费视频| 亚洲成人av在线免费| 亚洲av成人精品一区久久| 日韩在线高清观看一区二区三区| 在线天堂最新版资源| 久久99热这里只有精品18| 欧美成人a在线观看| 国产精品一二三区在线看| 大又大粗又爽又黄少妇毛片口| h日本视频在线播放| 最新中文字幕久久久久| 777米奇影视久久| 久久99热这里只频精品6学生| 一本久久精品| 国产精品久久久久久精品古装| 成人亚洲精品av一区二区| 激情五月婷婷亚洲| 亚洲av一区综合| 欧美日韩国产mv在线观看视频 | av免费观看日本| 日韩av免费高清视频| 国产精品久久久久久av不卡| 国产乱来视频区| 男女无遮挡免费网站观看| 成人亚洲欧美一区二区av| 91aial.com中文字幕在线观看| 国产亚洲91精品色在线| 少妇猛男粗大的猛烈进出视频 | 人体艺术视频欧美日本| 在现免费观看毛片| 性色av一级| 日韩欧美一区视频在线观看 | 国产精品国产三级国产av玫瑰| 精品午夜福利在线看| 国产亚洲91精品色在线| 国产69精品久久久久777片| 免费大片黄手机在线观看| 3wmmmm亚洲av在线观看| 国产精品久久久久久精品电影小说 | 2021少妇久久久久久久久久久| 一级毛片 在线播放| 91精品伊人久久大香线蕉| 又黄又爽又刺激的免费视频.| av网站免费在线观看视频| 久久影院123| 蜜桃亚洲精品一区二区三区| 亚洲欧美一区二区三区黑人 | 亚洲精品乱久久久久久| 亚洲精品中文字幕在线视频 | 国产精品麻豆人妻色哟哟久久| 日韩一本色道免费dvd| 超碰av人人做人人爽久久| 嘟嘟电影网在线观看| 三级男女做爰猛烈吃奶摸视频| 一边亲一边摸免费视频| 国产精品久久久久久精品古装| 少妇高潮的动态图| 午夜日本视频在线| 久久精品国产a三级三级三级| 日日啪夜夜爽| 高清视频免费观看一区二区| 黄片wwwwww| 插阴视频在线观看视频| 国产精品国产av在线观看| 一区二区三区精品91| 六月丁香七月| 午夜免费男女啪啪视频观看| 我的女老师完整版在线观看| 蜜桃亚洲精品一区二区三区| 永久免费av网站大全| 亚洲av一区综合| 在线观看人妻少妇| 亚洲欧美精品专区久久| 老司机影院成人| 99re6热这里在线精品视频| 99久久中文字幕三级久久日本| 欧美xxxx性猛交bbbb| 黄色配什么色好看| 久久久久久国产a免费观看| 国产精品蜜桃在线观看| 日韩人妻高清精品专区| 欧美成人a在线观看| 久久这里有精品视频免费| 国产黄频视频在线观看| 一本一本综合久久| 九草在线视频观看| 日本猛色少妇xxxxx猛交久久| 欧美成人精品欧美一级黄| 国产成人精品久久久久久| 美女xxoo啪啪120秒动态图| 97在线人人人人妻| 男的添女的下面高潮视频| 少妇的逼好多水| 观看美女的网站| 日韩伦理黄色片| 一级片'在线观看视频| 亚洲,一卡二卡三卡| 晚上一个人看的免费电影| 久久亚洲国产成人精品v| 日产精品乱码卡一卡2卡三| 久久久久久久午夜电影| 午夜激情福利司机影院| 国产欧美日韩精品一区二区| 1000部很黄的大片| 国产亚洲91精品色在线| 在线免费十八禁| 成年女人在线观看亚洲视频 | 亚洲精品日韩av片在线观看| 观看美女的网站| 99热全是精品| 男女那种视频在线观看| 免费在线观看成人毛片| 精品一区二区三卡| 身体一侧抽搐| 久久影院123| 免费在线观看成人毛片| 色5月婷婷丁香| 日韩av在线免费看完整版不卡| 一级二级三级毛片免费看| 亚洲成人一二三区av| 不卡视频在线观看欧美| 午夜福利高清视频| 亚洲精品国产成人久久av| 欧美xxxx性猛交bbbb| 日本一本二区三区精品| 成年女人看的毛片在线观看| 美女国产视频在线观看| 超碰av人人做人人爽久久| 在线免费十八禁| 免费黄频网站在线观看国产| 亚洲国产日韩一区二区| 国产成人精品久久久久久| 可以在线观看毛片的网站| 国产探花极品一区二区| 亚洲国产高清在线一区二区三| 亚洲国产欧美人成| 插阴视频在线观看视频| 欧美97在线视频| 久久久亚洲精品成人影院| av在线观看视频网站免费| 黄色怎么调成土黄色| 欧美成人午夜免费资源| 国产精品一区二区在线观看99| av网站免费在线观看视频| 免费电影在线观看免费观看| 国产成人aa在线观看| kizo精华| 中文欧美无线码| 国产精品久久久久久av不卡| 真实男女啪啪啪动态图| 色综合色国产| 各种免费的搞黄视频| 性插视频无遮挡在线免费观看| 直男gayav资源| 国产精品.久久久| 国产真实伦视频高清在线观看| 午夜亚洲福利在线播放| 少妇熟女欧美另类| 精品酒店卫生间| 精华霜和精华液先用哪个| 一区二区三区四区激情视频| 国内精品宾馆在线| 成人午夜精彩视频在线观看| 亚洲色图av天堂| 搡女人真爽免费视频火全软件| 美女xxoo啪啪120秒动态图| 久久99热6这里只有精品| 成人免费观看视频高清| 亚洲欧洲日产国产| 午夜福利视频1000在线观看| 日韩在线高清观看一区二区三区| 国产午夜精品久久久久久一区二区三区| 国产精品国产三级专区第一集| 色吧在线观看| 国产人妻一区二区三区在| 欧美成人精品欧美一级黄| 久久久久久久久久人人人人人人| 一区二区三区精品91| 日韩制服骚丝袜av| 九色成人免费人妻av| videossex国产| 晚上一个人看的免费电影| 亚洲精品一二三| 一区二区三区精品91| 国产精品久久久久久av不卡| 一级爰片在线观看| 热re99久久精品国产66热6| 午夜免费观看性视频| 日韩欧美精品v在线| 国国产精品蜜臀av免费| 51国产日韩欧美| 欧美日本视频| 欧美日韩综合久久久久久| 国产精品国产av在线观看| 欧美精品一区二区大全| 一级爰片在线观看| 少妇的逼水好多| av在线观看视频网站免费| 精品久久久精品久久久| 涩涩av久久男人的天堂| 身体一侧抽搐| 一本色道久久久久久精品综合| 日本爱情动作片www.在线观看| 国产欧美日韩一区二区三区在线 | 久久99蜜桃精品久久| 各种免费的搞黄视频| 91精品国产九色| 网址你懂的国产日韩在线| 亚洲国产欧美人成| 涩涩av久久男人的天堂| 午夜日本视频在线| 天美传媒精品一区二区| 高清在线视频一区二区三区| 日韩一本色道免费dvd| 国产免费一区二区三区四区乱码| 亚洲精品一二三| 亚洲激情五月婷婷啪啪| 亚洲怡红院男人天堂| 色播亚洲综合网| 久热这里只有精品99| 亚洲国产av新网站| 免费黄色在线免费观看| 精品国产一区二区三区久久久樱花 | 国产 精品1| 成人亚洲精品一区在线观看 | 免费高清在线观看视频在线观看| 欧美成人a在线观看| 精品一区二区三区视频在线| 亚洲欧美成人精品一区二区| 免费高清在线观看视频在线观看| 国产精品秋霞免费鲁丝片| 国产精品熟女久久久久浪| 色视频www国产| av在线亚洲专区| av免费在线看不卡| 99热6这里只有精品| 波野结衣二区三区在线| 精品久久久久久久久av| 亚洲精品国产av蜜桃| 日本黄大片高清| 青青草视频在线视频观看| 真实男女啪啪啪动态图| 日韩欧美精品v在线| 毛片女人毛片| 国产伦理片在线播放av一区| 色吧在线观看| 亚洲成人一二三区av| 天堂中文最新版在线下载 | 亚洲国产成人一精品久久久| 午夜福利高清视频| 在线免费观看不下载黄p国产| 男的添女的下面高潮视频| 亚洲欧美中文字幕日韩二区| 男的添女的下面高潮视频| 麻豆成人午夜福利视频| 日韩一本色道免费dvd| 日韩 亚洲 欧美在线| 午夜福利高清视频| 亚洲精品第二区| 亚洲精品自拍成人| 国产亚洲一区二区精品| 欧美日韩精品成人综合77777| 欧美成人精品欧美一级黄| av免费观看日本| 女的被弄到高潮叫床怎么办| 在线观看一区二区三区| 又爽又黄a免费视频| 亚洲国产精品专区欧美| 成人美女网站在线观看视频| 日韩av免费高清视频| av专区在线播放| 色哟哟·www| 亚洲丝袜综合中文字幕| 嫩草影院精品99| 91狼人影院| 亚洲av一区综合| 天天一区二区日本电影三级| 欧美成人a在线观看| 在线a可以看的网站| 欧美成人a在线观看| 亚洲成人av在线免费| 日韩免费高清中文字幕av| 日本免费在线观看一区| 亚洲欧美中文字幕日韩二区| 汤姆久久久久久久影院中文字幕| 欧美丝袜亚洲另类| 久久99蜜桃精品久久| 美女视频免费永久观看网站| 国产伦精品一区二区三区视频9| 成年av动漫网址| 久久久成人免费电影| 最近最新中文字幕大全电影3| 五月玫瑰六月丁香| 看黄色毛片网站| 美女xxoo啪啪120秒动态图| 内射极品少妇av片p| 亚洲色图综合在线观看| 在线a可以看的网站| 99久久精品国产国产毛片| 人妻夜夜爽99麻豆av| 国产黄色视频一区二区在线观看| 亚洲自拍偷在线| 中文字幕制服av| 精品一区二区免费观看| 亚洲三级黄色毛片| 亚洲av电影在线观看一区二区三区 | 亚洲一级一片aⅴ在线观看| 最近最新中文字幕大全电影3| xxx大片免费视频| 午夜免费观看性视频| 嫩草影院精品99| 国产探花极品一区二区| 欧美xxxx黑人xx丫x性爽| 亚洲电影在线观看av| 国产精品熟女久久久久浪| 国产爱豆传媒在线观看| 亚洲美女视频黄频| 精品亚洲乱码少妇综合久久| 精品午夜福利在线看| 80岁老熟妇乱子伦牲交| 黄色怎么调成土黄色| 白带黄色成豆腐渣| 国产在视频线精品| 在线观看av片永久免费下载| 婷婷色av中文字幕| 国产精品秋霞免费鲁丝片| 国产精品麻豆人妻色哟哟久久| 国产男女内射视频| 中国美白少妇内射xxxbb| 一区二区三区四区激情视频| 国产有黄有色有爽视频| av女优亚洲男人天堂| 欧美亚洲 丝袜 人妻 在线| 狂野欧美白嫩少妇大欣赏| 国产熟女欧美一区二区| 综合色av麻豆| 丝袜喷水一区| eeuss影院久久| 天堂俺去俺来也www色官网| 久久久久久九九精品二区国产| 男插女下体视频免费在线播放| 三级经典国产精品| 视频中文字幕在线观看| 赤兔流量卡办理| 国产男人的电影天堂91| 最近最新中文字幕大全电影3| 最新中文字幕久久久久| 国产精品伦人一区二区| 免费看a级黄色片| 少妇 在线观看| 久久久久久九九精品二区国产| 亚洲av免费高清在线观看| 精品国产露脸久久av麻豆| 精品熟女少妇av免费看| 亚洲av免费高清在线观看| 午夜日本视频在线| 成年av动漫网址| 丰满人妻一区二区三区视频av| av黄色大香蕉| 国产伦精品一区二区三区视频9| 久久久久国产精品人妻一区二区| 国产男女超爽视频在线观看| 深爱激情五月婷婷| 久久亚洲国产成人精品v| 久久久久久伊人网av| 国产精品嫩草影院av在线观看| 夜夜看夜夜爽夜夜摸| 在线观看国产h片| 美女高潮的动态| 久久99蜜桃精品久久| 婷婷色综合www| www.av在线官网国产| 中国国产av一级| 久久这里有精品视频免费| 中国三级夫妇交换| 国产精品久久久久久久久免| 日韩欧美精品v在线| 亚洲精品成人av观看孕妇| 亚洲av欧美aⅴ国产| 午夜爱爱视频在线播放| 我的老师免费观看完整版| 国产亚洲精品久久久com| 国产免费一级a男人的天堂| 99久国产av精品国产电影| 成人亚洲欧美一区二区av| 亚洲婷婷狠狠爱综合网| 免费观看av网站的网址| 亚洲国产最新在线播放| 欧美区成人在线视频| 欧美xxxx黑人xx丫x性爽| av福利片在线观看| 国语对白做爰xxxⅹ性视频网站| 国产av码专区亚洲av| 性色avwww在线观看| 亚洲av二区三区四区| 欧美xxxx黑人xx丫x性爽| 九九在线视频观看精品| 免费观看的影片在线观看| 三级男女做爰猛烈吃奶摸视频| 丝袜美腿在线中文| 国产白丝娇喘喷水9色精品| 国产高清国产精品国产三级 | 国产又色又爽无遮挡免| 亚洲av一区综合| 国产精品久久久久久精品古装| 三级经典国产精品| 国产综合精华液| 亚洲欧美成人精品一区二区| 美女xxoo啪啪120秒动态图| 91狼人影院| 国产视频首页在线观看| 国产成人午夜福利电影在线观看| 99九九线精品视频在线观看视频| 欧美日本视频| av在线播放精品| av天堂中文字幕网| 2021少妇久久久久久久久久久| 国产免费福利视频在线观看| 日本欧美国产在线视频| 纵有疾风起免费观看全集完整版| 狂野欧美激情性bbbbbb| 亚洲精品,欧美精品| 美女高潮的动态| 国产精品无大码| 18禁动态无遮挡网站| 99久久精品热视频| 国产一区有黄有色的免费视频| 有码 亚洲区| 亚洲av成人精品一区久久| 日日撸夜夜添| 99久久人妻综合| 亚洲国产精品专区欧美| 中文字幕人妻熟人妻熟丝袜美| 免费av毛片视频| 简卡轻食公司| 在线观看av片永久免费下载| 人妻少妇偷人精品九色| 三级国产精品欧美在线观看| 简卡轻食公司| 免费高清在线观看视频在线观看| 一级av片app| av免费观看日本| 国产精品蜜桃在线观看| 少妇人妻 视频| 涩涩av久久男人的天堂| 日韩不卡一区二区三区视频在线| 99热这里只有精品一区| 亚洲欧洲国产日韩| 尤物成人国产欧美一区二区三区| 亚洲av男天堂| 免费人成在线观看视频色| 亚洲欧洲日产国产| 久久这里有精品视频免费| 国产精品精品国产色婷婷| 午夜免费男女啪啪视频观看| 狂野欧美激情性xxxx在线观看| 特级一级黄色大片| 天天躁日日操中文字幕| 熟女电影av网| 国产成人aa在线观看| 久久99蜜桃精品久久| 又爽又黄无遮挡网站| 国产精品一及| 美女被艹到高潮喷水动态| 18禁在线播放成人免费| 久久精品国产亚洲av涩爱| 久久人人爽av亚洲精品天堂 | 97超碰精品成人国产| 综合色丁香网| 色网站视频免费| 美女内射精品一级片tv| 身体一侧抽搐| 免费看av在线观看网站| 99久久精品热视频| 久久这里有精品视频免费| 久久精品久久精品一区二区三区| 欧美区成人在线视频| 51国产日韩欧美| 久久精品久久久久久噜噜老黄| 黄色欧美视频在线观看| 免费大片18禁| 国产精品一区二区在线观看99| 人妻少妇偷人精品九色| 国产一区二区三区av在线| 伊人久久国产一区二区| 国产精品人妻久久久影院| 熟女av电影| 91精品伊人久久大香线蕉| 夜夜看夜夜爽夜夜摸| 亚洲精品日本国产第一区| 亚洲激情五月婷婷啪啪| 国产爽快片一区二区三区| 啦啦啦在线观看免费高清www| 人人妻人人爽人人添夜夜欢视频 | 久久这里有精品视频免费| 一级二级三级毛片免费看| 女人久久www免费人成看片| 精品久久久精品久久久| 久久ye,这里只有精品| 久久久久国产网址| 视频中文字幕在线观看| 国产高清不卡午夜福利| 在线免费十八禁| 久久国产乱子免费精品| 晚上一个人看的免费电影| 各种免费的搞黄视频| 欧美日本视频| 王馨瑶露胸无遮挡在线观看| av网站免费在线观看视频|