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

    基于切比雪夫最佳逼近意義下的CPI指數(shù)變化通道的原理與方法

    2015-11-03 15:06:16顧樂民
    經(jīng)濟(jì)數(shù)學(xué) 2015年3期

    摘要中國2000-2015年的186組CPI指數(shù)所體現(xiàn)的變化規(guī)律,可以用通道原理來描述.通道類似于交通圖中的道路,依據(jù)切比雪夫最佳逼近原理,以CPI函數(shù)為中心線,以極小化最大正負(fù)誤差為邊界線,包含了所有的CPI數(shù)據(jù),是一條具有最大安全范圍意義的數(shù)據(jù)通道.假設(shè)端點(diǎn)數(shù)據(jù)是建立在零誤差基礎(chǔ)上的最大權(quán)重數(shù)據(jù),未來短期內(nèi)CPI指數(shù)的變化不存在突變,或存在突變但其誤差不大于通道內(nèi)的最大正負(fù)誤差,則通道的延伸可以預(yù)測未來CPI指數(shù)的變化,以及變化的波動范圍.介紹了CPI通道的基本原理與方法、適用范圍、判別法則、具體算法等.通過多個數(shù)據(jù)處理實(shí)例,說明了直線通道具有簡單易懂、直觀性強(qiáng)、計(jì)算方便、適用范圍廣、符合性較好等特點(diǎn),驗(yàn)證了未來短期CPI指數(shù)變化基本是可知與可控的.

    關(guān)鍵詞切比雪夫;最佳逼近;零誤差型極小極大逼近;CPI指數(shù);CPI通道

    中圖分類號 文獻(xiàn)標(biāo)識碼A

    Principle and Method of CPI Index Change Channel

    Based on Chebyshev Optimal Approximation Significance

    GU Lemin

    (Tongji University,College of Material Science and Engineering, Shanghai200092,China)

    AbstractThe 186 sets of CPI index data from 2000 to 2015 years in China reflected some rules, which can be described by Channel principle and the corresponding method. Channel is a data channel with the largest security sense,which based on Chebyshev approximation theory,to the CPI function as its center line, to minimize the maximum positive and negative error as its boundary line,and contains all the CPI data.Assuming that the endpoint data is the largest weight data based on zero error, there is no mutation of CPI index in the short term, or mutation but the error is less than channel's maximum positive and negative error,then the future CPI index can be given by the channel through the forecast value and fluctuation value.This article introduced the CPI channel principle and its method, the applicable scope, the judging rules,the concrete algorithm, etc.Through some examples of data processing, this paper shows that the straight line channel is easy to understand, intuitively simple and strong, convenient in the computation, broad in the applicable scope, good in the compliance, and also verifies that the scope of future shortterm changes in the CPI index is knowable and controllable.

    Key wordsChebyshev; optimal approximation; zeroerror type of minimax approximation; CPI index; CPI channel

    1引言

    居民消費(fèi)價格指數(shù)(Consumer Price Index,簡稱CPI指數(shù)),是普通消費(fèi)者所購買的物品與勞務(wù)的總費(fèi)用的衡量標(biāo)準(zhǔn),反映了一定時期內(nèi)價格變動程度和趨勢的相對數(shù).CPI指數(shù)不僅受商品價格的影響,比如糧價[1]、房價[2]等,也有對其權(quán)重經(jīng)常進(jìn)行調(diào)整的一個動態(tài)過程,這使得CPI指數(shù)變化具有隨機(jī)性大、難以找到一般變化規(guī)律、難以進(jìn)行預(yù)測等特征.對CPI預(yù)測理論及方法成為許多學(xué)者關(guān)注的問題,目前主要有,基于小波分解自回歸模型分析法[3]、VAR 模型法[4]、ARIMA模型法[5]、神經(jīng)網(wǎng)絡(luò)法[6]、灰色GM(1,1)模型法[7]等.

    2000年以來的CPI指數(shù)積累了186組數(shù)據(jù),數(shù)據(jù)是離散的,孤立的,數(shù)據(jù)之間的關(guān)系是不明確的,數(shù)據(jù)的來源是有一定統(tǒng)計(jì)誤差的.這些看似雜亂無章的數(shù)據(jù)背后,仿佛總有一只無形的手在操縱著CPI數(shù)據(jù)的變化,稱這只無形的手為“隱函數(shù)”.或許這只無形的手根本就不存在,因?yàn)镃PI指數(shù)的波動含有大量的“人類因素”,或許存在但目前難以找到,但這不影響探討的本質(zhì).任何一個運(yùn)動、變化、發(fā)展的事物,都存在其本質(zhì)的內(nèi)在規(guī)律,都是可以從變化的過程中找到.

    構(gòu)建隱函數(shù)的目的,是要用數(shù)學(xué)的方法來探索CPI變化的某些規(guī)律.CPI指數(shù)屬于一種“近程有序,遠(yuǎn)程無序”的數(shù)據(jù)變化形式,在較短的局部范圍內(nèi),其變化具有一定的規(guī)律可循.從長期全局范圍看,其變化呈大波動狀失去規(guī)律.這也就是說,用具有軌跡特征的曲線是難以描述這種變化的,必須用其它的方法,一種既包含著曲線又不局限于曲線的方法去描述.endprint

    經(jīng)濟(jì)數(shù)學(xué)第 32卷第3期

    顧樂民:基于切比雪夫最佳逼近意義下的CPI指數(shù)變化通道的原理與方法

    切比雪夫(P.L.Chebyshev,1821~1894)創(chuàng)立的最佳一致逼近原理,最早源于19世紀(jì)對機(jī)器的機(jī)械運(yùn)動按理想設(shè)計(jì)運(yùn)動的研究.將該原理運(yùn)用于CPI指數(shù)變化,可以構(gòu)建一條切比雪夫最佳逼近意義下的CPI指數(shù)變化通道.CPI指數(shù)變化是有限的變化,可以用2條曲線,1條稱為上界限,另1條稱為下界限,將所有的數(shù)據(jù)都囊括其中,并形成一條延伸的通道.通道將雜亂無章的數(shù)據(jù)加以規(guī)范和約束,而隱函數(shù)必定在通道之內(nèi),通過數(shù)學(xué)的方法可以找

    到最佳逼近意義下的隱函數(shù),使CPI的變化成為可知與可控.由于通道具有連續(xù)性,變化具有慣性,所以通道的外延具有一定的預(yù)測效應(yīng),可以推斷出未來可能的變化趨勢,為決策提供有價值的參考.

    切比雪夫最佳逼近意義下的數(shù)據(jù)通道的建立與應(yīng)用,文獻(xiàn)[8]有較為詳細(xì)的描述.由于CPI指數(shù)具有變化莫測的特殊性,從最簡單的“直線通道”入手,通過建立通道而闡述其基本的原理與方法,并用186組數(shù)據(jù)按序做30個數(shù)據(jù)處理的實(shí)例,以檢驗(yàn)預(yù)測的效果.

    2 原理與方法

    切比雪夫最佳逼近的核心是最大絕對值誤差極小化,由此構(gòu)成了極小極大曲線擬合法,適用于一個封閉系統(tǒng)內(nèi)的描述,具有廣泛的應(yīng)用[9].由于最大誤差一般都是在端點(diǎn)處出現(xiàn),在定義區(qū)間外是發(fā)散的,這從切比雪夫多項(xiàng)式的所有圖形中可以看出,所以不適合用于預(yù)測.對預(yù)測而言,預(yù)測的誤差越小越好,而不是最大甚至發(fā)散.零誤差型極小極大逼近是切比雪夫最佳逼近原理的一個推廣,是通過若干零誤差點(diǎn)限制端點(diǎn)誤差為最大,達(dá)到端點(diǎn)外誤差不發(fā)散的目的,其理論基礎(chǔ)是零誤差型切比雪夫多項(xiàng)式[10],以及相應(yīng)的預(yù)測理論[11].對于CPI指數(shù)的預(yù)測,只需提供1個零誤差點(diǎn),在坐標(biāo)系上是指最右端(簡稱端點(diǎn),下同)的數(shù)據(jù)點(diǎn),這使復(fù)雜的問題可以簡化敘述.

    2.1零誤差型極小極大法(ZMinimax method)

    數(shù)據(jù)(xi,yi) i=1,2,…,m是隱函數(shù)y(x)在有定義的區(qū)間內(nèi)給出的m個離散點(diǎn)組,為找到隱函數(shù)y(x),設(shè)擬合函數(shù)f(x)=f(x,a),其中參數(shù)a=(a1,a2,…,an) ,n≤m,而a1,a2,…,an為n個不全為零的實(shí)數(shù).為使f(x)盡可能接近y(x),設(shè)誤差函數(shù)r(x)=y(x)-f(x),而誤差值ri是誤差函數(shù)r(x)上的具體數(shù)值:

    ri=yi-f(xi,a),i=1,2,…,m. (1)

    零誤差型極小極大法,是將端點(diǎn)數(shù)據(jù)誤差設(shè)定為零誤差rm=0條件下,依據(jù)最大絕對值誤差極小化的準(zhǔn)則來選擇參數(shù)a,即依據(jù)

    rm=ym-f(xm,a)=0,E=max1≤i≤m-1ri=max1≤i≤m-1yi-f(xi,a)→min(2)

    而構(gòu)成的一種方法,是切比雪夫最大絕對值誤差極小化基本準(zhǔn)則的一個推廣.

    2.2解的實(shí)現(xiàn)

    如果零誤差型極小極大解存在,即存在a=a使

    rm=ym-f(xm,a*)=0,E*=max1≤i≤m-1ri=max1≤i≤m-1yi-f(xi,a*)=min. (3)

    則至少存在1個零誤差點(diǎn)和n個切比雪夫交錯點(diǎn)x1,x2,…,xn使

    rm=ym-f(xm,a*)=0,E*=maxj=1,2,…,nrj=maxj=1,2,…,nyj-f(xj,a*)=min. (4)

    稱參數(shù)a為極小極大最佳擬合參數(shù),稱f(xj,a)為極小極大最佳擬合方程,稱E*為最佳逼近值,它們構(gòu)成了零誤差型極小極大逼近的一般解.

    2.3零誤差

    所謂零誤差就是沒有誤差或誤差為零.當(dāng)將端點(diǎn)數(shù)據(jù)的誤差設(shè)定為零時,可使短期的預(yù)測得到保證.從曲線變化的一般規(guī)律來看,零誤差的兩端是以“-,0,+”或“+,0,-”形式出現(xiàn),越接近零點(diǎn),誤差絕對值就越小.這就給出一個提示,將端點(diǎn)誤差設(shè)定為零誤差即ym-f(xm)=0,則對于端點(diǎn)外臨近的數(shù)據(jù),其預(yù)測誤差rm+1=ym+1-f(xm+1)的絕對值也必定不會大,這為預(yù)測的準(zhǔn)確性提供了理論依據(jù),僅在出現(xiàn)大隨機(jī)誤差的特例下才會無法成立[11].

    2.4CPI通道

    某個變化的過程和狀態(tài)可用“通道”來描述,例如,處于下降(或上升)的通道之中等.通道,來往暢通的道路,與交通圖中的道路相似,是擬人化的表達(dá).零誤差型極小極大逼近意義下的CPI指數(shù)變化通道,簡稱CPI通道(下同)有以下幾個特征:

    1)通道的構(gòu)成.通道由1條中心線,2條邊界線共同構(gòu)建;中心線是CPI函數(shù)轉(zhuǎn)化的曲線,也稱為路線是前行的指導(dǎo)線;邊界線是距中心線兩旁±E的曲線.通道用Channel(x)表示:

    Channel(x)=fx±E(5)

    2)通道的作用.通道將離散的、有隨機(jī)誤差的數(shù)據(jù)加以分類和規(guī)范,它依據(jù)最佳逼近值±E將全部的數(shù)據(jù)規(guī)范在通道內(nèi),指出了位于邊界線上的數(shù)據(jù),由于偏離中心線最遠(yuǎn),屬于波動最大的大誤差數(shù)據(jù);通道將端點(diǎn)數(shù)據(jù)的誤差設(shè)定為零誤差,廢除了所有數(shù)據(jù)是權(quán)重相等的慣例,使權(quán)重往端點(diǎn)數(shù)據(jù)傾斜,并使對未來的預(yù)測建立在零誤差的基礎(chǔ)上;通道包含了所有的數(shù)據(jù),所以隱函數(shù)必定在通道內(nèi),通過數(shù)學(xué)方法可以找到最佳逼近意義下的隱函數(shù),或近似隱函數(shù).

    3)通道的意義.將理論的指導(dǎo)路線與實(shí)際行走的軌跡聯(lián)系在一起,數(shù)據(jù)沿著中心線前行,但實(shí)際是在偏離和糾正偏離中前行的;通道指出了數(shù)據(jù)變化的最大范圍,限定了安全的最大界線;在最大安全范圍內(nèi)去探尋隱函數(shù),從而找到CPI變化的某些規(guī)律,用于解釋過于、指導(dǎo)現(xiàn)在、預(yù)測未來.

    3具體算法

    由于CPI數(shù)據(jù)的隨機(jī)性,波動性,難預(yù)測性,用通道的原理和方法尋找數(shù)據(jù)之間的關(guān)系有較好的效果.下面用圖示法加以介紹,用的是直線型通道.圖1橫坐標(biāo)x是月,縱坐標(biāo)y是CPI指數(shù)值(無量綱).圖中參與計(jì)算的CPI數(shù)據(jù)yi是12個,均在虛線之內(nèi),虛線外的數(shù)據(jù)有1個,不參與計(jì)算.endprint

    3.1通道的構(gòu)成

    通道的中心線是擬合函數(shù)描述f(x)=f(x,a),與2條邊界線共同構(gòu)建成通道f(x)±E.通道內(nèi)包含了12個數(shù)據(jù).由式(4)產(chǎn)生的最大正負(fù)誤差點(diǎn)位于邊界線上,如圖1中的點(diǎn)A和點(diǎn)B.通道內(nèi)有1個零誤差點(diǎn),位于數(shù)據(jù)的端點(diǎn),如圖1中的點(diǎn)C.

    3.2零誤差點(diǎn)

    零誤差點(diǎn)是人為設(shè)置的點(diǎn),目的有3個:首先是將相等權(quán)重的數(shù)據(jù)變?yōu)椴坏葯?quán)重,一般而言,距離現(xiàn)在最近的數(shù)據(jù)應(yīng)該有較大的權(quán)重,而較遠(yuǎn)的權(quán)重可以較小,這對于預(yù)測而言是合理的,所以端點(diǎn)數(shù)據(jù)是權(quán)重最大的數(shù)據(jù).其次是使未來的預(yù)測建立在誤差為零的基礎(chǔ)上,這對于未來的預(yù)測誤差,難以判斷是正還是負(fù)而言,是合理的.再次是一般在零點(diǎn)附近的數(shù)據(jù)其誤差絕對值一般都是較小的,這使短期預(yù)測準(zhǔn)確性有了理論上的依據(jù).

    3.3邊界線

    大誤差數(shù)據(jù)的出現(xiàn),會使邊界線外移,使通道變寬.為使通道收窄,必須將最大(正負(fù))誤差極小化.圖中由點(diǎn)A和點(diǎn)B的這2個最大誤差又都是極小化的,所以邊界線也是極小化的,使通道收窄.從安全意義上說,切比雪夫最佳逼近意義下的通道,是最窄的通道,也是數(shù)據(jù)變化的最大安全范圍,超出這個范圍就有可能是不安全或欠安全的.

    3.4預(yù)測假設(shè)條件

    若未來短期CPI指數(shù)的變化不存在突變,或存在突變但其最大絕對誤差不大于通道內(nèi)的最大絕對值誤差,則通道的外推延伸能較好的給出未來CPI指數(shù)變化的趨勢.預(yù)測是建立在預(yù)測值f(xm+1)與最大正負(fù)誤差±E基礎(chǔ)上的,描述了CPI指數(shù)未來可能的數(shù)值與最大的波動范圍:f(xm+1)±E.從統(tǒng)計(jì)概率角度出發(fā),大部分變化不會超出f(xm+1)±E,這樣就使得CPI指數(shù)未來的變化成為可知與可控.

    3.5判斷法則與具體算法

    判斷法則主要是判斷異常解是否存在,以及如何處理的問題.

    最大絕對值誤差與最小絕對值誤差(即零誤差)之間在本質(zhì)上是不相容的,強(qiáng)制將原本是最大誤差的端點(diǎn)改為零誤差點(diǎn),會導(dǎo)致方程結(jié)構(gòu)大的改變,甚至?xí)狗匠坛霈F(xiàn)一種病態(tài)狀,結(jié)果是預(yù)測準(zhǔn)確性變差.異常解出現(xiàn)的原因一般在于數(shù)據(jù)的隨機(jī)性偏大,數(shù)學(xué)模型選擇的不當(dāng)所致.具體算法包含判斷法則,主要步驟如下.

    1)判斷:先用最小二乘法對數(shù)據(jù)進(jìn)行預(yù)處理,進(jìn)行判斷.若用最小二乘法的數(shù)據(jù)處理,其最大絕對值誤差出現(xiàn)在端點(diǎn),且該誤差值較其他誤差值明顯放大,則該點(diǎn)就是異常數(shù)據(jù)點(diǎn),其解將可能會出現(xiàn)異常.最小二乘法是個簡單方便的數(shù)據(jù)處理法,它所獲得的最大絕對值誤差一般也是極小極大法的最大絕對值誤差,所以用最小二乘法進(jìn)行預(yù)處理,用的是其方便與有效.

    2)處理:對于異常數(shù)據(jù),可以通過增加或減少數(shù)據(jù)數(shù)目,或轉(zhuǎn)移零誤差點(diǎn),或改變數(shù)學(xué)模型等方式進(jìn)行處理;若異常數(shù)據(jù)雖然存在,但誤差在允許的范圍內(nèi),可不作處理.

    3)求解:取直線方程為f(xi)=a+bxi,對于式(4)設(shè)1≤j,k≤m-1,j≠k,可以通過

    a=ym-2ym-yj-yk2xm-xj-xkxm,

    b=2ym-yj-yk2xm-xj-xk,

    E=yj-ym+2ym-yj-yk2xm-xj-xkxm-xj.(6)

    獲得參數(shù)a,b及逼近值E.取不同的 j,k,使最大的E為極小minmaxE=E*,從而獲得最佳逼近值E*,此時獲得的參數(shù)即為最佳參數(shù)a*,b*.

    4數(shù)值分析

    2000年1月-2015年6月我國CPI指數(shù)來自國家統(tǒng)計(jì)局,共186個,每年的數(shù)據(jù)是12個,歸為1組(2015年除外),共有15組數(shù)據(jù),先以2003年的數(shù)據(jù)處理為例.

    4.1實(shí)例一

    2003年1月到12月的CPI數(shù)據(jù)有12個,由判斷法則進(jìn)行預(yù)處理,用最小二乘法得到的方程用P2003(xi)表示:P2003(xi)=99.41+0.26xi,經(jīng)判斷,端點(diǎn)i=12不是最大誤差點(diǎn),可以運(yùn)用零誤差型極小極大逼近,得到的方程用f2003(xi)表示,其中xi=1,2,…,12:

    f2003(xi)=99+0.35xi.(7)

    極小化的最大絕對值誤差出現(xiàn)在i=1,9處,為max|r|=1.05,由此構(gòu)建的2003年CPI通道為:

    Channel2003(x)=f2003x±1.05,其中f2003x是f2003(xi)在去掉下標(biāo)“i”后的函數(shù)表達(dá),定義區(qū)間為[1,12].將x13=13代入,得2004年1月CPI指數(shù)預(yù)測值及波動值:103.55±1.05.已知2014年1月CPI指數(shù)是103.2,在預(yù)測的范圍之內(nèi).將預(yù)測值與實(shí)際值進(jìn)行比較,預(yù)測的誤差為0.34%.

    文中的圖1就是2003年1月至12月的CPI變化,以及CPI通道,在虛線外的點(diǎn)就是2014年1月的CPI指數(shù),圖中有關(guān)說明可參見前文.

    4.2實(shí)例二

    如法炮制,按序?qū)?000,2001,…,2014年1月至12月的數(shù)據(jù)處理結(jié)果列于表1.

    從表1結(jié)果看,CPI實(shí)際值都落在預(yù)測值及波動范圍之內(nèi),所以預(yù)測的結(jié)果是有效的.其中2012年帶*的數(shù)據(jù)屬于異常的數(shù)據(jù)處理,該年11月的CPI指數(shù)為102.0,但12月增至104.6,而次年1月又回落到102.0,計(jì)算表明12月份的數(shù)據(jù)屬于異常數(shù)據(jù),通過零誤差點(diǎn)遷移,取絕對值誤差最小的點(diǎn)為零誤差點(diǎn),獲得表中的方程式.

    4.3實(shí)例三

    表2是按序2000,2001,…,2014年,當(dāng)年6月至次年5月的12個CPI數(shù)據(jù)處理以及對次年6月預(yù)測情況,共15組.其中帶*的2013/06 - 2014/05因出現(xiàn)3個零誤差點(diǎn),故斜率為0,屬于特殊方程.為了與表1區(qū)別,方程用符號g(xi)表示.

    5結(jié)論

    提供的通道原理與方法,具有簡單易懂、直觀性強(qiáng)、計(jì)算方便、適用范圍廣、符合性較好等特點(diǎn),是曲線擬合的一種推廣,目的是使預(yù)測及預(yù)測誤差成為可知與可控.由于直線通道是一個簡單的通道,在進(jìn)一步的探討中還需要逐步加以完善.

    參考文獻(xiàn)

    [1]鄒正方,黎智,李迪.國際糧價波動對我國CPI影響的實(shí)證分析:以玉米、大豆為例[J].數(shù)學(xué)的實(shí)踐與認(rèn)識,2012,42(17):41-46.

    [2]黃飛雪,金建東.金融危機(jī)前后中國房價指數(shù)對CPI的影響[J].經(jīng)濟(jì)數(shù)學(xué),2010,27(3):64-72.

    [3]陳升,李星野.基于小波分解自回歸模型的CPI預(yù)測[J].統(tǒng)計(jì)與決策,2012,349(1):18-20.

    [4]李慶華.基于VAR 模型的中國消費(fèi)價格指數(shù)分析[J].華中師范大學(xué)學(xué)報(人文社會科學(xué)版).2006,45(4): 56-61.

    [5]張本麗等.基于ARIMA模型的山東省居民消費(fèi)價格指數(shù)分析[J].魯東大學(xué)學(xué)報,2010,26(3):285-288.

    [6]劉海萍.神經(jīng)網(wǎng)絡(luò)在CPI預(yù)測中的應(yīng)用[C].第五屆(2010)中國管理學(xué)年會——市場營銷分會場論文集,2010.

    [7]曾波.居民消費(fèi)價格指數(shù)的GM(1,1)模型預(yù)測[J].統(tǒng)計(jì)與決策,2009,289(13):7-8.

    [8]顧樂民.基于切比雪夫最佳逼近原理的俄羅斯人口變化通道[J].俄羅斯研究,2015,191(2):178-192.

    [9]Gu LEMIN.Minimax curve fitting method in application of CD production function—— with the grain yield data in China as the example[J].Journal of Modern Agriculture.2013,2(3):43-55.

    [10]Lemin GU.ZeroError Type of chebyshev polynomials[J].International Journal of Modeling and Optimization.2013(3):272-277.

    [11]顧樂民.預(yù)測型切比雪夫多項(xiàng)式[J].計(jì)算機(jī)工程與應(yīng)用.2012,48(7):34-38.endprint

    免费搜索国产男女视频| 午夜福利在线观看吧| 香蕉国产在线看| 欧美日韩亚洲综合一区二区三区_| 免费女性裸体啪啪无遮挡网站| 99riav亚洲国产免费| 国内精品久久久久久久电影| 午夜日韩欧美国产| 少妇 在线观看| 精品电影一区二区在线| 亚洲精品美女久久久久99蜜臀| 香蕉丝袜av| 精品一区二区三区四区五区乱码| 高潮久久久久久久久久久不卡| 亚洲一区中文字幕在线| 国产又爽黄色视频| 亚洲九九香蕉| 精品久久蜜臀av无| 啦啦啦韩国在线观看视频| 少妇 在线观看| 国产精品香港三级国产av潘金莲| 最近最新中文字幕大全免费视频| 色av中文字幕| 老熟妇乱子伦视频在线观看| 婷婷精品国产亚洲av在线| 国产一区二区激情短视频| 中国美女看黄片| 精品欧美一区二区三区在线| 久久久久免费精品人妻一区二区 | 国内少妇人妻偷人精品xxx网站 | 久久精品国产亚洲av香蕉五月| 亚洲人成77777在线视频| 久久午夜亚洲精品久久| √禁漫天堂资源中文www| 无限看片的www在线观看| 亚洲全国av大片| 欧美午夜高清在线| 欧美成人一区二区免费高清观看 | 中国美女看黄片| 99久久国产精品久久久| www.熟女人妻精品国产| 亚洲精品国产区一区二| videosex国产| 免费看日本二区| 人人妻人人看人人澡| 久久精品91蜜桃| 久久国产精品男人的天堂亚洲| 久久精品国产亚洲av香蕉五月| 性色av乱码一区二区三区2| 国内精品久久久久久久电影| 美国免费a级毛片| 国产高清视频在线播放一区| 在线视频色国产色| 婷婷精品国产亚洲av在线| 久久久久国产精品人妻aⅴ院| 极品教师在线免费播放| 亚洲国产精品999在线| 黄色成人免费大全| 熟女电影av网| 亚洲国产精品合色在线| 黑人欧美特级aaaaaa片| 黄色a级毛片大全视频| 精品久久久久久久久久免费视频| 精品日产1卡2卡| 午夜成年电影在线免费观看| 欧美大码av| 人人澡人人妻人| 精品国产国语对白av| 757午夜福利合集在线观看| 老司机靠b影院| 亚洲自拍偷在线| 亚洲性夜色夜夜综合| 亚洲人成伊人成综合网2020| 亚洲va日本ⅴa欧美va伊人久久| av天堂在线播放| 亚洲狠狠婷婷综合久久图片| 最新在线观看一区二区三区| 手机成人av网站| 国内揄拍国产精品人妻在线 | 美女扒开内裤让男人捅视频| 在线十欧美十亚洲十日本专区| 女人高潮潮喷娇喘18禁视频| 久久中文字幕人妻熟女| 99久久国产精品久久久| 国产精品国产高清国产av| 亚洲第一av免费看| 久久天躁狠狠躁夜夜2o2o| 精品午夜福利视频在线观看一区| 天天一区二区日本电影三级| 两个人免费观看高清视频| 国产精品久久视频播放| 色播在线永久视频| 此物有八面人人有两片| 叶爱在线成人免费视频播放| 久久精品国产99精品国产亚洲性色| 一本大道久久a久久精品| 欧美av亚洲av综合av国产av| 亚洲精华国产精华精| 制服人妻中文乱码| 国产爱豆传媒在线观看 | 最近在线观看免费完整版| 欧美国产日韩亚洲一区| 一进一出抽搐gif免费好疼| 一个人免费在线观看的高清视频| 中文字幕精品亚洲无线码一区 | 国产成年人精品一区二区| 久久精品成人免费网站| 日日爽夜夜爽网站| 国产亚洲精品第一综合不卡| 国产又色又爽无遮挡免费看| 观看免费一级毛片| 成人午夜高清在线视频 | 两个人看的免费小视频| 中文资源天堂在线| 国产亚洲精品久久久久久毛片| 99久久久亚洲精品蜜臀av| 美女大奶头视频| 一a级毛片在线观看| 十八禁网站免费在线| 亚洲成av人片免费观看| 国产激情偷乱视频一区二区| 国产精品影院久久| 18美女黄网站色大片免费观看| 午夜久久久久精精品| 亚洲 欧美 日韩 在线 免费| 午夜福利免费观看在线| 校园春色视频在线观看| 亚洲成人久久性| 色综合站精品国产| 国产成人啪精品午夜网站| 女同久久另类99精品国产91| 欧美性猛交黑人性爽| 国产精品亚洲av一区麻豆| 激情在线观看视频在线高清| 一二三四社区在线视频社区8| 国产极品粉嫩免费观看在线| 搡老熟女国产l中国老女人| 日韩 欧美 亚洲 中文字幕| 中国美女看黄片| 欧美最黄视频在线播放免费| 亚洲欧美日韩无卡精品| 国产亚洲av嫩草精品影院| 国产1区2区3区精品| 精品人妻1区二区| 亚洲精品中文字幕一二三四区| 色尼玛亚洲综合影院| 亚洲精品中文字幕在线视频| 亚洲七黄色美女视频| 中文在线观看免费www的网站 | 欧美日韩中文字幕国产精品一区二区三区| 校园春色视频在线观看| 成人国产一区最新在线观看| 夜夜爽天天搞| 男男h啪啪无遮挡| 久久久国产欧美日韩av| 一区二区三区精品91| 久久久久久久久中文| 99re在线观看精品视频| 可以在线观看的亚洲视频| 欧美乱色亚洲激情| 国产一卡二卡三卡精品| 亚洲成a人片在线一区二区| www日本在线高清视频| 香蕉av资源在线| 一级作爱视频免费观看| 欧美不卡视频在线免费观看 | 欧美乱妇无乱码| 一区二区三区国产精品乱码| 美女午夜性视频免费| 免费看十八禁软件| 搞女人的毛片| 99精品欧美一区二区三区四区| 黄色片一级片一级黄色片| 欧美绝顶高潮抽搐喷水| 亚洲国产欧洲综合997久久, | 日韩 欧美 亚洲 中文字幕| 禁无遮挡网站| 成人国产综合亚洲| 视频区欧美日本亚洲| 老司机福利观看| 欧美日韩黄片免| 欧美日韩亚洲综合一区二区三区_| svipshipincom国产片| 国产成人欧美| 亚洲黑人精品在线| 日日爽夜夜爽网站| 一a级毛片在线观看| 99在线视频只有这里精品首页| 欧美日本视频| 亚洲av成人av| 色综合欧美亚洲国产小说| 999久久久精品免费观看国产| 黄色视频,在线免费观看| 精品久久久久久久久久久久久 | 日韩大码丰满熟妇| 99国产精品一区二区三区| 精品久久久久久,| 黑丝袜美女国产一区| 亚洲精品国产一区二区精华液| 亚洲国产日韩欧美精品在线观看 | 精品久久久久久久久久免费视频| 国产不卡一卡二| 国产欧美日韩一区二区精品| 亚洲专区国产一区二区| 久久精品影院6| av福利片在线| 久久国产乱子伦精品免费另类| 男人舔奶头视频| 欧美激情 高清一区二区三区| 性色av乱码一区二区三区2| 日本一区二区免费在线视频| 在线视频色国产色| 国产成人精品久久二区二区免费| 999久久久国产精品视频| 国产免费av片在线观看野外av| 国产精品 国内视频| 国产成人av教育| 精品久久久久久久久久免费视频| 激情在线观看视频在线高清| 大型黄色视频在线免费观看| 亚洲精品国产精品久久久不卡| 中文字幕另类日韩欧美亚洲嫩草| 不卡一级毛片| 精品人妻1区二区| 精品久久久久久久毛片微露脸| 欧美成狂野欧美在线观看| 欧美大码av| 一级a爱视频在线免费观看| 日韩精品青青久久久久久| 最好的美女福利视频网| 国产一区二区在线av高清观看| 午夜精品在线福利| 欧美久久黑人一区二区| 亚洲av五月六月丁香网| avwww免费| 成人三级黄色视频| 美女国产高潮福利片在线看| 久久婷婷成人综合色麻豆| 国产91精品成人一区二区三区| 亚洲一卡2卡3卡4卡5卡精品中文| 国产精品99久久99久久久不卡| 亚洲人成77777在线视频| 久久伊人香网站| 久久精品人妻少妇| 十八禁人妻一区二区| 精品卡一卡二卡四卡免费| 国产欧美日韩一区二区三| 欧美黄色淫秽网站| 人人妻,人人澡人人爽秒播| 嫁个100分男人电影在线观看| 99国产精品一区二区三区| 亚洲中文av在线| 丰满人妻熟妇乱又伦精品不卡| 2021天堂中文幕一二区在线观 | 欧美精品啪啪一区二区三区| 黑丝袜美女国产一区| 黄片播放在线免费| 宅男免费午夜| 黄色丝袜av网址大全| 国产av在哪里看| 亚洲五月婷婷丁香| 国产亚洲av高清不卡| 欧美日韩亚洲综合一区二区三区_| 看免费av毛片| 国产精品久久久久久精品电影 | 精品乱码久久久久久99久播| 欧美日韩亚洲国产一区二区在线观看| 亚洲国产精品sss在线观看| 在线观看免费午夜福利视频| 亚洲电影在线观看av| 日本a在线网址| 99国产精品99久久久久| 亚洲精品美女久久久久99蜜臀| 国产久久久一区二区三区| 午夜影院日韩av| 男人操女人黄网站| 亚洲熟女毛片儿| 久久精品国产亚洲av高清一级| 在线国产一区二区在线| 啦啦啦 在线观看视频| 高清在线国产一区| 久久婷婷人人爽人人干人人爱| 日韩一卡2卡3卡4卡2021年| 亚洲精品在线美女| 男人的好看免费观看在线视频 | 亚洲国产看品久久| 精品一区二区三区四区五区乱码| 大型黄色视频在线免费观看| 久99久视频精品免费| 久久精品91蜜桃| 一个人观看的视频www高清免费观看 | 97碰自拍视频| 18禁黄网站禁片午夜丰满| 亚洲欧美精品综合一区二区三区| 日日干狠狠操夜夜爽| 91在线观看av| 国产精品自产拍在线观看55亚洲| 日韩大尺度精品在线看网址| 1024香蕉在线观看| 成人18禁高潮啪啪吃奶动态图| 国产主播在线观看一区二区| 老司机福利观看| 亚洲美女黄片视频| 亚洲一区高清亚洲精品| 亚洲自拍偷在线| 热99re8久久精品国产| 色哟哟哟哟哟哟| 午夜福利在线观看吧| 色老头精品视频在线观看| 欧美日韩乱码在线| 欧美乱码精品一区二区三区| 亚洲性夜色夜夜综合| 日日爽夜夜爽网站| av有码第一页| 精品久久久久久久人妻蜜臀av| 99热只有精品国产| 亚洲人成网站高清观看| 亚洲第一青青草原| 人人澡人人妻人| 丰满的人妻完整版| 国产人伦9x9x在线观看| 国产免费av片在线观看野外av| 午夜福利成人在线免费观看| 男女下面进入的视频免费午夜 | 日韩精品青青久久久久久| 国产熟女午夜一区二区三区| 亚洲成人久久性| 午夜精品在线福利| 老司机靠b影院| 欧美成狂野欧美在线观看| 国产亚洲精品久久久久久毛片| 亚洲五月天丁香| 亚洲一区中文字幕在线| 99久久无色码亚洲精品果冻| 曰老女人黄片| 免费人成视频x8x8入口观看| 国产1区2区3区精品| av在线播放免费不卡| 十八禁网站免费在线| 一级a爱片免费观看的视频| 搡老妇女老女人老熟妇| 特大巨黑吊av在线直播 | 99热只有精品国产| 桃色一区二区三区在线观看| 国产精品亚洲一级av第二区| 国产一区二区在线av高清观看| 精品日产1卡2卡| 老熟妇仑乱视频hdxx| 亚洲七黄色美女视频| 中文字幕久久专区| 国产一区在线观看成人免费| 欧美成狂野欧美在线观看| 超碰成人久久| АⅤ资源中文在线天堂| 狂野欧美激情性xxxx| 夜夜看夜夜爽夜夜摸| 叶爱在线成人免费视频播放| 亚洲在线自拍视频| 国产在线观看jvid| 操出白浆在线播放| 一二三四在线观看免费中文在| 国产午夜精品久久久久久| 精品一区二区三区av网在线观看| 亚洲成av片中文字幕在线观看| 国产成人av激情在线播放| 黄片小视频在线播放| 国产aⅴ精品一区二区三区波| 国产蜜桃级精品一区二区三区| 日韩一卡2卡3卡4卡2021年| 韩国精品一区二区三区| 很黄的视频免费| 成人国产一区最新在线观看| 成年免费大片在线观看| 精华霜和精华液先用哪个| 最新在线观看一区二区三区| 美女午夜性视频免费| 成熟少妇高潮喷水视频| 在线永久观看黄色视频| 精品久久久久久久久久免费视频| 欧美激情极品国产一区二区三区| 久久国产亚洲av麻豆专区| 亚洲无线在线观看| 香蕉久久夜色| 亚洲av第一区精品v没综合| 欧美精品啪啪一区二区三区| 精品久久久久久久末码| 精品一区二区三区视频在线观看免费| 久久天躁狠狠躁夜夜2o2o| 久久久久免费精品人妻一区二区 | 国产亚洲精品综合一区在线观看 | 国内毛片毛片毛片毛片毛片| 桃色一区二区三区在线观看| 最近最新中文字幕大全电影3 | 大型av网站在线播放| 一级毛片高清免费大全| 婷婷六月久久综合丁香| 久久久久久亚洲精品国产蜜桃av| 一边摸一边做爽爽视频免费| 香蕉丝袜av| 怎么达到女性高潮| 国产高清激情床上av| 桃色一区二区三区在线观看| 丝袜在线中文字幕| 亚洲国产欧美一区二区综合| 久久久久免费精品人妻一区二区 | 日韩一卡2卡3卡4卡2021年| 亚洲 欧美 日韩 在线 免费| cao死你这个sao货| 精品少妇一区二区三区视频日本电影| 国产欧美日韩一区二区精品| 99国产极品粉嫩在线观看| e午夜精品久久久久久久| 国产野战对白在线观看| 亚洲男人的天堂狠狠| 久热这里只有精品99| 制服人妻中文乱码| 美女高潮到喷水免费观看| 黑人操中国人逼视频| a级毛片a级免费在线| 国产亚洲精品av在线| 叶爱在线成人免费视频播放| 不卡av一区二区三区| 51午夜福利影视在线观看| 欧美不卡视频在线免费观看 | 91老司机精品| 中文字幕人妻丝袜一区二区| 欧美一级a爱片免费观看看 | 热99re8久久精品国产| 人人妻人人澡欧美一区二区| 亚洲欧美精品综合久久99| 久久精品夜夜夜夜夜久久蜜豆 | 亚洲av片天天在线观看| 俄罗斯特黄特色一大片| 757午夜福利合集在线观看| 最近最新免费中文字幕在线| 亚洲电影在线观看av| 欧美日韩福利视频一区二区| 国产精品影院久久| 亚洲第一欧美日韩一区二区三区| 啦啦啦观看免费观看视频高清| 天堂√8在线中文| 一级毛片高清免费大全| 亚洲av第一区精品v没综合| 婷婷精品国产亚洲av| 最近最新免费中文字幕在线| 国产成人欧美| 成人国语在线视频| 国产欧美日韩一区二区三| 人妻丰满熟妇av一区二区三区| 午夜福利免费观看在线| 日韩精品中文字幕看吧| 给我免费播放毛片高清在线观看| 大型黄色视频在线免费观看| 婷婷丁香在线五月| 亚洲片人在线观看| 亚洲av电影不卡..在线观看| 久久婷婷成人综合色麻豆| 97碰自拍视频| 又大又爽又粗| 午夜激情av网站| 亚洲av第一区精品v没综合| 久久久久免费精品人妻一区二区 | 亚洲av电影在线进入| 亚洲在线自拍视频| 91麻豆精品激情在线观看国产| 欧美丝袜亚洲另类 | 视频区欧美日本亚洲| 日韩欧美国产在线观看| 两个人免费观看高清视频| 中出人妻视频一区二区| 亚洲片人在线观看| 久久草成人影院| 国产男靠女视频免费网站| 欧美色视频一区免费| 麻豆成人av在线观看| a级毛片在线看网站| 淫妇啪啪啪对白视频| 一区二区日韩欧美中文字幕| 国产日本99.免费观看| 人人妻人人澡欧美一区二区| bbb黄色大片| 少妇粗大呻吟视频| 一本一本综合久久| 日本免费一区二区三区高清不卡| 中文字幕精品亚洲无线码一区 | 麻豆av在线久日| 欧美 亚洲 国产 日韩一| 成在线人永久免费视频| 久久久久免费精品人妻一区二区 | 亚洲一区高清亚洲精品| 自线自在国产av| 日本一区二区免费在线视频| 久久欧美精品欧美久久欧美| 最新在线观看一区二区三区| 村上凉子中文字幕在线| x7x7x7水蜜桃| 特大巨黑吊av在线直播 | 成在线人永久免费视频| 亚洲精品av麻豆狂野| 亚洲av成人av| 看片在线看免费视频| 俄罗斯特黄特色一大片| 少妇熟女aⅴ在线视频| 99久久无色码亚洲精品果冻| 一级a爱视频在线免费观看| 熟女电影av网| 国内毛片毛片毛片毛片毛片| 色婷婷久久久亚洲欧美| 在线天堂中文资源库| 精品久久久久久久人妻蜜臀av| 日韩视频一区二区在线观看| 久久 成人 亚洲| 在线免费观看的www视频| 伦理电影免费视频| 亚洲av成人av| 亚洲av片天天在线观看| 国产亚洲精品久久久久5区| 白带黄色成豆腐渣| 最近最新免费中文字幕在线| 午夜两性在线视频| 无人区码免费观看不卡| 亚洲国产欧美日韩在线播放| 欧美成狂野欧美在线观看| 国内少妇人妻偷人精品xxx网站 | 老汉色∧v一级毛片| 中文字幕另类日韩欧美亚洲嫩草| www日本在线高清视频| 欧美成人一区二区免费高清观看 | 两个人看的免费小视频| 欧美黑人欧美精品刺激| 人人澡人人妻人| 美女高潮到喷水免费观看| 一级毛片高清免费大全| 俄罗斯特黄特色一大片| 亚洲专区字幕在线| 精品久久久久久久末码| 999久久久精品免费观看国产| 三级毛片av免费| 午夜影院日韩av| √禁漫天堂资源中文www| 国产成年人精品一区二区| 色综合婷婷激情| 日韩欧美免费精品| 国产伦一二天堂av在线观看| 午夜福利免费观看在线| 久久久久免费精品人妻一区二区 | 桃色一区二区三区在线观看| 久久国产乱子伦精品免费另类| 亚洲国产精品久久男人天堂| 亚洲七黄色美女视频| 日韩精品中文字幕看吧| 最近最新免费中文字幕在线| 黑人巨大精品欧美一区二区mp4| 午夜福利视频1000在线观看| 亚洲精品中文字幕在线视频| 久久久精品国产亚洲av高清涩受| 日本五十路高清| 熟女少妇亚洲综合色aaa.| 精品国产美女av久久久久小说| 亚洲欧美日韩高清在线视频| 成人亚洲精品一区在线观看| 老司机深夜福利视频在线观看| 亚洲黑人精品在线| 日韩中文字幕欧美一区二区| 高清毛片免费观看视频网站| 欧美亚洲日本最大视频资源| 欧美日韩一级在线毛片| 一区二区三区高清视频在线| 国产麻豆成人av免费视频| 人人妻,人人澡人人爽秒播| 在线av久久热| 亚洲第一av免费看| 熟女少妇亚洲综合色aaa.| 久久亚洲精品不卡| 亚洲av电影不卡..在线观看| 久久精品国产亚洲av香蕉五月| 一本精品99久久精品77| 国产蜜桃级精品一区二区三区| 一本精品99久久精品77| 天堂动漫精品| 波多野结衣巨乳人妻| 亚洲性夜色夜夜综合| 超碰成人久久| 夜夜爽天天搞| 女同久久另类99精品国产91| 桃红色精品国产亚洲av| 国产一级毛片七仙女欲春2 | 久久香蕉精品热| 午夜精品在线福利| 国产高清视频在线播放一区| 午夜激情福利司机影院| 国产av在哪里看| 可以免费在线观看a视频的电影网站| www日本黄色视频网| 色尼玛亚洲综合影院| 亚洲,欧美精品.| 中文字幕av电影在线播放| 黄色视频,在线免费观看| 日本熟妇午夜| 久久久久久久精品吃奶| 亚洲五月色婷婷综合| 最近在线观看免费完整版| 韩国精品一区二区三区| 最新在线观看一区二区三区| avwww免费| 久久久久精品国产欧美久久久| 亚洲久久久国产精品| 黄色毛片三级朝国网站| 亚洲九九香蕉| 免费无遮挡裸体视频| 午夜免费激情av| 一进一出抽搐gif免费好疼| 很黄的视频免费| 美女大奶头视频| 丰满人妻熟妇乱又伦精品不卡|