999精品在线视频,手机成人午夜在线视频,久久不卡国产精品无码,中日无码在线观看,成人av手机在线观看,日韩精品亚洲一区中文字幕,亚洲av无码人妻,四虎国产在线观看 ?

特征房價空間分析及連續型深度置信網絡預測

2021-07-30 03:29:34吳莞姝胡龍超趙凱
華僑大學學報(自然科學版) 2021年4期
關鍵詞:深度特征模型

吳莞姝, 胡龍超, 趙凱

(1. 華僑大學 建筑學院, 福建 廈門 361021;2. 華僑大學 數量經濟研究院, 福建 廈門 361021)

房地產行業在國民經濟中有著舉足輕重的地位,作為房地產價值體現的房價不僅關系到國民經濟的健康、平穩發展,其漲跌變動還關系到居民的財富及生活水平[1-2].影響房價的因素眾多,不僅受建筑特征的影響,還會因地區之間公共服務設施的資源配置不均衡而產生差異.人們根據對公共服務和設施的偏好選擇居住區域,特別是在物質生活水平日漸豐富的經濟社會環境下,居民往往愿意支付較高的價格以獲得優質的公共服務,而這部分優質公共資源的價值就資本化于房價之中[3].

針對房價的研究多基于特征價格模型.Ridker等[4]采用特征價格法分析空氣污染對房價的影響.Xiao等[5]利用特征向量空間濾波方法消除空間自相關性后,發現北京周圍的設施對北京房價的影響參差不齊.張驥[6]以北京市二手市場上的商品住宅和非商品住宅為研究對象,利用基于特征價格的配對回歸模型,研究北京學區房交易價格,發現北京市的學區房溢價已高出24.3%.文獻[7-8]以北京為例,通過特征價格模型探尋房價影響因素,研究證實地鐵、公交等交通基礎設施及優質教育資源、高水平的醫療機構、公園等公共服務設施對房價上漲皆具有明顯的正效應.Li等[9]則通過整合鏈家網站、Mobike網站及百度地圖興趣點的公開數據,分析上海公寓價格的空間模式及其與當地配套屬性的關聯,研究發現公園、學校、醫院和銀行等公共服務設施及娛樂、購物等私人服務設施推高了市中心地區的房價.

利用特征價格模型構建房價分析模型需要“先驗”地設定函數形式,這往往容易損失房價與其特征變量之間的深層次關系.近年來,國內外學者嘗試應用多種機器學習模型探討房價的變化趨勢和影響因素等問題.申瑞娜等[10]結合主成分分析和支持向量機,綜合考察影響上海住房價格的8種因素,并對上海房價進行預測.文獻[11-12]利用灰度GM(1,1)預測模型分別對福州市和周口市的房價走勢進行預測,并得到精準度較高的預測結果.張智鵬等[13]利用梯度提升樹(GBDT)算法對房價進行預測,實驗結果表明,公共設施、生活服務、學校、購物服務等是對房價產生明顯影響的因素.這些文獻均采用結構較為簡單的機器學習模型,并且分析數據的特征維度偏低.傳統機器學習方法難以全面且精確地挖掘特征因素和房價之間的聯系.

淺層BP神經網絡模型(BPNN)在預測上優于傳統機器學習模型,但仍存在學習速度慢、易陷入局部收斂等問題.而深度置信網絡采用無監督訓練方式,具有較好的降維性能.一方面,深度置信網絡能有效克服傳統人工神經網絡需要大量有監督信號和易陷入局部極小等缺點;另一方面,深度置信網絡可高效處理高維的數據并挖掘變量之間的深層關系.此外,深度置信網絡解決了大規模數據計算耗時問題且精度較高,并成功應用于多種人工智能問題的研究,尤其在圖像處理、聲音辨識和智能網絡分析等方面的應用中成效顯著[14-15].由于深度置信網絡在運算時使用數據離散化方法進行特征提取,隱藏層和可見層節點均為伯努利值(0或1),這使深度置信網絡不適用于對連續型變量的高精度預測[16-17].為提高模型的預測精度,學者們將深度置信網絡進行改進,使其能夠有效處理連續型的輸入變量[18-19].

盡管深度置信網絡在人工智能領域特別是模式識別任務中取得了較好的成果,但將其應用于現實經濟問題的研究仍較少見.基于此,本文嘗試將連續型深度置信網絡擴展至房價問題的研究中,依據特征房價理論并考慮到上海市二手房交易價格可能存在的空間相關性,構建空間計量模型以分析各特征變量對二手房交易價格的影響.在此基礎上,利用連續型深度置信網絡建立房價與多維影響因素之間的深度學習預測模型,深層挖掘其潛在規律.

1 研究介紹及數據探索

1.1 研究區域界定

由于上海是我國“超一線”城市,房地產市場發展較為成熟和完善,具有一定的代表性;同時,上海浦東和浦西在城市化建設和房屋價格上具有明顯的差異,這為探討建筑特征、區位特征和鄰里特征對二手房交易價格的影響效果提供了較好的素材,因此,選擇上海市作為研究區域.

1.2 數據來源

房屋交易數據源于“鏈家二手房交易平臺”(https:∥m.lianjia.com),鏈家的樓盤數據庫管理著160多個城市1.1億套真實的房產數據,依托互聯網對數據進行標準化管理,實現信息的無差別共享.

基于Python語言的爬蟲技術,按照不同行政區對鏈家上海市二手房交易平臺上的數據進行收集.利用Beautiful Soup對網頁返回結果進行重構,得到超文本標記語言(HTML)的樹狀結構,再使用正則表達式對所需信息進行提取,進而獲取變量數據.最終,所爬取的數據涉及房屋交易額、交易單價、百度地理坐標(BD08)、房屋戶型、所在樓層、建筑面積、戶型結構、建筑類型、建成年代、裝修情況、梯戶情況等變量信息,共計45 131條原始信息.同時,搜集整理上海市全部40家三甲醫院(及其分院)、上海市34所重點中小學及上海市所有地鐵站出口的地理坐標數據.上海市的主要三甲醫院及重點中小學的地理位置信息,如表1所示.

表1 上海市的主要三甲醫院及重點中小學的地理位置信息

1.3 數據清洗

上海市各行政區的原始數據交易時間跨度不一,為剔除時間跨度的影響,對各行政區數據進行切割,保留共同交易時間跨度的數據為研究樣本.另外,由于在房屋的交易權屬中個稅的收取不同,因而,刪除交易權屬中的售后公房.為避免極端價格對數據分析的影響,房屋用途中刪除別墅、車庫和商業辦公類房屋,只保留普通住宅.剔除含有缺失值的數據條目,最終獲得9 058個樣本,涉及房屋戶型、建筑面積、建成年代、所在樓層、裝修情況、配備電梯等變量.

房屋戶型變量的形式以字符數字組合為主,采用正則表達式將臥室、大廳、廚房和衛生間這幾個數值提取出來,并分別作為建筑特征的變量.建筑面積的原始數據中帶有面積單位m2,利用正則表達式剔除該單位,并把面積值變為浮點型數據.原始數據中建成年代為房屋建成年代,采用爬取數據的年份(2019年)減去建成年代的方式,計算房屋建成年數.所在樓層的原始數據為高樓層、中樓層和低樓層.對所在樓層進行數值化處理,把高樓層、中樓層和低樓層分別賦值為2,1和0.裝修情況和配備電梯為二值數據,裝修情況為已裝修或未裝修,配備電梯為有或無,利用1和0進行數值化處理.

鏈家官網的坐標數據來源于百度地圖,利用ArcGIS軟件將清洗后樣本數據的地理坐標轉化為WGS84坐標,其分布情況如圖1所示.

圖1 樣本數據的空間分布

1.4 特征變量的選取

借鑒以往研究,將二手房的特征梳理為建筑特征、區位特征和鄰里特征3類.二手房建筑特征為房屋本身的屬性,涉及的變量包括房屋戶型、所在樓層、建筑面積、建成年代、裝修情況及配備電梯.區位特征量化了二手房區位對整個城市的可達性,如出行成本等.

將二手房到城市中心的距離作為二手房的區位特征變量;以陸家嘴金融貿易中心區域的質心為城市中心點.鄰里特征通常指房屋周圍的環境及配套,如交通站點、學校、醫院等.選取到最近三甲醫院的距離、到最近重點中小學的距離及到最近地鐵站的距離體現鄰里特征.距離計算方式取大圓距離,大圓距離是將地球看作一個球形,計算球面上兩點的最短路徑.特征變量的相關說明,如表2所示.

表2 特征變量說明

1.5 描述性統計

清洗后數據的描述性統計,如表3所示.由表3可知:二手房成交單價最大值為147 668元·m-2,最小值為7 059元·m-2,均值為49 577元·m-2,偏度為1.011 452,屬于右偏數據,大多數交易價格集中在50 000元·m-2左右.

表3 數據的描述性統計

繪制QQPlot分布圖,將數據分布與正態分布進行對比,結果如圖2所示.圖2中:S為標準正態值;Q代表數據的分位數.由圖2可知:對數成交單價數據近似服從于正態分布.

圖2 對數化后的二手房成交單價分布

化后的利用空間趨勢分析將二手房成交單價投影到XZ和YZ平面上,繪制上海市二手房交易價格分布的空間趨勢,如圖3所示.圖3中:X為經度;Y為緯度;Z為單價.由圖3可知:不論在東西方向還是南北方向,上海市二手房交易單價都呈現由中心向兩頭遞減的趨勢.綜上可知,上海市二手房市場具有中心高價的特點.

圖3 二手房成交單價的空間趨勢

2 特征房價空間分析

2.1 空間自相關

(1)

由于二手房交易價格數據在空間上是點要素的形式,沒有多邊形的拓撲關系,在空間上的分布也較不均衡,故整體的空間關聯程度可以利用全局莫蘭指數判斷,莫蘭指數(I)的表達式為

(2)

計算得到全局空間自相關的檢驗結果如下:莫蘭指數I為0.236 862 3;統計量為-0.000 271;p近似為0.由檢驗結果可知:二手房交易價格的莫蘭指數顯著為正,說明上海市二手房交易價格具有空間集聚效應.

上海市二手房交易價格全局空間自相關散點圖,如圖4所示.通過莫蘭散點圖將空間自相關分為高-高集聚、高-低集聚、低-高集聚、低-低集聚這4種類型.圖4中:L為lnP的空間一階滯后;第1,3象限是高-高集聚、低-低集聚區域,即同質化明顯的區域;而第2,4象限是高-低集聚、低-高集聚區域,即異質性較強的區域.

圖4 上海市二手房交易價格全局空間自相關散點圖

由圖4可知:絕大多數樣本落入第1,3象限,少部分樣本落入第2象限,空間集聚特點較為明顯.

2.2 空間異質性

借助ArcGIS軟件繪制上海市二手房交易價格的LISA集聚狀況,如圖5所示.由圖5可知:中心城區的房價呈現高-高集聚的空間效應,且越靠近城市中心點,高-高集聚的特征越顯著;高-低集聚區域沿著高-高集聚區域的邊緣分布;而低-低集聚效應區域大多分布在上海周邊地區.

圖5 LISA集聚狀況 圖6 冷熱點分布

進一步,通過局部空間自相關檢驗探討分析上海市二手房交易價格的空間異質性.局部空間自相關水平的冷熱點分布,如圖6所示.由圖6可知:熱點區域和冷點區域均在99%的置信水平上顯著.上海市二手房交易價格呈現“中間高、四周低”的空間格局,相較于LISA集聚,冷熱點分布更寬,涉及更多邊緣樣本.城市中部的浦西七區、寶山區及閔行區的二手房交易價格為高-高集聚,環繞四周的嘉定區、青浦區、松江區、奉賢區和浦東新區外環城區的二手房交易價格為低-低集聚.

2.3 空間計量模型估計

上海市二手房交易價格具有空間關聯性,應選擇空間計量模型進行分析.拉格朗日乘數檢驗項L-Mlag,LMerr及其穩健值R-LMlag,R-LMerr的檢驗結果,如表4所示.由于LMlag,LMerr均顯著,需進一步比較R-LMlag和R-LMerr的顯著性,又因為R-LMerr顯著而R-LMlag不顯著,故選擇空間誤差模型(SEM)進行分析.

表4 拉格朗日乘數的檢驗結果

基于特征價格法建立SEM,探討影響上海市二手房交易價格的可能因素.SEM回歸結果,如表5所示.表5中:ES為標準誤差;λ為空間自相關系數;*,**,***分別表示在10%,5%,1%水平上影響有統計學意義.對數似然值為1 004.224;赤池信息準則AIC為-1 974.4.

由表5可知:除廚房數量外,其他特征變量對二手房交易價格的影響皆有統計學意義;已裝修、帶電梯、有客廳且洗手間數量較多的二手房交易價更高;臨近重點中小學、醫院和市中心的二手房交易價格較高;然而,臥室數量及建筑面積在一定程度上會對二手房交易價格產生一定的抑制作用,原因可能是上海市過高的單價抑制了人們對大面積住宅的需求;樓齡與樓層均在1%的顯著性水平下對房價有反向影響,但系數較小.

表5 SEM回歸結果

3 連續型深度置信網絡房價預測

3.1 連續型深度置信網絡簡介

連續型深度置信網絡(CDBNN)改造于深度置信網絡(DBN).DBN是由多個受限玻爾茲曼機 (RBM)逐層堆疊而成,其核心思想是自底向上每一層RBM對輸入數據進行提取、抽象,盡可能保留重要信息,訓練過程一般采用貪婪無監督方式,即逐層對DBN中的每一個RBM進行訓練.

RBM是一種基于能量的概率生成模型,生成模型是對特征和標簽之間的聯合分布進行建模.當可見層的狀態ν和隱藏層的狀態h確定后,RBM模型中的能量可以表示為

(3)

式(3)中:θ為參數向量;ai和bj分別為可見層第i個神經元上的偏置和隱藏層第j個神經元上的偏置;ωi,j為可見層神經元和隱藏層神經元之間的連接權重值.

基于能量函數,可得ν和h的聯合概率分布為

上式中:z(θ)為歸一化函數,使得概率之和為1.

依據聯合概率分布,可以得到在可見層狀態ν確定時,隱藏層每個神經元被激活的概率,以及在隱藏層狀態h確定時,可見層每個神經元被激活的概率分別為

(4)

(5)

一般采用梯度下降方法求取最優參數值,過程中涉及難以求解的歸一化函數z(θ),常用吉布斯(Gibbs)采樣方法近似計算[20].CDBNN是在DBN的基礎上改進,對式(4),(5)和激活函數σ(x)進行改進,使其適用于連續型數據,即

(6)

(7)

(8)

式(6)~(8)中:Ni(0,1),Nj(0,1)表示均值為0且方差為1的高斯隨機變量;φ為常量;θH和θL為漸近線,一般取樣本中的最大值和最小值.

由于連續型深度置信網絡是在深度置信網絡的基礎上衍生而來,因此,該方法同樣采用誤差反向傳播的算法進行網絡調優.CDBNN算法主要有以下8個步驟.

步驟1準備訓練數據D=(x1,x2,…,xn),共n個樣本,假設所有神經元的狀態使用狀態集{Si}表示,隨機初始化所有神經元的參數,設訓練的最大次數為K次.

步驟4根據步驟3所得的Si,同步驟2,計算隱藏層的重構神經元狀態Sj.

步驟5繼續隨機選擇下一個訓練樣本,返回步驟2,如果樣本集中的樣本都選完畢,則依據式(8)計算參數變化量,更新方式為wi,j(k+1)=wi,j(k)+Δw,ai(k+1)=ai(k)+Δai.

步驟6進行第k+1次訓練,當權重的變化量落入預定的范圍內,即|Δwi,j|<ε,其中,ε是預先設定的誤差范圍,或者訓練次數達到k次,則訓練停止.

步驟7將訓練好的RBM的輸出作為下一層RBM的輸入層輸入數據,按照步驟1~6進行訓練,直到訓練完DBN的所有RBM層.

步驟8網絡調優:完成DBN的訓練后,需進一步優化深度神經網絡權值.將訓練好的DBN網絡作為網絡的初始狀態,訓練得出的參數作為DBN的初始參數;然后,使用反向傳播的方法,運用梯度下降法對網絡的整體權值進行有監督的學習.

3.2 連續型深度置信網絡結構的確定

連續型深度置信網絡結構的確定實質上就是選擇深度置信網的超參數.待確定的神經網絡結構的超參數包括神經網絡的層數、神經網絡隱藏層的節點數、學習率的確定、高斯隨機變量中的方差值和樣本迭代次數的選擇及其他參數的選擇.

超參數調優即選擇超參數使網絡結構達到最優的效果,是訓練神經網絡的核心任務.目前,常用的超參數調優方法有網格搜索與隨機搜索.前者基于整個超參數空間進行搜索,速度較慢,但可獲得最優的超參數組合.后者速度快,但可能會錯過搜索空間中最優的超參數值.借鑒Snoek等[21]的思路,利用貝葉斯思想自動優化超參數,不僅能有效兼顧上述兩種方法的優點,還能借助Python的hyperopt模塊輕松實現優化超參數.主要超參數的估計值,如表6所示.

表6 主要超參數的估計值

3.3 預測模型性能的比較

在建立連續型深度置信網絡的過程中,將所有樣本按7∶3的比例隨機分成訓練集和測試集,先通過訓練集對模型進行訓練,再使用訓練后的模型對測試集數據進行預測.將文中的預測結果與現有文獻采用的支持向量機(SVM)、集成模型(采用Adaboost算法)和BP神經網絡模型的預測結果進行對比.連續型深度置信網絡、SVM、集成模型、BP神經網絡的預測誤差分別為0.006 67,0.007 61,0.008 42,0.029 03.BP神經網絡的預測誤差遠高于其他3個模型,這是由于隨機初始化使其難以達到全局最優值.而連續型深度置信網絡可預先對BP神經網絡進行預處理,有效緩解隨機初始化對最優預測的阻礙.此外,CDBNN的預測結果也略優于SVM和集成模型,表明CDBNN有更高的復雜度,能夠更加深入且全面地進行特征分析.

4種模型測試集樣本點的預測殘差絕對值,如圖7所示.圖7中:ε為殘差絕對值.由圖7可知:與其他模型相比,CDBNN的預測殘差總體情況更優,CDBNN能夠有效地解決BP神經網絡在預測模型上存在的不足.

圖7 測試集樣本點的預測殘差絕對值 圖8 殘差絕對值與二手房交易價格的關系

3.4 預測結果與討論

通過繪制殘差絕對值,可對CDBNN模型的預測結果進行評價.殘差絕對值與二手房交易價格的關系,如圖8所示.由圖8可知:對于房價偏低或偏高的區域,影響上海市二手房交易價格的因素較為復雜,不僅局限于房屋建筑特征變量;房價偏低的區域大部分偏離市中心,距上海市重點中小學、三甲醫院以及地鐵站距離較遠,利用特征價格法選取的變量對房價偏低區域的房價預測能力相對較差.這些區域的二手房交易價格可能會更多地受到其所在區域的亞中心及該區域所配套的基礎設施的影響.對于二手房交易價格偏高的區域,預測的殘差絕對值相對較大,可以認為當房價過高(P≥80 000元·m-2)時,二手房交易價格的影響因素更加復雜.這其中除特征變量之外,還可能與購房者的購房目的等因素有關.對上海市高房價區域的購買者來說,房價彈性相對較低,他們對高房價并不敏感;高房價的購買者對房屋的消費不僅在于其本身的價值,而可能是出于政策便利性和高房價周圍的鄰里交際環境.

各行政區預測結果的殘差絕對值平方,如圖9所示.由圖9可知:崇明區、黃浦區和靜安區的預測結果的殘差絕對值較大;崇明區的二手房交易價格偏低,而交易價格偏高的區域大多集中在黃浦區和靜安區這兩個市中心區域;浦東新區預測結果的殘差絕對值較小,基于特征房價探討的變量對浦東新區這樣的非市中心的二手房交易價格的預測效果較好;楊浦區和浦東新區的預測結果的殘差絕對值很相近,而靜安區和黃埔區的殘差絕對值相對較大.這可能是因為楊浦區和浦東新區隔海相望,楊浦區的經濟和浦東新區的經濟相互影響較大,黃浦區和靜安區作為上海一直以來的市中心,其二手房交易價格的影響因素較為復雜;而浦東新區是改革開放后繁榮的區域,受上個世紀90年代開放的房地產市場影響較大,所以,預測效果較好,特征價格法所選的建筑特征、鄰里特征和區位特征對新區房價的解釋力度更強;對于黃浦區和靜安區這樣的老中心區域,其房價的解釋力度相對較小;黃浦區和靜安區的二手房交易市場的影響因素已經超出特征價格變量的解釋范圍.

圖9 各行政區預測結果的殘差絕對值平方

4 結論

以上海市二手房交易市場為例,通過空間自相關分析,發現上海市二手房交易單價在空間上具有顯著的自相關效應.二手房交易價格在上海市核心區域存在高-高集聚效應,在周邊區域呈現低-低集聚效應,而在核心與周邊交界地區存在高-低集聚和低-高集聚的負向空間效應.與此同時,基于連續型深度置信網絡對特征二手房交易價格進行分析預測,發現特征變量對價格偏高區域的二手房交易價格解釋力度較小,價格偏高區域的二手房交易價格影響因素較為復雜.從區域角度分析,除中心區域外,基于深度置信網絡的特征變量對上海市二手房交易價格預測能力良好.連續型深度置信網絡不僅能有效地解決BP神經網絡在預測模型上存在的不足,而且與其他機器學習模型相比,連續型深度置信網絡能更精準地對房價進行預測,從而為政府部門進行房價預測提供理論支持和政策導引.

文中采用一種能夠處理大數據的深度學習模型,但由于獲取數據的難度大,僅選取十余個解釋變量.現實中,影響二手房交易價格的因素非常多,如加入更多的解釋變量,基于連續型深度置信網絡對特征二手房交易價格模型的預測將會更加精準.采用空間分析及深度學習技術對二手房交易價格進行研究,在各大城市均具有普遍適用性,可應用于其他城市的房價研究.

文中研究結果可為后續相關研究提供方法參考和模型借鑒.通過對上海市不同行政區及不同價格區間的房價預測模型效果進行差異性分析發現,在價格偏高的區域和上海市中心區域的預測效果較差.這為后續相關研究提供兩方面借鑒:一方面,在預測房價時,需要考慮到空間異質性的影響,應針對不同區域構建不同的預測模型;另一方面,為進一步提高房價預測的精度,需要在建筑、區位、鄰里等特征變量的基礎上納入更多相關的社會經濟要素,從而提升模型的預測能力.對房價走勢進行高精度預判,有助于政府制定調控政策.房地產市場調控一直是政府相關部門的工作重點,而穩定房價是調控的主要目標.對房價走勢進行高精度預測具有一定的現實意義,可為政府相關部門完善房地產市場、優化城市規劃設計提供一定的理論支持.

猜你喜歡
深度特征模型
一半模型
深度理解一元一次方程
重要模型『一線三等角』
重尾非線性自回歸模型自加權M-估計的漸近分布
如何表達“特征”
不忠誠的四個特征
當代陜西(2019年10期)2019-06-03 10:12:04
深度觀察
深度觀察
深度觀察
抓住特征巧觀察
主站蜘蛛池模板: 日韩av手机在线| 亚洲视频二| 国产理论一区| 久久综合结合久久狠狠狠97色| 亚洲成人网在线观看| 午夜精品久久久久久久无码软件| 小蝌蚪亚洲精品国产| 国产免费人成视频网| 99视频在线看| 中文字幕首页系列人妻| 狠狠色香婷婷久久亚洲精品| 免费在线视频a| 嫩草国产在线| 亚洲福利片无码最新在线播放 | 欧美翘臀一区二区三区| 久久精品中文字幕免费| 制服丝袜在线视频香蕉| 亚洲日韩精品综合在线一区二区| 国产成人综合久久精品下载| 亚洲乱亚洲乱妇24p| 少妇精品在线| 凹凸国产熟女精品视频| 久久综合丝袜长腿丝袜| 精品视频免费在线| 欧美无遮挡国产欧美另类| 亚洲美女AV免费一区| 91色综合综合热五月激情| 久久这里只有精品免费| 国产日韩精品一区在线不卡| 99精品高清在线播放| 久久久久国产一级毛片高清板| 91成人在线免费视频| 久久永久免费人妻精品| 欧美精品1区| 欧美一级夜夜爽www| 亚洲综合在线网| 日韩欧美国产区| 日韩欧美国产另类| 国产91av在线| 久久青草视频| www.亚洲一区| 精品天海翼一区二区| 日本高清免费一本在线观看| 国产成人av一区二区三区| 亚洲丝袜中文字幕| 91精品啪在线观看国产91| 久无码久无码av无码| 国产一区免费在线观看| 国内熟女少妇一线天| 9cao视频精品| 国产精品一线天| 国国产a国产片免费麻豆| 欧美日韩午夜视频在线观看| 91毛片网| 久久久91人妻无码精品蜜桃HD| 欧美精品v| 成人永久免费A∨一级在线播放| 香蕉伊思人视频| 久久久噜噜噜| 国产精品中文免费福利| 欧美激情第一欧美在线| 中文字幕无码制服中字| 午夜福利亚洲精品| 国产99欧美精品久久精品久久| 午夜无码一区二区三区在线app| 成人毛片在线播放| 欧美激情视频一区二区三区免费| 一级毛片网| 婷五月综合| 亚洲成在人线av品善网好看| 亚洲精品少妇熟女| 熟女视频91| 久久亚洲AⅤ无码精品午夜麻豆| 欧美劲爆第一页| 久久亚洲天堂| 福利国产微拍广场一区视频在线| 狠狠ⅴ日韩v欧美v天堂| 中文字幕在线观| 午夜视频免费一区二区在线看| 中文字幕2区| 亚洲婷婷在线视频| 亚洲av无码片一区二区三区|