• <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
    夜夜骑夜夜射夜夜干| 插逼视频在线观看| 欧美激情极品国产一区二区三区 | 日日摸夜夜添夜夜爱| 国产av精品麻豆| 中国美白少妇内射xxxbb| 在线观看人妻少妇| 热re99久久精品国产66热6| 久久久久网色| 亚洲av电影在线观看一区二区三区| 一本—道久久a久久精品蜜桃钙片| 97在线人人人人妻| 春色校园在线视频观看| 国产欧美日韩一区二区三区在线 | 国产亚洲午夜精品一区二区久久| 卡戴珊不雅视频在线播放| 麻豆精品久久久久久蜜桃| 精品亚洲乱码少妇综合久久| 91精品一卡2卡3卡4卡| 看十八女毛片水多多多| 看非洲黑人一级黄片| 久久午夜综合久久蜜桃| 国产片特级美女逼逼视频| 在线亚洲精品国产二区图片欧美 | 免费观看性生交大片5| 在线 av 中文字幕| 飞空精品影院首页| 国产精品久久久久成人av| 中文精品一卡2卡3卡4更新| 丰满乱子伦码专区| 久久 成人 亚洲| 午夜91福利影院| 美女内射精品一级片tv| 一区二区三区免费毛片| 亚洲欧洲日产国产| 亚洲美女黄色视频免费看| 日韩欧美一区视频在线观看| 中文字幕av电影在线播放| 午夜福利影视在线免费观看| 极品少妇高潮喷水抽搐| 插逼视频在线观看| 精品久久蜜臀av无| av不卡在线播放| 高清毛片免费看| 99国产精品免费福利视频| 久久免费观看电影| 国产片内射在线| av国产久精品久网站免费入址| 最近中文字幕高清免费大全6| av有码第一页| 伊人久久精品亚洲午夜| 亚洲精华国产精华液的使用体验| 欧美精品高潮呻吟av久久| 国产深夜福利视频在线观看| 亚洲av二区三区四区| 精品亚洲成a人片在线观看| a级毛色黄片| 中国国产av一级| 美女主播在线视频| 你懂的网址亚洲精品在线观看| 91久久精品电影网| 国产国语露脸激情在线看| 亚洲av成人精品一二三区| av在线老鸭窝| 国产成人一区二区在线| 男女无遮挡免费网站观看| 亚洲无线观看免费| 免费av不卡在线播放| 国产av码专区亚洲av| 免费人成在线观看视频色| 飞空精品影院首页| 国产精品久久久久久久久免| 国产69精品久久久久777片| 大又大粗又爽又黄少妇毛片口| 亚洲成色77777| 日本av免费视频播放| 男女边吃奶边做爰视频| 国产成人aa在线观看| av在线老鸭窝| 午夜日本视频在线| 啦啦啦中文免费视频观看日本| 性色av一级| 亚洲欧美成人精品一区二区| 日韩av在线免费看完整版不卡| 日韩av免费高清视频| 亚洲国产色片| 日韩成人伦理影院| 日本av免费视频播放| 精品少妇久久久久久888优播| 久久精品久久久久久久性| 亚洲久久久国产精品| av在线观看视频网站免费| 国产高清国产精品国产三级| 飞空精品影院首页| 你懂的网址亚洲精品在线观看| 成人国语在线视频| 久久国内精品自在自线图片| 人妻 亚洲 视频| 久久久久久久久久久久大奶| 黑人欧美特级aaaaaa片| 久久久久久久久久久丰满| 色网站视频免费| 亚洲av男天堂| 亚洲欧美清纯卡通| 秋霞伦理黄片| 亚洲精品视频女| 国产高清国产精品国产三级| 蜜桃国产av成人99| 高清毛片免费看| 午夜激情福利司机影院| 肉色欧美久久久久久久蜜桃| 亚洲av二区三区四区| 99热全是精品| 一级爰片在线观看| 精品视频人人做人人爽| 亚洲美女搞黄在线观看| 自线自在国产av| 亚洲欧美清纯卡通| 久久99热这里只频精品6学生| 大片免费播放器 马上看| 日本爱情动作片www.在线观看| 激情五月婷婷亚洲| 狂野欧美激情性xxxx在线观看| 在线观看免费日韩欧美大片 | 欧美老熟妇乱子伦牲交| 国产成人精品一,二区| 中文字幕亚洲精品专区| a级毛片黄视频| av国产久精品久网站免费入址| 高清av免费在线| 亚洲精品aⅴ在线观看| 亚洲欧美成人精品一区二区| 午夜视频国产福利| 九九在线视频观看精品| 人妻制服诱惑在线中文字幕| 黑人高潮一二区| 伦精品一区二区三区| 久热久热在线精品观看| 日韩成人av中文字幕在线观看| 日本av免费视频播放| 欧美日本中文国产一区发布| 乱人伦中国视频| 自线自在国产av| 欧美精品亚洲一区二区| videossex国产| 一边亲一边摸免费视频| 观看美女的网站| 免费大片18禁| 久久久久久久久久人人人人人人| 天美传媒精品一区二区| 人人妻人人澡人人看| 最新的欧美精品一区二区| 精品视频人人做人人爽| 国产视频内射| 国产高清国产精品国产三级| 视频区图区小说| 国产 一区精品| 国产亚洲av片在线观看秒播厂| 啦啦啦视频在线资源免费观看| 国产亚洲最大av| 亚洲精品久久久久久婷婷小说| 免费观看在线日韩| 人人妻人人澡人人爽人人夜夜| 欧美老熟妇乱子伦牲交| 久久久久久久亚洲中文字幕| h视频一区二区三区| 极品少妇高潮喷水抽搐| 狠狠精品人妻久久久久久综合| 2018国产大陆天天弄谢| 视频在线观看一区二区三区| 三级国产精品欧美在线观看| 久久狼人影院| 水蜜桃什么品种好| 少妇 在线观看| av卡一久久| av电影中文网址| 最后的刺客免费高清国语| www.色视频.com| 91精品三级在线观看| 欧美成人午夜免费资源| 国产精品久久久久久av不卡| 啦啦啦在线观看免费高清www| 五月伊人婷婷丁香| 欧美亚洲 丝袜 人妻 在线| 一区在线观看完整版| 国产男女内射视频| 一级片'在线观看视频| 菩萨蛮人人尽说江南好唐韦庄| 日韩av免费高清视频| 一级二级三级毛片免费看| 免费观看在线日韩| 国产精品久久久久久久久免| 国产精品久久久久久精品电影小说| 国产精品一二三区在线看| 插阴视频在线观看视频| 大片免费播放器 马上看| 一级,二级,三级黄色视频| av在线观看视频网站免费| 亚洲av男天堂| 成人影院久久| 一本大道久久a久久精品| 下体分泌物呈黄色| 国产老妇伦熟女老妇高清| av电影中文网址| 我要看黄色一级片免费的| 最近的中文字幕免费完整| 久久久精品94久久精品| 国产精品三级大全| 成人18禁高潮啪啪吃奶动态图 | 一本—道久久a久久精品蜜桃钙片| 夫妻午夜视频| 丰满迷人的少妇在线观看| 午夜福利影视在线免费观看| 91成人精品电影| 亚洲无线观看免费| 我的老师免费观看完整版| 日韩av免费高清视频| 搡老乐熟女国产| 久久精品久久精品一区二区三区| 一区二区三区四区激情视频| 夜夜骑夜夜射夜夜干| 国产精品国产三级专区第一集| 日韩av不卡免费在线播放| 国产极品粉嫩免费观看在线 | 午夜福利视频在线观看免费| 久久人人爽人人片av| 内地一区二区视频在线| 精品亚洲成a人片在线观看| 精品一区二区三卡| 人成视频在线观看免费观看| av卡一久久| 久久久久久久久久久丰满| 老司机影院毛片| 满18在线观看网站| 男女高潮啪啪啪动态图| 亚洲av成人精品一区久久| 满18在线观看网站| 丰满少妇做爰视频| 日韩视频在线欧美| 免费播放大片免费观看视频在线观看| 春色校园在线视频观看| 精品少妇黑人巨大在线播放| 久久久久精品久久久久真实原创| 欧美xxxx性猛交bbbb| 日日摸夜夜添夜夜添av毛片| 成人手机av| 亚洲综合精品二区| av视频免费观看在线观看| 国产国语露脸激情在线看| 久久午夜福利片| 赤兔流量卡办理| 国产精品偷伦视频观看了| 成人漫画全彩无遮挡| 精品人妻一区二区三区麻豆| 欧美精品人与动牲交sv欧美| 国模一区二区三区四区视频| 秋霞伦理黄片| 久久久久久久久大av| 男人爽女人下面视频在线观看| 伊人久久精品亚洲午夜| freevideosex欧美| 久久人人爽人人爽人人片va| 国产69精品久久久久777片| 大话2 男鬼变身卡| 久久综合国产亚洲精品| 国产淫语在线视频| 亚洲天堂av无毛| 成年人午夜在线观看视频| 两个人免费观看高清视频| 91久久精品电影网| 热99国产精品久久久久久7| 亚洲三级黄色毛片| 99久久中文字幕三级久久日本| 99九九在线精品视频| 免费黄频网站在线观看国产| 国产色爽女视频免费观看| 亚洲精品久久成人aⅴ小说 | 18禁在线无遮挡免费观看视频| 午夜福利影视在线免费观看| 欧美成人精品欧美一级黄| 精品99又大又爽又粗少妇毛片| 成人亚洲精品一区在线观看| 草草在线视频免费看| 一本久久精品| 天天影视国产精品| 九九久久精品国产亚洲av麻豆| 一本大道久久a久久精品| 免费大片黄手机在线观看| 国产伦精品一区二区三区视频9| 国产精品国产三级国产专区5o| 国产精品99久久久久久久久| 亚洲综合色网址| 九色亚洲精品在线播放| 日本黄大片高清| 欧美另类一区| 久久精品久久久久久久性| 国产成人精品福利久久| 最近最新中文字幕免费大全7| 超碰97精品在线观看| 国产探花极品一区二区| 97在线视频观看| 另类精品久久| freevideosex欧美| 欧美日韩国产mv在线观看视频| 十八禁网站网址无遮挡| 精品午夜福利在线看| 久久精品国产亚洲av天美| av.在线天堂| 在线观看免费高清a一片| 又黄又爽又刺激的免费视频.| 成人国产麻豆网| 亚洲综合精品二区| 色94色欧美一区二区| 国产男女内射视频| 精品视频人人做人人爽| 精品国产乱码久久久久久小说| 亚洲欧美日韩另类电影网站| 九草在线视频观看| 婷婷色综合大香蕉| 久久久久精品性色| 在线观看免费视频网站a站| 一本—道久久a久久精品蜜桃钙片| 国产黄色视频一区二区在线观看| 成年美女黄网站色视频大全免费 | 久久精品熟女亚洲av麻豆精品| 超色免费av| 欧美精品亚洲一区二区| 大香蕉久久网| 天美传媒精品一区二区| av在线老鸭窝| 又大又黄又爽视频免费| 下体分泌物呈黄色| 日本黄大片高清| 黄片播放在线免费| 日韩成人伦理影院| 国产精品人妻久久久影院| 新久久久久国产一级毛片| 一级,二级,三级黄色视频| 男人爽女人下面视频在线观看| 91精品三级在线观看| 青春草亚洲视频在线观看| 国产一区有黄有色的免费视频| 久久99精品国语久久久| 免费av中文字幕在线| 久久精品久久久久久久性| 夫妻性生交免费视频一级片| 天堂8中文在线网| 中文字幕久久专区| 蜜臀久久99精品久久宅男| 黄色视频在线播放观看不卡| 人妻制服诱惑在线中文字幕| 高清不卡的av网站| 91成人精品电影| 精品国产一区二区久久| a 毛片基地| 国产成人精品无人区| 久久午夜综合久久蜜桃| 97超碰精品成人国产| 22中文网久久字幕| 美女xxoo啪啪120秒动态图| 国产精品嫩草影院av在线观看| 97在线人人人人妻| 欧美最新免费一区二区三区| 中文欧美无线码| 久久国产亚洲av麻豆专区| 精品人妻熟女毛片av久久网站| 一本—道久久a久久精品蜜桃钙片| 老熟女久久久| 91国产中文字幕| 日日爽夜夜爽网站| 久久久欧美国产精品| 精品午夜福利在线看| 亚洲国产成人一精品久久久| 高清午夜精品一区二区三区| 免费日韩欧美在线观看| 久久久久国产网址| 在线观看免费高清a一片| 午夜影院在线不卡| 亚洲av.av天堂| 国产乱人偷精品视频| 亚洲精品视频女| 亚洲精品成人av观看孕妇| av有码第一页| 久久韩国三级中文字幕| 免费看不卡的av| 欧美少妇被猛烈插入视频| 亚洲av在线观看美女高潮| 一区二区av电影网| 精品人妻熟女av久视频| 日韩av在线免费看完整版不卡| 久久午夜福利片| 久久国产精品大桥未久av| 国国产精品蜜臀av免费| 国产成人午夜福利电影在线观看| 一区在线观看完整版| 18禁动态无遮挡网站| 我的老师免费观看完整版| 亚洲国产色片| 日韩成人伦理影院| 国产精品秋霞免费鲁丝片| 精品久久久久久电影网| 亚洲人成网站在线观看播放| 久久女婷五月综合色啪小说| 18在线观看网站| 久久人人爽av亚洲精品天堂| 国产精品一区www在线观看| 又黄又爽又刺激的免费视频.| 亚洲国产日韩一区二区| 婷婷色综合大香蕉| 久久精品久久精品一区二区三区| 91成人精品电影| .国产精品久久| 日本午夜av视频| 日韩强制内射视频| 中文字幕制服av| 女性生殖器流出的白浆| 欧美人与性动交α欧美精品济南到 | 国产极品粉嫩免费观看在线 | 在线亚洲精品国产二区图片欧美 | 我要看黄色一级片免费的| 插逼视频在线观看| 高清av免费在线| 少妇的逼水好多| 日韩三级伦理在线观看| 午夜视频国产福利| 中国三级夫妇交换| 久久精品熟女亚洲av麻豆精品| 婷婷色综合大香蕉| 婷婷成人精品国产| 又大又黄又爽视频免费| 亚洲av不卡在线观看| 精品国产乱码久久久久久小说| 91久久精品国产一区二区成人| 久久国产精品大桥未久av| 成人18禁高潮啪啪吃奶动态图 | 国产一区二区三区综合在线观看 | 亚洲国产最新在线播放| 赤兔流量卡办理| 一级毛片 在线播放| 韩国高清视频一区二区三区| 又大又黄又爽视频免费| 国产一级毛片在线| 菩萨蛮人人尽说江南好唐韦庄| 亚洲精品中文字幕在线视频| 久久精品久久久久久噜噜老黄| 亚洲欧美清纯卡通| 日韩av在线免费看完整版不卡| 国产视频内射| 婷婷色麻豆天堂久久| 97精品久久久久久久久久精品| 插逼视频在线观看| 亚洲欧美中文字幕日韩二区| 国产熟女欧美一区二区| a 毛片基地| 天堂中文最新版在线下载| 亚洲国产色片| 亚洲av日韩在线播放| 在线免费观看不下载黄p国产| 亚洲av.av天堂| 亚洲av免费高清在线观看| 看十八女毛片水多多多| 亚洲精品国产av成人精品| 日韩熟女老妇一区二区性免费视频| 啦啦啦视频在线资源免费观看| 久久av网站| 免费观看无遮挡的男女| 伊人久久精品亚洲午夜| 最近手机中文字幕大全| 免费看光身美女| 热re99久久国产66热| 午夜视频国产福利| 日日爽夜夜爽网站| 国产av码专区亚洲av| av国产久精品久网站免费入址| 我的老师免费观看完整版| 亚洲综合精品二区| 亚洲精品日韩在线中文字幕| 亚洲国产色片| 丝瓜视频免费看黄片| 久久热精品热| 一本一本综合久久| 日日摸夜夜添夜夜添av毛片| 亚洲天堂av无毛| 亚洲欧美一区二区三区国产| 国产成人精品在线电影| 在线观看美女被高潮喷水网站| 18禁观看日本| 亚洲欧美精品自产自拍| 美女内射精品一级片tv| 亚洲欧美成人精品一区二区| 精品一区在线观看国产| 久久久国产欧美日韩av| 精品久久久久久电影网| 麻豆乱淫一区二区| 日韩人妻高清精品专区| 免费观看性生交大片5| 99久久人妻综合| av福利片在线| 国内精品宾馆在线| 插逼视频在线观看| 激情五月婷婷亚洲| 美女视频免费永久观看网站| 汤姆久久久久久久影院中文字幕| 久久99热这里只频精品6学生| 亚洲人成网站在线观看播放| 国产精品成人在线| 91精品国产国语对白视频| 欧美日本中文国产一区发布| 欧美精品国产亚洲| 久久精品国产a三级三级三级| 视频在线观看一区二区三区| 日本欧美国产在线视频| 亚洲欧美成人精品一区二区| 亚洲内射少妇av| 欧美丝袜亚洲另类| 少妇精品久久久久久久| 久久久亚洲精品成人影院| 欧美日韩视频高清一区二区三区二| 久久毛片免费看一区二区三区| 精品国产乱码久久久久久小说| 欧美日韩国产mv在线观看视频| 免费大片黄手机在线观看| 国产男人的电影天堂91| av黄色大香蕉| 超色免费av| 国模一区二区三区四区视频| 大又大粗又爽又黄少妇毛片口| 国产精品无大码| 九色亚洲精品在线播放| 日韩,欧美,国产一区二区三区| 一区二区av电影网| 精品99又大又爽又粗少妇毛片| 久久免费观看电影| 日韩 亚洲 欧美在线| 久久鲁丝午夜福利片| 国产成人91sexporn| 水蜜桃什么品种好| 99视频精品全部免费 在线| 一本—道久久a久久精品蜜桃钙片| 午夜激情av网站| 伦理电影大哥的女人| 丰满乱子伦码专区| 国产精品国产三级国产专区5o| 欧美最新免费一区二区三区| 精品一区在线观看国产| 少妇被粗大的猛进出69影院 | 久热这里只有精品99| 如日韩欧美国产精品一区二区三区 | 肉色欧美久久久久久久蜜桃| 欧美性感艳星| 另类精品久久| 欧美精品人与动牲交sv欧美| 国产精品熟女久久久久浪| 精品酒店卫生间| 亚洲成人手机| 尾随美女入室| 日日撸夜夜添| 国产色爽女视频免费观看| 亚洲熟女精品中文字幕| 国产69精品久久久久777片| 狠狠婷婷综合久久久久久88av| 激情五月婷婷亚洲| 99久久人妻综合| 熟女电影av网| 老司机亚洲免费影院| 久久精品国产亚洲av涩爱| 日本爱情动作片www.在线观看| av视频免费观看在线观看| av一本久久久久| 乱人伦中国视频| 26uuu在线亚洲综合色| 少妇人妻精品综合一区二区| 岛国毛片在线播放| 久久韩国三级中文字幕| 麻豆成人av视频| 最近中文字幕高清免费大全6| 大又大粗又爽又黄少妇毛片口| 国国产精品蜜臀av免费| 乱码一卡2卡4卡精品| 91午夜精品亚洲一区二区三区| 少妇人妻精品综合一区二区| 九九在线视频观看精品| 国产亚洲精品久久久com| 人妻 亚洲 视频| 亚洲精品国产色婷婷电影| 国产日韩欧美在线精品| 日韩av免费高清视频| 精品亚洲成a人片在线观看| 十分钟在线观看高清视频www| 久久久久久久国产电影| 99久久中文字幕三级久久日本| 妹子高潮喷水视频| 国产熟女午夜一区二区三区 | h视频一区二区三区| av天堂久久9| 国产淫语在线视频| 狂野欧美激情性xxxx在线观看| h视频一区二区三区| 国产日韩一区二区三区精品不卡 | 亚洲精品中文字幕在线视频| 国产精品一区二区在线观看99| 91精品伊人久久大香线蕉| 亚洲精品一二三| 欧美日韩视频高清一区二区三区二| 午夜免费观看性视频| 欧美性感艳星| 日日摸夜夜添夜夜添av毛片| 啦啦啦中文免费视频观看日本| 赤兔流量卡办理| 欧美bdsm另类| 丁香六月天网| 精品久久久久久久久亚洲| 成年av动漫网址| 免费观看a级毛片全部| 高清在线视频一区二区三区|