• <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í)策略研究
    久久精品久久精品一区二区三区| 热99久久久久精品小说推荐| 伊人久久大香线蕉亚洲五| a级毛片黄视频| 一本久久精品| 黄色毛片三级朝国网站| 国产精品欧美亚洲77777| 在线观看免费午夜福利视频| 国产精品一国产av| 中文字幕另类日韩欧美亚洲嫩草| 夫妻性生交免费视频一级片| 人人妻人人澡人人爽人人夜夜| 老鸭窝网址在线观看| 亚洲人成电影免费在线| 午夜激情久久久久久久| 免费高清在线观看视频在线观看| 国产精品成人在线| 亚洲精品久久久久久婷婷小说| 一本色道久久久久久精品综合| 中文字幕最新亚洲高清| 中文字幕精品免费在线观看视频| 亚洲精品在线美女| 亚洲欧美一区二区三区国产| 国产黄频视频在线观看| 亚洲熟女精品中文字幕| 别揉我奶头~嗯~啊~动态视频 | cao死你这个sao货| 欧美在线一区亚洲| 欧美日韩亚洲综合一区二区三区_| 黄色一级大片看看| 久久性视频一级片| 国产视频首页在线观看| 少妇 在线观看| 在线观看免费视频网站a站| 亚洲七黄色美女视频| 亚洲免费av在线视频| 99久久综合免费| 好男人电影高清在线观看| 老汉色av国产亚洲站长工具| 黄色片一级片一级黄色片| 精品少妇久久久久久888优播| 午夜两性在线视频| 欧美变态另类bdsm刘玥| 只有这里有精品99| 丁香六月欧美| 精品亚洲乱码少妇综合久久| 国产精品av久久久久免费| 精品久久久久久电影网| 91九色精品人成在线观看| 99热全是精品| 狠狠精品人妻久久久久久综合| 中文字幕制服av| 欧美av亚洲av综合av国产av| 久久精品aⅴ一区二区三区四区| 成年动漫av网址| 亚洲第一av免费看| 欧美在线黄色| 久久精品aⅴ一区二区三区四区| 久久av网站| 欧美激情极品国产一区二区三区| 亚洲成人国产一区在线观看 | 欧美日韩亚洲国产一区二区在线观看 | 性色av乱码一区二区三区2| 午夜免费男女啪啪视频观看| 亚洲少妇的诱惑av| 午夜福利视频在线观看免费| 日本欧美视频一区| 老司机靠b影院| 激情视频va一区二区三区| 伊人久久大香线蕉亚洲五| 国产主播在线观看一区二区 | 在线看a的网站| 亚洲成色77777| 看十八女毛片水多多多| 欧美xxⅹ黑人| 欧美日韩亚洲综合一区二区三区_| 成人午夜精彩视频在线观看| 男女之事视频高清在线观看 | 免费久久久久久久精品成人欧美视频| 美女扒开内裤让男人捅视频| 国产高清videossex| 精品久久蜜臀av无| 久久九九热精品免费| 只有这里有精品99| 中文字幕人妻熟女乱码| 午夜老司机福利片| 各种免费的搞黄视频| 一本色道久久久久久精品综合| 一级毛片女人18水好多 | 九草在线视频观看| 91精品三级在线观看| 日本黄色日本黄色录像| 黄色一级大片看看| 国产男女超爽视频在线观看| 亚洲精品自拍成人| 欧美人与性动交α欧美精品济南到| 极品人妻少妇av视频| 亚洲一区二区三区欧美精品| 18在线观看网站| 亚洲美女黄色视频免费看| 婷婷色综合大香蕉| 高清av免费在线| 国精品久久久久久国模美| 亚洲国产欧美网| 国产1区2区3区精品| 丁香六月欧美| 99久久综合免费| 在现免费观看毛片| 国产精品免费视频内射| 999精品在线视频| 精品亚洲成国产av| 欧美黑人欧美精品刺激| 久久人妻熟女aⅴ| 久久免费观看电影| 国产精品av久久久久免费| 叶爱在线成人免费视频播放| 精品国产乱码久久久久久男人| 伊人亚洲综合成人网| 18在线观看网站| 一级毛片电影观看| 欧美成人精品欧美一级黄| 19禁男女啪啪无遮挡网站| 欧美日韩成人在线一区二区| 久久久精品国产亚洲av高清涩受| 亚洲精品久久久久久婷婷小说| 国产精品成人在线| 成人黄色视频免费在线看| 国产老妇伦熟女老妇高清| 欧美乱码精品一区二区三区| 午夜影院在线不卡| 男女边摸边吃奶| 久久久精品94久久精品| 亚洲中文日韩欧美视频| 免费在线观看视频国产中文字幕亚洲 | 欧美精品啪啪一区二区三区 | 久久99热这里只频精品6学生| 欧美 日韩 精品 国产| 人人妻人人澡人人爽人人夜夜| 亚洲av成人精品一二三区| 校园人妻丝袜中文字幕| 老司机深夜福利视频在线观看 | 国产一区二区 视频在线| 国产精品人妻久久久影院| 国产成人av激情在线播放| 精品卡一卡二卡四卡免费| 欧美人与性动交α欧美软件| 亚洲少妇的诱惑av| 电影成人av| 高潮久久久久久久久久久不卡| 中文字幕制服av| 人体艺术视频欧美日本| 久久国产精品人妻蜜桃| 建设人人有责人人尽责人人享有的| 男人添女人高潮全过程视频| 精品少妇黑人巨大在线播放| 久热这里只有精品99| 国产精品麻豆人妻色哟哟久久| 一个人免费看片子| 国产91精品成人一区二区三区 | 日本欧美视频一区| 男女边摸边吃奶| √禁漫天堂资源中文www| 一级毛片女人18水好多 | 精品一区二区三卡| 乱人伦中国视频| 久久精品aⅴ一区二区三区四区| 国产免费现黄频在线看| 天天躁夜夜躁狠狠久久av| 一本色道久久久久久精品综合| a级毛片黄视频| 国产欧美日韩一区二区三 | h视频一区二区三区| 国产精品国产av在线观看| 欧美黑人欧美精品刺激| 午夜免费男女啪啪视频观看| 亚洲人成77777在线视频| 亚洲九九香蕉| 国产福利在线免费观看视频| 午夜老司机福利片| 亚洲精品日本国产第一区| h视频一区二区三区| 人人妻人人澡人人爽人人夜夜| 亚洲精品一二三| 亚洲国产最新在线播放| 亚洲精品第二区| 中文字幕人妻丝袜一区二区| 中文精品一卡2卡3卡4更新| 又紧又爽又黄一区二区| 欧美日韩综合久久久久久| 亚洲综合色网址| 男人爽女人下面视频在线观看| 国产真人三级小视频在线观看| 亚洲国产日韩一区二区| 视频区欧美日本亚洲| 亚洲熟女精品中文字幕| 久久综合国产亚洲精品| 老司机亚洲免费影院| 欧美成人午夜精品| 999精品在线视频| avwww免费| 国产伦理片在线播放av一区| 51午夜福利影视在线观看| 美女扒开内裤让男人捅视频| 免费看av在线观看网站| 亚洲精品美女久久av网站| 嫁个100分男人电影在线观看 | 国产成人啪精品午夜网站| 日本猛色少妇xxxxx猛交久久| 久久青草综合色| 午夜视频精品福利| 国产欧美日韩一区二区三 | 国产午夜精品一二区理论片| 久久这里只有精品19| 一级黄色大片毛片| 国产熟女欧美一区二区| 欧美乱码精品一区二区三区| 午夜91福利影院| 侵犯人妻中文字幕一二三四区| 999久久久国产精品视频| 美女中出高潮动态图| 女性生殖器流出的白浆| 免费看十八禁软件| 欧美精品av麻豆av| 亚洲成人免费av在线播放| 婷婷色综合大香蕉| 午夜激情av网站| 欧美老熟妇乱子伦牲交| 在线观看免费日韩欧美大片| 欧美久久黑人一区二区| 婷婷色综合www| 丝袜在线中文字幕| 80岁老熟妇乱子伦牲交| 黄色一级大片看看| 欧美亚洲 丝袜 人妻 在线| 久久ye,这里只有精品| 在线 av 中文字幕| 精品一区二区三卡| 制服诱惑二区| 日韩 亚洲 欧美在线| 亚洲av电影在线进入| 国产高清国产精品国产三级| 精品视频人人做人人爽| 女人精品久久久久毛片| 三上悠亚av全集在线观看| 亚洲欧美精品自产自拍| www.熟女人妻精品国产| 久久久精品区二区三区| www.av在线官网国产| 免费人妻精品一区二区三区视频| 国产97色在线日韩免费| av国产精品久久久久影院| 亚洲欧美精品综合一区二区三区| 久久国产精品男人的天堂亚洲| 美女高潮到喷水免费观看| bbb黄色大片| 一本久久精品| 视频区欧美日本亚洲| 日本猛色少妇xxxxx猛交久久| 亚洲熟女毛片儿| 国产精品二区激情视频| 亚洲欧美精品自产自拍| 免费人妻精品一区二区三区视频| 最新在线观看一区二区三区 | 婷婷丁香在线五月| 亚洲av片天天在线观看| 亚洲久久久国产精品| 久久久久久亚洲精品国产蜜桃av| 久久人人爽人人片av| 国产成人系列免费观看| 国产午夜精品一二区理论片| 熟女av电影| 久久天堂一区二区三区四区| 久久精品国产亚洲av高清一级| 免费不卡黄色视频| 啦啦啦在线免费观看视频4| 侵犯人妻中文字幕一二三四区| 91精品三级在线观看| 99久久精品国产亚洲精品| 又黄又粗又硬又大视频| 少妇精品久久久久久久| 热99久久久久精品小说推荐| 女性生殖器流出的白浆| 电影成人av| 久久久久国产一级毛片高清牌| 国产免费一区二区三区四区乱码| 成人国产一区最新在线观看 | 9191精品国产免费久久| 亚洲第一青青草原| 夫妻午夜视频| 男人添女人高潮全过程视频| 国产成人影院久久av| 国产99久久九九免费精品| 日韩中文字幕欧美一区二区 | 亚洲五月婷婷丁香| 国产成人一区二区三区免费视频网站 | videosex国产| 这个男人来自地球电影免费观看| 欧美日韩亚洲综合一区二区三区_| 精品一区二区三区av网在线观看 | videosex国产| bbb黄色大片| 国产av国产精品国产| 欧美精品亚洲一区二区| 国产淫语在线视频| 久久99一区二区三区| 美女扒开内裤让男人捅视频| a 毛片基地| 国产精品久久久人人做人人爽| 久久天躁狠狠躁夜夜2o2o | 亚洲国产成人一精品久久久| 亚洲中文av在线| 亚洲成国产人片在线观看| 日韩一区二区三区影片| 久久久久精品人妻al黑| 亚洲人成电影观看| 大片电影免费在线观看免费| 精品第一国产精品| 久久国产精品影院| 国产淫语在线视频| 亚洲五月婷婷丁香| 天堂8中文在线网| 999精品在线视频| 亚洲九九香蕉| 一本久久精品| 夫妻性生交免费视频一级片| 国产免费视频播放在线视频| 久久国产精品男人的天堂亚洲| 如日韩欧美国产精品一区二区三区| 国产精品免费大片| 大片免费播放器 马上看| 十八禁高潮呻吟视频| 女性生殖器流出的白浆| 又大又黄又爽视频免费| 精品卡一卡二卡四卡免费| 午夜福利视频在线观看免费| 国产熟女欧美一区二区| 成年女人毛片免费观看观看9 | 国产精品av久久久久免费| 91精品伊人久久大香线蕉| 成年人免费黄色播放视频| 我的亚洲天堂| tube8黄色片| 亚洲熟女精品中文字幕| 丁香六月欧美| 亚洲国产欧美日韩在线播放| 男女床上黄色一级片免费看| 亚洲国产日韩一区二区| 男女免费视频国产| 成年女人毛片免费观看观看9 | 天天躁夜夜躁狠狠躁躁| 久久久亚洲精品成人影院| 精品卡一卡二卡四卡免费| 天天影视国产精品| 亚洲国产精品一区二区三区在线| 国产97色在线日韩免费| 啦啦啦视频在线资源免费观看| 午夜视频精品福利| 啦啦啦 在线观看视频| 在线观看免费日韩欧美大片| 日日夜夜操网爽| 宅男免费午夜| 精品久久蜜臀av无| 啦啦啦啦在线视频资源| 免费观看av网站的网址| 国产成人欧美在线观看 | 中文字幕高清在线视频| 97精品久久久久久久久久精品| kizo精华| 久久性视频一级片| 99久久99久久久精品蜜桃| av视频免费观看在线观看| 99热国产这里只有精品6| 亚洲专区国产一区二区| 久久人妻福利社区极品人妻图片 | 一区二区三区乱码不卡18| 蜜桃在线观看..| 久久久精品国产亚洲av高清涩受| 久热爱精品视频在线9| 99精国产麻豆久久婷婷| 久久久久精品国产欧美久久久 | 黄色怎么调成土黄色| 欧美日韩亚洲高清精品| 欧美黑人欧美精品刺激| 女人精品久久久久毛片| 大片免费播放器 马上看| 国产成人系列免费观看| 麻豆av在线久日| 亚洲欧美日韩高清在线视频 | 少妇人妻久久综合中文| 水蜜桃什么品种好| 成在线人永久免费视频| 两个人看的免费小视频| 巨乳人妻的诱惑在线观看| 少妇人妻 视频| 美女国产高潮福利片在线看| 国产真人三级小视频在线观看| 亚洲激情五月婷婷啪啪| 久久国产精品人妻蜜桃| 丁香六月天网| 汤姆久久久久久久影院中文字幕| 欧美国产精品一级二级三级| 五月天丁香电影| 男女床上黄色一级片免费看| 婷婷丁香在线五月| 精品一区二区三卡| 久久ye,这里只有精品| 亚洲av美国av| 国产成人a∨麻豆精品| 国产精品一区二区在线观看99| 9色porny在线观看| 又紧又爽又黄一区二区| 女人高潮潮喷娇喘18禁视频| 性色av乱码一区二区三区2| 欧美少妇被猛烈插入视频| 飞空精品影院首页| 日本a在线网址| 男女国产视频网站| 极品少妇高潮喷水抽搐| 欧美人与性动交α欧美精品济南到| 欧美 亚洲 国产 日韩一| av有码第一页| 天堂8中文在线网| 久久精品aⅴ一区二区三区四区| 操美女的视频在线观看| 亚洲精品中文字幕在线视频| 久久免费观看电影| 欧美人与性动交α欧美软件| 国产一卡二卡三卡精品| 亚洲国产av新网站| 又黄又粗又硬又大视频| kizo精华| 亚洲国产av影院在线观看| 一二三四在线观看免费中文在| 欧美激情极品国产一区二区三区| 大码成人一级视频| 超碰97精品在线观看| 国产成人精品无人区| 少妇精品久久久久久久| 国产黄色视频一区二区在线观看| 久久久国产精品麻豆| 制服人妻中文乱码| 丝袜脚勾引网站| 亚洲国产看品久久| 久久久久久久久久久久大奶| av国产久精品久网站免费入址| 久久鲁丝午夜福利片| 一本色道久久久久久精品综合| 欧美日韩黄片免| 大话2 男鬼变身卡| 亚洲精品成人av观看孕妇| 欧美精品一区二区免费开放| 老熟女久久久| 一级片'在线观看视频| 国产又色又爽无遮挡免| 国产1区2区3区精品| 少妇裸体淫交视频免费看高清 | 精品亚洲成国产av| 久久精品国产亚洲av涩爱| 又大又爽又粗| 国产在视频线精品| 人人妻人人爽人人添夜夜欢视频| 国产精品人妻久久久影院| 热re99久久精品国产66热6| 男女边摸边吃奶| 曰老女人黄片| 90打野战视频偷拍视频| 午夜久久久在线观看| 女人被躁到高潮嗷嗷叫费观| 看十八女毛片水多多多| 国产国语露脸激情在线看| 久久久久国产精品人妻一区二区| 国产在线免费精品| 国产男女内射视频| videos熟女内射| 菩萨蛮人人尽说江南好唐韦庄| 亚洲成人免费av在线播放| 女人精品久久久久毛片| 欧美+亚洲+日韩+国产| 嫁个100分男人电影在线观看 | 亚洲欧美色中文字幕在线| 大话2 男鬼变身卡| 欧美日韩亚洲国产一区二区在线观看 | 免费看不卡的av| 婷婷色综合大香蕉| 一个人免费看片子| 青春草亚洲视频在线观看| 久久毛片免费看一区二区三区| 亚洲成国产人片在线观看| 欧美+亚洲+日韩+国产| 久久九九热精品免费| 国产精品熟女久久久久浪| 国产淫语在线视频| 亚洲一码二码三码区别大吗| 在线观看免费高清a一片| 久久亚洲国产成人精品v| 激情五月婷婷亚洲| 无遮挡黄片免费观看| 国产一区二区三区av在线| 久久人人爽人人片av| 男女边摸边吃奶| 国产欧美日韩一区二区三区在线| 啦啦啦在线观看免费高清www| kizo精华| 我的亚洲天堂| 欧美人与善性xxx| 伦理电影免费视频| tube8黄色片| 中文字幕人妻熟女乱码| 我要看黄色一级片免费的| 国产成人一区二区在线| 爱豆传媒免费全集在线观看| 97精品久久久久久久久久精品| 永久免费av网站大全| 亚洲精品国产av成人精品| 一本一本久久a久久精品综合妖精| 成人18禁高潮啪啪吃奶动态图| 国产成人一区二区在线| 日韩免费高清中文字幕av| 免费高清在线观看视频在线观看| 新久久久久国产一级毛片| 久久久国产一区二区| 亚洲成av片中文字幕在线观看| 日韩免费高清中文字幕av| 一区二区三区精品91| 91字幕亚洲| 国产一区二区在线观看av| 另类精品久久| 只有这里有精品99| 免费日韩欧美在线观看| 久久精品aⅴ一区二区三区四区| 一区在线观看完整版| 亚洲国产欧美网| 国产淫语在线视频| 永久免费av网站大全| 国产1区2区3区精品| 精品一区二区三卡| 国产男女内射视频| 国产一区有黄有色的免费视频| 两个人看的免费小视频| 成年美女黄网站色视频大全免费| 日韩中文字幕欧美一区二区 | 免费日韩欧美在线观看| 久久鲁丝午夜福利片| 日韩av免费高清视频| 妹子高潮喷水视频| 色网站视频免费| 在线观看国产h片| 少妇精品久久久久久久| 最黄视频免费看| 欧美日韩黄片免| 欧美精品一区二区免费开放| 中文字幕亚洲精品专区| 777米奇影视久久| 黄色视频不卡| 热99国产精品久久久久久7| 男的添女的下面高潮视频| 一级黄色大片毛片| 日韩一本色道免费dvd| 手机成人av网站| 久久毛片免费看一区二区三区| 国产精品亚洲av一区麻豆| 999久久久国产精品视频| 制服诱惑二区| 亚洲精品国产区一区二| 麻豆国产av国片精品| 黄色 视频免费看| 久久性视频一级片| 国产欧美日韩一区二区三区在线| 大码成人一级视频| 久久毛片免费看一区二区三区| 亚洲情色 制服丝袜| 欧美日本中文国产一区发布| 亚洲欧美日韩高清在线视频 | 亚洲综合色网址| 极品少妇高潮喷水抽搐| 日本wwww免费看| 亚洲精品中文字幕在线视频| 日韩一本色道免费dvd| 免费日韩欧美在线观看| cao死你这个sao货| 亚洲欧美中文字幕日韩二区| 一级片免费观看大全| 一级黄片播放器| 嫁个100分男人电影在线观看 | 亚洲欧洲日产国产| 亚洲av日韩精品久久久久久密 | 高清黄色对白视频在线免费看| 国产精品久久久久久精品古装| 国产一区二区在线观看av| 亚洲第一av免费看| www.精华液| 另类亚洲欧美激情| 欧美久久黑人一区二区| 两个人免费观看高清视频| 日韩伦理黄色片| 婷婷色av中文字幕| 黄色片一级片一级黄色片| 成年av动漫网址| 99热网站在线观看| 亚洲国产看品久久| 黑人欧美特级aaaaaa片| 大香蕉久久网| 亚洲自偷自拍图片 自拍| 中文欧美无线码| 在线观看免费视频网站a站| 丰满少妇做爰视频| 久9热在线精品视频| 赤兔流量卡办理| 欧美 亚洲 国产 日韩一| bbb黄色大片| 人人妻人人爽人人添夜夜欢视频| 手机成人av网站| 女人精品久久久久毛片| 九草在线视频观看| 亚洲三区欧美一区|