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

    關(guān)于射電望遠鏡臺址航空信號處理方法的研究?

    2018-04-02 02:46:50許紅瑞陳卯蒸軍劉璇袁
    天文學報 2018年2期
    關(guān)鍵詞:射電航跡高峰

    許紅瑞陳卯蒸 劉 奇 尹 航 王 軍劉 璇袁 力

    (1中國科學院新疆天文臺烏魯木齊830011)

    (2中國科學院射電天文重點實驗室南京210008)

    (3中國科學院大學北京100049)

    1 引言

    H I(0.235–0.386 z)和1665 OH(0.448–0.624 z)是研究天體移動及規(guī)律預測的重要譜線,但由于紅移,譜線頻率與航空信號頻率發(fā)生重疊,如表1所示.民用航空信號(civil aviation signal,以下簡稱CAS)的工作波段為1–1.5 GHz之間約130 MHz帶寬的L波段.在正常的射電天文觀測中,CAS信號會從天線旁瓣進入接收機系統(tǒng),降低系統(tǒng)信噪比和觀測數(shù)據(jù)的有效性,增加數(shù)據(jù)處理的復雜性,嚴重影響了射電天文在此頻段內(nèi)的觀測效果,降低大型射電科學儀器的使用價值[1?2].鑒于民航航線的不斷增多和原有航線的調(diào)整,在新一代射電望遠鏡的選址工作中,雖然考慮了對地面無線電信號的電磁兼容和電磁干擾,但是對來自CAS的干擾卻難以避免.

    CAS主要來自民用航空飛機的機載廣播式自動相關(guān)監(jiān)視系統(tǒng)(Automatic Dependent Surveillance-Broadcast,以下簡稱ADS-B)和機載測距機系統(tǒng)(Distance Measuring Equipment,以下簡稱DME),ADS-B和DME信號占用L波段帶寬如表1所示.ADS-B的作用距離可達300 km以上,系統(tǒng)發(fā)射功率為51–57 dBm,信號編碼方式為脈沖編碼,比特率為1 Mb/s[2].DME系統(tǒng)機載詢問器發(fā)射頻率在1025–1150 MHz,帶寬占用125 MHz,信號帶寬1–3 MHz,響應(yīng)時間128μs[3].即ADS-B和機載DME產(chǎn)生的電磁干擾信號具有瞬態(tài)、極化和高輻射強度等特征.目前射電觀測數(shù)據(jù)的消干擾方法為屏蔽子積分頻譜中有干擾的通道,基于現(xiàn)有的消干擾方法,在干擾信號越多的情況下,有用的頻譜通道越少,獲得的天文數(shù)據(jù)的有用信息越少[4].國內(nèi)外針對由瞬時脈沖信號造成的射頻干擾多集中在后期處理,Willem Baan團隊提出避免干擾信號比抑制干擾信號更有意義,通過利用自然地形,提前獲知所觀測區(qū)域有哪些干擾頻譜存在,盡量避免或者消除干擾,從而有效建立觀測寧靜區(qū)[5].

    為避免射電望遠鏡觀測H I(0.235–0.386 z)和1665 OH(0.448–0.624 z)[6?8]紅移過程中出現(xiàn)的CAS,本文利用已開發(fā)的軟件無線電平臺,通過對ADS-B報文解析獲得航跡樣本點.基于航跡樣本點分布的時間特點,采用最小二乘多項式擬合,獲得航跡樣本點在時間上的分布高峰時段,結(jié)合獲得的高峰時段對航跡樣本點進行基于球面距離的二分k-均值聚類算法分析,從而獲得射電天文臺址周邊高峰時段飛機分布的主要區(qū)域.為安排射電天文觀測、避開CAS、建設(shè)射電天文觀測寧靜區(qū)及消干擾策略提供依據(jù).

    表1 飛機導航信號及紅移譜線在1025–1150 MHz頻段占用帶寬Table 1 Frequency occupation of the aircraft navigation signal and the redshift spectrum in 1025–1150 MHz

    2 數(shù)據(jù)存儲格式及解析過程

    本文基于已開發(fā)搭建的軟件無線電平臺,以新疆天文臺本部大樓為測試點,經(jīng)過測試解析,共獲得航跡信息12萬條,數(shù)據(jù)存儲格式如表2所示,其中第1列是飛機唯一的ICAO(國際民用航空組織)地址,Seen表示的是接收到飛機信號的實時時間,Altitude、Latitude和Longitude分別代表飛機此時的海拔高度、緯度和經(jīng)度.ADS-B信號經(jīng)天線接收,經(jīng)過射頻處理之后成為基帶信號,經(jīng)傳輸?shù)竭_PC端進行報文解析,獲得航跡信息包括接收時間、地址碼、航班號、經(jīng)緯度、高度等多個報文項.解析流程如圖1所示,其中,RF為射頻前端,AD為數(shù)模轉(zhuǎn)換子板,FPGA為現(xiàn)場可編程陣列,USRP是通用軟件無線電外設(shè).

    表2 ADS-B報文解析存儲格式Table 2 ADS-B message parse storage format

    圖1 ADS-B報文接收解析框圖Fig.1 Receiving and parse block diagram of ADS-B message

    3 飛機到測站距離及到達測站的干擾強度估算

    3.1 計算空間距離

    由于飛機在飛行過程中是動態(tài)的,接收到的點均為離散點,為了更直觀地表現(xiàn)飛機到測試點距離在時間上的分布,計算C、D兩點距離,其中C為飛機所處位置,其緯度、經(jīng)度、海拔分別為a1、c1、b1,D為測試點所處位置,其緯度、經(jīng)度、海拔分別為a2、c2、b2.利用(1)–(3)式,分別將緯度、經(jīng)度、海拔轉(zhuǎn)換為3維空間坐標x1、y1、z1和x2、y2、z2,利用(4)式計算C、D兩點距離d,根據(jù)距離畫出飛機在時間上的分布.

    3.2 估算干擾強度

    本文采用自由空間損耗模型估算信號強度,空間損耗公式為:

    其中f為頻率,單位為MHz,d為距離,單位為km.

    3.3 飛機到測站距離隨著時間分布及到測站信號強度估算結(jié)果

    根據(jù)(4)式計算了一天中飛機到測站距離隨著時間的分布圖,如圖2所示.由圖2可知:飛機到測站距離分布集中在20 km左右,分布范圍在10–80 km.相應(yīng)地,若飛機在視線范圍,到達測站的自由空間損耗在113–131 dB之間,若因地面等遮擋,飛機不在視線范圍內(nèi),損耗一般會比同距離自由空間更大些,由此可以估算航空信號在觀測站的干擾強度.

    圖2 飛機到測站距離隨時間分布圖Fig.2 The distribution of the distance between plane and station with time

    4 基于最小二乘法和二分k-均值聚類算法的航跡分析方法

    航線的主要分布時段和區(qū)域是眾多航線在時間空間上反映出來的主要特征,提取航線高峰時段分布區(qū)域的技術(shù)路線如圖3所示.獲取測試點航線高峰時段主要區(qū)域基本思路為:首先根據(jù)獲得的航跡點數(shù)量和時間的關(guān)系,分析每天接收的航跡點是否穩(wěn)定,針對航跡點統(tǒng)計量在24 h的分布情況,采用最小二乘法的多項式擬合獲得每天的高峰時段.在獲得高峰時段后,航跡點按照高峰時段篩選出經(jīng)緯度信息,之后應(yīng)用基于球面距離的二分k-均值聚類算法對高峰時段的經(jīng)緯度信息進行聚類,獲得主要航線分布范圍.

    4.1 高峰時段分析

    4.1.1 擬合原理

    在2017年3月5日到2017年3月17日之間,每天接收解析信息總次數(shù)如表3所示.在不考慮平臺性能和天氣原因的情況下,總的來看,每天接收信號數(shù)量大致是穩(wěn)定的.在此基礎(chǔ)上,分析一天的高峰時段.本文認為每天的高峰時段應(yīng)該是相同的,通過統(tǒng)計3月5日到3月17日之間每天24 h中每小時接收信號的總頻次,得到13組數(shù)據(jù)點,每組24個統(tǒng)計點.首先計算對應(yīng)各統(tǒng)計點每天平均樣本數(shù),得到一組24個平均數(shù)據(jù),再按這組平均數(shù)據(jù)進行基于最小二乘法的多項式擬合,得到航跡樣本點隨時間的分布趨勢.

    圖3 技術(shù)路線圖Fig.3 Technology roadmap

    表3 13組樣本的統(tǒng)計量Table 3 Statistics of 13 samples

    4.1.2 最小二乘多項式擬合

    給定離散數(shù)據(jù)點(xi,yi)(i=0,1,2,···,m),xi代表時間點,yi表示對應(yīng)時間點接收的統(tǒng)計量,Φ為次數(shù)不超過n(n≤m)的多項式構(gòu)成的函數(shù)類,通過給定的離散數(shù)據(jù)點建立擬合函數(shù)為:使得(6)式的平方和I最小.

    滿足(6)式的pn(x)為最小二乘擬合多項式.

    (8)式關(guān)于a0,a1,···,an的線性方程組,用矩陣表示為:

    方程組(9)的矩陣是一個對稱正定矩陣,故存在唯一解.從(9)式中解出ak(k=0,1,···,n),從而可得多項式:

    4.1.3 相關(guān)系數(shù)及其顯著性檢驗

    將平均樣本點(xi,yi)做多項式擬合時,還并不太了解x與y之間關(guān)系的密切程度.為此要用相關(guān)系數(shù)r進行判斷,其定義為:

    4.2 基于球面距離的二分k-均值航跡聚類算法

    本文目的是將目標航跡進行聚類[9]分析,將相似航跡點聚在一起,在眾多雜亂無章的航跡中找到出現(xiàn)航跡概率較大的區(qū)域范圍.雖然k均值方法也可以很好地將數(shù)據(jù)進行聚類分析,但是二分k-均值聚類算法是對k均值聚類算法的改進,一方面可以克服k均值聚類算法收斂于局部最小值的問題,另一方面可以加快運算速度.因本文數(shù)據(jù)為經(jīng)緯度信息,為了保證方法的普遍適用性,本文采用基于球面距離的二分k-均值聚類算法對航跡樣本點進行處理.

    4.2.1 計算球面距離

    由于我們的數(shù)據(jù)主要是經(jīng)緯度,但這些信息對于距離計算并不足夠.在北極附近每走幾米的經(jīng)度變化可能達到10?;而在赤道附近走相同的距離,帶來的變化可能僅零點幾米.為了統(tǒng)一規(guī)范化距離的計算,此處采用球面余弦定理來計算兩點之間的距離.A點緯度、經(jīng)度為(β1,α1),B點緯度、經(jīng)度為(β2,α2),R為地球半徑,則A、B兩點間的球面距離d1為:

    4.2.2 航跡聚類處理算法

    本文提出了基于聚類的航跡分析算法,其輸入包括,輸入1:航跡集合Q={q1,q2,q3,···,qn},其中qi(i=1,2,···,n)代表航跡點;輸入2: 指定k為希望得到的簇數(shù)目.輸出為航跡聚類結(jié)果,算法具體執(zhí)行步驟描述如下:利用二分k-均值算法對航跡進行分析,在算法中使用的距離為(12)式中的球面距離,得到航跡聚類結(jié)果,其中二分k-均值算法偽代碼描述為:

    (1)將所有點看成一個簇;

    (2)當簇數(shù)目小于k時,對于每一個簇,進行:

    (i)計算總誤差;

    (ii)在給定的簇上面進行k均值聚類;

    (iii)計算將該簇一分為二之后的總誤差;

    (3)選擇誤差最小的那個簇進行劃分操作.

    5 實驗結(jié)果分析與方法驗證

    5.1 基于最小二乘法的曲線擬合結(jié)果及分析

    基于2017年3月5日到2017年3月17日之間一天24 h的13組統(tǒng)計數(shù)據(jù)(每組數(shù)據(jù)有24個統(tǒng)計點),按一天24個統(tǒng)計點,先計算對應(yīng)各統(tǒng)計點每天平均樣本數(shù),再按這組平均數(shù)據(jù)進行基于最小二乘法的多項式擬合獲得的趨勢走向如圖4所示.由圖中連續(xù)的擬合曲線可以看出,一天之中的高峰時段分布在UTC時間的3點、9點和15點內(nèi),圖中離散點為2017年3月5日到2017年3月17日之間,24個統(tǒng)計點所對應(yīng)的實測統(tǒng)計值.

    經(jīng)計算對應(yīng)時間點的平均樣本數(shù)的相關(guān)性,以r=0.5為閾值,其中相關(guān)系數(shù)在閾值以上的占92%,因此我們認為得到的高峰時段的數(shù)據(jù)是可靠的,基于此,進行航跡聚類的實驗.

    圖4 2017年3月5日—2017年3月17日之間每日平均樣本數(shù)擬合分布趨勢Fig.4 Fitting of the daily average sample counting distribution between 2017-03-05—2017-03-17

    5.2 高峰時段航跡分布范圍聚類劃分

    采用實測數(shù)據(jù)對航跡聚類算法進行測試,數(shù)據(jù)分為6組,前3組作為測試劃分范圍,后3組驗證方法的有效性.第1組為3月5日–3月11日UTC時間3點內(nèi)所接收的全部航跡點,第2組為3月5日–3月11日UTC時間9點內(nèi)所接收的全部航跡點,第3組為3月5日–3月11日UTC時間15點內(nèi)所接收的全部航跡點,第4組為3月12日–3月17日UTC時間3點內(nèi)所接收的全部航跡點,第5組為3月12日–3月17日UTC時間9點內(nèi)所接收的全部航跡點,第6組為3月12日–3月17日UTC時間15點內(nèi)所接收的全部航跡點,由1、2、3組航跡數(shù)據(jù)聚類結(jié)果分別如圖5(左側(cè))、6(左側(cè))、7(左側(cè))所示.根據(jù)1、2、3組數(shù)據(jù)航跡聚類結(jié)果我們可以劃定高峰時段航跡出現(xiàn)的主要范圍區(qū)域.圖中紅色圓點代表新疆天文臺總部大樓,藍色圓點代表航跡點,紫色、紅色、綠色圈定的區(qū)域為在此時段內(nèi)航跡點主要聚集范圍.

    圖5 UTC 3點航跡聚類(左)和驗證航跡聚類(右)Fig.5 Track clustering during UTC 3(left)and track clustering verifying during UTC 3(right)

    5.3 方法驗證

    根據(jù)實驗劃定高峰時段范圍區(qū)域,將后3組高峰時段航跡點作為驗證數(shù)據(jù),聚類結(jié)果如圖5(右側(cè))、6(右側(cè))、7(右側(cè))所示.由圖中可以看出,后3組航跡點經(jīng)聚類后多數(shù)落入劃定區(qū)域范圍內(nèi).統(tǒng)計后3組高峰時段航跡點出現(xiàn)在劃定范圍區(qū)域的概率,我們將高峰時段接收的總的航跡點的集合稱作S={s1,s2···sn},落在區(qū)域范圍外的航跡點的集合稱為落在區(qū)域范圍內(nèi)概率公式為:驗證測試所得結(jié)果如表4.經(jīng)過聚類算法的處理結(jié)果表明:航跡在高峰時段的范圍比較確定,符合實際航跡分布特點.

    圖6 UTC 9點航跡聚類(左)和驗證航跡聚類(右)Fig.6 Track clustering during UTC 9(left)and track clustering verifying during UTC 9(right)

    圖7 UTC 15點航跡聚類(左)和驗證航跡聚類(右)Fig.7 Track clustering during UTC 15(left)and track clustering verifying during UTC 15(right)

    表4 概率統(tǒng)計Table 4 Probability statistics

    6 結(jié)論

    以新疆天文臺本部大樓為測試點,針對在測試點接收解析的時間、經(jīng)緯度的航跡樣本數(shù)據(jù),對數(shù)據(jù)按照時間節(jié)點初步統(tǒng)計篩選后,根據(jù)數(shù)據(jù)分布的時間特點,選用基于最小二乘法的多項式對數(shù)據(jù)進行擬合,獲得航跡樣本在一天的分布趨勢,經(jīng)計算相關(guān)系數(shù),驗證了用最小二乘法的多項式分析趨勢的可靠性.

    根據(jù)高峰時段的分布,將獲得的航跡樣本點分為兩組,一組用于測試,另一組用于驗證方法的有效性.在對第一組高峰時段數(shù)據(jù)進行處理之后獲得高峰時段航跡樣本主要分布范圍后,利用相同的方法對第二組數(shù)據(jù)進行航跡樣本的聚類處理,經(jīng)驗證航跡樣本點高達96%落在劃定的區(qū)域范圍內(nèi),驗證了該方法的有效性及實用性.通過數(shù)據(jù)擬合可有效地分析一天之中的高峰時段的分布,根據(jù)擬合結(jié)果對高峰時段航跡進行處理,劃定航跡出現(xiàn)的主要范圍區(qū)域,即可獲得一天之后高峰時段的航跡分布范圍,從而為射電天文觀測時段提供有效的參考,最大可能地避開射電天文觀測中CAS對射電天文觀測的影響和對射電望遠鏡的危害,為規(guī)劃建設(shè)中的射電天文臺址劃分寧靜區(qū)提供與相關(guān)部門交涉的依據(jù).

    [1]苗可可,王壯,程翥,等.天文研究與技術(shù),2015,12:433

    [2]Jones S R.Air Traffic Control Quaterly,2003,11:225

    [3]Robert J F,Zhang Q,Zheng Y,et al.ApJ,2005,129:2940

    [4]張群濤,劉奇,孫正文,等.科學技術(shù)與工程,2016,1:191

    [5]李耀華.基于LMS的射電天文大天線抗干擾技術(shù)的研究與實現(xiàn).長沙:國防科技大學,2011:8

    [6]吳忠祖,Haynes M P,Giovanlli R,等.天文學報,2015,56:112

    [7]Chang T C,Pen U L,Bandura K,et al.Nature,2010,466:463

    [8]Steven W E,Grant H.ApJS,2003,147:167

    [9]Maxime G,Ashok N S.IEEE Transactions on Intelligent Transportation Systems,2011,12:1511

    猜你喜歡
    射電航跡高峰
    誰能抓住“神秘天神”——快速射電暴?
    軍事文摘(2022年12期)2022-07-13 03:10:22
    射電星系
    美國的綠岸射電望遠鏡
    軍事文摘(2021年22期)2022-01-18 06:22:52
    病毒病將迎“小高峰”全方位布控巧應(yīng)對
    石慶云
    書香兩岸(2020年3期)2020-06-29 12:33:45
    夢的航跡
    青年歌聲(2019年12期)2019-12-17 06:32:32
    自適應(yīng)引導長度的無人機航跡跟蹤方法
    視覺導航下基于H2/H∞的航跡跟蹤
    世界著名射電望遠鏡
    太空探索(2016年9期)2016-07-12 10:00:01
    雨天早高峰,自在從容時。
    404 Not Found

    404 Not Found


    nginx
    一本—道久久a久久精品蜜桃钙片| 欧美成人午夜精品| 国产精品香港三级国产av潘金莲| 男女之事视频高清在线观看| 99热网站在线观看| 亚洲国产欧美网| 999久久久精品免费观看国产| 80岁老熟妇乱子伦牲交| 亚洲九九香蕉| 日韩中文字幕视频在线看片| 欧美乱码精品一区二区三区| 日本91视频免费播放| 欧美日本中文国产一区发布| 18禁黄网站禁片午夜丰满| 人人澡人人妻人| 丝袜脚勾引网站| 黄色 视频免费看| 精品亚洲乱码少妇综合久久| 美女视频免费永久观看网站| 老司机靠b影院| 精品人妻一区二区三区麻豆| 日韩 欧美 亚洲 中文字幕| 久久天堂一区二区三区四区| www.999成人在线观看| 黑丝袜美女国产一区| 亚洲黑人精品在线| 两性夫妻黄色片| 日本vs欧美在线观看视频| 国产97色在线日韩免费| 亚洲综合色网址| 国产99久久九九免费精品| 亚洲专区字幕在线| 一级a爱视频在线免费观看| 建设人人有责人人尽责人人享有的| 另类精品久久| 嫩草影视91久久| av在线老鸭窝| 蜜桃在线观看..| 国产在线免费精品| 肉色欧美久久久久久久蜜桃| 国产精品 国内视频| 国产av又大| 18禁黄网站禁片午夜丰满| 色综合欧美亚洲国产小说| 欧美久久黑人一区二区| 午夜福利影视在线免费观看| 脱女人内裤的视频| 蜜桃在线观看..| 999精品在线视频| 亚洲国产av影院在线观看| 久久精品亚洲av国产电影网| 日本av免费视频播放| 亚洲色图 男人天堂 中文字幕| 成人国语在线视频| 午夜老司机福利片| 国产97色在线日韩免费| 国产高清国产精品国产三级| 亚洲精品av麻豆狂野| 久久久久网色| 黄色 视频免费看| 日本黄色日本黄色录像| 国产av又大| 日本av手机在线免费观看| 999久久久精品免费观看国产| 美女高潮喷水抽搐中文字幕| 99热国产这里只有精品6| 窝窝影院91人妻| 久久久久久免费高清国产稀缺| 人人妻,人人澡人人爽秒播| 中文精品一卡2卡3卡4更新| 欧美在线一区亚洲| 老鸭窝网址在线观看| 久久久国产成人免费| 欧美久久黑人一区二区| 老司机在亚洲福利影院| 国产在线观看jvid| 日韩熟女老妇一区二区性免费视频| 母亲3免费完整高清在线观看| 男女高潮啪啪啪动态图| 午夜激情久久久久久久| 我的亚洲天堂| 久久久久国产一级毛片高清牌| 国产精品野战在线观看| 黑人巨大精品欧美一区二区mp4| 久久婷婷成人综合色麻豆| ponron亚洲| 一区福利在线观看| 国产三级黄色录像| 黄频高清免费视频| 免费在线观看黄色视频的| 国产探花在线观看一区二区| 久久国产精品影院| 草草在线视频免费看| 听说在线观看完整版免费高清| 色精品久久人妻99蜜桃| 女警被强在线播放| 国产精品久久电影中文字幕| 啦啦啦韩国在线观看视频| 婷婷亚洲欧美| 男女做爰动态图高潮gif福利片| 午夜福利在线观看吧| 日韩欧美国产在线观看| 久久久久久九九精品二区国产 | 国产av又大| 日日爽夜夜爽网站| 亚洲成人精品中文字幕电影| 中亚洲国语对白在线视频| 欧美久久黑人一区二区| 久久欧美精品欧美久久欧美| 亚洲国产精品成人综合色| 国产精品永久免费网站| 身体一侧抽搐| 在线观看美女被高潮喷水网站 | 在线观看日韩欧美| 18禁黄网站禁片午夜丰满| 亚洲在线自拍视频| 久久久久免费精品人妻一区二区| 久久精品91蜜桃| 啦啦啦免费观看视频1| 床上黄色一级片| 午夜福利在线观看吧| 国产精品亚洲av一区麻豆| 国产精品野战在线观看| 亚洲一区二区三区色噜噜| 精品人妻1区二区| 亚洲欧美日韩东京热| 久久国产精品影院| 成人av一区二区三区在线看| 变态另类丝袜制服| 麻豆国产97在线/欧美 | 国产精品久久电影中文字幕| 国产亚洲精品综合一区在线观看 | 桃红色精品国产亚洲av| 国产片内射在线| 777久久人妻少妇嫩草av网站| www日本黄色视频网| 久久中文看片网| 国产私拍福利视频在线观看| 国产精品一区二区三区四区久久| 亚洲av日韩精品久久久久久密| 两性夫妻黄色片| 精品第一国产精品| 不卡av一区二区三区| 波多野结衣高清无吗| 88av欧美| 亚洲五月天丁香| 成熟少妇高潮喷水视频| 超碰成人久久| 最好的美女福利视频网| 90打野战视频偷拍视频| 国产精品爽爽va在线观看网站| 国产高清videossex| 亚洲国产中文字幕在线视频| 成人永久免费在线观看视频| 国产精品亚洲一级av第二区| 五月伊人婷婷丁香| 最近视频中文字幕2019在线8| 亚洲一卡2卡3卡4卡5卡精品中文| 法律面前人人平等表现在哪些方面| 欧美成人午夜精品| 一个人免费在线观看电影 | www.熟女人妻精品国产| 99久久久亚洲精品蜜臀av| 嫩草影院精品99| 免费看日本二区| 一边摸一边做爽爽视频免费| 香蕉丝袜av| 国产高清videossex| 天天躁夜夜躁狠狠躁躁| 色精品久久人妻99蜜桃| 欧美乱码精品一区二区三区| 亚洲18禁久久av| 国产亚洲欧美98| 黄色成人免费大全| 舔av片在线| 日本在线视频免费播放| 村上凉子中文字幕在线| 好男人在线观看高清免费视频| 韩国av一区二区三区四区| 久久亚洲精品不卡| 夜夜夜夜夜久久久久| 又紧又爽又黄一区二区| 在线观看美女被高潮喷水网站 | 啦啦啦免费观看视频1| av超薄肉色丝袜交足视频| www国产在线视频色| 国产在线观看jvid| av欧美777| 欧美黑人巨大hd| 亚洲国产精品久久男人天堂| 亚洲熟女毛片儿| 亚洲熟妇中文字幕五十中出| 精品熟女少妇八av免费久了| netflix在线观看网站| 亚洲国产精品999在线| 亚洲 欧美 日韩 在线 免费| 亚洲,欧美精品.| 在线观看日韩欧美| 一本综合久久免费| 精品国产美女av久久久久小说| 一边摸一边抽搐一进一小说| 日韩大尺度精品在线看网址| 黄色 视频免费看| 久久久水蜜桃国产精品网| 非洲黑人性xxxx精品又粗又长| 国产亚洲精品一区二区www| 亚洲国产精品sss在线观看| 欧美zozozo另类| 久久久久精品国产欧美久久久| 日韩国内少妇激情av| 999久久久精品免费观看国产| 天天躁狠狠躁夜夜躁狠狠躁| 在线观看美女被高潮喷水网站 | 欧美成狂野欧美在线观看| 男女那种视频在线观看| 国产精品亚洲一级av第二区| 两个人视频免费观看高清| 精品一区二区三区av网在线观看| 久久这里只有精品中国| 又爽又黄无遮挡网站| 久久久久久免费高清国产稀缺| 欧美性猛交╳xxx乱大交人| 国产aⅴ精品一区二区三区波| 99国产极品粉嫩在线观看| 免费看十八禁软件| 在线十欧美十亚洲十日本专区| 深夜精品福利| 国产成人啪精品午夜网站| 日韩精品青青久久久久久| 最近最新免费中文字幕在线| 国产成+人综合+亚洲专区| 国产精品免费一区二区三区在线| 亚洲avbb在线观看| 天天躁狠狠躁夜夜躁狠狠躁| 午夜老司机福利片| 免费在线观看影片大全网站| 久久精品成人免费网站| 少妇熟女aⅴ在线视频| 久久久水蜜桃国产精品网| 国产久久久一区二区三区| 亚洲 国产 在线| 日韩高清综合在线| 精品福利观看| 久久久久久久午夜电影| 国产亚洲欧美98| 男插女下体视频免费在线播放| 99在线人妻在线中文字幕| 亚洲成av人片免费观看| 最近最新中文字幕大全电影3| 老司机深夜福利视频在线观看| 久久欧美精品欧美久久欧美| 国产精品精品国产色婷婷| 国产97色在线日韩免费| 久久精品亚洲精品国产色婷小说| 亚洲乱码一区二区免费版| 亚洲七黄色美女视频| 国产激情欧美一区二区| 亚洲精品色激情综合| av片东京热男人的天堂| 91老司机精品| 久久草成人影院| 老司机深夜福利视频在线观看| 成人三级黄色视频| 国产成人精品无人区| 精品一区二区三区av网在线观看| 神马国产精品三级电影在线观看 | 老司机福利观看| 成年版毛片免费区| 亚洲无线在线观看| 91麻豆av在线| 全区人妻精品视频| 久久天堂一区二区三区四区| 男插女下体视频免费在线播放| 一级毛片高清免费大全| 亚洲乱码一区二区免费版| 一级片免费观看大全| 欧美久久黑人一区二区| 日韩成人在线观看一区二区三区| 午夜亚洲福利在线播放| 久久久精品国产亚洲av高清涩受| 欧美高清成人免费视频www| 国产主播在线观看一区二区| 免费看十八禁软件| 看片在线看免费视频| 身体一侧抽搐| 91麻豆精品激情在线观看国产| 日韩有码中文字幕| 夜夜躁狠狠躁天天躁| 国产精品99久久99久久久不卡| 日韩三级视频一区二区三区| 一进一出抽搐动态| 色老头精品视频在线观看| 18禁国产床啪视频网站| 热99re8久久精品国产| 淫妇啪啪啪对白视频| 亚洲 欧美 日韩 在线 免费| 女生性感内裤真人,穿戴方法视频| 99久久精品国产亚洲精品| 成人精品一区二区免费| 无遮挡黄片免费观看| 国产伦人伦偷精品视频| 亚洲欧美精品综合久久99| 日韩中文字幕欧美一区二区| 美女高潮喷水抽搐中文字幕| 精品欧美一区二区三区在线| 精品一区二区三区视频在线观看免费| 久久香蕉国产精品| 午夜成年电影在线免费观看| 国产激情欧美一区二区| 巨乳人妻的诱惑在线观看| 俺也久久电影网| 真人做人爱边吃奶动态| av欧美777| 欧美黑人精品巨大| 一本久久中文字幕| 日日爽夜夜爽网站| 好男人电影高清在线观看| 国产精品av久久久久免费| 狂野欧美激情性xxxx| 国产精品久久久久久久电影 | 黄色成人免费大全| 哪里可以看免费的av片| 国产v大片淫在线免费观看| 欧美日韩黄片免| 天天添夜夜摸| 真人做人爱边吃奶动态| 国产精品日韩av在线免费观看| 日本一本二区三区精品| 99热只有精品国产| 国产精品一区二区免费欧美| 岛国视频午夜一区免费看| 久久久久精品国产欧美久久久| avwww免费| 亚洲第一欧美日韩一区二区三区| 欧美极品一区二区三区四区| 色在线成人网| 狂野欧美激情性xxxx| 国产三级在线视频| 夜夜夜夜夜久久久久| 亚洲精品在线观看二区| 国产97色在线日韩免费| 亚洲国产精品999在线| 夜夜夜夜夜久久久久| 不卡一级毛片| cao死你这个sao货| 日本撒尿小便嘘嘘汇集6| 中亚洲国语对白在线视频| www.熟女人妻精品国产| 日本黄色视频三级网站网址| 国产av不卡久久| 1024手机看黄色片| 亚洲av中文字字幕乱码综合| 男女下面进入的视频免费午夜| www日本黄色视频网| 在线观看66精品国产| 久久精品夜夜夜夜夜久久蜜豆 | 毛片女人毛片| 国产精品乱码一区二三区的特点| 长腿黑丝高跟| 色综合欧美亚洲国产小说| 亚洲成人久久性| 国产一区二区三区在线臀色熟女| 国产亚洲精品av在线| 欧美性猛交黑人性爽| 看片在线看免费视频| 五月玫瑰六月丁香| 精品一区二区三区av网在线观看| 变态另类丝袜制服| 国产成人精品久久二区二区免费| 欧美性猛交╳xxx乱大交人| 国产av麻豆久久久久久久| 在线观看一区二区三区| 久久久久久久午夜电影| 国产三级黄色录像| 国产区一区二久久| 亚洲va日本ⅴa欧美va伊人久久| 色在线成人网| 亚洲国产精品成人综合色| 国产av麻豆久久久久久久| 亚洲精品久久成人aⅴ小说| 国产精品美女特级片免费视频播放器 | 国产精品一及| 99久久国产精品久久久| 一边摸一边做爽爽视频免费| 亚洲狠狠婷婷综合久久图片| 精品国产乱码久久久久久男人| 在线观看美女被高潮喷水网站 | 午夜福利视频1000在线观看| 亚洲成av人片在线播放无| 免费在线观看完整版高清| 少妇的丰满在线观看| 又粗又爽又猛毛片免费看| 一本久久中文字幕| 天天一区二区日本电影三级| 女生性感内裤真人,穿戴方法视频| 成人欧美大片| 极品教师在线免费播放| 亚洲国产精品合色在线| 18美女黄网站色大片免费观看| 欧美成人午夜精品| 久久久国产成人免费| 三级国产精品欧美在线观看 | 一级毛片高清免费大全| 在线a可以看的网站| 国产精品免费一区二区三区在线| 亚洲人成伊人成综合网2020| 国产精品1区2区在线观看.| 中文字幕熟女人妻在线| 最新美女视频免费是黄的| 少妇人妻一区二区三区视频| 国产熟女午夜一区二区三区| 国产精品国产高清国产av| or卡值多少钱| 色综合婷婷激情| 久久中文字幕一级| 亚洲国产精品合色在线| 1024视频免费在线观看| 老司机靠b影院| 黄色片一级片一级黄色片| 大型黄色视频在线免费观看| 国产亚洲欧美98| 2021天堂中文幕一二区在线观| 国产精华一区二区三区| 曰老女人黄片| 美女 人体艺术 gogo| 精品午夜福利视频在线观看一区| 国产午夜福利久久久久久| 中亚洲国语对白在线视频| 国产免费男女视频| 国产精品久久久久久亚洲av鲁大| 午夜福利视频1000在线观看| 岛国在线观看网站| 亚洲狠狠婷婷综合久久图片| 黑人欧美特级aaaaaa片| av在线播放免费不卡| 亚洲成a人片在线一区二区| 每晚都被弄得嗷嗷叫到高潮| 午夜福利视频1000在线观看| 亚洲av五月六月丁香网| 性色av乱码一区二区三区2| 男女那种视频在线观看| 99热只有精品国产| 国产精品野战在线观看| 午夜亚洲福利在线播放| 小说图片视频综合网站| av福利片在线| 国产在线观看jvid| 国产又黄又爽又无遮挡在线| 亚洲18禁久久av| 一区二区三区国产精品乱码| 一个人免费在线观看电影 | 国产爱豆传媒在线观看 | 少妇被粗大的猛进出69影院| 国产99久久九九免费精品| 最近视频中文字幕2019在线8| 不卡av一区二区三区| 麻豆av在线久日| 精品国产超薄肉色丝袜足j| 99精品久久久久人妻精品| 神马国产精品三级电影在线观看 | 亚洲精品色激情综合| 国产99久久九九免费精品| av在线天堂中文字幕| 757午夜福利合集在线观看| 日韩大码丰满熟妇| 男人舔奶头视频| 精品午夜福利视频在线观看一区| 精品乱码久久久久久99久播| 国产激情欧美一区二区| а√天堂www在线а√下载| 精品不卡国产一区二区三区| 亚洲人与动物交配视频| 一a级毛片在线观看| 精品第一国产精品| 九九热线精品视视频播放| 18禁黄网站禁片免费观看直播| 久久久久精品国产欧美久久久| 欧美黑人欧美精品刺激| 亚洲av片天天在线观看| 日本 欧美在线| 免费在线观看影片大全网站| 久久精品国产99精品国产亚洲性色| av国产免费在线观看| 久久国产精品人妻蜜桃| 久久久精品欧美日韩精品| 成年版毛片免费区| av福利片在线| 国产一区二区在线观看日韩 | 麻豆国产97在线/欧美 | 1024视频免费在线观看| 91九色精品人成在线观看| 日韩三级视频一区二区三区| 九色国产91popny在线| 国产精品av视频在线免费观看| 精品日产1卡2卡| 亚洲成人中文字幕在线播放| 免费看日本二区| 亚洲国产精品sss在线观看| 免费在线观看成人毛片| 国产av在哪里看| 国产成人一区二区三区免费视频网站| 老熟妇仑乱视频hdxx| 狂野欧美激情性xxxx| 国产欧美日韩一区二区三| 成人一区二区视频在线观看| 一夜夜www| 国产精品久久久久久久电影 | 两性夫妻黄色片| 女人高潮潮喷娇喘18禁视频| 久久伊人香网站| 女生性感内裤真人,穿戴方法视频| 天天一区二区日本电影三级| 午夜精品在线福利| 国产在线观看jvid| 99热这里只有是精品50| 午夜亚洲福利在线播放| 不卡一级毛片| 99久久综合精品五月天人人| 老汉色av国产亚洲站长工具| 成人高潮视频无遮挡免费网站| 毛片女人毛片| 国产成人系列免费观看| 99热这里只有是精品50| 黄频高清免费视频| 三级国产精品欧美在线观看 | 老熟妇乱子伦视频在线观看| 日韩有码中文字幕| 国产高清激情床上av| 真人做人爱边吃奶动态| 色精品久久人妻99蜜桃| 国产精品影院久久| 中文在线观看免费www的网站 | 婷婷精品国产亚洲av| 午夜精品久久久久久毛片777| 丁香欧美五月| 啦啦啦韩国在线观看视频| 51午夜福利影视在线观看| 亚洲精品中文字幕在线视频| 人妻久久中文字幕网| 一二三四社区在线视频社区8| 亚洲欧美日韩东京热| 精品欧美国产一区二区三| 精华霜和精华液先用哪个| 91老司机精品| 国产黄色小视频在线观看| 999精品在线视频| 精品日产1卡2卡| 亚洲精品粉嫩美女一区| 国产精品一区二区精品视频观看| 国产精品美女特级片免费视频播放器 | 婷婷六月久久综合丁香| 成人精品一区二区免费| 亚洲精品久久国产高清桃花| 人人妻人人澡欧美一区二区| 亚洲av成人精品一区久久| 精品免费久久久久久久清纯| 少妇粗大呻吟视频| 欧洲精品卡2卡3卡4卡5卡区| 中文字幕熟女人妻在线| 亚洲电影在线观看av| 久久久久国产一级毛片高清牌| 哪里可以看免费的av片| 欧美精品亚洲一区二区| 日本一区二区免费在线视频| 丰满的人妻完整版| 亚洲国产日韩欧美精品在线观看 | 精品国产乱码久久久久久男人| 在线观看www视频免费| 91麻豆av在线| 国产亚洲欧美在线一区二区| 国产精品99久久99久久久不卡| 免费在线观看亚洲国产| 岛国视频午夜一区免费看| 欧美精品啪啪一区二区三区| 在线观看免费午夜福利视频| 两个人免费观看高清视频| 亚洲国产日韩欧美精品在线观看 | a级毛片在线看网站| 伊人久久大香线蕉亚洲五| www.熟女人妻精品国产| 日韩欧美免费精品| 999精品在线视频| www.精华液| 国产精品精品国产色婷婷| 亚洲国产高清在线一区二区三| 婷婷丁香在线五月| videosex国产| 中文亚洲av片在线观看爽| 久久久久国产一级毛片高清牌| 亚洲国产欧美人成| 久久国产精品影院| 少妇的丰满在线观看| 国产一区二区激情短视频| 黑人巨大精品欧美一区二区mp4| 国产单亲对白刺激| 亚洲精品美女久久av网站| 精品福利观看| 国产区一区二久久| 亚洲熟妇熟女久久| 国产精品免费视频内射| 欧美中文日本在线观看视频| 女生性感内裤真人,穿戴方法视频| 日本成人三级电影网站| 国产成人欧美在线观看| 麻豆国产97在线/欧美 | 免费人成视频x8x8入口观看| 亚洲精品在线观看二区| 怎么达到女性高潮| 午夜福利在线观看吧| 成人三级黄色视频| 一区二区三区国产精品乱码| 久久性视频一级片| 久久亚洲精品不卡| 久久精品亚洲精品国产色婷小说| 亚洲人成电影免费在线|