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

    基于改進(jìn)的遷移率模型的生物地理學(xué)優(yōu)化算法

    2019-10-31 09:21王雅萍張正軍顏?zhàn)雍?/span>金亞洲
    計(jì)算機(jī)應(yīng)用 2019年9期

    王雅萍 張正軍 顏?zhàn)雍? 金亞洲

    摘 要:生物地理學(xué)優(yōu)化(BBO)算法通過遷移和變異不斷更新棲息地,以尋找最優(yōu)解,其中遷移率模型的優(yōu)劣會(huì)直接影響算法的優(yōu)化性能。針對原始BBO算法采用線性遷移率模型適應(yīng)性不足的問題,基于Logistic函數(shù)、三次多項(xiàng)式函數(shù)以及雙曲正切函數(shù)提出了三種新的非線性遷移率模型,并應(yīng)用于原始BBO算法中。對17個(gè)典型的基準(zhǔn)函數(shù)進(jìn)行優(yōu)化性能測試,結(jié)果表明,基于雙曲正切函數(shù)的遷移率模型所得解更接近函數(shù)的全局最小值,總體表現(xiàn)優(yōu)于原始線性遷移率模型的BBO算法以及相關(guān)改進(jìn)算法中表現(xiàn)優(yōu)異的余弦遷移率模型。穩(wěn)定性測試結(jié)果表明,在不同的變異率下,基于雙曲正切函數(shù)的遷移率模型在多數(shù)測試函數(shù)上表現(xiàn)優(yōu)于原始線性遷移率模型。在滿足解多樣性的基礎(chǔ)上,該模型能夠較好地適應(yīng)非線性遷移問題,提高尋優(yōu)能力。

    關(guān)鍵詞:生物地理學(xué)優(yōu)化;遷移率模型;邏輯回歸函數(shù);三次多項(xiàng)式;雙曲正切函數(shù)

    中圖分類號(hào):TP301.6

    文獻(xiàn)標(biāo)志碼:A

    Biogeography-based optimization algorithms based on improved migration rate models

    WANG Yaping, ZHANG Zhengjun*, YAN Zihan, JIN Yazhou

    School of Science, Nanjing University of Science and Technology, Nanjing Jiangsu 210094, China

    Abstract:

    Biogeography-Based Optimization (BBO) algorithm updates habitats through migration and mutation continuously to find the optimal solution, and the migration model affects the performance of the algorithm significantly. In view of the problem of insufficient adaptability of the linear migration model used in the original BBO algorithm, three nonlinear migration models were proposed. These models are based on Logistic function, cubic polynomial function and hyperbolic tangent function respectively. Optimization experiments were carried out on 17 typical benchmark functions, and results show that the migration model based on hyperbolic tangent function performs better than the linear migration model of original BBO algorithm and cosine migration model with good performance of improved algorithm. Stability test shows that the migration model based on hyperbolic tangent function performs better than the original linear migration model with different mutation rates on most test functions. The model satisfies the diversity of the solutions, and better adapts to the nonlinear migration problem with improved search ability.

    Key words:

    Biogeography-Based Optimization (BBO); migration rate model; logistic regression function; cubic polynomial; hyperbolic tangent function

    0 引言

    受生物地理學(xué)[1-5]啟發(fā),Simon[6]于2008年提出了生物地理學(xué)優(yōu)化(Biogeography-Based Optimization, BBO)算法。BBO算法中的主要操作為遷移和變異。根據(jù)物種中每個(gè)個(gè)體的適應(yīng)度函數(shù),通過遷移和變異尋找全局最小值點(diǎn)。Simon在最初的BBO算法中采用線性遷移率模型,并實(shí)證證明與其他基于物種的優(yōu)化算法相比,BBO在多數(shù)特定測試函數(shù)上表現(xiàn)更好,具有一定的實(shí)際應(yīng)用價(jià)值。

    馬海平等[7]根據(jù)生物地理學(xué)物種分布情況,提出了指數(shù)遷移率模型、二次遷移率模型和余弦遷移率模型,在8個(gè)基準(zhǔn)函數(shù)上進(jìn)行測試分析;Ma[8]還提出了常數(shù)與線性混合遷移率模型和梯形遷移率模型[8]。Simon[9]提出了簡化版的BBO算法,并且利用概率論分析其優(yōu)化性能。吳斌等[10]提出了三種遷移操作模式,實(shí)驗(yàn)結(jié)果表明遷入主導(dǎo)的部分遷移式算子優(yōu)化性能最佳。Feng等[11]在更新棲息地的過程中引入了隨機(jī)棲息地,進(jìn)而豐富物種多樣性。Farswan等[12]在更新棲息地的過程中引入了遷入棲息地、隨機(jī)棲息地,增加了解的多樣性;他們還引入最佳棲息地,進(jìn)一步修正更新規(guī)則[13]。Guo等[14]提出均勻混合遷移算子、啟發(fā)式遷移算子以及拓展式遷移算子,提高了基本BBO算法的探索能力。唐繼勇等[15]提出了一種基于動(dòng)態(tài)選擇遷出地與混合自適應(yīng)遷入的優(yōu)化策略,提高了全局搜索能力。

    本文介紹了BBO算法,并針對BBO算法中線性遷移率模型適應(yīng)性不足的問題,分別基于Logistic函數(shù)、三次多項(xiàng)式以及雙曲正切函數(shù)提出了三種非線性遷移率模型應(yīng)用于BBO算法的遷移操作。實(shí)證表明,基于雙曲正切函數(shù)遷移率模型的BBO算法在多數(shù)測試函數(shù)下表現(xiàn)最佳,其解更接近函數(shù)的全局最小值。

    1 生物地理學(xué)優(yōu)化算法

    1.1 生物地理學(xué)

    生物地理學(xué)由自然學(xué)家Alfred Wallace和Charles Darwin于19世紀(jì)提出,其數(shù)學(xué)模型描述了物種在棲息地之間的遷移、新物種的產(chǎn)生以及物種的滅絕。棲息地適宜度指數(shù)(Habitat Suitability Index, HSI)用以衡量棲息地適合物種生存的程度,該值越大則越適宜物種生存。HSI由許多因素構(gòu)成,包括氣溫、降水、植被的多樣性以及土地面積等。這些影響HIS的因素稱為適宜度指數(shù)變量(Suitability Index Variable, SIV)。由于棲息地的資源有限,通常適宜度高的棲息地趨向于容納較多物種,且本地物種有較大的可能性遷出到其他棲息地,同時(shí)外來物種的遷入率較低;而適宜度低的棲息地趨向于包含較少的物種,且外來物種有較大的可能性遷入本地,同時(shí)遷出的可能性較小。此外,自然災(zāi)害等大變動(dòng)事件也可能會(huì)改變棲息地的狀態(tài),進(jìn)而影響HSI。

    為了形式化地描述BBO模型,假設(shè)某棲息地容納物種數(shù)量為S的概率為Ps。則由Ps在t時(shí)刻的取值,可得其在(t+Δt)時(shí)刻的取值:

    Ps(t+Δt)=Ps(t)(1-λsΔt-μsΔt)+Ps-1(t)λs-1Δt+Ps+1(t)λs+1Δt(1)

    其中:λs和μs分別為當(dāng)前棲息地物種數(shù)量為S時(shí)的遷入率和遷出率。式(1)是考慮以下三種情況所得的:

    1)當(dāng)前棲息地在t時(shí)刻的物種數(shù)量為S,并且在接下去的Δt時(shí)間內(nèi)不發(fā)生遷移。

    2)當(dāng)前棲息地在t時(shí)刻的物種數(shù)量為S-1,并且在接下去的Δt時(shí)間內(nèi)有一個(gè)物種遷入。

    3)當(dāng)前棲息地在t時(shí)刻的物種數(shù)量為S+1,并且在接下去的Δt時(shí)間內(nèi)有一個(gè)物種遷出。

    當(dāng)Δt足夠小時(shí),期間發(fā)生超過一次遷移的概率幾乎為0。令P=[P0,P1,…,Pn]T,n=Smax ,由此可得式(1)的極限形式:

    =AP(2)

    1.2 生物地理學(xué)優(yōu)化算法

    基于以上生物地理學(xué)理論,通過遷移和變異不斷更新棲息地的適應(yīng)度,構(gòu)成了生物地理學(xué)優(yōu)化算法。遷移步驟如算法1所示。

    算法1 BBO算法中的遷移步驟。

    程序前

    以正比于λi的概率選擇遷入棲息地Hi

    If Hi被選擇

    For每個(gè)棲息地

    以正比于μi的概率選擇遷出棲息地Hj

    If Hj被選擇

    從Hj中隨機(jī)選擇遷出SIV σ

    將Hi中隨機(jī)選取的SIV替換為σ

    End

    End

    End

    程序后

    變異步驟如算法2所示。

    算法2 BBO算法中的變異步驟。

    程序前

    對于棲息地Hi

    For 每個(gè)SIV

    根據(jù)λi和μi計(jì)算變異概率Pi

    以正比于Pi的概率選擇變異SIV為Hi(j)

    If Hi(j)被選擇

    將Hi(j)替換為隨機(jī)產(chǎn)生的SIV

    End

    End

    程序后

    BBO算法流程如算法3所示。

    算法3 BBO算法。

    程序前

    初始化物種數(shù)量、遷移率、變異率以及精英數(shù)。

    計(jì)算各棲息地的HSI

    While 停止條件未滿足 do

    將各棲息地按照HIS從高到低進(jìn)行排序

    利用算法1進(jìn)行遷移操作

    更新各棲息地的HIS

    利用算法2進(jìn)行變異操作

    保持當(dāng)前迭代所有棲息地中的精英

    End while

    程序后

    2 遷移率模型

    Simon采用線性模型來模擬自然界中物種遷移的規(guī)律,即認(rèn)為隨著物種數(shù)量的增加,遷入率呈線性遞減趨勢,同時(shí)遷出率呈線性遞增趨勢。在相關(guān)改進(jìn)算法中,Ma等[7]提出的余弦遷移率模型在其測試的多數(shù)函數(shù)下表現(xiàn)最佳,認(rèn)為接近自然規(guī)律的遷移率模型優(yōu)于簡單的線性模型。

    為了更好地適應(yīng)非線性遷移問題,本文提出三種不同變化趨勢的非線性遷移率模型,并與原始線性模型以及相關(guān)改進(jìn)算法中表現(xiàn)優(yōu)異的余弦遷移率模型[7]進(jìn)行對比分析。如圖1所示,在最大遷入率I和最大遷出率E下,遷入率λk和遷出率μk是關(guān)于物種數(shù)量k的函數(shù),k0為遷入率等于遷出率的平衡點(diǎn)。

    模型1 原始線性模型。

    線性模型的遷入率和遷出率如下:

    λk=I(1-k/n)

    μk=E(k/n)(9)

    如圖1(a)所示,在線性遷移率模型下,隨著物種數(shù)量的增加,遷入率從I開始線性下降至0,同時(shí)遷出率從0開始線性上升至E。

    模型2 余弦模型。

    余弦模型[7]的遷入率和遷出率如下。

    λk=I2coskπn+1

    μk=E2-coskπn+1(10)

    如圖1(b)所示,在余弦遷移率模型下,當(dāng)物種數(shù)量較多或較少時(shí),遷移率變化比較平穩(wěn);當(dāng)具有一定物種數(shù)量時(shí),遷移率相對變化較快。

    模型3 Logistic模型。

    Logistic模型的遷入率和遷出率如下:

    λk=I1-11+en2-k

    μk=E1+en2-k(11)

    如圖1(c)所示,在Logistic遷移率模型下,當(dāng)物種數(shù)量較多或較少時(shí),遷移率變化非常平穩(wěn);當(dāng)具有一定物種數(shù)量時(shí),遷移率變化很快。

    模型4 三次多項(xiàng)式模型。

    三次多項(xiàng)式模型的遷入率和遷出率如下:

    λk=I2-(2k-n)3n3+1

    μk=E2(2k-n)3n3+1(12)

    如圖1(d)所示,在三次多項(xiàng)式遷移率模型下,在物種數(shù)量較多或較少時(shí),遷移率變化較快;而具有一定物種數(shù)量時(shí),遷移率變化趨于平緩。

    模型5 雙曲正切變型模型。

    雙曲正切變型模型的遷入率和遷出率如下。

    λk=I2-ak-n2-a-k+n2ak-n2+a-k+n2+1

    μk=E2ak-n2-a-k+n2ak-n2+a-k+n2+1(13)

    如圖1(e)所示,在雙曲正切變型遷移率模型(本文取a=1.1)下,遷移率隨物種數(shù)量變化的趨勢類似于余弦模型,但是幅度比余弦模型更緩和些。

    3 實(shí)證分析

    在python3.6的環(huán)境下,針對17個(gè)典型的基準(zhǔn)函數(shù),分別利用提出的三種基于非線性遷移率模型的BBO算法進(jìn)行測試,并與原始線性遷移率模型以及相關(guān)改進(jìn)算法中表現(xiàn)優(yōu)異的余弦遷移率模型[7]進(jìn)行對比分析?;鶞?zhǔn)函數(shù)信息如下所示。

    1)Sphere:

    f1(x)=∑Di=1xi2; -10≤xi≤10, D=20

    2)Schwefel221:

    f2(x)=maxxi; 1≤i≤D, -10≤xi≤10, D=20

    3)Griewank:

    f3(x)=∑Di=1xi24000-∏Di=1cos(xi/i)+1;

    -600≤xi≤600, D=20

    4)Axis parallel hyper ellipsoid:

    f4(x)=∑Di=1ixi2; -5.12≤xi≤5.12, D=20

    5)Alpine:

    f5(x)=∑Di=1xi sin xi+0.1xi; -10≤xi≤10, D=20

    6)Exponential:

    f6(x)=-exp(-0.5∑Di=1xi2); -1≤xi≤1, D=20

    7)Cosine mixture:

    f7(x)=∑Di=1xi2-0.1∑Di=1cos(5πxi); -1≤xi≤1,D=20

    8)Zakharov function:

    f8(x)=∑Di=1xi2+(∑Di=10.5ixi)2+(∑Di=10.5ixi)4;

    -5≤xi≤10, D=20

    9)De jongs:

    f9(x)=∑Di=1ixi4; -5.12≤xi≤5.12, D=20

    10)Michalewicz:

    f10(x)=-∑Di=1sin(xi)sin2m(ix2i/π); 0≤xi≤π, D=10

    11)Easom:

    f11(x)=-cos(x1)cos(x2)e-(x1-π)2-(x2-π)2,

    -100≤xi≤100, D=2

    12)Ackley:

    f12(x)=-20 exp(-0.21D∑Di=1x2i)-

    exp(1D∑Di=1cos(2πxi))+20+e; -30≤xi≤30, D=20

    13)Schwefel222:

    f13(x)=∑Di=1xi+∏Di=1xi; -10≤xi≤10, D=20

    14)Step function

    f14(x)=∑Di=1(xi+0.5」)2; -100≤xi≤100, D=20

    15)Schwefel:

    f15(x)=418.9829D-∑Di=1xi sin(xi12);-500≤xi≤500, D=20

    16)Rastrigin:

    f16(x)=10D+∑Di=1[xi2-10 cos(2πxi)];-5.12≤xi≤5.12, D=20

    17)Salomon prol 3:

    f17(x)=1-cos(2π∑Di=1xi2)+0.1∑Di=1xi2;-100≤xi≤100, D=20

    所有測試函數(shù)有且僅有一個(gè)全局最小值點(diǎn),其中f3、 f5、 f7、 f12、 f15和f16有許多局部極小值點(diǎn),優(yōu)化算法在尋找最優(yōu)解時(shí)容易陷入局部極小值點(diǎn)而無法達(dá)到全局最小值點(diǎn)。

    f10的全局最小值為-9.66015, f6和f11的全局最小值為-1, f7的全局最小值為-2,其余函數(shù)全局最小值均為0。

    3.1 遷移率模型優(yōu)化性能分析

    為了平行比較基于不同遷移率模型的BBO算法的性能,有關(guān)參數(shù)設(shè)置如下:

    最大遷移率I=E=1,物種(棲息地)數(shù)量PopSize=30,迭代次數(shù)Iteration=50,變異率m=0.01。

    為了降低算法中隨機(jī)因素的影響,各模型獨(dú)立運(yùn)行50次,觀察平均值和最小值。5種模型測試結(jié)果的平均值如表1所示。

    由表1可知,從平均值的角度來看,不同遷移率模型優(yōu)化性能差別很大。總體來說,模型5在大多測試函數(shù)上表現(xiàn)最佳,模型2其次,而模型4表現(xiàn)差強(qiáng)人意,僅在3個(gè)測試函數(shù)上優(yōu)于模型1。

    17個(gè)函數(shù)在不同遷移率模型下求解的平均最小值誤差百分比變化趨勢如圖2所示。

    由圖2可知,就平均最小值誤差百分比而言,對于多數(shù)測試函數(shù)來說,模型4的誤差百分比最大,而模型2和模型5長期處于較低的誤差水平。在不考慮模型4的情況下,函數(shù)f5、 f6、 f8、 f9和f17受遷移率模型影響較大,其余函數(shù)受遷移率模型影響相對較小。

    5種模型測試結(jié)果的最佳最小值如表2所示。

    由表2可知,從最小值的角度來看,模型5的表現(xiàn)仍然優(yōu)異,并且模型4的優(yōu)化性能有所提升,而模型2的尋優(yōu)能力有所降低。

    17個(gè)函數(shù)在不同遷移率模型下求解的最佳最小值誤差百分比變化趨勢如圖3所示。

    由圖3可知,就最佳最小值誤差百分比而言,BBO算法尋找函數(shù)最小值時(shí)受遷移率模型影響很大。模型5的總體表現(xiàn)最穩(wěn)定,在14個(gè)測試函數(shù)上誤差小于模型1,并且長期處于低誤差水平。模型2的性能有所下降,僅在9個(gè)測試函數(shù)上誤差小于模型1。

    3.2 穩(wěn)定性分析

    改變物種最大變異率可以調(diào)節(jié)解的多樣性。為了研究不同

    遷移率模型的穩(wěn)定性,對比分析原始線性模型(模型1)和本文所

    提最佳模型(模型5)在不同遷移率下的表現(xiàn),如表3所示。

    由表3可知,不同變異率下所得平均值和最小值變化不一, f7、 f8、 f9受變異率的影響較大,其余函數(shù)受影響較小。總體上變異率為0.01時(shí)模型性能最佳,并且在多數(shù)測試函數(shù)上模型5的表現(xiàn)優(yōu)于模型1。

    4 結(jié)語

    生物地理學(xué)優(yōu)化算法性能受其遷移率模型影響顯著。本文提出三種新的非線性遷移率模型,并應(yīng)用于17個(gè)典型的基準(zhǔn)函數(shù)求最小值問題。結(jié)果顯示,在多數(shù)測試函數(shù)下,基于雙曲正切變體遷移率模型的BBO算法尋優(yōu)能力最佳,其解更接近函數(shù)的全局最小值。變異率作為影響B(tài)BO算法優(yōu)化性能的另一重要因素,有待進(jìn)一步研究。

    參考文獻(xiàn)

    [1]WALLACE A R. The geographical distribution of animals [J]. Nature, 1958, 182(4629): 140-141.

    [2]DARWIN C. The Origin of Species [M]. New York: Bantam Classics, 1995: 218-230.

    [3]HANSKI I, GILPIN M E. Metapopulation Biology [M]. New York: Academic Press, 1997: 63-67.

    [4]XING B, GAO W-J. Biogeography-based optimization algorithm [M]// Innovative Computational Intelligence: A Rough Guide to 134 Clever Algorithms, LSRL 62. Berlin: Springer, 2014: 81-91.

    [5]GUO W, WANG L, WU Q. Numerical comparisons of migration models for multi-objective-biogeography-based optimization [J]. Information Sciences, 2016, 328(C): 302-320.

    [6]SIMON D. Biogeography-based optimization [J]. IEEE Transactions on Evolutionary Computation, 2008, 12(6): 702-713.

    [7]馬海平,李雪,林升東.生物地理學(xué)優(yōu)化算法的遷移率模型分析[J].東南大學(xué)學(xué)報(bào)(自然科學(xué)版),2009,39(S1):16-21.(MA H P, LI X, LIN S D. Analysis of migration rate models for biogeography-based optimization [J]. Journal of Southeast University (Natural Science Edition), 2009, 39(S1): 16-21.)

    [8]MA H. An analysis of the equilibrium of migration models for biogeography-based optimization [J]. Information Sciences—Informatics and Computer Science, Intelligent Systems, Applications: An International Journal, 2010, 180(18): 3444-3464.

    [9]SIMON D. A probabilistic analysis of a simplified biogeography-based optimization algorithm [J]. Evolutionary Computation, 2011, 19(2): 167-188.

    [10]吳斌,林錦國,崔志勇.生物地理學(xué)優(yōu)化算法中遷移算子的比較[J].計(jì)算機(jī)工程與應(yīng)用,2012,48(25):61-64.(WU B, LIN J G, CUI Z Y. Compative study of migration operators in biogeography optimization [J]. Computer Engineering and Applications, 2012, 48(25): 61-64.)

    [11]FENG Q, LIU S, ZHANG J, et al. Biogeography-based optimization with improved migration operator and self-adaptive clear duplicate operator [J]. Applied Intelligence, 2014, 41(2): 563-581.

    [12]FSRSWAN P, BANSAL J C. Migration in biogeography-based optimization [C]// Proceedings of the 4th International Conference on Soft Computing for Problem Solving, AISC 336. Berlin: Springer, 2014: 389-401.

    [13]FARSWAN P, BANSAL J C, DEEP K. A modified biogeography based optimization [M]// Harmony Search Algorithm, AISC 382. Berlin: Springer, 2016: 227-238.

    [14]GUO W, WANG L, SI C, et al. Novel migration operators of biogeography-based optimization and Markov analysis [J]. Soft Computing, 2016, 21(22): 6605-6632.

    [15]唐繼勇,仲元昌,曾廣樸.基于遷出地動(dòng)態(tài)選擇與自適應(yīng)遷入策略的BBO算法[J].計(jì)算機(jī)科學(xué),2016,43(10):282-286.(TANG J Y, ZHONG Y C, ZENG G P. Biogeography-based optimization with adaptive immigration and dynamic selection emigration strategy [J]. Computer Science, 2016, 43(10): 282-286.)

    This work is partially supported by the National Natural Science Foundation of China (61773014).

    WANG Yaping, born in 1995, M. S. candidate. Her research interests include data mining.

    ZHANG Zhengjun, born in 1965, Ph. D., associate professor. His research interests include data mining, graphics and images.

    YAN Zihan, born in 1995, M. S. candidate. Her research interests include data mining.

    JIN Yazhou, born in 1993, M. S. candidate. His research interests include data mining.

    真人做人爱边吃奶动态| 国产亚洲精品久久久久久毛片| 91大片在线观看| 亚洲av电影不卡..在线观看| 美女 人体艺术 gogo| 成人亚洲精品一区在线观看| 99热6这里只有精品| 精品国产乱码久久久久久男人| 色综合亚洲欧美另类图片| 日韩 欧美 亚洲 中文字幕| 婷婷精品国产亚洲av| 国产免费男女视频| 国产熟女xx| 美国免费a级毛片| 欧美日韩瑟瑟在线播放| 日韩欧美 国产精品| 一区二区日韩欧美中文字幕| 亚洲 国产 在线| 国产精品香港三级国产av潘金莲| 欧美最黄视频在线播放免费| 亚洲 欧美一区二区三区| 不卡一级毛片| 久久国产亚洲av麻豆专区| 精品久久蜜臀av无| 欧美又色又爽又黄视频| 黑人巨大精品欧美一区二区mp4| 白带黄色成豆腐渣| 老司机深夜福利视频在线观看| 可以在线观看毛片的网站| 人人妻人人看人人澡| 精品国产乱子伦一区二区三区| 一卡2卡三卡四卡精品乱码亚洲| 99久久国产精品久久久| 一级作爱视频免费观看| 亚洲av电影在线进入| 看片在线看免费视频| 正在播放国产对白刺激| 欧美一级a爱片免费观看看 | 国产亚洲精品久久久久5区| 90打野战视频偷拍视频| 草草在线视频免费看| 久久久久亚洲av毛片大全| 亚洲 国产 在线| 久久久久亚洲av毛片大全| 亚洲真实伦在线观看| 日本黄色视频三级网站网址| 亚洲国产精品成人综合色| 亚洲精品粉嫩美女一区| 黄色视频,在线免费观看| 亚洲av电影在线进入| 亚洲一码二码三码区别大吗| 久久天躁狠狠躁夜夜2o2o| 丝袜人妻中文字幕| av在线播放免费不卡| av片东京热男人的天堂| 国产亚洲精品久久久久5区| 性色av乱码一区二区三区2| 国内毛片毛片毛片毛片毛片| 十八禁人妻一区二区| 两性夫妻黄色片| tocl精华| 国产精华一区二区三区| 色哟哟哟哟哟哟| 亚洲欧美日韩无卡精品| 两个人免费观看高清视频| 国产黄a三级三级三级人| 一个人观看的视频www高清免费观看 | 嫩草影视91久久| 午夜福利18| 欧美+亚洲+日韩+国产| 香蕉国产在线看| 亚洲黑人精品在线| 亚洲精品久久国产高清桃花| 中文字幕人成人乱码亚洲影| 精品国产国语对白av| 国产激情偷乱视频一区二区| 中文字幕精品亚洲无线码一区 | 国产精品一区二区三区四区久久 | 91国产中文字幕| cao死你这个sao货| 免费高清在线观看日韩| 亚洲七黄色美女视频| 亚洲av熟女| 成在线人永久免费视频| 叶爱在线成人免费视频播放| 亚洲欧美一区二区三区黑人| 高清毛片免费观看视频网站| 国产成人一区二区三区免费视频网站| 久久久久国产一级毛片高清牌| 久久久久久亚洲精品国产蜜桃av| 色av中文字幕| 中文字幕久久专区| 国产成+人综合+亚洲专区| 久久国产精品影院| 日日夜夜操网爽| 国产蜜桃级精品一区二区三区| 欧美日韩黄片免| 亚洲成人免费电影在线观看| 亚洲中文字幕一区二区三区有码在线看 | 18禁裸乳无遮挡免费网站照片 | 一进一出抽搐动态| 国产区一区二久久| 国产蜜桃级精品一区二区三区| 老汉色∧v一级毛片| 午夜福利18| 日本免费a在线| 成人欧美大片| 成人精品一区二区免费| 日日摸夜夜添夜夜添小说| 久久国产精品影院| 国产主播在线观看一区二区| 黄色视频,在线免费观看| 久久人妻福利社区极品人妻图片| 成人亚洲精品一区在线观看| 国产伦人伦偷精品视频| 色精品久久人妻99蜜桃| 亚洲成av人片免费观看| 18禁黄网站禁片午夜丰满| 成人国语在线视频| 成熟少妇高潮喷水视频| 波多野结衣巨乳人妻| 在线观看66精品国产| 最好的美女福利视频网| 久久精品国产清高在天天线| 一个人免费在线观看的高清视频| av超薄肉色丝袜交足视频| 免费在线观看亚洲国产| 国产人伦9x9x在线观看| 亚洲最大成人中文| 人人妻,人人澡人人爽秒播| 欧美av亚洲av综合av国产av| a级毛片a级免费在线| 午夜a级毛片| 黄色女人牲交| 变态另类丝袜制服| 可以在线观看毛片的网站| 91成人精品电影| 一个人观看的视频www高清免费观看 | 此物有八面人人有两片| 黑人欧美特级aaaaaa片| 一级作爱视频免费观看| 亚洲精品中文字幕一二三四区| 亚洲成人久久性| 老司机午夜十八禁免费视频| 日本五十路高清| 亚洲国产看品久久| 欧美三级亚洲精品| 黑人欧美特级aaaaaa片| 50天的宝宝边吃奶边哭怎么回事| 精品欧美一区二区三区在线| 亚洲九九香蕉| 亚洲成人精品中文字幕电影| 啦啦啦观看免费观看视频高清| 国产真实乱freesex| 国产日本99.免费观看| 可以在线观看的亚洲视频| 国产一卡二卡三卡精品| 国产精品综合久久久久久久免费| 精品国内亚洲2022精品成人| 亚洲第一欧美日韩一区二区三区| 国产成人一区二区三区免费视频网站| a级毛片在线看网站| 51午夜福利影视在线观看| 老汉色av国产亚洲站长工具| 久久精品国产清高在天天线| 亚洲国产看品久久| 巨乳人妻的诱惑在线观看| 亚洲精品在线观看二区| 成人国产一区最新在线观看| 搞女人的毛片| 午夜久久久久精精品| 一本一本综合久久| 欧美日韩中文字幕国产精品一区二区三区| 欧美日本亚洲视频在线播放| 免费在线观看视频国产中文字幕亚洲| 99热这里只有精品一区 | 一个人免费在线观看的高清视频| 男女床上黄色一级片免费看| 亚洲久久久国产精品| 国产欧美日韩一区二区三| 午夜福利欧美成人| 精华霜和精华液先用哪个| 神马国产精品三级电影在线观看 | 久久精品成人免费网站| 国产亚洲精品第一综合不卡| 国产精品国产高清国产av| 欧美一级a爱片免费观看看 | 久久中文看片网| 好男人电影高清在线观看| 亚洲一区中文字幕在线| 日日摸夜夜添夜夜添小说| svipshipincom国产片| 少妇的丰满在线观看| 可以在线观看的亚洲视频| 99在线视频只有这里精品首页| 国产亚洲av嫩草精品影院| tocl精华| 在线观看免费午夜福利视频| 亚洲欧美激情综合另类| 欧美大码av| 亚洲av电影不卡..在线观看| 欧美色欧美亚洲另类二区| 首页视频小说图片口味搜索| 在线播放国产精品三级| 国产99久久九九免费精品| 亚洲自偷自拍图片 自拍| 国产熟女午夜一区二区三区| 国产麻豆成人av免费视频| 日本一区二区免费在线视频| 国产单亲对白刺激| 91老司机精品| 亚洲第一欧美日韩一区二区三区| 亚洲成国产人片在线观看| 观看免费一级毛片| 久久国产亚洲av麻豆专区| 欧美不卡视频在线免费观看 | 丝袜人妻中文字幕| 国产精品久久久久久人妻精品电影| or卡值多少钱| 亚洲久久久国产精品| 国产一区二区在线av高清观看| 黄色女人牲交| 日韩欧美三级三区| 亚洲精品美女久久久久99蜜臀| 亚洲人成77777在线视频| 国产成人精品无人区| 最近在线观看免费完整版| 啪啪无遮挡十八禁网站| 嫩草影院精品99| 视频区欧美日本亚洲| 两个人视频免费观看高清| 久久久久久人人人人人| 国产成人系列免费观看| 久久99热这里只有精品18| 男女做爰动态图高潮gif福利片| 亚洲精品久久成人aⅴ小说| 天堂动漫精品| 俄罗斯特黄特色一大片| 俺也久久电影网| 免费在线观看亚洲国产| 成人18禁在线播放| 一区二区三区国产精品乱码| 又黄又粗又硬又大视频| 丝袜在线中文字幕| 国产一区在线观看成人免费| 亚洲熟妇熟女久久| 99国产综合亚洲精品| 久久久久久久精品吃奶| 长腿黑丝高跟| 日日干狠狠操夜夜爽| 在线天堂中文资源库| 国产蜜桃级精品一区二区三区| 国产欧美日韩一区二区三| 亚洲成人精品中文字幕电影| av欧美777| 久久久国产精品麻豆| 亚洲狠狠婷婷综合久久图片| 亚洲五月天丁香| 欧美日韩福利视频一区二区| 精品欧美一区二区三区在线| 久久 成人 亚洲| 亚洲 欧美 日韩 在线 免费| 精品国产一区二区三区四区第35| 久久中文字幕人妻熟女| 欧美色视频一区免费| 精品久久久久久成人av| 人人妻人人看人人澡| 国产成+人综合+亚洲专区| 亚洲va日本ⅴa欧美va伊人久久| 黑人操中国人逼视频| 亚洲专区字幕在线| 精品不卡国产一区二区三区| 精品国内亚洲2022精品成人| netflix在线观看网站| 欧美性长视频在线观看| 亚洲精品一卡2卡三卡4卡5卡| 亚洲精品国产区一区二| 男女午夜视频在线观看| 婷婷六月久久综合丁香| 欧美黑人欧美精品刺激| 90打野战视频偷拍视频| 色综合站精品国产| 亚洲av熟女| 麻豆国产av国片精品| 99久久国产精品久久久| 麻豆av在线久日| 久久久国产成人免费| 国产亚洲av嫩草精品影院| 精品午夜福利视频在线观看一区| 黄色毛片三级朝国网站| 国产精品一区二区精品视频观看| 日本五十路高清| 欧美黑人欧美精品刺激| a级毛片在线看网站| 夜夜夜夜夜久久久久| 成人三级做爰电影| 亚洲人成网站高清观看| 欧美乱码精品一区二区三区| 精品国产亚洲在线| 欧美av亚洲av综合av国产av| 色综合站精品国产| 别揉我奶头~嗯~啊~动态视频| 亚洲国产中文字幕在线视频| 久久久久久久久中文| 色老头精品视频在线观看| 国产精品久久久av美女十八| 久久久精品欧美日韩精品| 亚洲一卡2卡3卡4卡5卡精品中文| 亚洲天堂国产精品一区在线| 超碰成人久久| av中文乱码字幕在线| 亚洲精品美女久久av网站| av超薄肉色丝袜交足视频| 久久草成人影院| 9191精品国产免费久久| 亚洲av成人不卡在线观看播放网| 久久婷婷成人综合色麻豆| 亚洲成a人片在线一区二区| www.熟女人妻精品国产| 欧美大码av| 精品国产国语对白av| 天天一区二区日本电影三级| 又紧又爽又黄一区二区| 91成人精品电影| 国产亚洲欧美在线一区二区| 亚洲色图 男人天堂 中文字幕| 国产aⅴ精品一区二区三区波| 曰老女人黄片| 欧美黑人精品巨大| 中国美女看黄片| 国产人伦9x9x在线观看| 国产精品98久久久久久宅男小说| 一级a爱片免费观看的视频| 琪琪午夜伦伦电影理论片6080| 可以在线观看的亚洲视频| 国产aⅴ精品一区二区三区波| 欧美 亚洲 国产 日韩一| 国产精品乱码一区二三区的特点| 黑人巨大精品欧美一区二区mp4| 亚洲成人精品中文字幕电影| 国产精品,欧美在线| or卡值多少钱| 十八禁人妻一区二区| 国产成人av教育| 久久久久久久久免费视频了| 757午夜福利合集在线观看| 日韩免费av在线播放| 亚洲片人在线观看| 嫩草影院精品99| 欧美日韩福利视频一区二区| svipshipincom国产片| 欧美精品啪啪一区二区三区| 久久精品人妻少妇| 国产免费av片在线观看野外av| 免费观看人在逋| 久久精品91无色码中文字幕| 亚洲成国产人片在线观看| 日本在线视频免费播放| 日日夜夜操网爽| 一区二区三区激情视频| 宅男免费午夜| 好男人在线观看高清免费视频 | 亚洲成人免费电影在线观看| 国产一卡二卡三卡精品| 亚洲五月婷婷丁香| 精品国内亚洲2022精品成人| 熟女电影av网| 国产精品香港三级国产av潘金莲| 国内毛片毛片毛片毛片毛片| 欧美乱妇无乱码| 精品久久蜜臀av无| 成人永久免费在线观看视频| 免费看a级黄色片| 久久久久久大精品| 午夜精品在线福利| 黄色丝袜av网址大全| 久久精品91蜜桃| 19禁男女啪啪无遮挡网站| 国产极品粉嫩免费观看在线| 国产激情欧美一区二区| 伊人久久大香线蕉亚洲五| 露出奶头的视频| 欧美丝袜亚洲另类 | 国产黄色小视频在线观看| 嫁个100分男人电影在线观看| 日日摸夜夜添夜夜添小说| 久久热在线av| 麻豆成人av在线观看| 欧美黄色片欧美黄色片| 欧美人与性动交α欧美精品济南到| 99久久综合精品五月天人人| 一级毛片精品| 999精品在线视频| 国产黄a三级三级三级人| 一级作爱视频免费观看| 国产亚洲精品av在线| 一区福利在线观看| 给我免费播放毛片高清在线观看| 中文字幕久久专区| 亚洲精品一区av在线观看| 日韩av在线大香蕉| 高清在线国产一区| 在线观看66精品国产| 亚洲熟妇熟女久久| 免费在线观看亚洲国产| 国产亚洲精品综合一区在线观看 | 性欧美人与动物交配| 国产在线精品亚洲第一网站| 日韩中文字幕欧美一区二区| 宅男免费午夜| 1024手机看黄色片| 黑人欧美特级aaaaaa片| 一边摸一边做爽爽视频免费| 久久精品国产清高在天天线| 好男人在线观看高清免费视频 | 在线观看一区二区三区| 久久精品91无色码中文字幕| 窝窝影院91人妻| 欧美av亚洲av综合av国产av| 国产精品久久视频播放| 国产av又大| 久久国产亚洲av麻豆专区| 国内揄拍国产精品人妻在线 | 亚洲精品av麻豆狂野| 男女做爰动态图高潮gif福利片| 老司机午夜十八禁免费视频| 麻豆成人av在线观看| 免费人成视频x8x8入口观看| 激情在线观看视频在线高清| 精品久久久久久久末码| 精品国产超薄肉色丝袜足j| 岛国在线观看网站| а√天堂www在线а√下载| 搞女人的毛片| 亚洲黑人精品在线| 又紧又爽又黄一区二区| 久久午夜亚洲精品久久| 国产高清视频在线播放一区| 在线看三级毛片| 亚洲自拍偷在线| 国产成人欧美在线观看| 久久精品aⅴ一区二区三区四区| 成人18禁在线播放| 国产精品一区二区三区四区久久 | 国产主播在线观看一区二区| 成人18禁高潮啪啪吃奶动态图| 欧美午夜高清在线| 精品一区二区三区av网在线观看| 国产精品日韩av在线免费观看| 精品卡一卡二卡四卡免费| 香蕉丝袜av| 亚洲一区中文字幕在线| 久久久久九九精品影院| 国产成人欧美| xxxwww97欧美| 午夜a级毛片| 天堂影院成人在线观看| 99在线人妻在线中文字幕| 中文字幕高清在线视频| 51午夜福利影视在线观看| 欧美精品亚洲一区二区| 久久精品91无色码中文字幕| 熟女电影av网| 亚洲成人久久爱视频| 女人高潮潮喷娇喘18禁视频| 欧美另类亚洲清纯唯美| 免费在线观看完整版高清| 国产精品电影一区二区三区| 精品久久久久久久毛片微露脸| 一本大道久久a久久精品| 两个人视频免费观看高清| 一二三四在线观看免费中文在| 久久性视频一级片| 此物有八面人人有两片| 色综合欧美亚洲国产小说| 日韩精品免费视频一区二区三区| aaaaa片日本免费| 亚洲精品一区av在线观看| 欧美成人免费av一区二区三区| 久久99热这里只有精品18| 91老司机精品| 黄色视频,在线免费观看| 又黄又爽又免费观看的视频| 午夜福利视频1000在线观看| 男人舔奶头视频| 免费看日本二区| 男人舔女人下体高潮全视频| 黄片播放在线免费| 日韩视频一区二区在线观看| 白带黄色成豆腐渣| 亚洲五月婷婷丁香| 日韩欧美国产一区二区入口| 亚洲av电影不卡..在线观看| 岛国视频午夜一区免费看| 一区二区三区国产精品乱码| 最近最新中文字幕大全免费视频| 夜夜看夜夜爽夜夜摸| 国产极品粉嫩免费观看在线| 巨乳人妻的诱惑在线观看| 欧美大码av| 琪琪午夜伦伦电影理论片6080| 国产免费男女视频| 久久久久九九精品影院| www.自偷自拍.com| 亚洲男人的天堂狠狠| 男女床上黄色一级片免费看| 久久精品国产综合久久久| cao死你这个sao货| 欧美亚洲日本最大视频资源| 国产97色在线日韩免费| 一区二区日韩欧美中文字幕| 啦啦啦韩国在线观看视频| 久久人妻福利社区极品人妻图片| 在线观看66精品国产| 久久久国产成人免费| 亚洲精品一区av在线观看| 亚洲男人的天堂狠狠| 久久久久久久久久黄片| 日韩大尺度精品在线看网址| 88av欧美| 中文字幕久久专区| 99在线人妻在线中文字幕| 一级毛片女人18水好多| 亚洲熟妇熟女久久| 丁香六月欧美| 91成年电影在线观看| 亚洲精品av麻豆狂野| 19禁男女啪啪无遮挡网站| 国产精品电影一区二区三区| 国产一区二区激情短视频| 免费一级毛片在线播放高清视频| 国产精品日韩av在线免费观看| 免费一级毛片在线播放高清视频| 视频区欧美日本亚洲| 一区二区三区国产精品乱码| 国产一区二区三区在线臀色熟女| 中文资源天堂在线| 性色av乱码一区二区三区2| 不卡av一区二区三区| 亚洲av熟女| 欧美精品亚洲一区二区| 99riav亚洲国产免费| 国产亚洲av高清不卡| www.www免费av| 在线观看www视频免费| 国产av又大| ponron亚洲| 操出白浆在线播放| 午夜免费成人在线视频| 国产单亲对白刺激| 亚洲片人在线观看| 中文在线观看免费www的网站 | 久久久水蜜桃国产精品网| 一本大道久久a久久精品| 桃红色精品国产亚洲av| 亚洲成人精品中文字幕电影| 国产在线精品亚洲第一网站| 国产伦人伦偷精品视频| 动漫黄色视频在线观看| 久久草成人影院| 亚洲熟妇中文字幕五十中出| 亚洲专区字幕在线| 亚洲色图 男人天堂 中文字幕| 搡老岳熟女国产| 看免费av毛片| 99久久无色码亚洲精品果冻| 国产高清视频在线播放一区| 国产精品久久电影中文字幕| 欧美性猛交黑人性爽| 黄色视频不卡| 日本五十路高清| 亚洲国产欧美一区二区综合| 亚洲成国产人片在线观看| 亚洲午夜精品一区,二区,三区| 日韩视频一区二区在线观看| 欧美精品啪啪一区二区三区| 午夜精品在线福利| 日韩欧美一区二区三区在线观看| 国产精品一区二区免费欧美| 国产一区二区三区在线臀色熟女| 欧美国产精品va在线观看不卡| 桃红色精品国产亚洲av| 免费在线观看影片大全网站| 校园春色视频在线观看| 男女床上黄色一级片免费看| 又黄又粗又硬又大视频| 真人一进一出gif抽搐免费| 天天躁狠狠躁夜夜躁狠狠躁| 亚洲欧美日韩高清在线视频| 满18在线观看网站| 色哟哟哟哟哟哟| 欧美另类亚洲清纯唯美| 国产亚洲欧美在线一区二区| 久久久久久久久免费视频了| 日韩 欧美 亚洲 中文字幕| 两人在一起打扑克的视频| 少妇 在线观看| 制服诱惑二区| 亚洲五月色婷婷综合| 国产一区二区激情短视频| 最新美女视频免费是黄的| 精品国产美女av久久久久小说| 成人国语在线视频| 韩国精品一区二区三区| 麻豆国产av国片精品| 精品一区二区三区av网在线观看| 婷婷精品国产亚洲av| 亚洲国产精品成人综合色| 欧美日本视频| 亚洲九九香蕉| 亚洲精品av麻豆狂野| 欧美在线黄色| 性欧美人与动物交配| 午夜福利高清视频| 亚洲国产欧美一区二区综合|