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

    基于自適應(yīng)變鄰域搜索的火控雷達(dá)組網(wǎng)資源調(diào)度研究

    2025-03-20 00:00:00汪達(dá)旺陸志灃伍國華

    摘 要:在海上艦艇編隊(duì)防空作戰(zhàn)中,傳統(tǒng)火控雷達(dá)由于實(shí)行單目標(biāo)對單雷達(dá)跟蹤策略,難以滿足復(fù)雜的現(xiàn)代作戰(zhàn)需求。針對該問題,提出一種基于自適應(yīng)變鄰域搜索(adaptive variable neighborhood search, AVNS)的火控雷達(dá)組網(wǎng)資源調(diào)度方法。以海上艦艇編隊(duì)火控雷達(dá)資源為基礎(chǔ),考慮雷達(dá)剩余通道數(shù)、雷達(dá)探測范圍和預(yù)測航跡跟蹤覆蓋率約束,構(gòu)建以資源目標(biāo)距離、有效航跡覆蓋率和目標(biāo)威脅度為目標(biāo)函數(shù)的資源調(diào)度模型,設(shè)計(jì)基于AVNS的組網(wǎng)策略算法。在設(shè)置的兩種海上艦艇編隊(duì)模擬隊(duì)形場景中,與傳統(tǒng)規(guī)則算法和變鄰域搜索(variable neighborhood search, VNS)算法對比,驗(yàn)證所提方法的合理性和有效性。

    關(guān)鍵詞: 海上防空; 資源調(diào)度; 自適應(yīng)變鄰域搜索; 雷達(dá)組網(wǎng); 火控雷達(dá)

    中圖分類號(hào): E 955

    文獻(xiàn)標(biāo)志碼: ADOI:10.12305/j.issn.1001 506X.2025.02.16

    Research on resource scheduling of fire control radar networking based on

    adaptive variable neighborhood search

    WANG Dawang1, LU Zhifeng2, WU Guohua1,*

    (1. School of Automation Electro Mechanical Engineering Institute, Central South University, Changsha 410075, China;

    2. Shanghai Electro Mechanical Engineering Institute, Shanghai 201109, China)

    Abstract:In naval warship formation air defense combats, traditional fire control radars, due to its implemen tation of a single target to single radar tracking strategy, are difficult to meet complex modern combat needs." To address this issue, a fire control radar networking resource scheduling method based on adaptive variable neighborhood search (AVNS) is proposed. Based on the fire control radar resources of naval warship forma tions, considering the constraints of remaining radar channels, radar detection range, and predicted trajectory tracking coverage rate, a resource scheduling model is constructed with the objective functions of resource target distance, effective trajectory coverage rate, and target threat degree. A networking strategy algorithm based on AVNS is designed. The rationality and effectiveness of the proposed method is verified by comparing it with traditional rule algorithms and variable neighborhood search (VNS) algorithms in two naval warship formation simulated formation scenarios.

    Keywords:naval air defense; resource scheduling; adaptive variable neighborhood search (AVNS); radar networking; fire control radar

    0 引 言

    相控陣火控雷達(dá)(以下簡稱火控雷達(dá))是用于搜索、發(fā)現(xiàn)和跟蹤空中目標(biāo),連續(xù)測定目標(biāo)坐標(biāo),解算射擊諸元并控制火力攔截系統(tǒng)射擊的雷達(dá)1。由于艦船的空間有限、功能多、靈活性高,掃描速度快的火控雷達(dá)越來越廣泛地被各國海軍部署在艦艇上2?,F(xiàn)代軍艦中火控雷達(dá)和艦炮協(xié)同工作,由雷達(dá)發(fā)現(xiàn)目標(biāo)并進(jìn)行跟蹤,向控制系統(tǒng)傳送目標(biāo)信息,艦炮負(fù)責(zé)打擊目標(biāo)3。由此可見,火控雷達(dá)目標(biāo)探測跟蹤是現(xiàn)代與未來防空的首要任務(wù),也是火控系統(tǒng)的關(guān)鍵環(huán)節(jié)4,因此針對目標(biāo)探測跟蹤的火控雷達(dá)資源調(diào)度問題的研究非常重要。

    傳統(tǒng)的火控雷達(dá)調(diào)度方法由雷達(dá)工作方式的優(yōu)先級(jí)決定,并在調(diào)度過程中保持不變,自適應(yīng)能力差5。“宙斯盾”總設(shè)計(jì)師Baugh闡述相控陣?yán)走_(dá)資源調(diào)度的原理,分析調(diào)度策略的特點(diǎn)和影響調(diào)度算法設(shè)計(jì)的因素6。文獻(xiàn)[7]建立了一種搜索和跟蹤任務(wù)綜合的雙目標(biāo)約束優(yōu)化框架,并使用帕累托理論求解。文獻(xiàn)[8]基于工作指標(biāo)和任務(wù)影響提出一種滿載條件下的自適應(yīng)調(diào)度算法。文獻(xiàn)[9]采用插空法與帶有時(shí)間窗的一步回溯法相結(jié)合的自適應(yīng)算法求解。

    雖然大多數(shù)方法能夠滿足火控雷達(dá)資源調(diào)度需求,但是在海上防空作戰(zhàn)中,艦艇編隊(duì)往往需要協(xié)同作戰(zhàn),對編隊(duì)整體進(jìn)行雷達(dá)資源優(yōu)化調(diào)度10。艦艇編隊(duì)協(xié)同探測克服單艦艦載雷達(dá)系統(tǒng)因地球曲率影響而使探測視距受限的問題,通過采用不同艦船平臺(tái)雷達(dá)系統(tǒng)的協(xié)同工作機(jī)制,使其在復(fù)雜自然環(huán)境下的適應(yīng)能力、復(fù)雜電磁環(huán)境下的試驗(yàn)?zāi)芰σ约拔⑷跄繕?biāo)檢測能力等得到很大提升11。在協(xié)同作戰(zhàn)方式下,通過多平臺(tái)傳感器協(xié)同探測,擴(kuò)大探測范圍,可以更快、更可靠、更精確地發(fā)現(xiàn)目標(biāo),實(shí)現(xiàn)發(fā)現(xiàn)即打擊能力12-13;火控雷達(dá)數(shù)據(jù)共享,增強(qiáng)對作戰(zhàn)空間的一致性理解,支持火力協(xié)同打擊功能,網(wǎng)內(nèi)任意火力單元可利用其他單元火控雷達(dá)的數(shù)據(jù)進(jìn)行射擊,以實(shí)施最佳的打擊方案14。這種艦艇編隊(duì)火控雷達(dá)協(xié)同作戰(zhàn)的方式被稱為火控雷達(dá)組網(wǎng),火控雷達(dá)組網(wǎng)資源調(diào)度問題可以視為組合優(yōu)化問題15

    在海上防空作戰(zhàn)過程中,防空艦艇編隊(duì)必須使用火控雷達(dá)對來襲導(dǎo)彈進(jìn)行跟蹤,才能為下一步火力攔截提供實(shí)時(shí)的目標(biāo)狀態(tài)信息16。當(dāng)導(dǎo)彈來襲時(shí),導(dǎo)彈往往出現(xiàn)在距離艦艇很遠(yuǎn)的位置,而艦載火控雷達(dá)的探測范圍有限,無法跟蹤此時(shí)的來襲導(dǎo)彈17。為了彌補(bǔ)艦載雷達(dá)在探測距離上的“弱勢”,可以借助其他作戰(zhàn)單元的火控雷達(dá)進(jìn)行協(xié)同探測18-19。這里的作戰(zhàn)單元包含艦艇編隊(duì)的艦艇和預(yù)警機(jī)。預(yù)警機(jī)由于其有超遠(yuǎn)視距優(yōu)勢,可以協(xié)助艦載雷達(dá)共同完成對來襲導(dǎo)彈的長距離跟蹤20。如圖1所示,當(dāng)來襲目標(biāo)進(jìn)入機(jī)載雷達(dá)探測范圍后,機(jī)載雷達(dá)將對其進(jìn)行實(shí)時(shí)跟蹤,并占用機(jī)載雷達(dá)的一個(gè)通道數(shù)。當(dāng)其脫離機(jī)載雷達(dá)的探測范圍又進(jìn)入艦載雷達(dá)的探測范圍時(shí),機(jī)載雷達(dá)取消探測,釋放雷達(dá)通道,轉(zhuǎn)而由艦載雷達(dá)接力完成對目標(biāo)的跟蹤探測,并占用艦載雷達(dá)的一個(gè)通道數(shù)。這一過程被稱為雷達(dá)組網(wǎng)探測,簡稱雷達(dá)組網(wǎng)。利用雷達(dá)組網(wǎng),可以使艦艇編隊(duì)對目標(biāo)的探測范圍大大擴(kuò)展,進(jìn)而獲取準(zhǔn)確、充分的目標(biāo)信息,為火力攔截目標(biāo)提供信息支撐21。

    為應(yīng)對不斷變化的空中威脅,加強(qiáng)海上區(qū)域防空能力,美國海軍在“網(wǎng)絡(luò)中心戰(zhàn)”思想的指導(dǎo)下,提出并建設(shè)“海上一體化防空火控”體系,2015年在“羅斯福”號(hào)航母上形成初始作戰(zhàn)能力,現(xiàn)階段正在美航母戰(zhàn)斗群中廣泛推進(jìn)部署,不斷推陳出新,發(fā)展新構(gòu)型,以實(shí)現(xiàn)更大范圍和更具靈活性的一體化作戰(zhàn)22。該體系通過聯(lián)合E 2D預(yù)警機(jī)、宙斯盾防御系統(tǒng)和標(biāo)準(zhǔn) 6艦空導(dǎo)彈,建立起一個(gè)在預(yù)警機(jī)支援下的艦機(jī)協(xié)同反導(dǎo)打擊鏈。該打擊鏈的建立極大地縮短了導(dǎo)彈武器打擊目標(biāo)的反應(yīng)時(shí)間,使美國海軍形成了基于E 2D預(yù)警機(jī)的協(xié)同作戰(zhàn)能力23。因此,對將預(yù)警機(jī)資源加入火控雷達(dá)組網(wǎng)資源調(diào)度問題進(jìn)行研究非常必要。

    在現(xiàn)有研究中,火控雷達(dá)組網(wǎng)資源調(diào)度問題通常以跟蹤目標(biāo)總數(shù)為目標(biāo)函數(shù),以雷達(dá)資源和雷達(dá)探測范圍等為約束條件。文獻(xiàn)[24]以目標(biāo)跟蹤數(shù)目和雷達(dá)資源的最大化為優(yōu)化目標(biāo),利用隱枚舉法得到對雷達(dá)組網(wǎng)中多個(gè)雷達(dá)的資源調(diào)度方案。文獻(xiàn)[25]經(jīng)驗(yàn)證發(fā)現(xiàn)協(xié)調(diào)雷達(dá)資源管理相比獨(dú)立雷達(dá)資源管理提高了跟蹤性能。文獻(xiàn)[26]通過共享雷達(dá)資源減少目標(biāo)跟蹤負(fù)荷,增強(qiáng)目標(biāo)搜索能力和提高目標(biāo)跟蹤質(zhì)量。文獻(xiàn)[27]利用多種混合策略對基于分解的多目標(biāo)粒子群算法進(jìn)行改進(jìn),使用該算法對多功能雷達(dá)組網(wǎng)資源調(diào)度方案進(jìn)行優(yōu)化。文獻(xiàn)[28]使用目標(biāo)動(dòng)態(tài)威脅度驅(qū)動(dòng)的波束分配與駐留時(shí)間聯(lián)合優(yōu)化算法解決分布式組網(wǎng)相控陣?yán)走_(dá)多目標(biāo)跟蹤問題。此外,雷達(dá)組網(wǎng)還需要處理目標(biāo)預(yù)測航跡點(diǎn)29。雷達(dá)組網(wǎng)所需的點(diǎn)跡融合技術(shù)主要有點(diǎn)跡壓縮合并方法和點(diǎn)跡串行合并方法30。

    當(dāng)前對火控雷達(dá)協(xié)同作戰(zhàn)方式的算法研究尚淺,傳統(tǒng)規(guī)則算法仍為主流。但隨著資源和目標(biāo)數(shù)量的提升,傳統(tǒng)規(guī)則算法對雷達(dá)資源利用率考慮不足的短板逐漸顯露。因此,利用智能優(yōu)化算法實(shí)現(xiàn)火控雷達(dá)協(xié)同作戰(zhàn)很有必要31。針對海上艦艇編隊(duì)防空作戰(zhàn)場景下火控雷達(dá)資源調(diào)度問題,提出一種基于自適應(yīng)變鄰域搜索(adaptive variable neighborhood search, AVNS)的火控雷達(dá)組網(wǎng)資源調(diào)度方法,為艦艇編隊(duì)的來襲目標(biāo)跟蹤任務(wù)提供快速、高效的雷達(dá)組網(wǎng)調(diào)度方案。主要貢獻(xiàn)如下:

    (1) 根據(jù)雷達(dá)資源位置、雷達(dá)資源探測范圍、目標(biāo)預(yù)測航跡和預(yù)測航跡最短跟蹤覆蓋距離進(jìn)行可用航跡修正,將剩余可用航跡段縮放至同等大小,保證多目標(biāo)的數(shù)據(jù)結(jié)構(gòu)一致性,并提高計(jì)算效率。

    (2) 將預(yù)警機(jī)雷達(dá)資源同艦載雷達(dá)資源一起作為艦艇編隊(duì)整體雷達(dá)資源,更加符合現(xiàn)代作戰(zhàn)需求??紤]雷達(dá)資源剩余通道數(shù)量、修正航跡、跟蹤目標(biāo)總數(shù)、有效航跡覆蓋率、目標(biāo)威脅度和目標(biāo)跟蹤任務(wù)規(guī)劃需求,將火控雷達(dá)資源調(diào)度問題表述為組合優(yōu)化問題,建立火控雷達(dá)資源調(diào)度模型。

    (3) 設(shè)計(jì)3種鄰域結(jié)構(gòu)和AVNS策略,通過局部搜索算法提高模型尋優(yōu)效率,能在較短時(shí)間內(nèi)給出較優(yōu)方案,符合防空作戰(zhàn)對時(shí)效性的要求。

    (4) 以海上艦艇編隊(duì)防空作戰(zhàn)為背景,針對火控雷達(dá)組網(wǎng)資源調(diào)度問題,在10種來襲目標(biāo)跟蹤需求和兩種模擬海上艦艇編隊(duì)場景中,驗(yàn)證基于AVNS的火控雷達(dá)組網(wǎng)資源調(diào)度算法的有效性。

    1 火控雷達(dá)資源調(diào)度的數(shù)學(xué)模型

    將火控雷達(dá)資源調(diào)度問題的目標(biāo)函數(shù)構(gòu)建為最大有效航跡覆蓋率、最小目標(biāo)距離和最大目標(biāo)威脅度的3目標(biāo)函數(shù)。此外,考慮雷達(dá)資源探測范圍、有效航跡最短跟蹤覆蓋距離和雷達(dá)剩余通道數(shù)量3個(gè)約束。構(gòu)建一個(gè)3約束、3目標(biāo)的最大化多目標(biāo)問題模型。

    1.1 個(gè)體編碼與航跡點(diǎn)表示方法

    假設(shè)有一支由I個(gè)作戰(zhàn)單元組成的海上艦艇編隊(duì),其中有I1架預(yù)警機(jī)和I2艘水面艦艇,且在編隊(duì)組成中僅考慮有可用火控雷達(dá)資源且能夠參與防空目標(biāo)跟蹤任務(wù)的作戰(zhàn)單元。由于所有作戰(zhàn)單元均使用相控陣?yán)走_(dá),故可實(shí)現(xiàn)對周圍區(qū)域360°的跟蹤探測16。預(yù)警機(jī)與艦艇使用同一套艦載防空系統(tǒng),艦艇和預(yù)警機(jī)會(huì)根據(jù)自身位置和職責(zé),承擔(dān)不同區(qū)域的跟蹤任務(wù),而且預(yù)警機(jī)雷達(dá)的探測范圍大于艦艇雷達(dá)的探測范圍。

    假設(shè)有J個(gè)來襲目標(biāo),每個(gè)目標(biāo)的初始預(yù)測航跡點(diǎn)總數(shù)為M,給出第j個(gè)目標(biāo)的第m個(gè)預(yù)測三維航跡點(diǎn)(Xjm,Yjm,Zjm),j∈[1,J],m∈[0,M]。當(dāng)m=0時(shí),代表目標(biāo)當(dāng)前位置。對每一個(gè)目標(biāo),在同一時(shí)間只需要一部火控雷達(dá)進(jìn)行跟蹤。

    如圖2所示,用網(wǎng)格坐標(biāo)來表征艦機(jī)雷達(dá)位置,資源狀態(tài)信息中的艦船和預(yù)警機(jī)位置代表其裝備雷達(dá)的坐標(biāo)(Xi,Yi,Zi),i∈[1,I],(Xi,Yi,Zi),i∈[1,I]。假定使用的艦載雷達(dá)型號(hào)為SPY 1D,使用的機(jī)載雷達(dá)型號(hào)為AN/APY 9,網(wǎng)格內(nèi)通道數(shù)非空,即代表相應(yīng)位置的艦艇或預(yù)警機(jī)有可用火控雷達(dá)資源,形成帶有不同元素值的規(guī)則矩陣,以表示編隊(duì)各艦艇或者預(yù)警機(jī)裝備雷達(dá)的具體資源情況。同樣,使用目標(biāo)態(tài)勢信息中的目標(biāo)位置(Xj,Yj,Zj),j∈[1,J],(Xi,Yi,Zi),i∈[1,I]代表各目標(biāo)坐標(biāo)。

    此外,從接受跟蹤指令到計(jì)算完成輸出,導(dǎo)致仿真雷達(dá)系統(tǒng)實(shí)施存在一定時(shí)間延遲,需要將這段延遲時(shí)間Td也考慮在雷達(dá)組網(wǎng)跟蹤時(shí)間內(nèi)。

    1.2 目標(biāo)函數(shù)

    將火控雷達(dá)資源調(diào)度問題視為一個(gè)綜合考慮總距離優(yōu)化函數(shù)f(D)、總跟蹤有效航跡覆蓋率函數(shù)f(G)和總目標(biāo)威脅度判定函數(shù)f(W)最小的3目標(biāo)函數(shù)最小化問題,即

    min Z=f(D)+γf(G)+ρf(W)(1)

    式中:γ是平衡有效航跡覆蓋率函數(shù)f(G)與目標(biāo)函數(shù)Z的懲罰參數(shù),γ∈[1,2];ρ是平衡威脅度判定函數(shù)f(W)與目標(biāo)函數(shù)Z的懲罰參數(shù),ρ∈[1,2]。根據(jù)兩個(gè)懲罰參數(shù)調(diào)節(jié)各函數(shù)和目標(biāo)函數(shù)間的平衡關(guān)系,并根據(jù)多次仿真設(shè)定實(shí)驗(yàn)標(biāo)準(zhǔn)值,通過降低目標(biāo)函數(shù)Z來選擇更優(yōu)的雷達(dá)組網(wǎng)資源調(diào)度方案;f(D)代表雷達(dá)組網(wǎng)資源調(diào)度方案的距離優(yōu)化函數(shù)。Dj代表待跟蹤目標(biāo)j的距離優(yōu)化函數(shù);Di-jn代表第i部雷達(dá)與待跟蹤目標(biāo)j的第n個(gè)航跡點(diǎn)的距離。在對雷達(dá)組網(wǎng)調(diào)度方案進(jìn)行優(yōu)化選擇時(shí),假設(shè)此時(shí)對第j個(gè)目標(biāo)的雷達(dá)組網(wǎng)跟蹤方案如下:i1雷達(dá)跟蹤第j個(gè)目標(biāo)航跡第1至第n點(diǎn);i2雷達(dá)跟蹤目標(biāo)航跡第n+1點(diǎn)至第m點(diǎn);i1雷達(dá)和i2雷達(dá)的坐標(biāo)分別為(Xi1,Yi1,Zi1)和(Xi2,Yi2,Zi2);第j個(gè)目標(biāo)航跡第1點(diǎn)坐標(biāo)為(Xj1,Yj1,Zj1),第n點(diǎn)坐標(biāo)為(Xjn,Yjn,Zjn),第n+1點(diǎn)坐標(biāo)為(Xj(n+1),Yj(n+1),Zj(n+1)),第m點(diǎn)坐標(biāo)為(Xjm,Yjm,Zjm)。那么i1雷達(dá)對與第j個(gè)目標(biāo)需要跟蹤航跡點(diǎn)的起點(diǎn)和終點(diǎn)的距離分別為

    Di1-j1=2(Xi1-Xj1)2+(Yi1-Yj1)2+(Zi1-Zj1)2(2)

    Di1-jn=2(Xi1-Xjn)2+(Yi1-Yjn)2+(Zi1-Zjn)2(3)

    i2雷達(dá)與第j個(gè)目標(biāo)需要跟蹤航跡點(diǎn)的起點(diǎn)和終點(diǎn)的距離分別為

    Di2-j(n+1)=

    2(Xi2-Xj(n+1))2+(Yi2-Yj(n+1))2+(Zi2-Zj(n+1))2(4)

    Di2-jm=2(Xi2-Xjm)2+(Yi2-Yjm)2+(Zi2-Zjm)2(5)

    i1和i2雷達(dá)與第j個(gè)目標(biāo)的距離優(yōu)化函數(shù)定義如下:

    Di1-j=2D2i1-j1+D2i1-jn(6)

    Di2-j=2D2i2-j(n+1)+D2i2-jm(7)

    則對第j個(gè)目標(biāo)的雷達(dá)組網(wǎng)資源調(diào)度距離優(yōu)化函數(shù)為

    D(j)=Di1-j+Di2-j(8)

    對于所有目標(biāo)的整體雷達(dá)組網(wǎng)資源調(diào)度方案的距離優(yōu)化函數(shù)為

    f(D)=∑Jj=1D(j)(9)

    此外,還要考慮總跟蹤有效航跡覆蓋率最大,如下式所示:

    Gj=gjg′j(10)

    式(10)是第j個(gè)目標(biāo)有效航跡覆蓋率Gj的計(jì)算公式。gj表示對第j個(gè)目標(biāo)分配的跟蹤航跡點(diǎn)數(shù)量;g′j表示對第j個(gè)目標(biāo)的最長有效跟蹤航跡點(diǎn)數(shù)量。

    則總體方案跟蹤有效航跡覆蓋率函數(shù)如下式所示:

    f(G)=(J-J-)∑Ji=0(1-Gj)(11)

    式(11)表示在優(yōu)化過程中,期望雷達(dá)組網(wǎng)方案能夠跟蹤的目標(biāo)總數(shù)J-越大越好,對每個(gè)跟蹤目標(biāo)的有效跟蹤航跡覆蓋率越大越好。接下來,還需要考慮目標(biāo)的威脅值:

    f(W)=∑J-j=01Wj(12)

    式(12)是在雷達(dá)組網(wǎng)資源調(diào)度方案中待跟蹤目標(biāo)的威脅度判定函數(shù)。Wj表示第j個(gè)目標(biāo)的威脅度,Wj取1~10的自然數(shù),取值越大表示目標(biāo)威脅程度越高。式(12)的值越小,表示所有待跟蹤目標(biāo)的威脅度判定函數(shù)越小,雷達(dá)組網(wǎng)資源調(diào)度方案中威脅度高的待跟蹤目標(biāo)數(shù)量越多。

    1.3 約束條件

    在模型中,考慮雷達(dá)目標(biāo)探測距離、最短跟蹤航跡和雷達(dá)通道數(shù)量共3個(gè)約束。

    (1) 雷達(dá)目標(biāo)探測距離約束

    在編隊(duì)下達(dá)雷達(dá)目標(biāo)跟蹤指令時(shí),會(huì)為每個(gè)目標(biāo)分配火力攔截的目標(biāo)艦船,即每個(gè)目標(biāo)由一一對應(yīng)的艦艇負(fù)責(zé)火力攔截任務(wù)。已知每個(gè)目標(biāo)所對應(yīng)的任務(wù)艦艇,首先對預(yù)測目標(biāo)航跡點(diǎn)進(jìn)行篩選。由于每艘艦艇的各型導(dǎo)彈打擊范圍不同,取該艦艇所有類型導(dǎo)彈的打擊范圍的上下值作為艦艇打擊范圍的最大最小值

    Himinlt;Dij-mlt;Himax(13)

    式中:Himin為第i艘艦船導(dǎo)彈打擊范圍的最小值;Himax為第i艘艦船導(dǎo)彈打擊范圍的最大值;Dij-m為第j個(gè)目標(biāo)的第m個(gè)預(yù)測航跡點(diǎn)與對應(yīng)的第i艘艦艇的距離。通過式(13),可以得到滿足目標(biāo)艦艇打擊范圍的預(yù)測航跡點(diǎn),這些航跡點(diǎn)被稱為有效航跡點(diǎn)。

    (2) 最短跟蹤航跡約束

    考慮到雷達(dá)組網(wǎng)跟蹤是一個(gè)長時(shí)動(dòng)作,需要摒棄掉過短的航跡點(diǎn)連線,否則跟蹤時(shí)間過短,無法準(zhǔn)確制導(dǎo),將造成資源浪費(fèi)。

    ∑Mm=1ε(Dij-m-Himin)·ε(Himax-Dij-m)≥K

    (14)

    式中:i∈[0,I];j∈[0,J];ε(x)為階躍函數(shù);K為最少航跡點(diǎn)數(shù);M為預(yù)估航跡點(diǎn)總數(shù)。

    (3) 雷達(dá)通道數(shù)量約束

    對應(yīng)于每一艘艦艇或預(yù)警機(jī),要求在同一時(shí)刻,雷達(dá)組網(wǎng)方案分配給該艦機(jī)的跟蹤目標(biāo)數(shù)量不超過該艦機(jī)的雷達(dá)通道資源總數(shù)。由于單部雷達(dá)測定能力有限,如當(dāng)SPY 1D型號(hào)的火控跟蹤雷達(dá)的剩余通道數(shù)量為3時(shí),表示其能夠同時(shí)跟蹤的目標(biāo)數(shù)為3,故雷達(dá)組網(wǎng)資源調(diào)度方案需要考慮各雷達(dá)的剩余可用通道數(shù)量,即

    ∑Jj=0sij≤Si

    s.t. sij=0, 第i部雷達(dá)沒有跟蹤目標(biāo)j的任務(wù)

    sij=1, 第i部雷達(dá)有跟蹤目標(biāo)j的任務(wù)(15)

    式中:Si表示第i部雷達(dá)的通道資源總數(shù)。

    2 基于AVNS的火控雷達(dá)資源調(diào)度算法

    對輸入數(shù)據(jù)進(jìn)行目標(biāo)航跡篩選,刪除有效航跡點(diǎn)個(gè)數(shù)不滿足編隊(duì)雷達(dá)資源跟蹤時(shí)長的目標(biāo)。在篩選目標(biāo)預(yù)測航跡點(diǎn)時(shí),預(yù)測航跡點(diǎn)有400個(gè),數(shù)量多,影響計(jì)算速度。首先根據(jù)仿真雷達(dá)探測范圍,將不在任意艦機(jī)探測范圍內(nèi)的航跡點(diǎn)刪除。當(dāng)單個(gè)目標(biāo)預(yù)測航跡點(diǎn)中的有效航跡點(diǎn)個(gè)數(shù)少于20時(shí),認(rèn)為該航跡過短,不滿足雷達(dá)跟蹤條件,故放棄對該目標(biāo)的雷達(dá)跟蹤資源調(diào)度。再對篩選完的有效航跡點(diǎn)進(jìn)行取整縮小處理,將剩余航跡點(diǎn)個(gè)數(shù)取整縮小至20,用統(tǒng)一長度為20的列表表示,有利于后續(xù)鄰域變換和數(shù)據(jù)處理。對每個(gè)目標(biāo),列表中的數(shù)值為雷達(dá)編號(hào),連續(xù)的同等數(shù)值表示該段航跡均由相同的雷達(dá)跟蹤。對包含所有目標(biāo)的列表進(jìn)行AVNS,加快搜索效率。

    基于變鄰域搜索(variable neighborhood search, VNS)的火控雷達(dá)組網(wǎng)資源調(diào)度算法設(shè)計(jì)了3種鄰域變換結(jié)構(gòu)。第一種鄰域變換結(jié)構(gòu)將同列表中前后兩段不同數(shù)值的值段交換;第二種鄰域變換結(jié)構(gòu)將某一值段整體替換為其他值段;第三種鄰域變換結(jié)構(gòu)將某一值段的首端或末端延長格的值同化為該值段的值。當(dāng)目標(biāo)跟蹤方案中前后兩個(gè)值段的雷達(dá)編號(hào)不同時(shí),可由第一種鄰域變換結(jié)構(gòu)交換雷達(dá)編號(hào),得到新的雷達(dá)跟蹤方案。第二種鄰域變換結(jié)構(gòu)可以將目標(biāo)跟蹤方案中某一值段的雷達(dá)編號(hào)改變。第三種鄰域變換結(jié)構(gòu)可以將目標(biāo)跟蹤方案中處于值段末端的雷達(dá)編號(hào)改為與之相鄰的雷達(dá)編號(hào),并且當(dāng)整個(gè)列表均被同一雷達(dá)編號(hào)覆蓋時(shí),可以將值段末端隨機(jī)變化為其余雷達(dá)編號(hào)。3種鄰域結(jié)構(gòu)可以滿足對目標(biāo)列表的復(fù)雜變換,從而找到更優(yōu)的雷達(dá)組網(wǎng)跟蹤方案。3種鄰域結(jié)構(gòu)如圖3所示。

    由于VNS只能盲目地將所有鄰域結(jié)構(gòu)依次搜索一遍,缺乏搜索積累的經(jīng)驗(yàn)指導(dǎo),并且增加時(shí)間成本,不滿足防空作戰(zhàn)中火控雷達(dá)資源調(diào)度對時(shí)效性的要求,而自適應(yīng)機(jī)制則解決了上述問題。在上述3種鄰域結(jié)構(gòu)的基礎(chǔ)上,采用AVNS策略,可以提高搜索效率,在短時(shí)間內(nèi)得到更優(yōu)的雷達(dá)組網(wǎng)跟蹤方案。根據(jù)自適應(yīng)機(jī)制加入各種鄰域變換結(jié)構(gòu)的權(quán)重后,算法會(huì)根據(jù)各鄰域結(jié)構(gòu)的權(quán)重計(jì)算概率選擇下一次迭代時(shí)使用的鄰域,而不是簡單地重復(fù)迭代。各鄰域結(jié)構(gòu)的權(quán)重在搜索時(shí)會(huì)根據(jù)歷史表現(xiàn)和使用次數(shù)實(shí)現(xiàn)更新。在迭代計(jì)算時(shí),算法通過輪盤賭選擇策略生成當(dāng)前解的鄰域結(jié)構(gòu)。假設(shè)鄰域結(jié)構(gòu)集合Q包含3種鄰域結(jié)構(gòu){q1,q2,q3},對應(yīng)的權(quán)重集合ω包含3個(gè)權(quán)重{ω1,ω2,ω3}。h為迭代次數(shù),當(dāng)累計(jì)迭代10次以后且第r種鄰域結(jié)構(gòu)尋得更優(yōu)解時(shí),該鄰域的權(quán)重根據(jù)如下公式進(jìn)行更新。

    ωi=ωi+[h/10]h(16)

    在依概率選擇鄰域時(shí)遵循以下公式:

    Φr=ωr∑3i=1ωi, r=1,2,3(17)

    式中:Φr表示選擇第r種鄰域結(jié)構(gòu)的概率;ωr表示第r種鄰域結(jié)構(gòu)的權(quán)重。目標(biāo)跟蹤列表經(jīng)過AVNS尋優(yōu)后得到的最終方案需要轉(zhuǎn)換為時(shí)間窗輸出。由于使用統(tǒng)一長度為20的列表,只需要統(tǒng)計(jì)列表中每一值段的數(shù)值與長度即可。數(shù)值代表跟蹤火控雷達(dá)的編號(hào),長度為該雷達(dá)在20個(gè)航跡點(diǎn)中的跟蹤時(shí)段,對目標(biāo)列表中的值段進(jìn)行統(tǒng)計(jì)解碼,就可以得到最優(yōu)火控雷達(dá)組網(wǎng)跟蹤方案P′min。由于每部雷達(dá)都被裝備在對應(yīng)的艦船或者預(yù)警機(jī)上,故可以根據(jù)資源狀態(tài)信息I2找到對應(yīng)艦機(jī)編號(hào),并加入P′min。此時(shí)P′min中的值段長度并不能表示時(shí)間,根據(jù)前面取整縮小的操作,對時(shí)間進(jìn)行放大處理,可以得到每個(gè)值段長度代表的時(shí)間長度。該時(shí)間長度加上仿真輸入當(dāng)前時(shí)間和武器控制系統(tǒng)啟動(dòng)準(zhǔn)備時(shí)間,就能得到起始跟蹤時(shí)間和結(jié)束跟蹤時(shí)間。最后輸出到仿真系統(tǒng)的火控雷達(dá)資源調(diào)度來襲目標(biāo)組網(wǎng)跟蹤方案包括目標(biāo)編號(hào)、艦機(jī)編號(hào)、火控雷達(dá)編號(hào)、起始跟蹤時(shí)間和結(jié)束跟蹤時(shí)間。即對Pmin={T:[[J,L,ST,ET],…],…},T表示目標(biāo)編號(hào),J表示艦艇編號(hào);L表示雷達(dá)編號(hào);H表示最大迭代次數(shù);ST表示起始跟蹤時(shí)間;ET表示結(jié)束跟蹤時(shí)間。方案內(nèi)嵌列表數(shù)量大于1,表示有多部雷達(dá)對目標(biāo)進(jìn)行跟蹤。方案字典鍵值對數(shù)量大于1,表示對多個(gè)目標(biāo)進(jìn)行跟蹤。

    基于AVNS的火控雷達(dá)組網(wǎng)資源調(diào)度算法的基本流程如下:

    步驟 1 輸入目標(biāo)態(tài)勢信息I1、資源狀態(tài)信息I2、仿真系統(tǒng)信息I3、來襲目標(biāo)跟蹤任務(wù)規(guī)劃指令I(lǐng)4,并設(shè)定初始參數(shù)。

    步驟 2 篩選I1、I2、I3和I4中存在信息缺失的目標(biāo)及其信息。

    步驟 3 根據(jù)I4,篩除I1、I2和I3中與指令無關(guān)的目標(biāo)及其信息。

    步驟 4 根據(jù)I3中武器攔截范圍信息篩選I1目標(biāo)預(yù)測航跡點(diǎn),得到有效航跡點(diǎn)。

    步驟 5 基于規(guī)則生成初始火控雷達(dá)組網(wǎng)跟蹤方案P,計(jì)算其目標(biāo)函數(shù)值F(P),并設(shè)暫時(shí)最優(yōu)函數(shù)值F′min=F(P),暫時(shí)最優(yōu)方案P′min=P。

    步驟 6 判斷迭代次數(shù)h的終止條件,h小于最大迭代次數(shù)H,執(zhí)行步驟7,否則執(zhí)行步驟12。

    步驟 7 依概率Φ選擇鄰域結(jié)構(gòu)q,根據(jù)鄰域結(jié)構(gòu)q生成新的方案P′,并計(jì)算目標(biāo)函數(shù)值F(P′)。

    步驟 8 判斷當(dāng)前目標(biāo)函數(shù)值F(P′),F(xiàn)(P′)小于暫時(shí)最優(yōu)函數(shù)值F′min,執(zhí)行步驟9,否則執(zhí)行步驟6。

    步驟 9 令暫時(shí)最優(yōu)函數(shù)值F′min=F(P′),暫時(shí)最優(yōu)方案P′min=P′。

    步驟 10 判斷迭代次數(shù)h是否滿足權(quán)重更新條件,如果迭代次數(shù)h≥10,執(zhí)行步驟11,否則執(zhí)行步驟6。

    步驟 11 更新權(quán)重ω和概率Φ。

    步驟 12 根據(jù)當(dāng)前最優(yōu)方案P′min和目標(biāo)態(tài)勢信息I1計(jì)算雷達(dá)組網(wǎng)起始跟蹤時(shí)間ST和結(jié)束跟蹤時(shí)間ET。

    步驟 13 輸出包含目標(biāo)編號(hào)、艦艇編號(hào)、雷達(dá)起始跟蹤時(shí)間和結(jié)束跟蹤時(shí)間的最終方案P′min。

    具體的算法流程如圖4所示。

    3 仿真實(shí)驗(yàn)

    為了驗(yàn)證所提算法的性能,本節(jié)設(shè)置基于AVNS的雷達(dá)組網(wǎng)資源調(diào)度算法、基于VNS的雷達(dá)組網(wǎng)資源調(diào)度算法以及傳統(tǒng)規(guī)則算法的對比試驗(yàn)。

    3.1 實(shí)驗(yàn)場景

    圖5和圖6是兩種模擬的海上艦艇編隊(duì)場景,在每種場景下設(shè)置10種來襲目標(biāo)數(shù)量,分別為5,10,15,20,25,30,35,40,45,50,一共設(shè)置20個(gè)實(shí)驗(yàn)案例。

    第一種海上艦艇編隊(duì)模擬隊(duì)形如下。如圖5所示,以待跟蹤目標(biāo)數(shù)量為20為例,來襲導(dǎo)彈目標(biāo)從3個(gè)方向朝海上艦艇編隊(duì)襲來。其中編號(hào)1~8的目標(biāo)從場景圖左側(cè)出發(fā);編號(hào)9~14的目標(biāo)從場景圖左上方出發(fā);編號(hào)15~20的目標(biāo)從場景圖上方出發(fā)。編號(hào)1~4、9~11、15~17的目標(biāo)為高空目標(biāo),飛行速度為1 000 m/s;其余目標(biāo)為低空目標(biāo),飛行速度為300 m/s。海上艦艇編隊(duì)有2架預(yù)警機(jī)和2艘驅(qū)逐艦裝備火控雷達(dá),對來襲目標(biāo)進(jìn)行跟蹤。其中預(yù)警機(jī)1的火控雷達(dá)具有3個(gè)全時(shí)段可用的通道;預(yù)警機(jī)2的火控雷達(dá)具有5個(gè)0~800 s可用的通道;驅(qū)逐艦1的火控雷達(dá)具有4個(gè)全時(shí)段可用的通道和2個(gè)200 s以后可用的通道;驅(qū)逐艦2的火控雷達(dá)具有2個(gè)全時(shí)段可用的通道和1個(gè)220 s以后可用的通道。預(yù)警機(jī)對面目標(biāo)的最遠(yuǎn)跟蹤探測距離為270 km;艦載火控雷達(dá)對高空目標(biāo)的最遠(yuǎn)跟蹤探測距離為400 km;對低空目標(biāo)的最遠(yuǎn)跟蹤探測距離為30 km。

    第二種海上艦艇編隊(duì)模擬隊(duì)形如下。如圖6所示,以待跟蹤目標(biāo)數(shù)量為20為例,來襲導(dǎo)彈目標(biāo)從3個(gè)方向朝海上艦艇編隊(duì)襲來。其中編號(hào)1~8的目標(biāo)從場景圖左側(cè)出發(fā);編號(hào)9~14的目標(biāo)從場景圖左上方出發(fā);編號(hào)15~20的目標(biāo)從場景圖上方出發(fā)。編號(hào)1~4、9~11、15~17的目標(biāo)為高空目標(biāo),飛行速度為1 000 m/s;其余目標(biāo)為低空目標(biāo),飛行速度為300 m/s。海上艦艇編隊(duì)有2架預(yù)警機(jī)和2艘驅(qū)逐艦裝備火控雷達(dá),對來襲目標(biāo)進(jìn)行跟蹤。其中預(yù)警機(jī)1的火控雷達(dá)具有3個(gè)全時(shí)段可用的通道;預(yù)警機(jī)2的火控雷達(dá)具有5個(gè)0~800 s可用的通道;驅(qū)逐艦1的火控雷達(dá)具有4個(gè)全時(shí)段可用的通道和2個(gè)200 s以后可用的通道;驅(qū)逐艦2的火控雷達(dá)具有2個(gè)全時(shí)段可用的通道和1個(gè)220 s以后可用的通道。預(yù)警機(jī)對面目標(biāo)的最遠(yuǎn)跟蹤探測距離為270 km;艦載火控雷達(dá)對高空目標(biāo)的最遠(yuǎn)跟蹤探測距離為400 km;對低空目標(biāo)的最遠(yuǎn)跟蹤探測距離為30 km。

    3.2 參數(shù)設(shè)置

    所有算法獨(dú)立運(yùn)行30次,迭代搜索次數(shù)為H。

    (1) 所有案例Ei(1≤i≤20,i∈N)的迭代搜索次數(shù)H如表1所示。

    (2) 鄰域數(shù)量Qmax=3;目標(biāo)航跡點(diǎn)總數(shù)M=400。

    (3) 算法參數(shù):γ=1.1,ρ=1.2,ω1=0.33," ω2=0.33,ω3=0.33;低空目標(biāo)威脅值W取8,高空目標(biāo)威脅值W取6。

    (4) 規(guī)則算法參數(shù):γ=1.1,ρ=1.2;低空目標(biāo)威脅值W取8,高空目標(biāo)威脅值W取6。

    3.3 性能指標(biāo)

    選擇目標(biāo)跟蹤總時(shí)長和可攔目標(biāo)跟蹤覆蓋率作為算法性能測試指標(biāo)。

    目標(biāo)跟蹤總時(shí)長為跟蹤方案中所有目標(biāo)的跟蹤時(shí)長之和,其計(jì)算方式如下:

    Ta=∑jj=1tj(18)

    式中:tj是方案中目標(biāo)j的跟蹤時(shí)長。

    可攔目標(biāo)跟蹤覆蓋率的計(jì)算方法如下:

    Tw=∑jj=1t′j(19)

    式中:t′j是目標(biāo)j的待跟蹤時(shí)長;Tw是待跟蹤總時(shí)長。

    B=TaTw(20)

    式中:B是可攔跟蹤覆蓋率。

    3.4 仿真實(shí)驗(yàn)結(jié)果與分析

    3.4.1 模擬海上艦艇編隊(duì)場景1的仿真結(jié)果

    從表2和圖7可以看出,在所進(jìn)行的10組實(shí)驗(yàn)中,AVNS雷達(dá)組網(wǎng)資源調(diào)度算法所得到目標(biāo)雷達(dá)組網(wǎng)資源調(diào)度方案的可攔跟蹤覆蓋率相較于傳統(tǒng)基于規(guī)則的資源調(diào)度算法有明顯的優(yōu)勢,相較于VNS雷達(dá)組網(wǎng)資源調(diào)度算法也有提升。規(guī)則算法根據(jù)目標(biāo)所在的相對位置選擇并調(diào)度雷達(dá)資源,將會(huì)把距離最近的雷達(dá)資源分配給相應(yīng)的目標(biāo),而導(dǎo)致在目標(biāo)預(yù)測航跡后半段的覆蓋能力不足,使得總體可攔跟蹤覆蓋率不高,結(jié)果較差。VNS雷達(dá)組網(wǎng)資源調(diào)度算法通過對目標(biāo)威脅度度和整體雷達(dá)資源的考慮,會(huì)整合雷達(dá)資源對目標(biāo)實(shí)現(xiàn)接續(xù)跟蹤,從而提升目標(biāo)可攔跟蹤范圍覆蓋率,相比規(guī)則算法有明顯提升。所研究的AVNS雷達(dá)組網(wǎng)資源調(diào)度算法,在尋優(yōu)過程中對雷達(dá)組網(wǎng)的資源分配更加合理,使得目標(biāo)跟蹤范圍覆蓋率有進(jìn)一步提升,與規(guī)則算法相比平均提升14.06%,與VNS算法相比平均提升3.43%,有效證明所提出的AVNS雷達(dá)組網(wǎng)資源調(diào)度算法能夠?qū)崿F(xiàn)海上艦艇編隊(duì)作戰(zhàn)場景下的高效雷達(dá)資源調(diào)度。

    圖8~圖10是根據(jù)待跟蹤目標(biāo)數(shù)量為20的實(shí)驗(yàn)數(shù)據(jù)繪制的跟蹤效果甘特圖。為方便展示,數(shù)據(jù)均做取整處理。對于目標(biāo)數(shù)量為18的場景,在規(guī)則算法方案中,由于預(yù)警機(jī)2可用跟蹤時(shí)間有限,在到達(dá)800 s后該目標(biāo)沒有資源接續(xù)跟蹤。而在AVNS雷達(dá)組網(wǎng)資源調(diào)度算法方案中,由驅(qū)逐艦2實(shí)現(xiàn)對目標(biāo)的末端接續(xù)跟蹤,增加對單一目標(biāo)的可攔跟蹤時(shí)長。對于驅(qū)逐艦1,其6個(gè)可用通道數(shù)一開始全部指派完,但在高空目標(biāo)跟蹤完成后,由于規(guī)則算法不能組網(wǎng),通道釋放為空閑狀態(tài)后沒有其他跟蹤任務(wù)需要執(zhí)行,造成資源的浪費(fèi)。而在AVNS雷達(dá)組網(wǎng)資源調(diào)度算法方案中,驅(qū)逐艦1執(zhí)行完高空目標(biāo)的跟蹤任務(wù)后,可以繼續(xù)執(zhí)行對目標(biāo)12和目標(biāo)13兩個(gè)低空目標(biāo)的跟蹤任務(wù),提高了單一資源的利用效率。

    表3是對該次實(shí)驗(yàn)?zāi)繕?biāo)可攔跟蹤時(shí)段和可攔跟蹤覆蓋率的統(tǒng)計(jì)。在目標(biāo)數(shù)量較多、資源不足的情況下,規(guī)則算法對17個(gè)目標(biāo)實(shí)現(xiàn)跟蹤;VNS雷達(dá)組網(wǎng)資源調(diào)度算法對18個(gè)目標(biāo)實(shí)行跟蹤,其中對7個(gè)目標(biāo)接續(xù)跟蹤;AVNS雷達(dá)組網(wǎng)資源調(diào)度算法對18個(gè)目標(biāo)實(shí)行跟蹤,其中對8個(gè)目標(biāo)接續(xù)跟蹤。AVNS雷達(dá)組網(wǎng)資源調(diào)度算法相較于規(guī)則算法,在8個(gè)目標(biāo)的跟蹤時(shí)長上有提升,相較于VNS雷達(dá)組網(wǎng)資源調(diào)度算法在5個(gè)目標(biāo)的跟蹤時(shí)長上有提升。AVNS雷達(dá)組網(wǎng)資源調(diào)度算法的總體可攔跟蹤時(shí)長覆蓋率達(dá)到81.85%,相較于規(guī)則算法的71.87%,提升了9.98%,相較于VNS雷達(dá)組網(wǎng)資源調(diào)度算法的78.18%,提升了3.67%。

    圖11是對不同目標(biāo)數(shù)量下算法平均運(yùn)行時(shí)長的統(tǒng)計(jì)。由圖11可知3種算法的運(yùn)行時(shí)長均與待跟蹤目標(biāo)數(shù)量呈正相關(guān),均在16 ms至37 ms內(nèi)。兩種智能優(yōu)化算法使用規(guī)則算法計(jì)算初始解,因此運(yùn)行時(shí)長都大于規(guī)則算法,但相差都在8 ms以內(nèi)。且AVNS雷達(dá)組網(wǎng)資源調(diào)度算法的運(yùn)行時(shí)長整體小于VNS雷達(dá)組網(wǎng)資源調(diào)度算法。因此,可以認(rèn)為AVNS雷達(dá)組網(wǎng)資源調(diào)度算法能夠滿足防空作戰(zhàn)對火控雷達(dá)資源調(diào)度時(shí)效性的要求。

    3.4.2 模擬海上艦艇編隊(duì)場景2的仿真結(jié)果

    從表4和圖12可以看出,在所進(jìn)行的10組實(shí)驗(yàn)中,AVNS 雷達(dá)組網(wǎng)資源調(diào)度算法所得到目標(biāo)雷達(dá)組網(wǎng)資源調(diào)度方案的可攔跟蹤覆蓋率相較于規(guī)則算法仍然有明顯的優(yōu)勢,并且仍優(yōu)于VNS雷達(dá)組網(wǎng)資源調(diào)度算法。由于場景2的海上艦艇編隊(duì)均勻排布在來襲目標(biāo)預(yù)測航跡兩旁,更有利于雷達(dá)資源對目標(biāo)的跟蹤。相比與場景1,無論是規(guī)則算法、VNS雷達(dá)組網(wǎng)資源調(diào)度算法或是AVNS 雷達(dá)組網(wǎng)資源調(diào)度算法,場景2的整體可攔跟蹤覆蓋率在不同待跟蹤目標(biāo)數(shù)量下均有提升。同時(shí)與規(guī)則算法相比,AVNS雷達(dá)組網(wǎng)資源調(diào)度算法的目標(biāo)可攔跟蹤覆蓋率平均提升13.18%,與VNS雷達(dá)組網(wǎng)資源調(diào)度算法相比平均提升2.63%,有效證明所提出的AVNS雷達(dá)組網(wǎng)資源調(diào)度算法能夠?qū)崿F(xiàn)海上艦艇編隊(duì)作戰(zhàn)場景下的高效雷達(dá)資源調(diào)度。

    由圖13可知3種算法的運(yùn)行時(shí)長均與待跟蹤目標(biāo)數(shù)量呈正相關(guān),均在15 ms至41 ms內(nèi)。兩種智能優(yōu)化算法的運(yùn)行時(shí)長與規(guī)則算法均相差10 ms以內(nèi),且AVNS雷達(dá)組網(wǎng)資源調(diào)度算法運(yùn)行時(shí)長整體小于VNS雷達(dá)組網(wǎng)資源調(diào)度算法。因此,可以認(rèn)為AVNS雷達(dá)組網(wǎng)資源調(diào)度算法能夠滿足防空作戰(zhàn)對火控雷達(dá)資源調(diào)度時(shí)效性的要求。

    4 結(jié) 論

    本文以海上編隊(duì)防空反導(dǎo)作戰(zhàn)任務(wù)為背景,針對來襲目標(biāo)跟蹤問題,提出基于AVNS的雷達(dá)組網(wǎng)資源調(diào)度方法。首先,對初始目標(biāo)航跡點(diǎn)進(jìn)行篩選,以最大有效航跡覆蓋率、最小目標(biāo)距離和最大目標(biāo)威脅度為目標(biāo)函數(shù),對雷達(dá)資源探測范圍、有效航跡最小跟蹤覆蓋距離和雷達(dá)剩余通道數(shù)量進(jìn)行量化,將資源調(diào)度問題抽象為3約束、3目標(biāo)的多目標(biāo)優(yōu)化問題。設(shè)計(jì)3種鄰域變換算子并加入自適應(yīng)迭代策略,有效降低目標(biāo)函數(shù)值,最后通過時(shí)間窗值段統(tǒng)計(jì)得出雷達(dá)組網(wǎng)目標(biāo)跟蹤方案。

    為驗(yàn)證所提基于AVNS的雷達(dá)組網(wǎng)資源調(diào)度算法的有效性,在兩種模擬海上艦艇編隊(duì)場景的10種來襲目標(biāo)跟蹤需求中設(shè)置對比仿真實(shí)驗(yàn)。實(shí)驗(yàn)結(jié)果表明,在大多數(shù)場景中,所提算法均能得到目標(biāo)函數(shù)值更優(yōu)的雷達(dá)跟蹤方案,且能滿足防空作戰(zhàn)對火控雷達(dá)資源調(diào)度時(shí)效性的要求。

    參考文獻(xiàn)

    [1]XIAO B S. The development review of airborne fire control radar technology[C]∥Proc.of the China Automation Congress, 2021: 4502-4507.

    [2]JAMES J. A maritime conversation with America[J]. Orbis, 2022, 66(2): 166-183.

    [3]LIN D. Software and hardware implementation of signal processing for shipborne fire control radar[D]. Xi’an: Xidian University,2019.

    [4]LIEF L. Design considerations for a modern naval fire control radar[C]∥Proc.of the IEEE Radar Conference, 2010: 615-619.

    [5]DAI X R, SHI C G, WANG Z W, et al. Coalition game theore tic power allocation strategy for multi target detection in distribu ted radar networks[J]. Remote Sensing, 2023, 15(15): 3804.

    [6]Assistant Secretary of the Navy for Research. Updated responsibilities for management of naval integrated fire control counter air (NIFC CA)[R]. Washington, DC: Joint Memorandum, 2002.

    [7]YAN J K. Resource allocation for search and track application in phased array radar based on Pareto bi objective optimization[J]. IEEE Trans.on Vehicular Technology, 2019, 68(4): 3487-3499.

    [8]LIU D. Adaptive scheduling algorithm based on CPI and impact of tasks for multifunction radar[J]. IEEE Sensors Journal, 2019, 19(23): 11205-11212.

    [9]王巖松, 韓星, 陳春, 等. 綜合優(yōu)先級(jí)下相控陣火控雷達(dá)自適應(yīng)調(diào)度算法[J]. 火控雷達(dá)技術(shù), 2023, 52(1): 100-109.

    WANG Y S, HAN X, CHEN C, et al. A comprehensive priority based adaptive scheduling algorithm for phased array fire control radar[J]. Fire Control Radar Technology, 2023, 52(1): 100-109.

    [10]周克強(qiáng), 王海川. 基于相控陣?yán)走_(dá)的艦炮聯(lián)合反導(dǎo)技術(shù)[J]. 火力與指揮控制, 2019, 44(6): 159-163.

    ZHOU K Q, WANG H C. The anti missile research of combined naval gun system based on phased array radar[J]. Fire Control amp; Command Control, 2019, 44(6): 159-163.

    [11]閆明松, 魏旭鴻. 海上編隊(duì)艦載雷達(dá)協(xié)同探測技術(shù)研究[J]. 艦船電子對抗, 2021, 44(4): 7-11, 120.

    YAN M S, WEI X H. Research into the cooperative detection technology of shipborne radar for sea formation[J]. Shipboard Electronic Countermeasure, 2021, 44(4): 7-11, 120.

    [12]JIANG W, QI Z X, YE Z, et al. Research on cooperative detection technology of networked radar based on data fusion[C]∥Proc.of the 2nd China International SAR Symposium, 2021.

    [13]ZHOU S H, ZENG D X, LI H B, et al. Maneuvering target tracking in multi radar cooperative detection data fusion[C]∥Proc.of the 7th International Symposium on Advances in Electrical, Electronics, and Computer Engineering, 2022.

    [14]LI Z Q. Design and reliability evaluation of simulation system for fire control radar network[C]∥Proc.of the International Conference on Quality, Reliability, Risk, Maintenance, and Safety Engineering, 2012: 1314-1317.

    [15]GARCIA N, HAIMOVICH M A, COULON M, et al. Resource allocation in MIMO radar with multiple targets for non coherent localization[J]. IEEE Trans.on Signal Processing, 2014, 62(10): 2656-2666.

    [16]WANG K L. Research and implementation of shipboard phased array radar target tracking algorithm[D]. Xi’an: Xidian university, 2021.

    [17]LUO H, SHANG C X, HAN Z Z. HLA based anti reconnaissance fire control radar network data fusion and control simulation system development[J]. Advanced Materials Research, 2014, 3181(926-930): 2486-2489.

    [18]SHI C G, WANG Y J, SANA S. Joint transmit resource mana gement and waveform selection strategy for target tracking in distributed phased array radar network[J]. IEEE Trans.on Aerospace and Electronic Systems, 2022, 58(4): 2762-2778.

    [19]LI T T, MENG T K, MENG G L, et al. Formation optimization of airborne radar coordinated detection system using an improved artificial fish swarm algorithm[J]. Scientific Reports, 2024, 14(1): 248.

    [20]GRANT C J. Air and missile defense: transformations for 21st century warfighting[J]. Johns Hopkins APL Technical Digest, 2020, 35(2): 83-89.

    [21]YING Z. Research on dynamic anti jamming strategy of shipboard fire control radar[C]∥Proc.of the International Conference on Computer, Artificial Intelligence, and Control Engineering, 2022.

    [22]NICHOLAS A O. Distributed kill chain drawing insight for mosaic from the immune system and navy[R]. Santa Monica: Library of Congress Cataloging in Publication, 2021.

    [23]張財(cái)生, 司沈闖. 美軍預(yù)警機(jī)與艦載雷達(dá)協(xié)同預(yù)警探測[J]. 指揮控制與仿真, 2020, 42(3): 39-41.

    ZHANG C S, SI S C. Research on cooperative early warning detection of early warning aircraft and shipborne radar[J]. Command Control amp; Simulation, 2020, 42(3): 39-41.

    [24]CHENG Z L. Adaptive resource optimization management of phased array radar[D]. Chengdu: University of Electronic Science and Technology of China, 2022.

    [25]MOO P W, DING Z. Coordinated radar resource management for networked phased array radars[J]. IET Radar, Sonar amp; Navigation, 2015, 9(8): 1009-1020.

    [26]TUNCER O, CIRPAN H A. A networked radar resource mana gement approach utilizing target priority and maneuver for cooperative air defense fire control radars[J]. IEEE Access, 2023, 11(1): 136279-136291.

    [27]TANG J C. Research on resource scheduling method of multifunctional radar network[D]. Chengdu: University of Electronic Science and Technology of China, 2021.

    [28]宋曉程, 李陟, 任海偉, 等. 目標(biāo)動(dòng)態(tài)威脅度驅(qū)動(dòng)的分布式組網(wǎng)相控陣?yán)走_(dá)資源優(yōu)化分配算法[J]. 雷達(dá)學(xué)報(bào), 2023, 12(3): 629-641.

    SONG X C, LI Z, REN H W, et al. Threat driven resource allocation algorithm for distributed netted phased array radars[J]. Journal of Radars, 2023, 12(3): 629-641.

    [29]SHI C, JING D, SALOUS S, et al. Collaborative trajectory planning and resource allocation for multi target tracking in airborne radar networks under spectral coexistence[J]. Remote Sensing, 2023, 15(13): 3386.

    [30]WANG F J, DING X Y, KONG M. Application of multi sensor track fusion technology in fire control radar network[C]∥Proc.of the IEEE 11th International Conference on Electronic Measurement amp; Instruments, 2013: 563-567.

    [31]HE J, ZHU S R, SUN X Q. Research on flight scheduling of UAV in cooperation with passive radar for detection based on improved PSO algorithm[C]∥Proc.of the 4th International Conference on Computer Science and Communication Technology,2023.

    作者簡介

    汪達(dá)旺(1996—),男,碩士,主要研究方向?yàn)橹悄軟Q策任務(wù)規(guī)劃與資源調(diào)度。

    陸志灃(1981—),男,博士研究生,主要研究方向?yàn)槿斯ぶ悄芘c建模仿真。

    伍國華(1986—),男,教授,博士,主要研究方向?yàn)橛?jì)算智能、任務(wù)規(guī)劃與資源調(diào)度。

    一级爰片在线观看| 午夜福利网站1000一区二区三区| 最新中文字幕久久久久| 国产视频首页在线观看| 亚洲欧美中文字幕日韩二区| 直男gayav资源| 中文资源天堂在线| 亚洲精品aⅴ在线观看| 久久久成人免费电影| 综合色av麻豆| 色尼玛亚洲综合影院| 日韩一区二区视频免费看| 亚洲国产高清在线一区二区三| 美女大奶头视频| 激情五月婷婷亚洲| 国产麻豆成人av免费视频| 亚洲av成人精品一区久久| 亚洲色图av天堂| 91av网一区二区| 深爱激情五月婷婷| 国产精品伦人一区二区| 人妻制服诱惑在线中文字幕| 精品国产三级普通话版| 日韩av在线大香蕉| 免费看日本二区| 91精品一卡2卡3卡4卡| 色网站视频免费| 男人舔女人下体高潮全视频| 亚洲高清免费不卡视频| 日日啪夜夜爽| 肉色欧美久久久久久久蜜桃 | or卡值多少钱| 少妇熟女欧美另类| 午夜免费观看性视频| 熟女电影av网| 97超视频在线观看视频| 日韩亚洲欧美综合| 日韩一区二区视频免费看| 亚洲欧美一区二区三区黑人 | 麻豆乱淫一区二区| 波多野结衣巨乳人妻| 日韩欧美 国产精品| 亚洲综合色惰| 一本久久精品| 男女那种视频在线观看| 成人毛片60女人毛片免费| 国产成人91sexporn| 国产精品久久久久久av不卡| 99热这里只有精品一区| 三级经典国产精品| 亚洲av男天堂| 五月天丁香电影| 亚洲av成人精品一区久久| 欧美日韩一区二区视频在线观看视频在线 | 大又大粗又爽又黄少妇毛片口| 三级经典国产精品| 男的添女的下面高潮视频| 熟女电影av网| 亚洲精品日韩在线中文字幕| 婷婷色综合大香蕉| 日韩制服骚丝袜av| 国产 亚洲一区二区三区 | 99久久精品热视频| 国产精品人妻久久久影院| 欧美性感艳星| 亚洲精品影视一区二区三区av| 成人无遮挡网站| 美女cb高潮喷水在线观看| 久久精品夜色国产| 亚洲av电影不卡..在线观看| 黄色欧美视频在线观看| 99re6热这里在线精品视频| 久久久欧美国产精品| 欧美+日韩+精品| 亚洲性久久影院| 大陆偷拍与自拍| 午夜福利高清视频| 日韩欧美精品v在线| 亚洲av成人精品一二三区| 久久人人爽人人爽人人片va| 一级二级三级毛片免费看| av一本久久久久| 卡戴珊不雅视频在线播放| 国产亚洲av嫩草精品影院| 最近的中文字幕免费完整| 免费看日本二区| 精品酒店卫生间| 久久精品夜夜夜夜夜久久蜜豆| 色播亚洲综合网| 欧美高清性xxxxhd video| 国产黄片美女视频| 床上黄色一级片| 熟女人妻精品中文字幕| av专区在线播放| 亚洲激情五月婷婷啪啪| 国产精品av视频在线免费观看| 日本猛色少妇xxxxx猛交久久| 超碰av人人做人人爽久久| 午夜免费激情av| 久久99热6这里只有精品| 99久久人妻综合| 少妇的逼水好多| 国产精品1区2区在线观看.| 日产精品乱码卡一卡2卡三| 搡女人真爽免费视频火全软件| 日韩av免费高清视频| 一区二区三区四区激情视频| 五月天丁香电影| 中文字幕制服av| 国产午夜精品一二区理论片| 高清午夜精品一区二区三区| 国产视频首页在线观看| 国产91av在线免费观看| 身体一侧抽搐| 久久精品久久精品一区二区三区| 在线观看人妻少妇| 神马国产精品三级电影在线观看| 日本av手机在线免费观看| 欧美成人精品欧美一级黄| 三级经典国产精品| 久久久久久久午夜电影| 亚洲美女搞黄在线观看| 一区二区三区高清视频在线| 日本爱情动作片www.在线观看| 亚洲av电影不卡..在线观看| 国产老妇女一区| 国产精品一区二区三区四区久久| 国产精品伦人一区二区| 寂寞人妻少妇视频99o| 丰满乱子伦码专区| 精品久久久噜噜| 国产乱人视频| 亚洲aⅴ乱码一区二区在线播放| 久久久久久久久久久丰满| 你懂的网址亚洲精品在线观看| 午夜激情久久久久久久| 美女高潮的动态| 大香蕉久久网| 日日撸夜夜添| 国产精品国产三级国产专区5o| 久久鲁丝午夜福利片| 亚洲国产av新网站| 少妇人妻精品综合一区二区| 女人久久www免费人成看片| av免费在线看不卡| 乱系列少妇在线播放| 亚洲精品一二三| 亚洲自拍偷在线| 毛片一级片免费看久久久久| 99久久精品一区二区三区| 国产精品无大码| 天天躁日日操中文字幕| 99久久精品一区二区三区| 自拍偷自拍亚洲精品老妇| 久久97久久精品| 国产在视频线在精品| 我的老师免费观看完整版| 综合色av麻豆| 久久精品国产自在天天线| 亚洲aⅴ乱码一区二区在线播放| 亚洲国产精品成人综合色| 麻豆乱淫一区二区| 欧美一级a爱片免费观看看| 十八禁网站网址无遮挡 | av免费观看日本| 舔av片在线| 在现免费观看毛片| 日本三级黄在线观看| 久久久久久久久久成人| 永久免费av网站大全| 99久久精品热视频| 亚洲国产精品成人综合色| 久久久久性生活片| 久久久久久久亚洲中文字幕| 人体艺术视频欧美日本| 九九久久精品国产亚洲av麻豆| 日本一二三区视频观看| 午夜精品一区二区三区免费看| 在线天堂最新版资源| 男女国产视频网站| 伊人久久国产一区二区| 亚洲欧美精品专区久久| 午夜福利视频1000在线观看| 国产探花极品一区二区| 久久精品综合一区二区三区| 国产在线一区二区三区精| 国产女主播在线喷水免费视频网站 | 国产精品蜜桃在线观看| 中文字幕av成人在线电影| 日日摸夜夜添夜夜爱| 欧美一级a爱片免费观看看| 亚洲精品亚洲一区二区| 亚洲综合精品二区| 看黄色毛片网站| 亚洲av成人精品一二三区| 成人性生交大片免费视频hd| 精品不卡国产一区二区三区| 观看美女的网站| 精品人妻一区二区三区麻豆| 深夜a级毛片| 国产亚洲av嫩草精品影院| 尾随美女入室| 国产成人精品婷婷| 欧美 日韩 精品 国产| 日本黄大片高清| 青青草视频在线视频观看| 国产大屁股一区二区在线视频| 国产极品天堂在线| 毛片女人毛片| 欧美不卡视频在线免费观看| 最近视频中文字幕2019在线8| 亚洲美女搞黄在线观看| 日韩欧美精品v在线| 国产亚洲精品久久久com| 欧美+日韩+精品| 亚州av有码| 成人二区视频| 亚洲av国产av综合av卡| 春色校园在线视频观看| 女人久久www免费人成看片| 有码 亚洲区| 亚洲av不卡在线观看| 麻豆久久精品国产亚洲av| 久久精品久久精品一区二区三区| 免费av观看视频| www.色视频.com| 国产精品一二三区在线看| 人体艺术视频欧美日本| 女人十人毛片免费观看3o分钟| 联通29元200g的流量卡| 色综合亚洲欧美另类图片| 在线免费观看不下载黄p国产| 三级经典国产精品| 搡老妇女老女人老熟妇| 中文字幕久久专区| 男人舔女人下体高潮全视频| 女人被狂操c到高潮| 精品一区二区三区人妻视频| 欧美xxxx黑人xx丫x性爽| 日韩成人伦理影院| 爱豆传媒免费全集在线观看| 欧美高清性xxxxhd video| 亚洲精品久久久久久婷婷小说| 一区二区三区免费毛片| 亚洲最大成人中文| 舔av片在线| 亚洲av成人精品一区久久| 日韩中字成人| 婷婷色av中文字幕| 蜜臀久久99精品久久宅男| 男插女下体视频免费在线播放| 久久久久久久久久成人| 一级毛片电影观看| 国产精品麻豆人妻色哟哟久久 | 国产一区有黄有色的免费视频 | 在线免费观看不下载黄p国产| 精品不卡国产一区二区三区| 毛片一级片免费看久久久久| 久久鲁丝午夜福利片| 久久久久久国产a免费观看| 高清在线视频一区二区三区| 亚洲欧洲国产日韩| 如何舔出高潮| 国产 一区精品| 熟女人妻精品中文字幕| 亚洲成人av在线免费| 欧美日韩视频高清一区二区三区二| 国产成人a区在线观看| 看非洲黑人一级黄片| 熟女人妻精品中文字幕| 国产精品熟女久久久久浪| 国产伦在线观看视频一区| 免费av毛片视频| a级毛色黄片| 精品一区二区三区人妻视频| 18禁在线无遮挡免费观看视频| 国产精品不卡视频一区二区| 一级片'在线观看视频| 亚洲图色成人| 亚洲欧美日韩无卡精品| 少妇熟女欧美另类| 国产午夜福利久久久久久| 国精品久久久久久国模美| 少妇的逼水好多| av又黄又爽大尺度在线免费看| 亚洲久久久久久中文字幕| 久久99精品国语久久久| 午夜激情欧美在线| 人妻夜夜爽99麻豆av| 欧美xxxx黑人xx丫x性爽| 十八禁网站网址无遮挡 | 亚洲经典国产精华液单| 亚洲国产精品成人综合色| 久久久久免费精品人妻一区二区| 丰满乱子伦码专区| 成人午夜精彩视频在线观看| 又粗又硬又长又爽又黄的视频| av线在线观看网站| 国产不卡一卡二| 美女高潮的动态| 日韩一本色道免费dvd| 99久久精品国产国产毛片| 国产精品.久久久| 一级二级三级毛片免费看| 国产伦精品一区二区三区四那| 丰满人妻一区二区三区视频av| 18禁动态无遮挡网站| 免费看a级黄色片| 熟女人妻精品中文字幕| 男的添女的下面高潮视频| 99久久精品国产国产毛片| 欧美最新免费一区二区三区| 久久午夜福利片| 天天一区二区日本电影三级| 久久精品久久精品一区二区三区| 伊人久久精品亚洲午夜| 日日干狠狠操夜夜爽| 亚洲欧美成人精品一区二区| 狂野欧美激情性xxxx在线观看| 精品国内亚洲2022精品成人| av又黄又爽大尺度在线免费看| 高清在线视频一区二区三区| 99热网站在线观看| 中文在线观看免费www的网站| 少妇丰满av| 亚洲精品456在线播放app| 亚洲精品视频女| 在线免费十八禁| 成年女人看的毛片在线观看| 在现免费观看毛片| 午夜福利网站1000一区二区三区| 中文乱码字字幕精品一区二区三区 | 在线观看美女被高潮喷水网站| 欧美潮喷喷水| 久久久久久国产a免费观看| 青春草视频在线免费观看| 亚洲av免费在线观看| 国产在线一区二区三区精| 26uuu在线亚洲综合色| 久热久热在线精品观看| 国产一区二区在线观看日韩| 亚洲人成网站在线播| av免费观看日本| 超碰97精品在线观看| 亚洲真实伦在线观看| 99热6这里只有精品| 成人毛片a级毛片在线播放| 一本一本综合久久| 国产高潮美女av| 亚洲精品日本国产第一区| 欧美xxⅹ黑人| 精品欧美国产一区二区三| 午夜福利在线观看免费完整高清在| 亚洲婷婷狠狠爱综合网| 最近中文字幕2019免费版| 91精品伊人久久大香线蕉| 色综合色国产| 亚洲高清免费不卡视频| 在线天堂最新版资源| 插阴视频在线观看视频| 免费观看av网站的网址| 日本一二三区视频观看| 国产免费一级a男人的天堂| 毛片女人毛片| 成人毛片a级毛片在线播放| 久久久久网色| ponron亚洲| 国产人妻一区二区三区在| 寂寞人妻少妇视频99o| 日本欧美国产在线视频| 日韩欧美 国产精品| av女优亚洲男人天堂| 麻豆久久精品国产亚洲av| 日韩亚洲欧美综合| 男女啪啪激烈高潮av片| 久久人人爽人人片av| 男人舔奶头视频| 免费黄色在线免费观看| 街头女战士在线观看网站| 色视频www国产| 亚洲成人精品中文字幕电影| 特大巨黑吊av在线直播| 日韩制服骚丝袜av| 国产综合懂色| 777米奇影视久久| 青春草视频在线免费观看| 2018国产大陆天天弄谢| 国产成人午夜福利电影在线观看| 亚洲av一区综合| 99re6热这里在线精品视频| 日韩成人av中文字幕在线观看| 日韩国内少妇激情av| 一区二区三区四区激情视频| 人人妻人人澡人人爽人人夜夜 | 一个人看视频在线观看www免费| 日韩成人av中文字幕在线观看| 天天一区二区日本电影三级| 大片免费播放器 马上看| 在线a可以看的网站| 国产在视频线精品| 国产av码专区亚洲av| 99久久精品国产国产毛片| 人妻夜夜爽99麻豆av| 国产精品一区www在线观看| 中文天堂在线官网| 99热6这里只有精品| 国产黄频视频在线观看| 女的被弄到高潮叫床怎么办| 99热网站在线观看| 精品国内亚洲2022精品成人| 嫩草影院精品99| 高清视频免费观看一区二区 | 国产成人精品久久久久久| 五月天丁香电影| 精品人妻偷拍中文字幕| 国产精品一二三区在线看| 国产真实伦视频高清在线观看| 久久久久久久大尺度免费视频| 亚洲欧美成人综合另类久久久| 高清在线视频一区二区三区| 2018国产大陆天天弄谢| 国产伦理片在线播放av一区| 91精品伊人久久大香线蕉| 99视频精品全部免费 在线| 最近手机中文字幕大全| 国产精品综合久久久久久久免费| 麻豆乱淫一区二区| 26uuu在线亚洲综合色| 尤物成人国产欧美一区二区三区| 国产国拍精品亚洲av在线观看| 波多野结衣巨乳人妻| 国产精品女同一区二区软件| 九草在线视频观看| 男插女下体视频免费在线播放| 成人毛片60女人毛片免费| 午夜激情福利司机影院| 亚洲精品成人久久久久久| 亚洲图色成人| 一边亲一边摸免费视频| 亚洲国产日韩欧美精品在线观看| 91在线精品国自产拍蜜月| 国产精品一区二区在线观看99 | 成人亚洲精品一区在线观看 | av又黄又爽大尺度在线免费看| 可以在线观看毛片的网站| 亚洲最大成人av| 亚洲最大成人中文| 亚洲综合精品二区| 男女边摸边吃奶| 久久精品久久久久久久性| 精品久久久久久久久亚洲| 国产在线一区二区三区精| 亚洲熟妇中文字幕五十中出| 中文资源天堂在线| 街头女战士在线观看网站| 免费看av在线观看网站| 搡女人真爽免费视频火全软件| 国产色爽女视频免费观看| xxx大片免费视频| 可以在线观看毛片的网站| 18禁在线播放成人免费| 欧美最新免费一区二区三区| 简卡轻食公司| 亚洲成人精品中文字幕电影| 国产乱人偷精品视频| 久久久国产一区二区| 久久亚洲国产成人精品v| 91在线精品国自产拍蜜月| 亚洲成人av在线免费| av女优亚洲男人天堂| 国产黄频视频在线观看| 精品亚洲乱码少妇综合久久| 亚洲国产欧美人成| 免费看光身美女| 性色avwww在线观看| 国产女主播在线喷水免费视频网站 | 欧美三级亚洲精品| 国产中年淑女户外野战色| 婷婷色综合大香蕉| 成人美女网站在线观看视频| 舔av片在线| 国产精品.久久久| 午夜精品在线福利| 亚洲成人精品中文字幕电影| 国产综合懂色| 嫩草影院精品99| 看十八女毛片水多多多| 日韩亚洲欧美综合| 少妇人妻一区二区三区视频| 真实男女啪啪啪动态图| 搡老乐熟女国产| 色综合色国产| 国产老妇伦熟女老妇高清| 中文欧美无线码| 尤物成人国产欧美一区二区三区| 一级毛片电影观看| 国产精品.久久久| 久久久久网色| 五月伊人婷婷丁香| 成人综合一区亚洲| 69人妻影院| 永久网站在线| a级毛色黄片| 99久久精品热视频| 少妇的逼水好多| 亚洲国产精品成人综合色| 国产精品一区二区性色av| 成人鲁丝片一二三区免费| 一级毛片 在线播放| 国产成人午夜福利电影在线观看| 大片免费播放器 马上看| 男的添女的下面高潮视频| 久久草成人影院| 日韩三级伦理在线观看| 欧美+日韩+精品| 男插女下体视频免费在线播放| 汤姆久久久久久久影院中文字幕 | 亚洲乱码一区二区免费版| 精品久久久久久久久久久久久| 有码 亚洲区| 大香蕉久久网| 能在线免费观看的黄片| 久久精品夜夜夜夜夜久久蜜豆| 一级毛片aaaaaa免费看小| 成人午夜高清在线视频| 午夜福利网站1000一区二区三区| 国产精品蜜桃在线观看| 2018国产大陆天天弄谢| 精品酒店卫生间| 国产人妻一区二区三区在| 久久精品国产自在天天线| 联通29元200g的流量卡| 69av精品久久久久久| 91久久精品电影网| 免费黄网站久久成人精品| 亚洲国产欧美人成| 亚洲成色77777| 在线免费观看不下载黄p国产| 99久久精品热视频| 亚洲欧美成人精品一区二区| 亚洲va在线va天堂va国产| 国产三级在线视频| 久久久久久久亚洲中文字幕| 日本av手机在线免费观看| 大话2 男鬼变身卡| 午夜福利成人在线免费观看| 亚洲欧洲日产国产| 欧美最新免费一区二区三区| 国产精品国产三级国产专区5o| 国产成人免费观看mmmm| 欧美日本视频| 久久国内精品自在自线图片| 日本午夜av视频| 国产免费又黄又爽又色| 69av精品久久久久久| 国产乱人视频| 又爽又黄a免费视频| 美女xxoo啪啪120秒动态图| 插逼视频在线观看| 色视频www国产| 亚洲自偷自拍三级| 亚洲最大成人中文| 能在线免费观看的黄片| 久久99热这里只有精品18| 91精品一卡2卡3卡4卡| 最近手机中文字幕大全| 91av网一区二区| 可以在线观看毛片的网站| 嫩草影院入口| 国产精品国产三级专区第一集| 丝袜喷水一区| 国产高清有码在线观看视频| 中文字幕av成人在线电影| 又爽又黄a免费视频| 青青草视频在线视频观看| av在线蜜桃| 狠狠精品人妻久久久久久综合| a级毛片免费高清观看在线播放| 亚洲av二区三区四区| 在线免费观看不下载黄p国产| 日日撸夜夜添| 天天一区二区日本电影三级| 亚洲成人久久爱视频| 国产伦精品一区二区三区视频9| 91精品一卡2卡3卡4卡| 大香蕉97超碰在线| 男人舔女人下体高潮全视频| 免费av观看视频| 99久国产av精品国产电影| 黄片wwwwww| 成人午夜精彩视频在线观看| 亚洲综合精品二区| 美女主播在线视频| 国产成人午夜福利电影在线观看| 午夜福利视频1000在线观看| 少妇被粗大猛烈的视频| 国产精品人妻久久久久久| 国产亚洲av嫩草精品影院| 少妇被粗大猛烈的视频| 亚洲精品一二三| 在现免费观看毛片| 亚洲美女视频黄频| 国产伦在线观看视频一区| 亚洲欧美中文字幕日韩二区| 免费av毛片视频| 日本熟妇午夜| 51国产日韩欧美| 国产精品国产三级国产av玫瑰| 岛国毛片在线播放| 三级毛片av免费| 久久久久久久久久黄片| 久久鲁丝午夜福利片| 韩国av在线不卡| 国产国拍精品亚洲av在线观看| 久久久久久久久久人人人人人人|