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

    基于加權(quán)滑動(dòng)平均的磁盤使用率預(yù)測(cè)模型

    2016-02-23 07:19:08張宗華張海全李師航牛新征
    關(guān)鍵詞:計(jì)算技術(shù)

    張宗華,張海全, 李師航,牛新征

    1)南京南瑞集團(tuán)公司流程與信息管理中心,江蘇南京 211106;2)西南財(cái)經(jīng)大學(xué)經(jīng)濟(jì)信息工程學(xué)院,四川成都 611130;

    3)電子科技大學(xué)計(jì)算機(jī)科學(xué)與工程學(xué)院,四川成都 611731

    ?

    基于加權(quán)滑動(dòng)平均的磁盤使用率預(yù)測(cè)模型

    張宗華1,張海全1, 李師航2,牛新征3

    1)南京南瑞集團(tuán)公司流程與信息管理中心,江蘇南京 211106;2)西南財(cái)經(jīng)大學(xué)經(jīng)濟(jì)信息工程學(xué)院,四川成都 611130;

    3)電子科技大學(xué)計(jì)算機(jī)科學(xué)與工程學(xué)院,四川成都 611731

    摘要:為能提前做好擴(kuò)容準(zhǔn)備,提出一種改進(jìn)的加權(quán)滑動(dòng)平均(weighted moving average, WMA)模型,用以預(yù)測(cè)未來短期內(nèi)磁盤的使用率. 針對(duì)磁盤使用率序列變化較為平緩、要求滯后較小的特性,采用自相關(guān)和偏自相關(guān)系數(shù)法對(duì)模型定階,處理數(shù)據(jù)后,在不影響精度的前提下計(jì)算最小滯后值,并使用結(jié)合了拉依達(dá)準(zhǔn)則的權(quán)重轉(zhuǎn)移法來均衡權(quán)重,用多新息遞推最小二乘法對(duì)參數(shù)進(jìn)行更精確的估計(jì),以提高預(yù)測(cè)的準(zhǔn)確性. 通過Matlab仿真實(shí)驗(yàn)可知,該算法預(yù)測(cè)誤差小,滯后性弱,與原始WMA模型相比,具有更好的預(yù)測(cè)效果.

    關(guān)鍵詞:計(jì)算技術(shù);加權(quán)滑動(dòng)平均模型;磁盤使用率;自相關(guān)和偏自相關(guān)系數(shù)法;拉依達(dá)準(zhǔn)則;權(quán)重轉(zhuǎn)移;多新息遞推最小二乘法

    Disk usage prediction based on an improved

    weighted moving average method

    Zhang Zonghua1?, Zhang Haiquan1, Li Shihang2, and Niu Xinzheng3

    磁盤使用率的預(yù)測(cè)主要用于資源管理、故障管理、實(shí)時(shí)監(jiān)控以及對(duì)異常情況進(jìn)行警報(bào),防止磁盤寫滿造成數(shù)據(jù)丟失等問題. 目前有關(guān)磁盤使用率預(yù)測(cè)的研究不多.Murray等[1]提出了一種融合線性和指數(shù)回歸模型、自回歸積分滑動(dòng)平均模型及貝葉斯結(jié)構(gòu)時(shí)間序列模型的算法,而在磁盤使用率短期預(yù)測(cè)方面目前尚欠有效算法.

    短期預(yù)測(cè)要求精度較高,延遲較小,比較適合用時(shí)間序列模型[2]進(jìn)行預(yù)測(cè). 常見的時(shí)間序列預(yù)測(cè)模型包括用來對(duì)各種自然現(xiàn)象進(jìn)行預(yù)測(cè)的自回歸(auto regressive, AR)模型、通過一組時(shí)間序列逐次移動(dòng)完成計(jì)算的滑動(dòng)平均(moving average, MA)模型、由自回歸模型和滑動(dòng)平均模型迭加形成的自回歸滑動(dòng)平均(auto regressive moving average, ARMA)模型及經(jīng)過差分處理后轉(zhuǎn)化為ARMA并與之形式類似的自回歸積分滑動(dòng)平均模型(auto regressive integrated moving average, ARIMA). 磁盤使用率數(shù)據(jù)沒有明顯趨勢(shì)性,變化較平穩(wěn),而在時(shí)間序列模型中,MA模型最適合處理平穩(wěn)數(shù)據(jù),且運(yùn)行效率高,符合預(yù)測(cè)和警告的即時(shí)性要求,已在諸如股票交易、降雨風(fēng)速和金融時(shí)間序列等方面有大量應(yīng)用[3-4],取得了較好的預(yù)測(cè)效果. MA模型基本類型可以分為簡(jiǎn)單滑動(dòng)平均(simple moving average, SMA)模型和加權(quán)滑動(dòng)平均(weighting moving average, WMA)模型[5],SMA模型是用過去若干個(gè)時(shí)間點(diǎn)的平均值來預(yù)測(cè)當(dāng)前的時(shí)間值,但是忽略了各個(gè)時(shí)間點(diǎn)對(duì)當(dāng)前時(shí)間點(diǎn)的影響并不相同的問題,而實(shí)際上離當(dāng)前時(shí)間更接近的時(shí)間點(diǎn)通常都會(huì)對(duì)預(yù)測(cè)時(shí)間點(diǎn)造成更大的影響. 而WMA模型則可以在不同的時(shí)間點(diǎn)上分配不同的權(quán)值,讓那些對(duì)預(yù)測(cè)點(diǎn)影響更大的時(shí)間點(diǎn)擁有更大的權(quán)值,以減小過去較遠(yuǎn)時(shí)間點(diǎn)造成的干擾.

    針對(duì)以上問題,本研究選擇WMA模型對(duì)磁盤使用率進(jìn)行預(yù)測(cè),對(duì)傳統(tǒng)WMA模型進(jìn)行改進(jìn),以保證其在處理磁盤使用率數(shù)據(jù)時(shí)能有更好的預(yù)測(cè)效果. 在保證精度的前提下減小滯后值[6],并在傳統(tǒng)WMA模型上加入權(quán)重轉(zhuǎn)移[7],用來防止采集數(shù)據(jù)時(shí)的意外誤差對(duì)預(yù)測(cè)造成的影響,并結(jié)合拉依達(dá)準(zhǔn)則[8],將出現(xiàn)異常的時(shí)間點(diǎn)的權(quán)重值轉(zhuǎn)移到其他正常的時(shí)間點(diǎn)上,而不是像傳統(tǒng)的權(quán)重轉(zhuǎn)移法只對(duì)丟失值進(jìn)行處理,最后通過多新息遞推最小二乘算法[9]求出參數(shù)估計(jì),得到更好的擬合效果. 經(jīng)實(shí)驗(yàn)驗(yàn)證,與傳統(tǒng)的SMA和WMA比較,本算法可有效減小誤差,對(duì)磁盤使用率序列滯后現(xiàn)象的修補(bǔ)效果也十分良好.

    1磁盤使用率預(yù)測(cè)模型

    WMA模型屬于MA模型的一種. MA模型是一種常見的標(biāo)準(zhǔn)線性模型. 設(shè)εt是白噪聲序列, i為序列指標(biāo). 若時(shí)間序列yt滿足

    yt=εt-α1εt-1-α2εt-2-…-αqεt-q

    (1)

    則稱yt為q階滑動(dòng)平均序列,簡(jiǎn)稱MA(q)序列[10]. 式(1)中, αi是滑動(dòng)平均因子(moving average coefficients, MACs).

    本研究將WMA模型做出適當(dāng)改進(jìn),使其更符合磁盤使用率預(yù)測(cè)的要求,具體建模步驟如下:

    1) 通過磁盤使用率數(shù)據(jù)的自相關(guān)和偏自相關(guān)系數(shù)來對(duì)模型定階[11],確定的階數(shù)即為后面需要計(jì)算的參數(shù)個(gè)數(shù);

    2) 對(duì)時(shí)間序列進(jìn)行平穩(wěn)性檢驗(yàn)、消除奇異值和數(shù)據(jù)平滑[12]等處理,避免一些數(shù)據(jù)波動(dòng)造成的影響,減小隨機(jī)誤差,使后面要進(jìn)行的參數(shù)估計(jì)更加精確;

    3) 通過均方誤差公式計(jì)算得到最小誤差的最小滯后值,將得到的滯后值作為返回值參與計(jì)算,可在相同計(jì)算精度下減小滯后值;

    4) 使用結(jié)合了拉依達(dá)準(zhǔn)則的權(quán)重轉(zhuǎn)移法給各個(gè)時(shí)間點(diǎn)賦予權(quán)值,保證賦權(quán)的有效性;

    5) 用多新息遞推最小二乘算法進(jìn)行模型參數(shù)估計(jì)[13-21],將其值代替MA模型表達(dá)式中的參數(shù).

    至此模型建立完畢.依此模型,用采集得到的真實(shí)磁盤使用率數(shù)據(jù)進(jìn)行Matlab仿真,驗(yàn)證算法的有效性.

    1.1模型定階

    WMA模型的階數(shù)通過對(duì)應(yīng)MA模型的階數(shù)來確定. 本研究使用自相關(guān)和偏相關(guān)系數(shù)法來確定模型階數(shù),首先對(duì)序列的自相關(guān)系數(shù)(auto correlation function, ACF)和偏自相關(guān)系數(shù)(partial auto correlation function, PACF)進(jìn)行計(jì)算,然后通過對(duì)其拖尾性和截尾性的判斷來確定模型的類型及階數(shù). 拖尾性指這2個(gè)系數(shù)隨延遲k無限增長(zhǎng)以負(fù)指數(shù)的速度趨向于0;截尾性指它們?cè)趉大于模型階數(shù)后,其值變?yōu)?. 具體判斷規(guī)則如表1. 其中,AR(p)和MA(q)分別表示p階AR模型與q階MA模型.

    表1 模型及階數(shù)判斷標(biāo)準(zhǔn)

    設(shè)x1,x2,…,xT是平穩(wěn)的磁盤使用率時(shí)間序列XT的一個(gè)樣本,則其自協(xié)方差系數(shù)定義為

    (2)

    (3)

    樣本偏自相關(guān)系數(shù)定義為

    (4)

    其中,

    (5)

    由Bartlett公式,對(duì)于q>0, 自相關(guān)系數(shù)滿足正態(tài)分布

    (6)

    當(dāng)樣本容量充分大時(shí),有

    (7)

    (8)

    對(duì)于每個(gè)q>0, 檢查從第q+1個(gè)開始的ACF落入如式(9)范圍中的比例是否占總數(shù)的68.3%或95.5%左右.

    (9)

    如果在首階數(shù)q0之前,ACF都明顯不為0,而在q=q0時(shí),ACF中滿足式(9)的個(gè)數(shù)達(dá)到上述比例,則可判斷ACF序列在q0處截尾,即序列Xt為MA(q0)序列. 計(jì)算實(shí)際磁盤使用率序列得到的結(jié)果如表2.

    表2 自相關(guān)系數(shù)與偏自相關(guān)系數(shù)

    因表2數(shù)據(jù)滿足式(9)條件,則可判斷磁盤使用率數(shù)據(jù)在k=1處截尾,即q0=1, 故此時(shí)間序列是符合MA(1)模型的,即磁盤使用率時(shí)間序列是符合1階WMA模型的.

    1.2數(shù)據(jù)處理

    對(duì)磁盤使用率的數(shù)據(jù)處理包括平穩(wěn)性檢驗(yàn)、消除奇異值和平滑數(shù)據(jù)等方面. 表2自相關(guān)系數(shù)序列同時(shí)可以用作平穩(wěn)性檢驗(yàn)的標(biāo)準(zhǔn). 當(dāng)自相關(guān)系數(shù)序列快速衰減至0,則說明此時(shí)間序列是平穩(wěn)序列.由表2可見,磁盤使用率序列是平穩(wěn)序列.

    (10)

    (11)

    若序列中的值滿足式(12)條件則予以剔除,有

    (12)

    將與樣本均值的差值超過3倍標(biāo)準(zhǔn)差的數(shù)值確定為異常值,并將其用樣本均值代替.

    數(shù)據(jù)平滑的方法有滑動(dòng)窗口平均法、滑動(dòng)窗口擬合多項(xiàng)式平滑法、分量回歸平滑法、小波變換法和傅里葉變換法等. 磁盤使用率數(shù)據(jù)本身變化較為平緩,不需要過多的平滑處理,故為了節(jié)省運(yùn)算成本,并取得較好的平滑效果,本研究選取二次滑動(dòng)窗口平均法來進(jìn)行數(shù)據(jù)平滑,設(shè)定窗口大小為5個(gè)時(shí)間點(diǎn),原理如下:

    (13)

    其中, y(i)為第i個(gè)時(shí)間點(diǎn)的值; yy(j)為第j個(gè)時(shí)間點(diǎn)經(jīng)平滑處理后的值.

    對(duì)經(jīng)過1次平滑后的序列yy(j)用式(13)的方式進(jìn)行第2次平滑處理,經(jīng)2次處理后數(shù)據(jù)的隨機(jī)誤差減小. 用平滑后的序列進(jìn)行預(yù)測(cè),可以保證數(shù)據(jù)精度,提高預(yù)測(cè)準(zhǔn)確性.

    1.3滯后性的改進(jìn)

    由于滑動(dòng)平均因子的對(duì)稱性,WMA模型預(yù)測(cè)有大約(n+1)/2個(gè)時(shí)間點(diǎn)的滯后,而進(jìn)行短期的磁盤使用率預(yù)測(cè)時(shí)要求較小的滯后. 故在通過式(14)計(jì)算得到的一系列滯后值中,選取造成誤差最小的.

    (14)

    其中,

    (15)

    其中, α是滑動(dòng)平均因子. 根據(jù)原始輸入序列,可以得到輸出序列的滯后值為

    (16)

    式(16)使用了均方誤差作為衡量標(biāo)準(zhǔn),用于計(jì)算每個(gè)滯后值造成的誤差,其中造成最小誤差的滯后值就作為模型的因子來計(jì)算輸出,可有效減小磁盤使用率預(yù)測(cè)模型的滯后性問題.

    1.4改進(jìn)的權(quán)重轉(zhuǎn)移

    結(jié)合拉依達(dá)準(zhǔn)則法與權(quán)重轉(zhuǎn)移法,將不在假設(shè)檢驗(yàn)范圍內(nèi)的異常點(diǎn)的權(quán)重交移給當(dāng)前預(yù)測(cè)中權(quán)重最大的點(diǎn),針對(duì)磁盤使用率數(shù)據(jù)中出現(xiàn)的異常波動(dòng),加強(qiáng)優(yōu)化效果.

    傳統(tǒng)的權(quán)重轉(zhuǎn)移是基于WMA模型提出的,將模型中已分配權(quán)重但實(shí)際值不存在的時(shí)間點(diǎn)的權(quán)重,重新分配給對(duì)當(dāng)前預(yù)測(cè)目標(biāo)點(diǎn)影響最大的時(shí)間點(diǎn). 原始的權(quán)重轉(zhuǎn)移法只處理擁有權(quán)重但實(shí)際值為零的點(diǎn),即處理丟失值,但對(duì)其他異常數(shù)據(jù)并不敏感,只在極少數(shù)情況下才能起到作用,用途受到極大限制. 現(xiàn)將其與達(dá)依拉準(zhǔn)則結(jié)合,即在對(duì)每個(gè)當(dāng)前值進(jìn)行預(yù)測(cè)時(shí),將每個(gè)賦有權(quán)重的點(diǎn)依次排成序列,并對(duì)此權(quán)重序列進(jìn)行拉依達(dá)準(zhǔn)則的判斷,對(duì)于不滿足拉依達(dá)準(zhǔn)則的點(diǎn),將其權(quán)重值重新分配給對(duì)當(dāng)前預(yù)測(cè)目標(biāo)點(diǎn)影響最大的點(diǎn).

    如現(xiàn)有磁盤使用率數(shù)據(jù)中的5個(gè)時(shí)間點(diǎn)y1、y2、y3、y4和y5分別對(duì)應(yīng)權(quán)重值w1、w2、w3、w4和w5, 用以預(yù)測(cè)下一時(shí)間點(diǎn)y6, 根據(jù)加權(quán)滑動(dòng)平均公式有

    y6=y1w1+y2w2+y3w3+y4w4+y5w5

    (17)

    若通過計(jì)算發(fā)現(xiàn)y2滿足式(12),根據(jù)上述規(guī)則,將要除去的點(diǎn)的權(quán)值賦給對(duì)當(dāng)前目標(biāo)點(diǎn)影響最大的點(diǎn),而此例中y5距離目標(biāo)點(diǎn)y6最近,影響最大,故將y2的權(quán)值w2轉(zhuǎn)移到y(tǒng)5上,即此刻預(yù)測(cè)值滿足

    y6=y1w1+y3w3+y4w4+y5(w5+w2)

    (18)

    通過改進(jìn)的權(quán)重轉(zhuǎn)移,可有效避免因磁盤使用率數(shù)據(jù)的異常波動(dòng)而造成的預(yù)測(cè)誤差.

    1.5參數(shù)估計(jì)

    若將磁盤使用率預(yù)測(cè)模型描述為

    y(t)=D(z)v(t)

    D(z)=1+d1z-1+d2z-2+…+dvz-v

    (19)

    其中, y(t)為系統(tǒng)輸出觀測(cè)序列; v為白噪聲; z為單位后移算子; di為待計(jì)算參數(shù).且有

    z-1v(t)=v(t-1)

    (20)

    設(shè)整數(shù)p為新息長(zhǎng)度(innovation length), Y(t)為堆積輸出向量, V(t)為堆積誤差向量, Γ為信息矩陣, θ為對(duì)參數(shù)di的估計(jì)基于最小二乘優(yōu)化原理[21],有多新息遞推最小二乘法(multiple innovation recursive least squares,MILS)為

    (21)

    (22)

    (23)

    [φ(t),φ(t-1),…,φ(t-p+1)]T

    (24)

    具體步驟的偽代碼如下:

    Pseudocode of MILS

    set p0=106, P(0)=p0I,θ(0)=1n/p0;

    for t=1 to v do

    constuct P(t), V(p,t);

    for t=1 to n do

    choose a new p;

    end

    end

    該算法執(zhí)行了遞階計(jì)算過程,在設(shè)置初值后,用式(24)構(gòu)造Γ(p,t), 用式(22)構(gòu)造P(t), 用式(23)計(jì)算V(p,t), 用式(21)刷新估計(jì)值. 經(jīng)以上計(jì)算逐次刷新對(duì)di的估計(jì)值θ, 則可完成對(duì)磁盤使用率模型的所有參數(shù)估計(jì).

    通過第1節(jié)的各步驟,參數(shù)估計(jì)完畢,模型建立完成.

    2實(shí)驗(yàn)結(jié)果及分析

    為驗(yàn)證本算法的有效性,現(xiàn)對(duì)不同編號(hào)的10個(gè)磁盤進(jìn)行一段時(shí)間的使用率記錄,磁盤類型為jfs2,操作系統(tǒng)為aix. 本研究采取的記錄方式為每隔10 min對(duì)磁盤進(jìn)行1次采樣,每個(gè)磁盤的采樣點(diǎn)保證在300個(gè)以上. 得到每個(gè)磁盤使用率的時(shí)間序列,用于進(jìn)行Matlab仿真實(shí)驗(yàn).

    將傳統(tǒng)的WMA和SMA與本研究算法對(duì)比,用平均絕對(duì)誤差(mean absolute error, MAE)作為評(píng)判標(biāo)準(zhǔn). 平均絕對(duì)誤差是每個(gè)時(shí)間點(diǎn)上預(yù)測(cè)值與真實(shí)值之差絕對(duì)值的平均,相比平均誤差,能真實(shí)反映預(yù)測(cè)誤差,使結(jié)果具有良好的參考價(jià)值. MAE的計(jì)算公式為

    (25)

    其中, x(i)是原始時(shí)間序列; y(i)是預(yù)測(cè)時(shí)間序列.

    各算法的預(yù)測(cè)結(jié)果對(duì)比如表3.

    表3 算法平均絕對(duì)誤差對(duì)比

    由表3可以看到,本研究算法相對(duì)于原始的WMA,平均絕對(duì)誤差減小了24%~73%,在磁盤使用率序列上取得了良好的效果.

    本研究用Matlab實(shí)現(xiàn)仿真[22],對(duì)10組磁盤使用率數(shù)據(jù)進(jìn)行處理和預(yù)測(cè),部分實(shí)驗(yàn)結(jié)果如圖1.

    圖1 自相關(guān)系數(shù)Fig.1 Autocorrelation coefficients

    圖1中將表2的自相關(guān)系數(shù)用散點(diǎn)圖畫出,并將判斷是否落在范圍內(nèi)的標(biāo)準(zhǔn)差用與x軸平行的兩條實(shí)線標(biāo)記,可以更明確地看到Lag=1時(shí)的點(diǎn)落在實(shí)線外,而后面的值都迅速衰減到0左右,并在x軸附近進(jìn)行無規(guī)則的上下波動(dòng),根據(jù)以上情況可判斷該時(shí)間序列符合MA(1)模型.

    采集到的原始時(shí)間序列如圖2.

    圖2 采樣后的原始時(shí)間序列Fig.2 Original time series of sampling

    圖2是10組磁盤使用率序列中的一組,可以看到磁盤使用率時(shí)間序列極其平緩,甚至?xí)趲讉€(gè)連續(xù)時(shí)間點(diǎn)里數(shù)值保持不變,但在圖2中很明顯有2個(gè)點(diǎn)的數(shù)值比其他點(diǎn)高出很多,所以需要處理數(shù)據(jù),以避免一些非算法精度問題造成的影響. 圖3為處理后的時(shí)間序列.

    圖3 經(jīng)數(shù)據(jù)處理后的時(shí)間序列Fig.3 Time series after processing

    由圖3可以清晰看到磁盤使用率序列經(jīng)過600多次采樣,變化范圍都在94.8~95.0,明顯高于其他數(shù)據(jù)的值被判斷為異常值并予以剔除. 處理過后的序列消除了異常情況和隨機(jī)誤差的干擾,可以進(jìn)行預(yù)測(cè),其預(yù)測(cè)結(jié)果如圖4.

    圖4 預(yù)測(cè)結(jié)果擬合圖Fig.4 Fitting curve of prediction result

    圖4給出了本研究改進(jìn)的WMA模型的預(yù)測(cè)序列與原始時(shí)間序列的擬合圖,由于采樣點(diǎn)太多,此圖只能粗略地看到預(yù)測(cè)效果,所以截取了其中一小部分?jǐn)?shù)據(jù)并放大,與原始SMA和WMA算法對(duì)比.

    圖5截取了序列中變化較頻繁的一段,可以更好地看出擬合效果.由圖5可見,改進(jìn)后的WMA擬合效果優(yōu)于SMA和WMA,預(yù)測(cè)誤差明顯減小,能快速感知數(shù)據(jù)變化,且滯后也顯著減小,較好地解決了磁盤使用率預(yù)測(cè)中對(duì)高預(yù)測(cè)精度和較小滯后的需求.

    圖5 改進(jìn)算法與SMA和WMA的對(duì)比圖Fig.5 (Color online) Comparison diagram of prodiction by the proposed algorithm, SMA and WMA

    結(jié)語

    本研究提出了一種對(duì)WMA模型的改進(jìn)算法,對(duì)真實(shí)磁盤進(jìn)行使用率數(shù)據(jù)采樣,通過Matlab仿真驗(yàn)證了算法的有效性. 實(shí)驗(yàn)結(jié)果顯示本算法優(yōu)化效果良好,在減小預(yù)測(cè)誤差和滯后方面效果顯著.仿真過程中發(fā)現(xiàn),若數(shù)據(jù)出現(xiàn)連續(xù)的上下劇烈波動(dòng),則預(yù)測(cè)效果變差,即本算法對(duì)較尖銳的波峰波谷預(yù)測(cè)能力不強(qiáng),雖然這種數(shù)據(jù)出現(xiàn)劇烈波動(dòng)的情況在磁盤使用率數(shù)據(jù)中極為罕見,但也值得考慮并做出改進(jìn).

    引文:張宗華,張海全, 李師航,等. 基于加權(quán)滑動(dòng)平均的磁盤使用率預(yù)測(cè)模型[J]. 深圳大學(xué)學(xué)報(bào)理工版,2016,33(1):72-79.

    參考文獻(xiàn)/ References:

    [1] Murray S,Amaan M,Christoph A, et al. Projecting disk usage based on historical trends in cloud environment[C]// The 3rd ACM Workshop on Scientific Cloud Computing. Delft, Netherlands: Association for Computing Machinery, 2012.

    [2] Husna S H, Cui Lishan, Herny R H. Time series analysis of Web server logs for an online newspaper[C]// Proceedings of the 7th International Conference on Ubiquitous Information Management and Communication. Kota Kinabalu, Malaysia: Association for Computing Machinery, 2013.

    [3] Ahmad A S A,El-Shafie A,Naseri M, et al. Rainfall data analyzing using moving average (MA) model and wavelet multi-resolution intelligent model for noise evaluation to improve the forecasting accuracy[J]. Neural Computing and Applications, 2014, 25(7): 1853-1861.

    [4] Raudys A. Optimal negative weight moving average for stock price series smoothing[C]// Conference on Computational Intelligence for Financial Engineering. London: Institute of Electrical and Electronics Engineers, 2014.

    [5] Hansun Seng. A new approach of moving average method in time series analysis[C]// International Conference on New Media Studies.Tangerang,Indonesia: IEEE Computer Society, 2013.

    [6] Adrian L, Gao Junbin, Zheng Lihong. Optimizing the moving average[C]//World Congress on Computational Intelligence. Brisbane,Australia: Institute of Electrical and Electronics Engineers Inc, 2012.

    [7] 楊慧,許福栗. 基于權(quán)重轉(zhuǎn)移的加權(quán)滑動(dòng)平均模型改進(jìn)[J].計(jì)算機(jī)工程與應(yīng)用, 2014, 50(14):156-159.

    Yang Hui, Xu Fuli. Improvement of weighted moving average model based on transferring weights[J]. Computer Engineering and Applications, 2014, 50(14):156-159.(in Chinese)

    [8] Zhang Liang, Qin Yongyuan, Zhang Jinliang.Study of polynomial curve fitting algorithm for outlier elimination[C]// International Conference on Computer Science and Service System.Nanjing, China: IEEE Computer Society, 2011.

    [9] 周易,丁峰. 滑動(dòng)平均模型的最小二乘辨識(shí)方法比較研究[J].科學(xué)技術(shù)與工程, 2007, 18(7): 4570-4575.

    Zhou Yi, Ding Feng. Comparison of least squares identification for moving average models[J]. Science Technology and Engineering, 2007, 18(7): 4570-4575.(in Chinese)

    [10] 卜愛國,王炯. 基于MA(q)模型的動(dòng)態(tài)電源管理預(yù)測(cè)策略[J].計(jì)算機(jī)應(yīng)用研究,2011,28(7): 2516-2518.

    Bu Aiguo, Wang Jiong. Dynamic power management predictive policy based on MA(q) model[J]. Application Research of Computers, 2011,28(7): 2516-2518.(in Chinese)

    [11] 張賢達(dá). 用自相關(guān)確定MA模型階數(shù)的二種新方法[J].電子學(xué)報(bào), 1994, 22(7): 103-104.

    Zhang Xianda. Two new mothods for order determination of an MA model using autocorrelations[J]. Chinese Journal of Electronics, 1994, 22(7): 103-104.(in Chinese)

    [12] 王達(dá), 崔蕊.數(shù)據(jù)平滑技術(shù)綜述[J]. 電腦知識(shí)與技術(shù), 2009, 17(5):4507-4509.

    Wang Da, Cui Rui. Data smoothing technology summary[J]. Computer Knowledge and Technology, 2009, 17(5):4507-4509.(in Chinese)

    [13] 李守巨,霍軍周,曹麗娟. 盾構(gòu)機(jī)土壓平衡系統(tǒng)的ARMA模型及其參數(shù)估計(jì)[J]. 煤炭學(xué)報(bào), 2014, 39(11):2201-2205

    Li Shouju, Huo Junzhou, Cao Lijuan. Autoregressive moving average model and its parameter estimation for earth pressure balance system of shield[J]. Journal of China Coal Society, 2014, 39(11):2201-2205.(in Chinese)

    [14] 宋安超. 滑動(dòng)平均模型參數(shù)估計(jì)方法的仿真比較[J]. 大眾科技, 2011, 12(1):49-51.

    Song Anchao. Simulation and comparison of three methods of parametric estimation for moving average model[J]. Popular Science and Technology, 2011, 12(1):49-51.(in Chinese)

    [15] 汪海濱,龍俊波,查代奉. 分?jǐn)?shù)低階時(shí)頻滑動(dòng)平均模型參數(shù)估計(jì)[J]. 計(jì)算機(jī)工程與應(yīng)用, 2015, 50(20):178-182.

    Wang Haibin, Long Junbo, Zha Daifeng. Modeling and parameter estimation based on FLO-TFMA[J]. Computer Engineering and Applications, 2015, 50(20),178-182.(in Chinese)

    [16] 單銳,施蘇桐,劉文. 基于改進(jìn)共軛梯度思想的滑動(dòng)平均模型參數(shù)估計(jì)優(yōu)化方法[J]. 蘭州理工大學(xué)學(xué)報(bào), 2014, 40(1):144-147.

    Shan Rui, Shi Sutong, Liu Wen. The auto regression moving average model optimization method of parameter estimation based on the improved conjugate gradient thoughts[J]. Journal of Lanzhou University of Technology, 2014, 40(1):144-147.(in Chinese)

    [17] 鄧自立. 滑動(dòng)平均模型參數(shù)估計(jì)的Gevers-Wouters算法的指數(shù)收斂性[J]. 科學(xué)技術(shù)與工程, 2005, 20(5):1473-1484.

    Deng Zili. Exponential convergence of Gevers-Wouters algorithm for moving average model parameter estimation[J]. Science Technology and Engineering, 2005, 20(5):1473-1484.(in Chinese)

    [18] 單銳,劉雅寧,劉文. 改進(jìn)的差分自回歸移動(dòng)平均模型的共軛梯度參數(shù)估計(jì)法[J]. 河南科技大學(xué)學(xué)報(bào)自然科學(xué)版, 2015, 36(4):85-90.

    Shan Rui, Liu Yaning, Liu Wen. Improved conjugate gradient parameter estimation for autoregressive integrated moving average model[J]. Journal of Henan University of Science and Technology Natural Science, 2015, 36(4):85-90.(in Chinese)

    [19] 高艷普,王向東,王冬青. 多變量受控自回歸滑動(dòng)平均系統(tǒng)的極大似然辨識(shí)方法[J]. 山東大學(xué)學(xué)報(bào)工學(xué)版, 2015, 45(2):49-55.

    Gao Yanpu, Wang Xiangdong, Wang Dongqing. Maximum likelihood identification method for a multivariable controlled autoregressive moving average system[J]. Journal of Shandong University Engineering Science, 2015, 45(2):49-55.(in Chinese)

    [20] 李天一,鄭建榮.基于ARMAX模型的子空間辨識(shí)算法[J]. 計(jì)算機(jī)仿真, 2015, 32(1):310-313.

    Li Tianyi, Zheng Jianrong. A new subspace identification method based on ARMAX model[J]. Computer Simulation, 2015, 32(1):310-313.(in Chinese)

    [21] Ding Feng,Cheng Tongwen. Hierarchical least squares identification methods for multivariable ststems[J]. IEEE Transactions on Automatic Control,2005, 50(3): 397-402.

    [22] 王國峰,王子良,王太勇,等. Matlab在時(shí)間序列分析中的應(yīng)用[J]. 應(yīng)用科技, 2003, 30(5):36-38.

    Wang Guofeng, Wang Ziliang, Wang Taiyong, et al. Application of Matlab in time series analysis[J]. Applied Science and Technology, 2003, 30(5): 36-38.(in Chinese)

    【中文責(zé)編:坪梓;英文責(zé)編:遠(yuǎn)鵬】

    1) Process and Information Management Center, NARI Group Corporation, Nanjing 211106, Jiangsu Province, P.R.China

    2) College of Economics and Information Engineering, Southwestern University of Finance and Economics,

    Chengdu 611130,Sichuan Province, P.R.China

    3) College of Computer Science and Engineering, University of Electronic Science and Technology of China,

    Chengdu 611731, Sichuan Province, P.R.China

    Abstract:This paper proposes an improved weighted moving average (WMA) model to predict the usage of disks in the near future. Considering the characteristics of showing gentle change in disk usage and the requirements with small lag, we firstly utilize the autocorrelation and partial autocorrelation coefficient method to determine the order of the model. After processing the series, the minimum lag can be calculated on the premise without affecting the accuracy. Additionally, weights transferring combined with the Pauta criterion is used to balance the weight. At last, we estimate the parameters by using the multiple innovation recursive least squares to improve the result of prediction. According to the simulation result by Matlab, this algorithm is proved to have less result errors and a smaller lag. It provides a better prediction effect as compared with the original WMA model.

    Key words:computing technology; weighted moving average model; disk usage; autocorrelation and partial autocorrelation coefficient method; Pauta criterion; transferring weight; multiple innovation recursive least squares

    作者簡(jiǎn)介:張宗華(1977—),男,南京南瑞集團(tuán)公司工程師.研究方向:電力信息化. E-mail:hangzonghua@sgepri.sgcc.com.cn

    基金項(xiàng)目:國家自然科學(xué)基金資助項(xiàng)目(61300192) ;國家科技支撐計(jì)劃資助項(xiàng)目(2013BAH33F00);南京南瑞集團(tuán)公司研究基金資助項(xiàng)目(SGTYHT/14-XX-194)

    中圖分類號(hào):TP 301.6

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

    doi:10.3724/SP.J.1249.2016.01072

    猜你喜歡
    計(jì)算技術(shù)
    《計(jì)算技術(shù)與自動(dòng)化》2022年總目次索引
    基于云計(jì)算技術(shù)的FLAC3D軟件計(jì)算平臺(tái)的研發(fā)
    試論云計(jì)算技術(shù)在計(jì)算機(jī)網(wǎng)絡(luò)安全存儲(chǔ)中的應(yīng)用
    談云計(jì)算技術(shù)的安全性
    大數(shù)據(jù)、云計(jì)算技術(shù)在智慧城市中的應(yīng)用
    云計(jì)算技術(shù)在現(xiàn)代化辦公系統(tǒng)中的應(yīng)用
    電子制作(2017年8期)2017-06-05 09:36:15
    《物探化探計(jì)算技術(shù)》2016年1~6期總要目
    基于云計(jì)算技術(shù)的虛擬實(shí)訓(xùn)室設(shè)計(jì)與實(shí)現(xiàn)
    淺析云計(jì)算技術(shù)在高速公路收費(fèi)系統(tǒng)中的應(yīng)用
    淺談云計(jì)算技術(shù)下的互聯(lián)網(wǎng)金融應(yīng)用
    河南科技(2015年2期)2015-02-27 14:20:22
    国产1区2区3区精品| 久久精品亚洲精品国产色婷小说| 淫秽高清视频在线观看| 在线观看舔阴道视频| 日本a在线网址| 日本欧美视频一区| 亚洲国产精品合色在线| 国产色视频综合| 在线观看一区二区三区| 久久久国产成人免费| 他把我摸到了高潮在线观看| 国产精品免费视频内射| 免费搜索国产男女视频| 久久欧美精品欧美久久欧美| 亚洲人成电影免费在线| 黄色女人牲交| 国产精品亚洲一级av第二区| 国产高清videossex| 午夜免费激情av| 在线观看一区二区三区| 午夜两性在线视频| 亚洲avbb在线观看| 亚洲国产欧美一区二区综合| 国产精品自产拍在线观看55亚洲| 丁香欧美五月| 久久久久久人人人人人| 国产精品98久久久久久宅男小说| 窝窝影院91人妻| 久久久久久久久免费视频了| 性欧美人与动物交配| 后天国语完整版免费观看| 久久久精品欧美日韩精品| 国产成人欧美| 亚洲欧洲精品一区二区精品久久久| 国产伦人伦偷精品视频| 嫩草影视91久久| 免费在线观看日本一区| 一级片免费观看大全| 精品乱码久久久久久99久播| 亚洲一区中文字幕在线| 国产精品久久电影中文字幕| 十八禁人妻一区二区| 久久精品人人爽人人爽视色| 欧美一区二区精品小视频在线| 欧美日韩亚洲国产一区二区在线观看| 90打野战视频偷拍视频| 日本欧美视频一区| 亚洲熟女毛片儿| 91老司机精品| 视频区图区小说| 老汉色av国产亚洲站长工具| 少妇粗大呻吟视频| 9色porny在线观看| 一级作爱视频免费观看| 老汉色av国产亚洲站长工具| 女人高潮潮喷娇喘18禁视频| av网站在线播放免费| 嫁个100分男人电影在线观看| 美女高潮到喷水免费观看| 久久久国产欧美日韩av| 欧美日韩精品网址| 美女扒开内裤让男人捅视频| 亚洲国产欧美一区二区综合| 在线免费观看的www视频| 嫩草影视91久久| 欧美乱码精品一区二区三区| 老司机福利观看| 久久久久九九精品影院| 国产精品久久视频播放| 欧美一区二区精品小视频在线| 午夜免费鲁丝| 久久精品亚洲av国产电影网| 日韩欧美一区视频在线观看| 热99国产精品久久久久久7| 国产成人av激情在线播放| 久久久国产一区二区| videosex国产| 亚洲成人精品中文字幕电影 | 亚洲第一青青草原| 又黄又粗又硬又大视频| 精品熟女少妇八av免费久了| 亚洲中文av在线| 国产成人欧美| 一级毛片精品| 在线观看日韩欧美| 51午夜福利影视在线观看| 精品人妻在线不人妻| 999久久久国产精品视频| 老司机亚洲免费影院| 国产精品一区二区在线不卡| 成人手机av| 久久精品人人爽人人爽视色| 亚洲美女黄片视频| 欧美成人午夜精品| 法律面前人人平等表现在哪些方面| 国产精品久久久久久人妻精品电影| 午夜免费观看网址| 日韩免费av在线播放| 欧美精品亚洲一区二区| 级片在线观看| 一区二区三区精品91| 久久精品人人爽人人爽视色| 日韩大码丰满熟妇| 亚洲久久久国产精品| 黄色毛片三级朝国网站| 视频在线观看一区二区三区| 国产精品自产拍在线观看55亚洲| av在线播放免费不卡| 热re99久久精品国产66热6| 香蕉久久夜色| 黑人欧美特级aaaaaa片| 精品熟女少妇八av免费久了| 亚洲av熟女| 999精品在线视频| 久久九九热精品免费| 淫秽高清视频在线观看| 少妇裸体淫交视频免费看高清 | 精品一区二区三区四区五区乱码| 制服诱惑二区| 十八禁人妻一区二区| 欧美精品一区二区免费开放| 亚洲五月婷婷丁香| 一区在线观看完整版| 久久久国产成人免费| 亚洲人成伊人成综合网2020| 中文欧美无线码| 嫁个100分男人电影在线观看| 高清毛片免费观看视频网站 | 女人精品久久久久毛片| 亚洲欧美激情综合另类| 精品久久蜜臀av无| 欧美中文日本在线观看视频| 亚洲avbb在线观看| 黄色怎么调成土黄色| 免费一级毛片在线播放高清视频 | 中文字幕最新亚洲高清| 精品第一国产精品| 窝窝影院91人妻| 午夜福利影视在线免费观看| 黑人巨大精品欧美一区二区mp4| 纯流量卡能插随身wifi吗| 女人高潮潮喷娇喘18禁视频| 一级a爱片免费观看的视频| 日韩 欧美 亚洲 中文字幕| 亚洲第一青青草原| 欧美日韩中文字幕国产精品一区二区三区 | 男女午夜视频在线观看| 日韩欧美国产一区二区入口| 12—13女人毛片做爰片一| 国产男靠女视频免费网站| 精品久久久久久,| 黄色视频不卡| 少妇裸体淫交视频免费看高清 | 久热爱精品视频在线9| 一级a爱视频在线免费观看| 亚洲精品成人av观看孕妇| 少妇粗大呻吟视频| 国产三级在线视频| 国产在线观看jvid| 国产91精品成人一区二区三区| 19禁男女啪啪无遮挡网站| 国产三级黄色录像| 视频区欧美日本亚洲| av在线天堂中文字幕 | 欧美精品一区二区免费开放| 国产欧美日韩综合在线一区二区| 国产精品秋霞免费鲁丝片| 成人国语在线视频| 日本三级黄在线观看| 黄片播放在线免费| 午夜日韩欧美国产| 欧美精品啪啪一区二区三区| 麻豆久久精品国产亚洲av | 亚洲五月婷婷丁香| 精品国产一区二区久久| 国产成人欧美在线观看| aaaaa片日本免费| 高清av免费在线| 精品第一国产精品| 色播在线永久视频| 女同久久另类99精品国产91| 丝袜在线中文字幕| 日本欧美视频一区| 国产91精品成人一区二区三区| 91精品三级在线观看| 大陆偷拍与自拍| 久久久久久久久久久久大奶| 中文字幕人妻熟女乱码| 国产精品亚洲av一区麻豆| 国产精品一区二区在线不卡| 亚洲中文日韩欧美视频| 又黄又爽又免费观看的视频| 国产无遮挡羞羞视频在线观看| 日本免费一区二区三区高清不卡 | 亚洲成人免费电影在线观看| 精品电影一区二区在线| 一a级毛片在线观看| 亚洲精品美女久久av网站| 久久精品国产亚洲av高清一级| 中文欧美无线码| 国产精品美女特级片免费视频播放器 | 国产人伦9x9x在线观看| xxxhd国产人妻xxx| 激情在线观看视频在线高清| 亚洲欧美日韩另类电影网站| 激情视频va一区二区三区| 国产又色又爽无遮挡免费看| 久久久久久大精品| 超色免费av| 亚洲第一欧美日韩一区二区三区| 久久久久精品国产欧美久久久| netflix在线观看网站| 一级黄色大片毛片| 久久久国产欧美日韩av| 国产免费av片在线观看野外av| av欧美777| 国产精品二区激情视频| 99国产综合亚洲精品| 欧美成狂野欧美在线观看| 亚洲欧美精品综合久久99| 怎么达到女性高潮| 老司机深夜福利视频在线观看| 国产高清视频在线播放一区| 日韩精品免费视频一区二区三区| 岛国视频午夜一区免费看| 交换朋友夫妻互换小说| 国产成人一区二区三区免费视频网站| √禁漫天堂资源中文www| 免费一级毛片在线播放高清视频 | 黄色毛片三级朝国网站| av免费在线观看网站| 女人高潮潮喷娇喘18禁视频| 欧美黄色片欧美黄色片| 一夜夜www| 精品熟女少妇八av免费久了| 久久久国产一区二区| 精品高清国产在线一区| 国产精品1区2区在线观看.| 97碰自拍视频| 日韩欧美三级三区| 一个人观看的视频www高清免费观看 | 美女扒开内裤让男人捅视频| 免费在线观看日本一区| 欧美人与性动交α欧美软件| 999久久久精品免费观看国产| 女人高潮潮喷娇喘18禁视频| 亚洲精品美女久久久久99蜜臀| 母亲3免费完整高清在线观看| 色婷婷av一区二区三区视频| 国产1区2区3区精品| 亚洲午夜精品一区,二区,三区| 午夜精品国产一区二区电影| 免费av毛片视频| 欧美日韩乱码在线| 黄片大片在线免费观看| 法律面前人人平等表现在哪些方面| 在线国产一区二区在线| 国产精品亚洲一级av第二区| 波多野结衣av一区二区av| 国产一区二区三区综合在线观看| 精品国产乱码久久久久久男人| 中文字幕人妻熟女乱码| 人人澡人人妻人| 在线观看66精品国产| 国产在线观看jvid| 亚洲人成77777在线视频| 国产精华一区二区三区| 精品久久久久久成人av| 男男h啪啪无遮挡| 夜夜看夜夜爽夜夜摸 | 亚洲精品粉嫩美女一区| av网站免费在线观看视频| 中出人妻视频一区二区| 国产精品电影一区二区三区| 亚洲成人免费电影在线观看| 日韩欧美一区视频在线观看| 淫秽高清视频在线观看| 国产精品久久视频播放| 美女扒开内裤让男人捅视频| 黑人操中国人逼视频| 悠悠久久av| 亚洲久久久国产精品| 亚洲aⅴ乱码一区二区在线播放 | 99香蕉大伊视频| 精品日产1卡2卡| 淫妇啪啪啪对白视频| 水蜜桃什么品种好| 一级a爱片免费观看的视频| 嫩草影院精品99| 男人操女人黄网站| 久久久久精品国产欧美久久久| 热99国产精品久久久久久7| 日本五十路高清| 国产午夜精品久久久久久| 国产精品偷伦视频观看了| 亚洲精品国产色婷婷电影| 免费在线观看视频国产中文字幕亚洲| 亚洲国产欧美日韩在线播放| 亚洲视频免费观看视频| 男女下面插进去视频免费观看| 精品无人区乱码1区二区| 国产1区2区3区精品| 欧美日韩黄片免| 国产高清激情床上av| 亚洲av熟女| 精品欧美一区二区三区在线| 超碰成人久久| 日韩成人在线观看一区二区三区| 国产亚洲欧美在线一区二区| 桃红色精品国产亚洲av| www国产在线视频色| 精品国产国语对白av| 欧美性长视频在线观看| 新久久久久国产一级毛片| 成人亚洲精品一区在线观看| 99精品久久久久人妻精品| 怎么达到女性高潮| 丝袜在线中文字幕| 99热只有精品国产| 麻豆av在线久日| 免费久久久久久久精品成人欧美视频| 午夜激情av网站| 国产精品久久久人人做人人爽| 免费搜索国产男女视频| 亚洲第一欧美日韩一区二区三区| 亚洲一区中文字幕在线| 成人三级黄色视频| 啪啪无遮挡十八禁网站| 丰满人妻熟妇乱又伦精品不卡| 久热爱精品视频在线9| 久热这里只有精品99| 69av精品久久久久久| 熟女少妇亚洲综合色aaa.| 久久国产精品影院| 国产主播在线观看一区二区| 日日干狠狠操夜夜爽| 熟女少妇亚洲综合色aaa.| 亚洲五月色婷婷综合| 欧美午夜高清在线| 久久精品91蜜桃| 国产主播在线观看一区二区| a在线观看视频网站| 熟女少妇亚洲综合色aaa.| 亚洲中文日韩欧美视频| 伊人久久大香线蕉亚洲五| 女性生殖器流出的白浆| 中文欧美无线码| 日日爽夜夜爽网站| 久久人人97超碰香蕉20202| 女性生殖器流出的白浆| 午夜精品久久久久久毛片777| 久久久久久久久免费视频了| 老司机靠b影院| 女生性感内裤真人,穿戴方法视频| 18禁黄网站禁片午夜丰满| 精品卡一卡二卡四卡免费| 亚洲精品一区av在线观看| 国产成人精品无人区| 91成人精品电影| 亚洲男人天堂网一区| 久久性视频一级片| 久久国产亚洲av麻豆专区| 欧美+亚洲+日韩+国产| 久久中文字幕人妻熟女| 大码成人一级视频| 久热爱精品视频在线9| a级片在线免费高清观看视频| 在线观看免费午夜福利视频| 丝袜美足系列| 亚洲精品久久成人aⅴ小说| 99在线视频只有这里精品首页| 法律面前人人平等表现在哪些方面| 国产熟女xx| 国产精品久久久av美女十八| 大型av网站在线播放| 国产aⅴ精品一区二区三区波| 一区二区三区激情视频| 国产精品国产av在线观看| 在线观看午夜福利视频| 高潮久久久久久久久久久不卡| 悠悠久久av| 久久久精品欧美日韩精品| 国产97色在线日韩免费| 一二三四在线观看免费中文在| 国产一区二区三区综合在线观看| bbb黄色大片| 无限看片的www在线观看| 午夜a级毛片| 亚洲人成电影免费在线| 手机成人av网站| 久久精品影院6| 99久久精品国产亚洲精品| 757午夜福利合集在线观看| 亚洲精品一二三| 精品无人区乱码1区二区| 欧美 亚洲 国产 日韩一| 国产蜜桃级精品一区二区三区| 在线观看日韩欧美| 亚洲自偷自拍图片 自拍| av中文乱码字幕在线| 夜夜爽天天搞| 日日爽夜夜爽网站| 丁香欧美五月| 99riav亚洲国产免费| 国产亚洲精品综合一区在线观看 | ponron亚洲| 99国产综合亚洲精品| 亚洲avbb在线观看| 亚洲一区二区三区色噜噜 | 人人妻人人澡人人看| 99国产精品一区二区三区| 亚洲,欧美精品.| aaaaa片日本免费| 亚洲欧美精品综合一区二区三区| 国产一区二区激情短视频| 欧美另类亚洲清纯唯美| 亚洲国产欧美日韩在线播放| 亚洲欧美一区二区三区黑人| 国产aⅴ精品一区二区三区波| 91字幕亚洲| 夜夜爽天天搞| 精品一区二区三区av网在线观看| 99国产极品粉嫩在线观看| 日韩人妻精品一区2区三区| 精品高清国产在线一区| 一区福利在线观看| 日韩有码中文字幕| 亚洲色图av天堂| 久久久水蜜桃国产精品网| 亚洲欧美精品综合久久99| 国产成人欧美| 大型av网站在线播放| av国产精品久久久久影院| 日韩成人在线观看一区二区三区| 在线观看免费日韩欧美大片| 国产精品乱码一区二三区的特点 | 亚洲成人免费av在线播放| 日韩 欧美 亚洲 中文字幕| 久久午夜亚洲精品久久| 岛国在线观看网站| 老汉色av国产亚洲站长工具| 精品日产1卡2卡| 悠悠久久av| 亚洲黑人精品在线| 一区二区三区国产精品乱码| av在线播放免费不卡| 99久久国产精品久久久| 变态另类成人亚洲欧美熟女 | 欧美黄色片欧美黄色片| 日韩一卡2卡3卡4卡2021年| 精品一区二区三区视频在线观看免费 | 他把我摸到了高潮在线观看| 中出人妻视频一区二区| 12—13女人毛片做爰片一| 久久久久国内视频| 国产免费男女视频| 成年人免费黄色播放视频| 欧美中文综合在线视频| 97碰自拍视频| 亚洲国产欧美一区二区综合| 久久国产亚洲av麻豆专区| 亚洲第一青青草原| 亚洲免费av在线视频| 热re99久久精品国产66热6| 精品免费久久久久久久清纯| 国产激情欧美一区二区| 80岁老熟妇乱子伦牲交| 十分钟在线观看高清视频www| 女人被躁到高潮嗷嗷叫费观| 男女下面插进去视频免费观看| 一级毛片精品| 国产欧美日韩一区二区三| 久久久国产成人精品二区 | a级毛片在线看网站| 香蕉久久夜色| www.999成人在线观看| √禁漫天堂资源中文www| 涩涩av久久男人的天堂| 午夜精品在线福利| 激情在线观看视频在线高清| 日韩视频一区二区在线观看| 夜夜躁狠狠躁天天躁| 无限看片的www在线观看| 国产一区二区三区综合在线观看| 两个人看的免费小视频| 欧美中文综合在线视频| 久久亚洲精品不卡| 成人18禁在线播放| 99久久人妻综合| 热99re8久久精品国产| 大型黄色视频在线免费观看| 久久精品成人免费网站| 又黄又粗又硬又大视频| 99国产综合亚洲精品| 黑人操中国人逼视频| 亚洲av电影在线进入| 久久久久久久久中文| 久久香蕉国产精品| 国产1区2区3区精品| 国产三级黄色录像| 免费一级毛片在线播放高清视频 | 1024香蕉在线观看| 午夜成年电影在线免费观看| 亚洲片人在线观看| 午夜成年电影在线免费观看| 日韩欧美一区视频在线观看| 法律面前人人平等表现在哪些方面| av在线播放免费不卡| 国产亚洲欧美精品永久| 亚洲片人在线观看| 日本三级黄在线观看| 久久精品国产清高在天天线| 日韩大尺度精品在线看网址 | 18禁美女被吸乳视频| 日韩免费高清中文字幕av| 丝袜美腿诱惑在线| 麻豆成人av在线观看| 三级毛片av免费| 精品一品国产午夜福利视频| 日日摸夜夜添夜夜添小说| 免费在线观看完整版高清| 咕卡用的链子| 久久久久九九精品影院| 久久欧美精品欧美久久欧美| 丰满人妻熟妇乱又伦精品不卡| 国产一区二区激情短视频| 一边摸一边抽搐一进一出视频| 久久香蕉激情| 男女下面进入的视频免费午夜 | 国产色视频综合| 午夜影院日韩av| 国产日韩一区二区三区精品不卡| 久久中文字幕一级| 黑人巨大精品欧美一区二区mp4| 久久这里只有精品19| av视频免费观看在线观看| 男男h啪啪无遮挡| 亚洲五月天丁香| 自线自在国产av| 免费av中文字幕在线| 午夜福利免费观看在线| 精品少妇一区二区三区视频日本电影| 99热国产这里只有精品6| 精品一品国产午夜福利视频| 一级毛片高清免费大全| 1024视频免费在线观看| 女人高潮潮喷娇喘18禁视频| 欧美久久黑人一区二区| 久久人妻熟女aⅴ| 国产亚洲av高清不卡| 女人精品久久久久毛片| 亚洲成人国产一区在线观看| 咕卡用的链子| 两性夫妻黄色片| 搡老岳熟女国产| 免费在线观看日本一区| 亚洲aⅴ乱码一区二区在线播放 | 成人亚洲精品av一区二区 | 欧美日韩瑟瑟在线播放| 欧美日韩一级在线毛片| 在线观看一区二区三区| 99精品在免费线老司机午夜| 亚洲精品中文字幕一二三四区| 每晚都被弄得嗷嗷叫到高潮| 国产激情久久老熟女| 91成年电影在线观看| 手机成人av网站| 美女 人体艺术 gogo| 久久天堂一区二区三区四区| 伦理电影免费视频| 91av网站免费观看| 亚洲三区欧美一区| 中出人妻视频一区二区| 长腿黑丝高跟| 岛国在线观看网站| 久9热在线精品视频| 日本 av在线| 亚洲欧美一区二区三区久久| 十八禁人妻一区二区| 亚洲国产精品sss在线观看 | av视频免费观看在线观看| 韩国av一区二区三区四区| 香蕉国产在线看| 日本欧美视频一区| 国产高清国产精品国产三级| 免费女性裸体啪啪无遮挡网站| 99久久久亚洲精品蜜臀av| 日本三级黄在线观看| 女性生殖器流出的白浆| 国产亚洲精品久久久久久毛片| 成年女人毛片免费观看观看9| 19禁男女啪啪无遮挡网站| 亚洲专区国产一区二区| 成人国语在线视频| 黄色视频不卡| 女同久久另类99精品国产91| 成年女人毛片免费观看观看9| 久久香蕉国产精品| 在线观看免费视频日本深夜| 老司机午夜福利在线观看视频| 久久香蕉国产精品| 国产精品亚洲一级av第二区| 久久久久国产一级毛片高清牌| 女生性感内裤真人,穿戴方法视频| 777久久人妻少妇嫩草av网站| 精品一区二区三区四区五区乱码| 人人妻人人爽人人添夜夜欢视频| 亚洲av五月六月丁香网| 色婷婷av一区二区三区视频| 长腿黑丝高跟| 亚洲五月色婷婷综合| 91麻豆av在线| 免费高清在线观看日韩| 老汉色∧v一级毛片| 国产片内射在线|