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

    采用FFT-PSO策略的多約束稀疏陣天線方向圖綜合*

    2017-10-23 03:05:32潘紅光1薛紀(jì)康1高海南江能前
    電訊技術(shù) 2017年10期
    關(guān)鍵詞:波束寬度副瓣電平

    潘紅光1,薛紀(jì)康1,高海南**2,江能前

    (1.西安科技大學(xué) 電氣與控制工程學(xué)院,西安710054;2.西安衛(wèi)星測控中心,陜西 渭南714000;3.西安交通大學(xué) 經(jīng)濟(jì)金融學(xué)院,西安710049)

    采用FFT-PSO策略的多約束稀疏陣天線方向圖綜合*

    潘紅光1,薛紀(jì)康1,高海南**2,江能前3

    (1.西安科技大學(xué) 電氣與控制工程學(xué)院,西安710054;2.西安衛(wèi)星測控中心,陜西 渭南714000;3.西安交通大學(xué) 經(jīng)濟(jì)金融學(xué)院,西安710049)

    采用稀疏布陣的相控陣測控天線,陣元等幅激勵時的陣列峰值副瓣電平無法滿足指標(biāo)要求。在考慮天線增益、半功率波束寬度等指標(biāo)約束時設(shè)計了合適的適應(yīng)度函數(shù),并基于粒子群算法(PSO)對陣元激勵幅度進(jìn)行優(yōu)化以降低天線副瓣電平。同時,分析了等間距柵格平面陣列天線方向圖與離散傅里葉變換的關(guān)系,采用快速傅里葉變換(FFT) 降低了陣列方向圖計算復(fù)雜度。仿真結(jié)果表明,該方法可有效降低峰值副瓣電平和計算復(fù)雜度。

    航天測控;相控陣天線;稀疏平面陣列;粒子群優(yōu)化;方向圖綜合

    1 引 言

    隨著航天測控技術(shù)的發(fā)展和測控任務(wù)的增加,相控陣天線在測控設(shè)備天線的設(shè)計中應(yīng)用越來越廣泛。相控陣天線具備波束賦形和掃描靈活的優(yōu)勢,可用來提高陣列系統(tǒng)的整體性能,增加對輻射方向圖的控制能力[1-3]。實際中,需要天線陣列的優(yōu)化布陣技術(shù)以滿足其設(shè)計中成本、結(jié)構(gòu)和設(shè)備功能等方面的需求。一般來說,相控陣天線陣列設(shè)計分為稀疏和稀布兩種布陣方式[4]。稀疏陣列是將一定數(shù)目的天線單元從均勻間隔陣列中稀疏而形成的,其陣元間距是原均勻陣列的陣元間距的整數(shù)倍;而稀布陣列的天線單元在天線孔徑上隨機(jī)分布,其陣元間距是相互不可整除的[4-5]。相較而言,由于稀疏陣列的結(jié)構(gòu)設(shè)計和波束控制簡便易行,因此在工程實際中應(yīng)用更為常見。

    對稀疏陣列天線而言,工程中采用陣元等幅度激勵方法時,稀疏陣列天線方向圖的峰值副瓣電平(Peak SideLobe Level,PSLL) 很難滿足指標(biāo)要求,并會帶來一系列問題[6],如:較高的PSLL 會增大副瓣誤跟蹤的幾率,進(jìn)而增加任務(wù)風(fēng)險。因此,需要對稀疏平面陣列天線方向圖進(jìn)行低副瓣綜合以降低PSLL[7]。

    對測控(Telemetry,Track and Command,TT&C)天線而言,接收天線增益反映接收靈敏,半功率波束寬度BW(又稱主瓣寬度)對目標(biāo)測向精度有很大的影響。由于基于陣元激勵幅值的方向圖綜合時降低陣列天線峰值副瓣電平的同時,一般也會降低天線增益和展寬半功率波束寬度,導(dǎo)致接收信號電平門限的降低和測向數(shù)據(jù)誤差的增大,因此需要在設(shè)計適應(yīng)度函數(shù)時考慮考慮以上兩個指標(biāo)約束。

    稀疏平面陣列方向圖綜合問題中的目標(biāo)函數(shù)和約束條件往往是多參數(shù)、非線性、不可微的,因此傳統(tǒng)的解析方法和梯度尋優(yōu)數(shù)值方法均無法有效求解[9]。而遺傳算法(Genetic Algorithm,GA)、粒子群算法(Particle Swarm Optimization,PSO)、自適應(yīng)差分算法和模擬退火算法等智能算法[9-14]在求解多變量非線性優(yōu)化問題上具有很大的優(yōu)勢,算法簡單、易于實現(xiàn),因此在陣列天線方向圖綜合中得到了廣泛的應(yīng)用。譬如:文獻(xiàn)[14]中所采用的自適應(yīng)差分算法在處理分布式孔徑雷達(dá)的方向圖問題時,僅通過較少的仿真次數(shù)即可穩(wěn)定的獲得更低的旁瓣電平。

    粒子群算法源自對鳥群捕食行為的研究,由Kennedy 和Eberhart[15]提出。基于標(biāo)準(zhǔn)PSO算法及改進(jìn)的PSO算法優(yōu)化稀布陣列天線間距和稀疏陣列天線的稀疏位置和陣元激勵幅度、相位,可得到滿意的PSLL,但構(gòu)造適應(yīng)度函數(shù)時很少考慮工程中對天線增益和半功率波束寬度等指標(biāo)約束要求和柵瓣存在的影響。對于PSO等智能優(yōu)化算法,采用級數(shù)求和求解天線方向圖的方法也使得計算復(fù)雜度呈指數(shù)級增長[16]。當(dāng)陣元間距等于半波長的特殊情況下,可采用快速傅里葉變換(Fast Fourier Transform Algorithm,F(xiàn)FT)計算直接得到天線方向圖,但該條件需放寬至任意陣元間距情況方能滿足測控天線方向圖計算需求。

    本文建立了稀疏平面陣列天線數(shù)學(xué)模型,在此基礎(chǔ)上構(gòu)造了滿足TT&C天線性能指標(biāo)要求的適應(yīng)度函數(shù),并采用粒子群算法實現(xiàn)了稀疏陣列天線低副瓣方向圖綜合。另外,在分析等間距柵格平面陣列天線方向圖與FFT關(guān)系的基礎(chǔ)上,本文利用離散傅里葉變換(Discrete Fourier Transformation,DFT)的周期性拓展方向圖,得到了任意陣元間距情況陣列方向圖,降低了計算復(fù)雜度。仿真驗證了方法的有效性。

    2 稀疏平面陣列天線方向圖

    2.1稀疏陣列天線方向圖數(shù)學(xué)模型

    對于圖1所示的稀疏平面陣列天線,陣元只分布于x-y平面內(nèi)的M×N個等間距柵格上,dx為x方向的柵格間距,dy為y方向的柵格間距。定義矩陣A表示陣元稀疏信息,Am,n=1表示(m,n)坐標(biāo)柵格上有陣元,Am,n=0表示該位置陣元被稀疏掉,m∈[0,M-1],n∈[0,N-1]。平面陣列方向圖為

    (1)

    其中:xm=mdx;yn=ndy;θ為波束指向與陣列平面法線的夾角;φ為波束指向與陣列平面x正方向的夾角;Cm,n為每個陣元的口徑分布展開的復(fù)電壓,Cm,n=cm,nejΘm;n,cm,n為陣元激勵幅度,cm,n∈[0,1],Θm,n為歸一化陣元激勵相位,Θm,n∈[0,2π],均受陣元中的移相器控制[17]。

    圖1 稀疏平面陣列天線模型Fig.1 The model of sparse planar array antenna

    若僅考慮幅度激勵,則對?m∈[0,M-1]和?n∈[0,N-1],有Θm,n=0,式(1)可變?yōu)?/p>

    (2)

    由天線方向圖的乘積原理可知,相控陣天線的天線方向圖為陣元方向圖e(θ,φ)與陣列因子f(θ,φ)的乘積,即

    F(θ,φ)=e(θ,φ)f(θ,φ),

    (3)

    (4)

    2.2陣列天線方向圖的FFT計算方法

    陣元方向圖e(θ,φ)固定不變,因此計算陣列天線方向圖需首先計算陣列因子f(θ,φ)。 直接按照式(2)計算方向圖需要(MN)2次復(fù)數(shù)乘法和MN(MN-1)次復(fù)數(shù)加法,當(dāng)陣元個數(shù)較多時,計算量很大。通過變換將式(2)寫為

    (5)

    圖2 FFT 變換后的u-v空間示意圖Fig. 2 The u-v space diagram after FFT transformation

    對u和v進(jìn)行變換處理有

    (6)

    uL≤u≤uR,vL≤v≤vR。

    (7)

    通過FFT方法可將陣列因子方向圖的計算復(fù)雜度由(MN)2降低到(MN)lb(MN)數(shù)量級,提高了優(yōu)化算法的效率[16]。

    3 基于粒子群優(yōu)化的多約束方向圖綜合

    3.1粒子群優(yōu)化算法簡介

    粒子群算法首先在給定的解空間中隨機(jī)初始化粒子群,待優(yōu)化問題的變量數(shù)決定了解空間的維數(shù)。每個粒子均有初始位置、初始速度和優(yōu)化目標(biāo)構(gòu)造的適應(yīng)度函數(shù),在解空間中按照一定的速度迭代尋優(yōu)。每次迭代中,粒子通過跟蹤兩個極值來更新自己在解空間中的位置與速度,一個是單個粒子在迭代過程中找到的最優(yōu)個體極值,另一個為種群所有粒子在迭代過程中所找到的全局最優(yōu)極值。

    粒子群算法的速度、位置的迭代方程為

    (8a)

    (8b)

    3.2適應(yīng)度函數(shù)設(shè)計

    適應(yīng)度函數(shù)的設(shè)計是粒子群算法的關(guān)鍵,合適的適應(yīng)度函數(shù)能夠提高算法的優(yōu)化結(jié)果,因此,其設(shè)計必須根據(jù)所要優(yōu)化的問題具體分析。對于TT&C天線而言,接收天線增益可反映接收靈敏度,而半功率波束寬度BW(又稱主瓣寬度)也對目標(biāo)測向精度有很大的影響。基于陣元激勵幅值的方向圖綜合時,在降低陣列天線峰值副瓣電平的同時,一般也會降低天線增益和半功率波束寬度,導(dǎo)致接收信號電平門限的降低和測向數(shù)據(jù)誤差的增大,因此在設(shè)計適應(yīng)度函數(shù)時需要考慮以上兩個指標(biāo)約束。

    對于式(3)表示的方向圖函數(shù),若限定最大副瓣電平與主瓣電平之差(對數(shù))為KF,構(gòu)造適應(yīng)度函數(shù),即

    (9)

    其中:max函數(shù)為求極大值,且滿足

    (10)

    (11)

    其中:BW,S為設(shè)定的半功率波束寬度,BW不應(yīng)大于BW,S;FS為設(shè)定的天線增益,F(xiàn)max不應(yīng)小于FS;α、β、γ分別為對應(yīng)項的權(quán)值,權(quán)值越大表示該項指標(biāo)越重要。通過最小化該適應(yīng)度函數(shù),不斷迭代搜索滿足約束條件的最優(yōu)解,可實現(xiàn)方向圖的低副瓣綜合。

    3.3整體算法

    對于TT&C天線,離主瓣波束最近的若干峰值副瓣電平大小對目標(biāo)角度捕獲影響最大,因此,方向圖綜合的目的在于優(yōu)化陣元激勵幅度以降低一定角度范圍的峰值副瓣電平。為保證天線方向圖E面和H面的對稱,陣元激勵幅度也應(yīng)以陣元為中心一致。稀疏陣列天線的陣元稀疏位置由矩陣A控制,故被優(yōu)化變量為MN個等效陣元激勵幅度Im,n(Am,ncm,n)。 算法實施框圖如圖3所示。具體地,將等效陣元激勵幅度Im,n(m∈[0,M-1],n∈[0,N-1])視為粒子,并進(jìn)行初始化,可采用FFT算法計算得到陣元方向圖F(θ,ψ);在此基礎(chǔ)上,即可得到峰值副瓣電平、天線增益和半功率波束寬度等指標(biāo);綜合上述3個指標(biāo)得到適應(yīng)度函數(shù)值并進(jìn)行搜索優(yōu)化,完成一次迭代;最后,達(dá)到迭代次數(shù)或者適應(yīng)度函數(shù)值滿足誤差要求后優(yōu)化完成。

    圖3 粒子群優(yōu)化方向圖框圖Fig.3 The pattern diagram of particle swarm optimization

    算法1(整體算法):

    Step1 初始化等間距柵格節(jié)點(diǎn)上的陣元幅值及其速度值。

    Step2 由FFT計算陣列方向圖和初始種群適應(yīng)度函數(shù),并根據(jù)式(11)所示的適應(yīng)度函數(shù)得到全局最優(yōu)種群和個體最優(yōu)種群。

    Step3 按式(8)產(chǎn)生新的種群,并計算新種群的適應(yīng)度函數(shù),更新全局最優(yōu)和個體最優(yōu)值,然后更新種群速度。

    Step4 判斷是否達(dá)到終止條件,若是,終止迭代,否則轉(zhuǎn)向Step 3。

    注釋1:

    (1)粒子群算法作為一種常見的啟發(fā)式算法,其明顯優(yōu)勢是算法的有效性和程序?qū)崿F(xiàn)的便捷性。

    (2)由圖3及算法1的描述可知,該整體算法可作為解決此類問題的通用框架,能夠方便地采用遺傳算法、差分進(jìn)化算法等啟發(fā)式算法替換。

    (3)適應(yīng)度函數(shù)可以靈活的采用不同指標(biāo)進(jìn)行構(gòu)造,以滿足應(yīng)用中出現(xiàn)的不同需求。

    注釋2:

    (1)本文通過分析等間距柵格平面陣列天線方向圖與離散傅里葉變換的關(guān)系,采用FFT來降低陣列方向圖計算復(fù)雜度。在整體算法1中,每個粒子的每步迭代計算均采用FFT算法降低陣列方向圖的計算量,從而顯著降低了整體算法的計算量。

    (2)根據(jù)注釋1及注釋2(1),在采用其他啟發(fā)式算法(如遺傳算法、差分進(jìn)化算法等)后,整體算法的計算量主要取決于這些啟發(fā)式算法參數(shù)(如遺傳算法的種群和遺傳代數(shù)、差分進(jìn)化算法的個體數(shù)和迭代次數(shù)等)的設(shè)置。因此,已有文獻(xiàn)對于此類啟發(fā)式算法計算復(fù)雜度的比較結(jié)果也適用于本文的整體算法1。

    4 仿真與分析

    本仿真采用的臺式計算機(jī)參數(shù)如下:酷睿i5 4590系列處理器,CPU頻率為3.3 GHz,內(nèi)存容量為4 GB。

    測控陣列天線陣元稀疏情況如圖4所示,灰色方格為陣元,白色方格表示陣元被稀疏掉,陣列中陣元行數(shù)M=16,列數(shù)N=16。下行信號空間自由波長λ,陣元間距dx=dy=2λ。

    注釋3:

    工程應(yīng)用的設(shè)備工作在S頻段,信號波長為0.13~0.16 m,設(shè)備研制方選取陣元間距dx=dy=2λ,那么單個陣元尺寸為0.27~0.35 m,平面陣列尺寸約為5 m×5 m,在這個尺寸下,既方便結(jié)構(gòu)設(shè)計、降低成本,又能滿足性能要求。因此,選取陣元間距dx=dy=2λ是綜合衡量的結(jié)果。

    4.1波束未掃描時的仿真結(jié)果

    當(dāng)不存在波束掃描,即掃描角θ=0°時,若陣元為等幅激勵,則通過Matlab仿真得到陣列天線的三維空間歸一化方向圖如圖5 所示(坐標(biāo)u′、v′定義為u′=sinθcosφ、v′=sinθsinφ),對應(yīng)的E面歸一化方向圖如圖6所示,陣列天線理論增益為65.77 dB,方向圖峰值副瓣電平為-10.25 dB,半功率波束寬度為BW=1.47°。同時可以看到,在偏離主波束30° 存在第一柵瓣。在工程實際中,考慮到拼陣損失和數(shù)字移相器幅相誤差等因素影響,實際峰值幅瓣電平和半功率波束寬度一般均大于上述仿真值。

    圖4 陣元排布示意圖Fig.4 The diagram of array element configuration

    圖6 陣列天線E 面歸一化方向圖Fig.6 The E-plane normalized pattern of array antenna

    如表1所示,測控設(shè)備天線峰值副瓣電平要求小于-13 dB,采用粒子群算法進(jìn)行方向圖低副瓣綜合時,展寬后的半功率波束寬度不大于1.6°,天線增益下降不超過2 dB??紤]到測控設(shè)備天線測角的功能需求,方向圖關(guān)于陣面法線對稱分布,因此,圖4的每個子陣中的陣元激勵幅度關(guān)于陣面中心對稱,則粒子的維數(shù)為子陣柵格數(shù)64,粒子數(shù)目n=50,陣元激勵幅度范圍為[0,1],迭代次數(shù)N=200,c1=c2=2,慣性權(quán)重ω隨迭代次數(shù)從1到0.4線性遞減,ξ、η取[0,1]的隨機(jī)數(shù),適應(yīng)度函數(shù)各項加權(quán)系數(shù)α、β、γ分別為1、0.5、1。

    表1 優(yōu)化前后參數(shù)對比Tab.1 Comparison between the non-optimized and optimized parameters

    粒子群優(yōu)化算法性能如圖7所示,第49次迭代后適應(yīng)度函數(shù)值已基本接近于0。圖8所示為粒子群算法得到的陣元激勵幅值直方圖分布,圖9 所示為陣列天線E面方向圖,優(yōu)化后的峰值副瓣電平為-13.25 dB,天線增益下降1.89 dB,波束寬度為1.55°,均滿足設(shè)定的指標(biāo)要求。另外,與陣元疊加原理計算方向圖相比,F(xiàn)FT方法計算方向圖的每次平均迭代時間由4.58 s 降為0.23 s,因此FFT 方法的運(yùn)用提高了計算效率。

    圖7 粒子群優(yōu)化算法性能Fig.7 The algorithm performance of PSO

    圖8 陣元激勵幅值直方圖分布Fig.8 The amplitude histogram of array element excitation

    圖9 E面方向圖Fig.9 The E-plane pattern

    4.2進(jìn)一步討論

    4.2.1波束掃描時的仿真結(jié)果

    采用4.1節(jié)中θ=0°時優(yōu)化得到的陣元等效激勵幅度Im,n,并保持其他參數(shù)不變,通過改變波束掃描角θ的數(shù)值來分析波束掃描時的離線優(yōu)化結(jié)果,具體如表2所示。由表2可知,隨著掃描角度θ的增大,半功率波束寬度隨著掃描角度的增大逐漸展寬;同時,峰值副瓣電平逐漸增大、天線增益逐漸下降,這兩項指標(biāo)在θ=6°時均不滿足指標(biāo)要求。顯然,如果不能針對變化的掃描角逐個優(yōu)化陣元等效激勵幅度Im,n,那么,在掃描角度θ逐步增大的過程中,峰值副瓣電平、半功率波束寬度和天線增益這3項指標(biāo)均會逐步下降。

    表2 波束掃描時的優(yōu)化結(jié)果Tab.2 Optimization results with beam scanning

    4.2.2權(quán)重影響分析

    本節(jié)對適應(yīng)度函數(shù)式(11)中參數(shù)α、β和γ的選取及其對算法的影響展開討論。

    首先,在權(quán)重α、β和γ選取時,需要考慮峰值副瓣電平、半功率波束寬度和天線增益這3項指標(biāo)的重要程度。一般來說,若要體現(xiàn)某個指標(biāo)的重要性,其對應(yīng)的權(quán)重應(yīng)該較其他指標(biāo)的權(quán)重至少高出一個數(shù)量級,下面的仿真結(jié)果也驗證了該結(jié)論。

    其次,為了量化權(quán)重α、β和γ對式(11)的影響,采用4.1節(jié)中的仿真參數(shù),僅將3項指標(biāo)的參考值(分別為峰值副瓣電平≤16.00 dB、半功率波束寬度≤1.48°、天線增益≥65.76 dB)加以收緊,收緊后的各指標(biāo)如表3所示。我們共設(shè)計4組仿真:第一組,α、β和γ均取為1,作為后續(xù)3組仿真的比較基準(zhǔn);在第二、三、四組中,α、β和γ分別取為100,表示其對應(yīng)的某項指標(biāo)重要性高于其他指標(biāo)。以第四組為例,當(dāng)α、β和γ時,天線增益較第一組中增加2.1 dB,而峰值副瓣電平則距離指標(biāo)偏差增大;也即,增大的權(quán)重值明顯體現(xiàn)了其對應(yīng)指標(biāo)在優(yōu)化過程中的重要性。類似地,從第二、三組的結(jié)果中也可得出相似的結(jié)論。

    表3 不同權(quán)重取值對3項指標(biāo)的影響Tab.3 Influence of different weight on three indexes

    需要說明的是,4.1節(jié)中對于指標(biāo)參考值的選取主要考慮了實際工程中的影響,當(dāng)然,此種情況下,各權(quán)重對指標(biāo)重要性的區(qū)分度會有所降低,但其作用仍不容忽視,這是工程應(yīng)用中通常采用式(11)作為適應(yīng)度函數(shù)的原因之一。

    5 結(jié) 論

    為解決TT&C稀疏平面陣列天線峰值副瓣電平過高問題,本文在分析相控陣天線方向圖數(shù)學(xué)模型基礎(chǔ)上,將天線增益、半功率波束寬度等指標(biāo)加入適應(yīng)度函數(shù)設(shè)計中,采用粒子群算法優(yōu)化陣元激勵幅度,使峰值副瓣電平滿足指標(biāo)要求。FFT方法的運(yùn)用加快了陣列天線方向圖計算速度,顯著降低了計算復(fù)雜度。仿真試驗驗證了算法的有效性,算法可在陣列天線設(shè)計中進(jìn)行應(yīng)用。

    [1] YEPES L F,COVARRUBIAS D H,ALONSO M A,et al. Hybrid sparse linear array synthesis applied to phased antenna arrays[J]. IEEE Antennas & Wireless Propagation Letters,2014,13(13):185-188.

    [2] 王鵬毅,孔永飛,吳海洲. 靈活的全空域同時多波束測控技術(shù)[J]. 飛行器測控學(xué)報,2013,32(01):42-47.

    WANG Pengyi,KONG Yongfei,WU Haizhou. Flexible hemispherical simultaneous multi-beam TT&C technology[J]. Journal of Spacecraft TT&C Technology,2013,187:143-153. (in Chinese)

    [3] STUTZMAN W L,THIELE G A. Antenna theory and design[M]. 3rd ed. New Jersey:John Wiley & Sons,Inc.,2005.

    [4] BROOKNER E. Practical phased-array antenna systems[M]. Boston:Artech House,1991.

    [5] HAUPT R L. Thinned arrays using genetic algorithms[J]. IEEE Transactions on Antennas & Propagation,1994,42(7):993-999.

    [6] 陳客松. 稀布天線陣列的優(yōu)化布陣技術(shù)研究[D]. 成都:電子科技大學(xué),2006.

    CHEN Kesong. The research on the sparse antenna arrays optimization disposal technology[D]. Chengdu:University of Electronic Science and Technology of China,2006. (in Chinese)

    [7] CAORSI S,LOMMI A,MASSA A,et al. Peak sidelobe level reduction with a hybrid approach based on GAs and difference sets[J]. IEEE Transactions on Antennas and Propagation,2004,52(4):1116-1121.

    [8] GATTI M S,RAHMAT-SAMII Y. FFT applications to plane-polar near-field antenna measurements[J]. IEEE Transactions on Antennas and Propagation,1988,36(6):781-791.

    [9] 焦永昌,楊科,陳勝兵,等. 粒子群優(yōu)化算法用于陣列天線方向圖綜合設(shè)計[J]. 電波科學(xué)學(xué)報,2006,21(1):16-20.

    JIAO Yongchang,YANG Ke,CHEN Shengbing,et al. Application of particle swarm optimization in antenna array pattern synthesis[J]. Chinese Journal of Radio Science,2006,21(1):16-20. (in Chinese)

    [10] 史東湖,習(xí)靖,周志偉. 對稱平面稀疏陣列陣形優(yōu)化方法研究[J]. 飛行器測控學(xué)報,2013,32(04):302-305.

    SHI Donghu,XI Jing,ZHOU Zhiwei. Research on optimization method of symmetric thinned planar array[J]. Journal of Spacecraft TT&C Technology,2013,32(4):302-305. (in Chinese)

    [11] BOERINGER D,WERNER D. Particle swarm optimization versus genetic algorithms for phased array synthesis[J]. IEEE Transactions on Antennas and Propagation,2004,52(3):771-779.

    [12] 和潔,馮大政,孟超,等. 采用遺傳算法的MIMO雷達(dá)L型陣列稀布優(yōu)化[J]. 電訊技術(shù),2015,55(8):872-878.

    HE Jie,FENG Dazheng,MENG Chao,et.al. L-shape sparse array optimization based on genetic algorithm in MIMO radar[J].Telecommunication Engineering,2015,55(8):872-878. (in Chinese)

    [13] 劉恒,趙宏偉,李維梅,等. 平面稀疏陣列天線的約束優(yōu)化設(shè)計[J]. 電訊技術(shù),2016,56(2):166-170.

    LIU Heng,ZHAO Hongwei,LI Weimei,et al. Constrained optimization of planar thinned array antenna[J]. Telecommunication Engineering,2016,56(2):166-170. (in Chinese)

    [14] 湯小為,湯俊,萬爽,等. 改進(jìn)變異策略的自適應(yīng)差分進(jìn)化算法及其應(yīng)用[J]. 宇航學(xué)報,2013,34(7):1001-1007.

    TANG Xiaowei,TANG Jun,WAN Shuang,et al. Adaptive differential evolution algorithm with modified mutation strategy and its application[J]. Journal of Astronautics,2013,34(7):1001-1007. (in Chinese)

    [15] KENNEDY J,EBERHART R. Particle swarm optimization[C]//Proceedings of IEEE International Conference on Neural Networks.Perth:IEEE,1995:1942-1948.

    [16] 李建新,陳瑞. 基于FFT平面相控陣天線方向圖快速計算[J]. 微波學(xué)報,2009,25(6):12-17.

    LI Jianxin,CHEN Rui. Fast computation of planar phased array pattern based on FFT[J]. Journal of Microwaves,2009,25(6):12-17. (in Chinese)

    [17] BROWN A.Electronically scanned arrays MATLAB modeling and simulation[J].Florida:CRC Press,2012:33-35.

    [18] 張光義,趙玉潔. 相控陣?yán)走_(dá)技術(shù)[M]. 北京:電子工業(yè)出版社,2006:24-28.

    ThinnedArrayAntennaPatternSynthesiswithMultipleConstraintsUsingFFT-PSOStrategy

    PAN Hongguang1,XUE Jikang1,GAO Hainan2,JIANG Nengqian3
    (1.School of Electrical and Control Engineering,Xi′an University of Science and Technology,Xi′an 710054,China;2.Xi′an Satellite Control Center,Weinan 714000,China;3.School of Economics and Finance,Xi′an Jiaotong University,Xi′an 710049,China)

    For the thinned array antenna,the peak sidelobe level (PSLL) of the array with equal amplitude excitation can not meet the requirements of the designed performances. Considering the constraints of antenna gain and the half-power beamwidth,the appropriate fitness function is designed,and the excitation amplitude of the element based on particle swarm optimization (PSO) is optimized to reduce the antenna sidelobe level.Meanwhile,the relationship between the pattern of planar phased array with rectangular grid and the discrete Fourier transform (DFT) is analyzed,then the fast Fourier transform (FFT) is used to reduce the computational complexity. The simulation results show that the proposed approach can effectively reduce the maximum PSLL and the computational complexity.

    space TT&C;phased array antenna;thinned planar array;particle swarm optimization;pattern synthesis

    date:2017-03-03;Revised date:2017-07-14

    國家自然科學(xué)基金資助項目(61603295);陜西省科技廳工業(yè)科技公關(guān)項目(2016GY-040,2105GY020)

    **通信作者:gaohainan1505@163.com Corresponding author:gaohainan1505@163.com

    TN821

    A

    1001-893X(2017)10-1133-07

    潘紅光(1983—),男,山東臨沂人,2015年于西安交通大學(xué)獲博士學(xué)位,期間獲國家留學(xué)基金委資助赴美聯(lián)合培養(yǎng),現(xiàn)工作于西安科技大學(xué)電控學(xué)院,主要研究方向為預(yù)測控制、智能控制及其應(yīng)用;

    Email:hongguangpan@163.com

    薛紀(jì)康(1993—),男,陜西富平人,西安科技大學(xué)電控學(xué)院碩士研究生,主要研究方向為嵌入式控制、預(yù)測控制、智能控制;

    Email:jkxue1993@163.com

    高海南(1987—),男,安徽太和人,2013年于西安交通大學(xué)獲碩士學(xué)位,現(xiàn)工作于西安衛(wèi)星測控中心,目前研究方向為無線電測控與通信、預(yù)測控制、智能控制;

    Email:gaohainan1505@163.com

    江能前(1982—),男,江蘇蘇州人,西安交通大學(xué)經(jīng)濟(jì)與金融學(xué)院博士研究生,主要研究方向為應(yīng)用經(jīng)濟(jì)、環(huán)保能源、戰(zhàn)略管理等。

    Email:jnq999@hotmail.com

    10.3969/j.issn.1001-893x.2017.10.006

    潘紅光,薛紀(jì)康,高海南,等.采用FFT-PSO策略的多約束稀疏陣天線方向圖綜合[J].電訊技術(shù),2017,57(10):1133-1139.[PAN Hongguang,XUE Jikang,GAO Hainan,et al.Thinned array antenna pattern synthesis with multiple constraints using FFT-PSO strategy[J].Telecommunication Engineering,2017,57(10):1133-1139.]

    2017-03-03;

    2017-07-14

    猜你喜歡
    波束寬度副瓣電平
    毫米波信道中波束成形矢量的波束寬度
    C波段高增益低副瓣微帶陣列天線設(shè)計
    非均勻間距的低副瓣寬帶微帶陣列天線設(shè)計
    CINRAD/SA雷達(dá)天饋系統(tǒng)關(guān)鍵參數(shù)測量方法研究
    可編程超聲波測距系統(tǒng)的設(shè)計
    科技視界(2018年9期)2018-07-27 11:28:30
    雷達(dá)副瓣跟蹤機(jī)理分析
    基于調(diào)度間隔與波束寬度良好匹配的最優(yōu)V2V毫米波通信
    汽車文摘(2017年9期)2017-12-06 05:09:19
    NPC五電平Z源逆變器的設(shè)計研究
    S波段低副瓣平面陣列天線設(shè)計
    基于三電平光伏并網(wǎng)逆變器控制系統(tǒng)的研究與實踐
    亚洲国产欧美人成| 欧美又色又爽又黄视频| 亚洲在线自拍视频| xxxwww97欧美| 午夜精品久久久久久毛片777| 国产高清有码在线观看视频 | www.www免费av| 老鸭窝网址在线观看| 亚洲九九香蕉| 九色成人免费人妻av| 一二三四社区在线视频社区8| 香蕉av资源在线| 巨乳人妻的诱惑在线观看| 欧美av亚洲av综合av国产av| 亚洲美女黄片视频| 久久人人精品亚洲av| 久久久久亚洲av毛片大全| 全区人妻精品视频| 老司机靠b影院| 老司机深夜福利视频在线观看| 亚洲午夜理论影院| 黄色片一级片一级黄色片| 狠狠狠狠99中文字幕| 欧美成人性av电影在线观看| 神马国产精品三级电影在线观看 | 一进一出好大好爽视频| 日韩精品青青久久久久久| 神马国产精品三级电影在线观看 | xxxwww97欧美| 成人18禁在线播放| 最新美女视频免费是黄的| 伦理电影免费视频| 精品国内亚洲2022精品成人| 丁香六月欧美| 国产1区2区3区精品| 国产免费男女视频| 婷婷精品国产亚洲av| 成人特级黄色片久久久久久久| 妹子高潮喷水视频| 好男人电影高清在线观看| 国产精品久久久久久亚洲av鲁大| 亚洲精品粉嫩美女一区| 哪里可以看免费的av片| 亚洲狠狠婷婷综合久久图片| 国产精品精品国产色婷婷| 精品免费久久久久久久清纯| 欧美zozozo另类| 欧美成人一区二区免费高清观看 | 国产亚洲av嫩草精品影院| 禁无遮挡网站| 精品熟女少妇八av免费久了| 99精品久久久久人妻精品| 日韩高清综合在线| 黄片大片在线免费观看| 久久九九热精品免费| 久9热在线精品视频| 亚洲一区高清亚洲精品| 亚洲国产中文字幕在线视频| 亚洲在线自拍视频| 久久精品亚洲精品国产色婷小说| 在线a可以看的网站| 一区二区三区高清视频在线| 欧美一区二区精品小视频在线| 亚洲欧美日韩无卡精品| av欧美777| 亚洲精品色激情综合| 久久久久久国产a免费观看| 一夜夜www| 国产成人影院久久av| 日本一区二区免费在线视频| 欧美日韩亚洲国产一区二区在线观看| 亚洲av成人精品一区久久| 超碰成人久久| 一边摸一边做爽爽视频免费| 国产精品影院久久| 午夜福利欧美成人| 国产91精品成人一区二区三区| 国产91精品成人一区二区三区| 90打野战视频偷拍视频| 性色av乱码一区二区三区2| 夜夜爽天天搞| 最近在线观看免费完整版| 久久精品国产99精品国产亚洲性色| 亚洲天堂国产精品一区在线| 亚洲精品av麻豆狂野| 国产乱人伦免费视频| 日韩精品青青久久久久久| 99热只有精品国产| av在线天堂中文字幕| 1024视频免费在线观看| 亚洲人成伊人成综合网2020| 丰满人妻熟妇乱又伦精品不卡| 人人妻人人看人人澡| 久久国产精品影院| 国产av在哪里看| 哪里可以看免费的av片| 久久天躁狠狠躁夜夜2o2o| 99热6这里只有精品| 国产av一区在线观看免费| 亚洲国产日韩欧美精品在线观看 | 女人被狂操c到高潮| 亚洲人成网站高清观看| 日韩国内少妇激情av| 日韩欧美 国产精品| 国产精品久久久久久精品电影| 国产精品av久久久久免费| 欧美zozozo另类| 欧美日韩一级在线毛片| 久久久久精品国产欧美久久久| 一级毛片精品| 白带黄色成豆腐渣| 色哟哟哟哟哟哟| 色在线成人网| 亚洲精品粉嫩美女一区| 十八禁人妻一区二区| 精品一区二区三区四区五区乱码| 日韩欧美国产在线观看| 看黄色毛片网站| 亚洲av五月六月丁香网| 成人国产一区最新在线观看| 在线永久观看黄色视频| 俄罗斯特黄特色一大片| 18禁裸乳无遮挡免费网站照片| 国产69精品久久久久777片 | 高潮久久久久久久久久久不卡| 中文资源天堂在线| 岛国在线观看网站| 亚洲熟女毛片儿| 少妇人妻一区二区三区视频| 精品乱码久久久久久99久播| 制服丝袜大香蕉在线| 香蕉国产在线看| www.精华液| 欧美绝顶高潮抽搐喷水| 身体一侧抽搐| 久久精品国产清高在天天线| 欧美在线一区亚洲| 麻豆成人午夜福利视频| 欧美极品一区二区三区四区| 黑人操中国人逼视频| 亚洲 欧美 日韩 在线 免费| 男人舔女人下体高潮全视频| 天天添夜夜摸| 久久久久九九精品影院| 中出人妻视频一区二区| 看片在线看免费视频| 国产精品亚洲一级av第二区| 欧美在线黄色| 国产精品av久久久久免费| 91国产中文字幕| 午夜日韩欧美国产| 国产免费男女视频| 伊人久久大香线蕉亚洲五| 国产精品免费一区二区三区在线| а√天堂www在线а√下载| x7x7x7水蜜桃| 国产av麻豆久久久久久久| 国产精品一区二区精品视频观看| 国产精品久久久久久亚洲av鲁大| 日本精品一区二区三区蜜桃| 嫩草影院精品99| 69av精品久久久久久| 三级毛片av免费| 淫妇啪啪啪对白视频| 亚洲色图av天堂| 免费在线观看完整版高清| 三级毛片av免费| www日本在线高清视频| 久久天堂一区二区三区四区| 国产单亲对白刺激| 成人av一区二区三区在线看| 久久天躁狠狠躁夜夜2o2o| 国产av麻豆久久久久久久| 丁香欧美五月| 亚洲成人久久爱视频| 两性午夜刺激爽爽歪歪视频在线观看 | 免费观看人在逋| 日韩欧美国产在线观看| 少妇人妻一区二区三区视频| 欧美日韩一级在线毛片| 国内精品一区二区在线观看| 久久久国产精品麻豆| 成熟少妇高潮喷水视频| 好男人电影高清在线观看| 亚洲一码二码三码区别大吗| 成人欧美大片| 久久亚洲精品不卡| 久久亚洲真实| 很黄的视频免费| 99久久无色码亚洲精品果冻| 操出白浆在线播放| 日韩欧美国产在线观看| 欧美精品亚洲一区二区| 99re在线观看精品视频| 淫秽高清视频在线观看| 亚洲五月婷婷丁香| 精品久久久久久久久久久久久| 久久久久久久久中文| 成人特级黄色片久久久久久久| 日韩精品免费视频一区二区三区| 在线观看美女被高潮喷水网站 | 亚洲aⅴ乱码一区二区在线播放 | 女人高潮潮喷娇喘18禁视频| 黄色视频不卡| 国产麻豆成人av免费视频| 免费在线观看视频国产中文字幕亚洲| 深夜精品福利| 又爽又黄无遮挡网站| 在线观看www视频免费| 国内毛片毛片毛片毛片毛片| 脱女人内裤的视频| 色综合亚洲欧美另类图片| 久久这里只有精品19| 在线观看66精品国产| 精品久久久久久久毛片微露脸| 午夜成年电影在线免费观看| 在线视频色国产色| 18禁黄网站禁片午夜丰满| av福利片在线观看| 我要搜黄色片| 人人妻人人看人人澡| 精品国产美女av久久久久小说| bbb黄色大片| 国语自产精品视频在线第100页| 欧美一区二区国产精品久久精品 | 国产在线观看jvid| 非洲黑人性xxxx精品又粗又长| 亚洲色图av天堂| 露出奶头的视频| 久久久久精品国产欧美久久久| av中文乱码字幕在线| 91国产中文字幕| 成人永久免费在线观看视频| 久久香蕉国产精品| 亚洲avbb在线观看| 亚洲自拍偷在线| 香蕉国产在线看| 日韩精品中文字幕看吧| 欧美日韩乱码在线| 精品第一国产精品| 一进一出抽搐gif免费好疼| 国产片内射在线| 国产人伦9x9x在线观看| 777久久人妻少妇嫩草av网站| 日韩三级视频一区二区三区| 国产又黄又爽又无遮挡在线| 国产又色又爽无遮挡免费看| 国产亚洲欧美98| 麻豆久久精品国产亚洲av| 两个人视频免费观看高清| 一级毛片精品| 长腿黑丝高跟| 两个人看的免费小视频| 在线永久观看黄色视频| 亚洲avbb在线观看| 五月伊人婷婷丁香| 欧美黄色淫秽网站| 大型黄色视频在线免费观看| 久久香蕉激情| 三级国产精品欧美在线观看 | tocl精华| 国产亚洲av高清不卡| 一个人免费在线观看电影 | 999久久久国产精品视频| a级毛片在线看网站| 国产成+人综合+亚洲专区| 十八禁网站免费在线| 久久伊人香网站| 女人高潮潮喷娇喘18禁视频| 非洲黑人性xxxx精品又粗又长| 在线观看www视频免费| 欧美在线一区亚洲| 在线十欧美十亚洲十日本专区| 国产一区二区激情短视频| 亚洲第一电影网av| 日韩精品免费视频一区二区三区| 精品午夜福利视频在线观看一区| 国产三级在线视频| 非洲黑人性xxxx精品又粗又长| 亚洲精品在线美女| 丝袜人妻中文字幕| 国产精品一区二区精品视频观看| 国产精品久久久久久人妻精品电影| 亚洲av成人一区二区三| 亚洲国产欧美人成| 丁香六月欧美| 哪里可以看免费的av片| xxx96com| 88av欧美| 欧美久久黑人一区二区| 日韩欧美三级三区| 中文字幕最新亚洲高清| 国产精品av视频在线免费观看| 桃色一区二区三区在线观看| 精品人妻1区二区| 99精品久久久久人妻精品| 亚洲熟女毛片儿| 老司机午夜福利在线观看视频| 亚洲成人久久性| 国模一区二区三区四区视频 | 啦啦啦免费观看视频1| 老司机午夜十八禁免费视频| 日韩欧美 国产精品| 免费看十八禁软件| 一二三四在线观看免费中文在| 午夜亚洲福利在线播放| 亚洲自偷自拍图片 自拍| 国产精品日韩av在线免费观看| 久久久精品大字幕| 一二三四社区在线视频社区8| 一本精品99久久精品77| 日韩欧美一区二区三区在线观看| 亚洲男人的天堂狠狠| 美女大奶头视频| tocl精华| 三级国产精品欧美在线观看 | 高潮久久久久久久久久久不卡| 精品国产美女av久久久久小说| 美女扒开内裤让男人捅视频| 国产精品香港三级国产av潘金莲| 99热这里只有精品一区 | 18禁国产床啪视频网站| 国产aⅴ精品一区二区三区波| 午夜a级毛片| 狂野欧美激情性xxxx| 国产av在哪里看| 99久久国产精品久久久| 国产精品久久久久久亚洲av鲁大| 免费在线观看亚洲国产| 国产三级在线视频| 狂野欧美激情性xxxx| 亚洲精品一区av在线观看| 成人av在线播放网站| 久久精品91蜜桃| 久久久久九九精品影院| 在线观看www视频免费| 18禁裸乳无遮挡免费网站照片| 国产精品99久久99久久久不卡| 亚洲aⅴ乱码一区二区在线播放 | 国产熟女午夜一区二区三区| 身体一侧抽搐| 久久久国产成人免费| 不卡一级毛片| 99精品久久久久人妻精品| 黄色 视频免费看| tocl精华| 精品电影一区二区在线| 精品一区二区三区四区五区乱码| 免费看a级黄色片| 午夜老司机福利片| 欧美又色又爽又黄视频| 欧美丝袜亚洲另类 | 国产黄a三级三级三级人| 夜夜躁狠狠躁天天躁| 精品乱码久久久久久99久播| 怎么达到女性高潮| 国产黄a三级三级三级人| 国产精品电影一区二区三区| 久久久久久大精品| 中文资源天堂在线| 亚洲精品中文字幕一二三四区| 女人爽到高潮嗷嗷叫在线视频| 99国产综合亚洲精品| 好看av亚洲va欧美ⅴa在| 欧美性猛交╳xxx乱大交人| netflix在线观看网站| 天天添夜夜摸| 亚洲欧美日韩东京热| 国产黄色小视频在线观看| 特级一级黄色大片| 国产高清有码在线观看视频 | 欧美性猛交╳xxx乱大交人| 黑人操中国人逼视频| 极品教师在线免费播放| 超碰成人久久| 国产高清有码在线观看视频 | 99国产综合亚洲精品| 好看av亚洲va欧美ⅴa在| 免费看美女性在线毛片视频| 亚洲国产精品合色在线| 欧美成人一区二区免费高清观看 | 亚洲精品色激情综合| 亚洲在线自拍视频| 久久久久久国产a免费观看| 一级毛片高清免费大全| 九九热线精品视视频播放| 亚洲国产欧美网| 成在线人永久免费视频| 成人高潮视频无遮挡免费网站| 亚洲精品中文字幕一二三四区| 色老头精品视频在线观看| 夜夜爽天天搞| 亚洲国产高清在线一区二区三| 久久九九热精品免费| 免费在线观看影片大全网站| 亚洲人成伊人成综合网2020| 丁香欧美五月| 黄色丝袜av网址大全| 18禁国产床啪视频网站| 啦啦啦韩国在线观看视频| 国产一区二区在线av高清观看| 一区二区三区国产精品乱码| 制服丝袜大香蕉在线| 欧美日韩瑟瑟在线播放| www.www免费av| 久久久久久免费高清国产稀缺| 欧美成狂野欧美在线观看| 老汉色∧v一级毛片| 欧美日韩福利视频一区二区| 午夜福利在线在线| 成人三级做爰电影| 日韩欧美一区二区三区在线观看| 怎么达到女性高潮| 最近视频中文字幕2019在线8| 特级一级黄色大片| 亚洲五月天丁香| 午夜两性在线视频| 一级片免费观看大全| 国产成+人综合+亚洲专区| 九色成人免费人妻av| 成年女人毛片免费观看观看9| 亚洲一区二区三区色噜噜| 51午夜福利影视在线观看| 精品久久久久久久久久免费视频| 五月伊人婷婷丁香| 18禁黄网站禁片午夜丰满| 国产视频一区二区在线看| 成人午夜高清在线视频| 亚洲av成人一区二区三| 中文字幕人妻丝袜一区二区| 性欧美人与动物交配| 国产一区二区三区视频了| 一级片免费观看大全| 欧美乱色亚洲激情| 99精品欧美一区二区三区四区| 欧美一级毛片孕妇| 麻豆久久精品国产亚洲av| 欧美黑人巨大hd| 岛国在线观看网站| 在线观看免费日韩欧美大片| ponron亚洲| 国产视频内射| 在线观看一区二区三区| 观看免费一级毛片| 90打野战视频偷拍视频| 欧美性长视频在线观看| 18禁裸乳无遮挡免费网站照片| 人妻夜夜爽99麻豆av| x7x7x7水蜜桃| 亚洲片人在线观看| 看免费av毛片| av在线天堂中文字幕| 成人国产综合亚洲| 亚洲五月婷婷丁香| 亚洲精品在线观看二区| 日韩欧美在线二视频| 日本精品一区二区三区蜜桃| 欧美日本视频| 一级片免费观看大全| 亚洲中文字幕一区二区三区有码在线看 | 久久香蕉激情| 人人妻,人人澡人人爽秒播| 国产成人欧美在线观看| 日本 av在线| 久久久久久久久免费视频了| 成在线人永久免费视频| 亚洲精品一卡2卡三卡4卡5卡| 精品欧美国产一区二区三| 成人国产一区最新在线观看| 欧美日韩亚洲国产一区二区在线观看| 亚洲精品色激情综合| 欧美一级a爱片免费观看看 | 成人国产综合亚洲| 久久午夜亚洲精品久久| 国产亚洲精品第一综合不卡| 精品久久久久久久久久久久久| 欧美极品一区二区三区四区| 欧美黄色片欧美黄色片| 亚洲国产日韩欧美精品在线观看 | 免费看美女性在线毛片视频| 亚洲专区国产一区二区| 一进一出抽搐动态| 国产不卡一卡二| 最近最新中文字幕大全电影3| 床上黄色一级片| 欧美 亚洲 国产 日韩一| 国产精品久久视频播放| 大型黄色视频在线免费观看| 校园春色视频在线观看| 1024香蕉在线观看| 亚洲成av人片在线播放无| 亚洲黑人精品在线| 久久婷婷成人综合色麻豆| 日本撒尿小便嘘嘘汇集6| 亚洲片人在线观看| 国产av麻豆久久久久久久| 亚洲精品中文字幕一二三四区| 国产精品久久久久久久电影 | 久久久国产成人免费| 精品久久久久久,| 中文在线观看免费www的网站 | 亚洲色图av天堂| 久久香蕉激情| 国产真人三级小视频在线观看| 成人18禁高潮啪啪吃奶动态图| 淫秽高清视频在线观看| 国产精品98久久久久久宅男小说| 可以在线观看的亚洲视频| 国产一区在线观看成人免费| 久久久久久久久免费视频了| 我要搜黄色片| 精品国产超薄肉色丝袜足j| 国产区一区二久久| 日韩欧美三级三区| 欧美乱色亚洲激情| 欧美乱妇无乱码| 天天躁狠狠躁夜夜躁狠狠躁| 黑人欧美特级aaaaaa片| 制服丝袜大香蕉在线| 夜夜看夜夜爽夜夜摸| 99久久精品国产亚洲精品| 国产免费av片在线观看野外av| 波多野结衣巨乳人妻| 人成视频在线观看免费观看| 午夜福利高清视频| 亚洲avbb在线观看| 久久久久久国产a免费观看| 手机成人av网站| 99精品久久久久人妻精品| 中文字幕人成人乱码亚洲影| 窝窝影院91人妻| 欧美在线一区亚洲| 亚洲 国产 在线| 日韩欧美在线二视频| 亚洲av日韩精品久久久久久密| 亚洲一区中文字幕在线| 久99久视频精品免费| 国内揄拍国产精品人妻在线| 一级毛片高清免费大全| 国产精品久久视频播放| 97碰自拍视频| 一本精品99久久精品77| 久久人妻福利社区极品人妻图片| 精品一区二区三区视频在线观看免费| 最新在线观看一区二区三区| 免费在线观看完整版高清| 又黄又粗又硬又大视频| 好男人在线观看高清免费视频| 国产真人三级小视频在线观看| 成人一区二区视频在线观看| 欧美在线一区亚洲| 嫩草影视91久久| 国产精品 欧美亚洲| 90打野战视频偷拍视频| 一边摸一边做爽爽视频免费| 成人国语在线视频| 老汉色av国产亚洲站长工具| 成年版毛片免费区| 国内精品久久久久久久电影| 91麻豆精品激情在线观看国产| 99久久99久久久精品蜜桃| 十八禁人妻一区二区| 久久人人精品亚洲av| 国产成人av激情在线播放| 男女之事视频高清在线观看| 日韩欧美 国产精品| 国产久久久一区二区三区| 九九热线精品视视频播放| 精品少妇一区二区三区视频日本电影| 国产亚洲欧美98| 十八禁网站免费在线| 欧美不卡视频在线免费观看 | av欧美777| 日韩欧美在线乱码| 禁无遮挡网站| 国产又色又爽无遮挡免费看| 中文字幕高清在线视频| 一本久久中文字幕| 精华霜和精华液先用哪个| 欧美最黄视频在线播放免费| 最近最新免费中文字幕在线| 国产野战对白在线观看| 一级毛片女人18水好多| 老鸭窝网址在线观看| 免费看a级黄色片| 夜夜躁狠狠躁天天躁| 国产精品久久久久久人妻精品电影| 亚洲av成人不卡在线观看播放网| 欧美大码av| 狠狠狠狠99中文字幕| 成人国产一区最新在线观看| 最近最新中文字幕大全电影3| а√天堂www在线а√下载| 天天躁夜夜躁狠狠躁躁| 久久午夜综合久久蜜桃| 国产亚洲精品久久久久5区| 久久婷婷成人综合色麻豆| 国产久久久一区二区三区| 久久中文字幕一级| 国产精品综合久久久久久久免费| 色尼玛亚洲综合影院| 少妇裸体淫交视频免费看高清 | 久久久久国内视频| 亚洲国产精品999在线| 国产探花在线观看一区二区| 国产精品久久视频播放| 亚洲国产精品999在线| 国产激情偷乱视频一区二区| 久久精品成人免费网站| 搡老岳熟女国产| 免费观看精品视频网站| 麻豆av在线久日| 18禁黄网站禁片免费观看直播| 国产成人精品久久二区二区免费| 91国产中文字幕|