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

    多目標(biāo)優(yōu)化的艦船磁場(chǎng)建模方法*

    2021-09-03 08:26:52戴忠華周穗華張曉兵
    物理學(xué)報(bào) 2021年16期
    關(guān)鍵詞:磁偶極子磁矩艦船

    戴忠華 周穗華 張曉兵

    (海軍工程大學(xué)兵器工程學(xué)院504教研室, 武漢 430033)

    針對(duì)艦船磁場(chǎng)混合模型建模中存在的建模精度不高和穩(wěn)定性差的問(wèn)題, 提出一種高精度穩(wěn)定模型建立方法, 結(jié)合混合模型中磁偶極子參數(shù)與艦船結(jié)構(gòu)的相關(guān)性, 以建模精度和模型穩(wěn)定性為目標(biāo)構(gòu)造了多目標(biāo)函數(shù), 通過(guò)對(duì)多目標(biāo)函數(shù)優(yōu)化獲得合理的磁偶極子參數(shù), 間接地將建模求解問(wèn)題轉(zhuǎn)化為多目標(biāo)函數(shù)優(yōu)化問(wèn)題.利用多目標(biāo)粒子群優(yōu)化算法進(jìn)行求解, 得到了建模問(wèn)題求解結(jié)果的可選集, 以建模精度為基準(zhǔn)設(shè)計(jì)了從可選集中選取最佳結(jié)果的選擇規(guī)則.三種類(lèi)型的艦船船模實(shí)測(cè)數(shù)據(jù)建模結(jié)果表明: 本文方法所建模型相對(duì)誤差小于3%, 換算誤差小于6%, 能夠有效對(duì)艦船磁場(chǎng)進(jìn)行建模; 當(dāng)存在測(cè)量數(shù)據(jù)誤差時(shí), 本文方法建模求解結(jié)果穩(wěn)定,驗(yàn)證了文本方法建模具有較好的穩(wěn)定性.海上的某型艦船實(shí)測(cè)數(shù)據(jù)建模結(jié)果表明, 本文方法建模具有較高的建模精度和換算精度, 能夠有效地在相關(guān)的工程中應(yīng)用.

    1 引 言

    艦船目標(biāo)通常由鋼鐵材料制成, 受到地磁場(chǎng)的磁化會(huì)產(chǎn)生艦船磁場(chǎng), 當(dāng)艦船目標(biāo)進(jìn)入某一地磁場(chǎng)區(qū)域時(shí), 會(huì)造成該區(qū)域的磁場(chǎng)異常[1,2], 水中兵器正是利用該磁異常實(shí)現(xiàn)對(duì)目標(biāo)的探測(cè)和定位[3-11].為了提高艦船的防御能力, 艦船磁隱身技術(shù)成為各國(guó)研究的熱點(diǎn), 艦船磁場(chǎng)空間分布特性的獲取是艦船磁隱身技術(shù)研究的前提和關(guān)鍵.高精度的艦船磁場(chǎng)建??梢垣@取精確的艦船磁場(chǎng)空間分布特性, 通常, 艦船建模需要測(cè)量艦船某一平面上的完整磁場(chǎng)數(shù)據(jù), 運(yùn)用一定的建模技術(shù)完成對(duì)艦船磁場(chǎng)的建模.然而, 由于平面測(cè)量數(shù)據(jù)的不完整以及測(cè)量過(guò)程中存在的測(cè)量誤差, 導(dǎo)致所建模型存在建模精度低、穩(wěn)定性差、適用性不好等問(wèn)題, 如何從有限測(cè)量數(shù)據(jù)中獲得艦船的高精度穩(wěn)定模型是艦船磁場(chǎng)建模的研究重點(diǎn).

    目前, 主要的艦船磁場(chǎng)建模方法大致可以分為兩類(lèi): 一類(lèi)是依據(jù)嚴(yán)格的理論推導(dǎo)出來(lái)的, 如大平面法、邊界元法、有限元法等[12-14], 這類(lèi)方法計(jì)算精度高, 但對(duì)測(cè)量要求比較嚴(yán)格, 需要一個(gè)完整的測(cè)量包絡(luò)面, 在實(shí)際中很難滿(mǎn)足; 另一類(lèi)是磁體模擬法[15-17], 將艦船等效為位于艦船水線面的均勻磁化橢球體和均勻分布于吃水線上的磁偶極子陣列的混合模型, 均勻橢球體模擬艦船的宏觀磁場(chǎng),磁偶極子陣列模擬艦船的局部磁場(chǎng), 這類(lèi)方法可以在少量的磁場(chǎng)測(cè)量數(shù)據(jù)下完成精度較高的艦船磁場(chǎng)模型建立, 是工程中較為常用的方法.使用混合模型建模時(shí), 本質(zhì)是求解一個(gè)多維的超定方程, 該超定方程通常為一個(gè)病態(tài)程度較高的方程, 對(duì)輸入的敏感性極強(qiáng), 而磁場(chǎng)測(cè)量中不可避免地存在測(cè)量誤差, 導(dǎo)致建模求解結(jié)果的不穩(wěn)定, 如何求解建模病態(tài)超定方程是混合模型建模過(guò)程中的關(guān)鍵.

    對(duì)于建模病態(tài)超定方程的求解, 林春生和龔沈光[1]最早使用逐步回歸法進(jìn)行求解, 通過(guò)剔除超定方程系數(shù)矩陣中的不明顯列提高模型的穩(wěn)定性, 但當(dāng)系數(shù)矩陣病態(tài)程度過(guò)大時(shí), 該方法建模的精度較差; 楊明明等[18]使用正則化法來(lái)求解超定方程, 一定程度上解決了病態(tài)超定方程求解問(wèn)題, 但求解正確性往往依賴(lài)正則化參數(shù)的選取; 王金根等[19]首次考慮了測(cè)量數(shù)據(jù)誤差問(wèn)題, 先對(duì)誤差數(shù)據(jù)進(jìn)行剔除, 然后利用改進(jìn)的方法來(lái)對(duì)系數(shù)矩陣的冗余列進(jìn)行剔除, 一定程度上降低系數(shù)矩陣的條件數(shù).后續(xù),為了降低建模超定方程的病態(tài)程度, 多位學(xué)者通過(guò)優(yōu)化磁偶極子位置降低超定方程系數(shù)矩陣的條件數(shù)[20-25], 以增加模型穩(wěn)定性, 該類(lèi)方法在一定程度上解決了病態(tài)超定方程的穩(wěn)定求解問(wèn)題, 但為了降低系數(shù)矩陣條件數(shù)片面地優(yōu)化磁偶極子位置分布,忽略磁偶極子位置分布與艦船結(jié)構(gòu)的相關(guān)性, 使得所建模型存在精度不高、推算能力差等問(wèn)題.最近, 出現(xiàn)了考慮艦船實(shí)際結(jié)構(gòu)的艦船磁場(chǎng)反演技術(shù)[26-28], 具有計(jì)算精度高的優(yōu)點(diǎn), 然而, 這些技術(shù)需要對(duì)船體的結(jié)構(gòu)進(jìn)行較為準(zhǔn)確的三維建模, 在工程化應(yīng)用中較難實(shí)現(xiàn).同時(shí), 目前所有使用混合模型進(jìn)行建模的方法存在一個(gè)共同的問(wèn)題, 對(duì)磁偶極子個(gè)數(shù)的選取并未考慮, 而磁偶極子個(gè)數(shù)的多少對(duì)所建模型的有效性有較大的影響, 磁偶極子數(shù)過(guò)少, 模型精度較差, 偶極子數(shù)過(guò)多, 會(huì)造成局部擬合過(guò)度, 使得模型的推算性能不好.其實(shí), 混合模型中的磁偶極子個(gè)數(shù)和位置是與具體艦船的內(nèi)部結(jié)構(gòu)密切聯(lián)系的, 選擇合適的磁偶極子個(gè)數(shù)和正確分布的磁偶極子位置可以建立一個(gè)高精度穩(wěn)定的艦船磁場(chǎng)模型.

    本文將基于艦船磁場(chǎng)混合模型, 結(jié)合混合模型中磁偶極子參數(shù)與艦船結(jié)構(gòu)的相關(guān)性, 以建模精度和模型穩(wěn)定性為目標(biāo)構(gòu)造多目標(biāo)函數(shù), 通過(guò)多目標(biāo)函數(shù)優(yōu)化獲得合理的磁偶極子參數(shù), 間接地將建模求解問(wèn)題轉(zhuǎn)化為多目標(biāo)優(yōu)化問(wèn)題, 利用多目標(biāo)函數(shù)優(yōu)化算法對(duì)多目標(biāo)問(wèn)題進(jìn)行求解, 提出一種高精度穩(wěn)定艦船磁場(chǎng)建模方法.

    2 艦船磁場(chǎng)混合模型和多目標(biāo)函數(shù)構(gòu)造

    2.1 艦船磁場(chǎng)混合模型

    通常, 艦船磁場(chǎng)混合模型是指將艦船等效為一個(gè)均勻磁化旋轉(zhuǎn)橢球體和磁偶極子陣列的混合模型[1,16,17], 如圖1所示, 建立艦船磁場(chǎng)混合模型時(shí),均勻橢球體位于艦船中心, 其長(zhǎng)軸等于船長(zhǎng), 短軸等于船寬, 用于擬合艦船的宏觀磁場(chǎng), 磁偶極子陣列均勻分布于艦船吃水線上, 模擬艦船的局部不均勻磁場(chǎng).

    圖1 艦船磁場(chǎng)混合模型Fig.1.Ship magnetic field mixing model.

    假設(shè)艦船混合模型的磁偶極子個(gè)數(shù)為N, 坐標(biāo)分別為(αi, βi, γi), i =1,2,···,N .如圖1所示,以艦船中心建立坐標(biāo)系, 則在測(cè)量點(diǎn) Pj(xj,yj,zj)產(chǎn)生的磁場(chǎng)為[1]

    式中, Mxi, Myi, Mzi(i = 1, 2, ···, N)分 別 為 第i個(gè)磁偶極子x, y, z方向的磁矩, Mx(N + 1), My(N + 1),Mz(N + 1)為均勻旋轉(zhuǎn)橢球體在x, y, z方向磁矩,(axij, ayij, azij, bxij, byij, bzij, cxij, cyij, czij)分別為對(duì)應(yīng)的磁場(chǎng)計(jì)算系數(shù), 具體計(jì)算公式如下:

    1)當(dāng)i = 1, 2, ···, N時(shí), 為磁偶極子對(duì)應(yīng)的系數(shù)[1]:

    2) 當(dāng)i = N + 1時(shí), 為均勻橢球體對(duì)應(yīng)的系數(shù)[1]:

    將(1)式簡(jiǎn)寫(xiě)為矩陣形式:

    2.2 艦船磁場(chǎng)混合模型建模求解

    利用艦船混合模型建模時(shí), 通常測(cè)量艦船某一深度平面上的 n ×m 個(gè)測(cè)量點(diǎn)數(shù)據(jù), 如圖2所示.根據(jù)(1)式和(2)式可構(gòu)造求解方程組[1]:

    圖2 磁場(chǎng)測(cè)量Fig.2.Magnetic field measurement.

    式中, F = (F1, F1, ···, Fn × m)T, H = (H1, H1, ···,Hn × m)T.在方程組(3)式中, 測(cè)量點(diǎn)坐標(biāo)(xj, yj,zj)和測(cè)量點(diǎn)磁場(chǎng)(Hxj, Hyj, Hzj)T是已知的, 若混合模型中的磁偶極子均勻排布在吃水線上, 則可由艦船的長(zhǎng)度L和磁偶極子個(gè)數(shù)N計(jì)算得到磁偶極子坐標(biāo)(αi, βi, γi).此時(shí), (3)式為一線性方程, 磁矩參數(shù)M為唯一待求參數(shù), 若能準(zhǔn)確求解出磁矩參數(shù)M, 則可以根據(jù)(1)式計(jì)算艦船磁場(chǎng)空間分布,即完成艦船的磁場(chǎng)建模.

    若想求解方程組(3)式中磁矩參數(shù)M, 測(cè)量點(diǎn)數(shù)必須滿(mǎn)足 n ×m≥N+1 .通常, 為了獲得更為精確的艦船磁場(chǎng)模型, 測(cè)量點(diǎn)數(shù) n ×m 往往需要遠(yuǎn)大于等效源個(gè)數(shù) N +1 , 這使得方程組(3)式為方程階數(shù)很大的超定方程, 且其病態(tài)程度有時(shí)也會(huì)因測(cè)量點(diǎn)位置變得很大.因此, 對(duì)方程組(3)式的求解通常使用逐步回歸進(jìn)行求解.

    上述為使用混合模型進(jìn)行建模和求解的相關(guān)原理和方法, 其中艦船混合模型中的磁偶極子是均勻分布于艦船的吃水線上的.

    2.3 多目標(biāo)函數(shù)構(gòu)造

    2.2 節(jié)中使用混合模型進(jìn)行建模時(shí), 是將磁偶極子均勻分布于艦船的吃水線上, 磁偶極子的位置坐標(biāo)(αi, βi, γi)是已知量.然而, 混合模型中的磁偶極子代表艦船的局部磁場(chǎng)特性, 其數(shù)量和位置的分布與艦船內(nèi)部結(jié)構(gòu)密切相關(guān), 方程組(3)式實(shí)際上是對(duì)艦船磁場(chǎng)的一種擬合, 當(dāng)磁偶極子個(gè)數(shù)過(guò)多或者位置分布不正確時(shí), 對(duì)方程組(3)式的求解常常出現(xiàn)精度不高、穩(wěn)定性差、外推能力弱等問(wèn)題.

    為了建立更為合理的艦船磁場(chǎng)模型, 假設(shè)混合模型中的磁偶極子個(gè)數(shù)N和位置分布(αi, βi, γi)為未知量, 則由(1)式可知, 方程組(3)式中的系數(shù)矩陣F為磁偶極子個(gè)數(shù)N及其位置分布(αi, βi,γi)的函數(shù), 可通過(guò)優(yōu)化(4)式的相對(duì)誤差獲得磁偶極子個(gè)數(shù)及其分布, 建優(yōu)化立目標(biāo)函數(shù)f1如下:

    由前面2.2節(jié)可知, 方程組(3)式往往是一個(gè)病態(tài)的超定方程, 對(duì)輸入的敏感性極強(qiáng), 在解方程(4)過(guò)程中, 輕微的測(cè)量點(diǎn)數(shù)據(jù)擾動(dòng)都會(huì)嚴(yán)重影響模型求解結(jié)果.為了提高求解磁場(chǎng)模型的精度和穩(wěn)定性, 通常的做法是降低系數(shù)矩陣F的條件數(shù), 系數(shù)矩陣的條件數(shù)越小, 模型越穩(wěn)定, 對(duì)測(cè)量誤差越不敏感.由于系數(shù)矩陣F是磁偶極子個(gè)數(shù)及其位置坐標(biāo)的函數(shù), 因此, 可以引入系數(shù)矩陣條件數(shù)作為優(yōu)化目標(biāo)函數(shù), 通過(guò)對(duì)系數(shù)矩陣條件數(shù)的優(yōu)化而得到合理的磁偶極個(gè)數(shù)及其分布, 以系數(shù)矩陣條件數(shù)構(gòu)造的優(yōu)化目標(biāo)函數(shù)f2表示如下:

    式中, C ond(·) 為條件數(shù)計(jì)算.

    因此, 由上面引入的兩個(gè)優(yōu)化目標(biāo)函數(shù) f1和 f2,可構(gòu)造多目標(biāo)優(yōu)化函數(shù)如下:

    對(duì)多目標(biāo)函數(shù)(7)式進(jìn)行優(yōu)化求解, 獲得磁偶極子個(gè)數(shù)及其分布, 由于多目標(biāo)函數(shù)(7)式兼顧了建模精度和模型穩(wěn)定, 所求磁偶極子個(gè)數(shù)及其分布能夠保證一定模型精度下具備較好的模型穩(wěn)定性.同時(shí), 在目標(biāo)函數(shù)f1中的磁矩參數(shù)作為中間變量,間接地將建模問(wèn)題轉(zhuǎn)換為了對(duì)(7)式的多目標(biāo)優(yōu)化問(wèn)題.

    通常, 用于建模的磁偶極子個(gè)數(shù)N在2-15之間, 因此, (7)式中的變量N的取值約束為N ∈Ζ, N∈[2, 15].磁偶極子位置(αi, βi, γi)可由艦船的尺寸進(jìn)行約束, 如圖3所示.在混合模型中, 艦船等效橢球體, 長(zhǎng)軸為L(zhǎng)/2, 短軸為W/2, N個(gè)磁偶極子分布在橢球內(nèi), 將橢球體沿長(zhǎng)軸方向分為N等份, 則第i個(gè)磁偶極子的分布范圍為圖3中陰影部分所代表的區(qū)域, 即(7)式中, 磁偶極子位置(αi, βi, γi)分布的取值范圍為

    圖3 混合模型中磁偶極子分布范圍Fig.3.Distribution range of magnetic dipoles in the mixed model.

    式 中, xi=-L/2+(i-1)L/N , xi+1=-L/2-iL/N , i = 1, 2, ···, N + 1.

    3 基于多目標(biāo)粒子群優(yōu)化的建模求解算法

    由2.3節(jié)可知, 艦船磁場(chǎng)的高精度穩(wěn)定建??蓺w結(jié)為對(duì)(7)式多目標(biāo)函數(shù)的尋優(yōu)問(wèn)題, 聯(lián)合(8)式, 該多目標(biāo)優(yōu)化問(wèn)題可表示為:

    由(9)式可知, 待優(yōu)化參數(shù)的維數(shù)為3N + 1,與磁偶極子個(gè)數(shù)相關(guān), 且目標(biāo)函數(shù)是待優(yōu)化參數(shù)的非線性函數(shù), 對(duì)于多維度非線性目標(biāo)函數(shù)優(yōu)化問(wèn)題, 待求參數(shù)維數(shù)越高, 優(yōu)化求解越困難.混合模型中的磁偶極子個(gè)數(shù)的取值只能為整數(shù), 其取值范圍已經(jīng)給出, 故可將(9)式中的待求參數(shù)N作為已知量, 按照磁偶極子個(gè)數(shù)參數(shù)取值范圍依次取值,然后再進(jìn)行優(yōu)化求解, (9)式改寫(xiě)為

    式中的N為假設(shè)已知量.(10)式是關(guān)于磁偶極子位置分布的多目標(biāo)優(yōu)化問(wèn)題, 下面將研究該問(wèn)題的求解算法.

    3.1 多目標(biāo)粒子群優(yōu)化算法

    傳統(tǒng)多目標(biāo)優(yōu)化算法是將(10)式中的各目標(biāo)函數(shù)取加權(quán)和, 從而多目標(biāo)優(yōu)化問(wèn)題轉(zhuǎn)化為單目標(biāo)優(yōu)化問(wèn)題, 由前面可知, 目標(biāo)函數(shù) f1為建模相對(duì)誤差, f2為模型系數(shù)矩陣條件數(shù), 各目標(biāo)函數(shù)不具備可比性, 加權(quán)系數(shù)難于選取.多目標(biāo)問(wèn)題中, 各目標(biāo)之間往往是相互沖突的, 通常通過(guò)決策變量來(lái)協(xié)調(diào)權(quán)衡和折中處理, 使各子目標(biāo)函數(shù)盡可能達(dá)到最優(yōu).多目標(biāo)優(yōu)化問(wèn)題不存在單個(gè)最優(yōu)解, 而是Pareto最優(yōu)解集, 如何求得與真實(shí)Pareto前沿一致的Pareto最優(yōu)解集是多目標(biāo)函數(shù)優(yōu)化問(wèn)題的關(guān)鍵.

    目前, 通常用于解決多目標(biāo)函數(shù)優(yōu)化問(wèn)題的算法包括: 多目標(biāo)進(jìn)化算法(MOEA)[29]、多目標(biāo)蟻群算法(MASA)[30]、多目標(biāo)模擬退火算法(MSA)[31]、多目標(biāo)粒子群算法(MPSO)[32]等.MOEA對(duì)于低維優(yōu)化問(wèn)題具有較好的表現(xiàn), 在高維問(wèn)題中, 解集的多樣性不理想且耗時(shí)較長(zhǎng), 運(yùn)行效率不高;MASA在高維的多目標(biāo)函數(shù)優(yōu)化問(wèn)題上同樣需要較長(zhǎng)的搜索時(shí)間, 且易于出現(xiàn)早熟停滯現(xiàn)象; MSA對(duì)搜索空間(目標(biāo)函數(shù)的性質(zhì))不加任何限制, 可以是不連續(xù)的、不可微的, 并且也能求得Pareto邊界上多個(gè)不同方向的Pareto最優(yōu)解, 但是需要大量的迭代次數(shù), 因而收斂速度慢、優(yōu)化效率較低;MPSO具有高效的搜索能力, 有利于得到多目標(biāo)意義下的最優(yōu)解, 同時(shí), 其通用性比較好, 實(shí)現(xiàn)也較為容易, 適合處理多種類(lèi)型的目標(biāo)函數(shù)和約束,并且容易與傳統(tǒng)的優(yōu)化方法結(jié)合, 從而改進(jìn)自身的局限性, 更高效地解決問(wèn)題[33].本文所構(gòu)造的多目標(biāo)優(yōu)化問(wèn)題是一個(gè)高維的、復(fù)雜非線性問(wèn)題, 考慮到求解算法的效率和解的質(zhì)量, 本文將采用MPSO進(jìn)行Pareto最優(yōu)解集求解.

    多目標(biāo)粒子群優(yōu)化算法是在粒子群優(yōu)化算法(PSO)的基礎(chǔ)上進(jìn)行改進(jìn)的, 粒子群算法是由Kennedy[34]提出的基于群智能的優(yōu)化算法, 標(biāo)準(zhǔn)的粒子速度和位置更新公式如下:

    MPSO對(duì)PSO增加了外部?jī)?chǔ)備集保存算法迭代過(guò)程中出現(xiàn)的Pareto最優(yōu)解, 采取個(gè)體最優(yōu)全局最優(yōu) gt和外部?jī)?chǔ)備集的更新策略, 為防止算法早熟收斂, 增加種群多樣性策略.本文采用文獻(xiàn)[32]中的自適應(yīng)網(wǎng)格法的多目標(biāo)粒子群算法, 該算法根據(jù)Pareto支配關(guān)系更新個(gè)體最優(yōu); 采用自適應(yīng)網(wǎng)格法計(jì)算粒子的擁擠程度, 并從擁擠程度小的網(wǎng)格中隨機(jī)選擇一個(gè)粒子作為群體最優(yōu)粒子gt; 采用Pareto占優(yōu)準(zhǔn)則和自適應(yīng)網(wǎng)格法對(duì)外部?jī)?chǔ)備集進(jìn)行更新.

    自適應(yīng)網(wǎng)格法的多目標(biāo)粒子群算法的具體步驟如下:

    1) 輸入?yún)?shù): w, c1, c2, Niter(迭代次數(shù)), Nn(種群數(shù)), Nr(外部?jī)?chǔ)備集數(shù)), xrang(變量范圍),vrang(粒子速度范圍), Ngri(網(wǎng)格數(shù));

    2) 初始化種群粒子, 隨機(jī)選取可行域中的粒子初始位置和初始速度;

    3) 計(jì)算各粒子的適應(yīng)度值, 根據(jù)支配關(guān)系形成非支配解集;

    4) 更新外部?jī)?chǔ)備集;

    5) 更新個(gè)體最優(yōu)和全局最優(yōu);

    6) 按照(11)式和(12)式更新速度和粒子位置;

    7) 根據(jù)速度和位置的取值范圍調(diào)整速度和粒子位置的取值;

    8) 根據(jù)種群多樣性策略, 對(duì)粒子進(jìn)行多樣性處理;

    9) 檢查是否達(dá)到最大迭代次數(shù), 若是, 則終止程序; 如果未達(dá)到, 則繼續(xù)返回第二步.

    上述步驟中, 多目標(biāo)粒子群算法參數(shù)w, c1, c2,Niter, Nn, N r , Ngri可按常用的取值進(jìn)行設(shè)置, 粒子的取值范圍 xrang=[xmin,xmax] 可由(8)式計(jì)算,粒子速度范圍可設(shè)置為 vrang=[-0.4(xmax-xmin),0.4(xmax-xmin)] , 在算法初始化時(shí), 粒子取值和粒子速度取值在相應(yīng)的取值范圍內(nèi)選取.

    3.2 解的選取

    3.1 節(jié)中給出了求解多目標(biāo)優(yōu)化問(wèn)題(10)式的算法, 利用該算法能夠求解出優(yōu)化問(wèn)題(10)式的Pareto最優(yōu)解集, 磁偶極子位置的最優(yōu)分布可以從Pareto最優(yōu)解集中選取, 同時(shí), (10)式的優(yōu)化問(wèn)題是磁偶極子個(gè)數(shù)假設(shè)已知下的簡(jiǎn)化問(wèn)題, 而由(9)式可知, 磁偶極子個(gè)數(shù)也是需要優(yōu)化求解的.下面, 將對(duì)如何選取磁偶極子個(gè)數(shù)和如何從Pareto最優(yōu)解集中選取最佳的解進(jìn)行研究, 給出相應(yīng)的選擇規(guī)則.

    由前面可知, 磁偶極子個(gè)數(shù)的選取范圍為N∈[2,15] , 且都為整數(shù), 在求解(9)式時(shí), 依次選取 N =2,3,···15 , 則(9)式可簡(jiǎn)化為(10)式, 利用3.1節(jié)中的多目標(biāo)粒子群算法求解得到Pareto最優(yōu)解集 ΨN.所有磁偶極子個(gè)數(shù)下求得的Pareto最優(yōu)解集構(gòu)成可選集 Ψ =Ψ2∩Ψ3∩···∩Ψ15.如何從可選集中選取最佳的解, 需要考慮建模的實(shí)際情況, 對(duì)于艦船磁場(chǎng)建模, 最為重要的是模型精度,在保證模型精度的情況下應(yīng)盡可能使模型穩(wěn)定, 故可以引入從Ψ中選取的限制條件:

    式中, x為磁偶極子位置分布 (α1,β1,γ1,α2,β2,γ2,···,αN,βN,γN) , e0為建模最大允許相對(duì)誤差.當(dāng)選取了滿(mǎn)足(13)式的x, 則其對(duì)應(yīng)的磁偶極子個(gè)數(shù)為所選的磁偶極子數(shù).

    (13)式給出了磁偶極子位置和磁偶極子個(gè)數(shù)選擇規(guī)則, 按照該選擇規(guī)則可以獲得相應(yīng)的優(yōu)化磁偶極子參數(shù), 用該磁偶極子參數(shù)聯(lián)合(1)式-(3)式可以求解出艦船磁場(chǎng)的磁矩參數(shù), 利用磁矩參數(shù)根據(jù)(1)式便能獲取艦船的空間磁場(chǎng)分布, 即完成了艦船磁場(chǎng)的建模.

    4 試驗(yàn)結(jié)果分析

    4.1 試驗(yàn)設(shè)置

    試驗(yàn)時(shí), 對(duì)三種艦船船模深度平面 Z0和 Z1上三條測(cè)線上的磁場(chǎng)進(jìn)行了測(cè)量, 如圖4所示.三種艦船船模按照長(zhǎng)寬尺寸分為小、中、大三種船型,小型艦船船模參數(shù): 長(zhǎng)L = 57.2 m、寬W = 8.6 m、吃水深度D = 2.2 m; 中型艦船船模參數(shù): 長(zhǎng)L =76.5 m、寬W = 8.5, 吃水深度D = 5.3 m; 大型艦船船模參數(shù): 長(zhǎng)L = 153 m、寬W = 17.3 m, 吃水深度D = 5.1 m.測(cè)量參數(shù)和艦船參數(shù)如表1所列.

    圖4 測(cè)量模式Fig.4.Measurement mode.

    表1 試驗(yàn)參數(shù)Table 1.Test parameters.

    對(duì)上述三種船模利用第3節(jié)的方法進(jìn)行建模,建模性能的評(píng)估包括建模精度和換算精度, 建模精度可以用建模的相對(duì)誤差來(lái)評(píng)估, 見(jiàn)(4)式; 換算是指用一個(gè)平面的測(cè)量數(shù)據(jù)建模, 利用所建模型換算另一個(gè)深度平面上的磁場(chǎng), 換算精度可以用換算相對(duì)誤差來(lái)評(píng)估, 計(jì)算公式如下:

    建模使用的多目標(biāo)粒子群優(yōu)化算法的參數(shù)設(shè)置: 種群數(shù)目200, 迭代次數(shù)500, 慣性權(quán)重0.49,學(xué)習(xí)因子 c1=c2=2 , 網(wǎng)格數(shù)為20, 外部?jī)?chǔ)備集規(guī)模200.磁偶極子個(gè)數(shù)N的取值設(shè)定為: 2, 4, 6, 8,10, 12, 14; (13)式中小、中、大三種艦船船模的解選擇參數(shù) e0分別為0.03, 0.03, 0.03.

    4.2 不同艦船船模實(shí)測(cè)數(shù)據(jù)建模

    4.2.1 建模有效性分析

    按照4.1的試驗(yàn)參數(shù), 利用表1中三種艦船船模Z0深度平面上的三條航跡測(cè)量數(shù)據(jù)進(jìn)行建模,得到三種艦船建模的可選集分布, 分布結(jié)果如圖5-圖7所示, 圖中不同顏色的解代表不同磁偶極子數(shù)時(shí)的Pareto最優(yōu)解集, 橫坐標(biāo) f2和縱坐標(biāo)f1對(duì)應(yīng)(7)式中的多目標(biāo)函數(shù)值, 分別對(duì)應(yīng)模型系數(shù)矩陣條件數(shù)和建模相對(duì)誤差.

    從圖5-圖7可以看出, 隨著磁偶極子個(gè)數(shù)的增多, 對(duì)應(yīng)的Pareto最優(yōu)集的建模擬合相對(duì)誤差f1整體減小, 條件數(shù) f2整體增大, 說(shuō)明混合模型中的磁偶極子數(shù)越多, 建模精度越高, 但與此同時(shí),條件數(shù)會(huì)增大, 模型會(huì)越不穩(wěn)定.因此, 建模精度和模型穩(wěn)定是相互矛盾的, 只能折中選擇.

    按照(13)式的選擇規(guī)則, 分別對(duì)圖5-圖7中的結(jié)果進(jìn)行選擇, 選擇的結(jié)果如表2所列, 表2給出了三種艦船的建模選擇結(jié)果, 包括磁偶極子個(gè)數(shù)、對(duì)應(yīng)系數(shù)矩陣條件數(shù)和對(duì)應(yīng)的建模相對(duì)誤差.從表2可知, 小型艦船選擇磁偶極子個(gè)數(shù)為10時(shí), 系數(shù)矩陣條件數(shù)為85.13, 相對(duì)誤差為0.0295;中型艦船選擇磁偶極子個(gè)數(shù)為8時(shí), 系數(shù)矩陣條件數(shù)為76.5, 相對(duì)誤差為0.0290; 大型艦船選擇磁偶極子個(gè)數(shù)為14時(shí), 系數(shù)矩陣條件數(shù)為128.95, 相對(duì)誤差為0.0256.

    圖5 小型艦船建模結(jié)果可選集分布Fig.5.Selectable distribution of modeling results for small ships.

    圖6 中型艦船建模結(jié)果可選集分布Fig.6.Selectable distribution of modeling results for medium-sized ships.

    圖7 大型艦船建模結(jié)果可選集分布Fig.7.Selectable distribution of modeling results for large ships.

    表2 Z0深度平面上的三種艦船船模建模結(jié)果Table 2.Modeling results of three kinds of ships on the Z0 depth plane.

    圖8-圖10為表2中結(jié)果所對(duì)應(yīng)的磁偶極子分布情況, 從圖8可知, 小型艦船的磁偶極子主要分布于艦船尾部、中間靠后下方和中間靠前上方,與實(shí)際的艦船船模結(jié)構(gòu)分布一致, 艦船尾部磁偶極子對(duì)應(yīng)于艦船尾部的螺旋槳、舵、錨等結(jié)構(gòu), 中間靠后下方的磁偶極子對(duì)應(yīng)的是機(jī)艙部分(發(fā)動(dòng)機(jī))結(jié)構(gòu), 中間靠前上方的磁偶極子分布對(duì)應(yīng)的是艦橋艙室部分.同樣, 對(duì)比圖9和圖10中的磁偶極子分布與相應(yīng)艦船船模結(jié)構(gòu), 發(fā)現(xiàn)所求磁偶極子的分布與船模的實(shí)際結(jié)構(gòu)較為一致, 說(shuō)明利用本文方法建模可以得到更加合理的磁偶極子分布.

    圖8 小型艦船船模磁偶極子分布情況Fig.8.Distribution of magnetic dipole of small ship model.

    圖9 中型艦船船模磁偶極子分布情況Fig.9.Distribution of magnetic dipole of medium ship model.

    圖10 大型艦船船模磁偶極子分布情況Fig.10.Distribution of magnetic dipoles of large ship models.

    分別對(duì)表1中三種艦船船模的兩種深度 Z0和Z1進(jìn)行建模和換算, 得到相應(yīng)的建模相對(duì)誤差和換算誤差如表3所列.可以看出, 三種艦船船模的建模相對(duì)誤差都小于3%, 由淺到深的換算誤差低于5%, 由深換算到淺的換算誤差低于6%.

    表3 不同深度的建模相對(duì)誤差和換算相對(duì)誤差Table 3.Modeling relative errors and converted relative errors of different depths.

    上文對(duì)三種艦船船模的實(shí)測(cè)數(shù)據(jù)進(jìn)行了建模,從建模結(jié)果可知, 本文方法能夠有效對(duì)三種艦船船模進(jìn)行建模, 建模相對(duì)誤差低于3%, 由淺到深的換算誤差低于5%, 由深換算到淺的換算誤差低于6%, 求解的磁偶極位置能夠較好與真實(shí)艦船的結(jié)構(gòu)相符合, 驗(yàn)證了本文方法能夠有效對(duì)艦船磁場(chǎng)進(jìn)行高精度建模.

    4.2.2 不同方法建模結(jié)果對(duì)比

    下面分別用本文方法、文獻(xiàn)[1]、文獻(xiàn)[25]中的方法進(jìn)行建模對(duì)比, 文獻(xiàn)[1]使用磁偶極子均勻分布的混合模型進(jìn)行建模, 文獻(xiàn)[25]對(duì)混合模型建模進(jìn)行了改進(jìn), 為了提高建模穩(wěn)定性, 利用模擬退火算法優(yōu)化系數(shù)矩陣的條件數(shù), 以獲得磁偶極子分布.文獻(xiàn)[1]和文獻(xiàn)[25]中沒(méi)有給定磁偶極子個(gè)數(shù)的取值, 下面使用這兩種方法建模時(shí), 選擇與本文方法相同的磁偶極子個(gè)數(shù).

    利用不同的建模方法對(duì)表1中三種艦船船模的不同深度實(shí)測(cè)數(shù)據(jù)進(jìn)行建模和換算, 計(jì)算結(jié)果如表4所列, 表中“→”表示從一個(gè)深度換算到另一個(gè)深度, 從表4可知, 相對(duì)其他兩種方法, 文獻(xiàn)[1]中的磁偶極子均勻分布的混合模型系數(shù)矩陣條件數(shù)是最大的, 其建模相對(duì)誤差和換算相對(duì)誤差, 比文獻(xiàn)[25]小, 比本文方法大; 文獻(xiàn)[25]的方法系數(shù)矩陣條件數(shù)最小, 建模相對(duì)誤差和換算相對(duì)誤差最大, 導(dǎo)致這種結(jié)果原因是其僅僅通過(guò)優(yōu)化系數(shù)矩陣條件數(shù)來(lái)優(yōu)化磁偶極子分布, 雖然可以獲得最小的系數(shù)矩陣條件數(shù), 但忽略了磁偶極子位置與艦船結(jié)構(gòu)的關(guān)系, 使得模型與艦船實(shí)際模型不符, 導(dǎo)致建模相對(duì)誤差和換算誤差最大, 建模性能最差; 本文方法所建模型的系數(shù)矩陣條件數(shù)遠(yuǎn)遠(yuǎn)小于文獻(xiàn)[1],稍微大于文獻(xiàn)[25], 模型相對(duì)誤差最小, 換算誤差也是最小, 說(shuō)明本文方法建模性能要好于其他兩種方法.

    表4 不同方法的建模相對(duì)誤差和換算相對(duì)誤差Table 4.Modeling relative error and conversion relative error of different methods.

    4.2.3 建模穩(wěn)定性分析

    當(dāng)使用混合模型對(duì)艦船磁場(chǎng)進(jìn)行建模時(shí), 其本質(zhì)是求解一個(gè)病態(tài)超定方程, 病態(tài)超定方程對(duì)輸入極其敏感, 輸入有輕微的變化都會(huì)導(dǎo)致求解結(jié)果偏差很大, 而對(duì)于一個(gè)艦船磁場(chǎng)模型, 其磁矩參數(shù)應(yīng)當(dāng)是穩(wěn)定不變的, 由于測(cè)量誤差的存在, 常常會(huì)導(dǎo)致模型的求解結(jié)果(磁矩參數(shù))不穩(wěn)定.為了驗(yàn)證本文方法建模的穩(wěn)定性, 分別計(jì)算存在測(cè)量干擾下的建模精度和磁矩, 對(duì)比存在和不存在測(cè)量干擾下的建模精度和磁矩變化.

    現(xiàn)隨機(jī)選取三種艦船 Z0深度平面上的30個(gè)測(cè)量數(shù)據(jù)點(diǎn)增加20%的擾動(dòng), 利用本文方法、文獻(xiàn)[1]和文獻(xiàn)[25]中的方法進(jìn)行建模, 計(jì)算各種方法的建模相對(duì)誤差, 結(jié)果如表5所列.從表5中的計(jì)算結(jié)果可知, 文本方法在有干擾的情況下, 建模相對(duì)誤差不大于0.04, 相對(duì)于無(wú)干擾時(shí), 相對(duì)誤差波動(dòng)在0.01之內(nèi), 驗(yàn)證了本文方法對(duì)有干擾存在時(shí)的建模有效性; 對(duì)比其他兩種方法, 本文方法無(wú)論在有干擾和無(wú)干擾情況下, 建模相對(duì)誤差都是最小.文獻(xiàn)[25]在有干擾和無(wú)干擾下, 誤差波動(dòng)較小,是由于其通過(guò)優(yōu)化系數(shù)矩陣條件數(shù)提高了建模穩(wěn)定性, 但是, 由于忽略了磁偶極子分布與艦船結(jié)構(gòu)的關(guān)系, 使其建模相對(duì)誤差較大, 與前面4.2.2節(jié)中的不同深度上建模結(jié)果一致.綜上可知, 本文方法在測(cè)量數(shù)據(jù)存在干擾時(shí), 具有較高的建模精度, 能夠有效建模, 且相對(duì)其他兩種方法, 建模性能最好.

    表5 干擾下不同方法的建模相對(duì)誤差和換算誤差Table 5.Modeling relative errors and conversion errors of different methods under interference.

    圖11-圖13為不同艦船使用不同方法在有干擾和無(wú)干擾時(shí)磁矩求解的絕對(duì)誤差, 絕對(duì)誤差計(jì)算公式為

    圖11 干擾下不同方法求解的小型艦船磁矩絕對(duì)誤差 (a) x方向磁矩絕對(duì)誤差 Δ mx ; (b) y方向磁矩絕對(duì)誤差 Δ my ; (c) z方向磁矩絕對(duì)誤差 ΔmzFig.11.Absolute errors of magnetic moment of small ships solved by different methods under disturbance: (a) Absolute errors of magnetic moment in x-direction; (b) absolute errors of magnetic moment in y-direction; (c) absolute errors of magnetic moment in z-direction.

    圖13 干擾下不同方法求解的大型艦船磁矩絕對(duì)誤差 (a) x方向磁矩絕對(duì)誤差 Δ mx ; (b) y方向磁矩絕對(duì)誤差 Δ my ; (c) z方向磁矩絕對(duì)誤差 ΔmzFig.13.Absolute errors of magnetic moment of large ships solved by different methods under disturbance: (a) Absolute errors of magnetic moment in x-direction; (b) absolute errors of magnetic moment in y-direction; (c) absolute errors of magnetic moment in z-direction.

    式中, i = 1-N為對(duì)應(yīng)磁偶極子的序號(hào), i = N + 1對(duì)應(yīng)橢球體, j = x, y, z代表不同方向的磁矩,為無(wú)干擾下建模求解的磁矩,為存在干擾下建模求解的磁矩.

    圖12 干擾下不同方法求解的中型艦船磁矩絕對(duì)誤差 (a) x方向磁矩絕對(duì)誤差 Δ mx ; (b) y方向磁矩絕對(duì)誤差 Δ my ; (c) z方向磁矩絕對(duì)誤差 ΔmzFig.12.Absolute error of magnetic moment of medium ship solved by different methods under disturbance: (a) Absolute errors of magnetic moment in x-direction; (b) absolute errors of magnetic moment in y-direction; (c) absolute errors of magnetic moment in z-direction.

    從圖11-圖13中的結(jié)果可知, 本文方法在有干擾和無(wú)干擾下, 求解的磁矩絕對(duì)誤差小于 1 02,相對(duì)于艦船磁矩量級(jí)[1]( 1 04-107)來(lái)說(shuō)是很小的,可認(rèn)為求解的磁矩在有無(wú)干擾時(shí)幾乎是一致的, 說(shuō)明本文方法具有較好的穩(wěn)定性, 即本文方法能夠?qū)ε灤艌?chǎng)進(jìn)行穩(wěn)定建模.同時(shí), 對(duì)比圖11-圖13中三種方法所求磁矩在有干擾和無(wú)干擾時(shí)的磁矩絕對(duì)誤差可知, 文獻(xiàn)[1]中的方法求解磁矩相差最大, 文獻(xiàn)[25]中的方法次之, 本文方法相差最小,驗(yàn)證了在三種方法中, 本文方法建模具有最好的穩(wěn)定性.

    上文對(duì)三種艦船船模在有干擾和無(wú)干擾下進(jìn)行了建模, 對(duì)比了三種方法的建模相對(duì)誤差和磁矩絕對(duì)誤差, 驗(yàn)證了本文方法具有良好的穩(wěn)定性, 且相對(duì)其他兩種方法, 穩(wěn)定性是最好的.

    4.3 真實(shí)海況實(shí)測(cè)數(shù)據(jù)建模結(jié)果

    為了進(jìn)一步驗(yàn)證本文建模方法的現(xiàn)實(shí)有效性,利用實(shí)測(cè)的某型艦船通過(guò)數(shù)據(jù)進(jìn)行建模.實(shí)際測(cè)量是在消磁站中進(jìn)行的, 在海底13和22.3 m兩種深度平面上布置了5個(gè)三軸磁傳感器, 艦船的航行方向與傳感器列陣垂直, 分別采集了不同航向上的實(shí)測(cè)數(shù)據(jù), 測(cè)量系統(tǒng)為該消磁站自主研制的艦船水下磁測(cè)系統(tǒng), 能夠較為精確的測(cè)量艦船的磁場(chǎng), 并且能夠給出測(cè)量點(diǎn)相對(duì)于艦船坐標(biāo)系下的位置.如圖為測(cè)量場(chǎng)景示意圖14.

    圖14 測(cè)量場(chǎng)景示意圖Fig.14.The schematic diagram of measurement scene.

    利用本文方法進(jìn)行了建模, 計(jì)算不同深度平面上的不同航向的建模相對(duì)誤差和換算誤差(13 m → 22.3 m), 計(jì)算結(jié)果如表6所列.可以看出, 本文方法在現(xiàn)實(shí)艦船磁場(chǎng)建模中, 相對(duì)誤差小于0.03, 換算誤差小于0.07, 說(shuō)明本文方法能夠有效在實(shí)際工程中應(yīng)用.

    表6 真實(shí)艦船測(cè)量數(shù)據(jù)建模結(jié)果Table 6.Modeling results of real ship measurement data.

    5 結(jié) 論

    本文圍繞如何對(duì)艦船磁場(chǎng)高精度穩(wěn)定建模進(jìn)行了研究, 在艦船磁場(chǎng)混合模型建模方法的基礎(chǔ)上, 充分考慮了建模系數(shù)矩陣條件數(shù)和建模精度,構(gòu)造多目標(biāo)函數(shù)優(yōu)化建模求解問(wèn)題, 間接將模型求解轉(zhuǎn)換為多目標(biāo)函數(shù)優(yōu)化求解, 提出了艦船磁場(chǎng)高精度穩(wěn)定模型建模方法.實(shí)測(cè)船模數(shù)據(jù)建模試驗(yàn)結(jié)果表明, 本文方法能夠有效進(jìn)行建模, 建模相對(duì)誤差小于0.03, 換算誤差小于0.06, 和原有的建模方法相比具有較高的精度和較好的磁場(chǎng)換算性能; 通過(guò)分析有無(wú)干擾條件下的建模結(jié)果, 驗(yàn)證了本文建模方法受干擾影響小, 具有很好的建模穩(wěn)定性.利用本文方法對(duì)海上實(shí)測(cè)數(shù)據(jù)建模進(jìn)行建模, 建模結(jié)果表明本文方法具有較高的建模精度和換算精度,能夠有效地在實(shí)際工程中應(yīng)用.

    猜你喜歡
    磁偶極子磁矩艦船
    艦船通信中的噪聲消除研究
    艦船測(cè)風(fēng)傳感器安裝位置數(shù)值仿真
    基于遞推更新卡爾曼濾波的磁偶極子目標(biāo)跟蹤
    磁偶極子跟蹤的漸進(jìn)貝葉斯濾波方法
    艦船腐蝕預(yù)防與控制系統(tǒng)工程
    基于磁偶極子的磁場(chǎng)梯度張量縮并的試驗(yàn)驗(yàn)證及相關(guān)參數(shù)確定
    CoFeB/MgO磁隧道結(jié)的低電流密度磁矩翻轉(zhuǎn)特性
    水平磁偶極子電磁場(chǎng)特征研究
    兩種計(jì)算帶電輕子磁矩的嘗試
    河南科技(2014年23期)2014-02-27 14:18:52
    用微元法討論旋轉(zhuǎn)均勻帶電球的磁矩
    咕卡用的链子| 嫩草影院入口| 国产又爽黄色视频| 亚洲伊人色综图| 9热在线视频观看99| 女人爽到高潮嗷嗷叫在线视频| 一级片'在线观看视频| 精品国产乱码久久久久久男人| 欧美最新免费一区二区三区| 久久女婷五月综合色啪小说| 最近中文字幕2019免费版| 久久 成人 亚洲| 亚洲 欧美一区二区三区| 亚洲色图综合在线观看| 欧美日韩国产mv在线观看视频| 一本大道久久a久久精品| 五月天丁香电影| 日本色播在线视频| 黄片播放在线免费| 一级片'在线观看视频| 晚上一个人看的免费电影| 一区在线观看完整版| 免费日韩欧美在线观看| 高清黄色对白视频在线免费看| 亚洲欧美精品综合一区二区三区| 午夜av观看不卡| 久久久久久人妻| 三上悠亚av全集在线观看| 热re99久久精品国产66热6| 欧美国产精品va在线观看不卡| 妹子高潮喷水视频| 尾随美女入室| 丝袜人妻中文字幕| 丰满乱子伦码专区| 国产精品99久久99久久久不卡 | 国产亚洲精品第一综合不卡| 丝瓜视频免费看黄片| 欧美日本中文国产一区发布| 久久久久视频综合| 国产有黄有色有爽视频| 丝袜喷水一区| 咕卡用的链子| 免费久久久久久久精品成人欧美视频| 纵有疾风起免费观看全集完整版| 91成人精品电影| 亚洲成人手机| 国产男女内射视频| 1024视频免费在线观看| 亚洲国产毛片av蜜桃av| 午夜91福利影院| 香蕉国产在线看| 国产成人免费观看mmmm| 欧美xxⅹ黑人| 国产成人a∨麻豆精品| 国产精品 欧美亚洲| 亚洲国产av影院在线观看| 1024视频免费在线观看| 99国产精品免费福利视频| 国产淫语在线视频| 亚洲精品成人av观看孕妇| 两个人看的免费小视频| 日韩不卡一区二区三区视频在线| 国产男女超爽视频在线观看| 九九爱精品视频在线观看| 亚洲欧美清纯卡通| 日本vs欧美在线观看视频| 亚洲七黄色美女视频| 免费少妇av软件| 曰老女人黄片| 男女高潮啪啪啪动态图| 大香蕉久久成人网| 国产毛片在线视频| www.自偷自拍.com| 亚洲av综合色区一区| 亚洲国产精品成人久久小说| 一级毛片我不卡| 久久影院123| 少妇的丰满在线观看| 久久免费观看电影| 51午夜福利影视在线观看| www.av在线官网国产| 18禁国产床啪视频网站| 免费少妇av软件| 三上悠亚av全集在线观看| 九色亚洲精品在线播放| 无限看片的www在线观看| 久久99热这里只频精品6学生| 成人国产av品久久久| 女人爽到高潮嗷嗷叫在线视频| 青春草亚洲视频在线观看| 精品一区在线观看国产| 制服丝袜香蕉在线| 国产成人精品在线电影| 黑人巨大精品欧美一区二区蜜桃| 精品少妇黑人巨大在线播放| 十八禁高潮呻吟视频| 亚洲欧美激情在线| 国产视频首页在线观看| 久久久精品国产亚洲av高清涩受| 亚洲精品日韩在线中文字幕| 国产精品久久久人人做人人爽| 久久人人爽人人片av| 看免费av毛片| 熟女少妇亚洲综合色aaa.| 亚洲人成电影观看| 美女主播在线视频| 男人操女人黄网站| 欧美精品一区二区大全| 51午夜福利影视在线观看| 黑丝袜美女国产一区| 毛片一级片免费看久久久久| 中文天堂在线官网| 国产精品国产三级专区第一集| 亚洲欧美成人综合另类久久久| 精品久久蜜臀av无| 一级毛片电影观看| 中文字幕高清在线视频| 午夜激情av网站| 久久久久久久大尺度免费视频| 日韩精品免费视频一区二区三区| 久久精品国产亚洲av涩爱| av线在线观看网站| 国产黄频视频在线观看| 成人亚洲欧美一区二区av| 亚洲成人免费av在线播放| 久热这里只有精品99| 一级爰片在线观看| 国产午夜精品一二区理论片| 99国产精品免费福利视频| 欧美黄色片欧美黄色片| 99国产精品免费福利视频| 久久天躁狠狠躁夜夜2o2o | 久久99一区二区三区| 久久久久久久久久久久大奶| 看非洲黑人一级黄片| 亚洲色图 男人天堂 中文字幕| 自拍欧美九色日韩亚洲蝌蚪91| 国产激情久久老熟女| 女人精品久久久久毛片| 久久这里只有精品19| 在线看a的网站| 夫妻性生交免费视频一级片| 深夜精品福利| 亚洲精品一二三| 黄片小视频在线播放| 18禁动态无遮挡网站| 国产精品久久久久成人av| 69精品国产乱码久久久| 19禁男女啪啪无遮挡网站| 十八禁高潮呻吟视频| 麻豆精品久久久久久蜜桃| 老司机亚洲免费影院| 男男h啪啪无遮挡| 亚洲精品日本国产第一区| 热re99久久精品国产66热6| 国产亚洲一区二区精品| 在线观看人妻少妇| 国产在线视频一区二区| 精品国产露脸久久av麻豆| 男女边摸边吃奶| 国产一区二区三区av在线| e午夜精品久久久久久久| 久久ye,这里只有精品| 狂野欧美激情性xxxx| 啦啦啦中文免费视频观看日本| 一本大道久久a久久精品| 中文字幕亚洲精品专区| 亚洲欧美成人精品一区二区| 99国产精品免费福利视频| 国产福利在线免费观看视频| 国产极品粉嫩免费观看在线| 国产激情久久老熟女| 最近手机中文字幕大全| 飞空精品影院首页| 久久久精品免费免费高清| 精品国产一区二区三区久久久樱花| 老鸭窝网址在线观看| 多毛熟女@视频| 又黄又粗又硬又大视频| 在线亚洲精品国产二区图片欧美| 午夜福利视频在线观看免费| 久久久久精品人妻al黑| 国产免费福利视频在线观看| 18在线观看网站| 精品久久蜜臀av无| 黄色视频在线播放观看不卡| 亚洲图色成人| 亚洲,欧美,日韩| 日韩大码丰满熟妇| 久久久久视频综合| 日本爱情动作片www.在线观看| 在线免费观看不下载黄p国产| 亚洲美女搞黄在线观看| kizo精华| 免费人妻精品一区二区三区视频| 国产97色在线日韩免费| 日日摸夜夜添夜夜爱| 午夜激情av网站| 亚洲精品久久成人aⅴ小说| 国产高清不卡午夜福利| 日韩av在线免费看完整版不卡| 国产97色在线日韩免费| 国产一区二区 视频在线| 啦啦啦啦在线视频资源| 99国产精品免费福利视频| 午夜av观看不卡| 中国国产av一级| 亚洲国产成人一精品久久久| 午夜免费男女啪啪视频观看| 免费高清在线观看日韩| 女人精品久久久久毛片| 日本欧美国产在线视频| 在线观看www视频免费| 熟女av电影| 一级,二级,三级黄色视频| 日韩一卡2卡3卡4卡2021年| 青春草亚洲视频在线观看| 黄色视频在线播放观看不卡| 老熟女久久久| 久久毛片免费看一区二区三区| 精品国产一区二区三区久久久樱花| 亚洲国产av影院在线观看| 超碰成人久久| av线在线观看网站| 国产探花极品一区二区| 男女高潮啪啪啪动态图| 黑人猛操日本美女一级片| 天美传媒精品一区二区| 国产精品一二三区在线看| 男女床上黄色一级片免费看| 香蕉国产在线看| 9色porny在线观看| 人妻人人澡人人爽人人| 国产在线一区二区三区精| 男女床上黄色一级片免费看| 日韩av在线免费看完整版不卡| 一区二区av电影网| 亚洲欧洲精品一区二区精品久久久 | 天天躁日日躁夜夜躁夜夜| 老司机靠b影院| 十八禁网站网址无遮挡| 久久这里只有精品19| 欧美 亚洲 国产 日韩一| 自线自在国产av| 国产日韩欧美在线精品| 99国产精品免费福利视频| 亚洲国产av新网站| 精品久久久精品久久久| 在线观看一区二区三区激情| 黄频高清免费视频| 久久久精品国产亚洲av高清涩受| 亚洲一级一片aⅴ在线观看| 亚洲综合精品二区| 一级爰片在线观看| 国产成人精品久久久久久| 最近最新中文字幕免费大全7| 天美传媒精品一区二区| 欧美人与性动交α欧美软件| 人人妻,人人澡人人爽秒播 | 秋霞在线观看毛片| 久久精品熟女亚洲av麻豆精品| 欧美日韩一区二区视频在线观看视频在线| 国产黄频视频在线观看| 大片电影免费在线观看免费| 美女国产高潮福利片在线看| 巨乳人妻的诱惑在线观看| 亚洲av日韩在线播放| 9色porny在线观看| 晚上一个人看的免费电影| 一本一本久久a久久精品综合妖精| 一级片免费观看大全| 男女国产视频网站| 纯流量卡能插随身wifi吗| 久久久久久久久久久免费av| 可以免费在线观看a视频的电影网站 | 男人舔女人的私密视频| 午夜福利在线免费观看网站| 久久精品久久久久久噜噜老黄| 国产av国产精品国产| 纯流量卡能插随身wifi吗| 制服丝袜香蕉在线| 国产一区二区 视频在线| 国产精品一区二区精品视频观看| 久久久久久久国产电影| 国产福利在线免费观看视频| 久久午夜综合久久蜜桃| 亚洲国产成人一精品久久久| 国产一区二区激情短视频 | 亚洲图色成人| 一本一本久久a久久精品综合妖精| 久久精品国产亚洲av高清一级| 亚洲色图 男人天堂 中文字幕| 日本爱情动作片www.在线观看| 新久久久久国产一级毛片| 大陆偷拍与自拍| 亚洲精品乱久久久久久| 亚洲欧美一区二区三区久久| 黑人巨大精品欧美一区二区蜜桃| 久久人妻熟女aⅴ| 在线观看三级黄色| 99九九在线精品视频| 在线观看人妻少妇| 天天躁夜夜躁狠狠久久av| 啦啦啦中文免费视频观看日本| 极品少妇高潮喷水抽搐| 久久亚洲国产成人精品v| 精品免费久久久久久久清纯 | 久久久欧美国产精品| 久久性视频一级片| 啦啦啦在线免费观看视频4| 男女午夜视频在线观看| 老鸭窝网址在线观看| 免费av中文字幕在线| 在线观看一区二区三区激情| 成人国产av品久久久| 1024视频免费在线观看| 国产黄色视频一区二区在线观看| 国产成人精品久久二区二区91 | 19禁男女啪啪无遮挡网站| 黑人欧美特级aaaaaa片| 欧美黑人欧美精品刺激| svipshipincom国产片| 国产一区有黄有色的免费视频| 日韩制服骚丝袜av| 十八禁高潮呻吟视频| 国产极品粉嫩免费观看在线| 久久人人97超碰香蕉20202| 亚洲美女黄色视频免费看| 欧美成人午夜精品| 又大又爽又粗| 99热全是精品| 最新的欧美精品一区二区| 国产精品99久久99久久久不卡 | 蜜桃国产av成人99| 免费观看a级毛片全部| 天堂俺去俺来也www色官网| 人成视频在线观看免费观看| 亚洲一区中文字幕在线| 满18在线观看网站| 中文乱码字字幕精品一区二区三区| 亚洲成人av在线免费| 日日爽夜夜爽网站| 午夜福利免费观看在线| 欧美日韩国产mv在线观看视频| 成年人免费黄色播放视频| 亚洲av电影在线观看一区二区三区| avwww免费| 久久久久视频综合| 男人添女人高潮全过程视频| 欧美中文综合在线视频| 国产又色又爽无遮挡免| 国产人伦9x9x在线观看| 蜜桃在线观看..| 亚洲欧洲精品一区二区精品久久久 | 两个人看的免费小视频| 国产精品二区激情视频| 在线观看国产h片| 一区二区日韩欧美中文字幕| 亚洲精品国产色婷婷电影| tube8黄色片| 成人国语在线视频| 亚洲欧美色中文字幕在线| 18禁裸乳无遮挡动漫免费视频| 欧美 日韩 精品 国产| 免费观看a级毛片全部| 国产一级毛片在线| 制服诱惑二区| 国产爽快片一区二区三区| 少妇被粗大猛烈的视频| 在线亚洲精品国产二区图片欧美| 一区二区av电影网| 日本91视频免费播放| 精品一区二区免费观看| 国产精品麻豆人妻色哟哟久久| 操美女的视频在线观看| 日本91视频免费播放| 18禁裸乳无遮挡动漫免费视频| 无遮挡黄片免费观看| 菩萨蛮人人尽说江南好唐韦庄| 嫩草影视91久久| 日韩精品有码人妻一区| av在线app专区| 桃花免费在线播放| 国产精品久久久av美女十八| av网站在线播放免费| 日本欧美视频一区| 国产激情久久老熟女| 精品人妻在线不人妻| 午夜久久久在线观看| 国产淫语在线视频| 高清av免费在线| 午夜福利视频在线观看免费| 美女扒开内裤让男人捅视频| 亚洲欧美成人精品一区二区| 久久99热这里只频精品6学生| 成年人免费黄色播放视频| 久久国产亚洲av麻豆专区| 大香蕉久久网| 91精品国产国语对白视频| 夫妻性生交免费视频一级片| 热99国产精品久久久久久7| 国产免费福利视频在线观看| 国产精品免费视频内射| 一区二区三区精品91| 丰满饥渴人妻一区二区三| 国产黄频视频在线观看| 久久久久久久国产电影| 桃花免费在线播放| 一区福利在线观看| 中文字幕最新亚洲高清| 久久鲁丝午夜福利片| 啦啦啦在线免费观看视频4| 国产成人精品在线电影| 国产精品熟女久久久久浪| 伊人亚洲综合成人网| 欧美精品亚洲一区二区| 香蕉国产在线看| 三上悠亚av全集在线观看| 国产一区二区在线观看av| 无遮挡黄片免费观看| 一级片'在线观看视频| 国产av一区二区精品久久| 亚洲精品国产色婷婷电影| 欧美人与善性xxx| 亚洲精品久久午夜乱码| 黑丝袜美女国产一区| 不卡av一区二区三区| 伊人亚洲综合成人网| 黄网站色视频无遮挡免费观看| 一区在线观看完整版| 日韩制服丝袜自拍偷拍| 日本vs欧美在线观看视频| 少妇 在线观看| 叶爱在线成人免费视频播放| 欧美久久黑人一区二区| 男人添女人高潮全过程视频| 性少妇av在线| 午夜精品国产一区二区电影| 精品福利永久在线观看| 王馨瑶露胸无遮挡在线观看| 人妻 亚洲 视频| 国产精品秋霞免费鲁丝片| 热re99久久精品国产66热6| 中文天堂在线官网| 无限看片的www在线观看| videosex国产| 91老司机精品| 在线观看免费视频网站a站| 成人国语在线视频| 老熟女久久久| 少妇被粗大猛烈的视频| 黄网站色视频无遮挡免费观看| 国产精品久久久久久人妻精品电影 | 黑人欧美特级aaaaaa片| 赤兔流量卡办理| 欧美日韩亚洲国产一区二区在线观看 | 成人黄色视频免费在线看| 免费黄色在线免费观看| 日韩伦理黄色片| 女人高潮潮喷娇喘18禁视频| xxx大片免费视频| kizo精华| 亚洲婷婷狠狠爱综合网| 欧美精品人与动牲交sv欧美| 国产精品久久久久久精品古装| 中文字幕人妻熟女乱码| 国产野战对白在线观看| av一本久久久久| 99国产精品免费福利视频| 国产成人a∨麻豆精品| 久久精品国产亚洲av高清一级| 秋霞在线观看毛片| 亚洲美女视频黄频| 亚洲精品久久成人aⅴ小说| 两个人看的免费小视频| 色网站视频免费| 97人妻天天添夜夜摸| 91国产中文字幕| 国产在线一区二区三区精| 一边亲一边摸免费视频| 精品人妻熟女毛片av久久网站| 毛片一级片免费看久久久久| av网站免费在线观看视频| 国产乱来视频区| 精品国产乱码久久久久久小说| 少妇精品久久久久久久| 久久久精品94久久精品| 亚洲精品美女久久av网站| 又黄又粗又硬又大视频| 久久精品亚洲av国产电影网| 精品亚洲成a人片在线观看| 欧美国产精品一级二级三级| 亚洲精品美女久久av网站| 高清欧美精品videossex| 日韩av不卡免费在线播放| 久久久久久人人人人人| 自拍欧美九色日韩亚洲蝌蚪91| 欧美黄色片欧美黄色片| 亚洲综合精品二区| 亚洲精品,欧美精品| 日日啪夜夜爽| 国产精品一区二区在线不卡| 国产男女超爽视频在线观看| 少妇人妻精品综合一区二区| 亚洲av成人精品一二三区| avwww免费| 亚洲色图综合在线观看| 午夜福利网站1000一区二区三区| 亚洲三区欧美一区| 最新在线观看一区二区三区 | 亚洲国产日韩一区二区| 国产精品嫩草影院av在线观看| 两性夫妻黄色片| 日日啪夜夜爽| 中国国产av一级| 色网站视频免费| 亚洲精品第二区| 日本wwww免费看| av视频免费观看在线观看| 久久97久久精品| 男女下面插进去视频免费观看| 一区二区三区激情视频| 成人国产av品久久久| 国产在线视频一区二区| 色婷婷久久久亚洲欧美| 婷婷色综合大香蕉| 久久国产精品男人的天堂亚洲| 丰满迷人的少妇在线观看| 国产成人精品在线电影| 99热全是精品| 一级毛片 在线播放| 亚洲成人国产一区在线观看 | 日日啪夜夜爽| 久久久久人妻精品一区果冻| 丝袜美腿诱惑在线| 妹子高潮喷水视频| 亚洲国产欧美在线一区| 看免费成人av毛片| 中文乱码字字幕精品一区二区三区| 国产熟女午夜一区二区三区| 国产xxxxx性猛交| 啦啦啦中文免费视频观看日本| 亚洲一码二码三码区别大吗| 亚洲精品国产av蜜桃| www.精华液| 90打野战视频偷拍视频| 一边摸一边做爽爽视频免费| 又大又黄又爽视频免费| 亚洲人成网站在线观看播放| 亚洲国产最新在线播放| 大香蕉久久成人网| 在线观看一区二区三区激情| av在线播放精品| 伊人亚洲综合成人网| 自线自在国产av| 国产视频首页在线观看| 欧美激情高清一区二区三区 | 精品一区二区三区四区五区乱码 | 90打野战视频偷拍视频| 一边摸一边做爽爽视频免费| 亚洲精品乱久久久久久| 亚洲男人天堂网一区| 久久精品亚洲熟妇少妇任你| 在线观看免费午夜福利视频| 亚洲第一av免费看| 中文字幕av电影在线播放| 精品福利永久在线观看| 国产在线一区二区三区精| 久久影院123| 天堂中文最新版在线下载| 女的被弄到高潮叫床怎么办| 2018国产大陆天天弄谢| 成人亚洲精品一区在线观看| 少妇 在线观看| 中国三级夫妇交换| 中文字幕人妻丝袜一区二区 | 黄片无遮挡物在线观看| 国产一区有黄有色的免费视频| av国产精品久久久久影院| 男女下面插进去视频免费观看| 亚洲美女搞黄在线观看| 成人18禁高潮啪啪吃奶动态图| 亚洲欧美精品自产自拍| 成年av动漫网址| 男女国产视频网站| 免费观看av网站的网址| 亚洲国产成人一精品久久久| 精品一区在线观看国产| 你懂的网址亚洲精品在线观看| 欧美 日韩 精品 国产| 欧美另类一区| 中文字幕人妻丝袜制服| 老司机靠b影院| 激情五月婷婷亚洲| 叶爱在线成人免费视频播放| 晚上一个人看的免费电影| 激情五月婷婷亚洲| 久久人妻熟女aⅴ| 国产精品蜜桃在线观看| www.精华液| 欧美最新免费一区二区三区| 热re99久久精品国产66热6| 国产男女内射视频| 99国产综合亚洲精品| 国产无遮挡羞羞视频在线观看| 日本猛色少妇xxxxx猛交久久| 成人亚洲欧美一区二区av| 美女视频免费永久观看网站| 丰满乱子伦码专区| 男女之事视频高清在线观看 | 精品一区二区三卡| 久久国产亚洲av麻豆专区| 一区二区三区精品91| 久久精品熟女亚洲av麻豆精品| 午夜91福利影院|