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

    遺傳算法在通信衛(wèi)星軌道優(yōu)化設(shè)計(jì)中的應(yīng)用*

    2020-07-16 09:22:16鄧武東
    飛控與探測(cè) 2020年3期
    關(guān)鍵詞:星座適應(yīng)度算子

    鄧武東,陳 鋒,成 飛,賴 京

    (1.上海衛(wèi)星工程研究所·上?!?01109; 2.上海航天技術(shù)研究院·上?!?01109)

    0 引 言

    中低軌增強(qiáng)通信是對(duì)常規(guī)靜止軌道衛(wèi)星通信手段的補(bǔ)充,主要有針對(duì)目標(biāo)區(qū)域的大橢圓通信和全球低軌移動(dòng)星座通信。美國(guó)“戰(zhàn)術(shù)星”-4是一種大橢圓通信衛(wèi)星,它在頻段的選擇、通信體制的設(shè)計(jì)上,能夠與其現(xiàn)役或后續(xù)規(guī)劃的UHF通信星形成配合或者在能力上進(jìn)行增強(qiáng),并且可對(duì)特定目標(biāo)區(qū)域提供一定能力的連續(xù)通信支持。構(gòu)建基于區(qū)域型應(yīng)用的星座系統(tǒng),可為目標(biāo)地區(qū)提供準(zhǔn)實(shí)時(shí)的空間服務(wù)。為了解決大橢圓通信衛(wèi)星星座設(shè)計(jì)問(wèn)題,陳鋒等提出通過(guò)粒子群算法[1]對(duì)星座的軌道參數(shù)進(jìn)行優(yōu)化,但未尋找到全天時(shí)通信星座。相對(duì)于粒子群算法多用于連續(xù)問(wèn)題的求解[2],遺傳算法(Genetic Algorithm)還可用于對(duì)離散問(wèn)題的求解,且收斂性分析方法較為成熟。在科學(xué)和工程領(lǐng)域常用的遺傳算法可用于非線性、不可微和多峰值的復(fù)雜優(yōu)化,于亮等利用遺傳算法[3]實(shí)現(xiàn)對(duì)導(dǎo)航星座軌道參數(shù)的優(yōu)化設(shè)計(jì),優(yōu)化參數(shù)為全軌道參數(shù),模型計(jì)算復(fù)雜;ASGARIMEHR M以定位精度作為評(píng)價(jià)函數(shù)優(yōu)化了導(dǎo)航星座設(shè)計(jì)[4],評(píng)價(jià)函數(shù)的計(jì)算稍顯復(fù)雜。其他智能算法在航天器任務(wù)設(shè)計(jì)中也有著廣泛的應(yīng)用[5-7],本文以局部重點(diǎn)通信地區(qū)為例,基于改進(jìn)的遺傳算法提出了一種大橢圓通信衛(wèi)星星座軌道參數(shù)的優(yōu)化方法,模型計(jì)算量小,且逐代提升效能,收斂速度快。

    1 通信衛(wèi)星軌道的初始模型

    大橢圓通信衛(wèi)星遠(yuǎn)地點(diǎn)可選擇運(yùn)行在北極方向,可以實(shí)現(xiàn)軌道運(yùn)動(dòng)的大部分時(shí)間對(duì)北半球重點(diǎn)區(qū)域建鏈通信,以西太平洋作為局部重點(diǎn)通信地區(qū),進(jìn)行大橢圓通信衛(wèi)星軌道的設(shè)計(jì)優(yōu)化。

    1.1 軌道與星座設(shè)計(jì)

    軌道設(shè)計(jì)的任務(wù)是在載荷、星座規(guī)模一定的條件下實(shí)現(xiàn)對(duì)目標(biāo)區(qū)域的24小時(shí)實(shí)時(shí)通信支持。下文給出了設(shè)計(jì)輸入和設(shè)計(jì)目標(biāo),并闡述了關(guān)鍵軌道參數(shù)的選擇過(guò)程。

    (1)輸入條件

    衛(wèi)星軌道與星座設(shè)計(jì)的輸入條件為:

    ①目標(biāo)地區(qū):西太平洋區(qū)域?yàn)槟繕?biāo)通信站;

    ②通信仰角:10°;

    (2)設(shè)計(jì)目標(biāo)

    軌道與星座設(shè)計(jì)目標(biāo)如下:

    ①對(duì)目標(biāo)地區(qū)24小時(shí)不間斷訪問(wèn);

    ②衛(wèi)星的數(shù)量最少;

    ③載荷的作用距離最短;

    (3)軌道傾角選擇

    圖1給出了近地點(diǎn)幅角隨傾角在60天的變化情況,若采用30°低傾角,近地點(diǎn)幅角60天變化增加了113°,對(duì)于軌道覆蓋性能影響巨大,不可接受,因此須采用傾角為63.43°的臨界軌道凍結(jié)近地點(diǎn)幅角。

    圖1 近地點(diǎn)幅角隨傾角變化情況Fig.1 Peripheral angle change with different inclination

    (4)軌道高度選擇

    考慮大氣阻力的影響,軌道近地點(diǎn)高度暫定為658km。橢圓軌道的回歸圈數(shù)與軌道高度以及對(duì)應(yīng)的載荷作用距離如表1所示。在載荷作用距離約束下,可根據(jù)表1選擇合適的軌道高度和回歸天數(shù)。

    初步選擇回歸圈數(shù)為6圈的天回歸軌道、遠(yuǎn)地點(diǎn)高度在12093km左右的這條軌道。

    表1 橢圓軌道高度參數(shù)

    1.2 星座初步方案及效能

    在前述1.1節(jié)中,闡明擬采用大橢圓臨界傾角回歸軌道,利用衛(wèi)星在遠(yuǎn)地點(diǎn)附近運(yùn)行速度慢、運(yùn)行時(shí)間長(zhǎng)的特點(diǎn)以及回歸特性,再配合波束掃描能力或衛(wèi)星的姿態(tài)機(jī)動(dòng)能力,實(shí)現(xiàn)對(duì)目標(biāo)區(qū)域的長(zhǎng)時(shí)間通信覆蓋。結(jié)合通信載荷能力最遠(yuǎn)約20000km的通信能力,單星軌道具體參數(shù)如表2所示。

    表2 橢圓軌道參數(shù)

    (1)單星覆蓋性能

    衛(wèi)星軌道周期約為4小時(shí),每天運(yùn)行6圈,星下點(diǎn)軌跡每天重復(fù)。如下圖2和圖3所示,對(duì)目標(biāo)區(qū)域一天能形成四個(gè)時(shí)段的覆蓋,總時(shí)間可達(dá)約7.5小時(shí)。其中連續(xù)三軌每軌理論覆蓋時(shí)間超過(guò)2小時(shí)。因此,理論上只需要四顆星,通過(guò)相位和升交點(diǎn)赤經(jīng)的配合設(shè)計(jì),即可實(shí)現(xiàn)對(duì)目標(biāo)區(qū)域的24小時(shí)連續(xù)覆蓋。

    圖2 橢圓軌道單星目標(biāo)區(qū)域覆蓋特性Fig.2 Single star’s coverage characteristics of elliptical orbit for target area

    圖3 橢圓軌道單星目標(biāo)區(qū)域覆蓋時(shí)間Fig.3 Single star’s coverage time of elliptical orbit for target area

    (2)雙星共面橢圓軌道星座覆蓋性能分析

    采用上述軌道參數(shù),配置共軌道面但相位相差180°的兩顆大橢圓通信衛(wèi)星,軌道參數(shù)如表3所示。

    表3 雙星星座配置方案

    采用雙星接力,提供對(duì)地通信能力的方式,其覆蓋形式及特性如圖4所示:

    圖4 橢圓軌道雙星目標(biāo)區(qū)域覆蓋特性Fig.4 Double star’s coverage characteristics of elliptical orbit for target area

    圖5 橢圓軌道雙星目標(biāo)區(qū)域覆蓋時(shí)間Fig.5 Double star’s coverage time of elliptical orbit for target area

    由圖5可知,相位相差180°的兩顆大橢圓通信衛(wèi)星在上述橢圓軌道上接力工作,可對(duì)地提供近15小時(shí)的準(zhǔn)連續(xù)通信能力。

    (3)四星準(zhǔn)全天通信覆蓋星座設(shè)計(jì)

    進(jìn)一步,為了實(shí)現(xiàn)準(zhǔn)全天通信要求設(shè)計(jì)大橢圓通信衛(wèi)星星座方案。橢圓軌道參數(shù)不變的前提下,布置兩個(gè)軌道面,每個(gè)軌道面上部署兩顆相位相差180°的大橢圓通信衛(wèi)星,軌道參數(shù)如表4所示。覆蓋特性如圖6及圖7所示,可實(shí)現(xiàn)準(zhǔn)全天的通信覆蓋。

    表4 四星星座配置方案

    圖6 橢圓軌道四星目標(biāo)區(qū)域覆蓋特性Fig.6 Four star’s coverage characteristics of elliptical orbit for target area

    圖7 橢圓軌道四星目標(biāo)區(qū)域覆蓋時(shí)間Fig.7 Four star’s coverage time of elliptical orbit for target area

    根據(jù)上述軌道、星座仿真及覆蓋特性分析形成如下結(jié)論:

    ①單星橢圓軌道每天可提供3~4次通信覆蓋窗口,最多對(duì)目標(biāo)區(qū)域提供2.5小時(shí)的連續(xù)覆蓋通信時(shí)間;

    ②橢圓軌道兩星共面接力通信每天可提供超過(guò)12小時(shí)的準(zhǔn)連續(xù)通信覆蓋能力,在同一個(gè)軌道面上增加衛(wèi)星數(shù)量對(duì)提高連續(xù)覆蓋能力貢獻(xiàn)不大;

    ③采用兩個(gè)軌道面,每個(gè)軌道面上部署兩顆相位相差180°的大橢圓通信衛(wèi)星,可實(shí)現(xiàn)準(zhǔn)全天連續(xù)對(duì)目標(biāo)區(qū)域的通信覆蓋能力。

    2 方案選擇與優(yōu)化

    上述四星方案的升交點(diǎn)赤經(jīng)和相位的布置還不是最優(yōu)的,接下來(lái)通過(guò)增加了保優(yōu)算子的改進(jìn)的遺傳算法尋找到最優(yōu)布置,并且按照該布置方案設(shè)計(jì)24小時(shí)無(wú)縫覆蓋的通信星座。

    2.1 待優(yōu)化參數(shù)

    在高度、偏心率、傾角、近地點(diǎn)幅角給定的情況下,待優(yōu)化參數(shù)包括兩軌道面升交點(diǎn)赤經(jīng),四顆星相位,將待優(yōu)化參數(shù)作為遺傳算法編碼對(duì)象[8],如表5所示。

    表5 遺傳算法待優(yōu)化參數(shù)

    2.2 遺傳算法優(yōu)化方法

    遺傳算法是一種模擬自然界進(jìn)化過(guò)程的優(yōu)化算法。它利用編碼技術(shù)將對(duì)象轉(zhuǎn)換成自然界中個(gè)體的染色體形式,通過(guò)設(shè)計(jì)相應(yīng)的遺傳算子來(lái)操作由這些個(gè)體組成的種群,從而達(dá)到模擬生物種群進(jìn)化過(guò)程的目的。本文將采用遺傳算法對(duì)星座的布置進(jìn)行優(yōu)化求解。

    2.2.1 適應(yīng)度函數(shù)

    類似文獻(xiàn)[9]基于指標(biāo)的多目標(biāo)搜索方法,本文以適應(yīng)度函數(shù)作為評(píng)價(jià)指標(biāo)。適應(yīng)度函數(shù)是遺傳算法收斂性和穩(wěn)定性的關(guān)鍵因素之一,對(duì)最終的優(yōu)化結(jié)果也有重要影響。因此,本文以通信覆蓋的總時(shí)間作為適應(yīng)度函數(shù)。

    f=ttotal

    2.2.2 遺傳算子

    常規(guī)遺傳算子通常包括選擇算子、交叉算子和變異算子[10-13]。為提高算法收斂性和求解的全局最優(yōu)性,本文還提出了保優(yōu)算子。下面將針對(duì)這四個(gè)算子進(jìn)行設(shè)計(jì)。

    (1)選擇算子

    通過(guò)選擇算子選出種群中較優(yōu)的個(gè)體形成一個(gè)新的種群,可使得種群中的個(gè)體不斷逼近最優(yōu)解。選擇算子可采用常用的輪盤算法,設(shè)定個(gè)體被選中的概率如下,

    (1)

    式中fj為個(gè)體j的適應(yīng)度值,n為種群數(shù)目。

    在上述選擇方式下,個(gè)體的適應(yīng)度值越大,被選擇的概率也越大,則其基因構(gòu)造被遺傳到下一代的可能性越大。這種選擇方式使得種群中的個(gè)體都存在被選中的概率,保證了較小個(gè)體的遺傳性。

    (2)交叉算子

    交叉是指按交叉概率從種群中選擇部分個(gè)體,通過(guò)交換個(gè)體部分基因形成新的個(gè)體,從而形成新的種群。交叉算子可采用雙點(diǎn)交叉算子。由于個(gè)體采用實(shí)數(shù)編碼,所以交叉操作采用實(shí)數(shù)交叉法,如式(2)所示,第k個(gè)染色體ak和第l個(gè)染色體al在j位的交叉操作方法為:

    akj=akj(1-b)+aljbalj=alj(1-b)+akjb

    (2)

    其中akj代表第k個(gè)染色體ak在j位的值,b是[0,1]區(qū)間的隨機(jī)數(shù)。

    (3)變異算子

    變異算子模擬的是生物進(jìn)化過(guò)程中的基因突變現(xiàn)象,具體操作為:令變異概率為Pc,則變異的節(jié)點(diǎn)個(gè)數(shù)為[nSPc],S為單個(gè)個(gè)體的節(jié)點(diǎn)數(shù);根據(jù)生成的隨機(jī)數(shù)r決定變異方向,選中某些個(gè)體上的某些節(jié)點(diǎn),然后用新值替代這些節(jié)點(diǎn)上的原值。新值的產(chǎn)生過(guò)程如式3所示。變異算子維持了群體的多樣性,有利于防止早熟現(xiàn)象的出現(xiàn)。

    aij=aij+(amax-aij)kb,r>0aij=aij+(amin-aij)kb,r<0

    (3)

    其中aij代表第i個(gè)染色體ai在j位的值,amax和amin代表a的上下邊界值,k是變異系數(shù),b是[0,1]區(qū)間的隨機(jī)數(shù),r是[-1,1]區(qū)間的隨機(jī)數(shù)。

    (4)保優(yōu)算子

    通過(guò)上述三個(gè)算子的操作后,種群將由原來(lái)的父代變成性能更為優(yōu)異的子代。但是仍無(wú)法保證父代種群中最優(yōu)秀的個(gè)體能全部被遺傳到子代種群,從而降低了獲得最優(yōu)解的收斂速度,甚至導(dǎo)致最后尋求到的最優(yōu)解只是局部最優(yōu),而非全局最優(yōu)。文獻(xiàn)[14-15]的常規(guī)選擇算法都陷入了局部最優(yōu)的情況。因此,本文提出了一種具有“保優(yōu)”思想的保優(yōu)算子。其基本思想是:完全保留父代種群中適應(yīng)度值最大的部分個(gè)體至子代種群。具體操作方式為:分別對(duì)父代種群和子代種群中的所有個(gè)體按適應(yīng)度值的大小降序排列,若保優(yōu)率為Pa(大約為5%~10%),則父代種群中適應(yīng)度值排在前[nPa]個(gè)的個(gè)體將完全取代子代種群中適應(yīng)度值排在后[nPa]個(gè)的個(gè)體,形成一個(gè)性能更為優(yōu)異的子代,其中[]表示取整,不同于文獻(xiàn)中只保留父代最優(yōu)個(gè)體的保優(yōu)方法[16],本文的方法保優(yōu)率可調(diào),更為靈活。

    2.2.3 算法流程

    本文算法的計(jì)算步驟如下:

    (1)生成初始化種群;

    (2)根據(jù)設(shè)計(jì)的適應(yīng)度函數(shù),計(jì)算該代種群中所有個(gè)體的適應(yīng)度;

    (3)判斷是否滿足終止條件,若不滿足則到步驟(4),若滿足則到步驟(5);

    (4)根據(jù)交叉概率和變異概率,經(jīng)過(guò)選擇算子、交叉算子、變異算子和保優(yōu)算子產(chǎn)生下一代種群,然后返回步驟(3);

    (5)選出最后一代種群中的最優(yōu)個(gè)體,輸出相應(yīng)的最優(yōu)解,算法終止。

    圖8 遺傳算法流程圖Fig.8 Flowchart of genetic algorithm

    2.3 優(yōu)化結(jié)果

    為驗(yàn)證算法的正確性和有效性,本文運(yùn)用MATLAB軟件對(duì)所提出的算法進(jìn)行了仿真。對(duì)2.1節(jié)的星座中的8個(gè)待優(yōu)化參數(shù)進(jìn)行優(yōu)化求解,結(jié)果如表6所示。

    表6 遺傳算法優(yōu)化結(jié)果

    2.1節(jié)四星原始方案相對(duì)相位為0-180-0-180,遺傳算法優(yōu)化后的方案為0-180-90-270。從圖9的遺傳算法優(yōu)化過(guò)程發(fā)現(xiàn),適應(yīng)度函數(shù)逐代提升。優(yōu)化前,5天可提供約392894.850秒時(shí)長(zhǎng)的通信時(shí)間,可通信時(shí)間占比90.948%;優(yōu)化后,5天可提供約428721s的通信時(shí)長(zhǎng),可通信時(shí)間占比99.24%,且間隔次數(shù)較少。圖10為1天6圈回歸四星組網(wǎng)優(yōu)化前后的效能對(duì)比??傻贸鼋Y(jié)論:兩軌道面升交點(diǎn)赤經(jīng)差180°較為合適,相對(duì)相位0-180-90-270比原始設(shè)計(jì)0-180-0-180要優(yōu)。

    圖9 迭代過(guò)程適應(yīng)度變化曲線Fig.9 Fitness curve of iteration process

    (a) 優(yōu)化前通信時(shí)間圖

    (b) 優(yōu)化后通信時(shí)間圖 圖10 優(yōu)化前后通信時(shí)間對(duì)比Fig.10 Comparison of communication time before and after optimization

    3 尋找實(shí)時(shí)通信星座

    按照前述遺傳算法優(yōu)化得到的結(jié)論,擬采用相對(duì)相位0-180-90-270的四星組網(wǎng)方案;通過(guò)提高軌道高度的方法,選用表1中第二行的軌道參數(shù)方案,可增加通信時(shí)間。

    3.1 軌道及星座參數(shù)

    衛(wèi)星標(biāo)稱軌道設(shè)計(jì)為:

    軌道半長(zhǎng)軸:a=13517.69324 km

    近地點(diǎn)高度:hp=658 km

    遠(yuǎn)地點(diǎn)高度:ha=13621.1 km

    最大作用距離:18954.92939 km

    軌道傾角:i=63.43°

    偏心率:e= 0.479487

    近地點(diǎn)幅角:ω= 270°

    軌道平周期:T=260.68354 min

    回歸周期:2 天

    回歸總?cè)?shù):11

    按照遺傳算法優(yōu)化后的結(jié)論,即180°的升交點(diǎn)赤經(jīng)差和0-180-90-270的相對(duì)相位方案,配置2個(gè)軌道面,每個(gè)軌道面的參數(shù)及衛(wèi)星數(shù)量見(jiàn)表7,四星組網(wǎng)示意見(jiàn)圖11。

    表7 實(shí)時(shí)通信星座配置方案

    圖11 實(shí)時(shí)通信星座配置Fig.11 Real-time communication constellation configuration

    3.2 星座效能

    仿真了5天的通信情況如圖12所示,結(jié)果表明可為目標(biāo)區(qū)域提供24小時(shí)實(shí)時(shí)通信服務(wù)。

    圖12 實(shí)時(shí)通信星座可通信時(shí)間Fig.12 Communication time of Real-time communication constellation

    4 結(jié) 論

    本文介紹了國(guó)外戰(zhàn)術(shù)衛(wèi)星的發(fā)展情況、目標(biāo)區(qū)域?qū)崟r(shí)通信的實(shí)現(xiàn)方法、衛(wèi)星星座與軌道設(shè)計(jì)優(yōu)化方案,對(duì)類似任務(wù)衛(wèi)星軌道與星座設(shè)計(jì)具有一定的借鑒意義。通過(guò)改進(jìn)的遺傳算法對(duì)軌道六根數(shù)中的升交點(diǎn)赤經(jīng)、相位兩個(gè)參數(shù),共8個(gè)變量進(jìn)行了尋優(yōu)計(jì)算。經(jīng)多次獨(dú)立優(yōu)化過(guò)程驗(yàn)證,證明了優(yōu)化結(jié)果是一種較優(yōu)的設(shè)計(jì)方案;基于此進(jìn)一步通過(guò)軌道高度提高增加覆蓋能力,尋找到了可實(shí)時(shí)覆蓋的星座設(shè)計(jì),滿足目標(biāo)區(qū)域?qū)崟r(shí)通信的要求。

    猜你喜歡
    星座適應(yīng)度算子
    改進(jìn)的自適應(yīng)復(fù)制、交叉和突變遺傳算法
    擬微分算子在Hp(ω)上的有界性
    各向異性次Laplace算子和擬p-次Laplace算子的Picone恒等式及其應(yīng)用
    一類Markov模算子半群與相應(yīng)的算子值Dirichlet型刻畫
    星座
    健康女性(2017年3期)2017-04-27 21:34:31
    12星座之我愛(ài)洗澡
    星座
    都市麗人(2016年12期)2016-12-14 09:54:15
    Roper-Suffridge延拓算子與Loewner鏈
    星座
    都市麗人(2016年7期)2016-08-03 10:06:05
    基于空調(diào)導(dǎo)風(fēng)板成型工藝的Kriging模型適應(yīng)度研究
    最后的刺客免费高清国语| 内射极品少妇av片p| 亚洲第一区二区三区不卡| 日本色播在线视频| 最近视频中文字幕2019在线8| 丰满少妇做爰视频| 一卡2卡三卡四卡精品乱码亚洲| 亚洲人成网站在线观看播放| 男的添女的下面高潮视频| 国产伦精品一区二区三区视频9| 国产黄色小视频在线观看| 非洲黑人性xxxx精品又粗又长| 欧美97在线视频| 在线a可以看的网站| av免费观看日本| 欧美性猛交╳xxx乱大交人| 日韩 亚洲 欧美在线| 99久久成人亚洲精品观看| 中文资源天堂在线| 三级经典国产精品| 亚洲中文字幕日韩| 免费看光身美女| 亚洲欧美成人精品一区二区| 亚洲美女搞黄在线观看| 久久精品国产亚洲av天美| 在线播放国产精品三级| 看黄色毛片网站| 中文亚洲av片在线观看爽| 亚洲精品日韩在线中文字幕| 精品一区二区免费观看| 国产极品精品免费视频能看的| 日韩欧美精品v在线| www.色视频.com| 免费大片18禁| 在线免费十八禁| 99久久成人亚洲精品观看| 日韩一区二区视频免费看| 亚洲欧美成人精品一区二区| 搞女人的毛片| 熟女电影av网| 久久欧美精品欧美久久欧美| 国产精品av视频在线免费观看| 九色成人免费人妻av| 欧美性猛交╳xxx乱大交人| 天天躁日日操中文字幕| 99在线视频只有这里精品首页| 麻豆av噜噜一区二区三区| 精品99又大又爽又粗少妇毛片| 又爽又黄a免费视频| 国产又黄又爽又无遮挡在线| 久久久久久久久中文| 亚洲第一区二区三区不卡| 欧美日韩精品成人综合77777| 麻豆久久精品国产亚洲av| 午夜福利成人在线免费观看| 国产亚洲最大av| 国产三级中文精品| 网址你懂的国产日韩在线| 国产精品永久免费网站| 亚洲av一区综合| 婷婷色麻豆天堂久久 | 免费av不卡在线播放| 亚洲欧美精品综合久久99| 日本免费a在线| 免费看a级黄色片| 国产精品久久久久久久电影| 久久久久久久亚洲中文字幕| 久久精品夜色国产| 久久久久九九精品影院| 只有这里有精品99| 在线观看一区二区三区| or卡值多少钱| 国产精品嫩草影院av在线观看| 狂野欧美白嫩少妇大欣赏| 99久国产av精品| 国产久久久一区二区三区| 国产综合懂色| 欧美另类亚洲清纯唯美| 久久久a久久爽久久v久久| 精品久久久久久久末码| 午夜福利在线在线| 国产伦一二天堂av在线观看| 欧美三级亚洲精品| 亚洲精品成人久久久久久| 天堂√8在线中文| 卡戴珊不雅视频在线播放| 少妇人妻精品综合一区二区| 国产综合懂色| 欧美激情久久久久久爽电影| 亚洲av中文av极速乱| 久久精品91蜜桃| 国产色婷婷99| 性色avwww在线观看| 少妇猛男粗大的猛烈进出视频 | 麻豆成人午夜福利视频| 日韩欧美精品免费久久| 国产av不卡久久| 波野结衣二区三区在线| 国产乱人视频| 国产免费一级a男人的天堂| 久久久午夜欧美精品| 色噜噜av男人的天堂激情| 亚洲综合色惰| 国国产精品蜜臀av免费| 舔av片在线| 男女视频在线观看网站免费| 日本一二三区视频观看| 在现免费观看毛片| 午夜精品在线福利| 狂野欧美白嫩少妇大欣赏| 国产精品永久免费网站| 亚洲精华国产精华液的使用体验| 国产精品久久久久久久电影| 深夜a级毛片| av又黄又爽大尺度在线免费看 | 久久久国产成人精品二区| 十八禁国产超污无遮挡网站| 国产成人一区二区在线| 久久精品91蜜桃| 国内揄拍国产精品人妻在线| 日本wwww免费看| a级一级毛片免费在线观看| 国产大屁股一区二区在线视频| 国产一区二区在线av高清观看| 免费观看性生交大片5| 一本久久精品| 日韩亚洲欧美综合| 国产91av在线免费观看| 日本一本二区三区精品| 伦理电影大哥的女人| 国产三级中文精品| 亚州av有码| 91在线精品国自产拍蜜月| 国产麻豆成人av免费视频| 中文字幕人妻熟人妻熟丝袜美| 麻豆成人av视频| 一区二区三区高清视频在线| 热99在线观看视频| 我要搜黄色片| 高清在线视频一区二区三区 | 亚洲熟妇中文字幕五十中出| 午夜激情福利司机影院| 国产黄片美女视频| 在线免费观看的www视频| 在线免费十八禁| 国产伦一二天堂av在线观看| 男女下面进入的视频免费午夜| videossex国产| 国产色婷婷99| 亚洲欧美日韩无卡精品| 性色avwww在线观看| 免费一级毛片在线播放高清视频| 91aial.com中文字幕在线观看| 亚洲精品影视一区二区三区av| 夜夜爽夜夜爽视频| 97人妻精品一区二区三区麻豆| 日韩高清综合在线| 国产成人午夜福利电影在线观看| 久久精品久久久久久噜噜老黄 | 欧美又色又爽又黄视频| 免费搜索国产男女视频| 真实男女啪啪啪动态图| 久久久午夜欧美精品| 卡戴珊不雅视频在线播放| 久久久精品大字幕| 欧美日韩在线观看h| 女人十人毛片免费观看3o分钟| 国产久久久一区二区三区| 精品一区二区三区人妻视频| 亚洲av一区综合| videos熟女内射| 91午夜精品亚洲一区二区三区| av在线播放精品| 男女那种视频在线观看| 自拍偷自拍亚洲精品老妇| 久久午夜福利片| 99久久精品一区二区三区| 3wmmmm亚洲av在线观看| 精品一区二区三区视频在线| 亚洲欧美日韩东京热| 色吧在线观看| 99热这里只有是精品在线观看| 国语自产精品视频在线第100页| 中文精品一卡2卡3卡4更新| 亚洲av二区三区四区| 亚洲av.av天堂| 日本欧美国产在线视频| 最近视频中文字幕2019在线8| 欧美变态另类bdsm刘玥| 久久亚洲国产成人精品v| 国产色爽女视频免费观看| 久久久久久久亚洲中文字幕| 联通29元200g的流量卡| 精品久久久噜噜| 亚洲国产精品久久男人天堂| 国产乱来视频区| 女的被弄到高潮叫床怎么办| 熟女人妻精品中文字幕| 好男人视频免费观看在线| 综合色av麻豆| 国内揄拍国产精品人妻在线| 日韩精品青青久久久久久| 欧美日韩一区二区视频在线观看视频在线 | 一二三四中文在线观看免费高清| 最近2019中文字幕mv第一页| 欧美xxxx黑人xx丫x性爽| 国产精品国产三级国产专区5o | 最新中文字幕久久久久| 两个人视频免费观看高清| 高清毛片免费看| 一个人看的www免费观看视频| 日日啪夜夜撸| 高清av免费在线| 男插女下体视频免费在线播放| 女人被狂操c到高潮| 亚洲图色成人| 天堂av国产一区二区熟女人妻| 午夜久久久久精精品| 最近最新中文字幕大全电影3| 又爽又黄a免费视频| 亚洲天堂国产精品一区在线| 久久草成人影院| 亚洲精品影视一区二区三区av| 国产黄片视频在线免费观看| 免费观看的影片在线观看| 极品教师在线视频| 亚洲国产成人一精品久久久| 麻豆成人av视频| 少妇被粗大猛烈的视频| 亚洲高清免费不卡视频| 丝袜喷水一区| ponron亚洲| 丰满乱子伦码专区| 青春草国产在线视频| 国产精品乱码一区二三区的特点| 久久久久精品久久久久真实原创| 性插视频无遮挡在线免费观看| 老女人水多毛片| videos熟女内射| 午夜福利网站1000一区二区三区| 国产精品1区2区在线观看.| 国产免费男女视频| 伊人久久精品亚洲午夜| 神马国产精品三级电影在线观看| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 亚洲图色成人| 一级黄色大片毛片| 精品一区二区免费观看| 国产亚洲av嫩草精品影院| 中文字幕久久专区| 国产极品天堂在线| 18禁在线播放成人免费| 久久精品国产99精品国产亚洲性色| 看免费成人av毛片| 亚洲综合色惰| 久久久国产成人免费| 中文字幕av在线有码专区| 91aial.com中文字幕在线观看| 69av精品久久久久久| 国产亚洲av片在线观看秒播厂 | 内射极品少妇av片p| 免费播放大片免费观看视频在线观看 | 91精品国产九色| 国产爱豆传媒在线观看| 国产在视频线精品| h日本视频在线播放| 99九九线精品视频在线观看视频| 99久久精品一区二区三区| 日日撸夜夜添| 女的被弄到高潮叫床怎么办| 汤姆久久久久久久影院中文字幕 | 久久精品久久久久久久性| 亚洲欧美精品自产自拍| 亚洲高清免费不卡视频| 黄色配什么色好看| 在线观看美女被高潮喷水网站| 青春草视频在线免费观看| 色综合亚洲欧美另类图片| 日本黄色片子视频| 亚洲精品日韩在线中文字幕| 能在线免费观看的黄片| 亚洲av日韩在线播放| 久久午夜福利片| 亚洲欧美日韩东京热| 日韩强制内射视频| 大又大粗又爽又黄少妇毛片口| 草草在线视频免费看| 女人被狂操c到高潮| 亚洲精品乱码久久久久久按摩| 亚洲精品国产av成人精品| 黄色配什么色好看| 嫩草影院新地址| kizo精华| 久久久久久国产a免费观看| 日本一本二区三区精品| 国产伦精品一区二区三区视频9| 国产单亲对白刺激| 汤姆久久久久久久影院中文字幕 | 99久国产av精品| 日产精品乱码卡一卡2卡三| 狠狠狠狠99中文字幕| 国产综合懂色| 内地一区二区视频在线| 中文字幕人妻熟人妻熟丝袜美| 久久精品国产亚洲av涩爱| 久久精品人妻少妇| 偷拍熟女少妇极品色| 韩国高清视频一区二区三区| 一个人看视频在线观看www免费| 精品一区二区免费观看| 乱系列少妇在线播放| 亚洲欧美日韩东京热| 最近最新中文字幕免费大全7| 亚洲成人久久爱视频| 国内精品一区二区在线观看| 最近最新中文字幕免费大全7| 天天躁日日操中文字幕| 中国国产av一级| 国产人妻一区二区三区在| 亚洲欧美成人精品一区二区| 只有这里有精品99| 日日干狠狠操夜夜爽| 女的被弄到高潮叫床怎么办| 高清午夜精品一区二区三区| 亚洲精品国产成人久久av| 午夜福利成人在线免费观看| 欧美三级亚洲精品| 美女xxoo啪啪120秒动态图| 亚洲国产精品久久男人天堂| 成人毛片a级毛片在线播放| 日本免费a在线| 欧美一级a爱片免费观看看| 亚洲欧美精品专区久久| 一本久久精品| 91狼人影院| 美女黄网站色视频| 国产精品国产三级国产专区5o | 国产精品爽爽va在线观看网站| 精品99又大又爽又粗少妇毛片| 欧美精品一区二区大全| 亚洲国产精品成人综合色| 一级爰片在线观看| 亚洲美女视频黄频| 亚洲精品乱码久久久v下载方式| 欧美日韩精品成人综合77777| 久久久国产成人精品二区| 亚洲精品,欧美精品| 午夜激情福利司机影院| av在线观看视频网站免费| 精品无人区乱码1区二区| 男女啪啪激烈高潮av片| 中文字幕精品亚洲无线码一区| 一区二区三区高清视频在线| 国产麻豆成人av免费视频| 在线观看av片永久免费下载| 精品无人区乱码1区二区| 91午夜精品亚洲一区二区三区| 精品人妻熟女av久视频| 亚洲精品日韩在线中文字幕| 麻豆av噜噜一区二区三区| 国产三级在线视频| 国产精品久久久久久久久免| 看非洲黑人一级黄片| 中文字幕亚洲精品专区| 日本wwww免费看| av线在线观看网站| 免费电影在线观看免费观看| 精品久久久噜噜| 国产高清视频在线观看网站| 久久这里有精品视频免费| 国国产精品蜜臀av免费| 日韩成人伦理影院| 搡老妇女老女人老熟妇| 少妇高潮的动态图| 欧美性猛交╳xxx乱大交人| 久久精品国产鲁丝片午夜精品| 成人av在线播放网站| 99久久人妻综合| 亚洲三级黄色毛片| .国产精品久久| 日本午夜av视频| 淫秽高清视频在线观看| 18禁在线无遮挡免费观看视频| 久久人人爽人人片av| 别揉我奶头 嗯啊视频| 日韩一区二区视频免费看| 国产成人精品久久久久久| 亚洲精品日韩av片在线观看| 国产一区二区在线av高清观看| 97热精品久久久久久| 又爽又黄无遮挡网站| 中国美白少妇内射xxxbb| 国产三级中文精品| 99热网站在线观看| 男女视频在线观看网站免费| 非洲黑人性xxxx精品又粗又长| 免费一级毛片在线播放高清视频| 亚洲国产高清在线一区二区三| 久久这里只有精品中国| 99九九线精品视频在线观看视频| 人妻少妇偷人精品九色| 久久精品国产亚洲av涩爱| 国产成人精品婷婷| 亚洲国产精品国产精品| 村上凉子中文字幕在线| 乱系列少妇在线播放| 亚洲中文字幕日韩| 18禁在线无遮挡免费观看视频| 国产午夜精品一二区理论片| 久久亚洲精品不卡| 国产成人aa在线观看| 欧美日本视频| 97人妻精品一区二区三区麻豆| 最近中文字幕高清免费大全6| 99热6这里只有精品| 亚洲三级黄色毛片| 1000部很黄的大片| 日韩强制内射视频| 精品人妻熟女av久视频| 亚洲av成人精品一区久久| 国产老妇女一区| 九九久久精品国产亚洲av麻豆| 九九热线精品视视频播放| 久久久国产成人免费| 久久这里有精品视频免费| 久久婷婷人人爽人人干人人爱| 精品免费久久久久久久清纯| 久久久久久久国产电影| 波多野结衣高清无吗| 好男人在线观看高清免费视频| 在线a可以看的网站| 91久久精品国产一区二区成人| 麻豆成人av视频| 舔av片在线| 晚上一个人看的免费电影| 最近中文字幕2019免费版| 天美传媒精品一区二区| 午夜精品国产一区二区电影 | 国内揄拍国产精品人妻在线| 国产麻豆成人av免费视频| 亚洲经典国产精华液单| 免费看美女性在线毛片视频| 观看免费一级毛片| 搡女人真爽免费视频火全软件| 色视频www国产| 国产乱人视频| 欧美色视频一区免费| 久久久国产成人免费| 免费无遮挡裸体视频| 国产黄片视频在线免费观看| 免费av观看视频| 国产精品精品国产色婷婷| 亚洲av中文字字幕乱码综合| 99在线人妻在线中文字幕| 久久韩国三级中文字幕| 欧美成人精品欧美一级黄| 天天躁夜夜躁狠狠久久av| 亚洲高清免费不卡视频| 欧美激情国产日韩精品一区| 国产亚洲最大av| 精品久久久久久久人妻蜜臀av| 亚洲av成人av| 乱码一卡2卡4卡精品| 亚洲国产精品专区欧美| 亚洲国产成人一精品久久久| 国产午夜精品久久久久久一区二区三区| 亚洲人与动物交配视频| 亚洲av日韩在线播放| 国产在视频线精品| 精品国产一区二区三区久久久樱花 | 两个人的视频大全免费| 日本一二三区视频观看| 美女脱内裤让男人舔精品视频| 一个人观看的视频www高清免费观看| 国产精品99久久久久久久久| 91久久精品国产一区二区三区| 午夜视频国产福利| 国产精品日韩av在线免费观看| 久久这里有精品视频免费| 国产精品电影一区二区三区| 亚洲怡红院男人天堂| 亚洲av不卡在线观看| 村上凉子中文字幕在线| 色网站视频免费| 18禁在线播放成人免费| 搞女人的毛片| 国产黄片视频在线免费观看| 精品酒店卫生间| 卡戴珊不雅视频在线播放| .国产精品久久| 亚洲va在线va天堂va国产| 老司机影院毛片| 国产在视频线精品| 精品一区二区三区人妻视频| 麻豆成人午夜福利视频| 欧美一区二区精品小视频在线| 能在线免费观看的黄片| 日韩av不卡免费在线播放| 非洲黑人性xxxx精品又粗又长| 国产极品天堂在线| 欧美丝袜亚洲另类| ponron亚洲| 精品无人区乱码1区二区| 2021天堂中文幕一二区在线观| 亚洲国产成人一精品久久久| 99热网站在线观看| 精品久久久噜噜| 综合色av麻豆| 国产伦一二天堂av在线观看| 欧美激情在线99| 色网站视频免费| 精品一区二区三区视频在线| 亚洲精品亚洲一区二区| 午夜免费激情av| 国产一级毛片七仙女欲春2| 国产精品久久久久久精品电影小说 | 可以在线观看毛片的网站| 99在线人妻在线中文字幕| 欧美高清成人免费视频www| 一级毛片电影观看 | 免费观看在线日韩| 亚洲国产最新在线播放| 三级经典国产精品| 亚洲av熟女| 日韩国内少妇激情av| 免费大片18禁| 尾随美女入室| 中国国产av一级| 免费av观看视频| 久久久精品欧美日韩精品| 中文字幕av在线有码专区| 少妇的逼好多水| 成人综合一区亚洲| 亚洲国产欧洲综合997久久,| 男女那种视频在线观看| 久久久久精品久久久久真实原创| 亚洲图色成人| 国产伦精品一区二区三区视频9| 精品久久久久久久久久久久久| 久久久久精品久久久久真实原创| 久久精品国产99精品国产亚洲性色| 久久久亚洲精品成人影院| 免费黄色在线免费观看| 国产精品无大码| 中文字幕制服av| 国产精品无大码| 麻豆av噜噜一区二区三区| 国产 一区精品| 亚洲久久久久久中文字幕| 亚洲国产精品专区欧美| av在线老鸭窝| 又粗又硬又长又爽又黄的视频| 精品久久久久久久人妻蜜臀av| 最近视频中文字幕2019在线8| 内射极品少妇av片p| 国产成人a∨麻豆精品| 99视频精品全部免费 在线| 六月丁香七月| 寂寞人妻少妇视频99o| 最近2019中文字幕mv第一页| 免费看a级黄色片| 亚洲国产成人一精品久久久| 色综合亚洲欧美另类图片| 精品久久久噜噜| 久久久久久大精品| 高清视频免费观看一区二区 | 久久久久九九精品影院| 日韩欧美 国产精品| 青春草国产在线视频| 男女下面进入的视频免费午夜| videos熟女内射| 国产精品三级大全| 麻豆成人av视频| 国产精品久久久久久av不卡| 欧美日韩一区二区视频在线观看视频在线 | 久久亚洲精品不卡| 国产精品人妻久久久久久| 男人的好看免费观看在线视频| 日本色播在线视频| 美女被艹到高潮喷水动态| 一区二区三区免费毛片| 91aial.com中文字幕在线观看| 欧美97在线视频| 嫩草影院入口| 日本与韩国留学比较| 十八禁国产超污无遮挡网站| 女人十人毛片免费观看3o分钟| 亚洲精品乱码久久久久久按摩| 国产 一区 欧美 日韩| 国产高清有码在线观看视频| 又粗又爽又猛毛片免费看| 亚洲天堂国产精品一区在线| 日本爱情动作片www.在线观看| 22中文网久久字幕| 尤物成人国产欧美一区二区三区| 亚洲av一区综合| 日本免费一区二区三区高清不卡| 直男gayav资源| 日本免费a在线| 天天躁夜夜躁狠狠久久av| 午夜久久久久精精品| 国产高清有码在线观看视频| 亚洲av中文av极速乱| 亚洲在线观看片| 国产成人a∨麻豆精品| 日本与韩国留学比较| 可以在线观看毛片的网站| 最后的刺客免费高清国语| av国产免费在线观看| 国产又色又爽无遮挡免| 国产精品熟女久久久久浪| 精品国产三级普通话版| 国产伦在线观看视频一区| 熟女人妻精品中文字幕| 丰满人妻一区二区三区视频av|