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

    基于SOA優(yōu)化的電磁層析成像圖像重建研究

    2021-01-07 04:56:20王亞?wèn)|霍繼偉趙鵬飛王成飛
    計(jì)算機(jī)測(cè)量與控制 2020年12期
    關(guān)鍵詞:病態(tài)搜索算法正則

    王亞?wèn)|,劉 澤,霍繼偉,趙鵬飛,王成飛,袁 偉

    (1.北京交通大學(xué) 電子信息工程學(xué)院,北京 100044;2.中國(guó)鐵道科學(xué)研究院集團(tuán)有限公司 標(biāo)準(zhǔn)計(jì)量研究所,北京 100081)

    0 引言

    電磁層析成像(EMT,electromagnetic tomography)是一種近二、三十年發(fā)展起來(lái)的基于電磁感應(yīng)原理的過(guò)程層析成像技術(shù)[1],其目的是研究具有電磁特性的物質(zhì)在空間的分布。典型的EMT系統(tǒng)主要由4個(gè)部分組成:傳感器陣列(即檢測(cè)激勵(lì)線圈)、激勵(lì)信號(hào)產(chǎn)生和分配單元、數(shù)據(jù)采集和調(diào)制單元以及圖像重建和信息獲取單元。通過(guò)對(duì)激勵(lì)線圈通入交變電流,在被測(cè)物場(chǎng)中產(chǎn)生激勵(lì)磁場(chǎng);當(dāng)被測(cè)物場(chǎng)內(nèi)有具有導(dǎo)磁或?qū)щ姷奈矬w放入時(shí),原有磁場(chǎng)分布將會(huì)改變,線圈之間的感性耦合也將隨之改變[2-3]。通過(guò)切換不同的激勵(lì)線圈,可獲得多個(gè)投影方向的檢測(cè)信息,再運(yùn)用合適的圖像重建算法就可以重建出放入物場(chǎng)中物體的電導(dǎo)率或者磁導(dǎo)率的分布圖,從而獲得放入物體的信息。由此可見(jiàn),EMT具有非侵入、非接觸和快速性的特點(diǎn),目前在生物醫(yī)學(xué)、無(wú)損檢測(cè)、鋼軌探傷、異物檢測(cè)和工業(yè)煉鋼等領(lǐng)域均已展開(kāi)相應(yīng)研究,具有很好的應(yīng)用以及發(fā)展前景[4-6]。

    在EMT系統(tǒng)中,圖像重建算法是EMT技術(shù)研究中的重點(diǎn)和難點(diǎn),圖像重建算法的好壞將直接影響最終成像的準(zhǔn)確性。目前的圖像重建算法主要有線性反投影(LBP)[7]、Tikhonov正則化[8]、Landweber迭代[9]等算法。在傳統(tǒng)的EMT系統(tǒng)中,邊界測(cè)量值的微小變化就會(huì)導(dǎo)致求解出的物場(chǎng)內(nèi)的介質(zhì)分布有很大變化,導(dǎo)致求解過(guò)程不穩(wěn)定,因此EMT技術(shù)本身具有病態(tài)性,嚴(yán)重影響了重建圖像的質(zhì)量。目前常采用截?cái)嗥娈愔邓惴ê蚑ikhonov正則化算法來(lái)改善靈敏度矩陣的病態(tài)性[10-11],但都存在參數(shù)選取問(wèn)題。Tikhonov正則化算法主要依賴于正則化因子,但是正則化因子的選取較為復(fù)雜,在實(shí)際應(yīng)用中往往依據(jù)經(jīng)驗(yàn)進(jìn)行選取,具有較大的局限性。如果正則化參數(shù)選擇較小,原始的問(wèn)題就不能得到很好地近似,逆問(wèn)題的病態(tài)性的改進(jìn)效果不顯著;如果正則化參數(shù)選擇過(guò)大,雖然對(duì)誤差的敏感性會(huì)減少,但最終解會(huì)與正常值相差較大,達(dá)不到效果。

    電磁層析成像的逆問(wèn)題的病態(tài)性的根源是靈敏度矩陣的條件數(shù)較大,本文通過(guò)利用人群搜索算法(SOA,seeker optimization algorithm)對(duì)物場(chǎng)為空時(shí)的靈敏度矩陣進(jìn)行優(yōu)化,能夠明顯降低其條件數(shù)。將優(yōu)化過(guò)后的靈敏度矩陣應(yīng)用于Landweber迭代算法進(jìn)行圖像重建,仿真實(shí)驗(yàn)結(jié)果表明,該方法能夠明顯提高圖像重建的效果以及精度。

    1 理論基礎(chǔ)

    1.1 電磁層析成像數(shù)學(xué)模型

    電磁層析成像的研究包括正問(wèn)題以及逆問(wèn)題。正問(wèn)題是具有電磁特性的物質(zhì)在敏感場(chǎng)中的分布位置已知,求解其對(duì)激勵(lì)場(chǎng)的調(diào)制的結(jié)果;逆問(wèn)題則相反,是利用在正問(wèn)題中已經(jīng)得出的調(diào)制結(jié)果,來(lái)求解具有電磁特性的物質(zhì)的在物場(chǎng)的分布。在EMT中,其調(diào)制結(jié)果往往表現(xiàn)為邊界測(cè)量值,即線圈的感應(yīng)電壓。二者的關(guān)系可由圖1所示。

    圖1 正問(wèn)題與逆問(wèn)題關(guān)系圖

    EMT正問(wèn)題中的模型中線圈感應(yīng)電壓可用式(1)描述:

    Udec=f(Be(x,y),Bo(x,y),μ(x,y),σ(x,y),ε(x,y))

    (1)

    式中,Udec為測(cè)量電壓,Be為激勵(lì)主磁場(chǎng),Bo為由渦流產(chǎn)生的次級(jí)磁場(chǎng),σ為電導(dǎo)率,μ為磁導(dǎo)率,ε為介電常數(shù)。

    EMT的逆問(wèn)題就是利用已知的靈敏度矩陣S以及邊界測(cè)量值U來(lái)解出灰度矩陣G,再根據(jù)所求出的灰度矩陣G即可對(duì)靈敏場(chǎng)中被測(cè)物體進(jìn)行圖像重建。由于檢測(cè)線圈獲得的邊界測(cè)量值與物場(chǎng)空間之間的關(guān)系是非線性的,因此EMT的逆問(wèn)題可近似歸結(jié)為線性方程組式(2)的求解[12]。

    U=S·G

    (2)

    式中,U∈RM為歸一化之后的邊界測(cè)量值矩陣,S∈RM×N為歸一化后的靈敏度矩陣,G∈RN為歸一化之后的灰度矩陣。M為獨(dú)立測(cè)量值的個(gè)數(shù),N是逆問(wèn)題剖分單元數(shù)。

    1.2 EMT靈敏度矩陣分析

    在EMT逆問(wèn)題的成像計(jì)算中,靈敏度矩陣作為先驗(yàn)信息,反應(yīng)了特定的激勵(lì)檢測(cè)對(duì)于不同位置單元介質(zhì)變化的敏感程度,是圖像重建的基礎(chǔ)。在實(shí)際應(yīng)用中,獲取靈敏度矩陣主要有3種方法:模型擾動(dòng)法[13]、測(cè)量擾動(dòng)法[14]和場(chǎng)量提取法[15]。因模型擾動(dòng)法及測(cè)量擾動(dòng)法計(jì)算量大且較為繁瑣,因此在本文中,將采用簡(jiǎn)便易行,計(jì)算速度快的場(chǎng)量提取法。

    采用場(chǎng)量提取法進(jìn)行靈敏度矩陣的求取并畫(huà)出典型激勵(lì)檢測(cè)組合的靈敏度分布如圖2所示。

    圖2 典型靈敏度矩陣分布圖

    1.3 人群搜索算法應(yīng)用分析

    電磁層析成像中的圖像重建屬于第一類Fredholm積分方程問(wèn)題,其解具有不適定性,這主要表現(xiàn)在圖像重建的過(guò)程中靈敏度矩陣的病態(tài)性。條件數(shù)可以用來(lái)衡量靈敏度矩陣的病態(tài)程度,因此,可以采取一定的算法令靈敏度矩陣的條件數(shù)達(dá)到最小,從而改善其病態(tài)性,即轉(zhuǎn)化為一個(gè)優(yōu)化問(wèn)題。

    人群搜索算法[16]能夠直接模擬人的智能搜索行為,立足傳統(tǒng)的直接搜索算法,以搜尋隊(duì)伍為種群,搜尋者位置為候選解,通過(guò)模擬人類搜尋“經(jīng)驗(yàn)梯度”和不確定推理進(jìn)行搜索方向和搜索步長(zhǎng)的更新,并不斷更新位置,完成對(duì)問(wèn)題的最優(yōu)求解。

    人群搜索算法主要包含搜索步長(zhǎng)、搜索方向以及個(gè)體位置的更新[17],為了簡(jiǎn)化模型,提高運(yùn)算效率,在本文中,對(duì)目標(biāo)函數(shù)值進(jìn)行排序時(shí)采用線性隸屬函數(shù),定義如式(3)和式(4)所示:

    (3)

    Uij=Ui+(1-Ui)·rand,j=1,2,,D

    (4)

    式中,Ui是目標(biāo)函數(shù)i的隸屬度;Uij是j維搜索空間目標(biāo)函數(shù)值i的隸屬度;P為種群大?。籇為搜索空間維數(shù);Ki為種群函數(shù)值按照降序排列之后的序列編號(hào);rand()能夠生成0~1之間的隨機(jī)數(shù),用來(lái)模仿人的隨機(jī)行為。

    1.4 Landweber迭代算法

    由于邊界測(cè)量值與敏感場(chǎng)中的電導(dǎo)率分布是非線性的關(guān)系,因此通過(guò)任何非迭代的算法都幾乎不可能找到精確解。目前Landweber迭代算法具有良好的成像質(zhì)量以及精度,是目前應(yīng)用最廣泛的圖像重建算法之一。

    從原理上來(lái)講,Landweber算法是最速下降法的一種變形,其本質(zhì)是通過(guò)迭代來(lái)計(jì)算靈敏度矩陣的廣義逆矩陣[18]。一般地,Landweber迭代法更多的用來(lái)解決Fredholm積分方程的典型病態(tài)問(wèn)題。

    本文將采用LBP重建的圖像作為迭代的初值,Landweber迭代算法如式(5)所示:

    (5)

    式中,β是增益因子,用來(lái)控制收斂速度,gk+1為第k步迭代的圖像灰度值,g0為利用線性反投影算法計(jì)算得到的灰度初始值,S為靈敏度矩陣,U為測(cè)量電壓值。

    2 基于人群搜索算法優(yōu)化的EMT圖像重建

    基于人群搜索算法優(yōu)化的電磁層析圖像重建的流程主要分為兩個(gè)部分:第一部分是利用人群搜索算法對(duì)利用場(chǎng)量提取法得到的靈敏度矩陣進(jìn)行預(yù)處理,降低其條件數(shù);第二部分則是根據(jù)優(yōu)化之后靈敏度矩陣?yán)脠D像重建算法進(jìn)行成像。針對(duì)第二部分,本文所采用圖像重建算法為L(zhǎng)andweber迭代算法,為了便于表示,采用此方式進(jìn)行圖像重建結(jié)果的算法命名為O_Landweber迭代算法,流程如下:

    1)參數(shù)的初始化,在人群搜索算法中,設(shè)定種群規(guī)模為sizepop,迭代次數(shù)為maxgen,空間維數(shù)為D,最大隸屬度值Umax,最小隸屬度值Umin,權(quán)重最大值Wmax,權(quán)重最小值Wmin,種群個(gè)體最大值為popmax,種群個(gè)體最小值為popmin,EMT逆問(wèn)題剖分網(wǎng)格數(shù)為N,Landweber迭代算法的迭代次數(shù)為V。

    2)按照1.2節(jié)的方法計(jì)算靈敏度矩陣S,維度為D*N。

    3)按照參數(shù)初始化中定義的種群中個(gè)體的最大值和最小值,完全隨機(jī)化生成初始種群T,維度為1*D。

    4)根據(jù)人群搜索算法以適應(yīng)度值最小為最優(yōu)進(jìn)行搜索的特點(diǎn)[19],對(duì)文獻(xiàn)[19]中提出的適應(yīng)度值的求取公式進(jìn)行取負(fù)操作,并且以矩陣T中的元素為對(duì)角線元素生成對(duì)角陣T1,然后計(jì)算種群中的個(gè)體的適應(yīng)度值。定義如式(6)和式(7):

    T1=diag(T)

    (6)

    (7)

    其中:P是懲罰函數(shù),定義如式(8)所示:

    (8)

    式中,cond(T1)表示T1的條件數(shù)。

    5)按照1.3節(jié)的方法對(duì)搜索步長(zhǎng)、搜索方向進(jìn)行更新。

    6)如果達(dá)到迭代次數(shù),進(jìn)行下一步,否則返回第(4)步。

    7)我們依據(jù)所改善的靈敏度矩陣的條件數(shù)占初始靈敏度條件數(shù)的比例來(lái)最終衡量靈敏度矩陣的優(yōu)化效果,定義如式(9)所示,delta越大表明優(yōu)化效果越好。在人群搜索算法中,本文采取雙停止策略,即若連續(xù)循環(huán)3次delta的變化量都不超過(guò)0.1%,或者已經(jīng)循環(huán)優(yōu)化了20次,則進(jìn)入(8),否則進(jìn)入(4)繼續(xù)更新T。

    (9)

    式中,T1best表示該次循環(huán)后種群中的最佳個(gè)體。

    8)在得到最終優(yōu)化的T1best之后,根據(jù)式(10)更新靈敏度矩陣以及測(cè)量電壓。本文將按照線性反投影算法得到的初值進(jìn)行Landweber迭代算法的計(jì)算,將更新后的靈敏度矩陣以及測(cè)量電壓應(yīng)用到式(5)所定義的Landweber迭代算法中進(jìn)行迭代,直到達(dá)到迭代次數(shù)V時(shí) 停止,進(jìn)行圖像重建。

    (10)

    式中,U1和S1分別表示優(yōu)化之后的測(cè)量電壓以及靈敏度矩陣。

    基于SOA優(yōu)化的EMT圖像重建的流程圖如圖3所示。

    圖3 基于SOA優(yōu)化的EMT圖像重建流程圖

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

    為了驗(yàn)證本文所提出的利用人群搜索算法改善靈敏度矩陣的病態(tài)程度,然后進(jìn)行圖像重建的有效性,主要選取了6種物場(chǎng)分布,采用傳統(tǒng)的8線圈EMT系統(tǒng)結(jié)構(gòu)進(jìn)行仿真[20-21]。為了比較不同算法的仿真效果,本文采用LBP,Tikhonov正則化,Landweber迭代算法以及第2節(jié)提出的O_Landweber迭代算法進(jìn)行成像。

    3.1 仿真實(shí)驗(yàn)環(huán)境設(shè)定

    與二維仿真相比,三維仿真更能夠模擬實(shí)際空間的復(fù)雜情況,準(zhǔn)確地表達(dá)出各物體之間相互關(guān)系,得到的仿真結(jié)果也就更為準(zhǔn)確,故本文采用電磁場(chǎng)有限元仿真軟件對(duì)仿真模型進(jìn)行三維建模。利用Matlab來(lái)控制在仿真過(guò)程中激勵(lì)線圈的切換,激勵(lì)電流的設(shè)置以及敏感場(chǎng)中物理量的獲取。

    在建模時(shí)設(shè)定敏感場(chǎng)半徑為50 mm,線圈半徑為16 mm,激勵(lì)頻率為200 kHz,激勵(lì)電流為10 A,線圈材質(zhì)為銅。在進(jìn)行實(shí)驗(yàn)時(shí),所放入的圓柱形物質(zhì)為半徑為8 mm的銅棒,其余為空氣,本文所建立的模型如圖4所示。

    圖4 三維仿真模型

    設(shè)定電磁層析成像的剖分單元格數(shù)N為322,如圖5所示;在人群搜索算法中,設(shè)定種群規(guī)模為sizepop為100,迭代次數(shù)為maxgen為50,空間維數(shù)D為28,最大隸屬度值Umax為0.950 0,最小隸屬度值Umin為0.011 1,權(quán)重最大值Wmax為0.9,權(quán)重最小值Wmin為0.1,種群個(gè)體最大值為popmax為1,種群個(gè)體最小值popmin為0;設(shè)定Landweber迭代算法的迭代次數(shù)V為1 000次[22]。

    圖5 逆問(wèn)題網(wǎng)格剖分結(jié)果

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

    本文的仿真是在Matlab以及電磁場(chǎng)有限元仿真軟件中進(jìn)行的,利用Matlab來(lái)控制有限元仿真軟件的仿真進(jìn)程、仿真數(shù)據(jù)獲取以及圖像重建。仿真過(guò)程主要分為3個(gè)部分:第一部分為靈敏度矩陣的計(jì)算,第二部分為邊界測(cè)量值獲取,第三部分為圖像重建。首先獲取靈敏度矩陣,計(jì)算每個(gè)線圈在單位電流激勵(lì)下在被測(cè)物場(chǎng)中產(chǎn)生的磁場(chǎng),然后確定激勵(lì)和檢測(cè)線圈的位置,從場(chǎng)量中提取靈敏度[15],從而得到整個(gè)物場(chǎng)的靈敏度矩陣;獲取邊界測(cè)量值則采用單線圈激勵(lì)的策略,在電磁場(chǎng)有限元仿真軟件中,對(duì)3.1節(jié)中建立的模型中的8個(gè)銅線圈依次激勵(lì),共可得到8×8個(gè)測(cè)量值;最后依據(jù)靈敏度矩陣和邊界測(cè)量值采用不同的圖像重建算法在MATLAB中進(jìn)行圖像重建。

    本文通過(guò)采用循環(huán)利用人群搜索算法來(lái)改善靈敏度矩陣的病態(tài)程度,適應(yīng)度值和靈敏度矩陣條件數(shù)的收斂曲線如圖5所示,本次實(shí)驗(yàn)共循環(huán)了11次,每次循環(huán)次數(shù)maxgen為50,為了將所有的循環(huán)過(guò)程都呈現(xiàn)在一副圖片上,設(shè)定橫坐標(biāo)為總迭代次數(shù),共迭代11*50=550代,縱坐標(biāo)為每次迭代所對(duì)應(yīng)的適應(yīng)度值以及靈敏度矩陣條件數(shù)。

    圖6 適應(yīng)度值和靈敏度矩陣條件數(shù)收斂曲線

    由圖6可知,根據(jù)適應(yīng)度值收斂曲線可以看出,迭代次數(shù)小于300時(shí),隨著迭代次數(shù)的增加適應(yīng)度值下降明顯,當(dāng)?shù)螖?shù)大于300時(shí),適應(yīng)度值下降趨于平緩,逐漸收斂,其表現(xiàn)為靈敏度矩陣條件數(shù)也趨于收斂。根據(jù)靈敏度矩陣條件數(shù)收斂曲線可以看出,用人群搜索算法對(duì)靈敏度矩陣進(jìn)行優(yōu)化,能夠有效降低其條件數(shù)。在迭代次數(shù)小于100時(shí)靈敏度矩陣條件數(shù)下降明顯,當(dāng)?shù)螖?shù)大于300時(shí),適應(yīng)度值趨于穩(wěn)定,與適應(yīng)度值收斂曲線表現(xiàn)基本一致。

    靈敏度矩陣條件數(shù)的改善量收斂曲線如圖7所示。

    圖7 靈敏度矩陣條件數(shù)改善量收斂曲線

    由圖7可知,在人群搜索算法循環(huán)迭代的過(guò)程中,靈敏度矩陣條件數(shù)改善量隨循環(huán)次數(shù)的增加不斷增加;在增加的過(guò)程中,曲線的斜率逐漸減小,表明靈敏度矩陣條件數(shù)的改善速率不斷減小,數(shù)值差異也不斷縮小,即可以停止迭代優(yōu)化。結(jié)合圖6來(lái)看靈敏度矩陣的最初的條件數(shù)為632.621 8,優(yōu)化之后降至77.009 9,共降低了87.83%,有效地改善了靈敏度矩陣的病態(tài)程度仿真得到的圖像重建結(jié)果如表1所示。

    表1 仿真圖像重建結(jié)果

    如表1所示,第一列為仿真的物體的真實(shí)分布圖像,第二列為采用LBP算法重建所得到的圖像,第三列為采用Tikhonov正則化算法重建所得到的圖像,第四列為采用Landweber迭代算法重建所得到的圖像,第五列為本文提出O_Landweber迭代算法重建所得到的圖像。相對(duì)于LBP算法以及Tikhonov正則化算法,通過(guò)Landweber迭代算法以及本文所提出O_Landweber算法重建所得到的圖像能夠較為準(zhǔn)確的重建出物體所在的位置,并且能夠減少偽影。由表中第四列Landweber算法重建的圖像可知,在物體周圍存在較多的偽影,物體的位置還不夠精確,本文提出的O_Landweber迭代算法能夠進(jìn)一步減少物體周圍的偽影,還原物體的實(shí)際位置。

    為了以定量的方式來(lái)對(duì)以上4種算法的成像質(zhì)量以及成像精度進(jìn)行比較,本文采用圖像相對(duì)誤差(IE,image error)以及相關(guān)系數(shù)(CC,correlation coefficient)作為評(píng)價(jià)指標(biāo)進(jìn)行比較[23]。圖像相關(guān)系數(shù)是衡量變量之間線性相關(guān)程度的指標(biāo);圖像誤差是重建圖像與原始圖像之間的差異程度。若將采用兩種不同成像算法同一種物場(chǎng)分布的仿真結(jié)果進(jìn)行比較,如果其中一種仿真圖像的圖像相對(duì)誤差更小,而相關(guān)系數(shù)更大,那么就說(shuō)明該仿真結(jié)果的圖像重建的質(zhì)量最好。定義如式(11)和(12)所示:

    (11)

    (12)

    表2和表3給出了6中不同分布下4種不同的算法所得到的圖像重建結(jié)果的圖像誤差以及圖像相關(guān)系數(shù)。

    表2 不同圖像重建算法對(duì)各類型分布的圖像相關(guān)系數(shù)

    表3 不同圖像重建算法對(duì)各典型分布的圖像誤差模型

    由表2、表3可知,當(dāng)被測(cè)物場(chǎng)中含有一個(gè)物體時(shí),本文提出的O_Landweber迭代算法相比Landweber迭代算法的圖像相關(guān)系數(shù)能夠提升22%以上,圖像誤差能夠降低56%以上;當(dāng)被測(cè)物場(chǎng)中含有兩個(gè)物體時(shí),不論是圖像相關(guān)系數(shù)還是圖像誤差,改善幅度均沒(méi)有一個(gè)物體時(shí)明顯,但仍有提升,圖像相關(guān)系數(shù)能夠提升8%以上,圖像誤差能夠降低28%以上。本文提出的O_Landweber迭代算法相對(duì)于LBP、Tikhonov正則化以及Landweber迭代算法的圖像相關(guān)系數(shù)更大,圖像誤差更小,這表明其所成圖像與原始分布的圖像誤差較小且線性相關(guān)程度較高,與仿真圖像重建結(jié)果表現(xiàn)的基本一致。

    4 結(jié)束語(yǔ)

    由于EMT自身特性,導(dǎo)致靈敏度矩陣往往都具有病態(tài)性。為了進(jìn)一步提高電磁層析成像的成像質(zhì)量,本文提出先通過(guò)人群算法降低靈敏度矩陣病態(tài)程度,然后再進(jìn)行圖像重建來(lái)提高成像質(zhì)量。通過(guò)循環(huán)利用人群搜索算法對(duì)靈敏度矩陣進(jìn)行不斷優(yōu)化的方法,使靈敏度矩陣的條件數(shù)能夠在迭代300代之后趨于收斂,相對(duì)于初始靈敏度矩陣條件數(shù)能夠降低87.83%,靈敏度矩陣病態(tài)程度得到明顯改善。仿真實(shí)驗(yàn)表明:與傳統(tǒng)的直接用于成像的電磁層析成像算法相比,本文所提出的O_Landweber迭代算法能夠?qū)D像相關(guān)系數(shù)最高提升45%,降低圖像誤差最高75%;與LBP算法、Tikhonov正則化算法,Landweber迭代算法相比,本文所提出的算法能夠明顯消除偽影,提高成像質(zhì)量,進(jìn)一步精確物體的位置。

    猜你喜歡
    病態(tài)搜索算法正則
    改進(jìn)的和聲搜索算法求解凸二次規(guī)劃及線性規(guī)劃
    病態(tài)肥胖對(duì)門診全關(guān)節(jié)置換術(shù)一夜留院和早期并發(fā)癥的影響
    病態(tài)肥胖對(duì)門診關(guān)節(jié)置換術(shù)留夜觀察和早期并發(fā)癥的影響
    剩余有限Minimax可解群的4階正則自同構(gòu)
    類似于VNL環(huán)的環(huán)
    君子之道:能移而相天——王夫之《莊子解》對(duì)“社會(huì)病態(tài)”的氣論診療
    基于汽車接力的潮流轉(zhuǎn)移快速搜索算法
    基于逐維改進(jìn)的自適應(yīng)步長(zhǎng)布谷鳥(niǎo)搜索算法
    有限秩的可解群的正則自同構(gòu)
    基于跳點(diǎn)搜索算法的網(wǎng)格地圖尋路
    嫩草影院入口| 国产高清videossex| 黑人欧美特级aaaaaa片| 欧美不卡视频在线免费观看| 亚洲自拍偷在线| 亚洲av成人av| 男人的好看免费观看在线视频| 亚洲精品456在线播放app | 老司机午夜十八禁免费视频| 波多野结衣巨乳人妻| 欧美日韩黄片免| 日本 av在线| 18禁美女被吸乳视频| 亚洲成人久久爱视频| 亚洲av成人精品一区久久| 国产精品1区2区在线观看.| 一级a爱片免费观看的视频| 九九热线精品视视频播放| 日韩亚洲欧美综合| 国产精品野战在线观看| 国产爱豆传媒在线观看| 亚洲内射少妇av| 给我免费播放毛片高清在线观看| 国产精品一区二区免费欧美| 身体一侧抽搐| 亚洲熟妇中文字幕五十中出| 黄片大片在线免费观看| 久久人妻av系列| 国产精品1区2区在线观看.| 国产一区在线观看成人免费| 岛国在线免费视频观看| 亚洲第一欧美日韩一区二区三区| 日本熟妇午夜| 在线观看一区二区三区| 99久久无色码亚洲精品果冻| www.熟女人妻精品国产| 在线天堂最新版资源| 男女做爰动态图高潮gif福利片| 日本a在线网址| 成年免费大片在线观看| 亚洲国产精品合色在线| 国产免费av片在线观看野外av| 久久久久久久精品吃奶| 99热这里只有精品一区| 麻豆成人av在线观看| 久久精品夜夜夜夜夜久久蜜豆| 亚洲av第一区精品v没综合| 成年免费大片在线观看| 天天一区二区日本电影三级| 亚洲精品日韩av片在线观看 | 久久久久久久亚洲中文字幕 | 深夜精品福利| 亚洲成av人片免费观看| 免费在线观看日本一区| av天堂在线播放| 九九热线精品视视频播放| 三级毛片av免费| 国产伦精品一区二区三区四那| 老汉色∧v一级毛片| 国产欧美日韩精品亚洲av| 亚洲精品久久国产高清桃花| 亚洲人成网站在线播| 精品久久久久久久久久免费视频| 999久久久精品免费观看国产| 精品电影一区二区在线| 淫秽高清视频在线观看| 亚洲第一欧美日韩一区二区三区| 一个人看的www免费观看视频| 一a级毛片在线观看| 69人妻影院| 91麻豆av在线| 国产老妇女一区| 一进一出好大好爽视频| 亚洲av免费高清在线观看| 色综合亚洲欧美另类图片| 18禁美女被吸乳视频| 99久久综合精品五月天人人| 波野结衣二区三区在线 | 成年女人毛片免费观看观看9| 亚洲一区二区三区不卡视频| 亚洲精品成人久久久久久| 国产一区二区激情短视频| 成人性生交大片免费视频hd| 国产午夜精品论理片| 国内精品美女久久久久久| 国产精品一区二区三区四区久久| 97超级碰碰碰精品色视频在线观看| 小蜜桃在线观看免费完整版高清| av天堂在线播放| 欧美性猛交╳xxx乱大交人| 欧美3d第一页| 老熟妇乱子伦视频在线观看| 在线免费观看的www视频| 美女cb高潮喷水在线观看| 国产综合懂色| 国产乱人视频| 俄罗斯特黄特色一大片| 亚洲久久久久久中文字幕| 亚洲色图av天堂| 久久精品国产亚洲av香蕉五月| 国产精品久久视频播放| 天堂网av新在线| 两个人的视频大全免费| 深爱激情五月婷婷| 蜜桃久久精品国产亚洲av| 婷婷精品国产亚洲av| 国产精品久久久人人做人人爽| 午夜免费激情av| 精品国内亚洲2022精品成人| 人人妻,人人澡人人爽秒播| 国产精品精品国产色婷婷| 亚洲欧美精品综合久久99| 少妇高潮的动态图| 国产成年人精品一区二区| 99在线视频只有这里精品首页| 看片在线看免费视频| 波多野结衣高清无吗| 99国产综合亚洲精品| 狂野欧美激情性xxxx| 婷婷亚洲欧美| 一二三四社区在线视频社区8| 久99久视频精品免费| 亚洲国产欧美网| 午夜影院日韩av| 两个人视频免费观看高清| 内地一区二区视频在线| 亚洲国产高清在线一区二区三| 国产精品一区二区三区四区久久| 老司机午夜十八禁免费视频| 久久国产精品影院| 哪里可以看免费的av片| 日本黄色片子视频| 在线观看66精品国产| 久久久久久久久久黄片| 午夜激情欧美在线| 白带黄色成豆腐渣| 哪里可以看免费的av片| 精品午夜福利视频在线观看一区| 精品一区二区三区视频在线观看免费| 国产精品自产拍在线观看55亚洲| 高清在线国产一区| 男人舔奶头视频| 日韩欧美免费精品| 真实男女啪啪啪动态图| 国产高清视频在线播放一区| 欧美日韩乱码在线| 男女那种视频在线观看| 欧美+亚洲+日韩+国产| 国产精品女同一区二区软件 | 亚洲欧美激情综合另类| 亚洲专区国产一区二区| 国产一区二区亚洲精品在线观看| 99热这里只有精品一区| 久久欧美精品欧美久久欧美| 最近最新免费中文字幕在线| 成人精品一区二区免费| 一a级毛片在线观看| 桃红色精品国产亚洲av| 国产午夜精品论理片| 婷婷六月久久综合丁香| 三级毛片av免费| 亚洲国产精品合色在线| 琪琪午夜伦伦电影理论片6080| 欧美丝袜亚洲另类 | 成人特级黄色片久久久久久久| 波多野结衣高清作品| 美女 人体艺术 gogo| 天天一区二区日本电影三级| 免费av毛片视频| 国产淫片久久久久久久久 | 久久国产精品影院| 叶爱在线成人免费视频播放| 男女那种视频在线观看| www日本黄色视频网| 禁无遮挡网站| 制服人妻中文乱码| 亚洲电影在线观看av| 琪琪午夜伦伦电影理论片6080| 午夜激情福利司机影院| 欧美最新免费一区二区三区 | 国产精品三级大全| 一卡2卡三卡四卡精品乱码亚洲| 亚洲国产精品成人综合色| 男女那种视频在线观看| 精品久久久久久,| 俺也久久电影网| 真人做人爱边吃奶动态| 3wmmmm亚洲av在线观看| 真人一进一出gif抽搐免费| 久久精品亚洲精品国产色婷小说| 久久精品影院6| 亚洲精品456在线播放app | 欧美中文综合在线视频| 999久久久精品免费观看国产| 深爱激情五月婷婷| 人妻丰满熟妇av一区二区三区| 国产午夜精品久久久久久一区二区三区 | 久9热在线精品视频| 好看av亚洲va欧美ⅴa在| 九九久久精品国产亚洲av麻豆| 精品一区二区三区视频在线 | 高清毛片免费观看视频网站| 亚洲一区二区三区不卡视频| 国产淫片久久久久久久久 | 成人高潮视频无遮挡免费网站| 精品乱码久久久久久99久播| 久久草成人影院| 亚洲中文字幕日韩| 久久久久久大精品| 欧美国产日韩亚洲一区| 免费看a级黄色片| 欧美在线一区亚洲| 成人精品一区二区免费| 国产成年人精品一区二区| 女人被狂操c到高潮| 午夜福利在线观看免费完整高清在 | 中文字幕人妻丝袜一区二区| 中出人妻视频一区二区| 999久久久精品免费观看国产| 欧美日韩乱码在线| av专区在线播放| 搡女人真爽免费视频火全软件 | 岛国在线观看网站| 久久久国产精品麻豆| 在线播放国产精品三级| 久久国产乱子伦精品免费另类| 久久精品91无色码中文字幕| 欧美乱色亚洲激情| 久久久久久国产a免费观看| 久久伊人香网站| 91字幕亚洲| 免费看十八禁软件| 免费看光身美女| 女生性感内裤真人,穿戴方法视频| 亚洲av不卡在线观看| 国产免费一级a男人的天堂| 男女之事视频高清在线观看| 国产亚洲精品综合一区在线观看| 欧美黄色片欧美黄色片| 国产精品亚洲av一区麻豆| 色哟哟哟哟哟哟| 91在线精品国自产拍蜜月 | 亚洲欧美一区二区三区黑人| 久久久精品大字幕| 婷婷六月久久综合丁香| 亚洲精品456在线播放app | 在线天堂最新版资源| 日韩欧美免费精品| 亚洲国产欧美网| 黄片小视频在线播放| 国产精品亚洲av一区麻豆| 亚洲电影在线观看av| 深夜精品福利| 国产精品久久视频播放| 国产精品香港三级国产av潘金莲| 听说在线观看完整版免费高清| 国产精品久久视频播放| 欧美成人免费av一区二区三区| 国产伦精品一区二区三区四那| 精品一区二区三区视频在线观看免费| 亚洲,欧美精品.| 两人在一起打扑克的视频| 动漫黄色视频在线观看| 精品久久久久久久末码| 欧美+日韩+精品| 国产一区二区在线av高清观看| 国产成人啪精品午夜网站| 此物有八面人人有两片| 日韩亚洲欧美综合| 在线天堂最新版资源| 狠狠狠狠99中文字幕| 日日夜夜操网爽| 99久久精品热视频| 精品午夜福利视频在线观看一区| 2021天堂中文幕一二区在线观| 黄片大片在线免费观看| 亚洲av五月六月丁香网| 久久6这里有精品| 波多野结衣巨乳人妻| 88av欧美| 精品日产1卡2卡| 亚洲精品在线观看二区| 欧美色视频一区免费| 日本成人三级电影网站| 亚洲黑人精品在线| 亚洲国产日韩欧美精品在线观看 | 黄色丝袜av网址大全| 久久精品人妻少妇| 可以在线观看的亚洲视频| 国产精品久久久久久久电影 | 亚洲中文字幕一区二区三区有码在线看| 69av精品久久久久久| 好看av亚洲va欧美ⅴa在| 麻豆国产97在线/欧美| 亚洲天堂国产精品一区在线| 亚洲片人在线观看| 日韩人妻高清精品专区| 色吧在线观看| 黄色片一级片一级黄色片| 婷婷精品国产亚洲av| 天美传媒精品一区二区| 内射极品少妇av片p| or卡值多少钱| 怎么达到女性高潮| 国产精品一区二区三区四区久久| 国产成人a区在线观看| 国产黄片美女视频| 听说在线观看完整版免费高清| 欧美最黄视频在线播放免费| 少妇的逼好多水| 级片在线观看| 国产一区二区在线观看日韩 | 少妇熟女aⅴ在线视频| 淫秽高清视频在线观看| 少妇的逼水好多| 天堂√8在线中文| 久久久久精品国产欧美久久久| 亚洲精品日韩av片在线观看 | 一区二区三区高清视频在线| 午夜精品久久久久久毛片777| 欧美区成人在线视频| 国产精品久久久久久人妻精品电影| 91九色精品人成在线观看| 丰满乱子伦码专区| 可以在线观看毛片的网站| 国产精品一及| 日本在线视频免费播放| 色综合亚洲欧美另类图片| 蜜桃久久精品国产亚洲av| 一级a爱片免费观看的视频| 美女高潮的动态| 欧美又色又爽又黄视频| 午夜激情福利司机影院| av天堂在线播放| 美女 人体艺术 gogo| 久久久久久九九精品二区国产| 亚洲五月婷婷丁香| 日韩中文字幕欧美一区二区| 午夜免费激情av| 九九在线视频观看精品| 亚洲五月婷婷丁香| 身体一侧抽搐| 亚洲最大成人手机在线| 悠悠久久av| 国产亚洲精品久久久com| 欧美最新免费一区二区三区 | 国产精品嫩草影院av在线观看 | 亚洲在线自拍视频| 88av欧美| 国产在线精品亚洲第一网站| 美女高潮喷水抽搐中文字幕| 国产探花在线观看一区二区| 国产一区二区亚洲精品在线观看| 小蜜桃在线观看免费完整版高清| 舔av片在线| 麻豆一二三区av精品| 久久精品影院6| 怎么达到女性高潮| 久久精品综合一区二区三区| 9191精品国产免费久久| 日韩欧美免费精品| 久久久久国产精品人妻aⅴ院| 在线观看美女被高潮喷水网站 | 午夜两性在线视频| 男插女下体视频免费在线播放| 99视频精品全部免费 在线| 亚洲精华国产精华精| 白带黄色成豆腐渣| www.www免费av| 激情在线观看视频在线高清| 久久亚洲真实| 国产精品1区2区在线观看.| 最近最新免费中文字幕在线| 日韩中文字幕欧美一区二区| 欧美色视频一区免费| 两个人的视频大全免费| 少妇人妻精品综合一区二区 | 国产蜜桃级精品一区二区三区| 日本在线视频免费播放| 国产精品免费一区二区三区在线| 99热这里只有是精品50| 女生性感内裤真人,穿戴方法视频| 国产av麻豆久久久久久久| 哪里可以看免费的av片| 亚洲国产欧美网| 国产精品自产拍在线观看55亚洲| 一个人看视频在线观看www免费 | 人妻夜夜爽99麻豆av| 国产真实乱freesex| 两个人的视频大全免费| 精品无人区乱码1区二区| 伊人久久大香线蕉亚洲五| 一级黄片播放器| 欧美日本亚洲视频在线播放| 亚洲精品456在线播放app | 午夜影院日韩av| 国产成人欧美在线观看| 国产不卡一卡二| svipshipincom国产片| 国产精品一及| 久久久精品大字幕| 欧美性猛交黑人性爽| 亚洲成av人片在线播放无| 久久久精品欧美日韩精品| 成人永久免费在线观看视频| а√天堂www在线а√下载| 国产精品乱码一区二三区的特点| 91久久精品国产一区二区成人 | 日本a在线网址| 婷婷亚洲欧美| 亚洲精品久久国产高清桃花| 人妻夜夜爽99麻豆av| 一进一出抽搐动态| 国产av在哪里看| 日本黄色片子视频| 亚洲精品亚洲一区二区| 午夜免费激情av| 日韩精品中文字幕看吧| 90打野战视频偷拍视频| 中文字幕人妻熟人妻熟丝袜美 | 熟女人妻精品中文字幕| 夜夜看夜夜爽夜夜摸| 久久久久久久精品吃奶| 日韩精品中文字幕看吧| 深爱激情五月婷婷| 亚洲无线观看免费| 亚洲熟妇中文字幕五十中出| 国产在线精品亚洲第一网站| www日本黄色视频网| 美女 人体艺术 gogo| 不卡一级毛片| 亚洲黑人精品在线| 成人永久免费在线观看视频| 国产一区二区亚洲精品在线观看| 熟女少妇亚洲综合色aaa.| 国产成人系列免费观看| 欧美成人免费av一区二区三区| av福利片在线观看| 久久久国产精品麻豆| 最新美女视频免费是黄的| 男人舔女人下体高潮全视频| 99热6这里只有精品| 国产精品久久久久久亚洲av鲁大| 热99re8久久精品国产| 亚洲成人久久性| 最近在线观看免费完整版| 亚洲电影在线观看av| 日本a在线网址| 搞女人的毛片| 黄色视频,在线免费观看| 真实男女啪啪啪动态图| 国产真人三级小视频在线观看| 午夜精品久久久久久毛片777| 老司机在亚洲福利影院| 久久久精品欧美日韩精品| 一级毛片女人18水好多| 久久亚洲精品不卡| 亚洲成a人片在线一区二区| 午夜福利高清视频| 久久久久九九精品影院| h日本视频在线播放| 一二三四社区在线视频社区8| 精品久久久久久久久久久久久| 欧美av亚洲av综合av国产av| 国产成人av教育| 老汉色∧v一级毛片| 一边摸一边抽搐一进一小说| 国产探花极品一区二区| а√天堂www在线а√下载| 国产探花极品一区二区| 乱人视频在线观看| 国产亚洲精品av在线| 欧美zozozo另类| 亚洲中文字幕日韩| 亚洲人成伊人成综合网2020| 色噜噜av男人的天堂激情| 高清在线国产一区| 欧美日韩中文字幕国产精品一区二区三区| 欧美日韩国产亚洲二区| 国产精品 国内视频| 午夜两性在线视频| 国产真实乱freesex| 日韩成人在线观看一区二区三区| 国产aⅴ精品一区二区三区波| 久久精品综合一区二区三区| 日韩精品中文字幕看吧| 精品一区二区三区视频在线 | 成人亚洲精品av一区二区| 欧美黑人欧美精品刺激| 在线播放无遮挡| 亚洲aⅴ乱码一区二区在线播放| 全区人妻精品视频| 三级毛片av免费| 国产不卡一卡二| 高潮久久久久久久久久久不卡| 极品教师在线免费播放| 黄色日韩在线| 精品熟女少妇八av免费久了| 国产一区在线观看成人免费| 日韩大尺度精品在线看网址| 一夜夜www| 精品国产超薄肉色丝袜足j| xxxwww97欧美| 国产私拍福利视频在线观看| 亚洲国产日韩欧美精品在线观看 | ponron亚洲| 精品午夜福利视频在线观看一区| 国产高清视频在线观看网站| 在线观看舔阴道视频| 欧美成人一区二区免费高清观看| 午夜福利高清视频| 99热这里只有精品一区| 国产一区二区激情短视频| 久久久久国产精品人妻aⅴ院| 男人和女人高潮做爰伦理| 亚洲自拍偷在线| 久久国产精品人妻蜜桃| 久久久精品欧美日韩精品| 亚洲专区国产一区二区| 一级黄片播放器| 久久99热这里只有精品18| 麻豆一二三区av精品| 国产精品久久久久久人妻精品电影| 两个人的视频大全免费| 成人亚洲精品av一区二区| 99在线人妻在线中文字幕| 熟妇人妻久久中文字幕3abv| 久久久久九九精品影院| 午夜福利18| 国产在线精品亚洲第一网站| 观看美女的网站| 精品一区二区三区视频在线观看免费| 性欧美人与动物交配| 美女高潮喷水抽搐中文字幕| 国产免费男女视频| 精品欧美国产一区二区三| 一级a爱片免费观看的视频| 日韩大尺度精品在线看网址| 成人午夜高清在线视频| 精品一区二区三区人妻视频| 人妻久久中文字幕网| www.999成人在线观看| 香蕉久久夜色| 嫁个100分男人电影在线观看| 国产野战对白在线观看| 淫妇啪啪啪对白视频| 日本a在线网址| 99久久久亚洲精品蜜臀av| 欧美另类亚洲清纯唯美| 亚洲av日韩精品久久久久久密| 成年版毛片免费区| 99热6这里只有精品| 91麻豆精品激情在线观看国产| 香蕉久久夜色| 亚洲欧美日韩高清专用| 在线观看一区二区三区| 波多野结衣巨乳人妻| 国模一区二区三区四区视频| 久久久久亚洲av毛片大全| 国产av一区在线观看免费| 一个人免费在线观看的高清视频| 熟女电影av网| 国产69精品久久久久777片| 男女床上黄色一级片免费看| 国内精品一区二区在线观看| 18禁黄网站禁片午夜丰满| www.色视频.com| 免费av观看视频| 久久久久国内视频| 亚洲精品一区av在线观看| 国产在线精品亚洲第一网站| 最近最新中文字幕大全免费视频| 日韩人妻高清精品专区| 精品久久久久久久久久久久久| 欧美性感艳星| 精品一区二区三区视频在线观看免费| 欧美成人免费av一区二区三区| 国产一区二区亚洲精品在线观看| 麻豆成人av在线观看| 女人被狂操c到高潮| 床上黄色一级片| 在线观看66精品国产| 欧美3d第一页| 19禁男女啪啪无遮挡网站| av天堂中文字幕网| 欧美色欧美亚洲另类二区| 亚洲18禁久久av| 免费无遮挡裸体视频| 日本撒尿小便嘘嘘汇集6| 国产精品免费一区二区三区在线| 亚洲欧美日韩东京热| 成人性生交大片免费视频hd| 国产精品久久久人人做人人爽| 久久久久免费精品人妻一区二区| 精品欧美国产一区二区三| 国产精品av视频在线免费观看| 真人做人爱边吃奶动态| 一本久久中文字幕| 亚洲精品日韩av片在线观看 | 久久99热这里只有精品18| 国产伦在线观看视频一区| 最近在线观看免费完整版| 久久99热这里只有精品18| 少妇的逼水好多| 精品免费久久久久久久清纯| 亚洲av二区三区四区| 国产真实乱freesex| 女人高潮潮喷娇喘18禁视频| 亚洲精品456在线播放app | 欧美3d第一页| 精品国产三级普通话版| 99热这里只有是精品50| 淫秽高清视频在线观看| 岛国在线观看网站| 青草久久国产|