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

    海鷗算法的改進(jìn)及其在工程設(shè)計(jì)優(yōu)化問(wèn)題 中的應(yīng)用

    2022-12-22 12:28:48楊韜戴健
    關(guān)鍵詞:測(cè)試函數(shù)海鷗學(xué)習(xí)策略

    楊韜, 戴健

    遼寧工程技術(shù)大學(xué),工商管理學(xué)院,遼寧 葫蘆島 125100

    引 言

    智能優(yōu)化算法在解決復(fù)雜優(yōu)化的問(wèn)題上有著不俗的表現(xiàn),如灰狼優(yōu)化算法[1]、麻雀優(yōu)化算法[2]、蝗蟲(chóng)優(yōu)化算法[3]、粒子群優(yōu)化算法[4]、鯨魚(yú)優(yōu)化算 法[5]、蟻獅優(yōu)化算法[6]、海鷗優(yōu)化算法[7]等。海鷗優(yōu)化算法是Dhiman和Kumar在2019年提出的新型群智能優(yōu)化算法,與傳統(tǒng)PSO和GA算法相比,其結(jié)構(gòu)簡(jiǎn)單、易于實(shí)現(xiàn),參數(shù)的調(diào)整也很簡(jiǎn)便,這些優(yōu)點(diǎn)使其受到廣泛應(yīng)用。然而海鷗優(yōu)化算法仍然存在收斂速度慢、易陷入局部最優(yōu)的缺陷,導(dǎo)致其多樣性較差。

    目前相關(guān)專(zhuān)家關(guān)于群智能優(yōu)化算法的研究主要分為兩個(gè)方面:一是種群初始化方面,毛清華[8]、陳忠云[9]等用Logistics混沌映射對(duì)種群初始化,尹德鑫[10]用Fuch混沌映射對(duì)種群初始化,肖亞寧[11]用Tent混沌映射對(duì)種群初始化,馬馳等[12]用Tent混沌映射與對(duì)立學(xué)習(xí)策略融合的方式對(duì)種群初始化。二是種群搜索尋優(yōu)方面,Cao Y[13]、秦維娜[14]等通過(guò)加入Levy飛行機(jī)制增加搜索的隨機(jī)性,Chen X通過(guò)引入非線性搜索控制公式[15]提高算法的搜索速度和精度。

    受上述文獻(xiàn)的啟發(fā),在加入萊維飛行機(jī)制的基礎(chǔ)上,首先融合Fuch混沌映射[16]與精英反向?qū)W習(xí)策略[11,17]對(duì)海鷗種群初始化,提高種群質(zhì)量,其次根據(jù)余弦函數(shù)改進(jìn)自身行為特征參數(shù)A,將線性搜索非線性化,并將改進(jìn)的海鷗算法與原先算法和其他智能優(yōu)化算法在9個(gè)測(cè)試函數(shù)和3個(gè)工程設(shè)計(jì)優(yōu)化問(wèn)題上進(jìn)行性能對(duì)比。

    1 海鷗算法

    海鷗優(yōu)化算法[7,18-19]啟發(fā)于海鷗的遷徙和攻擊行為。為了避免碰撞,各海鷗初始位置不同,在這一機(jī)制中,海鷗試圖朝著最佳生存方向移動(dòng),以確定最優(yōu)方案。

    1.1 海鷗的遷徙模式

    海鷗群體需要在避免和其他海鷗碰撞、向最佳的鄰居飛行和移動(dòng)到最佳位置的前提下進(jìn)行遷徙:

    (1)避免碰撞

    由式(1)所示,變量A避免了海鷗與其它相鄰海鷗的碰撞,并確定自身新的位置:

    式(1)中,Cs為避免碰撞的新位置,t為迭代次數(shù),

    PS(t)為海鷗的當(dāng)前位置。A是海鷗在搜索空間中進(jìn)行全局搜索時(shí)代表了其自身行為的特征參數(shù),具體表述如下:

    式(2)中,Maxiteration為最大迭代次數(shù),fc為控制變量的頻率,fc值一般為2。

    (2)其他鄰居經(jīng)驗(yàn)

    海鷗中的候選個(gè)體會(huì)朝最佳鄰居的方向移動(dòng)探索優(yōu)值。

    式(3)中,Bs表示海鷗中的最優(yōu)候選位置,參數(shù)B為隨機(jī)參數(shù),用來(lái)權(quán)衡算法的探索和利用,Pbs(t)為海鷗個(gè)體向最優(yōu)位置移動(dòng)。

    式(4)中,R為[0,1]之間的隨機(jī)數(shù)。

    (3)移動(dòng)到最佳位置搜索代理

    更新步驟具體表述如下:

    式(5)中,DS表示海鷗個(gè)體搜索位置CS和全局最佳搜索位置BS之間的距離。

    1.2 海鷗的攻擊模式

    海鷗捕獵時(shí),呈螺旋狀移動(dòng),通過(guò)不斷改變攻擊的角度和速度改變螺旋的半徑,其飛行的高度可通過(guò)翅膀和體重維持。具體運(yùn)動(dòng)行為表述如下:

    式(6)中,r為海鷗在進(jìn)行螺旋運(yùn)動(dòng)時(shí)的半徑,a是定義的一個(gè)隨機(jī)角度,具體范圍是[0,2π],u和v是定值,一般情況下均取1。

    式(7)中,PS(t)為海鷗的攻擊位置,并更新其它海鷗搜索位置。

    2 改進(jìn)的海鷗算法

    2.1 加入混沌精英反向?qū)W習(xí)策略

    (1)Fuch映射

    海鷗優(yōu)化算法的初始種群存在分布不均勻的問(wèn)題,導(dǎo)致其相關(guān)性能不穩(wěn)定?;煦缬成淇衫闷浞蔷€性、普適性、隨機(jī)性、遍歷性和規(guī)律性的優(yōu)點(diǎn)增強(qiáng)種群多樣性和全局搜索能力。Fuch映射[10,16]在混沌特性、遍歷性和規(guī)律性等方面比傳統(tǒng)的Logistic映射和Tent映射更加優(yōu)越,具體表述如下(xn不為0):

    式(8)中xn不為0,n為當(dāng)前迭代次數(shù)。

    (2)精英反向?qū)W習(xí)策略

    精英反向?qū)W習(xí)策略[11,17]為反向?qū)W習(xí)策略的改進(jìn)版,其利用精英個(gè)體和反向種群使種群多樣性加強(qiáng),搜索空間擴(kuò)大。設(shè)Xij=(Xi,1,Xi,2,...Xi,d)為海鷗種群中的精英個(gè)體,d為優(yōu)化問(wèn)題的空間維度,精英反向解可表示為:

    式(9)中,i為種群個(gè)體,h為[0,1]上的動(dòng)態(tài)系數(shù),aj和bj為動(dòng)態(tài)邊界,分別為Xij第j維個(gè)體的最小、最大值。當(dāng)超出邊界范圍時(shí),則對(duì)其進(jìn)行重置,具體表述如下:

    (3)混沌精英反向?qū)W習(xí)策略

    受文獻(xiàn)[12]啟發(fā),將Fuch映射與精英反向?qū)W習(xí)策略進(jìn)行融合,根據(jù)適應(yīng)度大小對(duì)個(gè)體進(jìn)行排序,擇優(yōu)選取形成新的初始化種群,進(jìn)而提高初始化種群質(zhì)量。具體表述如下:

    式(11)中,ki為Fuch映射中的xn,?為逐項(xiàng)乘法運(yùn)算符號(hào),為融合Fuch映射的精英反向解。融合過(guò)后,用Fuch混沌均勻變化的優(yōu)勢(shì)來(lái)動(dòng)態(tài)壓縮原先初始種群的分布范圍,讓種群更加均勻化。

    2.2 改進(jìn)自身行為特征參數(shù)A

    針對(duì)海鷗算法全局搜索線性化的問(wèn)題,提出了一種非線性搜索方法[8,14-15]。通過(guò)引入余弦函數(shù),將原先線性變化的參數(shù)A能夠非線性地變化,能夠在海鷗搜索過(guò)程中提高算法的速度和精度。

    式(12)中,t為迭代次數(shù),Maxiteration為最大迭代次數(shù)。

    改進(jìn)的參數(shù)A呈非線性變化[8,14],在迭代前期參數(shù)A變化的幅度較小、減小的速度較慢,有效地?cái)U(kuò)大了海鷗的搜索范圍;在迭代后期,參數(shù)A變化幅度較大、減小的速度較快,從而提高了算法的收斂速度。改進(jìn)前后參數(shù)A的對(duì)比圖如圖1所示。

    圖1 改進(jìn)前后參數(shù)A的對(duì)比圖Fig.1 Comparison diagram of parameter A before and after improvement

    2.3 加入萊維飛行

    萊維飛行對(duì)于海鷗算法易陷入局部最優(yōu)的問(wèn)題有一定的改進(jìn),它是通過(guò)隨機(jī)游走機(jī)制去正確控制局部搜索[13-14,20],具體表述如下:

    式(13-15)中,0<θ<=2,E、F~N(0, σ2),Γ(x)是Gamma函數(shù),α表示步長(zhǎng),θ=1.5。

    在計(jì)算萊維飛行的步長(zhǎng)[21]時(shí),可根據(jù)目標(biāo)函數(shù)的實(shí)際情況自行調(diào)整參數(shù)k,在測(cè)試函數(shù)實(shí)驗(yàn)分析中,測(cè)試函數(shù)為f6時(shí),k取0.01,I-SOA(Improved-Seagull Optimization Algorithm)的性能優(yōu)勢(shì)更為明顯,測(cè)試函數(shù)為f1-f5、f7-f9時(shí),k取0.1,I-SOA的性能優(yōu)勢(shì)更為明顯。因此,I-SOA在面對(duì)不同類(lèi)型的問(wèn)題時(shí)可通過(guò)靈活調(diào)整k值去發(fā)揮自身最大的性能優(yōu)勢(shì)。具體表述如下:(測(cè)試f6時(shí)k取0.01,其余均取0.1)

    增加萊維飛行后,海鷗群體的更新位置公式為:

    2.4 改進(jìn)海鷗算法的執(zhí)行步驟

    (1)運(yùn)用Fuch混沌精英反向?qū)W習(xí)策略對(duì)海鷗種群進(jìn)行初始化。

    (2)設(shè)置算法的相關(guān)參數(shù):種群數(shù)量、最大迭代次數(shù)和參數(shù)A中的fc值等。

    (3)計(jì)算海鷗群體的初始適應(yīng)度值,互相比較后保留最佳適應(yīng)度值和最佳位置。

    (4)根據(jù)式(12)更新A。

    (5)計(jì)算海鷗新的位置并檢查是否越界。

    (6)計(jì)算更新位置后的海鷗個(gè)體適應(yīng)度值,并與步驟(3)的最佳適應(yīng)度值進(jìn)行比較,再次更新。

    (7)算法是否達(dá)到最大迭代次數(shù),如果是,則運(yùn)行停止;如果否,則跳轉(zhuǎn)至步驟(4)繼續(xù)運(yùn)行。

    2.5 時(shí)間復(fù)雜度計(jì)算

    標(biāo)準(zhǔn)SOA算法的時(shí)間復(fù)雜度主要包括3個(gè)部分:初始化種群O(1),計(jì)算初始適應(yīng)度值O(N),迭代計(jì)算最終適應(yīng)度值 O(NT)。綜上,SOA算法的時(shí)間復(fù)雜度為O(NT)。(N為種群數(shù)量,T為最大迭代次數(shù))

    I-SOA算法的時(shí)間復(fù)雜度主要包括:初始化種群O(1),計(jì)算初始適應(yīng)度值O(N),加入混沌精英反向?qū)W習(xí)策略O(shè)(NT),改進(jìn)自身行為的特征參數(shù)O(NT),加入萊維飛行O(NT),迭代計(jì)算最終適應(yīng)度值 O(NT)。綜上,I-SOA算法的時(shí)間復(fù)雜度為O(NT)。

    由以上可知,標(biāo)準(zhǔn)SOA算法和I-SOA算法的時(shí)間復(fù)雜度相同,I-SOA算法的改進(jìn)并不會(huì)降低算法的運(yùn)行效率。

    3 測(cè)試函數(shù)實(shí)驗(yàn)分析

    3.1 測(cè)試函數(shù)

    選取9個(gè)典型的測(cè)試函數(shù)進(jìn)行性能對(duì)比。9個(gè)測(cè)試函數(shù)中有單模態(tài)的基準(zhǔn)測(cè)試函數(shù)和多模態(tài)的基準(zhǔn)測(cè)試函數(shù),能夠比較全面地檢驗(yàn)算法的尋優(yōu)性能和收斂速度,具體表述如表1所示。

    表1 測(cè)試函數(shù)Table 1 Test function

    3.2 算法測(cè)試分析

    將I-SOA算法與標(biāo)準(zhǔn)SOA、PSO和GA算法進(jìn)行對(duì)比實(shí)驗(yàn)。種群規(guī)模為50,T=500。海鷗算法中A的起始值為2;粒子群算法Vmax為6,ω為0.6,c1和c2均為2;GA算法p1=0.7,p2=0.3。實(shí)驗(yàn)重復(fù)30次,并以30次最優(yōu)解的最優(yōu)值、最差值、平均值和標(biāo)準(zhǔn)差作為算法性能評(píng)價(jià)指標(biāo)。最優(yōu)值、最差值和平均值評(píng)價(jià)指標(biāo)能夠反映各算法尋優(yōu)的收斂情況和尋優(yōu)精度。(最優(yōu)值為30次最優(yōu)解中最接近于函數(shù)最優(yōu)解的值,最差值反之,平均值為30次最優(yōu)解的算術(shù)平均值)。標(biāo)準(zhǔn)差評(píng)價(jià)指標(biāo)能夠反映算法尋優(yōu)的穩(wěn)定性。(標(biāo)準(zhǔn)差為方差的算術(shù)平方根,方差為30次最優(yōu)解與其平均值差的平方的平均數(shù))。具體結(jié)果如表2所示。

    由表2的最優(yōu)值和平均值可知,I-SOA算法在f1-f9九個(gè)基準(zhǔn)測(cè)試函數(shù)中表現(xiàn)非常優(yōu)異。在f1-f5單模態(tài)函數(shù)中,I-SOA算法求解的函數(shù)值遠(yuǎn)優(yōu)于另外三種算法,尤其在求解函數(shù)f1時(shí),I-SOA算法比標(biāo)準(zhǔn)SOA算法提升了至少30個(gè)數(shù)量級(jí),求解函數(shù)f2時(shí),提升了20個(gè)數(shù)量級(jí)以上,求解函數(shù)f3和f4時(shí),提升了8個(gè)數(shù)量級(jí)以上。對(duì)于多模態(tài)函數(shù)f6-f9而言,I-SOA算法在求解函數(shù)f7和f9時(shí),均求得函數(shù)的理論最優(yōu)值,求解函數(shù)f8時(shí),I-SOA算法比標(biāo)準(zhǔn)SOA算法提升了10個(gè)數(shù)量級(jí)以上,雖然求解函數(shù)f6時(shí)I-SOA算法沒(méi)有體現(xiàn)出過(guò)大的優(yōu)勢(shì),但總體而言還是優(yōu)于其他三種標(biāo)準(zhǔn)算法。而表2中對(duì)4種算法的最差值進(jìn)行對(duì)比時(shí),I-SOA算法的優(yōu)勢(shì)更為明顯,其在9個(gè)基準(zhǔn)測(cè)試函數(shù)中的最差值均為最小值。將四種算法對(duì)9個(gè)測(cè)試函數(shù)求解30次最優(yōu)解的標(biāo)準(zhǔn)差進(jìn)行比較,除了I-SOA算法在函數(shù)f6上表現(xiàn)略差以外,在其余的8個(gè)測(cè)試函數(shù)中,標(biāo)準(zhǔn)差均比另外三種標(biāo)準(zhǔn)算法小,說(shuō)明了I-SOA算法的穩(wěn)定性較強(qiáng),不易陷入局部最優(yōu)。

    表2 各算法優(yōu)化結(jié)果對(duì)比Table 2 Comparison of optimization results of various algorithms

    綜上所述,I-SOA算法的各方面性能均優(yōu)于其他三種算法。

    圖6 f2收斂曲線對(duì)比Fig.6 Comparison of f2 convergence curves

    圖7 f3收斂曲線對(duì)比Fig.7 Comparison of f3 convergence curves

    圖8 f4收斂曲線對(duì)比Fig.8 Comparison of f4 convergence curves

    圖11 f7收斂曲線對(duì)比Fig.11 Comparison of f7 convergence curves

    圖12 f8收斂曲線對(duì)比Fig.12 Comparison of f8 convergence curves

    圖2-圖4為加入單個(gè)策略?xún)?yōu)化的I-SOA1、I-SOA2、I-SOA3與SOA在 多 模 態(tài) 函 數(shù)Schwefel p2.26中的求解表現(xiàn)。因海鷗優(yōu)化算法多樣性較差和篇幅所限,僅選取多模態(tài)函數(shù)Schwefel p2.26作為測(cè)試函數(shù)對(duì)僅加入混沌精英反向?qū)W習(xí)策略的I-SOA1和SOA進(jìn)行測(cè)試對(duì)比,收斂曲線對(duì)比圖如圖2所示(種群規(guī)模為50,T=500,與以下兩個(gè)改進(jìn)點(diǎn)測(cè)試對(duì)比參數(shù)相同),I-SOA1的收斂速度和尋優(yōu)精度均優(yōu)于SOA,證實(shí)了加入混沌精英反向?qū)W習(xí)策略的有效性。僅改進(jìn)參數(shù)A的I-SOA2和SOA在多模態(tài)函數(shù)Schwefel p2.26的測(cè)試收斂曲線對(duì)比圖如圖3所示,I-SOA2的收斂速度和尋優(yōu)精度均優(yōu)于SOA,證實(shí)了改進(jìn)參數(shù)A的有效性。僅加入萊維飛行的I-SOA3和SOA在多模態(tài)函數(shù)Schwefel p2.26的測(cè)試收斂曲線對(duì)比圖如圖4所示,I-SOA3的尋優(yōu)精度明顯優(yōu)于SOA,跳出了局部最優(yōu)的困境,證實(shí)了加入萊維飛行的有效性(k取0.01)。

    圖2 I-SOA1和SOA收斂曲線對(duì)比圖Fig.2 Comparison of I-SOA1 and SOA convergence curves

    圖3 I-SOA2和SOA收斂曲線對(duì)比圖Fig.3 Comparison of I-SOA2 and SOA convergence curves

    圖4 I-SOA3和SOA收斂曲線對(duì)比圖Fig.4 Comparison of I-SOA3 and SOA convergence curves

    圖5-圖9為I-SOA、SOA、PSO和GA算法在單模態(tài)函數(shù)中的求解表現(xiàn)。從圖中可知,I-SOA算法在函數(shù)f1-f5中收斂曲線向x軸下降的速度更快和幅度更大,這意味著I-SOA算法在這5個(gè)測(cè)試函數(shù)中尋求最優(yōu)解的精度更高。而另外3種標(biāo)準(zhǔn)算法的迭代曲線較為平緩,甚至出現(xiàn)停滯現(xiàn)象,尋優(yōu)精度也因此變得很差,尤其是f1-f4的收斂曲線對(duì)比圖,其與I-SOA算法形成鮮明的對(duì)比。I-SOA算法在單模態(tài)函數(shù)中的求解表現(xiàn)說(shuō)明了融合Fuch混沌映射和精英反向?qū)W習(xí)策略來(lái)初始化種群增加了海鷗種群的多樣性,從而提高了種群的質(zhì)量,加快了算法的收斂速度。而引入余弦函數(shù)改進(jìn)參數(shù)A,使線性搜索非線性化,從而改變海鷗步長(zhǎng),進(jìn)一步提高了收斂速度,尋優(yōu)精度也得到一定的提升。

    圖5 f1收斂曲線對(duì)比Fig.5 Comparison of f1 convergence curves

    圖9 f5收斂曲線對(duì)比Fig.9 Comparison of f5 convergence curves

    圖10-圖13為I-SOA、SOA、PSO和GA算法在多模態(tài)函數(shù)中的求解表現(xiàn),從這4張圖中可清晰看出I-SOA算法的優(yōu)越性,其收斂曲線下降的速度明顯快于其他三種標(biāo)準(zhǔn)算法,而尋優(yōu)精度也遠(yuǎn)勝于其他三種算法,尤其在求解f7和f9時(shí),I-SOA算法在迭代次數(shù)200左右均已收斂并求出了理論最優(yōu)值。I-SOA算法在多模態(tài)函數(shù)的求解表現(xiàn)進(jìn)一步說(shuō)明了融合Fuch混沌映射和精英反向?qū)W習(xí)策略的有效性,算法的收斂速度得到有效提高,改進(jìn)的參數(shù)A也在多模態(tài)函數(shù)的測(cè)試中體現(xiàn)了自身的優(yōu)勢(shì),加快了收斂速度并提高了尋優(yōu)精度。而在求解多模態(tài)函數(shù)時(shí),標(biāo)準(zhǔn)SOA算法表現(xiàn)較弱,萊維飛行機(jī)制的引入進(jìn)一步優(yōu)化了海鷗算法,擴(kuò)大了海鷗搜索解的范圍,擺脫了容易陷入局部最優(yōu)的困境,從而求得最優(yōu)解。

    圖10 f6收斂曲線對(duì)比Fig.10 Comparison of f6 convergence curves

    圖13 f9收斂曲線對(duì)比Fig.13 Comparison of f9 convergence curves

    綜上所述,I-SOA1、I-SOA2、I-SOA3算法在多模態(tài)函數(shù)Schwefel p2.26的收斂曲線對(duì)比圖中,其求解表現(xiàn)證實(shí)了加入單個(gè)策略?xún)?yōu)化的有效性。I-SOA算法在f1-f9的收斂曲線對(duì)比圖中,其收斂曲線的下降速度比標(biāo)準(zhǔn)SOA、PSO和GA算法更快,在面對(duì)易陷入局部最優(yōu)導(dǎo)致尋優(yōu)精度低的缺陷時(shí),能夠跳出局部最優(yōu)獲得全局最優(yōu)解。通過(guò)多種性能的對(duì)比證實(shí)了I-SOA算法的優(yōu)勢(shì)更為明顯。

    4 I-SOA算法在工程設(shè)計(jì)優(yōu)化問(wèn)題中的應(yīng)用

    壓力容器設(shè)計(jì)優(yōu)化問(wèn)題、三桿桁架設(shè)計(jì)優(yōu)化問(wèn)題和拉力彈簧設(shè)計(jì)優(yōu)化問(wèn)題為非線性規(guī)劃問(wèn)題中的不等式約束問(wèn)題,群智能優(yōu)化算法在解決這類(lèi)問(wèn)題時(shí)有一定的優(yōu)勢(shì),不過(guò)仍然存在最優(yōu)解質(zhì)量較低、問(wèn)題適應(yīng)性較弱和尋優(yōu)穩(wěn)定性較差的缺陷,因此將I-SOA算法用來(lái)求解3個(gè)工程設(shè)計(jì)優(yōu)化問(wèn)題,觀察其尋優(yōu)的精度、面對(duì)不同工程設(shè)計(jì)優(yōu)化問(wèn)題的求解能力和尋求最優(yōu)解的穩(wěn)定性能夠進(jìn)一步驗(yàn)證I-SOA算法的有效性和優(yōu)越性。

    4.1 壓力容器設(shè)計(jì)優(yōu)化問(wèn)題

    在壓力容器設(shè)計(jì)[22-23]優(yōu)化問(wèn)題中,主要目的是使壓力容器的總成本最低,其總成本主要包括材料、成型和焊接成本。壓力容器的結(jié)構(gòu)示意圖如圖14所示,其有4個(gè)優(yōu)化變量:壓力容器筒體的厚度(TS)、頭部半球的厚度(Th)、半球體的內(nèi)部半徑(R)、不考慮頭部半球的圓柱筒體的長(zhǎng)度L。令{TS、Th、R、L}為{x1、x2、x3、x4},x各維度的取值范圍和壓力容器設(shè)計(jì)的目標(biāo)函數(shù)及不等式約束條件如下所示。

    圖14 壓力容器結(jié)構(gòu)示意圖Fig.14 Structural diagram of pressure vessel

    x1、x2的取值范圍為[0,100];x3、x4的取值范圍為[10,200]。

    目標(biāo)函數(shù)為:

    不等式約束條件為:

    用標(biāo)準(zhǔn)SOA算法和I-SOA算法分別求解壓力容器設(shè)計(jì)問(wèn)題(種群規(guī)模為200、迭代次數(shù)為2000,fc為2,與以下兩個(gè)工程設(shè)計(jì)問(wèn)題參數(shù)相同),各自獨(dú)立運(yùn)行50次,并以50次尋優(yōu)結(jié)果的最優(yōu)值、平均值、最差值和標(biāo)準(zhǔn)差作為算法性能評(píng)價(jià)指標(biāo),具體數(shù)值如表3所示。將標(biāo)準(zhǔn)SOA算法、I-SOA算法和文 獻(xiàn)[22]中列出的19種算法的最優(yōu)值進(jìn)行比較,具體數(shù)值如表4所示。

    表3 SOA算法和I-SOA算法求解壓力容器設(shè)計(jì)問(wèn)題的尋優(yōu)結(jié)果對(duì)照表Table 3 Comparison of optimization results of SOA algorithm and I-SOA algorithm for solving pressure vessel design problem

    表4 與文獻(xiàn)[22]中其他算法求解壓力容器設(shè)計(jì)問(wèn)題的最優(yōu)值比較Table 4 Comparison with the optimal values of other algorithms in literature [22] for solving pressure vessel design problem

    由表3中的尋優(yōu)結(jié)果和圖15的收斂曲線可清晰的看出I-SOA算法在各個(gè)方面均優(yōu)于SOA算法。表4為21種智能優(yōu)化算法在壓力容器設(shè)計(jì)問(wèn)題上求得的最優(yōu)值,I-SOA算法除了微遜于4類(lèi)灰狼算法外(相差0.0232%以下),比其余的16種算法都要優(yōu)越。表3和表4表明了I-SOA算法求解壓力容器設(shè)計(jì)問(wèn)題的優(yōu)越性。

    圖15 SOA和I-SOA求解壓力容器設(shè)計(jì)問(wèn)題的收斂曲線Fig.15 Convergence curves of SOA and I-SOA for solving pressure vessel design problem

    4.2 三桿桁架設(shè)計(jì)優(yōu)化問(wèn)題

    在三桿桁架設(shè)計(jì)[22-23]的問(wèn)題中,主要目的是使三桿桁架的體積最小,其設(shè)計(jì)示意圖如圖16所示。因三桿桁架中桿x1和x3橫截面積相同,所以只需選取x1、x2兩個(gè)桿的橫截面積作為優(yōu)化變量。σ為三桿桁架在每個(gè)桁架構(gòu)件上受到的應(yīng)力。x1和x2的取值范圍、目標(biāo)函數(shù)和不等式約束條件如下所示。

    圖16 三桿桁架結(jié)構(gòu)示意圖Fig.16 Schematic diagram of three bar truss structure

    x1、x2的取值范圍為[0,1]。

    目標(biāo)函數(shù)為:

    不等式約束條件為:

    目標(biāo)函數(shù)中L=100cm,不等式約束條件中P= 2kN/cm2,σ=2kN/cm2。

    用標(biāo)準(zhǔn)SOA算法和I-SOA算法分別求解三桿桁架設(shè)計(jì)問(wèn)題,各自獨(dú)立運(yùn)行50次,并以50次尋優(yōu)結(jié)果的最優(yōu)值、平均值、最差值和標(biāo)準(zhǔn)差作為算法性能評(píng)價(jià)指標(biāo),具體數(shù)值如表5所示。將標(biāo)準(zhǔn)SOA算法、I-SOA算法和文獻(xiàn)[22]中列出的17種算法的最優(yōu)值進(jìn)行比較,具體數(shù)值如表6所示。

    表5 SOA算法和I-SOA算法求解三桿桁架設(shè)計(jì)問(wèn)題的尋優(yōu)結(jié)果對(duì)照表Table 5 Comparison of optimization results of SOA algorithm and I-SOA algorithm for solving three bar truss design problem

    表6 與文獻(xiàn)[22]中其他算法求解三桿桁架設(shè)計(jì)問(wèn)題的最優(yōu)值比較Table 6 Comparison with the optimal values of other algorithms in literature [22] for solving three bar truss design problem

    由表5中的尋優(yōu)結(jié)果與圖17中的收斂曲線可知,I-SOA算法在求解三桿桁架設(shè)計(jì)的問(wèn)題上雖然優(yōu)勢(shì)不明顯,但相比于比SOA算法依舊有微弱的優(yōu)勢(shì)。表6為SOA、I-SOA算法與文獻(xiàn)[22]中其他17種算法的最優(yōu)值進(jìn)行比較,在三桿桁架設(shè)計(jì)的問(wèn)題上,群智能優(yōu)化算法的表現(xiàn)都很不錯(cuò),I-SOA算法求解的最優(yōu)值在19種算法中位居前列,占據(jù)微弱的優(yōu)勢(shì)。表5、表6表明了I-SOA算法在求解三桿桁架設(shè)計(jì)問(wèn)題時(shí)有一定的優(yōu)越性。

    圖17 SOA和I-SOA求解三桿桁架設(shè)計(jì)問(wèn)題的收斂曲線Fig.17 Convergence curves of SOA and I-SOA for solving three bar truss design problem

    4.3 拉力彈簧設(shè)計(jì)優(yōu)化問(wèn)題

    在拉力彈簧的設(shè)計(jì)[24]問(wèn)題中,在滿(mǎn)足最小撓度、剪應(yīng)力和振動(dòng)頻率等約束的條件下使拉力彈簧的質(zhì)量最小,其結(jié)構(gòu)示意圖如圖18所示。設(shè)計(jì)變量主要有3個(gè),分別是彈簧線圈的直徑d(x1)、彈簧圈平均直徑D(x2)、有效繞圈的數(shù)量P(x3)。x1-x3的取值范圍、目標(biāo)函數(shù)和不等式約束條件如下所示。

    圖18 拉力彈簧結(jié)構(gòu)示意圖Fig.18 Structural diagram of tension spring

    x1的取值范圍為[0.25,1.30];x2的取值范圍為[0.05,2.00];x3的取值范圍為[2.00,15.00]。

    目標(biāo)函數(shù)為:

    不等式約束條件為:

    用標(biāo)準(zhǔn)SOA算法和I-SOA算法分別求解拉力彈簧設(shè)計(jì)問(wèn)題,各自獨(dú)立運(yùn)行30次,并將30次尋優(yōu)結(jié)果的最優(yōu)值、平均值、最差值和標(biāo)準(zhǔn)差與文獻(xiàn)[24] 4種算法進(jìn)行對(duì)比,具體數(shù)值如表7所示。

    表7 與文獻(xiàn)[24]中其他算法求解拉力彈簧設(shè)計(jì)問(wèn)題尋優(yōu)結(jié)果對(duì)照表Table 7 Comparison of optimization results of other algorithms in literature [24] for solving tension spring design problem

    由表7的尋優(yōu)結(jié)果可知,I-SOA算法在求解拉力彈簧設(shè)計(jì)問(wèn)題上優(yōu)勢(shì)明顯,其在與SOA算法和文獻(xiàn)[24]中其余4種算法的尋優(yōu)結(jié)果對(duì)比時(shí),最優(yōu)值最小,平均值和最差值僅次于CDE算法,標(biāo)準(zhǔn)差較CDE、AATM、IFA算法稍顯弱勢(shì)。綜合而言,I-SOA算法的性能優(yōu)勢(shì)最大。由圖19可知,相比于標(biāo)準(zhǔn)SOA算法,I-SOA算法在求解拉力彈簧問(wèn)題上收斂速度更快,尋優(yōu)結(jié)果更好。表7、圖19表明了I-SOA算法在求解拉力彈簧設(shè)計(jì)問(wèn)題上有較好的性能,占據(jù)一定的優(yōu)勢(shì)。

    圖19 SOA和I-SOA求解拉力彈簧設(shè)計(jì)問(wèn)題的收斂曲線Fig.19 Convergence curve of SOA and I-SOA for solving tension spring design problem

    以上3種工程設(shè)計(jì)優(yōu)化問(wèn)題的測(cè)試結(jié)果、文獻(xiàn)對(duì)比結(jié)果和收斂曲線對(duì)比圖顯示出SOA算法收斂速度慢和易陷入局部最優(yōu)的問(wèn)題,導(dǎo)致其尋優(yōu)精度較差。I-SOA算法相比SOA算法尋優(yōu)性能較好且收斂速度快,能夠跳出局部最優(yōu)從而找到更優(yōu)解,并且相比于其他智能優(yōu)化算法,I-SOA算法在3個(gè)不同的工程設(shè)計(jì)優(yōu)化問(wèn)題中均有不錯(cuò)的表現(xiàn),體現(xiàn)了I-SOA算法的普適性和尋優(yōu)穩(wěn)定性更好。I-SOA算法在工程設(shè)計(jì)優(yōu)化問(wèn)題中的表現(xiàn)體現(xiàn)了算法改進(jìn)的有效性和優(yōu)越性,也進(jìn)一步說(shuō)明了I-SOA算法在工程設(shè)計(jì)領(lǐng)域中有較強(qiáng)的應(yīng)用價(jià)值。

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

    針對(duì)海鷗算法收斂速度慢和易陷入局部最優(yōu)等問(wèn)題,提出了一種融合Fuch混沌精英反向?qū)W習(xí)策略和余弦函數(shù)的萊維飛行海鷗算法。并將改進(jìn)的海鷗算法與標(biāo)準(zhǔn)SOA、PSO和GA算法在9個(gè)基準(zhǔn)測(cè)試函數(shù)和3個(gè)工程設(shè)計(jì)優(yōu)化問(wèn)題上進(jìn)行性能對(duì)比,得出以下結(jié)論:

    (1)改進(jìn)的海鷗算法無(wú)論在尋求最優(yōu)解還是收斂速度上均有明顯的提升,說(shuō)明了改進(jìn)的有效性。首先Fuch混沌精英反向?qū)W習(xí)策略提高了種群的質(zhì)量,其次引入余弦函數(shù)對(duì)參數(shù)A進(jìn)行改進(jìn)提高了算法的收斂速度和全局尋優(yōu)能力,最后萊維飛行機(jī)制的加入使算法能夠跳出局部最優(yōu)。

    (2)I-SOA算法在9個(gè)基準(zhǔn)測(cè)試函數(shù)中收斂速度、尋優(yōu)精度和跳出局部最優(yōu)的能力均比標(biāo)準(zhǔn)SOA、PSO和GA算法優(yōu)越,證明了I-SOA算法改進(jìn)的有效性。

    (3)I-SOA算法在3個(gè)工程設(shè)計(jì)優(yōu)化問(wèn)題中,尋優(yōu)精度和跳出局部最優(yōu)的能力較強(qiáng),面對(duì)不同的工程設(shè)計(jì)優(yōu)化問(wèn)題均能求得更優(yōu)質(zhì)的解,體現(xiàn)了I-SOA算法的問(wèn)題適應(yīng)性和求解穩(wěn)定性更優(yōu),證明了I-SOA算法改進(jìn)的有效性。

    在之后的研究中,會(huì)繼續(xù)完善I-SOA算法,并將I-SOA算法應(yīng)用到更廣泛的實(shí)際問(wèn)題中去,提高I-SOA算法的應(yīng)用價(jià)值。

    利益沖突聲明

    所有作者聲明不存在利益沖突關(guān)系。

    猜你喜歡
    測(cè)試函數(shù)海鷗學(xué)習(xí)策略
    霸道海鷗誰(shuí)能治
    “海鷗”展翅 “美好”起飛
    具有收縮因子的自適應(yīng)鴿群算法用于函數(shù)優(yōu)化問(wèn)題
    高中生數(shù)學(xué)自主學(xué)習(xí)策略探討
    帶勢(shì)函數(shù)的雙調(diào)和不等式組的整體解的不存在性
    約束二進(jìn)制二次規(guī)劃測(cè)試函數(shù)的一個(gè)構(gòu)造方法
    約束二進(jìn)制二次規(guī)劃測(cè)試函數(shù)的一個(gè)構(gòu)造方法
    一種使用反向?qū)W習(xí)策略的改進(jìn)花粉授粉算法
    面向真實(shí)世界的測(cè)試函數(shù)Ⅱ
    基于微博的移動(dòng)學(xué)習(xí)策略研究
    国产一级毛片七仙女欲春2| 亚洲av电影在线进入| 天堂av国产一区二区熟女人妻| 老鸭窝网址在线观看| 女人十人毛片免费观看3o分钟| 国产极品精品免费视频能看的| 性欧美人与动物交配| 亚洲无线在线观看| 老女人水多毛片| 中文字幕免费在线视频6| 中文字幕熟女人妻在线| 国产高潮美女av| 国产久久久一区二区三区| 国产精品一及| 亚洲欧美日韩高清在线视频| 97热精品久久久久久| 一级a爱片免费观看的视频| 成人鲁丝片一二三区免费| 人人妻,人人澡人人爽秒播| 两人在一起打扑克的视频| 欧美一区二区国产精品久久精品| 九九在线视频观看精品| 亚洲第一电影网av| 两个人视频免费观看高清| 亚洲成人精品中文字幕电影| 精品国内亚洲2022精品成人| 免费人成视频x8x8入口观看| 婷婷精品国产亚洲av在线| 桃红色精品国产亚洲av| 国产精品乱码一区二三区的特点| 国产美女午夜福利| 国产亚洲欧美在线一区二区| 12—13女人毛片做爰片一| 淫妇啪啪啪对白视频| 日韩有码中文字幕| 琪琪午夜伦伦电影理论片6080| ponron亚洲| 午夜精品一区二区三区免费看| 亚洲av免费高清在线观看| 丰满人妻一区二区三区视频av| 一级黄色大片毛片| 麻豆国产av国片精品| 床上黄色一级片| 久久久久免费精品人妻一区二区| 非洲黑人性xxxx精品又粗又长| 99热只有精品国产| 国产精品久久久久久久电影| 亚洲精品乱码久久久v下载方式| 少妇裸体淫交视频免费看高清| 午夜福利高清视频| 亚洲在线观看片| 日本免费a在线| 亚洲国产精品999在线| 欧美+日韩+精品| 能在线免费观看的黄片| 成人特级av手机在线观看| 色吧在线观看| 中文字幕人成人乱码亚洲影| av在线蜜桃| 白带黄色成豆腐渣| 国产熟女xx| av专区在线播放| 中文在线观看免费www的网站| 久久99热6这里只有精品| av视频在线观看入口| 亚洲美女视频黄频| 久久精品国产亚洲av涩爱 | www.999成人在线观看| 成人av一区二区三区在线看| 99热这里只有是精品50| 欧美成狂野欧美在线观看| 亚洲精品粉嫩美女一区| 97碰自拍视频| 亚洲欧美日韩东京热| 成年女人毛片免费观看观看9| 免费搜索国产男女视频| 热99在线观看视频| 欧美日韩综合久久久久久 | 亚洲欧美精品综合久久99| 免费观看人在逋| 亚洲片人在线观看| 在线观看66精品国产| 99热这里只有是精品50| 可以在线观看的亚洲视频| a级毛片a级免费在线| 男人舔奶头视频| 日韩国内少妇激情av| 蜜桃亚洲精品一区二区三区| 中国美女看黄片| 久久精品影院6| 午夜影院日韩av| 波野结衣二区三区在线| 久久国产精品人妻蜜桃| 欧美成狂野欧美在线观看| 亚洲av电影在线进入| av国产免费在线观看| h日本视频在线播放| 极品教师在线视频| 哪里可以看免费的av片| 午夜日韩欧美国产| 国产精品伦人一区二区| 国产私拍福利视频在线观看| 色哟哟哟哟哟哟| 成人精品一区二区免费| 久久午夜福利片| av欧美777| 激情在线观看视频在线高清| 精品人妻一区二区三区麻豆 | 天堂动漫精品| 人妻夜夜爽99麻豆av| 午夜免费男女啪啪视频观看 | 熟女电影av网| 激情在线观看视频在线高清| 两性午夜刺激爽爽歪歪视频在线观看| 亚洲午夜理论影院| 小蜜桃在线观看免费完整版高清| 精品久久久久久久人妻蜜臀av| 国产激情偷乱视频一区二区| 免费人成在线观看视频色| 中文字幕久久专区| 国产成人a区在线观看| 十八禁国产超污无遮挡网站| 亚洲人成网站高清观看| 最近最新中文字幕大全电影3| 高清在线国产一区| 高清在线国产一区| 中文字幕人妻熟人妻熟丝袜美| 亚洲第一电影网av| 真实男女啪啪啪动态图| 久久国产精品人妻蜜桃| 久久久精品大字幕| 国产色婷婷99| 99精品在免费线老司机午夜| 窝窝影院91人妻| 国产伦精品一区二区三区四那| 色综合婷婷激情| 欧美乱色亚洲激情| 国产探花极品一区二区| 少妇熟女aⅴ在线视频| 在线十欧美十亚洲十日本专区| 在线十欧美十亚洲十日本专区| 看十八女毛片水多多多| 在线观看66精品国产| 亚洲国产高清在线一区二区三| 精品一区二区三区人妻视频| 国产精品综合久久久久久久免费| 国产爱豆传媒在线观看| av中文乱码字幕在线| 午夜福利18| 午夜福利欧美成人| 九九久久精品国产亚洲av麻豆| 久久婷婷人人爽人人干人人爱| 国产不卡一卡二| 久久精品国产清高在天天线| 色吧在线观看| 免费高清视频大片| 欧美性猛交黑人性爽| 99国产综合亚洲精品| 天美传媒精品一区二区| 国产高清有码在线观看视频| 中国美女看黄片| 久久久色成人| 中文字幕精品亚洲无线码一区| 国产成人欧美在线观看| 精品99又大又爽又粗少妇毛片 | 亚洲美女黄片视频| 久久精品91蜜桃| 国内精品一区二区在线观看| 一个人看的www免费观看视频| 国产亚洲精品久久久久久毛片| 国内精品美女久久久久久| 免费在线观看影片大全网站| 久久久久亚洲av毛片大全| 精品人妻视频免费看| 国产不卡一卡二| 日本与韩国留学比较| 国产不卡一卡二| 亚洲av美国av| 最近视频中文字幕2019在线8| 一个人免费在线观看电影| 99久久无色码亚洲精品果冻| 欧美黑人巨大hd| 麻豆国产av国片精品| 美女高潮的动态| 精品久久久久久久久久久久久| 脱女人内裤的视频| 亚洲,欧美,日韩| 亚洲av二区三区四区| 国产午夜福利久久久久久| 国产精品影院久久| 人人妻人人澡欧美一区二区| 午夜福利欧美成人| 精品午夜福利在线看| or卡值多少钱| 男女视频在线观看网站免费| 18禁裸乳无遮挡免费网站照片| 国产亚洲精品久久久com| 欧美色视频一区免费| 国产精品亚洲av一区麻豆| 高清在线国产一区| 国产伦精品一区二区三区四那| 欧美日韩瑟瑟在线播放| 成年女人看的毛片在线观看| 嫩草影视91久久| 日本 av在线| 女生性感内裤真人,穿戴方法视频| 午夜亚洲福利在线播放| 美女cb高潮喷水在线观看| 欧美黄色淫秽网站| www.熟女人妻精品国产| 国产在线男女| 精品福利观看| 91字幕亚洲| 久久精品国产亚洲av香蕉五月| 国产一区二区亚洲精品在线观看| 69av精品久久久久久| 久久午夜亚洲精品久久| 久久久国产成人精品二区| 国产精品一区二区免费欧美| 在线播放无遮挡| 日韩中文字幕欧美一区二区| 神马国产精品三级电影在线观看| 男女下面进入的视频免费午夜| 麻豆国产97在线/欧美| 欧美成人性av电影在线观看| 欧美色视频一区免费| 精品午夜福利视频在线观看一区| 午夜精品一区二区三区免费看| 国产午夜精品论理片| 国产午夜福利久久久久久| 美女免费视频网站| 欧美精品国产亚洲| 国产高清有码在线观看视频| av女优亚洲男人天堂| 欧洲精品卡2卡3卡4卡5卡区| 波多野结衣巨乳人妻| 小说图片视频综合网站| 国产三级中文精品| 国产私拍福利视频在线观看| 偷拍熟女少妇极品色| 亚洲av免费高清在线观看| 色综合欧美亚洲国产小说| 免费一级毛片在线播放高清视频| 国产精品一区二区三区四区免费观看 | 精品视频人人做人人爽| 午夜免费观看性视频| 成人免费观看视频高清| 国产黄色免费在线视频| 一本色道久久久久久精品综合| 精品少妇久久久久久888优播| 亚洲欧美清纯卡通| 欧美 日韩 精品 国产| 成人午夜精彩视频在线观看| 99精国产麻豆久久婷婷| 久久99蜜桃精品久久| 国产av不卡久久| 欧美xxxx黑人xx丫x性爽| 国产精品一二三区在线看| 国产色婷婷99| av福利片在线观看| 国产真实伦视频高清在线观看| 日韩av免费高清视频| 91精品伊人久久大香线蕉| 免费黄色在线免费观看| 亚洲人成网站高清观看| 看免费成人av毛片| 日韩av不卡免费在线播放| 高清av免费在线| 国产视频内射| 日本熟妇午夜| 只有这里有精品99| 69av精品久久久久久| 精品视频人人做人人爽| 精品久久久久久久久亚洲| 亚洲精品日本国产第一区| 欧美人与善性xxx| 大码成人一级视频| 干丝袜人妻中文字幕| www.色视频.com| 天天躁日日操中文字幕| 国产成人精品一,二区| 亚洲精品自拍成人| 日韩电影二区| 久久ye,这里只有精品| 久久韩国三级中文字幕| 51国产日韩欧美| 最近最新中文字幕大全电影3| 伊人久久精品亚洲午夜| 久久久久久九九精品二区国产| 国产精品久久久久久精品电影| 激情 狠狠 欧美| 久久久色成人| 又黄又爽又刺激的免费视频.| 十八禁网站网址无遮挡 | 国产乱来视频区| 精品久久久久久久人妻蜜臀av| 欧美三级亚洲精品| 色吧在线观看| 国产免费一区二区三区四区乱码| 干丝袜人妻中文字幕| 日韩一区二区视频免费看| 人体艺术视频欧美日本| 国产真实伦视频高清在线观看| 麻豆国产97在线/欧美| 六月丁香七月| 国产精品国产三级国产专区5o| 男女那种视频在线观看| 精品少妇黑人巨大在线播放| 日韩欧美 国产精品| 精品久久久久久电影网| 18禁在线播放成人免费| 九九爱精品视频在线观看| 人人妻人人爽人人添夜夜欢视频 | 日本黄大片高清| 少妇丰满av| 欧美日韩视频高清一区二区三区二| 久热久热在线精品观看| 麻豆成人av视频| 最新中文字幕久久久久| 嫩草影院入口| 国产综合懂色| 午夜福利在线观看免费完整高清在| 麻豆国产97在线/欧美| 国产精品一区www在线观看| 精品亚洲乱码少妇综合久久| 高清午夜精品一区二区三区| 又粗又硬又长又爽又黄的视频| 亚洲欧美日韩另类电影网站 | 极品教师在线视频| 亚洲经典国产精华液单| 97超视频在线观看视频| 麻豆成人午夜福利视频| 欧美变态另类bdsm刘玥| 只有这里有精品99| 国产高清国产精品国产三级 | av国产精品久久久久影院| 日韩一区二区三区影片| 国产一区有黄有色的免费视频| 高清午夜精品一区二区三区| 日本wwww免费看| 春色校园在线视频观看| 中文字幕亚洲精品专区| 成人免费观看视频高清| av免费观看日本| 18禁在线播放成人免费| 三级国产精品片| 国产成人一区二区在线| 色吧在线观看| 久久精品夜色国产| kizo精华| 亚洲成人中文字幕在线播放| 亚洲aⅴ乱码一区二区在线播放| 别揉我奶头 嗯啊视频| 久久精品人妻少妇| 精品久久久久久电影网| 国产精品熟女久久久久浪| 亚洲性久久影院| 国产精品国产三级国产av玫瑰| 丝袜脚勾引网站| 尾随美女入室| 精品一区二区免费观看| 人妻少妇偷人精品九色| 在线播放无遮挡| 国产精品一区二区性色av| 黄色怎么调成土黄色| 国产在线一区二区三区精| 大陆偷拍与自拍| 成人国产av品久久久| 久久久久久久大尺度免费视频| 亚洲性久久影院| 纵有疾风起免费观看全集完整版| 亚洲av.av天堂| 黄色怎么调成土黄色| 国产黄色视频一区二区在线观看| 免费少妇av软件| 熟女电影av网| 视频区图区小说| 国语对白做爰xxxⅹ性视频网站| 欧美成人午夜免费资源| 免费大片黄手机在线观看| 久久国产乱子免费精品| 深爱激情五月婷婷| 少妇熟女欧美另类| 久热这里只有精品99| 两个人的视频大全免费| 七月丁香在线播放| 国产精品久久久久久av不卡| 欧美 日韩 精品 国产| 色5月婷婷丁香| 亚洲av中文av极速乱| 少妇的逼好多水| 国产人妻一区二区三区在| 国产色婷婷99| av.在线天堂| 国产精品久久久久久精品电影| 又粗又硬又长又爽又黄的视频| 毛片女人毛片| 在线观看一区二区三区激情| 欧美最新免费一区二区三区| 男女下面进入的视频免费午夜| 在线观看美女被高潮喷水网站| 国产乱来视频区| 18+在线观看网站| 午夜福利在线在线| 深爱激情五月婷婷| 一级av片app| 看免费成人av毛片| 欧美日韩精品成人综合77777| 欧美日韩视频高清一区二区三区二| 女人被狂操c到高潮| 伦精品一区二区三区| 婷婷色综合大香蕉| 91久久精品国产一区二区成人| 蜜桃久久精品国产亚洲av| 亚洲经典国产精华液单| 精品久久久久久电影网| 一级毛片我不卡| 涩涩av久久男人的天堂| 国产午夜福利久久久久久| 亚洲av福利一区| 最近2019中文字幕mv第一页| 国产伦精品一区二区三区四那| 亚洲精品亚洲一区二区| 中国国产av一级| 男人狂女人下面高潮的视频| 国产一区二区三区综合在线观看 | 男男h啪啪无遮挡| 国产乱人偷精品视频| 插逼视频在线观看| 国产精品人妻久久久影院| 久久久色成人| 国产在线男女| 欧美亚洲 丝袜 人妻 在线| 欧美一区二区亚洲| 午夜精品国产一区二区电影 | 日韩精品有码人妻一区| 2021少妇久久久久久久久久久| 欧美变态另类bdsm刘玥| 国产真实伦视频高清在线观看| 国模一区二区三区四区视频| 男女无遮挡免费网站观看| 久久99蜜桃精品久久| 国产免费一区二区三区四区乱码| 国产成人精品一,二区| 欧美日韩在线观看h| 国语对白做爰xxxⅹ性视频网站| 内射极品少妇av片p| 黑人高潮一二区| 久久99热6这里只有精品| 人体艺术视频欧美日本| a级毛片免费高清观看在线播放| 又爽又黄无遮挡网站| 一级毛片黄色毛片免费观看视频| 亚洲欧洲日产国产| 亚洲久久久久久中文字幕| 国产午夜福利久久久久久| 91狼人影院| 神马国产精品三级电影在线观看| 国产免费视频播放在线视频| 久热这里只有精品99| 男女下面进入的视频免费午夜| 好男人在线观看高清免费视频| www.av在线官网国产| 国产精品三级大全| 国产爽快片一区二区三区| 日韩人妻高清精品专区| 一区二区三区免费毛片| 国产伦理片在线播放av一区| 少妇人妻久久综合中文| 狂野欧美白嫩少妇大欣赏| 国产毛片a区久久久久| 欧美少妇被猛烈插入视频| 亚洲久久久久久中文字幕| 男人舔奶头视频| 亚洲精品国产av成人精品| 天天一区二区日本电影三级| 麻豆国产97在线/欧美| 欧美高清性xxxxhd video| 日韩一本色道免费dvd| 六月丁香七月| 色5月婷婷丁香| 免费不卡的大黄色大毛片视频在线观看| 国产精品国产三级国产专区5o| 熟女电影av网| 91在线精品国自产拍蜜月| 久久久精品免费免费高清| av免费在线看不卡| 特大巨黑吊av在线直播| 中文天堂在线官网| 大片免费播放器 马上看| 欧美日本视频| 日韩制服骚丝袜av| 国产亚洲一区二区精品| 永久免费av网站大全| 毛片一级片免费看久久久久| av在线亚洲专区| 日韩三级伦理在线观看| 内地一区二区视频在线| 日本与韩国留学比较| 老司机影院成人| 一个人看的www免费观看视频| 久久ye,这里只有精品| 黄色欧美视频在线观看| 亚洲熟女精品中文字幕| 久久久久国产精品人妻一区二区| 精品国产一区二区三区久久久樱花 | 国产成人午夜福利电影在线观看| 黑人高潮一二区| 国产在线一区二区三区精| 国产 一区精品| 麻豆乱淫一区二区| 两个人的视频大全免费| 国语对白做爰xxxⅹ性视频网站| 岛国毛片在线播放| 超碰av人人做人人爽久久| 秋霞伦理黄片| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 国产成人精品一,二区| 少妇人妻 视频| 美女视频免费永久观看网站| 亚洲av中文av极速乱| 国产精品久久久久久精品电影小说 | 18禁动态无遮挡网站| 久久99蜜桃精品久久| 自拍欧美九色日韩亚洲蝌蚪91 | www.色视频.com| 18禁动态无遮挡网站| 黄片无遮挡物在线观看| 亚洲精品第二区| av网站免费在线观看视频| 亚洲av免费高清在线观看| 欧美日韩综合久久久久久| 最近中文字幕高清免费大全6| 午夜免费男女啪啪视频观看| 国产av不卡久久| 观看免费一级毛片| 亚洲综合色惰| 亚洲精品一区蜜桃| 亚洲美女视频黄频| 国产v大片淫在线免费观看| 国产精品国产三级专区第一集| 麻豆成人午夜福利视频| 亚洲欧美精品自产自拍| 欧美日韩国产mv在线观看视频 | 看免费成人av毛片| 校园人妻丝袜中文字幕| 欧美+日韩+精品| 久久综合国产亚洲精品| 亚洲国产精品999| 亚洲av电影在线观看一区二区三区 | 精品国产一区二区三区久久久樱花 | 夫妻性生交免费视频一级片| 91久久精品国产一区二区成人| 欧美+日韩+精品| 在线观看美女被高潮喷水网站| 人妻一区二区av| 久久综合国产亚洲精品| 免费看不卡的av| 久久影院123| 中国美白少妇内射xxxbb| 青春草亚洲视频在线观看| 国产精品不卡视频一区二区| 久久影院123| 亚洲性久久影院| 视频中文字幕在线观看| 色5月婷婷丁香| 午夜免费鲁丝| 久久鲁丝午夜福利片| 女人十人毛片免费观看3o分钟| www.av在线官网国产| 在线观看免费高清a一片| 国产黄片美女视频| 91精品伊人久久大香线蕉| 免费av不卡在线播放| 精品久久久精品久久久| 欧美日韩一区二区视频在线观看视频在线 | 午夜福利视频精品| 人人妻人人澡人人爽人人夜夜| 王馨瑶露胸无遮挡在线观看| 人体艺术视频欧美日本| 高清欧美精品videossex| 寂寞人妻少妇视频99o| 大香蕉97超碰在线| 99re6热这里在线精品视频| 91在线精品国自产拍蜜月| 国产精品一区二区三区四区免费观看| 欧美zozozo另类| 欧美最新免费一区二区三区| 精品久久久精品久久久| 日本-黄色视频高清免费观看| 2021少妇久久久久久久久久久| 1000部很黄的大片| 亚洲人与动物交配视频| 国产亚洲5aaaaa淫片| 亚洲精品中文字幕在线视频 | 国产精品嫩草影院av在线观看| 亚洲精品aⅴ在线观看| 午夜激情福利司机影院| 国产精品偷伦视频观看了| 欧美日韩亚洲高清精品| 亚洲国产av新网站| tube8黄色片| 夫妻午夜视频| 亚洲精品一二三| 最近中文字幕高清免费大全6| 久久久久久久久久久免费av| 国语对白做爰xxxⅹ性视频网站| 国产成人精品久久久久久| 亚洲人与动物交配视频| 街头女战士在线观看网站| 国产91av在线免费观看| 国产精品一及| 成人欧美大片| 超碰av人人做人人爽久久| 看十八女毛片水多多多|