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

    一種廣域稀布雷達低旁瓣波束形成方法

    2023-06-02 06:33:14李鑫宇王鵬毅夏雙志
    計算機測量與控制 2023年5期
    關鍵詞:方向優(yōu)化

    李鑫宇,王鵬毅,夏雙志,宋 晨

    (1.中國電子科技集團公司 第54研究所,石家莊 050081;2.中國衛(wèi)星 海上測控部,江蘇省 江陰 214400)

    0 引言

    廣域稀布雷達作為一種新體制雷達,其陣元廣域部署,在部分單元雷達受到干擾或破壞時,不會對整個雷達系統(tǒng)造成太大影響,同時可以從不同方向對目標進行探測,充分利用空域信息,提高目標的探測效率,很好的滿足了“三反一抗”(反偵察、反干擾、反隱身、抗反輻射打擊)的要求。對于由N個相同孔徑的陣元構成的廣域稀布雷達系統(tǒng),通過信號級合成處理,其信噪比增益可達到單元孔徑的N3倍,也滿足了對于雷達系統(tǒng)大威力的需求[1-2]。對于廣域稀布雷達陣列,方向圖中會產生較高的旁瓣。高旁瓣的干擾,會在探測時容易產生虛假的目標。對于降低方向圖旁瓣電平的方法,國內外對此已經進行了一定的研究。付云起等利用遺傳算法和模擬退火對不等間距稀布陣進行了綜合設計,使得陣列旁瓣電平得到了改善[3]。文獻[4]提出一種基于自監(jiān)督差分算子的自適應遺傳算法的天線位置優(yōu)化方法,并利用自監(jiān)督微分算子進行突變,使得方向圖可以產生較低的旁瓣電平。鄭志東等提出了一種利用接收波束合成的算法來對消發(fā)射陣元的柵瓣,有效地抑制了高的旁瓣電平[5]。文獻[6]提出了一種多階量化幅相加權的方法,優(yōu)化后的方向圖在所有頻段和掃描范圍內均存在較低且穩(wěn)定的低旁瓣,增益損失小,原理簡單易于工程實現(xiàn)。文獻[7]提出了一種改進的粒子群算法,該算法將MIMO雷達的聯(lián)合收發(fā)波束作為適應度函數,分別對發(fā)射陣列及接收陣列進行了優(yōu)化,該算法收斂速度較快,且可以達到主瓣不展寬的條件下,降低方向圖的旁瓣電平。文獻[8]提出了一種幅度加權的波束形成算法,該算法可以在小范圍展寬主瓣波束寬度的條件下,實現(xiàn)大幅度降低旁瓣電平的效果。文獻[9]提出了一種綜合算法,對于由迭代傅里葉算法得到的稀疏陣列,以陣元間距大于半波長的陣元作為優(yōu)化對象并采取差分進化算法,經優(yōu)化后的陣列旁瓣電平得到較好的降低效果文獻[10]提出一種改進的遺傳算法,通過用多個矩陣組合對MIMO陣列進行表示,并應用基于混沌序列的方法擾動種群,以避免優(yōu)化過程中局部收斂的情況,該優(yōu)化算法很好地解決了MIMO陣列的排布問題,且實用性強。然而,現(xiàn)有文獻的方法主要應用于非稀布陣列或陣元間距略大于的半倍波長的陣列,而廣域稀布陣列的陣元間距大多為百倍甚至千倍的半倍波長。其中以遺傳算法為代表的搜索算法,通過優(yōu)化布陣而降低旁瓣的作用十分有限,難以滿足雷達系統(tǒng)對目標探測的要求。且陣元間距過大會導致搜索范圍擴大,在提高計算量的同時,也難以避免陷入局部最優(yōu)解[13],所以需要進一步采取抑制旁瓣的綜合方法。

    對此,本文以稀布線陣為例,提出一種脈組內捷變頻聯(lián)合波束形成的方法,通過將不同頻點下的信號方向圖相參的形式進行聯(lián)合,并結合遺傳算法進行布陣優(yōu)化。仿真結果表明,二者相結合后的算法可以進一步降低稀布線陣的旁瓣電平。最后將其拓展至稀布面陣,驗證了本文方法在稀布面陣中同樣具有良好的效果。

    1 廣域稀布陣列模型

    如圖1所示是一個由N個陣元組成的稀布陣列模型,每個陣元的坐標為(xn,yn,zn)。

    圖1 稀布陣列模型

    由陣列天線的遠場方向圖乘積原理,陣列天線方向圖等于陣元因子乘陣列因子,得:

    (1)

    其中:每個陣元方向圖fn(φ,θ)相同,不失一般性,假設陣元的方向圖范圍足夠大,滿足全向性,即fn(φ,θ)=1;陣列因子為

    (2)

    其中:ΔRn=xncosθcosφ+yncosθsinφ+znsinθ,ΔR0=xncosθ0cosφ0+yncosθ0sinφ0+znsinθ0。f為工作頻率,(θ0,φ0)為波束指向,θ為俯仰角,φ為方位角。天線的照射孔徑函數是等幅分布的,即不進行幅度加權,則幅度加權系數為An=1。稀布陣列方向圖函數可表示為[6,10]:

    (3)

    2 算法分析

    本章以線陣為模型,提出了一種脈組內捷變頻聯(lián)合波束形成方法,并將該方法與遺傳算法相結合。

    一個陣列孔徑為L的直線陣列方向圖函數表達式可由式(3)簡化為:

    (4)

    其中:dn表示第n個陣元到第一個陣元的距離,且n≤N。

    2.1 遺傳算法

    遺傳算法最早由美國的J.H.Holland[17]教授提出,是通過對自然界生物進化機制的探究而發(fā)展起來的一種全局優(yōu)化方法,誕生于20世紀60年代;經過一系列后續(xù)的研究發(fā)展,在80年代,遺傳算法經D.E.Goldberg大量的研究總結歸納后形成。

    遺傳算法的本質是一種高效的全局搜索算法,它可以在搜索過程中獲得并積累搜索空間的信息,并通過自適應的方式控制搜索過程從而求得最優(yōu)解。遺傳算法在每一代中,以個體為單位在問題域中通過適應度值以及交叉、變異、選擇等遺傳操作,產生新的個體。這個過程產生的新的個體比原個體對環(huán)境的適應度更大,從而最終提升了整個種群的適應度。目前遺傳算法已廣泛應用于各個領域,并在布陣優(yōu)化中產生了較好的效果。

    簡單來說,遺傳算法通過把種群當做問題的一組解,對其采用選擇。交叉和變異等一系列遺傳學操作從而產生新的種群,并反復迭代使種群進化到近似最優(yōu)解的狀態(tài)。遺傳算法是自然遺傳學與計算機科學互相融合而形成的算法,下面介紹有關的生物學中的基本概念和術語[6]。

    1)群體和個體:群體是生物進化過程中所有個體的集合,種群的規(guī)模表示可行解集。個體是種群的基本單元,表示一個可行解。

    2)染色體和基因:染色體包含生物所有的遺傳信息,通常是可行解的一段編碼?;蚴沁z傳信息的基本單位,是染色體的一部分,即可行解編碼的一部分。

    3)遺傳編碼:遺傳編碼通常以二進制編碼和十進制編碼為主。

    4)實數編碼:因為二進制編碼是離散的,存在較大誤差,有時也不利于反應問題的特定知識,故本文采用實數編碼的方式,其特點是精確度較高,便于和連續(xù)優(yōu)化算法結合,更實用與數值優(yōu)化的問題。

    5)適應度:適應度是用來評估個體適應生存環(huán)境的能力。適應度函數則用來評價個體的優(yōu)劣。

    6)遺傳操作:遺傳操作是“選擇”“交叉”“變異”的統(tǒng)稱。優(yōu)勢個體的保留和劣勢個體的淘汰稱作“選擇”,個體間基因的隨機交換稱作“交叉”,個體間基因的隨機變換稱作“變異”。

    7)選擇:選擇即是通過計算個體的適應度值,通過特定的方法或規(guī)則,從第t代群體中選出優(yōu)良的個體遺傳到下一代?!拜啽P賭”的選擇方法是一種最早提出的選擇方法,它是一種基于比例的選擇方法,對于每個個體,適應度值所占比例大的后代更容易被保留;反之,適應度值所占比例小的后代則更容易被淘汰。為了選擇交叉及變異的個體需要進行多輪選擇,選擇的作用主要在于保留更優(yōu)良的基因,保證全局更好的收斂性能。

    8)交叉:交叉是指將經過選擇操作后的群體,對每個個體以特定的交叉概率,交換他們的部分同一位置的基因,從而產生新的個體。交叉運算是產生新的子代的主要方式,也在一定程度上可以起到避免群體的局部收斂的作用。

    9)變異:變異是指對群體的每一個個體,依據特定的變異概率對染色體上的基因改為其他等位基因值,從而形成新的個體。變異和交叉同樣產生新的個體,和交叉操作共同保證遺傳算法向更優(yōu)良的群體方向進化。

    如下為遺傳算法的具體流程:

    1)種群初始化。在值域范圍內用隨機數生成N個個體,同時設置最大進化代數G。設置進化代數計數器g=0。

    2)適應度值個體評價。計算群體中的每個個體的適應度值。

    3)選擇:把選擇算子作用到群體。計算個體的適應度,并根據適應度將優(yōu)良的個體遺傳到下一代。

    4)交叉:把交叉算子作用到群體,對每個個體依概率交換其部分基因,產生新的個體。

    5)變異:把變異算子作用到群體,對每個個體依概率改變某些基因為值域內的其他基因。

    6)終止條件判斷。當g≤G,則g=g+1,并重復步驟2);當g>G時,終止循環(huán),并將最大適應度值的個體輸出。

    如圖2為遺傳算法的流程圖,這里種群的個體數量設置為Np,每個個體都是一個實數值參數向量,設置其長度為L,表示為:

    圖2 遺傳算法的運算流程

    xi,g(i=1,2,…,NP)

    (5)

    式中,g為遺傳代數,i為每一代種群中個體的序號,Np為每一代種群中個體的數量。

    首先對種群的每個個體進行初始的編碼,從而創(chuàng)建優(yōu)化的起點。由優(yōu)化模型可知,搜索空間的范圍為[0,L]。為了保證陣列的孔徑大小不發(fā)生變化,需要使陣列的首尾都存在一個陣元。同時約束相鄰兩陣元的最小距離間隔為常數dc,則約束條件設置為 :

    (6)

    這里把di分成xi和(i-1)dc兩部分,則有:

    (7)

    其中:

    x1≤x2≤x3≤…≤xN

    (8)

    由上式可得:

    x1≤x2≤x3≤…≤xN∈[0,L-(N-1)dc]

    (9)

    經過上述的操作,就可以把個體的基因距離間隔dm間接地轉化成xi,搜索空間的范圍從[0,L]減小為[0,L-(N-1)dc],并通過隨機數生成器在搜索空間內生成N個隨機數。

    此外,要想使(8)式成立,需要將種群初始的中間變量進行從小到大的排序,通過式(7)求出真實距離間隔種群d,并令d1=0,dN=L。至此完成了種群的初始化。

    接下來是優(yōu)化變量適應度函數。對于波束指向為θi的陣列,由方向圖的最大旁瓣電平(MSLL)定義可以得到計算公式為:

    其中:S為線陣方向圖的旁瓣區(qū)間,F(xiàn)dB(θ)為歸一化的線陣方向圖函數。

    S={θ|θmin≤θ≤θi-φ0∪θi+φ0≤θ≤θmax},其中2φ0是方向圖的主瓣的零功率點。從而可以定義如下適應度函數為:

    (11)

    式中,θi為空域內指定的i個波束指向。遺傳算法通過該適應度函數優(yōu)化陣元的位置,使得陣列方向圖在空域內的指定指向中,所有方向圖中最大的旁瓣電平最低。

    然后進行選擇操作。采用“輪盤賭”的思想,利用每個個體適應度大小來決定保留其子代的概率,若種群大小是Np,個體的適應度是fiti,則該個體被選擇的概率為:

    (12)

    最后進行交叉和變異操作。將選擇操作后的個體分為奇數及偶數兩組個體,并根據設置的交叉概率Pc依概率對奇數組個體及偶數組個體進行部分基因的交換。具體操作為:選擇要交叉的一對個體,若個體的長度為L,則在[1,L-1]的范圍內選擇整數k作為交叉位置,并根據交叉概率進行交叉操作,雙方交換交叉位置的基因,從而形成一對新的個體。

    變異操作是為了使種群保證其多樣性,同時可以防止遺傳算法達到局部最優(yōu)解而收斂。具體操作為:對于交叉操作完成的種群,對每個個體的所有基因以[0,1]區(qū)間內的隨機數r進行變異操作,若r

    通過以上操作完成一次遺傳操作,完成后需要保證陣列的孔徑不變,即首尾需要保證有陣元,同時計算每個個體的適應度計算,保留最優(yōu)個體到下一代,并進行下一代的遺傳操作,當上述步驟重復至一定次數或者滿足一定條件時,遺傳算法完成,此時的最優(yōu)個體為最終的優(yōu)化結果。

    2.2 脈組內捷變頻聯(lián)合波束形成方法

    本節(jié)提出一種脈組內捷變頻聯(lián)合波束形成方法,該方法將若干頻點下的方向圖以相參的形式進行聯(lián)合,可以形成具有低旁瓣的方向圖,如下為具體實現(xiàn)過程。

    首先確定一個可用頻段范圍[fmin,fmax],并確定一個的頻點數M,并通過(13)式確定M個頻點中第i個頻點的頻率fiM:

    (13)

    其中:fmin表示可用頻段范圍內最低頻點,fmax表示可用頻段范圍內最高頻點。

    對選用的M個頻點按照(4)式進行波束形成,通過(14)式進行聯(lián)合波束形成,并計算最大旁瓣電平:

    (14)

    其中:fi為第i個頻點,c為光速,S表示方向圖的旁瓣區(qū)間,若方向圖主瓣的零功率點為2φ0,則S={θ|θmin≤θ≤θi-φ0∪θi+φ0≤θ≤θmax}。

    同理,通過改變選用的頻點數M,可以得到多個頻點組合下的方向圖函數,并根據(15)式選取使得MSLLM最小的頻點數目k:

    (15)

    最后按照下式將選定的k個頻點進行聯(lián)合波束形成:

    (16)

    由此可以實現(xiàn)對稀布陣列的脈組內捷變頻聯(lián)合波束形成。

    2.3 結合遺傳算法的脈組內捷變頻聯(lián)合波束形成

    對于均勻線陣,脈組內捷變頻聯(lián)合波束形成方法對于降低方向圖的旁瓣電平具有一定效果,但對于稀布線陣而言效果仍不足。在確定適用的頻段及頻點數后,結合遺傳算法對稀布線陣進一步優(yōu)化,優(yōu)化后的陣元布局,具有更低的旁瓣電平。方向圖的最大旁瓣電平(MSLLii)可以表示為:

    (17)

    其中:θii為波束指向,fik為確定的k個頻點中第i個頻點,S為線陣方向圖的旁瓣區(qū)間,F(xiàn)dB(θ)為歸一化的方向圖函數。

    遺傳算法通過設置適應度函數對種群進行多次的迭代優(yōu)化,從而可以得到更優(yōu)良的子代。本文中適應度函數定義為使陣列在指定的若干個指向中形成的方向圖的最大旁瓣最小。因此遺傳算法的適應度函數可定義為:

    ii= 1,2,3…

    (18)

    其中:x為陣元位置,設定遺傳算法的約束條件為:保證優(yōu)化后陣列孔徑不變,且相鄰陣元間距不小于指定距離。由此適應度函數進行優(yōu)化后的陣列,可得到具有更低旁瓣的方向圖。

    3 仿真與分析

    本章先后對稀布線陣及面陣進行仿真,結果驗證了本文方法對于降低廣域稀布雷達陣列的旁瓣具有較好的效果。

    3.1 廣域稀布線陣

    設定如下直線陣列,陣元數N=20,工作頻率f=2.3 GHz,孔徑為1 900 m,波束指向為0°,當陣列均勻排布時,即陣元間距d=100 m,陣列方向圖如圖3所示。

    圖3 均勻直線陣列方向圖

    根據仿真結果可以看出對于陣元間距遠大于半倍波長的均勻線陣,陣元間互耦效應十分明顯,方向圖中產生了大量高旁瓣。

    對此均勻線陣,分別選取2~2.2 GHz,2 ~2.4 GHz,2 ~2.6 GHz,2 ~2.8 GHz,2 ~3 GHz五個頻段下的2~10個頻點,進行脈組內捷變頻聯(lián)合波束形成仿真,部分仿真結果如下。

    1)在2~2.4 GHz 頻段下,3、7、10個頻點聯(lián)合波束形成的方向圖如圖4所示。

    圖4 2~2.4 GHz頻段下3、7、10個頻點聯(lián)合波束形成方向圖

    由仿真可以看到,聯(lián)合波束形成后方向圖中不存在和主瓣相同增益的旁瓣,經過3個頻點的聯(lián)合波束形成,陣列方向圖中最大旁瓣降低為-8.95 dB,經過7個頻點的聯(lián)合波束形成,陣列方向圖中最大旁瓣降低為-13.71 dB,經過10個頻點的聯(lián)合波束形成,陣列方向圖中最大旁瓣降低為-13.89 dB。由此可以初步得到結論:進行聯(lián)合波束形成所采用的頻點數越多,旁瓣降低的效果越好。

    2)不同頻段下的不同頻點聯(lián)合后的最大旁瓣電平如圖5所示。

    圖5 不同頻段不同頻點聯(lián)合后的最大旁瓣電平

    由仿真結果可以看到:

    1)聯(lián)合波束形成能達到的最好的降低旁瓣的效果會隨著頻帶寬度的增大而增大。

    2)對于所有頻段,當使用的頻點數較少時,其旁瓣降低的效果相近,這是因為各頻點方向圖的旁瓣間距離較大,難以產生較好的聯(lián)合效果。

    3)在每個頻段下,聯(lián)合波束形成后旁瓣降低的效果都隨著頻點數的增大而增大,并逐漸達到一個上限,當進行聯(lián)合的頻點達到一定數量后,再繼續(xù)增加頻點數,降低旁瓣電平的效果將不再明顯。這是因為當頻點數較多時,各頻點方向圖的旁瓣間距較小,且對于高頻部分將變得較為密集,使聯(lián)合效果降低。因此,對于頻帶寬度的選擇而言,小的頻帶寬度效果較差,大的頻帶寬度雖然效果優(yōu)良,但是實施難度也會增加;對于頻點數的選擇而言,頻點數少旁瓣降低的效果差;頻點數過多,會增加較大的計算量,故應選取能夠達到上限的最小頻點數。

    經比較,選擇2~2.4 GHz頻段下的7個頻點與遺傳算法相結合并進行仿真,優(yōu)化前的方向圖如圖6所示。

    圖6 7頻點聯(lián)合波束形成后的方向圖

    可以看到距主瓣較遠處仍存在較高的旁瓣,在[-9°,9°]的范圍內最大旁瓣電平為-7.76 dB。對于遺傳算法,優(yōu)化的約束為保證陣列孔徑不變,即陣列首尾都要有陣元,優(yōu)化后陣元最小間隔不小于 1 m,在[-20°,20°]的掃描范圍內,保證在[-9°,9°]的范圍內以3°的間隔,使所有指向中最大的旁瓣電平最小。優(yōu)化后的陣列方向圖如圖7所示。

    圖7 結合遺傳算法的脈組內捷變頻聯(lián)合波束形成的方向圖

    圖 8 均勻面陣方向圖

    圖9 不同頻段不同頻點聯(lián)合后的最大旁瓣電平(方位向)

    由仿真結果可以看到,兩種方法優(yōu)化后,在[-9°,9°]的空域范圍內,最大旁瓣電平從-7.76 dB進一步降低為-11.85 dB,降低了4.09 dB。表 1為前后兩種方法的進行的對比結果。

    表1 兩種方法優(yōu)化后最大旁瓣電平比較(線陣)

    3.2 廣域稀布面陣

    由3.1中仿真結果可以看到本文方法對于線陣具有較好的效果,由于面陣原理與線陣原理相似,在此將本文方法直接應用于稀布面陣中。

    設定如下稀布面陣,陣列孔徑90 m×90 m,工作頻率f=2.3 GHz,陣元數N=100,波束指向為(0°,0°)如圖8為10×10均勻布置時的方向圖。

    可以看到對于陣元間隔遠大于半倍波長的均勻面陣,方向圖中同樣產生大量高旁瓣。

    對于均勻面陣,同樣選取2~2.2 GHz,2~2.4 GHz,2~2.6 GHz,2~2.8 GHz,2~3 GHz五個頻段下的2~10個頻點,進行聯(lián)合波束形成仿真,仿真結果如下(對均勻面陣,方位向和俯仰向的結果相同,故只展示方位向的結果)。

    由仿真可以看到結論和3.1節(jié)中均勻線陣基本一致。

    本節(jié)選擇2~2.6 GHz頻段下的6個頻點,結合遺傳算法對平面陣列進行優(yōu)化,其中遺傳算法優(yōu)化的約束為保證陣列孔徑不變,即陣列4個角都要有陣元,優(yōu)化后陣元最小間隔不小于,在方位向及俯仰向的[-20°,20°]的掃描范圍內,保證在[-9°,9°]的范圍內以3°的間隔,使所有指向中最大的旁瓣電平最小。優(yōu)化后的陣列方向圖的仿真結果如圖10所示。

    圖10 結合遺傳算法的脈組內捷變頻聯(lián)合波束形成的方向圖

    由仿真結果可以看到優(yōu)化后方位向的最大旁瓣降低為-14.18 dB,俯仰向的最大旁瓣電平降低為-14.73 dB。

    由稀布線陣及面陣的仿真結果可以得出結論:通過將脈組內捷變頻聯(lián)合波束形成方法和遺傳算法相結合,可以有效的降低廣域稀布陣列方向圖旁瓣電平。

    4 結束語

    針對廣域稀布雷達陣列方向圖會產生大量高旁瓣的問題,本文提出一種結合遺傳算法的廣域稀布陣列脈組內捷變頻聯(lián)合波束形成方法。首先在一定帶寬,一定頻點范圍內,選擇合適的頻點組合,并結合遺傳算法,通過約束陣列孔徑,最小陣元間距,并通過多個指向對適應度函數進行設置進行優(yōu)化。仿真結果表明,本文方法可以保證廣域稀布雷達陣列在滿足自身的強抗干擾能力的同時,也能獲得低旁瓣的方向圖,對廣域稀布雷達低旁瓣的波束形成方法具有一定參考意義。對于如何選擇合適的頻點數及更好的聯(lián)合方式將是下一步的研究重點。

    猜你喜歡
    方向優(yōu)化
    超限高層建筑結構設計與優(yōu)化思考
    房地產導刊(2022年5期)2022-06-01 06:20:14
    2022年組稿方向
    計算機應用(2022年2期)2022-03-01 12:33:42
    2022年組稿方向
    計算機應用(2022年1期)2022-02-26 06:57:42
    民用建筑防煙排煙設計優(yōu)化探討
    關于優(yōu)化消防安全告知承諾的一些思考
    一道優(yōu)化題的幾何解法
    2021年組稿方向
    計算機應用(2021年4期)2021-04-20 14:06:36
    2021年組稿方向
    計算機應用(2021年3期)2021-03-18 13:44:48
    由“形”啟“數”優(yōu)化運算——以2021年解析幾何高考題為例
    2021年組稿方向
    計算機應用(2021年1期)2021-01-21 03:22:38
    亚洲av电影在线进入| 亚洲av五月六月丁香网| av片东京热男人的天堂| 19禁男女啪啪无遮挡网站| 草草在线视频免费看| 亚洲精品在线观看二区| www日本在线高清视频| 观看免费一级毛片| 精品午夜福利视频在线观看一区| 午夜精品在线福利| 国产精品久久久久久人妻精品电影| 黄色视频,在线免费观看| 狠狠狠狠99中文字幕| 亚洲狠狠婷婷综合久久图片| 国产成人系列免费观看| 国产激情偷乱视频一区二区| 日本在线视频免费播放| 国产单亲对白刺激| 法律面前人人平等表现在哪些方面| 美女高潮的动态| 狠狠狠狠99中文字幕| 99热只有精品国产| 欧美日本亚洲视频在线播放| 99视频精品全部免费 在线 | 午夜免费成人在线视频| 日韩精品中文字幕看吧| 国产精品亚洲美女久久久| 精品99又大又爽又粗少妇毛片 | 精品国产三级普通话版| 午夜两性在线视频| 国产亚洲精品av在线| 特大巨黑吊av在线直播| 给我免费播放毛片高清在线观看| 搡老妇女老女人老熟妇| 国产伦精品一区二区三区视频9 | 人人妻人人看人人澡| 一级毛片高清免费大全| 国产精品野战在线观看| 在线免费观看不下载黄p国产 | 啪啪无遮挡十八禁网站| 中亚洲国语对白在线视频| 欧美乱妇无乱码| 精品一区二区三区四区五区乱码| 黄片大片在线免费观看| 国产精品 国内视频| 欧美成人一区二区免费高清观看 | 国产精品一区二区精品视频观看| 1000部很黄的大片| 亚洲成人久久性| 三级男女做爰猛烈吃奶摸视频| 一本精品99久久精品77| 黑人欧美特级aaaaaa片| а√天堂www在线а√下载| 亚洲美女视频黄频| 男女做爰动态图高潮gif福利片| 精品乱码久久久久久99久播| 久久香蕉国产精品| 九九热线精品视视频播放| 国产主播在线观看一区二区| av天堂在线播放| 国产极品精品免费视频能看的| 精品久久久久久久人妻蜜臀av| 丁香六月欧美| 人妻丰满熟妇av一区二区三区| 小蜜桃在线观看免费完整版高清| 精品久久久久久久人妻蜜臀av| 久久久久久久久中文| 久久热在线av| 99热这里只有是精品50| 国产精品野战在线观看| 制服人妻中文乱码| svipshipincom国产片| av黄色大香蕉| 午夜精品一区二区三区免费看| 看黄色毛片网站| 日本熟妇午夜| 99热精品在线国产| 午夜精品在线福利| 1024手机看黄色片| 最近在线观看免费完整版| 国产美女午夜福利| 精品久久久久久久末码| 听说在线观看完整版免费高清| 国产欧美日韩精品亚洲av| 国产精品,欧美在线| 国产不卡一卡二| 亚洲成人精品中文字幕电影| 好男人电影高清在线观看| 岛国在线观看网站| 国产精品亚洲一级av第二区| 成人国产一区最新在线观看| 国产99白浆流出| 一级毛片高清免费大全| 国产在线精品亚洲第一网站| 亚洲国产精品久久男人天堂| 午夜精品在线福利| 免费在线观看亚洲国产| 中文字幕精品亚洲无线码一区| 久久精品国产99精品国产亚洲性色| 国内久久婷婷六月综合欲色啪| 国产真实乱freesex| 成人午夜高清在线视频| 亚洲无线在线观看| 亚洲欧美一区二区三区黑人| 美女cb高潮喷水在线观看 | 色综合婷婷激情| 91老司机精品| 午夜日韩欧美国产| 曰老女人黄片| 日韩人妻高清精品专区| 亚洲国产看品久久| 色综合站精品国产| 亚洲国产欧美人成| 日本 av在线| 亚洲精品在线观看二区| 最好的美女福利视频网| 观看免费一级毛片| 久久久久国产精品人妻aⅴ院| 一级黄色大片毛片| 国产野战对白在线观看| 欧美黄色片欧美黄色片| 岛国视频午夜一区免费看| 国产高潮美女av| 国产爱豆传媒在线观看| av在线蜜桃| av中文乱码字幕在线| 亚洲成人中文字幕在线播放| 久久天堂一区二区三区四区| 欧美性猛交╳xxx乱大交人| 精品久久久久久久末码| 两性夫妻黄色片| 精品一区二区三区视频在线 | 精品国产乱码久久久久久男人| 免费观看人在逋| 夜夜爽天天搞| 18禁黄网站禁片午夜丰满| 在线看三级毛片| 欧美3d第一页| 亚洲人成电影免费在线| 亚洲人成电影免费在线| 黄色片一级片一级黄色片| 欧美成人免费av一区二区三区| 免费观看精品视频网站| 免费一级毛片在线播放高清视频| 淫妇啪啪啪对白视频| 99精品在免费线老司机午夜| 国产精品电影一区二区三区| 亚洲一区二区三区不卡视频| 亚洲人成网站在线播放欧美日韩| 熟女电影av网| 国产一级毛片七仙女欲春2| 成人av一区二区三区在线看| 亚洲专区中文字幕在线| 国产三级在线视频| 国产主播在线观看一区二区| 国产亚洲欧美在线一区二区| 久久久久久大精品| 国产高清视频在线观看网站| 长腿黑丝高跟| 国产精品九九99| 国产精品九九99| 亚洲无线在线观看| 国产精品自产拍在线观看55亚洲| 天天添夜夜摸| 91久久精品国产一区二区成人 | 日本黄大片高清| 巨乳人妻的诱惑在线观看| 日韩三级视频一区二区三区| 久久久久久久久久黄片| 久久久水蜜桃国产精品网| 国内毛片毛片毛片毛片毛片| www.www免费av| 又爽又黄无遮挡网站| 国产综合懂色| 岛国在线观看网站| 999久久久精品免费观看国产| 日本五十路高清| 国产成年人精品一区二区| 51午夜福利影视在线观看| 精品熟女少妇八av免费久了| 在线播放国产精品三级| 国产伦精品一区二区三区四那| 麻豆一二三区av精品| 午夜精品一区二区三区免费看| 麻豆一二三区av精品| 亚洲人成网站在线播放欧美日韩| 婷婷精品国产亚洲av在线| 成年女人永久免费观看视频| 一本久久中文字幕| 日本三级黄在线观看| 国产伦一二天堂av在线观看| 国产精品亚洲一级av第二区| 亚洲自偷自拍图片 自拍| 搞女人的毛片| 一个人看的www免费观看视频| 亚洲电影在线观看av| 久久人人精品亚洲av| 中文字幕人妻丝袜一区二区| 亚洲国产精品成人综合色| 少妇的丰满在线观看| 中文字幕高清在线视频| 日本熟妇午夜| 12—13女人毛片做爰片一| 国产欧美日韩一区二区三| 免费人成视频x8x8入口观看| 午夜精品在线福利| 欧美另类亚洲清纯唯美| 天堂av国产一区二区熟女人妻| 1000部很黄的大片| avwww免费| 精品免费久久久久久久清纯| 欧美黑人巨大hd| 欧美黑人巨大hd| 99riav亚洲国产免费| 成年女人毛片免费观看观看9| 亚洲专区字幕在线| 日本免费一区二区三区高清不卡| 老汉色av国产亚洲站长工具| 少妇裸体淫交视频免费看高清| 69av精品久久久久久| 男人和女人高潮做爰伦理| 国产精品亚洲av一区麻豆| 国产又黄又爽又无遮挡在线| 美女免费视频网站| 99精品欧美一区二区三区四区| 听说在线观看完整版免费高清| 久久久久九九精品影院| 亚洲欧美精品综合久久99| 国产综合懂色| 中文亚洲av片在线观看爽| 男女床上黄色一级片免费看| 午夜精品在线福利| 国产成人精品久久二区二区免费| 一二三四在线观看免费中文在| 动漫黄色视频在线观看| 日韩精品青青久久久久久| 一a级毛片在线观看| 五月玫瑰六月丁香| 最近视频中文字幕2019在线8| 天天躁日日操中文字幕| 成人性生交大片免费视频hd| 色视频www国产| av黄色大香蕉| 国产亚洲精品一区二区www| 女同久久另类99精品国产91| 日韩人妻高清精品专区| 精品人妻1区二区| 他把我摸到了高潮在线观看| 色综合站精品国产| 男女视频在线观看网站免费| 久9热在线精品视频| 丰满人妻熟妇乱又伦精品不卡| 久久热在线av| 国产精品99久久99久久久不卡| 少妇人妻一区二区三区视频| 在线观看美女被高潮喷水网站 | 成人性生交大片免费视频hd| 午夜福利成人在线免费观看| 国产精品99久久久久久久久| 午夜福利在线在线| 丰满人妻一区二区三区视频av | 亚洲自拍偷在线| 18禁美女被吸乳视频| 成熟少妇高潮喷水视频| 国产毛片a区久久久久| 国产高清激情床上av| 精品久久久久久久毛片微露脸| 嫁个100分男人电影在线观看| 中文字幕最新亚洲高清| 女同久久另类99精品国产91| 亚洲美女视频黄频| 中国美女看黄片| 国产高清videossex| 国产精品亚洲美女久久久| 欧美三级亚洲精品| 神马国产精品三级电影在线观看| 色综合亚洲欧美另类图片| 欧美乱妇无乱码| 熟女人妻精品中文字幕| 少妇的逼水好多| 神马国产精品三级电影在线观看| 床上黄色一级片| 亚洲成人久久爱视频| 午夜日韩欧美国产| 亚洲专区中文字幕在线| 大型黄色视频在线免费观看| 国产三级中文精品| 国产真人三级小视频在线观看| 国产真人三级小视频在线观看| 最近视频中文字幕2019在线8| 很黄的视频免费| 在线国产一区二区在线| www.自偷自拍.com| 久久精品综合一区二区三区| 女警被强在线播放| 一级毛片女人18水好多| 999久久久精品免费观看国产| 久久久久久久久久黄片| 国产淫片久久久久久久久 | 日韩精品中文字幕看吧| 精华霜和精华液先用哪个| 狠狠狠狠99中文字幕| 午夜久久久久精精品| 亚洲国产高清在线一区二区三| 中文字幕熟女人妻在线| 非洲黑人性xxxx精品又粗又长| 欧美一级a爱片免费观看看| 欧美激情久久久久久爽电影| 久久国产乱子伦精品免费另类| 18禁美女被吸乳视频| 99热这里只有精品一区 | 久久精品夜夜夜夜夜久久蜜豆| 又大又爽又粗| 又爽又黄无遮挡网站| netflix在线观看网站| 国产精品综合久久久久久久免费| 欧美黄色片欧美黄色片| 中国美女看黄片| 国产成人啪精品午夜网站| 欧美精品啪啪一区二区三区| 村上凉子中文字幕在线| 噜噜噜噜噜久久久久久91| 在线观看免费午夜福利视频| 日日摸夜夜添夜夜添小说| 午夜福利18| 欧美zozozo另类| 久久久久精品国产欧美久久久| 成人国产综合亚洲| 精品人妻1区二区| 久久久久久国产a免费观看| 成年人黄色毛片网站| 天堂影院成人在线观看| 午夜日韩欧美国产| 国内少妇人妻偷人精品xxx网站 | e午夜精品久久久久久久| 高潮久久久久久久久久久不卡| 一级黄色大片毛片| 18禁黄网站禁片免费观看直播| 成在线人永久免费视频| 国产成人av教育| 免费av毛片视频| 黄频高清免费视频| av片东京热男人的天堂| 日日夜夜操网爽| 亚洲国产高清在线一区二区三| 国产主播在线观看一区二区| 狠狠狠狠99中文字幕| 国产成人欧美在线观看| 好看av亚洲va欧美ⅴa在| 日本黄大片高清| 国产高清视频在线观看网站| 国产亚洲av嫩草精品影院| 亚洲在线自拍视频| av天堂中文字幕网| 亚洲精品在线观看二区| 免费看美女性在线毛片视频| 三级毛片av免费| 国产免费男女视频| 啦啦啦免费观看视频1| 成在线人永久免费视频| xxxwww97欧美| 欧美大码av| 熟女电影av网| 亚洲国产欧美人成| 欧美丝袜亚洲另类 | 白带黄色成豆腐渣| 色综合亚洲欧美另类图片| 欧美日韩一级在线毛片| 在线国产一区二区在线| 一区福利在线观看| 久久性视频一级片| 久久久久久久久中文| 最新中文字幕久久久久 | 狂野欧美白嫩少妇大欣赏| 最新美女视频免费是黄的| 女同久久另类99精品国产91| 哪里可以看免费的av片| 精品一区二区三区四区五区乱码| 又紧又爽又黄一区二区| 精品乱码久久久久久99久播| 性欧美人与动物交配| 岛国在线观看网站| 一区二区三区国产精品乱码| 国产蜜桃级精品一区二区三区| 国产v大片淫在线免费观看| 国产爱豆传媒在线观看| 久久性视频一级片| 人人妻人人看人人澡| 一进一出好大好爽视频| 久久久国产成人精品二区| 国产极品精品免费视频能看的| 日韩欧美免费精品| 99热这里只有是精品50| 美女午夜性视频免费| 美女cb高潮喷水在线观看 | 巨乳人妻的诱惑在线观看| 可以在线观看的亚洲视频| 国产精品1区2区在线观看.| 嫩草影院精品99| 国产综合懂色| 搡老岳熟女国产| 熟女人妻精品中文字幕| 黄色日韩在线| 舔av片在线| 亚洲在线观看片| 男女下面进入的视频免费午夜| 亚洲人与动物交配视频| 亚洲成人免费电影在线观看| avwww免费| 99久久国产精品久久久| 在线观看免费视频日本深夜| 国产精品av久久久久免费| 亚洲av成人精品一区久久| 国产一级毛片七仙女欲春2| 老汉色av国产亚洲站长工具| 日本撒尿小便嘘嘘汇集6| 国产毛片a区久久久久| 亚洲中文av在线| 成年人黄色毛片网站| 美女被艹到高潮喷水动态| 国产成人系列免费观看| 午夜亚洲福利在线播放| 伦理电影免费视频| 亚洲国产精品sss在线观看| 欧美色视频一区免费| 国产精品98久久久久久宅男小说| 天天添夜夜摸| 久久久久精品国产欧美久久久| 最好的美女福利视频网| 欧美乱妇无乱码| 婷婷亚洲欧美| 听说在线观看完整版免费高清| 国产又黄又爽又无遮挡在线| 亚洲成av人片在线播放无| 欧美日本亚洲视频在线播放| 真人一进一出gif抽搐免费| 婷婷精品国产亚洲av| 欧美性猛交╳xxx乱大交人| 国产成人啪精品午夜网站| 伊人久久大香线蕉亚洲五| 最近最新免费中文字幕在线| 看黄色毛片网站| 一进一出好大好爽视频| 美女高潮的动态| 日本免费一区二区三区高清不卡| 国模一区二区三区四区视频 | 九色成人免费人妻av| 久久伊人香网站| 巨乳人妻的诱惑在线观看| 久久精品国产亚洲网站| 一区二区三区高清视频在线| 日日摸夜夜添夜夜爱| kizo精华| 爱豆传媒免费全集在线观看| 中文字幕免费在线视频6| 久久久久久九九精品二区国产| 一边摸一边抽搐一进一小说| 欧美日韩国产亚洲二区| 亚洲18禁久久av| 亚洲国产精品久久男人天堂| 黄片无遮挡物在线观看| 国产精品美女特级片免费视频播放器| 亚洲色图av天堂| 综合色av麻豆| 久久久精品大字幕| 久久精品国产自在天天线| 建设人人有责人人尽责人人享有的 | 精品一区二区免费观看| 毛片女人毛片| 九草在线视频观看| 日日干狠狠操夜夜爽| 亚洲天堂国产精品一区在线| av国产免费在线观看| 蜜桃亚洲精品一区二区三区| 蜜臀久久99精品久久宅男| 久久精品夜色国产| 国产精品久久久久久久电影| 欧美xxxx黑人xx丫x性爽| 岛国毛片在线播放| 亚洲av一区综合| 91精品一卡2卡3卡4卡| 男女国产视频网站| 夜夜爽夜夜爽视频| 亚洲熟妇中文字幕五十中出| 高清在线视频一区二区三区 | av视频在线观看入口| 欧美激情国产日韩精品一区| 一个人免费在线观看电影| 国产伦理片在线播放av一区| 久久久久久国产a免费观看| 天堂影院成人在线观看| 欧美日本亚洲视频在线播放| 人妻少妇偷人精品九色| 亚洲中文字幕一区二区三区有码在线看| 亚洲aⅴ乱码一区二区在线播放| 日日摸夜夜添夜夜添av毛片| 舔av片在线| 亚洲性久久影院| 男女国产视频网站| 免费一级毛片在线播放高清视频| 国产成人91sexporn| 欧美极品一区二区三区四区| 汤姆久久久久久久影院中文字幕 | 成年女人永久免费观看视频| 老师上课跳d突然被开到最大视频| 日韩精品有码人妻一区| 18禁动态无遮挡网站| 成人三级黄色视频| 人妻系列 视频| 亚洲国产精品sss在线观看| 久热久热在线精品观看| 亚洲精品,欧美精品| 日韩 亚洲 欧美在线| 日本免费在线观看一区| 成人无遮挡网站| 成人高潮视频无遮挡免费网站| 国产黄色小视频在线观看| 国产成人freesex在线| 美女xxoo啪啪120秒动态图| 天堂网av新在线| 97超视频在线观看视频| 成人亚洲精品av一区二区| 久久国产乱子免费精品| 秋霞在线观看毛片| 日韩欧美精品v在线| 特大巨黑吊av在线直播| 熟女电影av网| 久久久久久久久久黄片| 可以在线观看毛片的网站| 欧美一区二区国产精品久久精品| 免费观看精品视频网站| 久久久精品94久久精品| 中文资源天堂在线| 在线观看美女被高潮喷水网站| av在线亚洲专区| 视频中文字幕在线观看| 大又大粗又爽又黄少妇毛片口| 久久亚洲精品不卡| 蜜桃亚洲精品一区二区三区| 亚洲欧美精品专区久久| 麻豆乱淫一区二区| av在线观看视频网站免费| 午夜激情欧美在线| 久久久亚洲精品成人影院| 亚洲婷婷狠狠爱综合网| 女人久久www免费人成看片 | 看黄色毛片网站| 男人舔女人下体高潮全视频| 亚洲国产成人一精品久久久| 欧美日韩在线观看h| 婷婷色av中文字幕| 亚洲成人中文字幕在线播放| 亚洲精品乱码久久久v下载方式| 99久久精品热视频| 级片在线观看| 国产精品国产三级国产av玫瑰| 乱码一卡2卡4卡精品| 国产亚洲91精品色在线| 波多野结衣巨乳人妻| 久久综合国产亚洲精品| 日韩视频在线欧美| 亚洲国产精品专区欧美| 最近最新中文字幕大全电影3| 亚洲av免费在线观看| 欧美极品一区二区三区四区| 男人舔奶头视频| 国产精品伦人一区二区| 亚洲欧美成人综合另类久久久 | 免费播放大片免费观看视频在线观看 | 欧美高清性xxxxhd video| 爱豆传媒免费全集在线观看| 日本三级黄在线观看| 日韩制服骚丝袜av| 亚洲图色成人| 国产成人免费观看mmmm| 真实男女啪啪啪动态图| 国产一区亚洲一区在线观看| 久久久久久久亚洲中文字幕| 国产精品一及| 精品人妻一区二区三区麻豆| 久久久久久久久大av| 国产色爽女视频免费观看| 亚洲国产欧美在线一区| 真实男女啪啪啪动态图| 亚洲18禁久久av| 日韩欧美在线乱码| 搡老妇女老女人老熟妇| 女人久久www免费人成看片 | 国语自产精品视频在线第100页| 成人无遮挡网站| 亚洲欧美精品自产自拍| 人妻少妇偷人精品九色| 蜜臀久久99精品久久宅男| 国产精品伦人一区二区| 搡女人真爽免费视频火全软件| 老女人水多毛片| 国产精品一区二区三区四区久久| 国产精品野战在线观看| 成人漫画全彩无遮挡| 久久精品国产亚洲av涩爱| 少妇猛男粗大的猛烈进出视频 | 欧美色视频一区免费| 欧美精品一区二区大全| 在线免费十八禁| 国产高潮美女av| 精品人妻熟女av久视频| 精品一区二区三区人妻视频| 久久久色成人| 成人毛片60女人毛片免费| 国内精品宾馆在线| 精品一区二区三区人妻视频| 国产亚洲精品久久久com| 国产成人精品久久久久久| 精品一区二区三区视频在线| 成人毛片60女人毛片免费|