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

    基于NSGA-II的車載光學(xué)測量設(shè)備任務(wù)調(diào)度方案優(yōu)化

    2024-11-21 00:00:00許強強柴華
    關(guān)鍵詞:優(yōu)化設(shè)備

    摘" 要:

    針對車載光學(xué)測量設(shè)備任務(wù)調(diào)度方案優(yōu)化問題,提出了一種基于非支配排序的遺傳算法(non-dominated sorting genetic algorithm II, NSGA-II)的多目標(biāo)遺傳算法。首先,建立了包含約束、優(yōu)化指標(biāo)在內(nèi)的觀測任務(wù)調(diào)度問題的數(shù)學(xué)模型。其中,針對多優(yōu)化指標(biāo)進行巧妙處理,將某些不作為最優(yōu)指標(biāo)的優(yōu)化指標(biāo)作為指標(biāo)約束進行處理。其次,基于NSGA-II中的快速非優(yōu)超排序方法計算多目標(biāo)適應(yīng)度函數(shù)與選擇算子,多目標(biāo)優(yōu)化求解得到的Pareto最優(yōu)解集即為任務(wù)調(diào)度方案集。最后,通過仿真算例對所提算法進行了求解驗證。仿真結(jié)果表明,該算法能夠有效解決任務(wù)調(diào)度方案優(yōu)化問題,為車載光學(xué)測量設(shè)備的工程實踐提供了一定的參考。

    關(guān)鍵詞:

    車載光學(xué)測量設(shè)備; 任務(wù)調(diào)度; 多目標(biāo)優(yōu)化; 遺傳算法

    中圖分類號:

    V 11

    文獻標(biāo)志碼: A""" DOI:10.12305/j.issn.1001-506X.2024.07.21

    Optimization of task dispatch plan for vehicular optical observation

    equipment based on NSGA-II

    XU Qiangqiang*, CHAI Hua

    (Space Engineering University, Beijing 101416, China)

    Abstract:

    To improve the task dispatch plan for vehicular optical observation equipment, a multi-objective genetic algorithm based on non-dominated sorting genetic algorithm II (NSGA-II) is proposed. Firstly, the task dispatch problem is modeled with the constraints and optimum indexes. To deal with multiple optimum indexes, the optimum index which is not considered in the objective function is considered as a constraint. Secondly, the multi-objective fitness function and selection operator is calculated based on the fast non-dominated sorting method of NSGA-II. The Pareto solution set obtained by the multi-objective optimization is the task dispatch plan solution. Finally, the proposed algorithm is verified by a simulation example. The simulation results show that this method can solve the task dispatch plan problem effectively, which is valuable for the application of the vehicular optical observation equipment.

    Keywords:

    vehicular optical observation equipment; task dispatch; multi-objective optimization; genetic algorithm

    0" 引" 言

    空間目標(biāo)觀測任務(wù)具有隨機性強、約束條件多等突出特點[1],為了滿足不同任務(wù)需求,目前已經(jīng)形成了天基觀測、地基固定觀測、車載機動觀測、船載移動觀測[2]等觀測平臺。相比于天基觀測平臺的資源稀缺性,傳統(tǒng)的地基固定觀測平臺的觀測空間局限性,基于車載光學(xué)測量設(shè)備[34]的機動觀測平臺有效提高了目標(biāo)觀測的靈活性。以此為背景,為高效利用設(shè)備,文獻[5]以車載光學(xué)測量設(shè)備為研究對象,對其所需執(zhí)行的任務(wù)調(diào)度問題開展了研究,重點給出了相關(guān)數(shù)學(xué)模型。但是,在模型的求解方面,還有待進一步改進和優(yōu)化。本文在此基礎(chǔ)上,針對車載光學(xué)測量設(shè)備任務(wù)調(diào)度方案優(yōu)化問題,進一步研究該問題的求解算法,所提出的多目標(biāo)遺傳算法,可為工程應(yīng)用提供相關(guān)借鑒。

    1" 問題描述

    對本文所提出的任務(wù)調(diào)度方案優(yōu)化問題進行闡釋如下[5]。

    假設(shè)設(shè)備共有M臺,數(shù)學(xué)符號表示為Ei(i=1,2,…,M);將第i臺設(shè)備的起始位置用數(shù)學(xué)符號表示為Pi;出發(fā)時刻為t0;要求執(zhí)行任務(wù)的所有參與設(shè)備能夠在時間區(qū)間t1至t2內(nèi)完成對N個目標(biāo)Tj(j=1,2,…,N)的觀測任務(wù);設(shè)備可到達的觀測站有L個,用數(shù)學(xué)符號表示為Sk(k=1,2,…,L);求最優(yōu)設(shè)備調(diào)度計劃安排。

    上面闡釋的問題進一步分析可知,其內(nèi)涵為常見的武器和目標(biāo)之間如何進行對應(yīng)的問題,即武器目標(biāo)分配(weapon target assignment, WTA)問題[6]。目前,國內(nèi)外針對WTA問題的研究主要集中在兩方面:一是模型建立[712];二是算法求解[1322]。而模型建立又可以分為靜態(tài)WTA模型建立和動態(tài)WTA模型建立,其中前者研究較為深入。算法求解方面,國內(nèi)外學(xué)者也進行了大量的研究,主要分為傳統(tǒng)的非線性規(guī)劃問題求解算法和智能算法兩大類。傳統(tǒng)的非線性規(guī)劃問題求解算法包括分支定界法、隱枚舉法、動態(tài)規(guī)劃法和割平面法等[23],此類算法適用于小規(guī)模問題求解中,當(dāng)所求問題包含的約束條件和變量較多時,利用此類算法進行求解所消耗的時間也會大幅度增加。隨著計算機技術(shù)的發(fā)展,一些模擬自然界生物學(xué)現(xiàn)象或過程的啟發(fā)式算法得到了長足發(fā)展,比如典型的遺傳算法、粒子群優(yōu)化(particle swam optimization, PSO)算法、蟻群優(yōu)化(ant colony optimazation, ACO)算法、差分進化(differential evolution, DE)等算法,為解決大規(guī)模、復(fù)雜化的WTA問題提供了新途徑。

    針對多導(dǎo)彈對多目標(biāo)攔截的目標(biāo)分配問題,王儲等[24]在遺傳算法的基礎(chǔ)上,進一步做了優(yōu)化改進,以自適應(yīng)策略為基礎(chǔ),設(shè)計了一種多種群精英遺傳算法。以艦載聯(lián)合火力打擊為背景,陳曼等[25]以打擊失敗概率最小和消耗武器數(shù)量最少為優(yōu)化指標(biāo),設(shè)計了一種改進的多目標(biāo)PSO算法??紤]到基本的PSO算法有可能會陷入到局部解,陳曼等[26]將遺傳算法操作步驟和PSO的操作步驟相互融合,將遺傳算法的交叉算子與PSO相結(jié)合,提出了一種改進的PSO,并采用模擬退火的思路對種群中各粒子進行了更新,有效降低了算法陷入局部最優(yōu)的可能。針對靜態(tài)分配模型中存在的不足,張彥芳等[27]提出了一種基于ACO算法的準(zhǔn)動態(tài)防空武器分配算法,在一定程度上能夠適應(yīng)戰(zhàn)場態(tài)勢的動態(tài)變化。為提高WTA問題的求解速度和精度,吳文海等[28]提出了一種基于隨機鄰域的自適應(yīng)差分計算算法,充分結(jié)合了“最優(yōu)性”和“隨機性”,有效解決了動態(tài)WTA問題的尋優(yōu)精度和速度。為避免算法迭代過程中陷入局部最優(yōu)解,王力超等[29]引入云自適應(yīng)模型和鯰魚效應(yīng)思想,提出了一種改進的PSO算法。

    隨著國內(nèi)外學(xué)者對WTA問題研究的不斷深入,一些其他行之有效的方法也不斷涌現(xiàn)出來。如,針對現(xiàn)有智能算法在求解過程中容易陷入局部最優(yōu)的特點,邱少明等[30]在鯨魚優(yōu)化算法的基礎(chǔ)上,通過引入差分計算算法的變異策略和自適應(yīng)位置更新,提出了一種改進的鯨魚優(yōu)化算法來解決目標(biāo)分配問題。針對聯(lián)合遠程打擊WTA問題,劉雙雙等[31]提出了一種結(jié)合小生境淘汰思想的蝙蝠算法。針對動態(tài)WTA問題,張先劍等[32]構(gòu)建了一種綜合數(shù)學(xué)模型,該模型能夠有效考慮到雙方的博弈過程和策略,針對所建立的模型,并利用納什均衡和Pareto算法進行了逐步求解。

    雖然上述關(guān)于求解WTA問題的各種研究方法為解決WTA問題提供了一定的參考,但在模型的建立和問題求解過程中,仍存在一些假設(shè)和簡化條件。未來戰(zhàn)爭將是體系與體系之間的對抗,大規(guī)模、高動態(tài)、不確定性環(huán)境下的WTA問題仍是研究的重點。

    本文所提出的車載光學(xué)測量設(shè)備任務(wù)調(diào)度方案問題屬于一種動態(tài)、廣義的WTA問題。與傳統(tǒng)的WTA問題相比,該問題同時考慮了時間窗口和測站兩個方面因素。因此,該方案包含了4個元素,即觀測設(shè)備、觀測站、時間窗口和目標(biāo),四者相互關(guān)聯(lián)為建模和求解增加了困難。針對此問題,文獻[5]創(chuàng)新性地進行建模,采用一個一維數(shù)組來表示調(diào)度方案,并給出相關(guān)約束和評價指標(biāo)涵義。本文在此基礎(chǔ)上,進一步研究該問題的求解算法,提出了一種基于非支配排序的遺傳算法(non-dominated sorting genetic algorithm II, NSGA-II)的多目標(biāo)遺傳算法,可為車載光學(xué)測量設(shè)備的工程實踐提供一定的參考價值。

    2" 問題建模

    一般多目標(biāo)優(yōu)化問題的數(shù)學(xué)描述[33]如下:

    min F(xi)=[f1(xi),f2(xi),…,fm(xi)]T(1)

    x=(x1,x2,…,xn), y=F(x)(2)

    s.t.x=(x1,x2,…,xn)∈x∈Rm

    式中:x為決策空間;n為決策變量個數(shù);y表示觀測窗口所選用的裝備序號取值為[0,M];目標(biāo)函數(shù)向量F(x)包含m(m≥2)個優(yōu)化目標(biāo);Y∈Rm為目標(biāo)函數(shù)空間;f:Rn→Rm為目標(biāo)映射函數(shù)。

    針對于本文所研究的問題,建模過程分為以下3個步驟:① 對觀測方案進行數(shù)學(xué)描述;② 確定過程中所包含的約束;③ 確定方案優(yōu)化指標(biāo)。

    2.1" 觀測方案

    將M臺設(shè)備的一個觀測方案抽象為一個一維數(shù)列來表示,如圖1所示。

    其中,qtotal=∑Nj=1∑Lk=1qjk表示窗口數(shù)目,L表示測站個數(shù),N表示目標(biāo)個數(shù),qjk表示設(shè)備部署在Sk測站對Tj目標(biāo)的觀測窗口數(shù)目。

    從圖1中可以看出,該數(shù)列中數(shù)據(jù)個數(shù)等于窗口個數(shù),所有設(shè)備可用窗口中的第x個時間窗口可以由該數(shù)列中的第x個元素來表示。x為一個整數(shù),用y表示其取值,取值范圍為0~M。這樣就可以看出,所有設(shè)備中的第y臺設(shè)備利用第x個窗口對相應(yīng)的目標(biāo)進行觀測。

    2.2" 方案約束

    判斷一個方案是否可行,首先要判斷該觀測方案是否滿足相應(yīng)的約束條件。本文給出的約束包括設(shè)備轉(zhuǎn)場約束、測站進出約束、剩余觀測次數(shù)約束以及方案指標(biāo)約束。

    2.2.1" 設(shè)備轉(zhuǎn)場約束

    設(shè)備由初始位置轉(zhuǎn)場至第1個窗口對應(yīng)測站需滿足的時間約束可表示為

    t0+tmove(Pi,Sk1)+topen-tLi,j1,k1,q1≤0(3)

    式中:tmove(·)為裝備從初始位置到終點位置的運動時間。

    設(shè)備從第α個時間窗口對應(yīng)的觀測站轉(zhuǎn)至第α+1個窗口對應(yīng)的觀測站需滿足的時間約束為

    tRi,jα,kα,qα+tcool-tLi,jα+1,kα+1,qα+1≤0, kα=kα+1

    tRi,jα,kα,qα+tclose+tmove(Skα,Skα+1)+topen-tLi,jα+1,kα+1,qα+1≤0, kα≠kα+1(4)

    對每一臺設(shè)備,循環(huán)判斷是否滿足轉(zhuǎn)場約束,若不滿足,則停止計算,即可判定該方案不滿足該約束條件。

    2.2.2" 測站進出約束

    假設(shè)設(shè)備在第k個觀測站對目標(biāo)進行觀測時,滿足要求的時間窗口小于1個,則認(rèn)為該觀測站對設(shè)備進出無約束。

    但當(dāng)設(shè)備在第k個觀測站對目標(biāo)進行觀測時,滿足要求的時間窗口大于等于2個,則相鄰的兩個時間窗口β和β+1,需要滿足如下測站約束:

    tRiβ,jβ,k,qβ+tcool-tLiβ+1,jβ+1,k,qβ+1≤0, iβ=iβ+1

    tRiβ,jβ,k,qβ+tclose+topen-tLiβ+1,jβ+1,k,qβ+1≤0, iβ≠iβ+1(5)

    2.2.3" 剩余觀測次數(shù)約束

    剩余觀測次數(shù)約束是指在編排觀測任務(wù)時,分配給一臺設(shè)備的觀測任務(wù)數(shù)不應(yīng)超過其剩余觀測次數(shù)。具體的計算過程如下。

    步驟 1" 對于其中一個假設(shè)可行的方案,在第2.2.1節(jié)對觀測設(shè)備分析是否滿足在各場地轉(zhuǎn)換時間約束的時候,已經(jīng)可以同步計算出第i臺設(shè)備的所有時間窗口,把所有這些窗口看作一個集合,由該集合可知,分配給第i臺設(shè)備的觀測任務(wù)數(shù)為αi,則其對應(yīng)的剩余觀測次數(shù)約束表示為

    αi-Ri≤0(6)

    式中:Ri為第i臺設(shè)備總的剩余觀測次數(shù)。

    步驟 2" 重復(fù)步驟1,對選定觀測方案中對應(yīng)的設(shè)備進行判定,當(dāng)涉及的所有設(shè)備剩余觀測次數(shù)不足時,停止上述循環(huán)計算過程,即可判定該觀測方案不滿足要求。

    2.2.4" 方案指標(biāo)約束

    對多個觀測方案進行評價時可引入多種評價指標(biāo),針對不作為最優(yōu)指標(biāo)的評價指標(biāo)可作為指標(biāo)約束進行考慮,即使得不作為優(yōu)化指標(biāo)的評價指標(biāo)滿足設(shè)置的指標(biāo)范圍。優(yōu)化指標(biāo)約束的計算步驟如下。

    步驟 1" 針對優(yōu)化配置條件,自動識別約束指標(biāo)以及對應(yīng)的指標(biāo)范圍,為不等式向量構(gòu)建提供輸入。

    步驟 2" 構(gòu)建動態(tài)指標(biāo)約束不等式向量,并明確優(yōu)化求解模型的不等式規(guī)模。

    步驟 3" 根據(jù)設(shè)計的評價指標(biāo)體系計算各指標(biāo)值,并將納入不等式約束的評價指標(biāo)進行不等式約束計算,生成指標(biāo)不等式約束數(shù)組值,判斷是否滿足指標(biāo)約束。

    步驟1和步驟2均在初始化完成,只需進行一次操作,步驟3循環(huán)計算直至完成優(yōu)化計算生成優(yōu)化結(jié)果。

    上述4個約束即為一個合理的觀測方案所需滿足的約束條件。這里需要說明的是,設(shè)備轉(zhuǎn)場約束和測站進出約束的詳細(xì)描述,以及約束表達式中對應(yīng)的參數(shù)含義可具體參考文獻[5],此處不再贅述。

    2.3" 方案指標(biāo)

    方案指標(biāo)是判斷觀測方案優(yōu)劣的重要依據(jù),本文依照完備性、獨立性、定量性的選取原則,選取總機動距離、總觀測時長以及觀測目標(biāo)數(shù)目作為方案指標(biāo)。

    2.3.1" 總機動距離

    總機動距離描述了所選某一觀測方案中,所有涉及到的設(shè)備在從一個觀測站機動到另一個觀測站時,所設(shè)備行進的距離的加和。一般來說,總的距離越小,該方案的成本也就越小。因此,在選擇觀測方案時,應(yīng)該使總機動距離最小。其計算步驟如下。

    步驟 1" 對于其中一個假設(shè)可行的方案,在第2.2.1節(jié)對觀測設(shè)備分析是否滿足在各場地轉(zhuǎn)換時間約束的時候,已經(jīng)可以同步計算出第i臺設(shè)備的所有時間窗口,把所有這些窗口看作一個集合Wi,jα,kα,qα(α=1,2,…,αi),則該設(shè)備的機動距離可表示為

    di=dmove(Pi,Sk1)+∑αi-1α=1dmove(Skα,Skα+1)(7)

    式中:dmove(·)表示設(shè)備起點和終點之間的距離。

    步驟 2" 總機動距離可表示為

    dtotal=∑Mi=1di(8)

    這里進一步說明的是,如果在一個方案中,某一個設(shè)備沒有領(lǐng)取到相應(yīng)的任務(wù),則該設(shè)備的機動距離就可以記為0。

    2.3.2" 總觀測時長

    總觀測時長表示所選方案中,涉及到的所有時間窗口的長度和。為了能夠盡可能對目標(biāo)長時間觀測,以獲得較多數(shù)據(jù),總觀測時長越長越好。具體的計算過程如下。

    步驟 1" 對于其中一個假設(shè)可行的方案,在第2.2.1節(jié)對觀測設(shè)備分析是否滿足在各場地轉(zhuǎn)換時間約束的時候,已經(jīng)可以同步計算出第i臺設(shè)備的所有時間窗口,把所有這些窗口看作一個集合Wi,jα,kα,qα(α=1,2,…,αi),則該設(shè)備對目標(biāo)進行觀測的時長可表示為

    Δti=∑αiα=1(tRi,jα,kα,qα-tLi,jα,kα,qα)(9)

    步驟 2" 總觀測時長可表示為

    Δttotal=∑Mi=1Δti(10)

    這里進一步說明的是,如果在一個方案中,某一個設(shè)備沒有領(lǐng)取到相應(yīng)的任務(wù),則該設(shè)備的總觀測時長就可以記為0。

    2.3.3" 觀測目標(biāo)數(shù)目

    觀測目標(biāo)數(shù)目是指所有設(shè)備能夠看到的目標(biāo)的總數(shù),一般來說,觀測到的目標(biāo)越多,該方案所能獲取的效益越大。觀測到的目標(biāo)數(shù)目計算過程較為簡單,即針對某一個假設(shè)可行的方案,根據(jù)方案對應(yīng)的數(shù)列,計算觀測到的目標(biāo)的個數(shù)Ninvol即可。

    3" 問題求解

    前文已經(jīng)建立了車載測量設(shè)備觀測任務(wù)調(diào)度問題的數(shù)學(xué)模型,從約束條件和方案指標(biāo)可以看出,該問題屬于一個多目標(biāo)求解問題,本文擬采取一種多目標(biāo)遺傳算法進行求解。傳統(tǒng)的遺傳算法只有一個目標(biāo)函數(shù),遺傳操作過程中通常包括3個遺傳算子:選擇算子、交叉算子和變異算子。對于本文中有多個目標(biāo)函數(shù)待求解的問題,采用的多目標(biāo)遺傳算法可以和單目標(biāo)遺傳算法選取一樣的編碼策略。在遺傳操作過程中,交叉和變異時亦可采用相同的交叉算子和變異算子,而在選擇過程中使用到的選擇算子,和在進行適應(yīng)度即目標(biāo)函數(shù)計算時,就要根據(jù)Pareto最優(yōu)的概念來確定。因此,本文所采用的多目標(biāo)遺傳算法,在計算適應(yīng)度即優(yōu)化目標(biāo)函數(shù)的時候,主要是基于NSGA-II中的快速非優(yōu)超排序方法。

    3.1" NSGA-II原理

    和傳統(tǒng)的NSGA算法相比較,本文所提出的NSGA-II算法具有3點優(yōu)勢[34]。第一,基于劃分等級的思想,提出了一種新的快速非優(yōu)超排序算法,進一步簡化了計算過程;第二,引入了一個擁擠距離的概念和一個擁擠距離比較算子,如此便能夠?qū)ν粚蛹壦性氐哪繕?biāo)函數(shù)值進行標(biāo)定,替代了需要指定共享半徑的適應(yīng)度共享策略,保證了種群的多樣性;第三,引入了精英保留策略,將父代個體與子代個體進行合并,并進行飛支配排序,在一定程度上擴大了樣本空間,將排序優(yōu)先級較高的個體選入下一代種群,保證了個體的優(yōu)良率,提高了種群的整體水平。

    NSGA-II算法步驟簡要說明如下。

    步驟 1" 群體排序

    遺傳算法的一般順序為選擇、交叉、變異,在本文算法運用選擇操作之前,對種群進行劃分等級,劃分的標(biāo)準(zhǔn)主要是該種群中各個體的非劣解的水平。具體包括以下幾個劃分步驟。

    步驟 1.1" 把當(dāng)前種群中的所有非劣解的個體都劃分在同一個等級,且將該等級序號設(shè)為l。

    步驟 1.2" 將等級序號為1的個體從當(dāng)前的種群中刪除,這樣便可得到一個新的種群,在該新的種群中進一步找出非劣解,將該批次非劣解劃分為同一等級,等級設(shè)為2。

    步驟 1.3" 重復(fù)第一和第二步的分級排序過程,當(dāng)種群中所包含的個體全部被劃分為對應(yīng)的層級后,停止操作步驟。

    步驟 2" 確定擁擠距離

    在NSGA算法中,需要事先指定共享半徑,這對操作者經(jīng)驗要求較高。為了解決這一不足,NSGA-II引入了擁擠距離的概念,即每個個體與所在同一層級相鄰兩個個體之間的距離,用符號id表示。擁擠距離實質(zhì)上為一個最小長方形,在幾何圖形上,這個長方形表現(xiàn)出這樣的一個性質(zhì),即個體i位于長方形內(nèi)部,而且僅有該個體一個。同時,如果某個所在的長方形較大,即表現(xiàn)出的擁擠距離大,那么其就有更大的幾率去參與種群的繁殖和進化過程,這樣的話就會在一定程度上增加種群的多樣化。

    假設(shè)一個種群中有l(wèi)個非支配解,所有的解可以組成一個集合L,則有l(wèi)=|L|。則可采用以下步驟計算擁擠距離計算:

    步驟 2.1" 對給定種群中的每一個個體i,令L[i]d=0。

    步驟 2.2" 對給定種群中每一個目標(biāo)函數(shù)m,令L=sort(L,m),L[1]d=L[l]d=∞,其中,種群中的每個個體分層級排序操作用sort(L,m)表示。

    步驟 2.3" 從i=2到l-1循環(huán)計算判斷:

    L[i]d=L[i]d+(L[i+1]m-L[i-1]m)

    在上述操作步驟中,優(yōu)化問題的目標(biāo)函數(shù)包含多個,則該種群的非支配解集L中每一個個體對應(yīng)多個目標(biāo)函數(shù),其中第m個目標(biāo)函數(shù)的值可以用L[i]m計算得出。

    步驟 3" 擁擠距離比較

    在應(yīng)用NSGA-II算法進行求解時,希望最終得到的結(jié)果可以收斂至某個Pareto曲面,且分布記為均勻。因此,需進一步對個體進行擁擠距離比較。在完成步驟1中的群體排序和步驟2中的每個個體都確定了擁擠距離后。此時,種群所包含的每個個體i都具備了兩個性質(zhì):

    (1) 每個個體都具備了非支配序?qū)蛹壭蛱杋r;

    (2) 每個個體都具備了擁擠距離id。

    這里定義如下偏序關(guān)系:≥n為i≥nj,如果irgt;jr,或者ir=jr且idgt;jd。該偏序關(guān)系表示為:如果群體中的兩個個體經(jīng)過分層級排序后,層級序號不同,那么就選取層級序號小的個體;如果兩個個體在同一層級,即層級序號一樣,則根據(jù)擁擠距離進行判斷,如果該個體周圍較為稀疏,則選取該個體。

    步驟 4" 精英保留

    第一,需要根據(jù)父代種群和子代種群生成一個新的種群Rt=Pt∪Qt,種群Rt內(nèi)包含2N個個體。第二,對新種群Rt進行群體排序,確定個體層級號,將當(dāng)前得到的第一層級個體作為新的父代種群個體,當(dāng)新種群個體數(shù)量為N后,得到新的父代種群Pt+1。根據(jù)種群Pt+1,繼續(xù)采用遺傳操作,形成子代種群Qt+1。

    算法主體流程操作如下。

    步驟 4.1" 對給定種群中每一個個體進行循環(huán)計算,同時標(biāo)記出種群中各個個體之間的支配關(guān)系,為每個個體的Pareto秩賦值為0。

    步驟 4.2" 剔除種群中的Pareto最優(yōu)解。

    步驟 4.3" 種群經(jīng)過步驟4.2操作后,將種群中其他剩下的個體的Pareto秩都加1,然后進一步剔除剩余種群中的Pareto最優(yōu)解。

    步驟 4.4" 重復(fù)進行步驟4.3,直到所有個體都被剔除。

    步驟 4.5" 用一個大的實數(shù)減去經(jīng)過上述步驟后,每個個體得到的Pareto秩,將得到的結(jié)果作為每一個個體的適應(yīng)度值。

    步驟 4.6" 對每個個體采用錦標(biāo)賽選擇操作,根據(jù)適應(yīng)度進行判斷,適應(yīng)度大的個體保留下來,如果存在個體之間的適應(yīng)度相同,則進行步驟4.7。

    步驟 4.7" 計算種群中所有Pareto秩一樣的個體的目標(biāo)函數(shù),比較目標(biāo)函數(shù)之間的差,以目標(biāo)函數(shù)差的最小值表示該種群中個體之間的擁擠密度。在步驟4.6中,采用錦標(biāo)賽選擇操作時,擁擠密度小的個體,即種群個體對應(yīng)目標(biāo)函數(shù)差異最小值大的個體生存。

    3.2" 約束處理方法

    觀測方案優(yōu)選問題具有較強的非線性約束,若采用常規(guī)的遺傳算法流程將很難獲取可行觀測方案。結(jié)合問題屬性,以獲取可行觀測方案為第一準(zhǔn)則,將約束處理分為兩部分進行處理:

    (1) 在初始化、交叉、變異過程中進行約束處理,以保證每次優(yōu)化計算均可產(chǎn)生可行觀測方案。

    (2) 在計算目標(biāo)函數(shù)與約束條件值時進行約束懲罰處理,以保證所有優(yōu)化結(jié)果均滿足約束條件。

    具體處理策略如下。

    (1) 初始化

    針對執(zhí)行觀測任務(wù)的所有觀測設(shè)備逐個進行初始化操作,對每一個設(shè)備進行初始化后,判斷是否存在重復(fù)窗口,如果存在重復(fù)窗口,則需要重新進行初始化。當(dāng)所有設(shè)備都被初始化完畢后,根據(jù)所設(shè)定的約束條件進行是否滿足約束判斷,如果不滿足約束條件,則需要對觀測設(shè)備重新初始化。

    (2) 交叉與變異

    當(dāng)種群執(zhí)行交叉和變異操作時,對不同個體的同一窗口進行重復(fù)性判斷,即觀察同一窗口是否對應(yīng)于不同的觀察設(shè)備。如果存在此種情況,則對該窗口對應(yīng)的后一個觀測設(shè)備進行重新隨機初始化。重復(fù)進行上述窗口重復(fù)性計算,當(dāng)重復(fù)次數(shù)大于設(shè)定的某一個閾值時,將該窗口對應(yīng)的后一個觀測設(shè)備值設(shè)為0。為減少計算時間,此過程不判斷其他約束條件是否滿足。

    (3) 約束懲罰

    根據(jù)優(yōu)化過程中的優(yōu)化變量計算所有約束值,若不滿足約束,對目標(biāo)函數(shù)進行懲罰處理,剔除不滿足約束的打擊方案。

    (4) 優(yōu)化自適應(yīng)設(shè)計

    為適應(yīng)測量設(shè)備數(shù)量的任意變化輸入,開展優(yōu)化自適應(yīng)設(shè)計。優(yōu)化自適應(yīng)設(shè)計主要體現(xiàn)在各段優(yōu)化變量根據(jù)用戶輸入數(shù)據(jù)自動構(gòu)建優(yōu)化變量結(jié)構(gòu),確定優(yōu)化模型,進行優(yōu)化計算。優(yōu)化自適應(yīng)設(shè)計根據(jù)測量設(shè)備與窗口對應(yīng)測量設(shè)備類型進行自動匹配,確定優(yōu)化設(shè)計變量維數(shù)和段數(shù),根據(jù)用戶輸入的優(yōu)化目標(biāo)自動判斷優(yōu)化模型類型。

    3.3" 算法流程設(shè)計

    多目標(biāo)優(yōu)化問題的求解算法流程,如圖2所示。

    4" 仿真算例

    設(shè)定某一場景,確定相應(yīng)的初始參數(shù),按照國際標(biāo)準(zhǔn)時間生成格式(universal time code generator, UTCG),對時刻t0取值為:2018年5月24日0時0分0秒;對時刻t1的取值為2018年6月1日0時0分0秒; 對時刻t2的取值為2018年6月8日0時0分0秒。

    假設(shè)觀測車的最大行駛速度為60 km/h,到達觀測站后,工作人員將車載觀測設(shè)備進行展開耗時為topen=3 h,當(dāng)觀測任務(wù)結(jié)束后,工作人員將車載觀測設(shè)備進行收攏耗時為tclose=1 h,當(dāng)設(shè)備完成某一次觀測任務(wù)后,進行下一次觀測任務(wù)前,設(shè)備冷卻耗時為tcool=1 h。

    設(shè)定某一次觀測任務(wù)需要6臺車載觀測設(shè)備,其基本參數(shù)如表1所示。假設(shè)該批次車載觀測設(shè)備需要對4個空間衛(wèi)星目標(biāo)實施觀測,以收集目標(biāo)數(shù)據(jù)??臻g衛(wèi)星在t1時刻的經(jīng)典軌道六根數(shù)如表2所示。同時,設(shè)定完成該次觀測任務(wù)時,有4個觀測站可以供車載觀測設(shè)備使用,其參數(shù)如表3所示。

    對于設(shè)備觀測時間窗口的獲取,本文借助了商業(yè)軟件衛(wèi)星工具箱(satellite tool kit, STK)中的衛(wèi)星過境計算工具箱Access來計算。具體步驟如下。

    步驟 1" 在STK建立一個場景,場景中增添4個觀測站模塊,在每一個觀測站上增設(shè)一個傳感器對象Sensor模塊,在參數(shù)設(shè)置中,把Sensor類型簡化為一個圓錐形,該圓錐形的錐角取20°。

    步驟 2" 利用STK中的Access模塊,計算4個觀測站對不同目標(biāo)進行觀測時所具備的時間窗口。

    根據(jù)上述步驟,本文利用Acess模塊便可獲取的4個觀測站對4個空間目標(biāo)進行觀測時所具備的所有時間窗口,如表4所示。

    由表4可知,在本算例中,觀測窗口全集共有31個元素。因此,本算例建立的最優(yōu)化模型的自變量為長度為31的一維數(shù)組。

    以總觀測時長最長和觀測目標(biāo)次數(shù)最少為優(yōu)化指標(biāo),并將總機動距離限制在6 000 km以內(nèi)。采用本文所給出的多目標(biāo)優(yōu)化求解算法,得到的觀測方案共計1 358個。

    以總觀測時長最長降序排列的最優(yōu)觀測方案如表5所示。該方案對應(yīng)的優(yōu)化指標(biāo)總觀測次數(shù)為17次,總觀測時長為399 s,總的機動距離為5 459.63 km,滿足總距離限制約束。

    以總觀測次數(shù)最少升序排列的最優(yōu)觀測方案如表6所示。該方案對應(yīng)的優(yōu)化指標(biāo)總觀測次數(shù)為4次,總觀測時長為119 s,總的機動距離為2 223 km,滿足總距離限制約束。

    5" 結(jié)" 論

    本文針對車載光學(xué)測量設(shè)備任務(wù)調(diào)度方案優(yōu)化問題,提出了一種基于NSGA-II的多目標(biāo)遺傳算法。首先,建立了包含約束、優(yōu)化指標(biāo)在內(nèi)的數(shù)學(xué)模型,通過將某些不作為本次最優(yōu)指標(biāo)的優(yōu)化指標(biāo)作為約束進行巧妙化處理,有效降低了求解難度;其次,采用基于NSGA-II中的快速非優(yōu)超排序方法來計算多目標(biāo)適應(yīng)度函數(shù)與選擇算子,多目標(biāo)優(yōu)化求解得到的Pareto最優(yōu)解集即為任務(wù)調(diào)度方案集。最后,通過仿真算例對所提出的算法進行了求解驗證,仿真結(jié)果表明該算法能夠有效解決任務(wù)調(diào)度方案優(yōu)化問題。

    參考文獻

    [1] 羅劍, 于小紅, 王杰娟. 多平臺空間目標(biāo)觀測協(xié)同調(diào)度問題[J]. 火力指揮與控制, 2021, 46(12): 8187.

    LUO J, YU X H, WANG J J. Cooperative scheduling problem space target observation of multi-platform[J]. Fire Control amp; Command Control, 2021, 46(12): 8187.

    [2] 崔嵩. 船載空間目標(biāo)觀測系統(tǒng)的任務(wù)規(guī)劃[J]. 電訊技術(shù), 2021, 61(12): 15141520.

    CUI S. Mission planning of shipborne space target observation system[J]. Telecommunication Engineering, 2021, 61(12): 15141520.

    [3] 劉瑩奇. 輕型車載大口徑光電跟測系統(tǒng)光學(xué)設(shè)計[J]. 光學(xué)技術(shù), 2012, 38(5): 583587.

    LIU Y Q. Optical design of electro-optical vehicle-based tracking system with large aperture and lightweight[J]. Optical Technique, 2012, 38(5): 583587.

    [4] 劉瑩奇, 王志, 劉欣悅, 等. 米級車載高分辨光電成像系統(tǒng)光學(xué)設(shè)計[J]. 紅外與激光工程, 2011, 40(8): 15121516.

    LIU Y Q, WANG Z, LIU X Y, et al. Optical design of vehicle-based high resolution E-O imaging system using meter class telescope[J]. Infrared and Laser Engineering, 2011, 40(8): 15121516.

    [5] 柴華, 宋旭民, 趙乾, 等. 車載光學(xué)測量設(shè)備任務(wù)調(diào)度問題建模與求解[J]. 國防科技大學(xué)學(xué)報, 2021, 43(5): 19.

    CHAI H, SONG X M, ZHAO Q, et al. Modeling and solving for task dispatch problem of vehicular optical observation equipment[J]. Journal of National University of Defense Technology, 2021, 43(5): 19.

    [6] SAHIN M A, LEBLEBICIOGLU K. A standard expert system for weapon target assignment problem[C]∥Proc.of the International Symposium on Performance Evaluation of Computer amp; Telecommunication Systems, 2009: 221224.

    [7] LU Y P, CHEN D Z. A new exact algorithm for the weapon-target assignment problem[J]. Omega, 2021, 98: 102138.

    [8] HAMZEHI S, BOGENBERGER K, FR ANECK P, et al. Combinatorial reinforcement learning of linear assignment problems[C]∥Proc.of the IEEE Intelligent Transportation Systems Conference, 2019: 33143321.

    [9] SUMMERS D S, ROBBINS M J, LUNDAY B J. An approximate dynamic progmming approach for comparing firing policies in a networked air defense environment[J]. Computers and Operations Research, 2020, 117: 104890.

    [10] LEBOUCHER C, SHIN H S, MENEC S L, et al. Optimal weapon target assignment based on an geometric approach[J]. IFAC Proceedings Volumes, 2013, 46(19): 341346.

    [11] SHALUMOV V, SHIMA T. Weapon-target-allocation strategies in multi-agent target-missile-defender engagement[J]. Journal of Guidance, Control, and Dynamics, 2017, 40(10): 24522464.

    [12] HU X W, LUO P C, ZHANG X N. IACO algorithm for weapon-target assignment problem in air combat[C]∥Proc.of the 2nd International Conference on Intelligent Systems, Metaheuristics amp; Swarm Intelligence, 2018: 3540.

    [13] KLINE A G, AHNER D K, LUNDAY B J, et al. Real-rime heuristic algorithms for the static weapon target assignment problem[J]. Journal of Heuristics, 2019, 25(3): 377397.

    [14] AHNER D K, PARSON C R. Optimal multi-stage allocation of weapons to targets using adaptive dynamic programming[J]. Optimization Letters, 2015, 9(8): 16891701.

    [15] CHANG T Q, KONG D P, HAO N, et al. Solving the dynamic weapon target assignment problem by an improved artificial bee colony algorithm with heuristic factor initialization[J]. Applied Soft Computing, 2018, 70: 845863.

    [16] XIN B, WANG Y P, CHEN J. An efficient marginal-return-based constructive heuristic to solve the sensor-weapon-target assignment problem[J]. IEEE Trans.on System, Man, and Cybernetics: Systems, 2019, 49(12): 25362547.

    [17] LI X Y, ZHOU D Y, YANG Z, et al. A novel genetic algorithm for the synthetical sensor-weapon-target assignment problem[J]. Applied Sciences, 2019, 9(18): 3803.

    [18] ZHANG K, ZHOU D Y, YANG Z, et al. Constrained multi-objective weapon target assignment for area targets by efficient evolutionary algorithm[J]. IEEE Access, 2019, 7: 176339176360.

    [19] LI J, XIN B, PARDALOS P M, et al. Solving bi-objective uncertain stochastic resource allocation problems by the CVaR-based risk measure and decomposition-based multi-objective evolutionary algorithms[J]. Annals of Operations Research, 2021, 296: 639666.

    [20] XU W Q, CHEN C, DING S X, et al. A bi-objective dynamic collaborative task assignment under uncertainty using modified MOEA/D with heuristic initialization[J]. Expert Systems with Applications, 2020, 140: 112844.

    [21] FU G Y, WANG C, ZHANG D Q, et al. A multi-objective particle swam optimization algorithm based on multi-population coevolution for weapon-target assignment[J]. Mathematical Problems in Engineering, 2019, 2019(1): 111.

    [22] ZHANG K, ZHOU D Y, YANG Z, et al. Efficient decision approaches for asset-based dynamic weapon target assignment by a receding horizon and marginal return heuristic[J]. Electronics, 2020, 9(9): 1511.

    [23] KLINE A, AHNER D, HILL R. The weapon-target assignment problem[J]. Computers and Operations Research, 2019, 105: 226236.

    [24] 王儲, 南英, 許航. 一種新的精英遺傳算法及在多彈攔截分配策略的應(yīng)用[J]. 航天控制, 2021, 39(4): 5965.

    WANG C, NAN Y, XU H. A new elite genetic algorithm and its application in multi-missile interception assignment strategy[J]. Aerospace Control, 2021, 39(4): 5965.

    [25] 陳曼, 周鳳星, 張成堯. 改進MOPSO的聯(lián)合火力打擊目標(biāo)分配[J]. 火力指揮與控制, 2019, 44(9): 125129.

    CHEN M, ZHOU F X, ZHANG C Y. Improved MOPSO joint fire strike target assignment[J]. Fire Control amp; Command Control, 2019, 44(9): 125129.

    [26] 陳曼, 周鳳星. 改進粒子群算法的艦載武器目標(biāo)分配[J]. 火力指揮與控制, 2018, 43(11): 7276.

    CHEN M, ZHOU F X. Shipborne weapon target assignment based on improved particle swarm optimization[J]. Fire Control amp; Command Control, 2018, 43(11): 7276.

    [27] 張彥芳, 閆德恒, 王冀揚, 等. 一種基于蟻群算法的準(zhǔn)動態(tài)防空武器分配算法[J]. 火力指揮與控制, 2016, 41(9): 112116.

    ZHANG Y F, YAN D H, WANG J Y, et al. A quasi-dynamic defense weapon assignment based on ant-colony optimization[J]. Fire Control amp; Command Control, 2016, 41(9): 112116.

    [28] 吳文海, 郭曉峰, 周思羽, 等. 改進差分進化算法求解武器目標(biāo)分配問題[J]. 系統(tǒng)工程與電子技術(shù), 2021, 43(4): 10121021.

    WU W H, GUO X F, ZHOU S Y, et al. Improved differential evolution algorithm for solving weapon target assignment pro-blem[J]. Systems Engineering and Electronics, 2021, 43(4): 10121021.

    [29] 王力超, 喬勇軍, 李永勝. 基于CE-CAPSO武器目標(biāo)分配優(yōu)化算法[J]. 火力指揮與控制, 2020, 45(11): 8287.

    WANG L C, QIAO Y J, LI Y S. Optimization algorithm for weapon target assignment based on CE-CAPSO[J]. Fire Control amp; Command Control, 2020, 45(11): 8287.

    [30] 邱少明, 劉良成, 張學(xué)翠, 等. 基于改進鯨魚優(yōu)化算法的武器目標(biāo)分配[J]. 火力指揮與控制, 2021, 46(7): 2731.

    QIU S M, LIU L C, ZHANG X C, et al. Research on weapon target based on improved whale optimization algorithm[J]. Fire Control amp; Command Control, 2021, 46(7): 2731.

    [31] 劉雙雙, 許瑞明, 潘俊杰. 基于小生境蝙蝠算法的聯(lián)合遠程打擊武器目標(biāo)分配問題建模與求解[J]. 裝備學(xué)院學(xué)報, 2017, 28(2): 9398.

    LIU S S, XU R M, PAN J J. Research on the modeling and solving of the joint long-range strike weapon target assignment problem based on the niche bat algorithm[J]. Journal of Equipment Academy, 2017, 28(2): 9398.

    [32] 張先劍. 空陸攻防博弈的動態(tài)武器目標(biāo)分配[J]. 國防科技大學(xué)學(xué)報, 2019, 41(2): 185190.

    ZHANG X J. Land defense weapon versus target assignment against air attack[J]. Journal of National University of Defense Technology, 2019, 41(2): 185190.

    [33] GUO X H, JI M J, ZHAO Z W, et al. Global path planning and multi-objective path control for unmanned surface vehicle based on modified particle swarm optimization (PSO) algorithm[J]. Ocean Engineering, 2020, 216: 107693.

    [34] 趙瑞. 多目標(biāo)遺傳算法應(yīng)用的研究[D]. 天津: 天津大學(xué), 2005.

    ZHAO R. Analysis about the application of the multi-objective genetic algorithm[D]. Tianjin: Tianjin University, 2005.

    作者簡介

    許強強(1990—),男,講師,博士,主要研究方向為任務(wù)規(guī)劃與評估。

    柴" 華(1988—),男,副教授,博士,主要研究方向為任務(wù)規(guī)劃與評估。

    猜你喜歡
    優(yōu)化設(shè)備
    諧響應(yīng)分析在設(shè)備減振中的應(yīng)用
    超限高層建筑結(jié)構(gòu)設(shè)計與優(yōu)化思考
    民用建筑防煙排煙設(shè)計優(yōu)化探討
    關(guān)于優(yōu)化消防安全告知承諾的一些思考
    一道優(yōu)化題的幾何解法
    由“形”啟“數(shù)”優(yōu)化運算——以2021年解析幾何高考題為例
    基于VB6.0+Access2010開發(fā)的設(shè)備管理信息系統(tǒng)
    基于MPU6050簡單控制設(shè)備
    電子制作(2018年11期)2018-08-04 03:26:08
    500kV輸變電設(shè)備運行維護探討
    如何在設(shè)備采購中節(jié)省成本
    国产一区二区在线观看av| 少妇熟女欧美另类| 少妇被粗大的猛进出69影院 | 伊人久久国产一区二区| 国产毛片在线视频| 国产精品无大码| 婷婷色麻豆天堂久久| 国产成人91sexporn| 日本wwww免费看| av天堂久久9| av女优亚洲男人天堂| 狠狠精品人妻久久久久久综合| tube8黄色片| 欧美最新免费一区二区三区| 国产一区二区在线观看av| 综合色丁香网| 大码成人一级视频| 黄色欧美视频在线观看| 亚洲图色成人| 久久精品国产亚洲av天美| 黑人巨大精品欧美一区二区蜜桃 | 亚洲精品国产色婷婷电影| 国产成人91sexporn| 制服诱惑二区| 欧美人与性动交α欧美精品济南到 | 成人黄色视频免费在线看| 亚洲国产av影院在线观看| 中国三级夫妇交换| 久久久精品94久久精品| 青春草亚洲视频在线观看| 国产成人免费无遮挡视频| 最近2019中文字幕mv第一页| 欧美激情 高清一区二区三区| av.在线天堂| 久久久久久伊人网av| 嫩草影院入口| 亚洲av成人精品一区久久| 国产精品久久久久久久电影| 成人国语在线视频| 欧美日韩av久久| 美女中出高潮动态图| 国产免费福利视频在线观看| 又粗又硬又长又爽又黄的视频| 午夜91福利影院| 精品亚洲乱码少妇综合久久| 欧美三级亚洲精品| 精品久久国产蜜桃| 麻豆精品久久久久久蜜桃| 在线播放无遮挡| 这个男人来自地球电影免费观看 | 亚洲,欧美,日韩| 蜜臀久久99精品久久宅男| 亚洲精品国产色婷婷电影| 免费黄色在线免费观看| 久久精品国产亚洲av涩爱| 欧美精品一区二区免费开放| 日日摸夜夜添夜夜爱| 国产欧美日韩一区二区三区在线 | 午夜91福利影院| 涩涩av久久男人的天堂| 交换朋友夫妻互换小说| 亚洲av电影在线观看一区二区三区| 热99国产精品久久久久久7| tube8黄色片| 欧美xxⅹ黑人| 激情五月婷婷亚洲| 18禁动态无遮挡网站| 欧美激情极品国产一区二区三区 | 黄色怎么调成土黄色| 婷婷成人精品国产| 欧美最新免费一区二区三区| 啦啦啦在线观看免费高清www| 日韩免费高清中文字幕av| 午夜福利视频在线观看免费| 男的添女的下面高潮视频| 一级a做视频免费观看| a 毛片基地| 考比视频在线观看| 国产午夜精品久久久久久一区二区三区| 欧美少妇被猛烈插入视频| 免费黄色在线免费观看| 免费大片黄手机在线观看| 91aial.com中文字幕在线观看| 亚洲欧洲日产国产| 性色av一级| 精品少妇黑人巨大在线播放| www.色视频.com| 最后的刺客免费高清国语| 又黄又爽又刺激的免费视频.| 免费高清在线观看日韩| 2018国产大陆天天弄谢| 青春草国产在线视频| 伦精品一区二区三区| 国产精品一国产av| 一本久久精品| 考比视频在线观看| 又大又黄又爽视频免费| av电影中文网址| 久久久久久久亚洲中文字幕| 男女边摸边吃奶| 在线观看国产h片| 婷婷色av中文字幕| 国产有黄有色有爽视频| 激情五月婷婷亚洲| 亚洲不卡免费看| 人体艺术视频欧美日本| 欧美另类一区| 成年av动漫网址| 大陆偷拍与自拍| 国产欧美另类精品又又久久亚洲欧美| 水蜜桃什么品种好| 国产毛片在线视频| 日韩大片免费观看网站| 国产 精品1| a级片在线免费高清观看视频| 蜜臀久久99精品久久宅男| xxxhd国产人妻xxx| 97在线人人人人妻| 汤姆久久久久久久影院中文字幕| 成人漫画全彩无遮挡| 欧美3d第一页| 久久国产精品大桥未久av| 久久久久视频综合| 亚洲av综合色区一区| 女人久久www免费人成看片| 亚洲综合色惰| 51国产日韩欧美| 人妻一区二区av| 欧美日韩视频精品一区| av不卡在线播放| 另类亚洲欧美激情| 日本wwww免费看| 色婷婷久久久亚洲欧美| 黄色配什么色好看| 免费久久久久久久精品成人欧美视频 | 卡戴珊不雅视频在线播放| 最近最新中文字幕免费大全7| 日本91视频免费播放| 精品亚洲成a人片在线观看| 亚洲精品久久成人aⅴ小说 | 国产免费现黄频在线看| 狠狠婷婷综合久久久久久88av| 国产精品无大码| 一级二级三级毛片免费看| 亚洲四区av| 亚洲图色成人| 人妻制服诱惑在线中文字幕| 久久久午夜欧美精品| 国产精品三级大全| 人妻少妇偷人精品九色| 国产不卡av网站在线观看| 三级国产精品欧美在线观看| 亚洲av欧美aⅴ国产| 欧美日韩精品成人综合77777| 桃花免费在线播放| 卡戴珊不雅视频在线播放| 国产精品久久久久久精品古装| 国产在线一区二区三区精| 人妻制服诱惑在线中文字幕| 免费高清在线观看视频在线观看| 国产成人91sexporn| 亚洲精品久久午夜乱码| 欧美精品亚洲一区二区| 免费少妇av软件| 国产av精品麻豆| 视频区图区小说| 日本av手机在线免费观看| 国产视频内射| 亚洲国产精品999| 2022亚洲国产成人精品| av天堂久久9| 午夜激情福利司机影院| 日韩在线高清观看一区二区三区| 免费观看av网站的网址| 精品午夜福利在线看| 一级毛片我不卡| 国产精品嫩草影院av在线观看| 亚洲av免费高清在线观看| 国产亚洲精品第一综合不卡 | 18禁观看日本| 精品一区二区三卡| 亚洲av二区三区四区| 国产黄色免费在线视频| 国产国拍精品亚洲av在线观看| 精品少妇内射三级| 久久久久久久国产电影| 国产一区二区三区综合在线观看 | 国产黄频视频在线观看| 精品久久久精品久久久| 国产精品久久久久久精品古装| 日韩电影二区| 国产成人精品无人区| 国产成人精品婷婷| 亚洲精品日本国产第一区| 国产色婷婷99| 这个男人来自地球电影免费观看 | 欧美精品亚洲一区二区| 午夜免费男女啪啪视频观看| 精品人妻偷拍中文字幕| 国产亚洲午夜精品一区二区久久| 国产精品久久久久成人av| 丝袜喷水一区| 永久免费av网站大全| 亚洲国产精品一区三区| av线在线观看网站| 交换朋友夫妻互换小说| 亚洲欧洲精品一区二区精品久久久 | 天堂8中文在线网| 成人二区视频| 黄色一级大片看看| 日韩中文字幕视频在线看片| 久久久久久人妻| 欧美精品人与动牲交sv欧美| av在线播放精品| 欧美一级a爱片免费观看看| 亚洲精品日韩av片在线观看| 亚洲精品日本国产第一区| 久久鲁丝午夜福利片| 久久久国产一区二区| 纯流量卡能插随身wifi吗| 人妻夜夜爽99麻豆av| av国产精品久久久久影院| 久久国产精品男人的天堂亚洲 | 18禁在线播放成人免费| 一级毛片aaaaaa免费看小| 日韩伦理黄色片| 人妻少妇偷人精品九色| 亚洲第一区二区三区不卡| xxxhd国产人妻xxx| a级片在线免费高清观看视频| 日本wwww免费看| 亚洲精品第二区| 成人亚洲欧美一区二区av| 成人18禁高潮啪啪吃奶动态图 | 国产一区二区三区av在线| 99久国产av精品国产电影| 欧美+日韩+精品| 日韩三级伦理在线观看| 亚洲成人一二三区av| 亚洲av日韩在线播放| 高清不卡的av网站| 99re6热这里在线精品视频| 中文精品一卡2卡3卡4更新| 亚洲国产成人一精品久久久| 久久99蜜桃精品久久| 亚洲av免费高清在线观看| 欧美激情极品国产一区二区三区 | 国产精品一区www在线观看| 成人亚洲精品一区在线观看| 综合色丁香网| 欧美激情 高清一区二区三区| 国产成人91sexporn| 久久精品久久久久久久性| 亚洲欧美精品自产自拍| 国语对白做爰xxxⅹ性视频网站| 午夜免费男女啪啪视频观看| 亚洲内射少妇av| 免费黄网站久久成人精品| 亚洲精品自拍成人| av女优亚洲男人天堂| 国产高清不卡午夜福利| 99热国产这里只有精品6| 一级黄片播放器| 狠狠婷婷综合久久久久久88av| 亚洲av电影在线观看一区二区三区| 国产在线一区二区三区精| 久久99蜜桃精品久久| 久久久久网色| 夜夜看夜夜爽夜夜摸| 久久99热这里只频精品6学生| 欧美97在线视频| 亚洲av成人精品一二三区| 三级国产精品片| 亚洲欧美一区二区三区黑人 | 91久久精品国产一区二区三区| 特大巨黑吊av在线直播| 国产午夜精品久久久久久一区二区三区| 丰满少妇做爰视频| 人人妻人人澡人人爽人人夜夜| 精品久久蜜臀av无| 精品人妻一区二区三区麻豆| 国产精品偷伦视频观看了| 国产精品熟女久久久久浪| 国产免费现黄频在线看| 久久青草综合色| 午夜免费男女啪啪视频观看| 夫妻午夜视频| 日韩精品免费视频一区二区三区 | 日本vs欧美在线观看视频| 天堂俺去俺来也www色官网| 高清黄色对白视频在线免费看| 观看美女的网站| 日韩精品免费视频一区二区三区 | 大片电影免费在线观看免费| 天天影视国产精品| 亚洲伊人久久精品综合| 高清毛片免费看| 如日韩欧美国产精品一区二区三区 | 中文字幕人妻熟人妻熟丝袜美| 蜜臀久久99精品久久宅男| av女优亚洲男人天堂| 成人亚洲欧美一区二区av| 一区二区三区四区激情视频| 欧美+日韩+精品| 伦理电影免费视频| 精品人妻在线不人妻| 黑人高潮一二区| 蜜桃久久精品国产亚洲av| 18在线观看网站| 亚洲四区av| 九九在线视频观看精品| 亚洲国产最新在线播放| 亚洲成人av在线免费| 免费观看无遮挡的男女| 亚洲在久久综合| 中文字幕最新亚洲高清| 亚洲精华国产精华液的使用体验| 美女国产视频在线观看| 亚洲图色成人| 九九久久精品国产亚洲av麻豆| 欧美丝袜亚洲另类| 街头女战士在线观看网站| 少妇 在线观看| 91久久精品电影网| 色婷婷久久久亚洲欧美| 纵有疾风起免费观看全集完整版| 一二三四中文在线观看免费高清| 亚洲精品一二三| 自拍欧美九色日韩亚洲蝌蚪91| 嘟嘟电影网在线观看| 这个男人来自地球电影免费观看 | 亚洲精品,欧美精品| 最近的中文字幕免费完整| 成人毛片60女人毛片免费| 我要看黄色一级片免费的| 国产免费现黄频在线看| 亚洲成人手机| 色网站视频免费| 亚洲高清免费不卡视频| 免费高清在线观看日韩| 国语对白做爰xxxⅹ性视频网站| 伊人久久精品亚洲午夜| 黄片播放在线免费| 免费看光身美女| 成人漫画全彩无遮挡| 国产成人午夜福利电影在线观看| 国产精品欧美亚洲77777| 精品一区二区三卡| 国产国拍精品亚洲av在线观看| 18禁观看日本| 日韩不卡一区二区三区视频在线| 婷婷成人精品国产| 国产一区二区在线观看日韩| 久热久热在线精品观看| 最新中文字幕久久久久| 九九爱精品视频在线观看| 国产一区二区在线观看日韩| 久久精品国产亚洲av天美| 国产极品粉嫩免费观看在线 | 在线看a的网站| 国语对白做爰xxxⅹ性视频网站| 精品酒店卫生间| 久久青草综合色| 伊人久久精品亚洲午夜| 久久影院123| 欧美激情 高清一区二区三区| 2018国产大陆天天弄谢| 国产午夜精品一二区理论片| 9色porny在线观看| 亚洲精品美女久久av网站| av国产精品久久久久影院| 亚洲在久久综合| 国产综合精华液| 国产成人一区二区在线| 午夜免费鲁丝| 亚洲中文av在线| 91精品国产国语对白视频| 伦理电影大哥的女人| 国产成人a∨麻豆精品| 国产成人精品福利久久| 能在线免费看毛片的网站| 99国产综合亚洲精品| 99九九在线精品视频| 国产男女超爽视频在线观看| 高清欧美精品videossex| 永久网站在线| 熟女人妻精品中文字幕| 亚洲精品,欧美精品| 久久久精品区二区三区| 欧美bdsm另类| 新久久久久国产一级毛片| 我的女老师完整版在线观看| 男人添女人高潮全过程视频| 亚洲精品国产色婷婷电影| 国产精品不卡视频一区二区| 丰满迷人的少妇在线观看| 欧美人与善性xxx| 国产av国产精品国产| 少妇高潮的动态图| 丰满少妇做爰视频| 久久99一区二区三区| 在现免费观看毛片| 亚洲av中文av极速乱| 三级国产精品片| 午夜av观看不卡| 亚洲av不卡在线观看| 亚洲欧美成人综合另类久久久| 18禁观看日本| 97超视频在线观看视频| 日韩电影二区| 国产日韩一区二区三区精品不卡 | 各种免费的搞黄视频| 免费少妇av软件| 亚洲欧洲日产国产| 亚洲无线观看免费| 免费观看性生交大片5| 国产成人精品福利久久| 亚洲欧美清纯卡通| 欧美精品人与动牲交sv欧美| 在线观看www视频免费| 最黄视频免费看| 亚洲精品乱码久久久v下载方式| 日韩伦理黄色片| av黄色大香蕉| 777米奇影视久久| 精品一品国产午夜福利视频| 成人18禁高潮啪啪吃奶动态图 | 男女边吃奶边做爰视频| 人妻一区二区av| 日日爽夜夜爽网站| av.在线天堂| 精品人妻熟女av久视频| 国产精品一国产av| 精品人妻一区二区三区麻豆| 日本猛色少妇xxxxx猛交久久| 久久精品熟女亚洲av麻豆精品| 2018国产大陆天天弄谢| 久久av网站| 中文字幕av电影在线播放| 老司机影院毛片| 免费日韩欧美在线观看| freevideosex欧美| 纯流量卡能插随身wifi吗| 久久久久久伊人网av| 免费少妇av软件| 最近最新中文字幕免费大全7| 午夜福利视频在线观看免费| 丰满乱子伦码专区| 成年人免费黄色播放视频| 精品国产国语对白av| 高清黄色对白视频在线免费看| 曰老女人黄片| 亚洲美女黄色视频免费看| 三上悠亚av全集在线观看| 成人国产av品久久久| 亚洲av福利一区| 午夜视频国产福利| 久久久久久久精品精品| 精品国产一区二区三区久久久樱花| 一本大道久久a久久精品| 91久久精品电影网| 成人二区视频| 最近最新中文字幕免费大全7| 亚洲第一av免费看| 久久ye,这里只有精品| 黄色配什么色好看| 2021少妇久久久久久久久久久| 人妻夜夜爽99麻豆av| 国产乱来视频区| 看非洲黑人一级黄片| 少妇精品久久久久久久| 亚洲欧美日韩另类电影网站| 国产精品无大码| 99热这里只有精品一区| 成人国语在线视频| 在线播放无遮挡| 91精品国产国语对白视频| 亚洲欧美精品自产自拍| 国产69精品久久久久777片| 日本wwww免费看| 免费人成在线观看视频色| 久久影院123| 亚洲经典国产精华液单| 国产片特级美女逼逼视频| 在线观看免费日韩欧美大片 | 国产男人的电影天堂91| 亚洲国产精品999| 大片免费播放器 马上看| 人成视频在线观看免费观看| 国产一级毛片在线| 亚州av有码| 免费看不卡的av| 午夜福利视频精品| 91久久精品电影网| 纯流量卡能插随身wifi吗| av又黄又爽大尺度在线免费看| 人妻 亚洲 视频| 久久精品久久精品一区二区三区| 大陆偷拍与自拍| 国产精品熟女久久久久浪| 寂寞人妻少妇视频99o| 美女福利国产在线| 女性被躁到高潮视频| 人体艺术视频欧美日本| 美女内射精品一级片tv| 亚洲国产日韩一区二区| 综合色丁香网| www.av在线官网国产| 我要看黄色一级片免费的| 亚洲av中文av极速乱| 亚洲不卡免费看| 国产成人精品婷婷| 久久鲁丝午夜福利片| 国产精品嫩草影院av在线观看| 美女大奶头黄色视频| 91精品国产九色| 国产精品国产三级国产av玫瑰| 午夜视频国产福利| 国产精品久久久久久av不卡| 母亲3免费完整高清在线观看 | 久久久久视频综合| 久久久久久久久久成人| 国产精品久久久久成人av| 街头女战士在线观看网站| 午夜福利视频在线观看免费| 天堂8中文在线网| 久久人妻熟女aⅴ| 国产无遮挡羞羞视频在线观看| 精品一区二区免费观看| av视频免费观看在线观看| av卡一久久| a级毛片在线看网站| 亚洲不卡免费看| 欧美另类一区| 老司机影院毛片| 少妇熟女欧美另类| 我的女老师完整版在线观看| 亚洲性久久影院| av不卡在线播放| 免费观看a级毛片全部| 成人综合一区亚洲| 国精品久久久久久国模美| 中文天堂在线官网| 亚洲精品乱码久久久v下载方式| 一边亲一边摸免费视频| 亚洲中文av在线| 狂野欧美白嫩少妇大欣赏| 亚洲国产av影院在线观看| 热99久久久久精品小说推荐| 欧美日韩国产mv在线观看视频| 日韩人妻高清精品专区| 成年av动漫网址| 激情五月婷婷亚洲| 99视频精品全部免费 在线| 丰满饥渴人妻一区二区三| 亚洲精品色激情综合| 久久久久久久久久成人| 国产片特级美女逼逼视频| 久久久久精品久久久久真实原创| 国产免费又黄又爽又色| 久久97久久精品| 日韩人妻高清精品专区| 久久精品国产鲁丝片午夜精品| 黑人猛操日本美女一级片| 国产成人精品在线电影| 秋霞在线观看毛片| 免费不卡的大黄色大毛片视频在线观看| 日韩在线高清观看一区二区三区| 美女内射精品一级片tv| 一区二区三区精品91| 天天影视国产精品| 91久久精品国产一区二区成人| 成人毛片60女人毛片免费| 九九爱精品视频在线观看| 亚洲精品国产色婷婷电影| 高清午夜精品一区二区三区| 久热久热在线精品观看| √禁漫天堂资源中文www| 亚洲成人av在线免费| 超碰97精品在线观看| 国产国拍精品亚洲av在线观看| 亚洲国产日韩一区二区| 九色亚洲精品在线播放| av免费在线看不卡| 成年人免费黄色播放视频| 51国产日韩欧美| 女人精品久久久久毛片| 不卡视频在线观看欧美| 亚洲人与动物交配视频| 一级毛片 在线播放| 欧美日本中文国产一区发布| 国产成人aa在线观看| av免费在线看不卡| 久久鲁丝午夜福利片| 午夜激情福利司机影院| 黑人高潮一二区| 日韩熟女老妇一区二区性免费视频| 国产黄色视频一区二区在线观看| 亚洲av成人精品一二三区| 尾随美女入室| 成年人午夜在线观看视频| 国产一区二区在线观看av| 欧美人与善性xxx| 欧美精品亚洲一区二区| 日韩大片免费观看网站| 黄色一级大片看看| 久久精品国产亚洲网站| 啦啦啦在线观看免费高清www| www.av在线官网国产| 国产乱来视频区| 久久久国产一区二区| 欧美日韩综合久久久久久| 晚上一个人看的免费电影| 中文字幕最新亚洲高清| 国产免费一级a男人的天堂|