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

    基于Pisarenko 譜估計(jì)的GPS 多普勒頻偏快速捕獲

    2012-08-29 06:35:40楊曉波
    科技視界 2012年20期
    關(guān)鍵詞:譜估計(jì)虛警環(huán)路

    楊曉波

    (石家莊職業(yè)技術(shù)學(xué)院電子工程系 河北 石家莊 050081)

    0 引言

    提高GPS 接收機(jī)多普勒頻移的捕獲時(shí)間指標(biāo)具有重要的實(shí)際意義[1].GPS 接收機(jī)偽碼的捕獲是載波多普勒和碼相位的二維搜索,如圖1 所示.

    圖1 偽碼相位和載波多普勒二維搜索

    多普勒搜索步長(zhǎng)一般取2/3T[2],T 為積分時(shí)間,所以一般多普勒頻偏得估計(jì)精度為幾百赫茲, 這對(duì)于PLL 來(lái)講過(guò)于粗糙,并且PLL 對(duì)動(dòng)態(tài)非常敏感,所以一般接收機(jī)先由FLL 將多普勒頻差估計(jì)到一個(gè)比較小的范圍, 然后再切換到PLL 進(jìn)入跟蹤狀態(tài), 這樣就在多普勒的捕獲時(shí)間中增加了FLL 的收斂時(shí)間,同時(shí)捕獲虛警后的滯留時(shí)間也會(huì)加大,對(duì)此文獻(xiàn)[3]提出對(duì)解擴(kuò)后的信號(hào)利用FFT 來(lái)估計(jì)載波多普勒,由于FFT 頻率分辨率受到信號(hào)長(zhǎng)度的限制,該方法也不能取得很好的效果.例如, 在5MHz 的采樣率下要達(dá)到100Hz 的分辨率就需要10ms 數(shù)據(jù),50000 個(gè)采樣點(diǎn),這樣不但增加了捕獲時(shí)間同時(shí)還增大了計(jì)算量和存儲(chǔ)空間.文獻(xiàn)[4]提出了將多普勒搜索單元設(shè)窄一些,如250Hz、125Hz 提高多普勒頻率估計(jì)的分辨率.但該方法會(huì)加大搜索空間,非常耗時(shí).本文在研究解擴(kuò)后GPS 信號(hào)特點(diǎn)的基礎(chǔ)上提出利用Pisarenko 譜估計(jì)的方法利用很少的數(shù)據(jù)來(lái)精確估計(jì)多普勒以減少載波多普勒的平均捕獲時(shí)間.

    1 FFT 快速捕獲原理與多普勒平均捕獲時(shí)間

    1.1 多普勒FFT 快速捕獲原理

    圖1 為PN 碼FFT 快捕的中頻數(shù)字化結(jié)構(gòu), 由數(shù)字下變頻、FFT 相關(guān)、檢測(cè)判決、解調(diào)、FLL 環(huán)路、符號(hào)時(shí)鐘恢復(fù)與保護(hù)等幾部分組成. 其中FLL 環(huán)路的工作流程是從接收數(shù)據(jù)流中提取頻偏估計(jì)值, 經(jīng)環(huán)路濾波后去控制本地NCO 的頻率,使其跟蹤輸入中頻信號(hào)的頻率變化, 從而使下變頻的輸出始終為零中頻數(shù)字基帶信號(hào).

    接收機(jī)接收到的信號(hào)模型為:

    其中M:可見(jiàn)衛(wèi)星個(gè)數(shù);Di:電文;Pi:偽隨機(jī)碼;ωi:GPS 信號(hào)載波頻率加多普勒頻率;τi:信號(hào)延遲;φi:載波初始相位;n(t):雙邊功率譜密度為N0/2 加性高斯白噪聲.

    圖2 多普勒頻偏FFT 快捕的中頻數(shù)字化結(jié)構(gòu)

    FFT 相關(guān)器同相、正交輸出:

    其中,ε 為接收的偽隨機(jī)碼與本地偽隨機(jī)碼未完全對(duì)齊的基碼數(shù),ωd為載波的多普勒頻移,T 為相關(guān)積分時(shí)間,φ 為相位差,nI和nQ是相互獨(dú)立的基帶高斯白噪聲,相關(guān)輸出后計(jì)算同相、正交包絡(luò):

    一個(gè)頻點(diǎn)上包絡(luò)的最大值與恒虛警門(mén)限比較如果包絡(luò)大于門(mén)限,認(rèn)為捕獲成功,在沒(méi)有信號(hào)時(shí),包絡(luò)是服從瑞利分布的隨機(jī)過(guò)程[2],經(jīng)本文研究FFT 快捕恒虛警門(mén)限為:

    其中,n 為FFT 點(diǎn)數(shù),σ 為1-sigma 噪聲功率, 但對(duì)于多普勒來(lái)講這只是實(shí)現(xiàn)了粗捕, 然后將粗捕后的多普勒傳遞給FLL, 經(jīng)過(guò)一定時(shí)間后環(huán)路鎖定則系統(tǒng)轉(zhuǎn)入頻率跟蹤狀態(tài),即捕獲成功,否則將重新回到粗略估計(jì)狀態(tài).

    1.2 多普勒頻偏平均捕獲時(shí)間

    圖3 多普勒頻偏捕獲狀態(tài)轉(zhuǎn)移圖

    根據(jù)以上分析可知,多普勒頻移捕獲算法分為兩個(gè)步驟:粗略估計(jì)和捕獲驗(yàn)證.首先,把接收機(jī)的頻域搜索范圍劃分為若干相等的頻率區(qū)間,各區(qū)間的中心頻率點(diǎn)設(shè)為fn.粗略估計(jì)就是對(duì)頻點(diǎn)fn上的最大相關(guān)值與門(mén)限比較, 如果大于門(mén)限則將NCO 固定在此頻率點(diǎn)上.系統(tǒng)隨即進(jìn)入捕獲驗(yàn)證狀態(tài),FLL環(huán)路對(duì)殘留頻偏進(jìn)行捕獲,若經(jīng)過(guò)時(shí)間Ts(FLL 收斂時(shí)間)后環(huán)路鎖定(殘留頻偏小于某個(gè)規(guī)定值),則系統(tǒng)轉(zhuǎn)入頻率跟蹤狀態(tài)(即捕獲成功),否則將重新回到粗略估計(jì)狀態(tài).

    所以多普勒頻移的捕獲過(guò)程具有馬爾科夫性, 圖3 中狀態(tài)n(1,2,3,…)對(duì)應(yīng)系統(tǒng)的粗略估計(jì)狀態(tài),另外兩種情況分別是捕獲狀態(tài)和虛警狀態(tài),捕獲狀態(tài)是馬爾科夫過(guò)程的穩(wěn)態(tài),設(shè)系統(tǒng)初始時(shí)刻處于狀態(tài)1, 經(jīng)FFT 快捕后判斷進(jìn)入其它狀態(tài)的轉(zhuǎn)移路徑:a.進(jìn)入捕獲狀態(tài);b.進(jìn)入狀態(tài)2,即系統(tǒng)沒(méi)有滿(mǎn)足判決條件仍處于粗捕狀態(tài);c.進(jìn)入虛警狀態(tài),經(jīng)驗(yàn)證后再回到粗捕狀態(tài)如此循環(huán),直到系統(tǒng)進(jìn)入捕獲狀態(tài),該過(guò)程轉(zhuǎn)移圖中各支路的增益表達(dá)式:

    其中G 是頻率單元數(shù),Tc是一個(gè)頻點(diǎn)的駐留時(shí)間,Ts是虛警滯留時(shí)間,于是得到多普勒平均捕獲時(shí)間的生成函數(shù):

    式中PD,PFA,PM分別表示多普勒頻移捕獲的檢測(cè)概率、虛警概率和漏檢概率,由于PD+PFA+PM=1,H1(1)=PD,H2(1)=PFA,H3(1)=PM,H4(1)=1 平均捕獲時(shí)間計(jì)算得:

    其中G 為頻點(diǎn)的個(gè)數(shù),Ts是虛警滯留時(shí)間,也是多普勒虛警驗(yàn)證時(shí)間,對(duì)于有FLL 輔助的環(huán)路,就是FLL 的收斂時(shí)間,如果能夠?qū)⑦@個(gè)時(shí)間縮短,將縮短多普勒的平均捕獲時(shí)間.

    2 Pisarenko 估計(jì)多普勒頻偏具體算法

    由式(1)如果第j 個(gè)信號(hào)被捕獲,則解擴(kuò)后的信號(hào)模型為:

    dsj(t)=Dj(t-τj)sin(ωjt+φj)

    式(9)等號(hào)右邊第二項(xiàng)中Rij為第j 個(gè)偽隨機(jī)碼和第i 個(gè)偽隨機(jī)碼的互相關(guān), 研究表明互相關(guān)部分最大可能功率為-21dB[5],第三項(xiàng)為偽隨機(jī)碼和加性白噪聲的相關(guān).這兩部分可以看作高斯白噪聲, 第一項(xiàng)為載波部分, 其中載波相位φj是[-π,π ]之間均勻分布的隨機(jī)數(shù).所以,解擴(kuò)后的信號(hào)可以認(rèn)為是平穩(wěn)諧波過(guò)程與高斯白噪聲的和,這樣的觀測(cè)信號(hào),可以利用Pisarenko 譜估計(jì)對(duì)載波的多普勒頻率進(jìn)行高分辨率的估計(jì)[6]. Pisarenko 譜估計(jì)由于使用了奇異值分解和總體最小二乘方法,整個(gè)計(jì)算具有好的數(shù)值穩(wěn)定性,而且AR 階數(shù)和參數(shù)估計(jì)具有非常高的精確度[7].所以本文采用這種方法估計(jì)載波多普勒.

    當(dāng)諧波信號(hào)在加性白噪聲中被觀測(cè)時(shí), 它的AR 參數(shù)和MA 參數(shù)完全相同一個(gè)特殊的隨機(jī)過(guò)程, 這一特殊的ARMA過(guò)程所服從的法方程為:

    與修正的Yule-Walker 方程類(lèi)似, 該方程可以構(gòu)造超定的方程組,并使用SVD-TLS[6]算法求解.

    具體算法:

    步驟1 利用觀測(cè)數(shù)據(jù)的樣本自相關(guān)函數(shù)估計(jì)R?y(k)構(gòu)造法方程的擴(kuò)展階自相關(guān)矩陣

    式中pe>2p,并且M>p;

    步驟2 將矩陣Re當(dāng)作增廣矩陣,利用算法SVD-TLS 確定AR 的階數(shù)2p 和系數(shù)向量a 的總體最小二乘估計(jì);

    步驟3 計(jì)算特征多項(xiàng)式

    步驟4 利用

    計(jì)算載波的頻率.

    3 計(jì)算機(jī)仿真及分析

    3.1 載波多普勒估計(jì)精度

    載波頻率為1.25MHz,多普勒頻率設(shè)為300Hz,采樣頻率5MHz,若多普勒搜索步長(zhǎng)為500Hz 則粗捕后多普勒最大誤差250Hz,通過(guò)一個(gè)以捕獲頻率為中心頻率帶寬1K 的帶通濾波器,此時(shí)擴(kuò)頻增益為33dB,考慮實(shí)際的較弱GPS 信號(hào)信噪比-20dB,解擴(kuò)后信噪比為13dB,仿真利用1ms 數(shù)據(jù).

    由于估計(jì)的諧波頻率只有1 個(gè),所以式(10)中p=1,取M=12,pe=10,當(dāng)然也可以取其它值,只要滿(mǎn)足擴(kuò)展階自相關(guān)矩陣的要求即可.圖5 為信噪比13dB 時(shí)20 次仿真的譜估計(jì)圖,圖中實(shí)線為估計(jì)諧波功率譜密度,虛線所示為真實(shí)頻率,可以看出每次的估計(jì)結(jié)果非常穩(wěn)定.

    圖4 本文方法估計(jì)載波多普勒

    同時(shí)本文應(yīng)用三級(jí)最小二乘ARMA 參數(shù)模型法和FFT估計(jì)載波多普勒, 數(shù)據(jù)長(zhǎng)度1ms, 其它條件相同, 下表是Pisarenko、ARMA 參數(shù)模型法及FFT 三種估計(jì)方法的對(duì)比,表1 為仿真進(jìn)行50 次的頻率均方誤差(RMS),由于FFT 的頻率分辨率對(duì)于1ms 的數(shù)據(jù)長(zhǎng)度為1kHz, 其輸出頻率為1kHz 的整數(shù)倍,所以其誤差不隨信噪比變化.

    表1 1ms 數(shù)據(jù)頻率估計(jì)RMS(Hz)

    考慮偽碼未對(duì)齊情況. 碼片誤差最差情況0.5 個(gè)碼片,如圖5 所示,圖中實(shí)線表示本地碼和接收信號(hào)偽碼對(duì)其情況,虛線為差0.5 碼片的相關(guān)輸出,輸出功率降低6dB,所以會(huì)造成估計(jì)精度降低.表2 為仿真進(jìn)行50 次的頻率均方誤差(RMS).

    圖5 本地碼和接收碼未對(duì)齊對(duì)相關(guān)輸出峰值功率的影響

    從表1, 表2 可以看出Pisarenko 的多普勒估計(jì)精度比ARMA 及FFT 高,雖然其估計(jì)誤差受信噪比的影響,但是在考慮低信噪比GPS 信號(hào)情況下,其估計(jì)誤差完全能夠滿(mǎn)足實(shí)際情況.

    表2 1ms 數(shù)據(jù)考慮碼片未對(duì)齊情況頻率估計(jì)RMS(Hz)

    3.2 運(yùn)算量的研究及對(duì)平均多普勒捕獲時(shí)間的對(duì)比

    仿真條件: 估計(jì)的多普勒范圍為±10kHz, 多普勒步長(zhǎng)500Hz,多普勒粗捕誤差250 為Hz,則頻率單元為21 個(gè),積分時(shí)間為1ms,在MATLAB 仿真環(huán)境下,進(jìn)行100 次仿真,FLL 平均需50ms 達(dá)到穩(wěn)定(頻率誤差達(dá)到20Hz 內(nèi)),每次需要40 e6次浮點(diǎn)運(yùn)算, 而本文采用的頻譜估計(jì)方法要20 e6 次浮點(diǎn)運(yùn)算,是FLL 收斂時(shí)間的一半.

    由式(8),若Gc=21,PFA=0.01,PD=0.95,Tc=0.001 情況下,若Ts=0.05 得平均捕獲時(shí)間0.0616s, 若Ts=0.025 得平均捕獲時(shí)間0.0363s,可以減少42%的時(shí)間,可見(jiàn),采用高分辨率的載波多普勒估計(jì)替代FLL 對(duì)多普勒剩余殘差估計(jì)可以大大減少?gòu)牟东@到跟蹤的時(shí)間.

    4 結(jié)論

    以上研究表明,對(duì)于解擴(kuò)后的GPS 信號(hào),視為平穩(wěn)隨機(jī)過(guò)程和高斯白噪聲的和, 通過(guò)諧波恢復(fù)Pisarenko 譜估計(jì)法,對(duì)GPS 多普勒進(jìn)行高分辨率估計(jì)是完全可行的, 可以取得非常好的的精度,以及穩(wěn)定的結(jié)果,而且計(jì)算量也相對(duì)FLL 對(duì)多普勒剩余殘差估計(jì)小得多,可以減少40%左右的GPS 接收機(jī)從捕獲到跟蹤的時(shí)間。S

    [1]黃振,陸建華,楊士中.衛(wèi)星通信中多普勒頻移的快速捕獲[J].電子學(xué)報(bào),2003,31(7)1052-1056.

    [2]Kaplan,E.,Understanding GPS:Principles and Applications,Artech House,1996.

    [3]Tsui,J.,Fundamentals of Global Positioning System Receivers: A Software Approach, Wiley Inter-Science,2000.

    [4]Lin,D.,Tsui,J.,”Comparison of Acquisition Methods for Software GPS Receiver”. Proceedings of ION, pp.2385-2391, September 2000.

    [5]Parkinson,B., et al.(editors),Global Positioning System. Volume 1.American Institute of Astronautics and Aeronautics, 1996.

    [6]張賢達(dá).現(xiàn)代信號(hào)處理[M].2 版.北京:清華大學(xué)出版社,2003.

    [7]Zhang X D .Estimation of frequencies of sinusoids in colored ARMA noise via singular value decomposition.Proc.IEEE ISCAS’89,Portland,OH.1989:1315-1318.

    猜你喜歡
    譜估計(jì)虛警環(huán)路
    頻率步進(jìn)連續(xù)波雷達(dá)電磁輻射2階互調(diào)虛警干擾效應(yīng)規(guī)律
    一種電阻式應(yīng)變傳感器的數(shù)據(jù)處理方法
    上海市中環(huán)路標(biāo)線調(diào)整研究
    上海公路(2018年4期)2018-03-21 05:57:46
    基于MATLAB實(shí)現(xiàn)的AR模型功率譜估計(jì)
    空管自動(dòng)化系統(tǒng)二次代碼劫機(jī)虛警分析
    BIT虛警影響因素分析和解決方案
    經(jīng)典功率譜估計(jì)方法的研究
    Welch譜估計(jì)的隨機(jī)誤差與置信度
    Buck-Boost變換器的環(huán)路補(bǔ)償及仿真
    單脈沖雷達(dá)導(dǎo)引頭角度跟蹤環(huán)路半實(shí)物仿真
    午夜老司机福利剧场| 中文字幕亚洲精品专区| 黄色怎么调成土黄色| 国产精品人妻久久久影院| 三级国产精品片| 国产成人免费无遮挡视频| 久久人妻熟女aⅴ| 日韩一区二区视频免费看| 欧美日韩亚洲高清精品| 精品一区二区三区视频在线| 亚洲av中文av极速乱| 日本欧美国产在线视频| 午夜激情久久久久久久| 久久久久久久久大av| 亚洲真实伦在线观看| 亚洲av成人精品一区久久| 啦啦啦中文免费视频观看日本| 亚洲精品色激情综合| 午夜福利在线观看免费完整高清在| 在线观看www视频免费| 久久人人爽av亚洲精品天堂| 深夜a级毛片| 色94色欧美一区二区| a 毛片基地| 国产日韩一区二区三区精品不卡 | 啦啦啦视频在线资源免费观看| 日本午夜av视频| 性高湖久久久久久久久免费观看| 久久精品国产亚洲av天美| 老司机亚洲免费影院| 国产亚洲最大av| 一区在线观看完整版| 岛国毛片在线播放| 91成人精品电影| 国产黄色视频一区二区在线观看| 成人18禁高潮啪啪吃奶动态图 | 国产成人精品一,二区| 久久久久精品性色| 亚洲真实伦在线观看| 国产淫语在线视频| 一区二区三区免费毛片| 日韩不卡一区二区三区视频在线| 欧美+日韩+精品| 欧美一级a爱片免费观看看| 中文字幕制服av| 99九九线精品视频在线观看视频| 亚洲熟女精品中文字幕| 我的老师免费观看完整版| 伊人久久精品亚洲午夜| 成年人免费黄色播放视频 | 久久精品国产鲁丝片午夜精品| 97精品久久久久久久久久精品| 热99国产精品久久久久久7| 春色校园在线视频观看| 一本色道久久久久久精品综合| 高清在线视频一区二区三区| 综合色丁香网| 一级毛片电影观看| 色视频在线一区二区三区| 色视频在线一区二区三区| 亚洲一级一片aⅴ在线观看| av有码第一页| 午夜视频国产福利| 一区二区三区四区激情视频| 亚洲av福利一区| 夫妻午夜视频| 午夜福利视频精品| 大香蕉97超碰在线| 国产精品久久久久成人av| 免费观看在线日韩| 涩涩av久久男人的天堂| 特大巨黑吊av在线直播| 丰满饥渴人妻一区二区三| 免费观看无遮挡的男女| 精品国产一区二区久久| 2018国产大陆天天弄谢| 午夜福利在线观看免费完整高清在| 精品人妻熟女毛片av久久网站| 男女无遮挡免费网站观看| 在线观看人妻少妇| av在线app专区| 日日啪夜夜爽| 亚洲第一av免费看| 一级爰片在线观看| 美女视频免费永久观看网站| 日本黄色片子视频| 久久人人爽av亚洲精品天堂| 国产精品国产三级国产av玫瑰| 婷婷色麻豆天堂久久| 国产高清不卡午夜福利| 亚洲人成网站在线观看播放| 高清在线视频一区二区三区| 久久精品国产亚洲网站| 亚洲国产最新在线播放| 我要看日韩黄色一级片| 最近最新中文字幕免费大全7| 久久99精品国语久久久| 你懂的网址亚洲精品在线观看| 国产毛片在线视频| 男人添女人高潮全过程视频| 91精品国产九色| 成年人午夜在线观看视频| 丰满少妇做爰视频| 婷婷色综合www| 精品久久国产蜜桃| 少妇被粗大猛烈的视频| 一级,二级,三级黄色视频| 欧美精品一区二区大全| 国产亚洲5aaaaa淫片| 国产极品粉嫩免费观看在线 | 精品国产国语对白av| 最近手机中文字幕大全| 曰老女人黄片| 国产乱来视频区| 久久久久久久大尺度免费视频| 卡戴珊不雅视频在线播放| 我的老师免费观看完整版| 日韩成人伦理影院| 99re6热这里在线精品视频| 国产精品99久久久久久久久| 伊人久久精品亚洲午夜| 久久国产精品男人的天堂亚洲 | 国产精品.久久久| 精品一区二区三卡| 日韩一区二区视频免费看| 久久毛片免费看一区二区三区| 色94色欧美一区二区| 永久网站在线| 欧美日本中文国产一区发布| 丰满少妇做爰视频| 亚洲精品国产色婷婷电影| 在线观看美女被高潮喷水网站| 日日摸夜夜添夜夜爱| 日本猛色少妇xxxxx猛交久久| av免费观看日本| 黄片无遮挡物在线观看| 亚洲精品久久久久久婷婷小说| 国产女主播在线喷水免费视频网站| 亚洲第一区二区三区不卡| av.在线天堂| 高清午夜精品一区二区三区| 内射极品少妇av片p| 国产淫语在线视频| 免费少妇av软件| 香蕉精品网在线| 国产 精品1| 女性被躁到高潮视频| 人妻人人澡人人爽人人| 91久久精品国产一区二区成人| 亚洲情色 制服丝袜| 激情五月婷婷亚洲| 有码 亚洲区| 日韩三级伦理在线观看| 亚洲伊人久久精品综合| 欧美日韩精品成人综合77777| 一区二区三区精品91| a级毛片在线看网站| 国产免费一级a男人的天堂| 日本免费在线观看一区| 亚洲成人av在线免费| 成人亚洲精品一区在线观看| 岛国毛片在线播放| 国产成人精品福利久久| 成人综合一区亚洲| 嘟嘟电影网在线观看| 欧美成人午夜免费资源| 少妇丰满av| 婷婷色麻豆天堂久久| 久久久亚洲精品成人影院| 亚洲第一区二区三区不卡| 中文天堂在线官网| 午夜日本视频在线| 日韩伦理黄色片| 一级a做视频免费观看| 中文字幕av电影在线播放| videossex国产| 午夜91福利影院| 亚洲av免费高清在线观看| 欧美xxxx性猛交bbbb| 久久精品国产亚洲av天美| 国产又色又爽无遮挡免| 另类亚洲欧美激情| 亚洲美女搞黄在线观看| 国产91av在线免费观看| 免费黄网站久久成人精品| 亚洲国产毛片av蜜桃av| 国产成人午夜福利电影在线观看| h日本视频在线播放| 欧美xxⅹ黑人| a级毛片在线看网站| 欧美bdsm另类| 久久ye,这里只有精品| 99九九线精品视频在线观看视频| 久久韩国三级中文字幕| 亚洲人与动物交配视频| 一级av片app| 国产精品伦人一区二区| 精品亚洲成a人片在线观看| 色5月婷婷丁香| 99精国产麻豆久久婷婷| 成人毛片a级毛片在线播放| 成人亚洲欧美一区二区av| 波野结衣二区三区在线| 精品人妻偷拍中文字幕| 免费播放大片免费观看视频在线观看| 亚洲精品,欧美精品| 久久热精品热| 亚洲精品,欧美精品| 久久久久精品久久久久真实原创| 搡女人真爽免费视频火全软件| 七月丁香在线播放| 欧美另类一区| 丰满迷人的少妇在线观看| 91在线精品国自产拍蜜月| 国精品久久久久久国模美| a 毛片基地| 秋霞在线观看毛片| 国产精品久久久久成人av| 国产精品久久久久久精品电影小说| 男人爽女人下面视频在线观看| 水蜜桃什么品种好| 国产69精品久久久久777片| 中文字幕制服av| 精品亚洲成国产av| 国产欧美日韩精品一区二区| 亚洲精华国产精华液的使用体验| 十八禁网站网址无遮挡 | 一边亲一边摸免费视频| 少妇的逼好多水| 噜噜噜噜噜久久久久久91| 日本爱情动作片www.在线观看| 六月丁香七月| 男人和女人高潮做爰伦理| 久久精品夜色国产| 免费看日本二区| 自线自在国产av| 欧美一级a爱片免费观看看| 亚洲,欧美,日韩| 亚洲经典国产精华液单| 夜夜看夜夜爽夜夜摸| av免费观看日本| av天堂久久9| 毛片一级片免费看久久久久| 美女内射精品一级片tv| 好男人视频免费观看在线| 国产精品久久久久成人av| 老熟女久久久| 亚洲在久久综合| 2022亚洲国产成人精品| 日韩制服骚丝袜av| 黑人猛操日本美女一级片| 能在线免费看毛片的网站| 在线观看www视频免费| 秋霞伦理黄片| 99精国产麻豆久久婷婷| 亚洲成人av在线免费| 国产精品99久久久久久久久| 国产午夜精品一二区理论片| 国产伦理片在线播放av一区| 在线看a的网站| 中国三级夫妇交换| 亚洲精品日本国产第一区| 中文字幕久久专区| 欧美激情国产日韩精品一区| 久久久久久久久久久免费av| 国产精品成人在线| 18禁在线播放成人免费| 美女内射精品一级片tv| 男人爽女人下面视频在线观看| 午夜精品国产一区二区电影| 伊人久久精品亚洲午夜| 国产精品成人在线| 美女xxoo啪啪120秒动态图| 久久精品夜色国产| 欧美+日韩+精品| 特大巨黑吊av在线直播| 99热这里只有是精品50| 国产伦精品一区二区三区四那| 国产成人一区二区在线| 男的添女的下面高潮视频| www.色视频.com| 一区二区三区四区激情视频| 亚洲国产色片| 国产视频内射| 久久精品国产a三级三级三级| 国产又色又爽无遮挡免| 国产熟女欧美一区二区| 国产真实伦视频高清在线观看| 中文精品一卡2卡3卡4更新| 在线观看免费高清a一片| 国产成人aa在线观看| 精品人妻一区二区三区麻豆| 日韩欧美一区视频在线观看 | 免费不卡的大黄色大毛片视频在线观看| 午夜激情久久久久久久| 国产av国产精品国产| 熟女人妻精品中文字幕| 蜜桃久久精品国产亚洲av| av视频免费观看在线观看| 久久久国产精品麻豆| 亚洲av中文av极速乱| 亚洲国产精品专区欧美| 青春草亚洲视频在线观看| 亚洲综合精品二区| 成人美女网站在线观看视频| 国产成人一区二区在线| 全区人妻精品视频| 亚洲欧美精品专区久久| 欧美精品国产亚洲| 91精品国产国语对白视频| 亚洲欧美日韩另类电影网站| 久久久a久久爽久久v久久| 日本黄色日本黄色录像| a级毛片在线看网站| 日日撸夜夜添| 久久99热这里只频精品6学生| 精品人妻熟女av久视频| av黄色大香蕉| 国产成人一区二区在线| 两个人的视频大全免费| 九九久久精品国产亚洲av麻豆| 最近中文字幕高清免费大全6| 国产精品99久久99久久久不卡 | 天天操日日干夜夜撸| 大码成人一级视频| 一区二区三区四区激情视频| 成年人午夜在线观看视频| 亚洲不卡免费看| 国产毛片在线视频| 亚洲欧美精品专区久久| 美女xxoo啪啪120秒动态图| 国产无遮挡羞羞视频在线观看| 亚洲三级黄色毛片| 全区人妻精品视频| av福利片在线| 啦啦啦中文免费视频观看日本| 美女cb高潮喷水在线观看| 久久久a久久爽久久v久久| 久久久精品免费免费高清| 久久亚洲国产成人精品v| 成年人免费黄色播放视频 | 国产成人精品无人区| 国产一区二区在线观看av| 国精品久久久久久国模美| 高清毛片免费看| 久久鲁丝午夜福利片| 男女边摸边吃奶| 亚洲自偷自拍三级| 色婷婷久久久亚洲欧美| 内地一区二区视频在线| 伦精品一区二区三区| 2022亚洲国产成人精品| 成人影院久久| a 毛片基地| 中文字幕久久专区| 天堂中文最新版在线下载| av在线观看视频网站免费| 国产精品无大码| 国产精品一区二区在线不卡| 老司机影院成人| 中文字幕人妻熟人妻熟丝袜美| 国产成人午夜福利电影在线观看| 两个人的视频大全免费| 精品午夜福利在线看| 国产精品熟女久久久久浪| 国产日韩欧美亚洲二区| 色网站视频免费| av免费在线看不卡| 精品久久久精品久久久| 亚洲国产精品国产精品| 欧美最新免费一区二区三区| 日本wwww免费看| 18禁在线播放成人免费| 国产精品久久久久久久久免| 一个人看视频在线观看www免费| 只有这里有精品99| 国产精品欧美亚洲77777| 中文字幕精品免费在线观看视频 | 大话2 男鬼变身卡| 免费不卡的大黄色大毛片视频在线观看| videossex国产| 99久久人妻综合| 成人亚洲精品一区在线观看| 精品国产乱码久久久久久小说| 日韩成人av中文字幕在线观看| a级毛色黄片| 亚洲欧美成人精品一区二区| av网站免费在线观看视频| 日韩成人av中文字幕在线观看| 六月丁香七月| 国产深夜福利视频在线观看| 五月开心婷婷网| 搡女人真爽免费视频火全软件| 国产成人免费无遮挡视频| 狠狠精品人妻久久久久久综合| 国产一区二区三区综合在线观看 | 国产亚洲91精品色在线| 国产精品熟女久久久久浪| 赤兔流量卡办理| 一区二区三区精品91| 少妇的逼好多水| 97精品久久久久久久久久精品| 欧美97在线视频| 寂寞人妻少妇视频99o| 精品久久国产蜜桃| 99热全是精品| 亚洲国产日韩一区二区| 色吧在线观看| 国产精品嫩草影院av在线观看| 久久久a久久爽久久v久久| 国产精品久久久久久久电影| 丰满人妻一区二区三区视频av| 狂野欧美白嫩少妇大欣赏| 在线精品无人区一区二区三| 一区在线观看完整版| 日日爽夜夜爽网站| 久久精品国产亚洲av天美| 成人亚洲精品一区在线观看| 亚洲av.av天堂| 人妻人人澡人人爽人人| 亚洲丝袜综合中文字幕| 男人舔奶头视频| 一二三四中文在线观看免费高清| 久久久久久久亚洲中文字幕| 国产永久视频网站| 黄色怎么调成土黄色| 最近中文字幕高清免费大全6| 久久国产亚洲av麻豆专区| 高清毛片免费看| 国产免费一级a男人的天堂| 黄色毛片三级朝国网站 | 91aial.com中文字幕在线观看| 国产精品嫩草影院av在线观看| 成人亚洲精品一区在线观看| 男女边吃奶边做爰视频| 又黄又爽又刺激的免费视频.| av福利片在线| 亚洲精品第二区| 国产白丝娇喘喷水9色精品| 国产精品免费大片| 亚洲av在线观看美女高潮| 男女啪啪激烈高潮av片| 国产黄片视频在线免费观看| 午夜视频国产福利| 亚洲精品乱久久久久久| 人妻系列 视频| 精品国产露脸久久av麻豆| 嫩草影院新地址| 亚洲国产精品国产精品| 国产一区二区三区av在线| 乱系列少妇在线播放| 亚洲伊人久久精品综合| 欧美日韩综合久久久久久| 成人国产av品久久久| 日韩亚洲欧美综合| 麻豆成人av视频| 大片电影免费在线观看免费| 日本黄色片子视频| 美女脱内裤让男人舔精品视频| 免费不卡的大黄色大毛片视频在线观看| 亚洲婷婷狠狠爱综合网| 亚洲精华国产精华液的使用体验| 草草在线视频免费看| 一级黄片播放器| 五月开心婷婷网| 精品久久久久久久久av| 秋霞在线观看毛片| 高清午夜精品一区二区三区| 国产亚洲午夜精品一区二区久久| 一级黄片播放器| 成人亚洲精品一区在线观看| 91久久精品国产一区二区三区| 国产深夜福利视频在线观看| 汤姆久久久久久久影院中文字幕| 亚洲精品国产av成人精品| 欧美老熟妇乱子伦牲交| 七月丁香在线播放| 最近的中文字幕免费完整| 亚洲,欧美,日韩| 亚洲精品色激情综合| 亚洲精品乱码久久久v下载方式| 纵有疾风起免费观看全集完整版| 久久青草综合色| 免费人成在线观看视频色| 一区二区av电影网| 天堂俺去俺来也www色官网| av天堂中文字幕网| 国产精品免费大片| 精品99又大又爽又粗少妇毛片| 国产女主播在线喷水免费视频网站| 黄色怎么调成土黄色| 成人综合一区亚洲| 国产黄频视频在线观看| 国语对白做爰xxxⅹ性视频网站| 国产精品不卡视频一区二区| av免费在线看不卡| 亚洲精品中文字幕在线视频 | 国产真实伦视频高清在线观看| 麻豆成人午夜福利视频| 国产精品.久久久| 丰满乱子伦码专区| 欧美精品高潮呻吟av久久| 国产午夜精品久久久久久一区二区三区| 亚洲va在线va天堂va国产| 亚洲欧美成人精品一区二区| 精品一区二区免费观看| 少妇人妻精品综合一区二区| 久久久久久久久久人人人人人人| 亚洲精品国产色婷婷电影| 丝瓜视频免费看黄片| 欧美精品一区二区大全| 91在线精品国自产拍蜜月| 久久精品熟女亚洲av麻豆精品| 午夜激情福利司机影院| 色94色欧美一区二区| 亚洲欧美日韩卡通动漫| 不卡视频在线观看欧美| 亚洲av电影在线观看一区二区三区| 中文资源天堂在线| 日韩成人av中文字幕在线观看| 人人妻人人爽人人添夜夜欢视频 | 麻豆精品久久久久久蜜桃| 国产极品天堂在线| 秋霞在线观看毛片| 日产精品乱码卡一卡2卡三| 久久人人爽人人片av| 免费高清在线观看视频在线观看| 亚洲欧美一区二区三区国产| 男的添女的下面高潮视频| 国产淫片久久久久久久久| 高清黄色对白视频在线免费看 | 亚洲欧美精品自产自拍| 亚洲四区av| 国产熟女欧美一区二区| 人妻人人澡人人爽人人| 插逼视频在线观看| 我要看黄色一级片免费的| 九九在线视频观看精品| 高清欧美精品videossex| 欧美xxⅹ黑人| 亚洲国产色片| 欧美日韩综合久久久久久| av在线播放精品| 久久久欧美国产精品| 亚洲av男天堂| 久久精品久久久久久久性| 欧美日韩在线观看h| 18禁动态无遮挡网站| 26uuu在线亚洲综合色| 久久精品国产亚洲av天美| 久久久久久久久大av| 老女人水多毛片| 自拍欧美九色日韩亚洲蝌蚪91 | xxx大片免费视频| 久久精品熟女亚洲av麻豆精品| 成人特级av手机在线观看| 色5月婷婷丁香| 国产高清有码在线观看视频| h日本视频在线播放| 插阴视频在线观看视频| 亚洲国产精品一区二区三区在线| 黄色视频在线播放观看不卡| 在线观看免费日韩欧美大片 | 国产片特级美女逼逼视频| 久久久亚洲精品成人影院| 丰满饥渴人妻一区二区三| 婷婷色综合大香蕉| 精品亚洲乱码少妇综合久久| 十分钟在线观看高清视频www | 欧美bdsm另类| 久热这里只有精品99| 亚洲高清免费不卡视频| 国产极品天堂在线| 又黄又爽又刺激的免费视频.| 午夜免费鲁丝| 亚洲av欧美aⅴ国产| 99久久人妻综合| 久久久久人妻精品一区果冻| 人体艺术视频欧美日本| 啦啦啦啦在线视频资源| 国产成人精品婷婷| 性高湖久久久久久久久免费观看| 亚洲人与动物交配视频| 九九久久精品国产亚洲av麻豆| 在线播放无遮挡| 性色avwww在线观看| 欧美日韩视频精品一区| 日韩精品免费视频一区二区三区 | 国产欧美日韩综合在线一区二区 | 色吧在线观看| 99国产精品免费福利视频| 五月天丁香电影| 成人国产麻豆网| 日韩在线高清观看一区二区三区| 最近手机中文字幕大全| 美女主播在线视频| 少妇的逼水好多| 午夜福利影视在线免费观看| 国产一区二区三区av在线| 一级毛片aaaaaa免费看小| 六月丁香七月| 黄片无遮挡物在线观看| 黄色欧美视频在线观看| 国产成人a∨麻豆精品| 亚洲精品,欧美精品| 国产色婷婷99| 精品卡一卡二卡四卡免费| 亚洲欧美精品专区久久| 国产淫语在线视频| av国产久精品久网站免费入址| 亚洲真实伦在线观看| 亚洲三级黄色毛片| 91精品国产国语对白视频| 欧美亚洲 丝袜 人妻 在线| 成人美女网站在线观看视频|