• <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.

    欧美高清成人免费视频www| 又黄又爽又刺激的免费视频.| 简卡轻食公司| 老司机影院成人| 欧美日韩视频高清一区二区三区二| 亚洲av电影在线观看一区二区三区 | 国模一区二区三区四区视频| 日韩成人av中文字幕在线观看| 欧美精品国产亚洲| 日韩大片免费观看网站| 国产精品久久久久久精品电影| 亚洲av成人精品一二三区| 国产精品无大码| 2022亚洲国产成人精品| 中文天堂在线官网| 日韩一本色道免费dvd| 日本免费在线观看一区| 国产成人精品久久久久久| 国产高潮美女av| 午夜福利网站1000一区二区三区| 国产女主播在线喷水免费视频网站 | 黄片无遮挡物在线观看| 久久精品人妻少妇| 神马国产精品三级电影在线观看| 男女边摸边吃奶| 亚洲自拍偷在线| 日本色播在线视频| 街头女战士在线观看网站| 又大又黄又爽视频免费| 欧美bdsm另类| 看非洲黑人一级黄片| 五月伊人婷婷丁香| 成人国产麻豆网| 26uuu在线亚洲综合色| eeuss影院久久| 男人舔女人下体高潮全视频| 热99在线观看视频| 七月丁香在线播放| 男女下面进入的视频免费午夜| 99久久精品热视频| 午夜激情久久久久久久| 综合色av麻豆| 一夜夜www| 91精品伊人久久大香线蕉| 午夜老司机福利剧场| 一本久久精品| 国内精品宾馆在线| 午夜激情福利司机影院| 国产美女午夜福利| 99热这里只有是精品在线观看| 欧美极品一区二区三区四区| 精品亚洲乱码少妇综合久久| 国产v大片淫在线免费观看| 久久久久久久久久久丰满| 一级毛片我不卡| 成人午夜高清在线视频| 国产成年人精品一区二区| 一级a做视频免费观看| 久久久国产一区二区| 国产精品麻豆人妻色哟哟久久 | 大又大粗又爽又黄少妇毛片口| 日韩视频在线欧美| 秋霞伦理黄片| 日韩电影二区| a级一级毛片免费在线观看| 人妻一区二区av| 熟妇人妻久久中文字幕3abv| 精品少妇黑人巨大在线播放| 国产成人freesex在线| 国产探花极品一区二区| 国内精品一区二区在线观看| 国产成人freesex在线| 成人亚洲欧美一区二区av| 在线a可以看的网站| 日本黄大片高清| 中文天堂在线官网| 男女下面进入的视频免费午夜| 日韩大片免费观看网站| 国产成人免费观看mmmm| 少妇的逼好多水| 国产精品女同一区二区软件| 久久6这里有精品| av女优亚洲男人天堂| 嘟嘟电影网在线观看| 成人高潮视频无遮挡免费网站| 夜夜爽夜夜爽视频| 久久亚洲国产成人精品v| 看免费成人av毛片| 亚洲av二区三区四区| 国产女主播在线喷水免费视频网站 | 26uuu在线亚洲综合色| 精品午夜福利在线看| 观看美女的网站| 听说在线观看完整版免费高清| 嫩草影院新地址| 禁无遮挡网站| 亚洲欧美中文字幕日韩二区| 一级黄片播放器| 亚洲国产成人一精品久久久| 精品久久久噜噜| 麻豆成人av视频| www.av在线官网国产| 国产有黄有色有爽视频| 亚洲精品乱码久久久久久按摩| 成人国产麻豆网| 国产不卡一卡二| 丰满人妻一区二区三区视频av| 国国产精品蜜臀av免费| 偷拍熟女少妇极品色| 久久国内精品自在自线图片| 韩国av在线不卡| 亚洲国产精品国产精品| 免费人成在线观看视频色| 亚洲欧美日韩东京热| 欧美区成人在线视频| 久久久午夜欧美精品| 日韩伦理黄色片| 久久久久久久久久黄片| av在线天堂中文字幕| 高清欧美精品videossex| 亚洲成人一二三区av| 中文在线观看免费www的网站| 看黄色毛片网站| 欧美最新免费一区二区三区| 婷婷色麻豆天堂久久| 高清日韩中文字幕在线| 日韩强制内射视频| 日韩av在线大香蕉| 毛片女人毛片| 亚洲欧美一区二区三区国产| 日本黄色片子视频| 日本三级黄在线观看| 国产伦精品一区二区三区四那| 日本熟妇午夜| 1000部很黄的大片| 久久久久精品久久久久真实原创| 国产av不卡久久| 午夜福利视频精品| 免费在线观看成人毛片| 久久久久久久大尺度免费视频| 天堂影院成人在线观看| av专区在线播放| 少妇熟女aⅴ在线视频| 国产精品久久久久久久电影| 午夜福利在线在线| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 九九在线视频观看精品| 男女啪啪激烈高潮av片| 久久精品久久久久久噜噜老黄| 日韩欧美三级三区| 久久久国产一区二区| 欧美精品国产亚洲| 听说在线观看完整版免费高清| 97在线视频观看| 97热精品久久久久久| 婷婷色av中文字幕| 亚洲国产精品国产精品| 亚洲在久久综合| ponron亚洲| 我的女老师完整版在线观看| 中国国产av一级| 高清欧美精品videossex| 超碰97精品在线观看| 国产高清三级在线| 久久久久国产网址| 欧美日韩综合久久久久久| 久久这里有精品视频免费| 亚洲最大成人手机在线| 美女主播在线视频| 精品久久久久久久人妻蜜臀av| 国产精品一区www在线观看| 美女xxoo啪啪120秒动态图| 日韩制服骚丝袜av| 你懂的网址亚洲精品在线观看| 国产一级毛片七仙女欲春2| 麻豆久久精品国产亚洲av| 婷婷六月久久综合丁香| 视频中文字幕在线观看| 亚洲国产色片| 亚洲精品国产av蜜桃| 看黄色毛片网站| 亚洲欧美一区二区三区黑人 | 在线观看一区二区三区| 男女边吃奶边做爰视频| 最近手机中文字幕大全| 又爽又黄a免费视频| 欧美不卡视频在线免费观看| 国产精品福利在线免费观看| 国国产精品蜜臀av免费| 国产人妻一区二区三区在| 一级毛片aaaaaa免费看小| 国产免费一级a男人的天堂| 青青草视频在线视频观看| 别揉我奶头 嗯啊视频| 精品人妻偷拍中文字幕| videossex国产| 国产免费又黄又爽又色| 日韩不卡一区二区三区视频在线| 99热这里只有是精品在线观看| 亚洲成人av在线免费| av一本久久久久| 亚洲美女搞黄在线观看| 日产精品乱码卡一卡2卡三| 大片免费播放器 马上看| 99久久精品一区二区三区| 中文欧美无线码| 午夜福利视频1000在线观看| 亚洲精华国产精华液的使用体验| 一级毛片 在线播放| 麻豆乱淫一区二区| 麻豆成人av视频| 美女黄网站色视频| 国产精品蜜桃在线观看| 日韩av在线免费看完整版不卡| 日日摸夜夜添夜夜爱| 欧美潮喷喷水| 久久精品夜夜夜夜夜久久蜜豆| 国产精品久久久久久精品电影小说 | 亚洲精品日本国产第一区| 欧美激情国产日韩精品一区| 有码 亚洲区| 99热这里只有精品一区| 国产精品国产三级专区第一集| 大片免费播放器 马上看| 人人妻人人看人人澡| 九九久久精品国产亚洲av麻豆| 狂野欧美白嫩少妇大欣赏| 国产高清不卡午夜福利| 午夜激情久久久久久久| 国产乱人偷精品视频| 日产精品乱码卡一卡2卡三| 一区二区三区乱码不卡18| av在线亚洲专区| 亚州av有码| 一级av片app| 在线天堂最新版资源| 99九九线精品视频在线观看视频| 秋霞伦理黄片| 精品午夜福利在线看| 女的被弄到高潮叫床怎么办| 久久99精品国语久久久| 亚洲精品日本国产第一区| 少妇高潮的动态图| 超碰av人人做人人爽久久| 国产成人a区在线观看| 亚洲国产精品sss在线观看| 国产老妇伦熟女老妇高清| 国产成人精品福利久久| 国产综合精华液| 国产色婷婷99| 精品国产三级普通话版| 韩国高清视频一区二区三区| 久久热精品热| 国内精品宾馆在线| 熟女电影av网| 午夜福利在线观看吧| 国产成人一区二区在线| 久久久精品94久久精品| 日韩电影二区| 亚洲精品色激情综合| 欧美 日韩 精品 国产| 国产永久视频网站| 亚洲最大成人手机在线| 欧美区成人在线视频| 搡老妇女老女人老熟妇| 欧美日本视频| 精品久久久久久久久av| 国产成人a区在线观看| 久久精品久久久久久噜噜老黄| 韩国av在线不卡| 天堂中文最新版在线下载 | 午夜福利视频1000在线观看| 麻豆成人av视频| 亚洲欧洲日产国产| 青春草视频在线免费观看| 成人国产麻豆网| 麻豆av噜噜一区二区三区| 亚洲第一区二区三区不卡| 成人午夜高清在线视频| 内射极品少妇av片p| 色哟哟·www| 中文在线观看免费www的网站| av国产免费在线观看| 最近的中文字幕免费完整| 波野结衣二区三区在线| 床上黄色一级片| 亚洲精品第二区| 国产成人精品久久久久久| 久久精品综合一区二区三区| 欧美日韩亚洲高清精品| 欧美日韩在线观看h| 水蜜桃什么品种好| av网站免费在线观看视频 | 岛国毛片在线播放| 两个人视频免费观看高清| 少妇的逼水好多| 亚洲国产欧美在线一区| 久久久久久久亚洲中文字幕| 在线免费观看的www视频| 免费黄网站久久成人精品| 亚洲欧美精品专区久久| 亚洲精品成人av观看孕妇| 永久网站在线| 婷婷色综合www| 99久久九九国产精品国产免费| 乱人视频在线观看| 在线播放无遮挡| 日本猛色少妇xxxxx猛交久久| videos熟女内射| 成人国产麻豆网| 干丝袜人妻中文字幕| av网站免费在线观看视频 | 久久久久久久久久久免费av| 欧美一区二区亚洲| 欧美日本视频| 免费看美女性在线毛片视频| 国产伦理片在线播放av一区| 深爱激情五月婷婷| 亚洲国产日韩欧美精品在线观看| 我要看日韩黄色一级片| 国产成人午夜福利电影在线观看| 欧美另类一区| 亚洲人成网站在线观看播放| 久久久久久久午夜电影| 80岁老熟妇乱子伦牲交| 久久精品综合一区二区三区| 欧美变态另类bdsm刘玥| 丰满乱子伦码专区| 国产精品人妻久久久久久| 日韩av免费高清视频| 免费黄频网站在线观看国产| 免费电影在线观看免费观看| 色尼玛亚洲综合影院| 国产不卡一卡二| 91精品国产九色| 91aial.com中文字幕在线观看| 国产黄色视频一区二区在线观看| 日韩av在线免费看完整版不卡| av女优亚洲男人天堂| 国产一级毛片七仙女欲春2| 欧美日韩国产mv在线观看视频 | 小蜜桃在线观看免费完整版高清| 国产亚洲午夜精品一区二区久久 | 免费看光身美女| 一级二级三级毛片免费看| 波多野结衣巨乳人妻| 精品久久国产蜜桃| 精品久久久久久久人妻蜜臀av| 亚洲内射少妇av| 联通29元200g的流量卡| 少妇的逼好多水| 免费观看精品视频网站| 国产精品三级大全| 人人妻人人看人人澡| 一级毛片电影观看| 欧美日本视频| 日本免费a在线| 亚洲性久久影院| 欧美性猛交╳xxx乱大交人| 国产又色又爽无遮挡免| 乱人视频在线观看| 精品少妇黑人巨大在线播放| 亚洲精品久久午夜乱码| 欧美不卡视频在线免费观看| 亚洲aⅴ乱码一区二区在线播放| 日韩,欧美,国产一区二区三区| 青青草视频在线视频观看| 国产午夜精品论理片| 综合色av麻豆| 国产伦一二天堂av在线观看| 国产精品女同一区二区软件| 不卡视频在线观看欧美| 少妇猛男粗大的猛烈进出视频 | 精品国产三级普通话版| 精品久久久久久久久久久久久| 亚洲av成人精品一二三区| 亚洲精品乱码久久久v下载方式| 在线观看一区二区三区| 国产乱人偷精品视频| 美女国产视频在线观看| 国产一区有黄有色的免费视频 | 狠狠精品人妻久久久久久综合| 日本色播在线视频| 亚洲美女搞黄在线观看| 国产精品国产三级国产av玫瑰| 男女那种视频在线观看| 777米奇影视久久| 国产探花在线观看一区二区| a级毛片免费高清观看在线播放| 久久这里有精品视频免费| 久久久精品94久久精品| 高清视频免费观看一区二区 | 夜夜看夜夜爽夜夜摸| 高清欧美精品videossex| 日本猛色少妇xxxxx猛交久久| 最近最新中文字幕大全电影3| 亚洲精品乱久久久久久| 成年女人在线观看亚洲视频 | 久久精品夜色国产| 人人妻人人澡欧美一区二区| 嘟嘟电影网在线观看| 国产一区二区在线观看日韩| 高清在线视频一区二区三区| 国产一区亚洲一区在线观看| 国产伦精品一区二区三区四那| a级毛片免费高清观看在线播放| 精品久久久久久久人妻蜜臀av| 国产视频首页在线观看| 免费看光身美女| 69av精品久久久久久| 成人亚洲精品一区在线观看 | 精品熟女少妇av免费看| 久久国内精品自在自线图片| 精品人妻偷拍中文字幕| 亚洲一区高清亚洲精品| 男人舔女人下体高潮全视频| 国产亚洲av片在线观看秒播厂 | av线在线观看网站| 97超视频在线观看视频| 爱豆传媒免费全集在线观看| 欧美潮喷喷水| 深爱激情五月婷婷| 成人高潮视频无遮挡免费网站| 男女视频在线观看网站免费| 最近中文字幕2019免费版| 午夜福利视频1000在线观看| 男女啪啪激烈高潮av片| 久久草成人影院| 热99在线观看视频| 精品人妻一区二区三区麻豆| 日韩三级伦理在线观看| 80岁老熟妇乱子伦牲交| 国产黄色小视频在线观看| 欧美激情久久久久久爽电影| 麻豆国产97在线/欧美| 小蜜桃在线观看免费完整版高清| 校园人妻丝袜中文字幕| 在线a可以看的网站| 日本欧美国产在线视频| 亚洲精品久久久久久婷婷小说| 亚洲精品456在线播放app| 男女下面进入的视频免费午夜| 午夜久久久久精精品| 国产成人精品婷婷| 亚洲美女搞黄在线观看| 色综合亚洲欧美另类图片| 丝袜喷水一区| 国产亚洲最大av| 内射极品少妇av片p| 大片免费播放器 马上看| 天堂中文最新版在线下载 | 亚洲真实伦在线观看| 久久草成人影院| 亚洲伊人久久精品综合| 免费少妇av软件| 亚洲精品国产成人久久av| 99热6这里只有精品| 国产毛片a区久久久久| 欧美成人一区二区免费高清观看| av在线蜜桃| 久久精品国产亚洲网站| 日韩不卡一区二区三区视频在线| 免费少妇av软件| 日韩强制内射视频| 免费无遮挡裸体视频| 高清欧美精品videossex| 天天躁夜夜躁狠狠久久av| 美女高潮的动态| 婷婷六月久久综合丁香| 成年女人在线观看亚洲视频 | 看非洲黑人一级黄片| 久久久久久久久久久丰满| 国产成人精品婷婷| 麻豆成人av视频| 精品久久久久久久末码| 午夜福利高清视频| 女的被弄到高潮叫床怎么办| 日韩三级伦理在线观看| 中文资源天堂在线| 午夜福利在线观看吧| 一区二区三区高清视频在线| 免费看日本二区| 高清日韩中文字幕在线| 午夜免费男女啪啪视频观看| 免费观看无遮挡的男女| 午夜福利在线观看免费完整高清在| 亚洲最大成人手机在线| 国产精品福利在线免费观看| av黄色大香蕉| 国产亚洲最大av| 精品酒店卫生间| 男人和女人高潮做爰伦理| 亚洲精品日韩在线中文字幕| 欧美成人午夜免费资源| 一级a做视频免费观看| 国产成人freesex在线| 最近手机中文字幕大全| 亚洲真实伦在线观看| 亚洲图色成人| 国产成人福利小说| 亚洲无线观看免费| 亚洲欧美日韩卡通动漫| 插逼视频在线观看| 国产精品美女特级片免费视频播放器| 午夜日本视频在线| 国产91av在线免费观看| 国产成人福利小说| 大香蕉久久网| 亚洲欧美中文字幕日韩二区| 精品久久久久久久久亚洲| 男女那种视频在线观看| 久久综合国产亚洲精品| 全区人妻精品视频| 国产亚洲一区二区精品| 亚洲怡红院男人天堂| 中文字幕久久专区| 日本爱情动作片www.在线观看| av.在线天堂| 国产中年淑女户外野战色| 91av网一区二区| 久久久久久久国产电影| 精品久久久久久久久av| 亚洲精品一二三| 国产黄片视频在线免费观看| 亚洲第一区二区三区不卡| 国产69精品久久久久777片| 99久久中文字幕三级久久日本| kizo精华| 亚洲国产最新在线播放| 2022亚洲国产成人精品| 中文字幕人妻熟人妻熟丝袜美| a级毛片免费高清观看在线播放| 女人久久www免费人成看片| 国产人妻一区二区三区在| 欧美bdsm另类| 国产伦精品一区二区三区四那| 丝袜美腿在线中文| 大陆偷拍与自拍| 麻豆乱淫一区二区| 我的女老师完整版在线观看| 日韩一区二区视频免费看| 少妇猛男粗大的猛烈进出视频 | 久久精品综合一区二区三区| 成人亚洲精品一区在线观看 | 97在线视频观看| 成人亚洲精品一区在线观看 | 日韩强制内射视频| 欧美区成人在线视频| 韩国高清视频一区二区三区| 97超视频在线观看视频| 不卡视频在线观看欧美| 亚洲国产色片| 看非洲黑人一级黄片| 午夜免费观看性视频| 欧美成人a在线观看| 亚洲精品日韩在线中文字幕| 亚洲天堂国产精品一区在线| 一级二级三级毛片免费看| 亚洲国产高清在线一区二区三| 麻豆av噜噜一区二区三区| 黑人高潮一二区| 亚洲欧美日韩东京热| 简卡轻食公司| 日韩制服骚丝袜av| 中文资源天堂在线| 99久国产av精品国产电影| 亚洲精品成人久久久久久| 美女被艹到高潮喷水动态| 国产精品久久久久久av不卡| freevideosex欧美| 在线观看人妻少妇| 嫩草影院新地址| 免费观看在线日韩| 18禁动态无遮挡网站| 一边亲一边摸免费视频| 欧美不卡视频在线免费观看| 午夜激情欧美在线| 国产成人免费观看mmmm| 欧美+日韩+精品| 在线 av 中文字幕| 久久久久精品久久久久真实原创| 久久久久久久久中文| 亚洲国产最新在线播放| 深夜a级毛片| 欧美潮喷喷水| 天天躁日日操中文字幕| 久久久午夜欧美精品| 老师上课跳d突然被开到最大视频| 免费看美女性在线毛片视频| 天堂中文最新版在线下载 | 国产大屁股一区二区在线视频| 一个人免费在线观看电影| 午夜福利网站1000一区二区三区| 中文字幕亚洲精品专区| 亚洲国产欧美在线一区| 看十八女毛片水多多多| 国产精品一二三区在线看| 欧美精品国产亚洲| .国产精品久久| 久久99热6这里只有精品| 五月玫瑰六月丁香| 女人久久www免费人成看片| 免费在线观看成人毛片| 丰满人妻一区二区三区视频av| 伊人久久精品亚洲午夜| 51国产日韩欧美| 国产 亚洲一区二区三区 | 欧美人与善性xxx| 国内精品宾馆在线| 日韩不卡一区二区三区视频在线| 国产伦理片在线播放av一区| 久久国内精品自在自线图片| 欧美日韩精品成人综合77777| 国产午夜精品一二区理论片| 一本久久精品|