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

    考慮采購成本與組合風(fēng)險的構(gòu)件優(yōu)化選擇方法

    2017-09-15 07:26:47吳志樵盧祥遠(yuǎn)牟立峰唐加福
    中國管理科學(xué) 2017年8期
    關(guān)鍵詞:排序構(gòu)件成本

    吳志樵,盧祥遠(yuǎn),牟立峰,唐加福

    (1.東北財經(jīng)大學(xué)管理科學(xué)與工程學(xué)院,遼寧 大連 116025;2.上海大學(xué)悉尼工商學(xué)院,上海 200444)

    考慮采購成本與組合風(fēng)險的構(gòu)件優(yōu)化選擇方法

    吳志樵1,盧祥遠(yuǎn)1,牟立峰2,唐加福1

    (1.東北財經(jīng)大學(xué)管理科學(xué)與工程學(xué)院,遼寧 大連 116025;2.上海大學(xué)悉尼工商學(xué)院,上海 200444)

    本文圍繞著采購成本和組合風(fēng)險這兩個重要的系統(tǒng)非功能需求,將構(gòu)件選擇問題定義為在滿足系統(tǒng)需求的情況下,為各模塊選擇具體的實(shí)施構(gòu)件,實(shí)現(xiàn)系統(tǒng)總體設(shè)計最優(yōu)。針對所建立的雙目標(biāo)整數(shù)規(guī)劃模型,結(jié)合麥克勞林展開式二階拉格朗日余項,給出了精確的關(guān)于系統(tǒng)組合風(fēng)險預(yù)測的誤差分析方法,據(jù)此得出應(yīng)用該模型估計系統(tǒng)總體風(fēng)險的適用條件。進(jìn)一步設(shè)計了求解全部有效解集的啟發(fā)式算法。并對由小到大四種規(guī)模的算例,采用衡量多目標(biāo)優(yōu)化算法的重要指標(biāo)——生成有效解向量的比率(ONVGR),比較了本文提出算法與NSGA-II算法。

    系統(tǒng)構(gòu)件選擇;優(yōu)化模型;多目標(biāo)優(yōu)化

    1 引言

    基于構(gòu)件的產(chǎn)品開發(fā)(CBD)是以構(gòu)件為中心組織整個產(chǎn)品開發(fā)過程,主要包括構(gòu)件設(shè)計、構(gòu)件選擇、構(gòu)件測試與適配、構(gòu)件更新、構(gòu)件集成及產(chǎn)品規(guī)劃設(shè)計等多階段。區(qū)別于傳統(tǒng)的從零開始開發(fā)方式,構(gòu)件選擇是通過采購市場中已有構(gòu)件來實(shí)現(xiàn)。構(gòu)件選擇既是CBD開發(fā)的基礎(chǔ)工作,也是貫穿整個開發(fā)過程的核心工作。目前針對這部分的研究多數(shù)都是關(guān)于構(gòu)件技術(shù)的實(shí)現(xiàn)細(xì)節(jié)[1], 而缺少宏觀上的決策指導(dǎo)構(gòu)件的選擇。當(dāng)建立一個由若干模塊(或稱子系統(tǒng))構(gòu)成的系統(tǒng)時,為每一模塊從一組具有相似功能屬性,而具有不同非功能屬性(如成本與組合風(fēng)險)的備選構(gòu)件集合中選擇一個構(gòu)件實(shí)施,達(dá)到系統(tǒng)總體最優(yōu)是統(tǒng)籌學(xué)研究中的關(guān)鍵問題。

    自上世紀(jì)九十年代以來,許多的優(yōu)化方法被引入構(gòu)件選擇問題中,按照優(yōu)化目標(biāo)的不同,主要分為風(fēng)險最優(yōu)問題[2]或成本最優(yōu)問題[3-6]。然而,在實(shí)際的系統(tǒng)開發(fā)過程中對構(gòu)件選擇問題,僅單獨(dú)考慮成本或者風(fēng)險的最優(yōu)解往往不是決策者真正需要的,有時決策者也并不需要單純最優(yōu)解。而協(xié)同考慮這兩個目標(biāo)的有效解卻顯得更為重要并有意義[7-8]。由于這兩個目標(biāo)間復(fù)雜的非線性關(guān)系,有學(xué)者通過設(shè)計非支配排序算法(Nondominated Sorting Genetic Algorithm, NSGA)來求解近似最優(yōu)解[9]。但亞啟發(fā)式算法NSGA本身有兩個缺陷:1)不能保證求出全部的有效解;2)求出的有效解中沒有平衡成本與風(fēng)險兩個目標(biāo)的決策信息,使決策者面臨從有效解集中選擇適合組織的有效解的困難[9-10]。

    本文通過對構(gòu)件的成本與組合風(fēng)險兩項指標(biāo)關(guān)系進(jìn)行研究,建立一個雙目標(biāo)0-1整數(shù)規(guī)劃模型,為開發(fā)者提供滿足系統(tǒng)需求情況下考慮成本與組合風(fēng)險最優(yōu)的構(gòu)件選擇方法。同時,根據(jù)模型特點(diǎn)設(shè)計基于支持有效解集的啟發(fā)式算法求解全部有效解集。采用衡量多目標(biāo)優(yōu)化算法的重要指標(biāo)——生成有效解向量的比率(ONVGR)來比較本文提出算法與NSGA-II算法。

    1.1 考慮成本與組合風(fēng)險的構(gòu)件選擇問題描述

    文中使用的參數(shù)分別是:

    n為完成目標(biāo)系統(tǒng)需要開發(fā)的模塊個數(shù);

    mj為模塊j可以選擇的構(gòu)件個數(shù);

    cij表示第j個模塊選擇第i種構(gòu)件實(shí)施時的采購成本($);

    ηij表示第j個模塊選擇第i種構(gòu)件實(shí)施時的失敗率(%);

    sj為第j個模塊在此產(chǎn)品架構(gòu)實(shí)現(xiàn)時被調(diào)用的平均次數(shù)(次);

    Satij(k)表示當(dāng)?shù)趈個模塊選擇第i種構(gòu)件開發(fā)時,第k個目標(biāo)的可滿足性水平。

    文中使用的變量為:

    不同于簡單的產(chǎn)品開發(fā),CBD可以系統(tǒng)地組織一組產(chǎn)品同時開發(fā)。從產(chǎn)品線管理角度來說,產(chǎn)品線可視為在一組共享資源基礎(chǔ)上建立的產(chǎn)品組合,產(chǎn)品組合中各產(chǎn)品分別滿足特定領(lǐng)域市場的需求。不失一般性地,為了便于理解產(chǎn)品線中的共享資源通過抽象的模塊進(jìn)行表示,并且每一模塊可以通過使用一個構(gòu)件進(jìn)行實(shí)施,而且每一模塊需要實(shí)現(xiàn)至少一個功能目標(biāo)。根據(jù)面向目標(biāo)的需求工程(The Goal-Oriented Requirement Engineering Approach,GORE)[11-12],需求經(jīng)過抽象、細(xì)分均轉(zhuǎn)化為可以量化考核的目標(biāo),這些目標(biāo)可以根據(jù)組織的需要被劃分為若干評分檔次,決策者可以通過測試、分析和匹配各構(gòu)件資產(chǎn)執(zhí)行時的效果,使用可滿足性函數(shù)對資產(chǎn)進(jìn)行度量評分。圖1表示在面向目標(biāo)的需求工程下構(gòu)件選擇的框架,以模塊為基本單元來衡量質(zhì)量目標(biāo)和構(gòu)件。

    圖1 面向目標(biāo)的需求工程下構(gòu)件選擇框架

    在此框架下,圍繞著成本和組合風(fēng)險這兩個重要的非功能需求,構(gòu)件選擇決策問題可以定義為在滿足系統(tǒng)需求的情況下,為各模塊選擇具體實(shí)施的構(gòu)件,最終實(shí)現(xiàn)總體開發(fā)成本和組合風(fēng)險最優(yōu)的目標(biāo)。

    1.2 組合風(fēng)險的描述與誤差分析

    根據(jù)Berman[2]的研究,模塊的錯誤發(fā)生概率服從指數(shù)分布。那么,假設(shè)以pj表示模塊j在系統(tǒng)執(zhí)行時至少發(fā)生一次錯誤的概率,則模塊j沒有錯誤發(fā)生的概率可以表示為:

    (1)

    (2)

    最優(yōu)化系統(tǒng)的組合風(fēng)險也就可以表示成為最小化至少一次發(fā)生錯誤的概率。這樣,目標(biāo)函數(shù)表示為:

    (3)

    假設(shè)系統(tǒng)中任意模塊出現(xiàn)錯誤,便導(dǎo)致整個系統(tǒng)錯誤。根據(jù)麥克勞林前兩階展開式,可以將原組合風(fēng)險目標(biāo)方程(3)近似表示為:

    (4)

    其中,將(3)式展開為(4)式的麥克勞林展開式的二階拉格朗日余項為:

    (5)

    下面考慮該拉格朗日余項,即(5)式趨于零的條件。

    故此,(5)式在完成目標(biāo)系統(tǒng)需要開發(fā)的模塊個數(shù)n趨于無窮的時候接近于零。而事實(shí)上為完成目標(biāo)系統(tǒng)需要開發(fā)的模塊個數(shù)n總會是有限個,并且n越大,系統(tǒng)出錯的概率就會越高,進(jìn)而導(dǎo)致組合風(fēng)險存在近似誤差。

    (6)

    假設(shè)系統(tǒng)對組合風(fēng)險估計誤差的可接受波動大小為δ,則由(5)式導(dǎo)致的誤差在δ范圍內(nèi)時,所提出的求解模型是有效的。根據(jù)上述,顯然t滿足(6)式時,(5)式取得最大值。若已知系統(tǒng)的誤差可接受波動大小δ,則:

    θ∈(0,e-1/2)

    (7)

    當(dāng)系統(tǒng)組合風(fēng)險估計誤差的可接受波動滿足(7)式時,我們認(rèn)為模型是有效的。

    1.3 系統(tǒng)需求約束的描述

    本文采用目前被業(yè)界廣泛使用的面向目標(biāo)的需求工程(GORE)[11-12]來衡量系統(tǒng)的系統(tǒng)需求。GORE主要是通過抽象和細(xì)分將需求轉(zhuǎn)化成為可以量化考核的目標(biāo),再把這些目標(biāo)根據(jù)組織的需要劃分為若干評分檔次,決策者可以通過測試、分析和匹配各構(gòu)件執(zhí)行時的效果,使用可滿足性水平反映各構(gòu)件(功能模塊)的度量評分。假設(shè)滿意度水平參數(shù)(Satisfaction Level,Satij(k))在[0,1]區(qū)間上,并且已經(jīng)通過匹配各目標(biāo)檔次得到。通過對各目標(biāo)的加權(quán)求和,可以得到第j個模塊選擇第i種構(gòu)件開發(fā)時的總體可滿足性水平Satij:

    (8)

    其中,τk表示各目標(biāo)k的權(quán)重。

    則整個系統(tǒng)的需求約束可以描述為:

    (9)

    表示整個系統(tǒng)的運(yùn)行至少需要達(dá)到某一可接受水平R。式中,ωj表示模塊j的權(quán)重。

    1.4 數(shù)學(xué)模型

    通過以上對成本、組合風(fēng)險和系統(tǒng)需求的描述,可以得到滿足系統(tǒng)需求下考慮成本和組合風(fēng)險最優(yōu)的優(yōu)化模型(Model of Cost and combinational Risk Optimization under a goal satisfaction constraint,C&R):

    (10)

    (11)

    (12)

    (13)

    xij=0/1,?i,j

    (14)

    其中,式(10)表示整個系統(tǒng)的總成本。式(13)表示對于各模塊而言,無論選擇何種構(gòu)件,只使用一種構(gòu)件實(shí)施。

    2 啟發(fā)式算法

    可以看出,C&R模型為雙目標(biāo)整數(shù)規(guī)劃模型(The Bi-objective 0-1 Integer Programming),決策者使用目前如ILOG Cplex 10、Lindo和Lingo12.0等優(yōu)化軟件不能求解,而使用如遺傳算法、模擬退火算法和分散搜索算法等智能優(yōu)化算法僅能求出部分有效解。為此,本文根據(jù)該模型的特點(diǎn)設(shè)計了三階段啟發(fā)式算法對其進(jìn)行求解,提供全部非劣解集 (Non-dominated Solution),以方便模型使用。

    2.1 啟發(fā)式算法的假設(shè)條件與相關(guān)概念

    假設(shè)系統(tǒng)中任意模塊出現(xiàn)錯誤,便導(dǎo)致整個系統(tǒng)錯誤。根據(jù)麥克勞林前兩階展開式(The First Two Terms of Maclaurin Series Expansion),可以將原組合風(fēng)險目標(biāo)方程(4)近似表示為:

    (15)

    不同于單目標(biāo)優(yōu)化模型具有唯一的最優(yōu)解,一般來講,多目標(biāo)優(yōu)化問題并不存在唯一最優(yōu)解,而是往往通過一個解的集合來表示求解結(jié)果。該集合描述的是與集合之外的解相比較,集合中的解至少有一個目標(biāo)值優(yōu)于集合之外的一個解,且其他目標(biāo)值均不比這個集合外的解劣,這種解集中的解我們稱之為非支配解(Non-dominated Solution)或者Pareto解,并將這種解集稱為有效解集(Efficient Solution)[13]。

    為了更有效地描述本算法,下面給出C&R模型中非支配解的具體定義。

    假設(shè)引入一個在 [0,1]之間取值的偏好值λ(Preferences),通過對兩個原目標(biāo)函數(shù)加權(quán)求和,構(gòu)成一個新的目標(biāo)函數(shù):

    z(x,λ)=λzCOS(x)+(1-λ)zREL(x)

    (16)

    這個偏好值反應(yīng)決策者對成本與組合風(fēng)險的偏好信息,通過將原來的兩個目標(biāo)加權(quán)成為一個目標(biāo),再求解該優(yōu)化模型,可以得到原問題的一個有效解,并且這一有效解反應(yīng)了決策者事先定義好的對成本與組合風(fēng)險的偏好信息。

    從理論上講,如果可以遍歷偏好值λ在[0,1]之間的取值,并以各加權(quán)后的目標(biāo)函數(shù)作為新的目標(biāo),在滿足約束(12)~(14)的條件下對模型進(jìn)行求解,在凸可行域上就可以得到全部非支配解。然而這種方法存在兩個局限性:第一,遍歷所有偏好值λ的取值在實(shí)際中是不能實(shí)現(xiàn)的;第二,對于非凸可行域的問題,使用該方法不能求解出全部非支配解。為此,本文提出一種可以獲得全部非支配解的算法。通過在滿足約束(12)~(14)下,最小化加權(quán)后的目標(biāo)函數(shù)(16)可以獲得部分有效解,這部分有效解我們將其稱之為支持的有效解(Supported Efficient Solutions,Supported-ES);同時,將有效解集中的其他有效解稱之為非支持的有效解(Non-supported Efficient Solutions,Non- supported-ES)。此算法應(yīng)用非支持的有效解必存在于兩個連續(xù)的支持有效解之間的理論[14],基于此理論,在本文提出的三階段算法中:第一階段對全部可能的聯(lián)合效用值(Aggregated Efficiency)進(jìn)行排序,由于排序的變化取決于偏好λ的取值,得到排序穩(wěn)定的λ的區(qū)間范圍,以此有效地解決了遍歷全部偏好λ的取值的困境;第二階段通過引入貪婪啟發(fā)式規(guī)則,遍歷得出全部支持的有效解;第三階段中,通過遍歷任意兩個連續(xù)的支持有效解構(gòu)成的區(qū)間,得到全部非支持的有效解。至此,全部的有效解就得到了。

    為了更清楚地說明算法的求解過程,引入如下定義:

    通過上述定義可以進(jìn)一步得到考慮帶偏好值的聯(lián)合效率值(Aggregated Efficiency)eij(λ)為:

    eij(λ)=αijλ+βij(1-λ)=βij+(αij-βij)λ

    (17)

    它表示第j個模塊選擇第i種構(gòu)件實(shí)施時,每單位滿足度增加所付出的成本與組合風(fēng)險聯(lián)合作用。

    因此,對于一個給定的λ′,可以得到其聯(lián)合效用值排序:

    el1j(λ′)≤…≤elij(λ′)≤eli+1j(λ′)≤…≤elmjj(λ′),

    where=<1,2,…,mj>

    (18)

    遍歷所有可能存在的聯(lián)合效用值的排序,就得到所有的確定支持有效解,相應(yīng)的每一對應(yīng)有效解的排序,我們稱之為有效排序(Efficient Order,EO)。

    2.2 第一階段:確定聯(lián)合效用值的有效排序

    本小節(jié)中將確定所有可能的有效排序(EO),以及各EO的對應(yīng)偏好值λ的區(qū)間。隨著偏好值λ在可行域[0, 1]的變化,EO不是固定不變的,首先介紹確定各模塊有效排序的λ子區(qū)間([0,u1j],…,[utj,ut+1j],…,[uTj,1]))的算法。

    第一階段算法的流程如下:

    開始

    第1步:初始化t←0,λt←0,utj←0;

    第2步:以λt的取值,將所有構(gòu)件按照式(18)進(jìn)行排序;

    第3步:通過求解下面優(yōu)化問題確定有效的偏好值的上界(upper bounds)

    Maxλ

    s.t.elij(λ)≤eli+1j(λ)

    i=1,…,mj-1,0≤λ≤1

    (19)

    假設(shè)用λmax來表示(19)模型的最優(yōu)解,如果λmax≥1,則令t←t+1,utj←1

    即表示模塊j的全部有效的偏好值區(qū)間已經(jīng)確定完畢,進(jìn)行第4步;

    第4步: 令t←t+1,utj←λmax,λt←λmax+ε(其中ε為干擾項(perturbation factor),ε>0)。進(jìn)行第2步;

    結(jié)束

    此外,式(19)可以使用αij和βij參數(shù),通過如下優(yōu)化模型(20)進(jìn)行求解:

    (αi+1j-βi+1j-(αij-βij))<0,i=1,…,mj}

    (20)

    推論1 各模塊的有效排序個數(shù)等于模型(20)在從0遍歷到1的過程中所取得的偏好值λmax的個數(shù)。

    根據(jù)第一階段啟發(fā)式(Procedure-1)可以得出各模塊的有效排序及其對應(yīng)的子區(qū)間,并通過合并全部子區(qū)間可以得出有效排序不改變的偏好值區(qū)間([0,v1],…, [vg,vg+1],…, [vG,1])。

    2.3 第二階段:求解支持有效解

    對于每一個有效排序不變的偏好值λ的區(qū)間,通過貪婪啟發(fā)式過程求解出支持有效解,遍歷全部區(qū)間便可以得到全部的支持有效解。

    第二階段算法的流程如下:

    開始

    第1步:初始化h←0,g←0,λh←0,EfficientSet←{φ},ExploreSet←{φ};

    第3步:循環(huán)直到滿足系統(tǒng)需求約束(12)為止:

    如果xli+1j=non,則將xlij鎖定為必選解;

    否則對ExploreSet∪{xli+1j},EfficientSet←{xli+1j}使用構(gòu)件xli+1j更新}

    第4步:如果滿足Dantzig檢驗(yàn)條件 (17);則保存支持到支持有效解集中,

    并初始化EfficientSet←{φ},ExploreSet←{φ}。

    第5步:λh=vg+1+ε(其中ε為干擾項(perturbation factor),ε>0);

    如果λh≥1, 則終止算法(說明所有空間已經(jīng)搜索完畢);,

    否則返回第2步;

    結(jié)束

    在此過程中需要用到的Dantzig規(guī)則如下:

    (21)

    s.t.ρ1j≥βij+(αij-βij)λ,i∈EfficientSet

    ρ1j≥0,ρ2j≥0,1≥λ≥0

    使用z*和λ*表示優(yōu)化模型(21)的最優(yōu)解。根據(jù)Dantzig 規(guī)則,可以進(jìn)一步得到推論2。

    通過對每一模塊使用(Procedure-2),可以獲得全部支持有效解,并且每一支持有效解都會帶有相應(yīng)的成本與組合風(fēng)險的偏好信息(Trade-off Information ofλ),以輔助決策者進(jìn)行有效解選擇。

    2.4 第三階段:求解非支持有效解

    而對于剩余的有效解,也就是之前提到的非支持有效解必存在于兩個連續(xù)的支持有效解之間的理論,將通過本階段啟發(fā)式算法進(jìn)一步搜索,以確保全部有效解被遍歷。

    第三階段算法的流程如下:

    開始

    第1步:初始化支持有效解集(Supported-ES set)

    第2步:構(gòu)造搜索空間,空間中的下界(Lower bound)由有效解的連線構(gòu)成,而上界(Upper bound)為將C&R模型松弛為0-1連續(xù)背包問題后有效解的連線。

    第3步:通過建立典型搜索技術(shù)分支定界算法(Branch-and-bound)[14]或者動態(tài)規(guī)劃算法(Dynamic programming)[15]搜索各目標(biāo)值下降空間(Reduced objective Space),當(dāng)遍歷過所有解空間之后,全部的非支持有效解也就確定完畢。

    結(jié)束

    此外,加權(quán)目標(biāo)函數(shù)(Aggregated Functions)可以進(jìn)一步用來確定分支定界搜索(Branch-and-Bound)下降空間。

    特別地,為了方便大規(guī)模問題的求解,由于各模塊的處理過程間具有平行關(guān)系(Procedure-1),因此可以在不同的處理器上進(jìn)行運(yùn)算,最后通過合成得到有效排序穩(wěn)定的區(qū)間。

    3 案例分析

    本節(jié)以郵箱系統(tǒng)為例來驗(yàn)證前文所提優(yōu)化模型及其啟發(fā)式算法在實(shí)際應(yīng)用中的有效性。

    一個郵件數(shù)據(jù)請求從發(fā)件人所處的用戶代理服務(wù)器(Mail User Agent,MUA),可以由四個模塊構(gòu)成:Web Mail Servers、Wap Mail Servers、Secure Mail、Network Communication(為方便表述,后文將其分別縮寫為模塊1、模塊2、模塊3、模塊4)??紤]到產(chǎn)品的多樣性,產(chǎn)品1、產(chǎn)品2、產(chǎn)品3三類可如圖2進(jìn)行開發(fā)。

    圖中的連線表示不同產(chǎn)品開發(fā)時需要用到的模塊。例如開發(fā)產(chǎn)品1,將用到Web Mail Servers和 Wap Mail Servers模塊作為產(chǎn)品的構(gòu)成模塊。

    各構(gòu)件所需操作成本均使用成本表示,其可以通過平均的工時成本(Man-hours)來計算。求解出各模塊在不同的可選擇構(gòu)件下所對應(yīng)的成本參數(shù)cij、及組合風(fēng)險值ηij,結(jié)果見表1。并根據(jù)各模塊的使用頻率得到各模塊的權(quán)重值為(0.472, 0.311, 0.087, 0.066)。

    圖2 郵件系統(tǒng)產(chǎn)品線架構(gòu)圖

    構(gòu)件序號模塊1模塊2模塊3模塊4成本組合風(fēng)險滿意度成本組合風(fēng)險滿意度成本組合風(fēng)險滿意度成本組合風(fēng)險滿意度110005.10.7118007.20.6110004.60.671155.10.5124207.10.734208.70.724206.50.74855.30.5834606.80.634608.10.664106.80.68954.90.6244106.90.645107.10.683606.90.791004.60.7358504.80.6514505.60.558005.50.751553.70.5765455.20.8310156.10.746425.80.8312150.710.9676805.10.9118645.50.913997.10.785293.40.88

    表2描述了當(dāng)R=0.8時得到的支持有效解,每個解由編碼給出。編碼“7716”表示四個模塊分別使用的是構(gòu)件7、構(gòu)件7、構(gòu)件1、構(gòu)件6。圖3反映的是支持有效解在成本與組合風(fēng)險這兩個目標(biāo)空間上的情況。

    表2 啟發(fā)式算法獲得的支持有效解

    圖3 R=0.8時得到的支持有效解

    4 算法性能的評估

    算法性能的評估是多目標(biāo)優(yōu)化問題重點(diǎn)研究的內(nèi)容之一。本文針對從小到大4個規(guī)模的問題, 采用衡量多目標(biāo)優(yōu)化算法的重要指標(biāo)——生成有效解向量的比率(Overall Nondominated Vector Generation Ratio, ONVGR)來比較本文提出算法與NSGA-II算法。本文用PFt表示本文提出算法獲得的有效解的個數(shù),其值越大說明決策者有越多的選擇。類似地,用PFn來表示NSGA-II算法獲得的有效解的個數(shù);則生成效解數(shù)比率可表示為ONVGR=PFn/PFt。

    首先隨機(jī)生成四組算例,算例規(guī)模分別為10×10、 25×25、25×50和50×100(模塊數(shù)×備選構(gòu)件數(shù)), 縮寫為P-1、P-2、P-3和P-4。它們各自的成本, 失敗率和可滿足性水平均在區(qū)間[100, 5000]、[0.0001, 0.001]和[0.5, 1]內(nèi)隨機(jī)生成, 并將可接受可滿足性水平設(shè)置為0.6。NSGA-II算法的種群數(shù)設(shè)置為100, 迭代次數(shù)設(shè)置為200。本文所列出結(jié)果為NSGA-II算法執(zhí)行10次的平均值。為了客觀地比較上述兩種算法,本文采用兩種NSGA-II參數(shù)設(shè)置策略: 1) 交叉率設(shè)置為0.5、0.7、0.9,同時變異率設(shè)置為0.03;2)變異率設(shè)置為 0.01、0.03、0.05,同時交叉率設(shè)置為0.9。

    上述兩個算法比較結(jié)果如表3和表4所示,依據(jù)ONVGR指標(biāo),NSGA-II算法有效性僅相當(dāng)于本文提出算法71%~96%。而且NSGA-II的不同變異率和交叉率下的ONVGR比率均劣于本文提出算法。

    表3 NSGA-II.不同變異率下ONVGR比率

    表4 NSGA-II.不同交叉率下ONVGR比率

    5 結(jié)語

    本文通過建立了一個雙目標(biāo)0-1整數(shù)規(guī)劃模型,提供滿足系統(tǒng)需求情況下的成本最小且組合風(fēng)險最低的解決方案。通過分析麥克勞林展開式的二階拉格朗日余項,給出了該解決方案的適用條件,指出在系統(tǒng)組合風(fēng)險估計誤差的波動在一定可接受范圍內(nèi)時,該估計模型是有效的。并設(shè)計三階段啟發(fā)式算法對其進(jìn)行求解,前兩階段確定支持有效解集、第三階段得到非支持有效解集,以提供全部非劣解集,用于從非劣解集中選擇最適合方案實(shí)施。最后,通過算例分析驗(yàn)證了該模型和算法的適用性。

    [1] 邵良杉, 張照. 組件、COM與Windows DNA技術(shù)在MIS中應(yīng)用研究[J]. 中國管理科學(xué),2000,8(S1): 131-135.

    [2] Berman O,Cutler M. Optimal software implementation considering reliability and cost[J]. Computers & Operations Research,1997,25(10): 857-868.

    [3] Cortellessa V,Marinelli F,Potena P. Automated selection of software components based on cost/reliability tradeoff[C]//Proceedings of European Workshop on Software Architectwre,2006,4344:66-81.

    [4] Cortellessa V,Marinelli F, Potena P. An optimization framework for "build-or-buy" decisions in software architecture[J]. Computers & Operations Research,2008,35(10): 3090-3106.

    [5] Boehm B W,Valerdi R. Achievements and challenges in Cocomo-based software resource estimation[J]. IEEE Software,2008,25(5):74-83.

    [6] Wu Zhiqiao, Kwong C K. Integrated model for software component selection with simultaneous consideration of implementation and verification[J]. Computers & Operations Research, 2012,39(12):3376-3393.

    [7] Wang Zai, Tang Ke, Yao Xin. Multi-objective approaches to optimal testing resource allocation in modular software systems[J]. IEEE Transactions on reliablity,2010, 59 (3):563-575.

    [8] Zhang Yuanyuan, Harman M, Mansouri A. The multi-objective next release problem[C]//Proceedings of the 9th Annual Conference on Genetic and Evolutionary Computation,ACM,London,England,July 7-11,2007.

    [9] Stummer C,Kiesling E. A multicriteria decision support system for competence-driven project portfolio selection[J]. International Journal of Information Technology & Decision Making, 2009,8(2):379-401.

    [10] Wu Zhiqiao, Tang Jiafu,Kwong C K, et al. A model and its algorithm for software reuse optimization problem with simultaneous reliability and cost consideration[J]. International Journal of Innovative Computing, Information and Control,2011,7(5):2611-2622.

    [11] Lamsweerde A V,Letier E. Handling obstacles in goal-oriented requirements engineering[J]. IEEE Transactions on software engineering,2000,26(10):978-1005.

    [12] Letier E,Kramer J,Magee J,et al. Deriving event-based transition systems from goal-oriented requirements models[J]. Automated Software Engineering,2008,15(2):175-206.

    [13] Chen Jiangyong, Lin Qiuzhen, Hu Qingbin. Application of novel clonal algorithm in multiobjective optimization[J]. International Journal of Information Technology & Decision Making,2010, 9(2):239-266.

    [14] Visee M,Teghem J,Pirlot M, et al. Two-phases method and branch and bound procedures to solve the bi-objective knapsack problem[J]. Journal of Global Optimization,1998,12(2):139-155.

    [15] Bazgan C, Hugot H, Vanderpooten D. Solving efficiently the 0-1 multi-objective knapsack problem[J]. Computers & Operations Research,2009,36(1):260-279.

    An Optimization Model for System Component Selection to Minimize Cost and Combinational Risk

    WU Zhi-qiao1,LU Xiang-yuan1, MU Li-feng2, TANG Jia-fu1

    (1.College of Management Science and Engineering,Dongbei University of Finance and Economics,Dalian 116025, China;2.SHU-UTS SILC Business School, Shanghai University, Shanghai 200444,China)

    In contrast to the traditional process of product development, Component-Based Development (CBD) focuses on building products by integrating previously-existing components. So to start with, an available set of components should be identified. In this paper, this major problem of component-based system development involves the effective evaluation and selection alternative system components is addressed by considered the cost and combinational risk factors. Based on a bi-objective 0-1 integer programming, an optimization model is proposed that can assist decision-makers in selecting system components for minimizing cost and combinational risk, and satisfying system requirements. The condition of application of the proposed model is further proposed based on the Maclaurin expansion with second order Lagrange remained term. To solve the model efficiently, a supported efficient solution based algorithm is presented that can find the entire set of efficient solutions. Computational experience also describes in solving the problems using the Metaheuristics and the proposed algorithm.

    component selection;optimization model;multiple objective programming

    1003-207(2017)08-0158-08

    10.16381/j.cnki.issn1003-207x.2017.08.017

    2015-05-25;

    2015-11-24

    國家自然科學(xué)基金資助項目(71301107,71671028)

    吳志樵(1981-),男(錫伯族),遼寧沈陽人,東北財經(jīng)大學(xué)科學(xué)與工程學(xué)院副教授,博士生導(dǎo)師,博士,研究方向:服務(wù)系統(tǒng)工程,E-mail: wuzhiqiao@dufe.edu.cn.

    TP311.5

    A

    猜你喜歡
    排序構(gòu)件成本
    排序不等式
    2021年最新酒駕成本清單
    河南電力(2021年5期)2021-05-29 02:10:00
    恐怖排序
    節(jié)日排序
    溫子仁,你還是適合拍小成本
    電影(2018年12期)2018-12-23 02:18:48
    刻舟求劍
    兒童繪本(2018年5期)2018-04-12 16:45:32
    建筑構(gòu)件
    建筑構(gòu)件
    建筑構(gòu)件
    建筑構(gòu)件
    亚洲 国产 在线| 一级黄片播放器| 深爱激情五月婷婷| 日韩欧美三级三区| 五月伊人婷婷丁香| 日本黄色视频三级网站网址| 久久久久久久久久成人| 久久人人精品亚洲av| 亚洲av二区三区四区| 小蜜桃在线观看免费完整版高清| 一本精品99久久精品77| 99在线视频只有这里精品首页| 亚洲熟妇熟女久久| 中文亚洲av片在线观看爽| 国产精品乱码一区二三区的特点| 波多野结衣巨乳人妻| 怎么达到女性高潮| 国产单亲对白刺激| 免费在线观看亚洲国产| 成年人黄色毛片网站| 99riav亚洲国产免费| 成人午夜高清在线视频| 午夜福利免费观看在线| 麻豆国产97在线/欧美| 国产国拍精品亚洲av在线观看| 国产欧美日韩一区二区三| 亚洲精品成人久久久久久| 国产三级中文精品| 国产高清视频在线播放一区| 亚洲专区国产一区二区| 久久久久国产精品人妻aⅴ院| 中文资源天堂在线| 国产激情偷乱视频一区二区| 精品人妻一区二区三区麻豆 | 有码 亚洲区| 国产精品,欧美在线| 亚洲精品成人久久久久久| av国产免费在线观看| 日韩精品青青久久久久久| 日韩欧美 国产精品| 国产精品久久久久久久电影| 免费av不卡在线播放| 激情在线观看视频在线高清| 在线国产一区二区在线| 精品久久久久久久人妻蜜臀av| 国产av麻豆久久久久久久| 美女被艹到高潮喷水动态| 日韩大尺度精品在线看网址| 一本精品99久久精品77| 观看免费一级毛片| av在线老鸭窝| 午夜久久久久精精品| 直男gayav资源| 日本三级黄在线观看| 欧洲精品卡2卡3卡4卡5卡区| 亚洲不卡免费看| 婷婷丁香在线五月| 婷婷丁香在线五月| 国产 一区 欧美 日韩| 在线播放国产精品三级| 亚洲av电影不卡..在线观看| 校园春色视频在线观看| 亚洲精品456在线播放app | 中亚洲国语对白在线视频| 欧美zozozo另类| 日韩av在线大香蕉| 久久6这里有精品| 亚洲av熟女| 欧美激情在线99| 夜夜爽天天搞| 中国美女看黄片| 久久久久国内视频| 午夜免费激情av| 久9热在线精品视频| 国产淫片久久久久久久久 | 黄色丝袜av网址大全| 丰满人妻熟妇乱又伦精品不卡| 国产在视频线在精品| 国产在视频线在精品| 乱码一卡2卡4卡精品| 人妻丰满熟妇av一区二区三区| 高潮久久久久久久久久久不卡| 成人美女网站在线观看视频| 欧美黑人巨大hd| a级毛片a级免费在线| 亚洲狠狠婷婷综合久久图片| 51午夜福利影视在线观看| 一区二区三区四区激情视频 | 国产主播在线观看一区二区| 亚洲自拍偷在线| 白带黄色成豆腐渣| 国产黄片美女视频| 女同久久另类99精品国产91| 在线a可以看的网站| 夜夜看夜夜爽夜夜摸| 免费在线观看亚洲国产| 很黄的视频免费| 桃色一区二区三区在线观看| 老司机午夜福利在线观看视频| 夜夜夜夜夜久久久久| 亚洲人成网站在线播| 日韩欧美在线二视频| 国产精品永久免费网站| 天堂影院成人在线观看| 亚洲在线自拍视频| 美女被艹到高潮喷水动态| 日韩欧美在线乱码| 两性午夜刺激爽爽歪歪视频在线观看| 亚洲va日本ⅴa欧美va伊人久久| 白带黄色成豆腐渣| 动漫黄色视频在线观看| 午夜福利高清视频| av天堂在线播放| 搡老岳熟女国产| 久久婷婷人人爽人人干人人爱| 宅男免费午夜| 欧美性猛交黑人性爽| 露出奶头的视频| 观看美女的网站| 久久6这里有精品| 日韩高清综合在线| av女优亚洲男人天堂| 老司机午夜十八禁免费视频| 亚洲一区高清亚洲精品| 色av中文字幕| 国产精品一区二区免费欧美| 欧美最黄视频在线播放免费| 国产在线男女| 色5月婷婷丁香| 国产精品国产高清国产av| 麻豆成人av在线观看| 亚洲第一区二区三区不卡| 国产精品一区二区性色av| 此物有八面人人有两片| 国产男靠女视频免费网站| 少妇被粗大猛烈的视频| 国产免费av片在线观看野外av| 国产成年人精品一区二区| 国产三级黄色录像| 欧美性感艳星| 啪啪无遮挡十八禁网站| 禁无遮挡网站| 99热这里只有是精品在线观看 | 永久网站在线| 国产午夜精品久久久久久一区二区三区 | 内射极品少妇av片p| 又紧又爽又黄一区二区| 欧美日韩黄片免| 嫩草影视91久久| 欧美bdsm另类| 欧美精品国产亚洲| 国内毛片毛片毛片毛片毛片| 女人十人毛片免费观看3o分钟| 十八禁人妻一区二区| 久99久视频精品免费| 久久草成人影院| 国产在线精品亚洲第一网站| 午夜视频国产福利| 天天一区二区日本电影三级| 欧美三级亚洲精品| 午夜免费成人在线视频| 久久精品国产亚洲av香蕉五月| 国产一区二区亚洲精品在线观看| 女人十人毛片免费观看3o分钟| 给我免费播放毛片高清在线观看| 一进一出抽搐动态| 免费看日本二区| 国产成人a区在线观看| 色5月婷婷丁香| 一区二区三区激情视频| 97碰自拍视频| 国产精品一及| 我要看日韩黄色一级片| 天天躁日日操中文字幕| 18禁黄网站禁片午夜丰满| 欧美日韩亚洲国产一区二区在线观看| av黄色大香蕉| netflix在线观看网站| 激情在线观看视频在线高清| ponron亚洲| 草草在线视频免费看| 亚洲av成人av| 天美传媒精品一区二区| 亚洲欧美日韩高清专用| 色视频www国产| 最近最新中文字幕大全电影3| 99久久99久久久精品蜜桃| 欧美日韩综合久久久久久 | 啦啦啦观看免费观看视频高清| 91狼人影院| 国产午夜福利久久久久久| 搡老岳熟女国产| 国产国拍精品亚洲av在线观看| 性色av乱码一区二区三区2| 欧美xxxx黑人xx丫x性爽| 免费在线观看成人毛片| 最近最新中文字幕大全电影3| 国产亚洲精品综合一区在线观看| 免费无遮挡裸体视频| 国产精品1区2区在线观看.| 国产不卡一卡二| 97热精品久久久久久| av在线老鸭窝| 一进一出抽搐gif免费好疼| 少妇的逼好多水| 一区二区三区免费毛片| 十八禁国产超污无遮挡网站| 老熟妇仑乱视频hdxx| 国产一区二区亚洲精品在线观看| 99热6这里只有精品| 国产成+人综合+亚洲专区| 日韩欧美精品v在线| av视频在线观看入口| 亚洲av不卡在线观看| 国产成人a区在线观看| 一区二区三区高清视频在线| 亚洲国产欧美人成| 简卡轻食公司| 久久久久久久久中文| av在线观看视频网站免费| 亚洲av熟女| 亚洲精品乱码久久久v下载方式| 免费av观看视频| 在线天堂最新版资源| 亚洲国产精品sss在线观看| www.999成人在线观看| av视频在线观看入口| 欧美日韩福利视频一区二区| 亚洲人成网站在线播| 真人做人爱边吃奶动态| 在线天堂最新版资源| 一级a爱片免费观看的视频| 男女床上黄色一级片免费看| 丝袜美腿在线中文| 亚洲人成伊人成综合网2020| 国产在视频线在精品| 国产免费男女视频| 欧美又色又爽又黄视频| 欧美午夜高清在线| 欧美zozozo另类| x7x7x7水蜜桃| 看免费av毛片| 少妇高潮的动态图| 久久人人精品亚洲av| 长腿黑丝高跟| 在线观看午夜福利视频| 国产69精品久久久久777片| 99精品在免费线老司机午夜| 精品福利观看| 亚洲国产精品999在线| 噜噜噜噜噜久久久久久91| 欧美色视频一区免费| 色尼玛亚洲综合影院| 成人永久免费在线观看视频| 成人av一区二区三区在线看| 免费av毛片视频| 91在线精品国自产拍蜜月| 国产中年淑女户外野战色| 久久国产精品人妻蜜桃| 国产黄片美女视频| 久久国产乱子伦精品免费另类| 色哟哟·www| 午夜免费激情av| 日韩欧美在线二视频| 一本综合久久免费| 日韩人妻高清精品专区| 哪里可以看免费的av片| 18美女黄网站色大片免费观看| 老司机午夜福利在线观看视频| 中文资源天堂在线| 久久热精品热| 一二三四社区在线视频社区8| 国模一区二区三区四区视频| 老司机福利观看| 色5月婷婷丁香| 99国产极品粉嫩在线观看| 哪里可以看免费的av片| 嫁个100分男人电影在线观看| 丰满人妻熟妇乱又伦精品不卡| 国产伦在线观看视频一区| 国产男靠女视频免费网站| 国产成人福利小说| 美女cb高潮喷水在线观看| 亚洲男人的天堂狠狠| 亚洲无线观看免费| 91午夜精品亚洲一区二区三区 | 久久这里只有精品中国| 在线观看美女被高潮喷水网站 | 精品一区二区三区av网在线观看| 亚洲精品一卡2卡三卡4卡5卡| 少妇被粗大猛烈的视频| 亚洲av一区综合| 99久久精品国产亚洲精品| 淫妇啪啪啪对白视频| 成人特级黄色片久久久久久久| 男女做爰动态图高潮gif福利片| 超碰av人人做人人爽久久| 中文字幕精品亚洲无线码一区| 亚洲av成人精品一区久久| 国产极品精品免费视频能看的| 国产精品,欧美在线| 男插女下体视频免费在线播放| 日韩欧美国产在线观看| 婷婷精品国产亚洲av| 亚洲精品成人久久久久久| 嫩草影院新地址| 国产日本99.免费观看| 久久久久久久久中文| 97碰自拍视频| 日韩欧美免费精品| 色视频www国产| 搡老熟女国产l中国老女人| 午夜福利在线在线| 一级毛片久久久久久久久女| 97超级碰碰碰精品色视频在线观看| 波多野结衣巨乳人妻| 午夜福利高清视频| 亚洲av日韩精品久久久久久密| 三级毛片av免费| 色综合婷婷激情| 午夜影院日韩av| 午夜亚洲福利在线播放| 欧美中文日本在线观看视频| 在线观看舔阴道视频| 欧美日本亚洲视频在线播放| 国内精品一区二区在线观看| 9191精品国产免费久久| 亚洲精品亚洲一区二区| 九色成人免费人妻av| 1024手机看黄色片| 免费人成在线观看视频色| 两个人的视频大全免费| 久久精品影院6| 美女 人体艺术 gogo| a级一级毛片免费在线观看| 好男人电影高清在线观看| 九九热线精品视视频播放| 能在线免费观看的黄片| 看黄色毛片网站| 久久久久亚洲av毛片大全| 亚洲性夜色夜夜综合| 在线看三级毛片| 欧美日韩国产亚洲二区| 国产av在哪里看| 免费看美女性在线毛片视频| 国产v大片淫在线免费观看| 黄色配什么色好看| avwww免费| 老鸭窝网址在线观看| 高清在线国产一区| 日本成人三级电影网站| 88av欧美| 欧美日韩福利视频一区二区| 日韩欧美在线乱码| 国产午夜精品久久久久久一区二区三区 | 精品乱码久久久久久99久播| 久久九九热精品免费| 脱女人内裤的视频| 国产爱豆传媒在线观看| 国产精品爽爽va在线观看网站| 久久99热6这里只有精品| 欧美绝顶高潮抽搐喷水| 国产不卡一卡二| 久久亚洲真实| 国产综合懂色| av在线蜜桃| 噜噜噜噜噜久久久久久91| 国产亚洲欧美98| 欧美精品啪啪一区二区三区| 国产主播在线观看一区二区| 国模一区二区三区四区视频| 99热这里只有是精品50| 午夜福利高清视频| 制服丝袜大香蕉在线| 精品日产1卡2卡| 亚洲av二区三区四区| 88av欧美| 51国产日韩欧美| 欧美日韩国产亚洲二区| 国内揄拍国产精品人妻在线| 最后的刺客免费高清国语| 中文字幕久久专区| 99国产极品粉嫩在线观看| 亚洲av熟女| 乱码一卡2卡4卡精品| 亚洲一区二区三区色噜噜| 午夜福利免费观看在线| 免费看日本二区| 亚洲 欧美 日韩 在线 免费| 精品久久久久久久久av| 国产视频内射| 国产精品国产高清国产av| 久久久久久久久大av| 国产欧美日韩一区二区精品| 搡老熟女国产l中国老女人| av欧美777| 最近在线观看免费完整版| 韩国av一区二区三区四区| 久久精品国产自在天天线| 高清毛片免费观看视频网站| 啦啦啦观看免费观看视频高清| 久久久久亚洲av毛片大全| 国产精品久久久久久人妻精品电影| 亚洲avbb在线观看| 全区人妻精品视频| 一个人免费在线观看的高清视频| 亚洲欧美精品综合久久99| 全区人妻精品视频| 90打野战视频偷拍视频| 久久99热6这里只有精品| 网址你懂的国产日韩在线| 波多野结衣高清作品| 最新在线观看一区二区三区| 亚洲色图av天堂| 一二三四社区在线视频社区8| 日韩亚洲欧美综合| 91狼人影院| 欧美午夜高清在线| 非洲黑人性xxxx精品又粗又长| 成人国产综合亚洲| 亚洲18禁久久av| 国产精品亚洲av一区麻豆| 成人国产综合亚洲| 十八禁网站免费在线| 国产中年淑女户外野战色| 亚洲av五月六月丁香网| 亚洲熟妇熟女久久| 久久久国产成人精品二区| 免费av不卡在线播放| 国产精品三级大全| 成年女人毛片免费观看观看9| 成人特级黄色片久久久久久久| 舔av片在线| 久久伊人香网站| 91久久精品电影网| 99国产精品一区二区蜜桃av| 国产成人a区在线观看| 亚洲自拍偷在线| 桃红色精品国产亚洲av| 中文资源天堂在线| 香蕉av资源在线| 午夜激情福利司机影院| av福利片在线观看| 亚洲无线在线观看| 欧美国产日韩亚洲一区| 欧美+亚洲+日韩+国产| 亚洲七黄色美女视频| 欧美最黄视频在线播放免费| 亚洲国产精品999在线| 日韩亚洲欧美综合| 国内精品一区二区在线观看| 深夜精品福利| 又爽又黄无遮挡网站| 欧美三级亚洲精品| 国产伦人伦偷精品视频| 国产大屁股一区二区在线视频| 91午夜精品亚洲一区二区三区 | 综合色av麻豆| 一本精品99久久精品77| 人人妻人人澡欧美一区二区| a在线观看视频网站| 在线播放无遮挡| 深爱激情五月婷婷| 免费在线观看成人毛片| av天堂中文字幕网| 亚洲av电影在线进入| 亚洲av电影在线进入| 国产人妻一区二区三区在| 亚洲国产精品合色在线| 国产av不卡久久| 国产精品久久电影中文字幕| 亚洲精品亚洲一区二区| 嫩草影院入口| 成人精品一区二区免费| 国产蜜桃级精品一区二区三区| 久久99热6这里只有精品| 亚洲人成网站在线播放欧美日韩| 综合色av麻豆| 老女人水多毛片| 精华霜和精华液先用哪个| 1024手机看黄色片| h日本视频在线播放| 亚洲av电影在线进入| av在线老鸭窝| 国产一区二区在线观看日韩| 国产美女午夜福利| 中文亚洲av片在线观看爽| 国产一区二区三区在线臀色熟女| 久久久久久久亚洲中文字幕 | x7x7x7水蜜桃| 欧美日韩黄片免| 男女视频在线观看网站免费| 人人妻人人看人人澡| 精品乱码久久久久久99久播| 老司机午夜福利在线观看视频| 波多野结衣高清无吗| 亚洲五月婷婷丁香| 欧美精品啪啪一区二区三区| 哪里可以看免费的av片| 久久久久精品国产欧美久久久| 久久人人精品亚洲av| 能在线免费观看的黄片| 久久国产精品影院| 精品免费久久久久久久清纯| 天堂av国产一区二区熟女人妻| 亚洲综合色惰| 国产精品久久久久久久久免 | 一进一出好大好爽视频| 蜜桃亚洲精品一区二区三区| 久久热精品热| 免费观看的影片在线观看| 99久久99久久久精品蜜桃| 亚洲精品粉嫩美女一区| 99久久成人亚洲精品观看| 深夜a级毛片| 日韩亚洲欧美综合| 在线天堂最新版资源| 久久6这里有精品| 免费av观看视频| 成年女人永久免费观看视频| 国产成人a区在线观看| 草草在线视频免费看| 老司机午夜福利在线观看视频| av在线观看视频网站免费| 亚洲国产色片| 亚洲人成网站在线播放欧美日韩| 一区二区三区四区激情视频 | 精品久久久久久久人妻蜜臀av| 免费电影在线观看免费观看| 国产精品久久电影中文字幕| 精品久久久久久,| 熟女人妻精品中文字幕| 永久网站在线| 高清毛片免费观看视频网站| 婷婷亚洲欧美| 欧美日韩乱码在线| 91九色精品人成在线观看| 亚洲av美国av| 日本在线视频免费播放| 午夜福利18| 非洲黑人性xxxx精品又粗又长| 国产爱豆传媒在线观看| 久久九九热精品免费| 国产老妇女一区| 欧美黑人巨大hd| 欧美绝顶高潮抽搐喷水| 免费看美女性在线毛片视频| 国产精品人妻久久久久久| АⅤ资源中文在线天堂| 日本一本二区三区精品| 久久久成人免费电影| 亚洲人成网站在线播| 一级a爱片免费观看的视频| 女人被狂操c到高潮| 国产一区二区在线av高清观看| av视频在线观看入口| 日本熟妇午夜| 国产欧美日韩一区二区精品| 三级男女做爰猛烈吃奶摸视频| 国产精品自产拍在线观看55亚洲| 国产av麻豆久久久久久久| 精品午夜福利在线看| 男女那种视频在线观看| 久久久久久久久久黄片| 男女那种视频在线观看| 国产91精品成人一区二区三区| 别揉我奶头 嗯啊视频| 国产精品av视频在线免费观看| 如何舔出高潮| 亚洲精品成人久久久久久| 如何舔出高潮| 欧美在线黄色| 免费看日本二区| 97人妻精品一区二区三区麻豆| 国产精品伦人一区二区| 亚洲精品一区av在线观看| 亚洲人成伊人成综合网2020| 欧美一区二区精品小视频在线| av女优亚洲男人天堂| 91久久精品电影网| 可以在线观看的亚洲视频| 99久久久亚洲精品蜜臀av| 久久久色成人| 一区二区三区激情视频| 高清毛片免费观看视频网站| 成人特级黄色片久久久久久久| 日韩高清综合在线| 国产精品1区2区在线观看.| 69av精品久久久久久| 国产精品一区二区三区四区免费观看 | av中文乱码字幕在线| 国产成人av教育| 日韩精品青青久久久久久| 美女 人体艺术 gogo| av在线观看视频网站免费| 99热这里只有是精品在线观看 | 久9热在线精品视频| 亚洲成人久久性| 成人av在线播放网站| 国产乱人视频| 无人区码免费观看不卡| 99国产极品粉嫩在线观看| 最新在线观看一区二区三区| 亚洲国产精品999在线| 蜜桃亚洲精品一区二区三区| 丁香欧美五月| 亚洲人成网站在线播| 国产伦精品一区二区三区视频9| 亚洲avbb在线观看| 国产高清激情床上av| 尤物成人国产欧美一区二区三区| 国产精品,欧美在线| 天堂网av新在线| 精品久久久久久久久av| 99久久精品一区二区三区|