亚洲免费av电影一区二区三区,日韩爱爱视频,51精品视频一区二区三区,91视频爱爱,日韩欧美在线播放视频,中文字幕少妇AV,亚洲电影中文字幕,久久久久亚洲av成人网址,久久综合视频网站,国产在线不卡免费播放

        ?

        基于模糊近似支持向量回歸的股價(jià)預(yù)測研究

        2021-04-06 10:53:08王冰玉劉勇軍
        關(guān)鍵詞:模型

        王冰玉,劉勇軍

        (華南理工大學(xué),廣東 廣州 510640)

        0 引 言

        股價(jià)預(yù)測主要基于相關(guān)預(yù)測模型預(yù)測股價(jià)的未來變化趨勢以捕捉相應(yīng)的市場行情,從而促進(jìn)股票選擇。近年來,股價(jià)預(yù)測是熱點(diǎn)問題。許多學(xué)者和業(yè)界人士都從理論和實(shí)踐層面對(duì)其股價(jià)預(yù)測模型進(jìn)行相關(guān)研究,其研究方法包括GARCH[1-3]、模糊時(shí)間序列[4-5]、ARIMA[6-8]等。由于計(jì)算機(jī)技術(shù)的發(fā)展和廣泛應(yīng)用,各種機(jī)器學(xué)習(xí)和量化投資模型也被逐漸地應(yīng)用到股票預(yù)測領(lǐng)域,例如支持向量機(jī)、人工神經(jīng)網(wǎng)絡(luò)、遺傳算法和差分進(jìn)化算法等股價(jià)預(yù)測模型。

        Yu和Yan基于PSR方法和DL的長短期記憶網(wǎng)絡(luò)(LSTMs),設(shè)計(jì)了一種基于DNN的股價(jià)預(yù)測模型,并對(duì)不同時(shí)期的多個(gè)股票指數(shù)進(jìn)行預(yù)測[9]。Onoh等利用協(xié)調(diào)搜索(HS)和遺傳算法,建立混合人工神經(jīng)網(wǎng)絡(luò)模型,并根據(jù)數(shù)據(jù)集的統(tǒng)計(jì)和財(cái)務(wù)表現(xiàn)進(jìn)行實(shí)證分析[10]。Moghaddam等研究人工神經(jīng)網(wǎng)絡(luò)對(duì)納斯達(dá)克股票日匯率的預(yù)測能力,從而對(duì)用反向傳播算法訓(xùn)練的幾種前饋神經(jīng)網(wǎng)絡(luò)進(jìn)行了評(píng)估[11]。Hjek等將情緒與財(cái)務(wù)指標(biāo)相結(jié)合,使用一個(gè)多層感知器神經(jīng)網(wǎng)絡(luò)來預(yù)測異常股票回報(bào)率[12]。綦方中等提出一種基于PCA-IFOA-BP神經(jīng)網(wǎng)絡(luò)的股票價(jià)格預(yù)測模型,并借助上證指數(shù)進(jìn)行預(yù)測[13]。孟葉等選取K-近鄰、梯度提升和自適應(yīng)提升這3個(gè)分類器,通過改進(jìn)的投票算法聚合成一個(gè)新的分類器模型,對(duì)指數(shù)行情數(shù)據(jù)進(jìn)行學(xué)習(xí)分類[14]。

        雖然基于神經(jīng)網(wǎng)絡(luò)的混合模型是預(yù)測股票市場指數(shù)的有效預(yù)測方法,但是卻存在著諸如黑箱技術(shù)、過擬合、收斂速度慢、陷入局部極小等局限性。為了克服這些局限性,Cortes和Vapnik提出的支持向量機(jī)(SVM)方法已經(jīng)成為股票市場指數(shù)預(yù)測領(lǐng)域的一種流行研究方法,它采用結(jié)構(gòu)風(fēng)險(xiǎn)最小化的原則,以最小化泛化誤差的上界。通過應(yīng)用支持向量機(jī),不太可能出現(xiàn)過擬合,而且最優(yōu)解也可能是全局的[15]。Nayak等構(gòu)建一種支持向量機(jī)與K近鄰法相結(jié)合的混合框架,用于印度股市指數(shù)的短期、中期和長期預(yù)測[16]。Chen等引入特征加權(quán)向量,進(jìn)一步提出特征加權(quán)支持向量機(jī)和特征加權(quán)K最近鄰的基本混合框架,以有效地預(yù)測股票市場指數(shù)[17]。Xiao等提出一種將SVM和SSA相結(jié)合的組合模型,并建立基準(zhǔn)模型進(jìn)行對(duì)比分析[18]。Lee創(chuàng)建一種基于支持向量機(jī)(SVM)和混合特征選擇方法的股市趨勢預(yù)測模型[19]。Fung等提出一種新的近端支持向量機(jī)(PSVM)模型,并使用公開的數(shù)據(jù)集上的計(jì)算結(jié)果表明,PSVM不僅具有與SVM相當(dāng)?shù)臏y試集正確性,而且具有相當(dāng)快的計(jì)算時(shí)間[20]。姚瀟和余樂安在PSVM的基礎(chǔ)上,引入模糊隸屬度的思想,提出模糊近似支持向量機(jī)(FPSVM),并利用兩個(gè)公開的信用數(shù)據(jù)集進(jìn)行實(shí)證研究驗(yàn)證該模型的有效性[21]。張貴生等針對(duì)股價(jià)的非線性特點(diǎn),提出基于近鄰互信息特征選擇SVM-GARCH的股價(jià)預(yù)測模型[22]。張冰等提出具有局部信息挖掘功能的DNN加權(quán)算法對(duì)eplion-TSVR模型進(jìn)行改進(jìn),并借助上證A股的高頻數(shù)據(jù)進(jìn)行實(shí)證預(yù)測[23]。

        綜上所述,首先雖然SVR能夠在股價(jià)預(yù)測中發(fā)揮更好的表現(xiàn)力,但是都未考慮股票數(shù)據(jù)中所含的噪聲信息對(duì)股價(jià)預(yù)測的影響,其次現(xiàn)有研究未曾考慮投資者對(duì)兩種預(yù)測誤差((1)預(yù)測值>實(shí)際值;(2)預(yù)測值<實(shí)際值)的不同偏好。有的投資者更注重收益,而有的投資者更關(guān)注損失,因此,不能簡單直接地使用SVR模型進(jìn)行股價(jià)預(yù)測,而是應(yīng)該根據(jù)實(shí)際情況,針對(duì)投資者的不同目標(biāo)賦予不同的偏好值,構(gòu)建有效的股價(jià)預(yù)測模型進(jìn)行股票價(jià)格的準(zhǔn)確預(yù)測。

        綜合考慮以上研究現(xiàn)狀和存在的不足之處,該文構(gòu)建信噪比特征向量,并借用現(xiàn)有文獻(xiàn)表明的相關(guān)股價(jià)預(yù)測指標(biāo),選取歷史數(shù)據(jù)、趨向、反趨向、能量、量價(jià)、波動(dòng)和信噪比等其他七個(gè)方面的指標(biāo)作為輸入變量,考慮到投資者對(duì)預(yù)測誤差的不同偏好,引入模糊隸屬度和雙邊權(quán)重測量方法,構(gòu)建基于信噪比的模糊近似支持向量回歸模型進(jìn)行股價(jià)預(yù)測。

        1 模糊近似支持向量機(jī)

        模糊近似支持向量機(jī)(FPSVM)是由姚瀟和余樂安在2012年提出的[21]。該模型是為了減小訓(xùn)練樣本的奇異點(diǎn)和噪聲對(duì)模型的干擾,在PSVM的基礎(chǔ)上,引入了模糊隸屬度。設(shè)含有N個(gè)訓(xùn)練樣本的訓(xùn)練集樣本對(duì){(x1,y1),(x2,y2),…,(xN,yN)},其中xi∈Rd為第i個(gè)訓(xùn)練樣本的輸入向量,d為樣本空間維度,yi∈{+1,-1}為對(duì)應(yīng)輸出值,則FPSVM模型的具體形式如下:

        (1)

        其中,C>0為錯(cuò)誤項(xiàng)的懲罰參數(shù),εi為松弛變量,Φ(x)為非線性映射函數(shù),wT為特征空間維數(shù),b為待定的標(biāo)量參數(shù),mi為隸屬度,表示第i個(gè)樣本點(diǎn)對(duì)超平面的貢獻(xiàn)率。mi越小,則誤差項(xiàng)所占的比例越小,在整個(gè)模型的影響就越小。當(dāng)所有的mi=1時(shí),F(xiàn)PSVM就退化為PSVM。

        2 股價(jià)預(yù)測輸入指標(biāo)的獲取

        由于買賣反彈,價(jià)格變化的離散性、交易規(guī)模的差異和訂單流的戰(zhàn)略組成部分等因素導(dǎo)致所觀察到的價(jià)格過程是一個(gè)包含噪聲的過程[24]?,F(xiàn)實(shí)生活中,一般所觀察的對(duì)數(shù)價(jià)格過程包含對(duì)數(shù)有效價(jià)格和噪聲兩部分,相應(yīng)的股價(jià)對(duì)數(shù)形式為:

        (2)

        (3)

        噪聲的存在會(huì)影響股價(jià)預(yù)測。若該文能在SVM方程的構(gòu)建中進(jìn)一步消除噪聲,則可提高股價(jià)預(yù)測精度。信噪比(SNR)是描述信號(hào)中有效成分與噪聲成分的比例關(guān)系參數(shù),可以有效地對(duì)噪聲進(jìn)行處理[25]。該文考慮引入SNR降低噪聲干擾。根據(jù)現(xiàn)有文獻(xiàn),可得到信噪比(SNR)的計(jì)算公式[26]:

        (4)

        此外,股票市場中,開盤價(jià)、最低價(jià)、最高價(jià)和總交易量歷史數(shù)據(jù)通常被用做輸入指標(biāo)。最近相關(guān)學(xué)者研究表明一些技術(shù)指標(biāo)有助于更好地預(yù)測股價(jià)[17]。故該文同時(shí)借助相關(guān)的技術(shù)指標(biāo)作為輸入變量,具體將這些指標(biāo)分為:

        (1)趨向指標(biāo):移動(dòng)平均線(MA)、指數(shù)移動(dòng)平均線(EMA)、異同移動(dòng)平均線(MACD)和動(dòng)量指標(biāo)(MTM);

        (2)反趨向指標(biāo):相對(duì)強(qiáng)度指數(shù)(RSI);

        (3)能量指標(biāo):AR、BR和成交量變異率(VR);

        (4)量價(jià)指標(biāo):多空比率凈額(DK);

        (5)波動(dòng)指標(biāo):平均真實(shí)范圍(ATR)。

        3 基于FPSVR的股價(jià)預(yù)測模型構(gòu)建

        3.1 FPSVR模型構(gòu)建

        直接采用SVR模型進(jìn)行股票預(yù)測會(huì)存在以下不足:首先原有的SVR模型都未考慮投資者對(duì)股價(jià)預(yù)測誤差的不同偏好情況,默認(rèn)投資者對(duì)預(yù)測誤差的偏好是一致的,使得預(yù)測結(jié)果不能準(zhǔn)確反映投資者的投資策略,其次SVR模型無法處理噪聲和奇異點(diǎn)對(duì)模型的干擾狀況,影響模型的預(yù)測準(zhǔn)確度。

        故而,該文得到基于模糊近似支持向量回歸的股價(jià)預(yù)測模型,其具體形式如下:

        (5)

        其中,i為所選取股票的第i個(gè)樣本,N為訓(xùn)練樣本總數(shù),xi為第i個(gè)輸入向量,xi=(xi1,xi2,…,xi15)分別代表開盤價(jià)、最高價(jià)、最低價(jià)、總交易量、移動(dòng)平均線、指數(shù)移動(dòng)平均線、異同移動(dòng)平均線、動(dòng)量指標(biāo)、相對(duì)強(qiáng)弱指數(shù)、人氣指標(biāo)、意愿指標(biāo)、成交量比率、多空比率額、真實(shí)波動(dòng)幅度均值和信噪比,yi為實(shí)際股票收盤價(jià)。

        通過構(gòu)造拉格朗日函數(shù)對(duì)式(5)進(jìn)行求解,得到:

        (6)

        分別對(duì)w,b,mi,ni求偏導(dǎo),得:

        (7)

        對(duì)式(7)求解得到:

        (8)

        故上述規(guī)劃問題變?yōu)椋?/p>

        (9)

        f(x)=w*Φ(x)+b=

        (10)

        3.2 隸屬度的獲取

        (11)

        4 實(shí)證分析

        4.1 數(shù)據(jù)的選取

        選擇滬深300成份股的股票日數(shù)據(jù)進(jìn)行股票收盤價(jià)的實(shí)證分析,相應(yīng)時(shí)間序列的日期為2008年1月1日至2019年12月31日,數(shù)據(jù)來源于東方財(cái)富數(shù)據(jù)庫。滬深300成份股共包括300只股票,故該文總共選取300只股票作為總樣本,并從中隨機(jī)選取30只股票進(jìn)行股票收盤價(jià)預(yù)測。所有的實(shí)證均在同一系統(tǒng)環(huán)境下運(yùn)行,系統(tǒng)運(yùn)行環(huán)境為PC(CPU 2.60 GHz,4.00 GB RAM),操作系統(tǒng)為Windows 10,仿真軟件為Matlab R2016b。此外,采用libsvm處理SVR。并將樣本數(shù)據(jù)集分為訓(xùn)練樣本數(shù)據(jù)集和測試樣本數(shù)據(jù)集,選擇樣本數(shù)據(jù)的前80%作為訓(xùn)練樣本,后20%作為測試樣本。

        滬深300成份股是滬深證券交易所于2005年4月8日聯(lián)合發(fā)布的反映A股市場整體走勢的指數(shù),由上海和深圳證券市場中選取300只A股作為樣本,其中滬市有179只,深市121只,綜合反映深交所上市A股的股價(jià)走勢。它具有業(yè)績優(yōu)于整體、對(duì)家方案較優(yōu)、股改行情明顯等優(yōu)勢,經(jīng)常被研究者用于股價(jià)預(yù)測。

        為了更好地反映所提出的基于信噪比的FPSVR股價(jià)預(yù)測模型的有效性,將所選取的股票時(shí)間序列按照股市的波動(dòng)情況(牛市、震蕩市和熊市)分為三個(gè)階段,進(jìn)行階段式預(yù)測。為了更清晰地展示股市狀況,以招商銀行為例,對(duì)階段式預(yù)測進(jìn)行說明,如圖1所示。

        圖1 招商銀行2008-2019年日收盤價(jià)的變化趨勢

        從圖1可看出,招商銀行的收盤價(jià)變化趨勢可分為三個(gè)階段,其中H代表熊市階段,為階段1;N代表牛市階段,為階段2;其余時(shí)間的指數(shù)變化趨勢大致代表震蕩市階段,為階段3。

        本研究為確保所有的輸入特征位于相同參數(shù)范圍內(nèi),以防止大范圍的輸入特征壓倒其他輸入特征,采用式(12)對(duì)數(shù)據(jù)進(jìn)行標(biāo)準(zhǔn)化處理:

        (12)

        其中,xij為第i個(gè)樣本的第j個(gè)輸入特征值,xminj為樣本數(shù)據(jù)第j個(gè)輸入特征的最小值,xmaxj為樣本數(shù)據(jù)第j個(gè)輸入特征的最大值。

        4.2 相關(guān)性分析

        Spearman相關(guān)系數(shù)不但可以衡量存在非線性關(guān)系的相關(guān)變量之間的相關(guān)程度,而且未對(duì)數(shù)據(jù)有嚴(yán)格的假設(shè)要求,故本節(jié)采用Spearman相關(guān)系數(shù)進(jìn)行輸入特征變量與股票收盤價(jià)之間的相關(guān)性分析。為了更清晰地說明特征變量與股票收盤價(jià)之間的相關(guān)性,進(jìn)一步地以招商銀行為例展示輸入特征與股票收盤價(jià)的相關(guān)性分析結(jié)果,如表1所示。

        由表1可看出:

        (1)在置信度為0.01時(shí),開盤價(jià)、最高價(jià)、最低價(jià)、MA和EMA都與股票收盤價(jià)之間的Spearman相關(guān)系數(shù)>0.95,說明這些指標(biāo)都與股票收盤價(jià)存在很強(qiáng)的正相關(guān)性,這些輸入指標(biāo)數(shù)值增加時(shí),股票收盤價(jià)也增加。

        表1 各指標(biāo)相關(guān)性分析

        (2)在置信度為0.01時(shí),SNR與股票收盤價(jià)的Spearman相關(guān)系數(shù)>0.6,說明SVR與股票收盤價(jià)之間存在中等程度的正相關(guān)性,當(dāng)SNR增大時(shí),股票收盤價(jià)隨之增大。

        (3)不管是在置信度為0.01時(shí),還是在置信度為0.05時(shí),MACD、BR、VR和OBV都與股票收盤價(jià)存在非相關(guān)關(guān)系,但是這些指標(biāo)與其他指標(biāo)之間存在相關(guān)關(guān)系,說明這些指標(biāo)可以影響其他指標(biāo),從而間接影響股票收盤價(jià)的變化。

        從上述發(fā)現(xiàn)中可得到以下結(jié)論:①SNR與預(yù)測變量股票收盤價(jià)之間存在直接的正相關(guān)性,進(jìn)而說明所構(gòu)建的SNR特征變量是有效的;②這些輸入特征之間存在緊密的相關(guān)性,且都與預(yù)測變量股票收盤價(jià)之間存在直接或間接的相關(guān)關(guān)系。

        4.3 評(píng)估準(zhǔn)則

        眾所周知,誤差是指測量值與真實(shí)值或?qū)嶋H值之間的差值,主要用來衡量測量結(jié)果的準(zhǔn)確度,其中平均絕對(duì)百分比誤差因能避免誤差相互抵消的問題,可以準(zhǔn)確反映實(shí)際預(yù)測誤差的大小,均方根誤差不僅能對(duì)一組測量中的特大或特小誤差反映非常敏感,而且能夠很好地反映出測量的精密度,故經(jīng)常被用于預(yù)測模型的準(zhǔn)確度測量。該文將平均絕對(duì)百分比誤差(MAPE)和均方根誤差(RMSE)作為股價(jià)預(yù)測模型的評(píng)估準(zhǔn)則,其表達(dá)式如下所示:

        (13)

        (14)

        考慮到上述實(shí)驗(yàn)條件,對(duì)于所有RMSE和MAPE,它們的值越小,所構(gòu)建的股票收盤價(jià)預(yù)測模型的預(yù)測性能越好。

        4.4 結(jié)果分析

        為了驗(yàn)證所提出的基于信噪比的FPSVR股價(jià)預(yù)測模型的有效性,在收盤價(jià)預(yù)測的過程中,選用三種基準(zhǔn)模型預(yù)測方法,分別為模型1:未加入信噪比的支持向量回歸(未加入SNR的SVR);模型2:加入信噪比的支持向量回歸(加入SNR的SVR);模型3:未加入信噪比的FPSVR,并將這三種預(yù)測方法與文中所提的模型4—基于信噪比的FPSVR模型進(jìn)行對(duì)比分析。

        由于金融時(shí)間序列的動(dòng)態(tài)特性是非線性的,文中將FPSVR模型應(yīng)用于收盤價(jià)預(yù)測時(shí),使用高斯核函數(shù)作為核函數(shù),因?yàn)楦咚购撕瘮?shù)在一般的平滑假設(shè)下往往會(huì)有很好的性能[27]。此外,本節(jié)采用基于十折交叉驗(yàn)證的網(wǎng)格搜索方法,對(duì)傳統(tǒng)的支持向量回歸參數(shù)進(jìn)行了選擇,并對(duì)FWSVR模型也使用同樣的參數(shù),所使用的參數(shù)如表2所示。

        表2 基于信噪比的FPSVR模型所使用的最優(yōu)參數(shù)

        此外,采用遺傳算法(GA)對(duì)FPSVR股價(jià)預(yù)測模型進(jìn)行求解,GA中所使用的相關(guān)參數(shù)為G=2 000,popsize=50,其中G為進(jìn)化代數(shù),popsize為種群規(guī)模,相應(yīng)的適應(yīng)度函數(shù)變化情況如圖2所示。

        圖2 GA的適應(yīng)度函數(shù)的變化情況

        從圖2中可看出,在進(jìn)化代數(shù)為2 000時(shí),隨著GA算法的進(jìn)化代數(shù)增加,適應(yīng)度函數(shù)趨于穩(wěn)定,從而說明這種情況下,GA算法可以很好地對(duì)FPSVR進(jìn)行求解。

        下面根據(jù)上述所建立的四種不同模型對(duì)滬深300成份股中隨機(jī)選擇的30只股票時(shí)間序列進(jìn)行收盤價(jià)預(yù)測,并給出了四種模型在不同階段的預(yù)測誤差,如表3所示。

        從上述實(shí)驗(yàn)結(jié)果中可以看出:

        (1)與模型1和模型2相比,模型3和模型4的預(yù)測誤差更低,說明加入模糊隸屬度和雙邊權(quán)重的FPSVR模型可以更好地實(shí)現(xiàn)股價(jià)預(yù)測。

        (2)模型2和模型4的預(yù)測誤差分別低于模型1和模型3的預(yù)測誤差,從而表明加入信噪比特征變量后的股價(jià)預(yù)測模型準(zhǔn)確度更高。

        (3)分階段來看,不管是SVR股價(jià)預(yù)測模型還是FPSVR的股價(jià)預(yù)測模型,在階段1和階段2的預(yù)測誤差均較高于階段3的預(yù)測誤差,說明在震蕩時(shí)期,這兩種模型更適用于股價(jià)預(yù)測。

        (4)分階段來看,在階段1和階段2時(shí)期,F(xiàn)PSVR的股價(jià)預(yù)測模型的預(yù)測誤差要遠(yuǎn)低于SVR股價(jià)預(yù)測模型的預(yù)測誤差,從而表明,該文所構(gòu)建的FPSVR模型可以彌補(bǔ)SVR模型在階段1和階段2時(shí)期的股價(jià)預(yù)測誤差較大的不足,可以更好地實(shí)現(xiàn)股價(jià)預(yù)測。

        表3 模型對(duì)比分析結(jié)果

        5 結(jié)束語

        為了更好地實(shí)現(xiàn)股價(jià)的精準(zhǔn)預(yù)測,提出了基于信噪比的模糊近似支持向量回歸(FPSVR)的股價(jià)預(yù)測模型,并通過實(shí)證研究、相關(guān)性分析和對(duì)比分析三個(gè)方面實(shí)現(xiàn)和驗(yàn)證所提模型的準(zhǔn)確性。研究結(jié)果表明,與現(xiàn)有模型相比,該模型不僅在震蕩期,而且在牛市和熊市期均可實(shí)現(xiàn)股價(jià)的精準(zhǔn)預(yù)測。該模型通過對(duì)滬深成份股2008年至2019年的股票日數(shù)據(jù)進(jìn)行研究和分析,實(shí)現(xiàn)對(duì)股票的準(zhǔn)確預(yù)測,為投資者提供一種更有效的預(yù)測方法,從而有助于投資者的投資決策。

        猜你喜歡
        模型
        一半模型
        一種去中心化的域名服務(wù)本地化模型
        適用于BDS-3 PPP的隨機(jī)模型
        提煉模型 突破難點(diǎn)
        函數(shù)模型及應(yīng)用
        p150Glued在帕金森病模型中的表達(dá)及分布
        函數(shù)模型及應(yīng)用
        重要模型『一線三等角』
        重尾非線性自回歸模型自加權(quán)M-估計(jì)的漸近分布
        3D打印中的模型分割與打包
        美女国产毛片a区内射| 国产精品一区二区午夜久久| 日日噜噜夜夜久久密挑| 日本最新一区二区三区在线视频| 亚洲av无码片vr一区二区三区| 野外性史欧美k8播放| 天堂在线观看av一区二区三区| 久久久噜噜噜噜久久熟女m| 嫩草伊人久久精品少妇av| 国产成人综合久久亚洲精品 | 亚洲日韩中文字幕在线播放| 日本高清色一区二区三区| 日本在线观看一区二区三| 免费中文熟妇在线影片| 中文字幕久久久久人妻无码| 久久亚洲一区二区三区四区五| 日本道色综合久久影院| 91久久综合精品久久久综合 | 蜜臀性色av免费| 免费特级黄毛片| 亚洲精品国产熟女久久| 国产av天堂亚洲av刚刚碰| 日本黄网站三级三级三级| 日本午夜福利| 在线看亚洲十八禁网站| 青青草手机视频免费在线播放| 国产av无码专区亚洲avjulia| 亚洲爆乳无码专区| 国产三级精品三级男人的天堂| 色婷婷综合久久久中文字幕| 国产欧美日韩a片免费软件| 亚洲精品无码不卡在线播放he| 99热这里只有精品3| 97国产精品麻豆性色| 手机免费在线观看av网址| 亚洲中国精品精华液| 老师翘臀高潮流白浆| 仙女白丝jk小脚夹得我好爽| 成人免费播放视频影院| 中文字幕在线日亚洲9| 色综合久久丁香婷婷|