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

    變采樣間隔自適應(yīng)EWMAt控制圖統(tǒng)計(jì)設(shè)計(jì)

    2015-12-02 01:27:12常志遠(yuǎn)孫金生
    關(guān)鍵詞:警戒線鏈長(zhǎng)間隔

    常志遠(yuǎn),孫金生

    (南京理工大學(xué) 自動(dòng)化學(xué)院,江蘇 南京 210094)

    0 引言

    1924年,Shewhart博士設(shè)計(jì)了均值控制圖,作為一種有效的統(tǒng)計(jì)質(zhì)量控制方法,控制圖在科學(xué)研究與工業(yè)生產(chǎn)中得到了廣泛的應(yīng)用。隨著科技的發(fā)展,生產(chǎn)過程越來越復(fù)雜,早期的控制圖漸漸不能滿足對(duì)生產(chǎn)過程監(jiān)控的需求。因此,學(xué)者們也在不斷地改進(jìn)控制圖的設(shè)計(jì)方法。

    針對(duì)Shewhart控制圖對(duì)小偏移不敏感的問題,1959年Roberts[1]提出指數(shù)加權(quán)移動(dòng)平均(Exponentially Weighted Moving Average,EWMA)控制圖,這種控制圖具有設(shè)計(jì)簡(jiǎn)單、對(duì)小偏移敏感的特性。EWMA 控制圖提出后受到了廣泛的關(guān)注,學(xué)者們對(duì)此進(jìn)行了深入研究。例如:Lucas等[2]深入研究了EWMA 控制圖的特性,給出了計(jì)算EWMA控制圖平均鏈長(zhǎng)ARL的Markov 鏈方法;Capizzi等[3]將自適應(yīng)EWMA 濾波方法應(yīng)用于EWMA 控制圖,提出自適應(yīng)EWMA(Adaptive EWMA,AEW-MA)控制圖,該控制圖方法可以克服常規(guī)EWMA控制圖只對(duì)固定偏移敏感的問題;Woodall等[4]給出了表征控制圖“慣性(inertia)”的量化指標(biāo),并且指出AEWMA 控制圖要優(yōu)于常規(guī)EWMA 控制圖;Saccucci等[5]提出變采樣間隔EWMA(Variable Sampling Intervals EWMA,VSI-EWMA)控制圖,將變采樣間隔技術(shù)引入EWMA 控制圖,提高了EWMA 控制圖的檢測(cè)速度;Lee等[6]深入研究了多變量VSI-EWMA 控制圖,并給出了控制圖最優(yōu)設(shè)計(jì)算法;薛麗[7]給出了同時(shí)檢測(cè)均值和標(biāo)準(zhǔn)差的EWMA 控制圖的經(jīng)濟(jì)設(shè)計(jì)模型,并優(yōu)化了控制圖參數(shù)。此外,胡雪龍等[8]考慮帶有輔助樣本的X(Supplematary Sample X,SSX)控制圖的經(jīng)濟(jì)設(shè)計(jì),通過模式搜索算法對(duì)費(fèi)用函數(shù)進(jìn)行優(yōu)化,得到了SSX 控制圖的最優(yōu)決策變量。

    以上控制圖設(shè)計(jì)方法均需假設(shè)第一階段的參數(shù)估計(jì)是準(zhǔn)確的,但是這種假設(shè)在實(shí)際中很難做到,當(dāng)方差估計(jì)不準(zhǔn)確時(shí),會(huì)造成控制圖的誤報(bào)漏報(bào)率大大增加。針對(duì)該問題,Zhang等[9]提出t控制圖,證明當(dāng)服從標(biāo)準(zhǔn)正態(tài)分布的檢測(cè)變量出現(xiàn)均值漂移時(shí),均值與標(biāo)準(zhǔn)差的比值服從非中心t分布,并指出t控制圖的優(yōu)勢(shì)在于控制限與方差沒有直接關(guān)系,而是樣本大小的函數(shù);Caleno 等給出了短周期下Shewhartt控制圖的經(jīng)濟(jì)設(shè)計(jì)模型[10];Castagliola等給出了短周期下變樣本大小的t控制圖的設(shè)計(jì)方法以及性能分析[11];Kazemzadeh 等研究 了VSIEWMAt控制圖的特性,并給出了VSI-EWMAt控制圖的優(yōu)化設(shè)計(jì)方法[12]。VSI-EWMAt控制圖的缺點(diǎn)在于固定的參數(shù)只對(duì)固定的偏移敏感,小的平滑系數(shù)會(huì)導(dǎo)致控制圖對(duì)大偏移檢測(cè)的嚴(yán)重滯后。

    針對(duì)文獻(xiàn)[12]提出的VSI-EWMAt控制圖存在的問題,本文提出VSI-AEWMAt控制圖。在Capizzi等[3]提出的AEWMA 控制圖的基礎(chǔ)上,本文首先給出了AEWMAt控制圖的設(shè)計(jì)方法。結(jié)合變采樣間隔技術(shù),在AEWMAt控制圖中加入警戒線,從而給出了根據(jù)統(tǒng)計(jì)量的大小調(diào)節(jié)采樣間隔的VSI-AEWMAt控制圖。變采樣間隔提高了AEWMA 控制圖的檢測(cè)速度,t統(tǒng)計(jì)量可以使控制圖對(duì)方差的估計(jì)具有魯棒性。本文還給出了計(jì)算VSI-AEWMAt控制圖平均報(bào)警時(shí)間ATS的Markov鏈方法,在此基礎(chǔ)上提出一種VSI-AEWMAt控制圖的參數(shù)優(yōu)化算法?;谠搩?yōu)化算法,計(jì)算了特定采樣大小下的控制圖最優(yōu)設(shè)計(jì)參數(shù)。最后,以ATS作為指標(biāo),評(píng)價(jià)了VSI-AEWMAt控制圖與VSI-EWMAt控制圖的性能,結(jié)果表明相同情況下VSI-AEWMAt控制圖具有更小的ATS。

    1 AEWMAt控制圖

    假設(shè){Xi,1,Xi,2,…,Xi,n}來自 同一正 態(tài)分布 的總體,其中i=1,2,…表示采樣時(shí)刻,n表示每個(gè)樣本的大小。Xij~N(μ0+aσ0,,i=1,2,…,1≤j≤n,其中μ0 和分別為正態(tài)分布的均值和方差。a=0,b=1時(shí)過程處于受控狀態(tài),其他情況表明過程出現(xiàn)了異常。每個(gè)采樣時(shí)刻樣本均值和標(biāo)準(zhǔn)差Si分別為:

    過程受控時(shí),通過式(1)中的兩項(xiàng)可以構(gòu)造一個(gè)服從t分布的統(tǒng)計(jì)量Ti,

    Zhang[9]證明了當(dāng)a≠0或者b≠1時(shí),Ti服從非中心t分布,

    式中:Yi表示第i時(shí)刻的AEWMA 統(tǒng)計(jì)量,λ為平滑系數(shù),ei=Ti-Yi-1,φ(ei)為估計(jì)問題中的誤差項(xiàng)。令ω(ei)=φ(ei)/ei,UCL和LCL分別為上下控制限。通過整理,式(4)可重寫為

    由式(5)可知,ω(ei)作為等效的平滑系數(shù),可根據(jù)誤差ei的變化而變化,從而實(shí)現(xiàn)EWMA 算法的自適應(yīng)。本文取φ(ei)為如下形式:

    式中γ為一個(gè)閾值,當(dāng)控制圖設(shè)計(jì)完成后γ為常數(shù)。

    統(tǒng)計(jì)量Ti方差的大小只與樣本大小n有關(guān),第一階段的方差估計(jì)準(zhǔn)確與否并不影響受控狀態(tài)下統(tǒng)計(jì)量Ti的分布形式。

    2 VSI-AEWMAt控制圖

    因?yàn)椴捎檬剑?)的自適應(yīng)算法在偏移加大時(shí)可以采用更大的平滑系數(shù),所以AEWMA 控制圖可以提高EWMA 控制圖對(duì)大偏移的檢測(cè)能力。然而,固定采樣間隔(Fixed Sampling Interval,F(xiàn)SI)AEWMAt控制圖在異常發(fā)生時(shí)仍不能及時(shí)報(bào)警,導(dǎo)致生產(chǎn)損失加大,變采樣間隔技術(shù)可以有效解決這一問題。

    2.1 VSI-AEWMAt控制圖構(gòu)造

    變采樣間隔控制圖是在固定采樣間隔控制圖的基礎(chǔ)上加入上下兩條警戒線(warning line),當(dāng)監(jiān)測(cè)變量落在警戒線以內(nèi)時(shí)采用大的采樣間隔,反之采用小的采樣間隔。設(shè)定AEWMA 的上警戒線UWL=ρ和下警戒線LWL=-ρ,則VSI-AEWMAt控制圖的工作原理如下:

    (1)當(dāng)被監(jiān)測(cè)變量落在警戒線之間,即-ρ≤Yi≤ρ時(shí),VSI-AEWMAt控制圖采用大的采樣間隔hl。

    (2)當(dāng)被監(jiān)測(cè)變量落在警戒線與控制限之間,即LCL≤Yi<-ρ或ρ<Yi≤UCL時(shí),采用小的采樣間隔hs。

    (3)當(dāng)被監(jiān)測(cè)變量落在控制限以外,即Yi<LCL或UCL<Yi時(shí),控制圖報(bào)警。

    2.2 VSI-AEWMAt控制圖ATS的計(jì)算

    平均鏈長(zhǎng)ARL是評(píng)價(jià)固定采樣間隔控制圖性能的一項(xiàng)重要指標(biāo),但是不適用于變采樣間隔控制圖,變采樣間隔控制圖需用ATS來評(píng)價(jià)其性能。根據(jù)文獻(xiàn)[13],VSI-AEWMAt控制圖的ATS計(jì)算過程如下:

    將控制限以內(nèi)的區(qū)間[UCL,LCL]等分為m(m為奇數(shù))個(gè)小區(qū)間,每個(gè)小區(qū)間Iu的寬度為D=(UCL-LCL)/m,每個(gè)小區(qū)間的中心為Su=LCL+(u-0.5)D,u=1,2,…,m。控制限以內(nèi)的各個(gè)狀態(tài)為轉(zhuǎn)移態(tài),控制限以外的區(qū)間為吸收態(tài)。當(dāng)統(tǒng)計(jì)量Yi落入某一小區(qū)間Iu內(nèi)時(shí),認(rèn)為Yi處在狀態(tài)Su。因此可以求得狀態(tài)之間的一步轉(zhuǎn)移概率矩陣

    式中:R為一個(gè)m×m的實(shí)值矩陣,U為一個(gè)全1列向量,0為全0行向量。R中的每一個(gè)元ruv代表i時(shí)刻狀態(tài)u轉(zhuǎn)移到狀態(tài)v的概率,由式(4)得

    由式(6)可知φ(e)是可逆的,其逆函數(shù)為

    將式(9)代入式(8),可得

    AEWMAt控制圖的平均鏈長(zhǎng)可以按照下式求?。?/p>

    式中Pm為一個(gè)1×m的行向量,除了第(m+1)/2個(gè)元為1外其他均為0。

    ARL的精度與m的選取有很大關(guān)系,m越大,ARL的精度越高。當(dāng)m達(dá)到一定值時(shí)再增大m對(duì)提高ARL計(jì)算精度的作用不大,本文選擇m=151。

    VSI-AEWMAt控制圖ATS的計(jì)算與平均鏈長(zhǎng)的計(jì)算類似,只需將式(11)的ARL計(jì)算公式稍作修改即可,VSI-AEWMAt控制圖的ATS可以按照下式來求取:

    式中G為一個(gè)m×1的列向量,其元素由采樣間隔hl與hs組成。G中元素gu的取值規(guī)則為,當(dāng)統(tǒng)計(jì)量落入警戒線內(nèi)時(shí)采用大的采樣間隔,反之采用小的采樣間隔。又因?yàn)榻y(tǒng)計(jì)量的狀態(tài)由Su表征,所以有:

    (1)當(dāng)統(tǒng)計(jì)量Su包含在警戒線以內(nèi),即Su∈[UWL,LWL]時(shí),gu=hl。

    (2)當(dāng)統(tǒng)計(jì)量Su在警戒線以外,即Su∈[-∞,LWL]∪[UWL,∞]時(shí),gu=hs。

    3 VSI-AEWMAt控制圖參數(shù)優(yōu)化

    對(duì)于固定采樣間隔AEWMAt控制圖,平均報(bào)警時(shí)間ATSFSI是通過其平均鏈長(zhǎng)ARLFSI與固定采樣間隔hFSI相乘,即

    VSI-AEWMAt控制圖的平均報(bào)警時(shí)間ATSVSI求取可等價(jià)為平均鏈長(zhǎng)ARLVSI與平均采樣間隔E(h)相乘,即

    變采樣間隔控制圖與固定采樣間隔控制圖作比較的原則是其受控狀態(tài)下的平均報(bào)警時(shí)間相等,失控狀態(tài)下的平均報(bào)警時(shí)間越短越好。令式(13)中的hFSI=E(h),則可以在相同穩(wěn)態(tài)平均鏈長(zhǎng)指標(biāo)下比較固定采樣間隔控制圖與變采樣間隔控制圖的性能。因此,本文將VSI-AEWMAt控制圖的優(yōu)化分為兩個(gè)階段:①優(yōu)化AEWMAt控制圖的ARL,使AEWMAt控制圖受控狀態(tài)平均鏈長(zhǎng)ARL0為固定值,失控狀態(tài)下的平均鏈長(zhǎng)ARL1盡量?。虎谇笕‘?dāng)前情況下的最優(yōu)警戒線,使VSIAEWMAt控制圖的ATS達(dá)到最優(yōu)。

    3.1 VSI-AEWMAt控制圖ARL 的優(yōu)化方法

    對(duì)于AEWMA 控制圖的ARL優(yōu)化,Capizzi等[3]給出了一種利用分層、加權(quán)思想將多目標(biāo)優(yōu)化問題轉(zhuǎn)化為兩個(gè)單目標(biāo)優(yōu)化問題的求解算法。Shu[14]則首先將AEWMA 控制圖簡(jiǎn)化為常規(guī)EWMA 控制圖來處理,求出最優(yōu)的控制限UCL與平滑系數(shù)λ,然后尋找最優(yōu)的γ。Shu[14]的優(yōu)化算法相對(duì)簡(jiǎn)單,但是結(jié)果的保守性較大。本文采用Capizzi等[3]設(shè)計(jì)的優(yōu)化算法,算法描述如下:

    (1)選取需要達(dá)到的穩(wěn)態(tài)平均鏈長(zhǎng)ARL0,設(shè)為B;選取小偏移μ1 與大偏移μ2 的具體值。

    (2)針對(duì)大偏移μ2 計(jì)算參數(shù)組合θ1={UCL,λ,γ},使在ARL0=B前提下失控狀態(tài)下的平均鏈長(zhǎng)最小。

    (3)選擇一個(gè)小的正數(shù)α(本文取α=0.05),針對(duì)小偏移μ1 優(yōu)化參數(shù)組合θ*={UCL,λ,γ},使下式成立:

    由此得到的參數(shù)組合θ*={UCL,λ,γ}可以保證AEWMAt控制圖在小偏移情況下的失控狀態(tài)平均鏈長(zhǎng)接近最小值,并且在過程出現(xiàn)大偏移時(shí)控制圖也能具有理想的性能。

    對(duì)于式(15)和式(16)的求解,Capizzi等[3]使用了200 000次迭代的模擬退火算法,將模擬算法的結(jié)果作為Nelder-Mead單純型法的初值來求取最終結(jié)果。這種求解過程需要耗費(fèi)大量的計(jì)算資源和計(jì)算時(shí)間,結(jié)果并非總是最優(yōu)。粒子群優(yōu)化算法(Particle Swarm Optimization,PSO)作為一種高效的全局優(yōu)化算法得到了廣泛的關(guān)注和研究,傳統(tǒng)的PSO 算法性能已經(jīng)得到了很大提高。本文使用Wang等[15]提出的DNSPSO 算法對(duì)式(15)和式(16)進(jìn)行尋優(yōu),該算法在傳統(tǒng)PSO 的基礎(chǔ)上加入相鄰兩代粒子之間交叉的多樣性促進(jìn)機(jī)制(diversity enhancing mechanism)與鄰域搜索策略(neighborhood search strategies),這些改進(jìn)使DNSPSO 算法能夠以更大的概率搜索到全局最優(yōu)解。相比于Capizzi等[3]采用的模擬退火算法與Nelder-Mead單純型法相結(jié)合的算法,采用DNSPSO 對(duì)AEWMA 控制圖進(jìn)行優(yōu)化可以節(jié)省大量計(jì)算資源。

    3.2 VSI-AEWMAt控制圖ATS 的優(yōu)化方法

    在求出VSI-AEWMAt控制圖對(duì)應(yīng)的最優(yōu)設(shè)計(jì)參數(shù)θ*={UCL,λ,γ}的基礎(chǔ)上,對(duì)VSI-AEWMAt控制圖的ATS進(jìn)行優(yōu)化,此時(shí)只需求取最優(yōu)的警戒線ρ即可。

    由于該模型下的ATS優(yōu)化是一個(gè)多目標(biāo)優(yōu)化問題,常規(guī)的ATS優(yōu)化算法不再適用。類似于AEWMA 控制圖ARL優(yōu)化算法,本文提出一種新的ATS優(yōu)化算法。在受控狀態(tài)下平均報(bào)警時(shí)間ATS0盡量接近其最大值A(chǔ)RL0×hl的前提下,要求失控狀態(tài)下的平均報(bào)警時(shí)間盡量小。

    結(jié)合3.1節(jié)的內(nèi)容,可以得到平均報(bào)警時(shí)間的優(yōu)化步驟:

    (1)選取需要達(dá)到的受控狀態(tài)平均鏈長(zhǎng)ARL0,設(shè)為B;選取小偏移μ1 與大偏移μ2 的具體值。

    (2)選擇最大采樣間隔hl和最小采樣間隔hs,根據(jù)(1)設(shè)定的ARL0得到受控狀態(tài)下的最大平均報(bào)警時(shí)間ATS0=B×hl。

    (3)針對(duì)大偏移μ2 優(yōu)化警戒線ρ′,取一個(gè)接近1的正數(shù)η0(本文取η0=0.95),根據(jù)3.1節(jié)的最優(yōu)設(shè)計(jì)參數(shù)θ*={UCL,λ,γ},使得下式成立:

    (4)針對(duì)小偏移μ1 優(yōu)化警戒線ρ*,選取一個(gè)接近0的正數(shù)η1(本文取η1=0.05),使下式成立:

    由此得到最優(yōu)警戒線為ρ*,即為最終結(jié)果。

    這種優(yōu)化方法的優(yōu)點(diǎn)在于將四維優(yōu)化問題做降維處理,得到一個(gè)三維優(yōu)化問題與一個(gè)一維優(yōu)化問題,從而大大降低原優(yōu)化問題的復(fù)雜性。求解式(17)和式(18)優(yōu)化問題的方法與式(15)和式(16)類似,也是通過使用DNSPSO 優(yōu)化方法求得數(shù)值解。

    4 仿真研究

    對(duì)受控狀態(tài)情況下服從標(biāo)準(zhǔn)正態(tài)分布N(0,1)的過程進(jìn)行仿真研究。選取小偏移μ1 分別為0.25,0.5,1,對(duì)應(yīng)的大偏移μ2 分別為4,5,6。首先,計(jì)算樣本大小n=5~10時(shí)VSI-AEWMAt控制圖的最優(yōu)設(shè)計(jì)參數(shù)λ,h,γ及其對(duì)應(yīng)的ARL,結(jié)果如表1所示。

    表1 方差估計(jì)準(zhǔn)確時(shí)VSI-AEWMAt控制圖的最優(yōu)參數(shù)及對(duì)應(yīng)的ARL

    基于表1的結(jié)果,計(jì)算采樣間隔{hl,hs}分別為{2.5,0.2},{2.0,0.2},{1.7,0.3}時(shí)ATS最優(yōu)警戒線的值,如表2所示。

    表2 樣本大小為5,不同采樣采樣間隔的下的警戒線

    由表2可知,警戒線隨著偏移的加大而變大,隨著控制限的增加而增加。此外,可以發(fā)現(xiàn)警戒線的選取與均值偏離實(shí)際值的大小以及穩(wěn)態(tài)ARL的大小有更顯著的關(guān)系,與采樣間隔的關(guān)系則不是很明顯。警戒線表征了質(zhì)量特性發(fā)生偏移的概率,因此警戒線的取值將受被監(jiān)測(cè)變量的分布形式及分布參數(shù)的影響。本文還計(jì)算了樣本大小n=5、采樣間隔{hl,hs}為{2.0,0.2}和{1.7,0.3}時(shí)的ATS,結(jié)果如表3和表4所示。

    表3 樣本大小為5,采樣間隔hl=2.0,hs=0.2時(shí)的ATS

    表4 樣本大小為5,采樣間隔hl=1.7,hs=0.3時(shí)的ATS

    由表3和表4可知,受控狀態(tài)下ATS達(dá)到了最大化的要求,失控狀態(tài)下的ATS隨著偏移的增大迅速減小。選取文獻(xiàn)[9]中w=0.3對(duì)應(yīng)的情況,與本文偏移{μ1,μ2}={0.5,5}對(duì)應(yīng)的情況進(jìn)行比較,以失控狀態(tài)下的平均報(bào)警時(shí)間ATS1與受控狀態(tài)下的平均報(bào)警時(shí)間ATS0的比值K=ATS1/ATS0作為目標(biāo),比較結(jié)果如圖1所示。

    由圖1可知,VSI-AEWMAt控制圖的性能優(yōu)于針對(duì)不同偏移分別優(yōu)化的VSI-EWMAt控制圖的性能。選取文獻(xiàn)[9]中{hl,hs}={1.7,0.3},w=0.3時(shí)針對(duì)偏移為0.1 與2.0 設(shè)計(jì)的控制圖,與本文{hl,hs}={1.7,0.3},ATS0=599 時(shí),偏 移{μ1,μ2}={0.25,4}的情況進(jìn)行比較,仍然選用失控狀態(tài)下平均報(bào)警時(shí)間與受控狀態(tài)下平均報(bào)警時(shí)間的比值作為目標(biāo),結(jié)果如圖2所示。

    由圖2可知,當(dāng)VSI-EWMAt控制圖參數(shù)針對(duì)固定偏移設(shè)計(jì)時(shí),會(huì)出現(xiàn)對(duì)其他偏移不敏感的問題,實(shí)際偏移與控制圖設(shè)計(jì)的偏移距離越遠(yuǎn),這種不敏感性越明顯。而VSI-AEWMAt控制圖可以較好地避免這種情況。因此在實(shí)際應(yīng)用中VSI-AEWMAt控制圖比VSI-EWMAt控制圖更加簡(jiǎn)便有效。

    t控制圖對(duì)方差估計(jì)具有魯棒性,因此本文研究了b=0.8,0.9,1.0,1.1,1.2 五種 情況下VSIAEWMAt控制圖的ATS隨a變化的情況。在樣本大小n=5,hl=2.0,hs=0.2,μ1=1,μ2=6,ARL=370的情況下進(jìn)行仿真分析,結(jié)果如圖3所示。

    從圖3可以看出,方差估計(jì)不準(zhǔn)確會(huì)在a≠0時(shí)對(duì)t控制圖的ATS帶來影響,但不會(huì)影響到穩(wěn)態(tài)ATS,體現(xiàn)了t控制圖具有對(duì)方差估計(jì)魯棒性的特點(diǎn)。當(dāng)方差估計(jì)偏小時(shí),會(huì)使失控狀態(tài)下的ATS減??;當(dāng)方差估計(jì)偏大時(shí),會(huì)使失控狀態(tài)下的ATS增大。由非中心參數(shù)的表達(dá)式可知,這種變化趨勢(shì)是正確的。

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

    固定平滑系數(shù)的EWMAt控制圖只對(duì)優(yōu)化設(shè)計(jì)時(shí)所選用的偏移敏感,且異常發(fā)生時(shí)固定采樣間隔不能及時(shí)報(bào)警。本文針對(duì)以上問題構(gòu)造了VSIAEWMAt控制圖,并給出了VSI-AEWMAt控制圖ATS計(jì)算的Markov鏈方法;運(yùn)用多目標(biāo)優(yōu)化理論,優(yōu)化了VSI-AEWMAt控制圖的決策參數(shù);以ATS作為評(píng)價(jià)指標(biāo),比較了VSI-AEWMAt控制圖與VSI-EWMAt控制圖。結(jié)果表明,VSI-AEWMAt控制圖可以在不損失小偏移檢測(cè)能力的前提下提高對(duì)大偏移的檢測(cè)能力,且VSI特性明顯縮短了失控狀態(tài)下的平均報(bào)警時(shí)間。

    同時(shí),VSI-AEWMAt控制圖的優(yōu)化過程比較復(fù)雜,如何提高優(yōu)化算法的有效性與減小結(jié)果的保守性值得進(jìn)一步研究。另外,VSI-AEWMAt控制圖的經(jīng)濟(jì)設(shè)計(jì)也是值得下一步研究的課題。

    [1]ROBERTS S W.Control chart tests based on geometric moving averages[J].Technometrics,1959,3(1):239-250.

    [2]LUCAS J M,SACCUCCI M S.Exponentially weighted moving average control schemes:properties and enhancements[J].Technometrics,1990,32(1):1-12.

    [3]CAPIZZI G,MASAROTTO G.An adaptive exponentially weighted moving average control chart[J].Technometics,2003,45(3):199-207.

    [4]WOODALL W H,MAHMOUND A M.The inertial properties of quality control charts[J].Technometric,2005,47(4):425-436.

    [5]SACCUCCI M S,AMIN R W,LUCAS J M.Exponentially weighted moving average control schemes with variable sampling intervals[J].Communications in Statistics-Simulation and Computation,1992,21(3):627-657.

    [6]LEE M H,KHOO M B C.Design of a multivariable exponentially weighted moving average control chart with variable sampling intervals[J].Computational Statistics,2014,29(1/2):189-214.

    [7]XUE Li.Economic design of variable sampling intervals EWMA chart for monitoring process mean and standard deviation[J].Computer Integrated Manufacturing Systems,2013,19(6):1369-1376(in Chinese).[薛 麗.監(jiān)控均值標(biāo)準(zhǔn)差的可變抽樣區(qū)間EWMA 圖經(jīng)濟(jì)設(shè)計(jì)[J].計(jì)算機(jī)集成制造系統(tǒng),2013,19(6):1369-1376.]

    [8]HU Xuelong,SUN Jinsheng.Economic design of X control chart with supplementary samples[J].Computer Integrated Manufacturing Systems,2015,21(1):160-168(in Chinese).[胡雪龍,孫金生.帶有輔助樣本的X 控制圖經(jīng)濟(jì)設(shè)計(jì)[J].計(jì)算機(jī)集成制造系統(tǒng),2015,21(1):160-168.]

    [9]ZHANG L,CHEN G,CASTAGLIOLA P.On t and EWMA t charts for monitoring changes in the process mean[J].Quality and Reliability Engineering International,2009,25(8):933-945.

    [10]CELANO G,CASTAGLIOLA P,TROVATO E,et al.The economic performance of the shewhart t chart[J].Quality and Reliability Engineering International,2012,28(2):159-180.

    [11]CASTAGLIOLA P,CELANO G,F(xiàn)ICHERA S,et al.The variable sample size t control chart for monitoring short production runs[J].The International Journal of Advanced Manufacturing Technology,2013,66(9-12):1353-1366.

    [12]KAZEMZADEH R B,KARBASIAN M,BABAKHANI M A.An EWMA t chart with variable sampling intervals for monitoring the process mean[J].The International Journal of Advanced Manufacturing Technology,2013,66(1-4):125-139.

    [13]REYNOLDS M R Jr,AMIN R W,AMOLD J C.CUSUM chart with variable sampling intervals[J].Techmometrics,1990,32(4):371-384.

    [14]SHU L.An adaptive exponentially weighted moving average control chart for monitoring process variances[J].Journal of Statistical Computation and Simulation,2008,78(4):367-384.

    [15]WANG H,SUN H,LI C H,et al.Diversity enhanced particle swarm optimization with neighborhood search[J].Information Sciences,2013,223:119-135.

    猜你喜歡
    警戒線鏈長(zhǎng)間隔
    黨建賦能“鏈長(zhǎng)制”落實(shí)落地
    “警戒線”
    安全藝苑
    中泰紡織集團(tuán):做最強(qiáng)“鏈長(zhǎng)”,引領(lǐng)新疆紡織邁向新高度
    間隔問題
    健康月歷
    健康之家(2020年7期)2020-11-02 02:18:28
    間隔之謎
    休哈特控制圖的改進(jìn)
    烷基鏈長(zhǎng)及肽鏈電荷分布對(duì)脂肽雙親分子自組裝及水凝膠化的影響
    上樓梯的學(xué)問
    www.自偷自拍.com| 亚洲五月色婷婷综合| 中文字幕高清在线视频| 亚洲成国产人片在线观看| 欧美在线一区亚洲| 91麻豆精品激情在线观看国产 | 久久中文字幕人妻熟女| 日本欧美视频一区| 亚洲熟女毛片儿| 热re99久久国产66热| 黄色丝袜av网址大全| 99在线视频只有这里精品首页| 男女之事视频高清在线观看| 亚洲狠狠婷婷综合久久图片| 亚洲片人在线观看| 成人国产一区最新在线观看| 色婷婷久久久亚洲欧美| 亚洲中文日韩欧美视频| 精品欧美一区二区三区在线| 亚洲一卡2卡3卡4卡5卡精品中文| 免费女性裸体啪啪无遮挡网站| 1024香蕉在线观看| 涩涩av久久男人的天堂| 亚洲欧美激情在线| 国产亚洲精品综合一区在线观看 | 中文字幕人妻熟女乱码| 这个男人来自地球电影免费观看| 嫩草影院精品99| 国产一区在线观看成人免费| 99久久久亚洲精品蜜臀av| 久久久久国产精品人妻aⅴ院| 国产成人av激情在线播放| 一进一出抽搐动态| 中国美女看黄片| 欧美激情高清一区二区三区| 19禁男女啪啪无遮挡网站| 精品国产美女av久久久久小说| 久久国产乱子伦精品免费另类| 国产成人欧美| 一级a爱视频在线免费观看| 亚洲专区字幕在线| 亚洲欧美一区二区三区久久| 丁香六月欧美| 黑人欧美特级aaaaaa片| 欧美中文综合在线视频| 热99re8久久精品国产| 在线观看免费日韩欧美大片| 黄片大片在线免费观看| 人成视频在线观看免费观看| 一夜夜www| 国产一区二区三区在线臀色熟女 | 欧美日韩福利视频一区二区| 在线看a的网站| 大陆偷拍与自拍| 亚洲熟妇熟女久久| 精品乱码久久久久久99久播| 亚洲精品一二三| 日本 av在线| 天天躁狠狠躁夜夜躁狠狠躁| 老司机靠b影院| 热re99久久精品国产66热6| 亚洲av五月六月丁香网| 欧美激情久久久久久爽电影 | 91老司机精品| 欧美中文综合在线视频| 99国产精品一区二区蜜桃av| 日本a在线网址| 欧美日韩乱码在线| 国产成人精品在线电影| 国产高清激情床上av| 久久国产精品男人的天堂亚洲| 国产精品影院久久| 香蕉国产在线看| 国产免费av片在线观看野外av| 午夜精品国产一区二区电影| 精品国产亚洲在线| 亚洲欧美日韩另类电影网站| 日韩国内少妇激情av| 久久精品国产综合久久久| 黄色毛片三级朝国网站| 正在播放国产对白刺激| 欧美另类亚洲清纯唯美| 人妻久久中文字幕网| 首页视频小说图片口味搜索| 很黄的视频免费| 国产熟女午夜一区二区三区| 丝袜人妻中文字幕| 91国产中文字幕| 久久精品91无色码中文字幕| 母亲3免费完整高清在线观看| 国产一区二区激情短视频| 在线天堂中文资源库| 两个人免费观看高清视频| 在线观看免费视频网站a站| 午夜影院日韩av| 亚洲伊人色综图| 精品国产超薄肉色丝袜足j| 天堂俺去俺来也www色官网| 999精品在线视频| 黄色女人牲交| 国产xxxxx性猛交| 国产亚洲av高清不卡| 水蜜桃什么品种好| 午夜91福利影院| netflix在线观看网站| 亚洲av片天天在线观看| 一级毛片女人18水好多| 久久 成人 亚洲| 国产1区2区3区精品| 韩国av一区二区三区四区| 亚洲国产欧美网| 国产精品电影一区二区三区| 精品第一国产精品| 国内久久婷婷六月综合欲色啪| 很黄的视频免费| 午夜精品久久久久久毛片777| 一区在线观看完整版| 国产黄色免费在线视频| av超薄肉色丝袜交足视频| 亚洲一码二码三码区别大吗| 欧美在线一区亚洲| 51午夜福利影视在线观看| 俄罗斯特黄特色一大片| 日韩欧美免费精品| 大型av网站在线播放| 夜夜爽天天搞| 99国产综合亚洲精品| 极品教师在线免费播放| 午夜亚洲福利在线播放| 国产一区二区激情短视频| 久久精品人人爽人人爽视色| 国产亚洲欧美在线一区二区| 宅男免费午夜| 国产精品免费一区二区三区在线| 国产精品免费一区二区三区在线| 国产欧美日韩精品亚洲av| av有码第一页| 51午夜福利影视在线观看| 最近最新免费中文字幕在线| 久久天躁狠狠躁夜夜2o2o| 久久久久国内视频| 国产精品 国内视频| 中亚洲国语对白在线视频| 免费在线观看黄色视频的| 91国产中文字幕| 精品国产一区二区久久| 夜夜夜夜夜久久久久| 黑人猛操日本美女一级片| 两个人看的免费小视频| 啦啦啦在线免费观看视频4| 黄色毛片三级朝国网站| 好看av亚洲va欧美ⅴa在| 看片在线看免费视频| 亚洲欧美一区二区三区久久| 看黄色毛片网站| 手机成人av网站| 国产精品综合久久久久久久免费 | 久久人人爽av亚洲精品天堂| 精品国产乱子伦一区二区三区| 色综合婷婷激情| 国产精品乱码一区二三区的特点 | 中文字幕av电影在线播放| 精品高清国产在线一区| 国产亚洲精品一区二区www| 激情在线观看视频在线高清| 少妇 在线观看| av天堂在线播放| 精品一区二区三卡| 午夜免费观看网址| 中文字幕人妻丝袜一区二区| 夜夜躁狠狠躁天天躁| 脱女人内裤的视频| 黑人操中国人逼视频| 色老头精品视频在线观看| 妹子高潮喷水视频| 亚洲专区中文字幕在线| 神马国产精品三级电影在线观看 | 悠悠久久av| 人人妻,人人澡人人爽秒播| 精品一区二区三卡| 岛国视频午夜一区免费看| 国产蜜桃级精品一区二区三区| 国产男靠女视频免费网站| 两个人免费观看高清视频| 国产亚洲欧美在线一区二区| 免费女性裸体啪啪无遮挡网站| 又大又爽又粗| 久久久国产成人免费| 一级作爱视频免费观看| 精品欧美一区二区三区在线| 黑人欧美特级aaaaaa片| 9191精品国产免费久久| 在线观看www视频免费| 男女午夜视频在线观看| 亚洲视频免费观看视频| 极品教师在线免费播放| 俄罗斯特黄特色一大片| 久久99一区二区三区| 国产精品 国内视频| 91国产中文字幕| 18禁黄网站禁片午夜丰满| 欧美成人免费av一区二区三区| 久久人人爽av亚洲精品天堂| av电影中文网址| 亚洲专区字幕在线| 久热爱精品视频在线9| 亚洲成人久久性| 最好的美女福利视频网| 国产精品国产av在线观看| 成人免费观看视频高清| 亚洲一码二码三码区别大吗| 成人黄色视频免费在线看| av片东京热男人的天堂| 国内久久婷婷六月综合欲色啪| 首页视频小说图片口味搜索| 淫妇啪啪啪对白视频| 国产免费男女视频| 国产高清视频在线播放一区| 美女午夜性视频免费| 长腿黑丝高跟| 夜夜爽天天搞| 美女 人体艺术 gogo| 亚洲avbb在线观看| 麻豆一二三区av精品| 国产精品免费视频内射| 99久久人妻综合| 国产精品二区激情视频| 高清黄色对白视频在线免费看| 亚洲一区二区三区欧美精品| 久久草成人影院| 成人18禁在线播放| 日日夜夜操网爽| 国产熟女xx| 精品少妇一区二区三区视频日本电影| 国产欧美日韩一区二区三区在线| 男人舔女人下体高潮全视频| 十八禁网站免费在线| 一夜夜www| 午夜老司机福利片| 国产成人啪精品午夜网站| 中文字幕高清在线视频| 午夜两性在线视频| 午夜影院日韩av| 日韩三级视频一区二区三区| av网站免费在线观看视频| 欧美日韩福利视频一区二区| 国产精品一区二区在线不卡| 两性午夜刺激爽爽歪歪视频在线观看 | 女生性感内裤真人,穿戴方法视频| 黄频高清免费视频| 午夜免费激情av| 桃色一区二区三区在线观看| 中文字幕av电影在线播放| 久热爱精品视频在线9| 亚洲少妇的诱惑av| 一二三四社区在线视频社区8| 精品免费久久久久久久清纯| 一级作爱视频免费观看| www.熟女人妻精品国产| 无限看片的www在线观看| 亚洲精品中文字幕在线视频| 国产又爽黄色视频| 在线天堂中文资源库| 在线观看午夜福利视频| 亚洲免费av在线视频| 婷婷丁香在线五月| 人妻丰满熟妇av一区二区三区| www日本在线高清视频| 嫩草影视91久久| 男女做爰动态图高潮gif福利片 | 老汉色∧v一级毛片| 精品国产美女av久久久久小说| 丝袜美腿诱惑在线| 女警被强在线播放| 窝窝影院91人妻| 亚洲色图av天堂| 99国产精品一区二区三区| 999精品在线视频| 日韩欧美一区视频在线观看| 国产伦人伦偷精品视频| aaaaa片日本免费| 窝窝影院91人妻| 亚洲午夜精品一区,二区,三区| 99久久人妻综合| 一二三四社区在线视频社区8| 久久久久久久久久久久大奶| 狂野欧美激情性xxxx| 免费一级毛片在线播放高清视频 | 亚洲自偷自拍图片 自拍| 97碰自拍视频| 一二三四社区在线视频社区8| 久久久久久久久久久久大奶| 91麻豆av在线| 日韩免费高清中文字幕av| 久久久久亚洲av毛片大全| 亚洲熟妇中文字幕五十中出 | 国产熟女xx| 午夜精品国产一区二区电影| 亚洲av美国av| 日韩人妻精品一区2区三区| 欧美日韩亚洲国产一区二区在线观看| 久久精品亚洲精品国产色婷小说| 91成人精品电影| 国产一区二区激情短视频| 久久久久九九精品影院| 伦理电影免费视频| 欧美精品亚洲一区二区| 在线观看一区二区三区激情| av有码第一页| 国产成人一区二区三区免费视频网站| 91大片在线观看| 侵犯人妻中文字幕一二三四区| 久久国产精品男人的天堂亚洲| 少妇裸体淫交视频免费看高清 | 亚洲人成77777在线视频| 女警被强在线播放| 亚洲av片天天在线观看| 国产欧美日韩一区二区三区在线| 亚洲精品美女久久久久99蜜臀| 久久精品91无色码中文字幕| 黄色丝袜av网址大全| 亚洲欧美日韩另类电影网站| 男女做爰动态图高潮gif福利片 | 亚洲国产精品sss在线观看 | 极品人妻少妇av视频| 精品久久久精品久久久| 嫩草影视91久久| 国产成人精品在线电影| 日本vs欧美在线观看视频| 99久久99久久久精品蜜桃| 亚洲国产欧美日韩在线播放| 亚洲熟妇中文字幕五十中出 | 欧美精品亚洲一区二区| 午夜福利一区二区在线看| 国产午夜精品久久久久久| 69精品国产乱码久久久| 欧美久久黑人一区二区| 亚洲成人久久性| 成年人免费黄色播放视频| 欧美日韩一级在线毛片| 99精品在免费线老司机午夜| 亚洲成人久久性| 美女高潮到喷水免费观看| 国产精品美女特级片免费视频播放器 | 成人特级黄色片久久久久久久| 男女床上黄色一级片免费看| 成年人黄色毛片网站| 正在播放国产对白刺激| 亚洲 欧美 日韩 在线 免费| 亚洲av成人av| 99久久精品国产亚洲精品| 香蕉久久夜色| 午夜免费成人在线视频| 91av网站免费观看| 人妻久久中文字幕网| 亚洲一区二区三区不卡视频| 亚洲国产精品sss在线观看 | 亚洲一区高清亚洲精品| 国产精品乱码一区二三区的特点 | svipshipincom国产片| 免费久久久久久久精品成人欧美视频| 国产免费男女视频| 久久人人精品亚洲av| 啪啪无遮挡十八禁网站| 国产人伦9x9x在线观看| 丝袜在线中文字幕| 亚洲人成77777在线视频| 精品午夜福利视频在线观看一区| 亚洲色图av天堂| 美女高潮喷水抽搐中文字幕| 中文字幕人妻熟女乱码| 女人爽到高潮嗷嗷叫在线视频| 亚洲少妇的诱惑av| 国产亚洲精品久久久久5区| 精品午夜福利视频在线观看一区| 好看av亚洲va欧美ⅴa在| 欧美日本亚洲视频在线播放| 久久久久久久精品吃奶| 久久人妻av系列| 亚洲 国产 在线| 国产精品久久久久久人妻精品电影| 久久影院123| 久久婷婷成人综合色麻豆| 极品教师在线免费播放| 精品日产1卡2卡| 一级a爱视频在线免费观看| 欧美激情 高清一区二区三区| 女警被强在线播放| 老司机深夜福利视频在线观看| 精品国产一区二区三区四区第35| 久久久久九九精品影院| 一级毛片高清免费大全| 午夜影院日韩av| 久久精品91无色码中文字幕| 日韩大尺度精品在线看网址 | 中文亚洲av片在线观看爽| 成年人免费黄色播放视频| 最新在线观看一区二区三区| 亚洲专区国产一区二区| 亚洲精品美女久久久久99蜜臀| 日韩一卡2卡3卡4卡2021年| 国产精品99久久99久久久不卡| 亚洲男人天堂网一区| 久久这里只有精品19| 757午夜福利合集在线观看| 宅男免费午夜| 日韩有码中文字幕| av电影中文网址| 丰满的人妻完整版| 精品熟女少妇八av免费久了| 高清欧美精品videossex| 国产精品国产av在线观看| 如日韩欧美国产精品一区二区三区| 91精品国产国语对白视频| 新久久久久国产一级毛片| 久久久久久久久久久久大奶| 精品一区二区三区av网在线观看| 超碰97精品在线观看| videosex国产| 精品日产1卡2卡| 亚洲精品中文字幕在线视频| 香蕉丝袜av| 国产高清国产精品国产三级| 久久人人97超碰香蕉20202| 啪啪无遮挡十八禁网站| 淫妇啪啪啪对白视频| 欧美日本亚洲视频在线播放| 亚洲一码二码三码区别大吗| 国产日韩一区二区三区精品不卡| 久久久久亚洲av毛片大全| 色精品久久人妻99蜜桃| 亚洲 欧美一区二区三区| 久久久国产欧美日韩av| a在线观看视频网站| 日日爽夜夜爽网站| 亚洲第一av免费看| 青草久久国产| 亚洲午夜理论影院| 久久亚洲精品不卡| 中文字幕人妻熟女乱码| 99精国产麻豆久久婷婷| 成熟少妇高潮喷水视频| 麻豆久久精品国产亚洲av | 91大片在线观看| 最好的美女福利视频网| 十八禁人妻一区二区| 一区二区三区国产精品乱码| 国产不卡一卡二| 久久久久九九精品影院| 91在线观看av| 亚洲视频免费观看视频| 日本一区二区免费在线视频| 欧美 亚洲 国产 日韩一| 久热这里只有精品99| 国产av又大| 热99国产精品久久久久久7| 啦啦啦 在线观看视频| svipshipincom国产片| 91成人精品电影| 欧美人与性动交α欧美精品济南到| 久久久久国产精品人妻aⅴ院| 亚洲av成人av| √禁漫天堂资源中文www| 午夜精品国产一区二区电影| 亚洲精华国产精华精| 精品国产国语对白av| 精品乱码久久久久久99久播| 黑人猛操日本美女一级片| 亚洲片人在线观看| 精品久久久久久成人av| 亚洲精品一区av在线观看| 欧美精品一区二区免费开放| 久久久久国产精品人妻aⅴ院| 日韩三级视频一区二区三区| 亚洲男人的天堂狠狠| 好男人电影高清在线观看| 免费看十八禁软件| 日韩免费高清中文字幕av| 一二三四社区在线视频社区8| 精品日产1卡2卡| 97碰自拍视频| 超色免费av| 岛国在线观看网站| 国产一区二区三区在线臀色熟女 | 男人舔女人下体高潮全视频| 日韩一卡2卡3卡4卡2021年| 久久精品91蜜桃| 久久99一区二区三区| 成人精品一区二区免费| 91麻豆精品激情在线观看国产 | 黄网站色视频无遮挡免费观看| 免费观看人在逋| 国产精品野战在线观看 | 精品日产1卡2卡| 国产精品综合久久久久久久免费 | 精品国产美女av久久久久小说| 欧美黄色淫秽网站| 亚洲专区中文字幕在线| 91国产中文字幕| 手机成人av网站| 日韩欧美免费精品| 精品国内亚洲2022精品成人| 免费少妇av软件| 一进一出抽搐动态| 午夜两性在线视频| 最近最新中文字幕大全免费视频| 热re99久久国产66热| 亚洲精品一二三| 欧美日韩亚洲国产一区二区在线观看| 国产三级在线视频| 人人妻人人澡人人看| 国产99白浆流出| 18美女黄网站色大片免费观看| 成人永久免费在线观看视频| av有码第一页| 国产精品一区二区免费欧美| 亚洲片人在线观看| 99久久综合精品五月天人人| 精品人妻1区二区| 日韩有码中文字幕| 不卡一级毛片| 久久久久国内视频| 亚洲一区中文字幕在线| 99国产综合亚洲精品| 久久久久久久久中文| 亚洲熟妇熟女久久| 亚洲情色 制服丝袜| 丰满的人妻完整版| av网站在线播放免费| 日本五十路高清| 国产av又大| 亚洲熟女毛片儿| 中亚洲国语对白在线视频| 757午夜福利合集在线观看| 丝袜美腿诱惑在线| 亚洲av成人一区二区三| 亚洲精品国产精品久久久不卡| 人人澡人人妻人| 另类亚洲欧美激情| 免费久久久久久久精品成人欧美视频| 亚洲欧美日韩无卡精品| 亚洲成a人片在线一区二区| www.自偷自拍.com| 女警被强在线播放| 亚洲精品中文字幕在线视频| 这个男人来自地球电影免费观看| 成年版毛片免费区| 亚洲av第一区精品v没综合| av电影中文网址| 久久中文字幕人妻熟女| 91在线观看av| 亚洲人成77777在线视频| 看片在线看免费视频| 欧美日本中文国产一区发布| 精品久久久久久久久久免费视频 | 久久精品国产亚洲av高清一级| 亚洲av成人一区二区三| 免费高清视频大片| 久久人妻av系列| 免费在线观看完整版高清| 中文字幕最新亚洲高清| 久热这里只有精品99| 成人三级做爰电影| 亚洲第一av免费看| 亚洲成国产人片在线观看| 无遮挡黄片免费观看| 欧美丝袜亚洲另类 | 午夜福利一区二区在线看| 精品午夜福利视频在线观看一区| 欧美日本亚洲视频在线播放| 桃色一区二区三区在线观看| 亚洲av成人不卡在线观看播放网| 操出白浆在线播放| ponron亚洲| 99在线视频只有这里精品首页| 麻豆久久精品国产亚洲av | 欧美精品亚洲一区二区| av网站免费在线观看视频| 激情在线观看视频在线高清| 亚洲精品中文字幕一二三四区| 国产精品久久久久成人av| av在线天堂中文字幕 | 欧美一区二区精品小视频在线| 日韩有码中文字幕| 亚洲精品一二三| 美女 人体艺术 gogo| 女人高潮潮喷娇喘18禁视频| 精品无人区乱码1区二区| 久久欧美精品欧美久久欧美| 欧美成狂野欧美在线观看| 激情视频va一区二区三区| 欧美乱妇无乱码| 妹子高潮喷水视频| 男人操女人黄网站| 久久婷婷成人综合色麻豆| 久久香蕉激情| 天天躁夜夜躁狠狠躁躁| 热99re8久久精品国产| 国产精品影院久久| 成人国语在线视频| 国产又爽黄色视频| www.精华液| 欧美激情 高清一区二区三区| 淫妇啪啪啪对白视频| 88av欧美| 国产1区2区3区精品| 国产免费现黄频在线看| 欧美乱码精品一区二区三区| 中文字幕另类日韩欧美亚洲嫩草| 久久久久国产精品人妻aⅴ院| 日韩人妻精品一区2区三区| 亚洲 欧美 日韩 在线 免费| 国产熟女午夜一区二区三区| 午夜精品国产一区二区电影| 免费不卡黄色视频|