期刊大全 雜志訂閱 SCI期刊 投稿指導 期刊服務 文秘服務 出版社 登錄/注冊 購物車(0)

首頁 > 精品范文 > 主成分分析論文

主成分分析論文精品(七篇)

時間:2022-09-28 02:20:43

序論:寫作是一種深度的自我表達。它要求我們深入探索自己的思想和情感,挖掘那些隱藏在內心深處的真相,好投稿為您帶來了七篇主成分分析論文范文,愿它們成為您寫作過程中的靈感催化劑,助力您的創作。

主成分分析論文

篇(1)

成分分析法,又稱主分量分析法是指相關的經濟變量間通常存在著起主導作用的決定性因素,通過對原始變量的相關矩陣內部結構進行分析,找出幾個不相關的綜合指標來線性表示原來的變量,主成分之間既互不相關,又盡可能多的包含了原指標集合。這種方法首先由Hotelling提出,其主要思想是降維。Stone(1947)對美國1929-1938年間的17項國民經濟統計指標進行主成分分析,發現完全可以用三個經濟指標來概括原來的17項指標,大大簡化了數據分析。M.Scott(1961)對英國157個城鎮的發展水平進行主成分分析,發現原57個測度指標完全可以由5個綜合變量替代,既解決了原指標間的信息重疊問題,又簡化了原指標體系的指標結構,主成分分析由此推廣。邱東(1990)系統闡述了主成分分析法的定義、基本思想、基本步驟和特點,認為主成分分析法可以消除評價指標間的相關影響,并且伴隨數學變換過程生成信息量權數和系統效應權數,保證了客觀性。同時也指出了主城分析法在計算綜合評價值未充分考慮指標的重要程度等不足,主要適用于被評價對象較多的綜合評價。隨后,眾多學者對此提出了改進:孟生旺(1992)針對原始數據的標準化處理和主成分個數的選擇問題,認為標準化不如均值化的無量綱處理方法,提出了非標準化主成分分析法。陳述云等(1995)通過對原始數據作對數—中心化轉換,用原始變量的非線性組合表示主成分,同時重點分析樣本協方差矩陣而非相關系數矩陣,提出了非線性主成分法。朱泰英等(2004)提出了加權主成分分析法,認為可以將主成分分析法的客觀分析和層次分析法的主觀分析有機結合。王璐等(2006)在對主成分分析法的權數、降維等問題的研究上,提出了首先要按主成分分析法對指標體系進行分類,得到各方面的評價值后再進行主成分分析,最終得到綜合評價值的二重主成分分析法。段力誌等(2009)在傳統主成分分析法基礎上,首先對原始指標值進行預處理,再借助軟件,將原始數據轉化為少數幾個主成分的線性組合,并進行加權變換,得到改進的主成分綜合值。白雪梅等(1995)則分析了“均值化”、“標準化”、“極差正規化”三種方法的選擇條件是保證方差損失最小。陳衍泰等(2004)認為主成分分析法具有全面性、可比性和客觀合理性等優點,比較適合對評價對象進行分類,但需要大量數據,函數意義不夠明顯,不能反映客觀發展水平。蘇為華(2012)提出經典的R型主成分本質是單項指標標準化結果的加權算術平均值,比當量平均法復雜。趙利等(2013)通過主成分分析法對宏觀經濟中影響城鎮勞動就業因素分析時,提出主成分為宏觀經濟和技術進步,通過VAR模型對主成分進行分析,得出宏觀成分中對城鎮勞動就業影響最大的是消費、產業結構和城市化水平,而技術進步成分中影響最大的是技術進步的結論。黃利文(2013)針對主成分分析中存在的未考慮負向因子的影響,以及采用線性加權法時確定權重方法不統一,評價結果非唯一等缺陷,提出了逼近理想點的主成分分析法,更好地反映了原始數據信息,并較為客觀地給出了綜合評價結果。林海明等(2013)認為主成分分析因缺乏應用條件的考慮而導致評價結果不具合理性甚至錯誤,通過分析因子分析法因子載荷陣的簡單結構、加權算術平均數的合理性,得出主成分分析的應用條件是:指標是正向、標準化的;主成分載荷陣達到更好的簡單結構時,主成分正向,且主成分與變量顯著相關。

二、因子分析

法因子分析法是指從被評對象的觀察變量的相關度出發,利用降維的思想,把繁雜的變量盡可能歸納為幾個綜合因子進行分析的的一種多變量統計分析方法。其基本思想是:將觀察變量按相關度的高低或聯系的緊密程度進行分類,類別內部變量相關性高,聯系緊密,而類別之間的變量則相關度較低,聯系稀疏,每一類變量則代表一個公共因子。具體步驟為:

三、逼近理想解的排序法

篇(2)

關鍵詞 期刊評價;區間數據;主成分分析;相關分析

中圖分類號 文獻標識碼 A

Analyzing the Current Development of Periodicals

Based on Interval-Censored Data

LI Jing-bo

(Periodical Agency, Hunan University, Changsha, Hunan 410082, China)

Abstract The article uses all periodicals' data collected in CSSCI database to build a model by the ways of relative analysis and interval-censored data's principal component analysis. It analyzes all periodicals in the database objectively and systematically, and explores its characters, regularity and reasons, which provides some ways to manage academic journal better.

Key words periodicals' comment; interval-censored data; principal component analysis; relative analysis

1 引 言

學術期刊作為科學技術事業的重要組成部分,具有傳播知識、傳播思想、傳播信息的重要功能,與科技創新和進步有著密不可分的關系.另一方面,學術期刊作為衡量一個國家的科技發展水平的尺度,標志著該國家的生產力水平高低及科學文化事業的興旺發達程度,會產生巨大的社會影響和經濟效益.因此,辦好學術期刊是建立我國科學和技術創新體系的重要環節[1].

自20 世紀80 年代以來,中國的學術期刊出現了一個“黃金時代”,數量急劇增長,似乎呈現出一片繁榮的景象.據統計,目前,全世界共有近24萬種期刊,85%在發達國家出版發行;學術性期刊約10 萬種,其中重要科技期刊8 萬種,幾乎全為發達國家所擁有.與發達國家相比較而言,我國所擁有的學術期刊并不多,2008年全國共出版期刊95 49種,平均期印數16 767萬冊,總印數31.05億冊,總印張157.98億印張,其中,學術期刊有6 000余種(哲學、社會科學類2 339種),在整個期刊中所占比例在40%左右[2].

從數量上講,我國僅次于美國,成為世界學術期刊第二大國.雖然從學術期刊和學術論文的絕對數量上來看,我國算得上是一個學術期刊大國,但是我國學術期刊發展仍面臨著發展水平參差不齊、優秀稿源外流等諸多問題.因此,建立科學的評價模型,分析我國學術期刊的發展現狀和特點,以期對我國學術期刊的健康發展發揮有效的導向作用,并促進其管理的標準化、規范化.

期刊評價是文獻計量學研究的重要組成部分,它通過對學術期刊的發展規律和增長趨勢進行量化分析,揭示學科文獻數量在期刊中的分布規律.

評價學術期刊的文獻計量指標多種多樣,各有側重.在評價期刊時,單一使用任何一種指標都可能造成某些方面不同程度的片面性,而選擇許多個指標時,又常因個數太多而增加對問題分析的復雜性,且人為確定指標權重會產生主觀偏差,加上對于指標之間的相關性未給充分的考慮,造成所確定的權重并不是實際計算過程中所體現的真實權重[3].

隨著人們對期刊評價工作重要性認識的不斷加深,評價方法有了很大的進展,出現了諸如綜合評分、專家討論、費用效益分析、層次分析、模糊聚類、模糊評判、灰色統計、熵值法以及改進熵值法,灰色關聯分析、人工神經網絡方法、主成分分析方法等評價方法.這些方法在評價工作中展示了各自的優點,同時也暴露出各自的弱點和缺陷,在此不一一贅述.

現有期刊評價的研究大多聚焦于通過建立指標體系,對某一類學術期刊進行綜合排序,并且所使用的數據樣本規模較小.本文將采用CSSCI數據庫中收錄的全部期刊數據,通過相關分析、區間數據主成分分析等方法建立模型,實現數據降維和可視化,客觀地對整個數據庫中各學科期刊的整體特征進行全面、系統的分析,探尋其中的特征、規律以及原因.

2 數據描述與指標選擇

本文的數據來自于中文社會科學引文索引(Chinese Social Sciences Citation Index,簡稱CSSCI)是由南京大學中國社會科學研究評價中心開發研制的引文數據庫,用來檢索中文人文社會科學領域的論文收錄和被引用情況.CSSCI中收錄期刊均為學術性期刊,以反映我國社會科學研究和科研成果應用為主要方向,充分顯示我國的科研水平,因此利用其數據可以一定程度反映我國學術期刊的整體情況.

2007-2008年,CSSCI來源期刊共收錄528種期刊,CSSCI擴展版收錄期刊的數量為152種,CSSCI來源集刊共收錄86種期刊,共計收錄766種期刊.

CSSCI使用的指標分為兩類:(一類反映情況,包括發文量、基金論文數、發文機構數、篇均參考文獻量、自引率、引用半衰期等6個指標,

另一類屬于引證指標,包括影響因子、即年指標、總被引頻次、自被引率、被引半衰期等5個指標[4].

下面對上述11個指標的含義進行簡要介紹[5].

1)發文量.絕對數量指標,指某一期刊在一定時期內(一年)所刊登的全部論文數,反映了期刊的信息含量.發文量越高,說明期刊的信息含量就越大.

2)基金論文數.

絕對數量指標,指期刊在當年中受基金資助的論文數量可測度期刊在學術交流中的地位,表征論文產出性質的重要指標,可以衡量期刊論文學術質量.基金論文數越多,說明期刊論文的學術質量就越高.

3)發文機構數.絕對數量指標,指期刊當年所中所涉及的不同機構的數量,可測度期刊論文的機構分布情況,衡量期刊科學生產能力,機構分布越廣,說明期刊具有開放性和作者隊伍具有廣泛性.發文機構數越多,說明期刊的影響范圍就越廣.

4)篇均參考文獻量.相對數量指標,指在給定的時間內(一年),期刊中的平均參考文獻量,在數值上等于期刊參考文獻數/期刊論文總數,通??梢苑从称诳招畔⒌哪芰σ约翱茖W交流程度的高低.篇均參考文獻量越多,說明期刊吸收的信息越多.

5)自引率.相對數量指標,指該刊在當年的自引次數占該刊當年參考文獻總數的比例.自引率越高,說明期刊引用論文中的大部分是發表在自身期刊上的,學術交流程度則越低,且若自引率過高,往往說明該期刊有不良的非正常自引現象.

6)引用半衰期.絕對數量指標,指該期刊當年引用類型為期刊的參考文獻中較新的一半是在多長時間中發表的,可測度期刊文獻老化的速度[6].

一般來說,半衰期長的期刊比短的期刊影響更深遠一些.但該指標受學科的內容、性質等因素的制約,比如說發展穩定的學科比發展較快、較活躍的學科長,基礎理論學科的比技術學科的長,歷史悠久的比新興學科的長.因此該指標在判斷期刊學術質量時往往不易解釋和操作.

7)影響因子.相對數量指標,指期刊前2 年的被引次數占前2 年該刊所總數的比例.

表征論文被引強度的重要指標,反映學術影響力.影響因子可克服由于發文量不同所造成的對期刊被引率的偏差,使期刊總被引頻次這一絕對指標變成了相對指標,所以能較好地反映期刊被使用的真實客觀情況,從而可以對不同期刊的引用次數和質量進行比較,以說明其利用率在科學上的重要性.通常影響因子越大,可以認為期刊在科學發展和文獻交流過程中的作用和影響較大,其學術水平也較高.影響因子是期刊評價中最重要的指標之一.

8)即年指標.相對數量指標,指某刊當年的被引次數占該刊當年總數的比例,是表征論文被引速度的主要指標.即年指標越大,說明期刊論文被引的數量多且速度快.

9)總被引頻次.絕對數量指標,指某刊自創刊以來所刊登的全部論文在某一年被其他期刊(包括本刊)引用的總次數,可測度期刊自創刊以來的學術影響力,是從信息反饋的角度評價期刊的基本指標之一,表征期刊學術質量的重要指標[4].

體現了被引用過程的總體規模.總被引頻次也是期刊評價中最重要的指標之一.

10)自被引率.相對數量指標,指某期刊當年被本刊引用的次數占該刊當年被引用總次數的比例.與自引率類似,自被引率越高反而不好.

11)被引半衰期.絕對數量指標,指某期刊在某年被引用的全部論文中較新的一半是在多長時間中發表的.

與引用半衰期類似,被引半衰期在期刊評價中的受各種因素影響大且不易直接判斷好壞.

以上11個指標中,自引率和自被引率兩個指標屬于反映期刊水平的逆向指標,而引用半衰期和被引半衰期兩個指標不能簡單地用數值大小來解釋期刊水平的優劣,故只選取余下的7個指標(發文量、基金論文數、發文機構數、篇均參考文獻量、影響因子、即年指標、總被引頻次)進行以下的建模分析.

3 模型建立及結果分析

本節采用相關分析和區間數據主成分分析等方法建立模型,從期刊計量指標進行篩選和大規模期刊數據的“打包”處理兩方面,評價期刊的發展水平.

3.1 相關分析

相關分析是研究隨機變量之間是否存在某種依存關系,并對具體有依存關系的變量探討其相關方向以及線性相關程度的一種統計方法.

采用SPSS 15.0軟件對7個指標的原始數據進行相關分析,得到見表1.

表1顯示了7個指標之間的相關系數,可以看出,發文量與基金論文數、發文機構數的相關性很強,影響因子與即年指標、總被引頻次中度相關[4].

由于基金論文數、發文機構數和發文量的相關程度非常高,僅用其中一個指標就能夠代表三者所包含的絕大部分信息.而發文量的應用更具普遍性,這里篩掉基金論文數、發文機構數兩個指標,用余下的5個指標繼續對全部期刊進行數據分析.

3.2 區間數據主成分分析

區間數據是一種非常重要的數據概念.與傳統的分析方法不同,區間數據分析研究的數據表單元不再是一般意義下的定量和定性數值,而是一個實數域區間.經打包處理后的數據集合在幾何上表現為一個高維的超矩形,這種處理方式不僅能夠反映大規模數據的集中趨勢和離散程度,同時有助于從整體上把握數據集合的表現特征和內在規律[7].

設有一個N個樣本點和p個變量的數據表XN×p=xijN×p.根據數據系統中的某些重要屬性,將這N個樣本點劃分成n類子集合:S1,S2,…,Sn.如果用區間數據ek來概括Sk,常用的方法有:

ek=x-kj,kj=min i∈Sk(xij),max i∈Sk(xij),

j=1,2,…,p. (1)

在實際應用中,還可以采用適當的分位點作為每個區間的上下限.

一般地, n個數據集合Sk,k=1,…,n,可以被概括成一個n×p維的“樣本點×定量區間變量”類型的數據表,其形式為.

n×p=[x11,11][x12,x12]…[x1p,x1p][x21,x21][x22,x22]…[x2p,x2p]……

[xn1,xn1][xn2,xn2]…[xnp,xnp]

=(e1)′(e2)′(en)′.(2)

這里,ei=[xi1,xi1][xi2,xi2]…[xip,xip]′,被稱為區間數據表的樣本點.

由于區間數據的特殊性,在對以區間數據為單元的數據表進行主成分分析時,首先要按照一定的展平算法,將區間數據表展平為普通數據表.本文采用因素區間數據展平算法.該方法首先對每一個數據集合求主成分,然后分別用各個主成分的最大值和最小值來概括區間數據,并用高維超矩形的頂點作為樣本點生成普通數據表.該方法能夠以很高的精度代表原始數據的取值范圍,計算出的主軸誤差較小[8].

本文采用發文量、總被引頻次、篇均參考文獻量、影響因子、即年指標等5個文獻計量指標對CSSCI收錄的全部期刊進行區間數據的主成分分析.

采用Matlab軟件對區間數據進行主成分分析,得到結果見表2.

由表2的分析結果,前兩個主成分累計貢獻率接近70%,能夠較好地代表5個原始變量的大部分信息,因此模型提取第一和第二主成分對期刊進行綜合評價,其結果是合理和有效的.

表3的因子載荷矩陣反映了第一和第二主成分和5個原始變量的相關性.可以看到,第一主成分與5個原始變量均正相關,與影響因子、總被引頻次呈強正相關,說明第一主成分主要反映的是期刊被引用的情況,體現其學術影響力,以影響因子這一指標為代表.

第二主成分與發文量、總被引頻次正相關,特別是和發文量強正相關,與剩下其他三個變量負相關,說明第二主成分反映期刊的載文信息含量,以發文量這一指標為代表.

根據主成分分析的結果,可將符號對象在第一、第二主成分上的得分情況繪制成圖,觀察期刊的特征.如圖1所示,主平面圖的橫軸和縱軸的含義與主成分分析的第一主軸和第二主軸相同,描述期刊的區間數據在主平面上均以“十”字表示.其中,每一個“十”字的交叉點表示相應學科期刊的中心位置,而十字的長短則分別表示該類期刊在第一與第二主軸的離散范圍.通過觀察期刊在主平面圖上的位置,可以把握它們各自在兩方面指標上的表現特征,進而了解各學科期刊的整體特點.

圖1 因子載荷圖4 結論與建議

近年來,我國科技水平快速發展,學術論文的數量和質量也在不斷提高,然而相較而言,我國學術期刊與國際的差距卻非常明顯.因此必須正確認識我國學術期刊的發展現狀,采用客觀、科學、合理的手段評價學術期刊,指導學術期刊的管理工作和發展方向.

本文的研究對CSSCI數據庫中收錄的數百種學術期刊進行了全面、系統描述,并且通過建立模型,得出了一些有益的結論.

首先,學術期刊的發展差異很大,一方面反映在不同學科之間,另一方面反映在同一學科內水平參差不齊,因此在評價期刊時,不宜用同一套指標對全部期刊大排名,而用學科內排名.

其次,就評價指標而言,各種指標反映的信息主要集中在兩個維度.其一是體現被引用情況的一類指標,如影響因子、總被引頻次等,其次是反映信息含量的一類指標,如發文量等.在對學術期刊進行綜合評價時,要挖掘系統的內部特性,重點把握其主要差異所在,便于期刊的管理和遴選.

近年來,隨著經濟的發展和社會的變革,學術界的浮燥之風,也影響到學術期刊的質量.這一方面表現在有重大影響的社會科學成果和具有深遠意義的原創性成果比較匱乏,缺少創新,缺少健康和必要的學術探討,有水準的爭鳴和書評少之又少.另一方面則表現為低水平重復的社會科學成果大量剩余,大部分學術論文質量平平,甚至偏低[2].

學術期刊所面臨的問題已經引起了廣泛的重視.隨著科研體制的進一步完善,學術期刊也將與時俱進,順應時代和社會的發展,尋找自身的出路.

首先要提高期刊的綜合質量,引領學術創新.提高刊物的質量首先就是要爭取好的稿源,通過擴大交流,舉辦研討會等方式,主動爭取好的稿源.

其次提高編輯隊伍的專業素質.應由掌握除編輯學之外的某一領域的專業知識,且學有所長、有一定影響力的專家擔任學術期刊的領導,而其領導下的編輯工作者隊伍是由一批具備一定科研能力、在某一領域已有所成就或有較強科研潛力的中青年學者組成.學術期刊編輯自身的科研能力對于學術期刊發展具有十分重要的作用.

還有就是加快學術期刊國際化步伐.學術期刊承載著重要的社會責任,要努力推進合作中的學術交流,躋身于世界知名期刊行列.提高學術意識,反映我國學術研究的狀況,尤其是科技期刊要依托我國學術研究的專業優勢,真正在國際上代表這一領域的較高水平.學術期刊應與時俱進,打造亞洲一流乃至世界一流的學術期刊.參考文獻

[1] 杜玲.我國學術期刊數字化面臨的問題與對策[J].天水行政學院學報,2010,(6):114-117.

[2] 冷曉玲.學術期刊現狀和發展的思考[J]. 濰坊學院學報,2008,(5):119-122.

[3] 邱均平. 文獻計量學[M ]. 北京:科學技術文獻出版社,1988:2-34.

[4] 杜飛,房濤,葛學銘,等.1999~2008年中國應用生理學雜志載文統計分析[J].中國應用生理學雜志,2010,(2):252-255.

[5] 張建勇.中國科學計量指標:期刊引證報告[M].北京:中國科學院文獻情報中心出版社,2006.

[6] 譚果林.科技情報信息出版物有效性評價指標分析[J].太原科技,2009(10):42-44.

篇(3)

【關鍵詞】人臉識別 判別過程 歐氏距離 馬氏距離

隨著技術的不斷進步,人臉識別系統正在逐漸進入日常應用,在企業、住宅安全、刑偵、自助服務、信息安全等領域有著廣泛的市場前景。目前人臉識別方法主要有:幾何特征的人臉識別方法,基于主成分分析(PCA)的人臉識別方法等。

1 識別方法簡介

幾何特征的人臉識別方法研究開始于上世紀60年代末,通過考量眼、鼻、嘴等的形狀和它們之間的幾何關系實現人臉幾何特征的提取。此方法識別速度快,所占內存小,但識別率較低。

基于主成分分析的人臉識別方法主要是基于KL(Karhunen-Loeve)變換的識別方法,KL變換以矢量信號的協方差矩陣的歸一化正交特征矢量所構成的正交矩陣來對該矢量信號正交變換。在圖像處理中,高維的圖像空間經過KL變換后得到一組新的正交基,保留其中重要的正交基,由這些基可以張成低維線性空間。如果人臉在這些低維線性空間的投影具有可分性,將這些投影作為識別的特征矢量。PCA具有良好的去相關特性,KL變換后的矢量信號分量互不相關。同時KL變換是在均方誤差測度下,失真最小的一種變換。PCA也是目前人臉識別的主流方法。

PCA特征值提取時首先計算矩陣樣本的協方差矩陣,然后計算協方差矩陣的本征向量e1,e2,…,eN的本征值。通過PCA訓練多個樣本的特征值即可作為此人的特征,在識別過程中,用已有的特征值和提取出的本征值計算距離,作為人臉識別的判別依據。

2 主要算法分析

距離計算主流的方法有歐式距離和馬氏距離。

歐氏距離(Euclidean distance)可以描述為在n維空間內,最短的線的長度。歐氏距離可以表示為:

歐氏距離將臉部特征的之間的差別等同看待,并沒有對于主要特征重點考量,而馬氏距離可以解決此問題。

馬氏距離也可以定義為兩個服從同一分布并且其協方差矩陣為Σ的隨機變量與的差異程度。對于樣本集Xi,馬氏距離可以表示為:

其中為樣本平均值,S表示樣本的協方差。

馬氏距離不受量綱的影響,兩點之間的馬氏距離與原始數據的測量單位無關,馬氏距離還可以排除變量之間的相關性的干擾。

3 系統實現方法

在人臉識別系統的判別方法設計中,參考馬氏距離的加權思想,對特征值的歐氏距離進行加權。在提取人臉特征值時,取(樣本數-1)作為提取的特征值的數量。查看特征值所對應的圖片效果,以29個特征值為例,效果如下:

對于權值的選取,考慮到相近特征值所表示的特征之間差異性很小,并且其表示的實質性參量的不確定性,所以將一組固定數量且相鄰的特征值取同一權值。得到經過加權的距離后,計算置信度,置信度定義為:

其中 為訓練樣本值,aj為樣本權值。

如訓練樣本取30張人臉圖片,其特征值選30-1=29個,固定若干個特征值一組,每組用同一權值aj。特征值越靠前,表示其可參考性越大,所賦予的權值也應該越高。

對于非加權算法的歐氏距離,通常認為置信度高于0.5表示兩張圖片相似程度高,低于0.5表示相似程度低。但在人臉識別過程中發現,在訓練樣本不是特別多的情況下判決門限會隨著訓練樣本的增加而增加,而且判決門限值也會高于0.5,實驗中判決門限一般定在0.8以上,判決門限的設定也和前期圖片預處理的方法有關。

4 實驗結果分析

實驗訓練樣本數分別取30,40,50,60,70,綜合10個不同的人在同一位置分別計算加權算法置信度和非加權算法置信度來作對比。分析實驗結果,得到4條平均置信度隨訓練樣本數變化的曲線。訓練樣本數較小的時候置信度隨訓練樣本數的增加而增加比較明顯。但在訓練樣本數較大的情況下,置信度隨訓練樣本數變化不明顯。對于正確的樣本,加權算法所得到的置信度高于非加權算法得到的值,而對于錯誤的樣本,加權算法所得到的置信度低于非加權算法得到的值,所以采用加權算法可以加強判別過程的區分度,提高樣本的識別度。

參考文獻

[1]吳迪.基于特征臉改進算法的人臉識別技術的研究[D].哈爾濱:哈爾濱工業大學碩士論文,2007.

[2]李.基于主成分分析的人臉識別[D].濟南:山東大學碩士學位論文,2008.

[3]柯曉華.人臉識別算法研究及實現[D].哈爾濱:哈爾濱工程大學碩士學位論文,2008.

[4]孫偉,李曉飛.基于PCA的實時人臉識別系統[J].中國多媒體通信,2013(04).

[5]何振學,張貴倉,譙鈞,楊林英.對稱核主成分分析及其在人臉識別中的應用[J].計算機工程,2013(03).

[6]郭浩,王國宇.一種基于中值思想的改進人臉識別方法[J].現代電子技術,2013(02).

[7]李冠楠,李強.一種基于人臉核心特征的PCA人臉識別算法及應用[J].電子器件,2012(05).

篇(4)

關鍵詞:主成分分析 證券公司 競爭力

中圖分類號:F830.91 文獻標志碼:A文章編號:1673-291X(2011)28-0074-03

一、研究背景

隨著中國經濟持續穩定的發展,證券化程度的不斷提高,中國資本市場的規模不斷擴大,這為證券公司創新發展提供了良好的機遇,券商特別是已上市的證券公司的綜合競爭力也得到了長足的發展。在當前激烈的市場競爭中,如何對上市類證券公司的綜合競爭力進行衡量和分析,以達到科學比較,相互借鑒的良性市場機制,從而提升公司的競爭力,已成為證券監管部門、證券公司和投資者關注的問題。

關于證券公司綜合競爭力評價的研究,具有代表性的有單項指標評價法和綜合指標評價法。單項指標評價法以證券業協會公布的證券公司業績排名為代表,包括資產、利潤、承銷收入等單項指標排名,單項指標排名只能反映證券公司某一方面的能力,不能反映證券公司的綜合實力,具有一定的局限性。王曉芳、王學偉等(2008)選取中國50家證券公司為研究對象,選擇經營中的18個參考指標,運用因子分析法對證券公司2006年度的經營狀況做出了客觀全面的評價。運用綜合指標評價法對中國證券公司競爭力進行研究具有積極的意義,但是很少有學者運用綜合評價法對中國上市類證券公司股票競爭力進行評價和分析。

本文以15家上市證券公司為研究對象,選擇其2008―2010年度反映上市公司獲利能力和經營發展能力的六個指標,運用主成分分析法提取影響這些上市證券公司競爭力的綜合指標,并對這些公司的綜合競爭力進行了排名和分析,希望能夠給證券公司和投資者有所幫助。

二、主成分分析法的基本原理

主成分分析是在降低維度思想下產生的處理高維數據的統計方法,即用為數較少的、信息互不重疊的新變量來反映原變量提供的大部分信息,從而通過對為數較少的新變量的分析達到解決問題的目的。

主成分分析可廣泛的應用于經濟指標的評價,其具體的步驟如下:

首先,對原始指標數據矩陣進行標準化處理,并基于標準化的數據計算相關系數矩陣R。計算相關系數矩陣R的目的是通過其求出主成分。

其次,求出相關系數矩陣R的特征值λk(k=1,2,…n)以及各樣本主成分的方差貢獻率、累計方差貢獻率,并根據累計方差貢獻率選取主成分Zj(j=1,2,…r)。

最后,以每個主成分的方差貢獻率為權數,構造綜合評價函數F:

F=■αiZi,其中,αi=■,1≤i≤r。

并根據綜合評價函數對上市證券公司綜合競爭力進行排名。

三、實證結果與分析

1.樣本數據的選取及指標體系的構建。為了較為全面了解中國上市證券公司的綜合競爭力,特選取如下六個指標進行分析:X1:每股凈收益;X2:凈資產收益率;X3:主營業務收益率;X4:主營業務增長率;X5:凈資產增長率;X6:總資產增長率。其中前三個變量反映了上市證券公司的獲利能力,后三個變量反映了公司的經營發展能力。下頁表1給出了15家上市證券公司關于以上六個指標在2008―2010年三年取值的加權平均,權數分別取0.2、0.3和0.5。

由于廣發證券在2010年2月借殼巨虧的S延邊路上市,在2008―2010年度的相關財務指標異常,因此,本文選取其余14家上市證券公司為研究對象。

2.用統計數據進行主成分分析。首先,對原始指標數據矩陣進行標準化處理,并求出其相關系數矩陣:

RR=1 0.661-0.06-0.596 -0.152-0.361 0.661 1-0.591 -0.513 0.0700.013-0.06-0.5911 0.340 -0.309 -0.412-0.596-0.513 0.340 10.0290.343-0.1520.070 -0.309 0.029 1 0.742-0.3610.013-0.4120.343 0.742 1

其次,確定主成分分量。其特征值、方差貢獻率和方差累計貢獻率(如表2所示)。

表2各主成分特征值、方差貢獻率和累計方差貢獻率

由表2可知,前兩個主成分的累計貢獻率已達到76.033%(大于75%),這說明提取的兩個主成分能夠解釋這六個變量75%以上的信息。因此選取前兩個主成分做進一步分析,這兩個主成分分量的計算公式是:

Z1=-0.844X1-0.835X2+0.397X3+0.835X4+0.176X5+

0.381X6

Z2=-0.162X1+0.363X2-0.739X3-0.017X4+0.819X5+

0.864X6

分析上述數據和公式,發現:第一主成分Z1主要反映在X3,X4指標上,說明主營業務利潤高,主營業務增長快。第二主成分Z2主要反映在X5,X6指標上,說明凈資產和總資產增長快。在第一主成分和第二主成分中,X1指標的系數都是負數,表明每股凈收益指標在描述上市證券公司獲利能力上是不足的。上市證券公司獲利能力主要體現在反映主營業務指標的第一主成分上,經營發展能力主要反映在凈資產和總資產的增長速度上。

3.各主成分排序。基于上面提取的兩個主成分對14家上市證券公司排名,每個公司在各個主成分上的得分是按照主成分公式計算而得。然后根據提取的兩個主成分計算出綜合評價函數F,F=0.4067Z1+0.3536Z2。

該函數使用各主成分的貢獻率為權重。利用F函數對14家上市證券公司進行綜合排名。所有排名(如表3所示)。

四、結論與分析

從表3可以看出,在主成分Z1中,太平洋、西南證券和山西證券排名靠前,說明這三家公司獲利能力較強,在主成分Z2中,西南證券,山西證券和華泰證券排名靠前,說明這三家公司經營發展能力較強。在綜合評價函數F中,西南證券、山西證券和太平洋證券排名靠前,說明這三家上市證券公司的綜合競爭能力較強。在證券業協會的行業排名中居前的中信證券、東北證券在這里排名卻靠后,說明這幾家公司的綜合競爭力在證券類上市公司中不具優勢。無論是第一主成分、第二主成分還是綜合評價函數的得分比較接近,差距不大,說明中國的這些上市類證券公司處于平穩發展期,綜合競爭力差別不大。

上市類證券公司的排名情況與公司財務狀況有直接的關系,然而,選取不同的財務指標可能會產生不同的結果,并且主成分分析法對指標的選擇帶有一定的主觀性。同時,衡量公司競爭力的因素還有很多。

參考文獻:

[1]王曉芳,王學偉.基于因子分析的中國證券公司競爭力研究[J].現代商貿工業,2008,(1):139-141.

[2]哈達.基于主成分分析的上市公司核心競爭力評價研究[J].內蒙古科技與經濟,2009,(11):10-13.

[3]諸藝方,章亞男.制酒行業上市公司股票競爭力分析[J].經濟研究導刊,2010,(36):72-73.

篇(5)

關鍵詞:生態城市;主成分分析;綜合評價

“生態城市”是20世紀70年代聯合國教科文組織發起的“人與生物圈(MAB)”計劃首先提出了“生態城市”這一概念。隨著人口的迅速增長,工業化的深入發展,諸如城市的大氣污染、水污染、垃圾污染、地面沉降、噪音污染;城市的基礎設施落后、水資源短缺、能源緊張;城市的人口膨脹、交通擁擠、住宅短缺、土地緊張,以及城市的風景旅游資源被污染、名城特色被破壞等一系列問題日益突出。如何實現城市經濟社會發展與生態環境建設的協調統一,就成為國內外城市建設共同面臨的一個重大理論和實際問題。因此,建設生態城市已成為城市之間競爭的焦點,許多城市紛紛加快了城市轉型發展的步伐,把建設“生態城市”、“花園城市”、“山水城市”、“綠色城市”作為奮斗目標和發展模式。近年來,關于生態城市建設的研究成為國內外的熱點話題。曾芳芳介紹了國內外“生態城市”的研究背景以及中外學者對其內涵的認識,重點闡述了國內外學者對其的研究理論,并提出生態城市建設的理論研究在生態城市評估方面的不足,為后期的生態城市發展理論研究起著重要的作用;李文君等以西安市為例,對西安市生態城市建設進行分析,為以后具體的生態城市的建設研究提供了借鑒等。商洛作為國家“南水北調”中線工程的主要水源地之一,加強對商洛生態環境的保護,顯得尤為重要。本文以商洛市為研究對象,結合商洛市城市建設實際,構建商洛生態城市評價指標體系,運用主成分分析等方法,深入探討商洛市生態城市發展水平及其影響因素,為打造中國西部最美的山水園林城市規劃提供建設性的意見。

一、生態城市建設評價指標體系的構建與評價方法

(一)生態城市綜合評價指標的選取本文指標選取主要是通過頻度統計方法,參考關于生態城市評價的相關研究論文中的指標體系,結合商洛城市建設的實際,堅持科學合理、針對性強,反映生態城市發展本質的指標作為商洛市生態城市綜合評價指標體系進行科學評價,具體指標體系構建詳見表1。(二)數據來源及評價方法指標體系所用原始數據來源于《商洛市統計年鑒》(2007-2016),城市環境質量報告書、統計公報、政府部門的調研數據等資料和相關文獻資料、書籍、學術論文、研究報告等。本文擬采用主成分分析方法,通過對城市經濟、社會和環境3個子系統39個指標進行分析確定權重,綜合評價商洛市生態城市發展水平。

二、商洛市生態城市建設水平的實證分析

生態城市建設的主要目標就是實現可持續發展,是一項涉及到經濟、社會、人口、科技、資源與環境等子系統組成的時空尺度高度耦合的復雜動態開放巨系統的系統工程。因此,在生態城市建設過程當中,要以生態學理論為指導,充分考慮生態城市的自身特點,不僅關注生態城市的整體發展,還要關注城市社會、經濟、環境子系統在“關系”上的協調。文章從經濟、社會與環境三個子系統方面入手,通過SPSS19.0對這39個指標10年來的各具體數據進行主成分分析,提取影響生態城市發展的因子,更好地為商洛生態城市建設規劃提供參考。

(一)經濟子系統生態化建設水平的影響因素

對經濟子系統的12個指標通過主成分分析,確定兩個主成分F1、F2,由表2分析可以看出各指標與各主成分之間的相關程度。第一主成分與城鎮人均可支配收入、農村人均純收入、人均GDP、財政一般預算收入、固定資產投資占GDP比重等指標相關程度較高,均達75%以上,且都是衡量經濟發展水平的重要指標,主要反映了經濟發展的整體水平。F1中尤其是城鎮人均可支配收入、農村人均純收入、人均GDP、財政一般預算收入的相關程度高達97%以上,說明這幾個指標在反映商洛經濟發展水平上,對其影響更大。第二主成分主要與GDP增長率、一般財政預算收入增長率、旅游業增加值占GDP比重等指標有一定的相關性,這些指標相比第一主成分指標,相關程度所占百分比較低,但是對經濟發展的影響又不可忽視。經濟被視為一個城市、區域、乃至國家發展的“硬實力”,其發展水平的高低將直接決定社會發展的各個方面。因此,在建設生態城市的同時,要統籌好影響經濟發展的基礎指標和潛在影響力指標因子,才能為生態城市的建設提供更好的發展基礎和動力。

(二)社會子系統生態化建設水平的影響因素

通過對社會子系統15個指標進行主成分分析,確定三個主成分F1、F2、F3,由表3可以看出:第一主成分(F1)與人均保險費、在校大學生人數、城市化水平、衛生技術人員、醫院床位數等指標相關程度較高,其相關程度都在90%以上,這些指標既是對城市功能建設、社會保障、教育醫療等基礎設施方面的成績的反映,同時也是影響社會生態化發展的主要因素。因此,商洛市生態城市建設中社會子系統生態化的發展要對上述指標不斷優化和升級,有利于商洛生態城市的總體建設。與第二主成分(F2)相關程度較高的指標有:城鎮恩格爾系數、農村恩格爾系數、人均城市道路面積等指標,主要反映了商洛人民生活質量和城市社保障水平,這是對第一主成分的有力補充,其發展程度必然對社會生態化產生重要影響。第三主成分(F3)主要是常住人口、人均生活用水量、授權專利數指標,其中關于人口的發展又是對經濟、社會和生態發展等方面都起著巨大的影響。較多的人口可以為區域的發展提供充足的勞動力和消費市場;但是過多的人口也會對資源、環境保護等方面產生巨大的壓力。所以,可以將第三主成分解釋為人口質量與數量對商洛生態城市建設的影響因素。

(三)環境子系統生態化建設水平的影響因素

由表4可以看出,第一主成分(F1)與廢氣治理設施處置能力、城鎮人均公園綠地面積、工業廢水排放量、工業固體廢物綜合處置量、工業固體廢物排放量等指標其相關程度達到75%以上,這些指標主要說明工業三廢的處理水平及環境質量狀況。同時說明在生態城市的建設中對生態環境的改善應重點放在工業三廢的處理方面,不斷提高生態環境質量。第二主成分(F2)與污水處理率、工業固體廢物利用率、人均水資源量等指標有較大的相關性,說明商洛在生態城市建設過程中對水資源的保護及固體廢棄物的利用程度。第三主成分(F3)與全年城市空氣質量達標率等指標有關,主要反映了城市環境治理后的天氣狀況。通過分析說明,商洛生態化城市建設中對環境的保護不僅要注重工業三廢的處理,同時還要兼顧水環境污染和空氣污染保護防治,進而使商洛環境質量整體水平得到有效提高。

三、結論

生態城市建設的主要目標就是實現可持續發展,本文通過經濟、社會與環境系統的實證分析結果來看,生態城市建設涉及到各個方面,商洛市生態城市建設不僅要注重整體發展水平,還要重視經濟、社會、環境各子系統之間的協調發展,共同促進商洛市生態城市建設的步伐,促進商洛經濟又好又快發展。通過經濟子系統12個指標主成分分析得出:在生態城市建設過程中,經濟生態化要把重點放在與第一主成分相關程度較高的指標(城鎮人均可支配收入、農村人均純收入、人均GDP、財政一般預算收入、固定資產投資占GDP比重)方面,主要反映在經濟發展水平方面。但由于商洛地處秦嶺腹地,自然環境復雜,基礎薄弱,經濟發展水平低。因此在生態城市化建設中應不斷加強對外交流學習,借鑒成功經驗,大力發展特色區域經濟,不斷提高經濟發展水平,促進商洛經濟生態化建設步伐。但與第二主成分相關的指標也不容忽視。通過對社會子系統的15個指標主成分分析得出:在生態城市建設過程中,社會生態化的重點放在與第一主成分(F1)相關程度較高的指標(人均保險費、在校大學生人數、城市化水平、衛生技術人員、醫院床位數)等方面,不斷完善城市功能、社會保障、教育醫療等基礎設施的建設,促進商洛市社會生態化進程。通過對社會子系統的12個指標進行主成分分析得出:環境生態化重點要放在與第一主成分相關程度高的指標(廢氣治理設施處置能力、城鎮人均公園綠地面積、工業廢水排放量、工業固體廢物綜合處置量、工業固體廢物排放量)等方面。生態環境是經濟和社會協調發展的前提,隨著城市人口的不斷增多,工業化發展的不斷深入,商洛市環境生態化要不斷提高對工業三廢的處理及利用水平,合理開發利用各種資源,防治水污染和空氣污染,優化生態環境,不斷提高環境質量,確保經濟、社會和環境的協調發展,促進商洛市生態城市建設整體水平的提高。

作者:王高建 單位:商洛學院

參考文獻:

[1]曾芳芳.國內外生態城市建設的理論研究[J].科技廣場,2011(10)

[2]李文君;田巧娣.基于主成分分析法的生態城市建設綜合評價———以陜西省西安市為例.經濟視角(上旬刊),2015(1)

[3]周曉唯,王拓.陜西省生態城市建設評價指標體系[J].城市環境與城市生態,2009(6)

[4]趙維良,紀曉嵐,柳中權.主成分分析在城市生態安全評價中的應用——以上海為例[J].科技進步與對策,2009(3)

篇(6)

關鍵詞 PM2.5 線性插值 Box-Cox變換 主成分分析 預測分析

中圖分類號:X823 文獻標識碼:A DOI:10.16400/ki.kjdkz.2016.11.071

0 引言

從長春市2014年10、11、12月份月平均PM2.5濃度以及優良級天數監測結果(數據來自吉林省環境保護廳)可以看出:這三個月長春市優良級天數比例均小于50.0%,幾乎整個冬季都處于連續污染狀態,霧霾問題嚴重,再加上嚴寒的天氣,使得市民易于感染呼吸道疾病,甚至加重比如敏感、哮喘病等類患者的病情。因此準確預測并及時公布大氣環境狀況變得越來越重要。

近年來,專家學者開展了一些相關的研究工作。劉小生等①提出了一種基于基因表達式編程的PM2.5濃度預測研究;彭斯俊等②提出了一種基于ARIMA模型的PM2.5預測模型;陳俏等③提出了一種基于支持向量機和回歸法的大氣污染物濃度預測模型。這些文獻通過研究個別因素對大氣環境狀況提出了預測方法,但是氣象因素對PM2.5的影響是十分復雜的,實際情況中往往是不同氣象因素相互影響的結果。尤其是ARIMA模型只突出了時間因素在預測中的作用,沒有考慮到外界具體因素的影響,因而存在著預測誤差的缺陷,當遇到外界發生較大變化往往會有較大偏差。

本文旨在用與PM2.5濃度相關性強的因素,綜合考慮PM10、CO、NO2、SO2四項指標對PM2.5濃度的影響,并基于主成分分析理論提取幾個互不相關的主成分進行回歸分析,最終得到準確度較高的大氣環境預測模型。

1 材料與方法

1.1 數據預處理

1.1.1 補充缺失數據

本文數據來自天氣后報網,共研究PM2.5、PM10、CO、NO2、SO2五項指標,個別日期(共4天)的數據缺失。這時,我們使用SAS統計分析軟件,運用插值法補全缺失值。

1.1.2 Box-Cox非線性變換④⑤

將原始數據中PM2.5、PM10、CO、NO2、SO2等因子依次記為、、、、,因事先由散點圖分析可得,PM2.5與PM10、CO具有良好的線性關系,考慮到變換的簡便性最終選定對PM2.5、PM10和CO做變換%d的值為0,記變換后的PM2.5、PM10和CO為、和;對NO2和SO2進行Box-Cox變換的過程中,最優%d的取值是依據最大似然估計的方法原理來確定,由SAS統計軟件計算得到,最終選擇NO2和SO2的最優%d值依次為0.5、0,經過Box-Cox變換后的NO2和SO2依次用下列符號標記:和。

1.2 主成分分析原理

主成分分析⑥是將多指標化為少數幾個綜合指標的一種統計分析方法。這些主成分能夠反映原始變量的絕大部分信息,它們通常表示為原始變量的線性組合,且各個主成分之間互不相關。這樣在研究復雜問題時就可以只考慮少數幾個主成分且不止于損失太多信息,從而更容易抓住主要矛盾,解釋事物內部變量之間的規律性,同時使問題得以簡化,提高分析效率。

2 結果

2.1 主成分分析

本過程主要通過SAS軟件⑦實現,詳細程序參照附件。輸出結果(表1)給出了各變量之間的相關系數矩陣。可以看出:與之間的相關系數為0.8341,呈現非常強的相關性;

與,與之間的相關系數均為0.6800以上,有較強的相關性,其他變量之間相關性則相對較弱。不滿足多元線性回歸的基本假定條件,這也是本文選擇主成分回歸的主要原因之一。

輸出結果(表2)給出了相關系數矩陣的特征值、上下特征值之差、各主成分的方差貢獻率以及累計貢獻率??梢钥闯?,第一主成分的方差貢獻率為71.34%,前兩個主成分的累計貢獻率已達87.69%,因此,只需前面兩個主成分就可以概括這組數據。根據相關系數矩陣的各個特征值的特征向量,可以寫出前三主成分得分:

2.2 主成分回歸

現在用對前兩個主成分和做普通最小二乘回歸,得到主成分回歸方程為:

=3.93932+0.36567Ha0.01363 (3)

但是斜率的t檢驗p值0.4123>0.05,未通過顯著性檢驗,即認為與之間的線性回歸關系不顯著,需對模型進一步的調整。

經過多次對u和主成分、、之間進行不同模型模擬對比分析,最終采用逐步回歸法選取與、做最小二乘回歸,輸出結果如表3。雖然信息量從原來的87.69%降到了79.78%,但截距和斜率的t檢驗p值

為了得到和、、、之間的關系,運用R軟件⑧編寫計算系數的函數代碼作變換,得到還原后的主成分回歸方程為:

3 結果分析

3.1 殘差分析

回歸方程通過了t檢驗,只是表明變量之間的線性關系是顯著的,但不能保證數據擬合得很好,也不能排除由于意外原因而導致的數據不完全可靠,比如異常值出現、周期性因素干擾等。借助SAS統計分析軟件以回歸預測值作橫軸,以殘差 =Ha作縱軸,將相應的殘差點畫在直角坐標系上,得到殘差圖如圖1。

從殘差圖上看出,殘差是圍繞隨機波動的,表明模型中不存在異方差、相關性問題,模型的基本假定是滿足的。

3.2 預測

為了更好地反映模型的擬合效果,我們對2015年1月2日到3月27日的PM2.5進行預測,并與實際值對比。預測曲線圖如圖2,除了少數幾天PM2.5預測值偏差較大外,其它模擬效果比較樂觀,而且相對誤差相對較小,很大程度上接近了PM2.5的實際觀測值。進一步地證明該模型能夠較準確地預測未來短期的PM2.5趨勢和水平。

4 結論

本文利用長春市2014年空氣質量歷史數據,基于主成分分析理論建立了大氣環境預測模型,并預測2015年1月至3月的污染物PM2.5濃度,通過將其與實際值檢測值進行對比分析,表明運用此模型在PM10、CO、NO2、SO2各項指標已知的情況下,能夠對PM2.5進行準確性預測。

模型的判定系數和修正分別為0.9225和0.9221,均方殘差平方根為0.1866。這說明回歸方程的擬合效果非常好,能夠很好地反映并及時公_大氣環境狀況,讓公眾更加精確地感知到本地空氣質量的真實情況,還可以有的放矢地治理大氣污染,通過數據分析找出污染源頭。

注釋

① 劉小生,李勝,趙相博.基于基因表達式編程的PM2.5濃度預測研究[J].江西理工大學學報,2013.34(5):1-5.

② 彭斯俊,沈加超,朱雪.基于ARIMA模型的PM2.5預測[J].安全與環境工程,2014.21(6):125-128.

③ 陳俏,曹根牛,陳柳.支持向量機應用于大氣污染濃度預測[J].計算機技術與發展,2010.20(1):250-252.

④ 張誠.基于Box-Cox變換的城市火災起數的模型研究[D].合肥工業大學碩士學位論文,2013.

⑤ 胡宏昌,樊獻花.廣義Box-Cox變換[J].周口師范學院學報,2006.23(5):17-19.

⑥ 何曉群.應用回歸分析(第三版)[M].北京:中國人民大學出版社,2011.

篇(7)

[論文摘要]根據相似性度量理論中的χ2統計量,在專家評級的基礎上,構造出定性指標的相關系數矩陣,利用主成分分析法選取主成分,并作因子分析,進而根據因子載荷矩陣,得到主成分與各原始指標間的相關系數,最后在一定閾值標準下,舍掉相關系數絕對值較小的指標,從而達到用客觀合理的方法對定性指標進行篩選的目的。該方法可被借鑒應用于物流領域中供應商定性評估指標的篩選。

一、引言

供應商評估和選擇是企業的一個重要決策,一個好的供應商是指供應商擁有制造高質量產品的加工技術,擁有足夠的生產能力,以及能夠在獲得利潤的同時提供有競爭力的產品。同一產品在市場上的供應商數目越來越多,供應商的多樣性更使得供應商的評估和選擇工作變得復雜,需要一個規范的標準來操作。供應商評估首先要解決的是供應商評估指標體系的確立問題。

一個具體的綜合評價問題,評價指標的選取是否合適,直接影響到綜合評價的結論[1]。指標并不是選得越多越好,太多重復性的指標會增加指標權重確定的難度。綜合各學者的觀點,對供應商評估的指標可以包括技術水平、產品質量、生產能力、供應能力、價格、地理位置、可靠性(信譽)、售后服務、柔性等[2],這些都是定性指標。對定性指標的篩選,以前用得較多的是模糊數學的方法[3],主觀性較強,最后得出的結果可能差異非常大。這兩年,有一些學者引進粗糙集理論對指標進行屬性約簡[4],該方法實用性強,不失為定性指標篩選的一個好方法。然而,需要事先知道各樣本對應的決策屬性的值,使得該方法在很多的綜合評價問題得不到應用。

對定量指標的篩選有一些統計的方法,如:條件廣義方差極小[5]、極大不相關[6]、灰關聯分析[7]等,這些方法客觀性較強,評價結果比較穩定。定性指標與定量指標區別的關鍵在于定量指標是可度量的而定性指標不可度量,如果把定性指標定量化,定量指標的篩選方法就可以為定性指標所用。然而,對定性指標的定量化研究多采用模糊數學的方法[8],該方法應用結果并不能讓人滿意。那么,是否可以考慮避開定性指標定量化這一問題,另辟蹊徑,嫁接定量指標篩選的統計分析方法對定性指標進行篩選呢?

在定量的指標分析中,相關系數是一個重要的量,它能相當好地反映變量之間的線性相關的程度。主成分分析法就是建立在相關系數矩陣基礎上的,在主成分分析法基礎上作因子分析,根據因子載荷矩陣中主成分與各原始指標的相關系數的絕對值大小,就可以實現對原始指標的取舍。因此,對于定性指標,如果能求出兩兩指標間的相關系數,而不必對指標一一進行定量化處理,也能對指標進行篩選?,F在問題的關鍵在于,對定性指標是否能找出較好地反映變量間線性相關程度的量?Pearson的χ2統計量為這一問題的解決帶來了曙光。

二、用Pearson的χ2統計量求定性指標的相關系數

設x,y都是離散的隨機變量,x可以是r個狀態x,x,…,x之一,x不一定是數,但可以用指定的數或字母等來表示這一狀態,這樣的離散變量是定性資料的一種概率性的描述方法。類似地,y可以處于y,y,…,y這c個狀態之一。如果觀察了(x,y)的n個樣本,就可以計算出x=x,y=y的頻數、x=x的頻數以及y=y的頻數i Pearson據此定義了χ2統計量[9]

三、供應商定性指標篩選模型的構建及其求解

作為x與y關聯性的度量值。文獻[1]指出,χ2有下列性質:0≤x≤n(min(r,c)-1)(2)由此,如果引入θ=x/n(min(r,c)-1)θ≥0(3)θ就在[0,1]之內,它與相關系數ρ的平方在[0,1]之內相仿,那么θ就可以作為x與y相關系數的近似值。

相關系數問題的解決,為定性指標篩選方法的尋找提供了一塊墊腳石。主成分分析是建立在相關系數矩陣基礎上的一種多元統計方法,對該方法的進一步分析可以達到對指標篩選的目的。因此,結合Pearson的χ2統計量和因子分析法,可以構造定性指標篩選模型及其求解步驟如下:

(一)利用χ2統計量求定性指標的相關系數為避免符號過多引起歧義,對定性指標仍用x和y來表示,只是對x和y的描述稍加改變即可。為簡明起見,下面列了一個r×c雙向列聯表。指標x有r個水平x,x,…,x,指標y有c個水平y,y,…,y。觀察了(x,y)的n個樣本,其中指標x處在水平x、指標y處在水平y的個體數為表中之

分別是x處在水平x的個體數和y處在水平四、計算示例y的個體數。于是,首先計算式(1),然后計算式(3),把得到的數開方并取非負值,即可求得x與y的相關系數。同理,可以求出其他兩兩指標的相關系數,并構造出關于所有指標的相關系數矩陣。

(二)在相關系數矩陣基礎上用主成分分析法確定主成分求相關系數矩陣的特征值和特征向量,對特征值按照從大到小的順序進行排序,并一一計算貢獻率和累積貢獻率,按照累積貢獻率的值對綜合指標進行取舍。一般的,取臨界值為85%。一旦累積貢獻率超過85%,剩下的特征值對應的綜合指標即被舍掉了,前面的幾個綜合指標被確定為主成分。

(三)在主成分分析基礎上作因子分析,求得因子載荷矩陣對主成分作因子分析,即用各主元對應的特征值的平方根乘以該特征值對應的特征向量各分量,組合而成的矩陣即為因子載荷矩陣,該矩陣的各分量反映的是各原始指標與主成分的相關系數。

(四)確定閾值,對定性指標進行篩選在綜合評價中,主成分是最后確定為對樣本進行綜合評價的量。因此,原始指標與主成分相關系數越小,對綜合評價的結果影響力也就越小。對原始指標與主成分的相關系數進行排序,按一定標準確定閾值,把小于該閾值的相關系數對應的原始指標篩選掉,就可以達到對定性指標篩選的目的。

某生產企業針對供應商的評估列出下列指標:技術水平(X1)、產品質量(X2)、生產能力(X3)、供應能力(X4)、價格(X5)、地理位置(X6)、可靠性(信譽)(X7)、信息共享和相互信任(X8)、類似的產品經驗(X9)、地理范圍和提供的服務范圍(X10)、客戶的兼容性(X11)、售后服務(X12)、供貨期(X13)、交貨準確率(X14)、柔性(X15)[10]。現在要對某供應商進行評估,首先需要對評價指標進行篩選??梢钥紤]在供應商評估和選擇小組評級的基礎上用上述模型實現對上述15個指標進行篩選。

供應商評估和選擇小組選擇A、B、C、D、E五個等級之一分別對被評價對象的上述15個指標進行評價,其中,A、B、C、D、E分別代表“很好”、“較好”、“一般”、“差”和“很差”。9份評級結果如下表2:

現在要根據前述理論、模型及求解方計算出這15個指標兩兩間的相關系數。以標X1和X2為例,X1和X2有同樣的5個水平AB、C、D、E,分別對應腳標1、2、3、4、5比如

如此類推,窮極X1和X2的5個水平的率,即可計算出(1)式,然后計算(3式,最后開方取非負值,就可以得到X1X2的相關系數約為0.14。類似地,可以求其他指標間的相關系數,進而可以構建出下的相關系數矩陣:

由R的特征值我們計算出第一主成分的差貢獻率為79.6%,于是,我們只用第一成分這一綜合評價指標就能較好地反映原數據的變換情況。對第一主成分的特征值方再與第一主成分對應的特征向量各分量法指乘,這樣得到的向量的每個分量分別對應第一主成分與各原始指標的相關系數,相關系數絕對值的大小反映了綜合評價指標(第一主成分)與各原始指標的關系密切程度,因此,對相關系數先取絕對值,再按照絕對值的大小從大到小進行排列。排列結果如下表3:

由表3可以看出,X4、X9、X11這三項指標與綜合評價指標的密切度明顯低于其他指標。經過分析發現,供貨能力(X4)與供貨期(X13)和交貨準確率(X14)內容重合度很高,用供貨期和交貨準確率這兩項指標就可以代表供貨能力,因此供貨能力屬重復指標,應該刪掉。同時,類似的產品經驗(X9)和客戶的兼容性(X11)與綜合評價指標的密切度低很可能是因為對這兩項指標的界定沒有一個很明確的標準,從表2也可以看出9位評價者對這兩項指標評級差異度非常大。因此,X9和X11對整個評價中干擾度較大,也應該考慮刪除。

五、結束語

本文根據相似性度量理論中的χ2統計量,在專家評級的基礎上,構造出相關系數矩陣,進而用主成分分析法和因子分析法,對一個評價對象的多評價因素指標進行篩選。而且,為提高計算的精度,可以采取增加樣本的方法,使頻率盡可能地接近于概率。從前面的示例可以看出,基于Pearson的χ2統計量的定性指標篩選方法能取得比較科學、合理的結果,因此可以為涉及主觀因素眾多的供應商評估和選擇問題所借鑒。

參考文獻

[1]胡永宏,賀思輝,綜合評價方法[M],北京:科學出版社,2000:21

[2]魏際剛,鄭志軍,企業物流管理基礎[M],深圳:海天出版社,2004:42-43

相關文章
成年女人永久免费观看视频,国产成人精品亚洲线观看,国产91精品成人一区二区三区,激情久久男人天堂
日韩亚洲区欧美精品| 宅男午夜激情一区二区免费播放| 成人在线免费观看一区二区三区| 久久精品亚洲久久久| 久久精品av| 国内精品人妻久久毛片αpp| 3344成人在线免费视频| 宝宝好涨水快流出来免费视频| 亚洲国产精品久久免费观看| 超碰国产精品97| av网站免费在线观看| 国产黑色丝袜高跟鞋在线网站| 天堂а√在线官网| 国产精品黄片播放| 国产乱码精品一区二区三区四| 最新国产精品精品| 午夜久久福利影院| AV无码免费看| 亚洲av欧美aⅴ日本综合| 国产午夜免费啪视频观看视频| 日韩亚洲av电影在线观看| 久久精品一区二区三区国产av| 熟女人妻一区二区免费 | 亚洲乱码一区二区中文字幕| 免费色av天堂| 麻豆成人观看| 波多野结衣久久中文字幕| 再深一点太爽了舒服死了| av天堂伊人网| 一级毛片免费不卡观看| 亚洲专区12p| 自拍偷区色综合| 国产偷国产偷高清精品| 欧美久久激情图| 不卡久久av| 日韩电影一区二区在线看| 99久久精品午夜| 人妻久久中文字幕网| 国产做a爰片久久毛| 国产在线永久免费视频m| 人人妻人人澡人人爽的公开视频| 夫妇交换性4中文字幕| 精品久久久久久,| av亚洲aⅴ久久精品| 亚洲av免费毛片| 三级黄a毛片| 国产日韩一区二区三区精品不卡| 亚洲成人网av| 有码人妻中文字幕| 有码av在线一区二区三区合集| 女人18毛片免费看看| 高清一区二区三区四区乱码| 精品国产一区二区三区四区导航| 亚洲夫妻在线观看视频| 婷婷激情麻豆五月| 嘿咻男女动态图| 男人边吃奶边做细节描述| 午夜日韩色视频| 国内一级毛片免费观看| 又色又激情的免费视频| 噗嗤噗嗤好涨好爽太深了| 国产avxxx| 男人的天堂av最新版本| 亚洲综合最大av网站| 熟女国产专区| 无遮挡十八禁国产免费网站| 被十几个男人扒开腿猛戳| 最近国产中文字幕视频| 鲁死你AV资源站| 国产床戏无遮挡大尺度网站| 国产午夜视频在线免费观看| 日本a√视频在线| 最近2019免费中文字幕视频亚洲熟妇| 亚洲区自拍视频在线观看| 久久婷婷人人爽人人干人人爱| 老司机免费午夜福利视频| 国产床上激情在线观看视频| 99国产精品欲88av| 欧美第一页色| 欧美一级做a爰片久久毛片| 国产丝袜制服诱惑av| 日韩岛国大片av久久久久久久| 不卡av一区二区三区| 欧美亚洲精品91| 白袜男高中生GAY网站| 愿我如星卿如月,夜夜流光相皎洁 久久久久久久精品国产亚洲蜜臀 人妻中文制服中文 | 亚洲欧美日韩综合精品| 久久精品国产亚洲av香:| 日本蜜桃免费观看mv| 神马久久中文字幕| 亚洲国模私拍人体GOGO| 久久久久人妻高清| 上司部长出轨漂亮人妻| 国产动漫一区二区三区在线观看| 欧美aa在线观看| 99国产精品久久久久久另类| 西西人体大胆啪啪私拍| av免费在线观看全网站| 久久精品国产亚洲av香蕉明| 黄色片网站久久久久| 成人午夜特黄AAAAA片男男| 国产一区二区三区播放视频 | a级毛色黄片| 黑人妇女毛片| 日本少妇特殊按摩2| 国产av午夜福利电影| 国产精品嫩草影院一二三区入口| 亚洲中文字幕乱码视频| 色欲人妻综合网| 狠狠躁夜夜躁中文字幕| 国产欧美成人一区二区网站| 19禁免费视频无码网站| 91香蕉在线精品视频| 阿娇被陈冠希亲下面无遮挡| 97国语自产视频在线免费观看| 暖暖日本中文免费| 中文字幕av一区乱码| 老太交CHINESEBBW| 成熟女人牲交视频| 福利午夜久久| 羞羞影院午夜男女爽爽应用| 久久综合精品国产二区小说| 国产精品午夜性视频入口| 日本公妇被公侵犯中文字幕2| 久久人妻中文字幕av| 国产av色黄网站| 国产一级二级三级女| 疯狂边爱边吃奶视频| 永久免费看啪啪网站入口| 欧美性色欧美精品视频69| 99精品免费视频久久久久久久久| 真实宾馆乱子伦| 成年黄色视频在线播放| 精品国内亚洲2022精品成人 | 欧美老妇与ZOZOZO交| 国产粉嫩av高清久久| 国产午夜亚洲精品AⅤ | 人妻av无码专区久久| 免费看在线黄色a级毛片| 久久久久久久久久人人人人人人| 午夜免费一级片| 久久精品国产精品亚洲下载| 91丝袜美腿在线| 国产精品一线视频| 一本综合久久| 中国XXXX真实自拍| 在线观看免费亚洲黄片| 美女黄频视频网站大全| 宝贝把腿张开我要添你下边| 欧美一区二区三区免费的网站| 一级,二级,三级黄色视频| av在线播放亚洲精品| 亚洲精品成人va| 男女边摸边吃奶边做爽动态图| 午夜无遮挡免费视频| 欧美美性大战久久久久久| 亚洲三级中文字幕在线| 中国少妇高潮喷水免费可以看| 最新中文字幕亚洲区| 国产91精品小伙和老熟女泻火| 欧美的第一页| 久久精品夜夜夜夜夜久久蜜豆| 精品麻豆婷婷久久99| 日韩欧美一区二区三区国产| 朋友的丰满人妻| 最新2019国产天天| 一区二区三区国产免费电影| 色婷婷婷亚洲| 扒开未发育完全的小缝视频| 老司机午夜福利看片| 久久精品麻豆国产| 国产一级特黄高清在线大片| 91av看片在线| 免费高清欧美一区二区三区| 午夜网址观看| 久久av影院精品| 亚洲欧美精品综合一区二区三区| 综合激情丁香久久狠狠| 欧美视频一区二区三区在线| 成人男女在线播放| 粉嫩蜜臀久久精品久久久久酒店 | 做床爱免费视频在线观看| 欧美亚洲日本视频在线| 日韩国产欧美a| 婬荡的秘书呻吟波多野结衣| 亚洲三级中文字幕在线| 试看美女被日| 一区二区三区国产福利| 黄片av高清不卡免费看| 久久久一本精品| 十分钟在线观看免费观看完整| 黄网站色视频免费观看无下载| 成年男女免费视频网站无 | av天堂亚洲图| 亚洲高清免费不卡视频| 初学生免费自慰网站高清| 国内精品九九在线视频| 亚洲深深色噜噜狠狠爱网站| 久久激情影院久久| 男女啪啪无遮挡| 国产一区二区高清视频女| 新不卡在线播放av| 无码国产福利AV私拍| 国产一级a爱做片免费观看| 欧美日韩精品电影一区二区三区| 在线视频自拍久| 国产精品久久久久毛片真精品| 欧美人与动牲交片免费播放| 国产卡一卡二卡三卡10| 客厅乱H伦亲女| 亚洲av在线观看一区在线| 国产超黄a级视频免费看| 可以免费看黄片的视频| 九九久久精品国产亚洲av麻豆| 东京热人妻无码人AV| 欧美日韩一区二区高清在线播放| 大香蕉一区二区在线| 成人免费激情在线视频| 女同毛片免费播放| 亚洲欧美日韩精品人妻| 69一区二区三区在线观看| 国产精品乱码一区二区三区软件| 欧洲av福利大| 免费观看日韩欧美黄片| 绯色av一本一道道久久精品| 人妻中文字幕小视频| 国产成人福利免费视频打野战| 国产精品亚洲片av在线播放卡| 亚洲成a人ⅴ片| 欧美aaaaa级黑人又粗又长| 精品国产专区| 天堂av亚洲av国产av电影| 欧美日韩精品国产| 女人18毛片xxⅹ水真多|| 噗嗤噗嗤太深了啊快停下学长| 18日本学生无套高潮片| 天堂中文资源在线8| 精品一区二区三区四区在线| 伦精品一区二区三区| 久久久久久久亚洲中文字幕| 91色欧美视频| 老司机性感三级视频| 天天躁夜夜躁狠狠躁躁| 宅男噜66国产精品观看| 午夜福利久久久影院| 婷婷色综合www| 疯狂做受DVD播放免费| 亚洲精品色婷婷在线影院 | 久久精品国产亚洲av高另类| 国产精品尤物久久久| 免费看不卡的av| 熟女电影av网| 亚洲伊人精品网| 日本熟妇午夜| 美女露出100%奶头的视频| 欧美多毛熟女视频| 久久久久久精品一级毛片蜜| 亚洲中文字幕人妻系列| 麻豆国产麻豆| 欧美成人一区在线观看| 久久久久77777人人人人人| a级毛片在线| 国产精品视频大全一区二区| 少妇被黑人整得嗷嗷叫| 欧美日韩国产高清在线观看| 久久在热精品| 日本特别黄的视频| 极品粉嫩国产免费观看| 国产内射在线激情一区| 未发育学生的女A片在线观看| 嫩模被啪的呻吟不断| 欧美日本国产在线不卡视频| 男人操女人黄网站| 少妇熟女图亚洲| 色噜噜日韩精品欧美一区二区| 国产综合精品亚洲| 免费看成人国产一区二区三区 | 亚洲人午夜射精精品日韩| 久久人人爽人人爽从片a| 小h片在线播放免费毛片| 一个人在线看片网址| 欧洲毛片av片一级毛片| 色欲香天天天综合网站无码| 22中文网久久字幕| 国产精品一区二区三区免费观看| 国产福利精品99| 成人久久18秋免费网站麻豆| 福利国产在线视频| 99精品人妻少妇一| 精品人人爱人人做| 一区,二区,三区免费观看| 亚洲综合欧美日韩一区| 天堂av日韩av在线| 97国语自产视频在线免费观看| 欧美国产日韩一二三区| 少妇特殊按摩高潮不止| 黄色看点欧美成人黄色看点 | av不卡免费在线看| 嫩草影院福利版| 天天狠天天透天干天天| Al—娜塔莉·波特曼| 国产乱码精品一区二区16| 久久国产精品2020| 国产乱码在线观看一区二区| v电影v亚洲v欧美v国产| 八个少妇沟厕小便漂亮各种大屁股| ass拉美富婆大屁股pics| 亚洲av综合色区av| 男女人搞黄色视频在线观看| a级毛片在线网站| 亚洲精品色婷婷在线影院| 5D肉蒲团之性战奶水| 亚洲一区二区三区在线播放_| 精品国产乱码久久久久久站| 欧美人与拘性视交免费| 欧美成人看片a免费| 无码熟妇人妻AV在线电影| 14萝自慰专用网站| 97久久精品人妻人人玩| 国产suv精品一区二区在线| 亚洲大片av一区| 女人18毛片免费a级毛片| 国产成年人精品一区二区| 精品久久精品色综合| 亚洲中文字幕久久乱码| 久久精品亚洲日本桃乃木香奈| 疼死了大粗了放不进去视频| 深爱激情五月婷婷| 我要看免费的国产黄视频| 亚洲av片子| 男女爱爱好爽视频免费看| 久久天躁狠狠躁夜夜2o2o| 亚洲最大av电影网| 骚妇视频在线免费观看| 免费99久久| 亚洲中文字幕久久999| 午夜福利亚洲自偷自拍| 成年女人黄小视频| 午夜精品一区二区三区在免费看| 男同GAY片AV网站| 欧美精品啪啪一区二区三区| 91人人爽人人澡人人妻| 四虎国产精品永久入口| 亚洲国产精品久久久久久久简单 | 人人人妻人人澡人人爽欧美| 国产精品V日韩精品V欧美精品 | 免费观看精品久久| 18禁网站成人 | 禁无遮挡网站| 欧美成人色站在线视频| 国产超黄a级视频免费看| 一区二区三区视频日韩| 精品国产福利一区二区三区| 欧美精品第六区| 国产精品国产三级国产专不?| 久久特黄免费视频| 人妻系列中文字幕在线| 日本XXXXX片免费播放| 中文字幕无线码一区2015| 午夜福利在线观看免费高清| 人妻激情乱人伦| 日本欧美一区二区视频在线观看| 亚洲第一青青自偷自拍| 99久久精品国产精品久久| 亚洲午夜av电影在线观看| 黄色日韩一级二级| 不卡的av观看| 无国产精品内射视频| 精品日产一区二区三区在线观看 | 经典人妻一区二区| 成人av网站看| 亚洲国产色在线日韩| 苍井空被躁2小时在线观看| 我和岳交换夫妇爽完整版| 午夜福利在线观看免费高清| 七妺福利精品导航大全| 亚洲免费观看国产| 毛片a级毛片免费观看精品| 日韩高清动作片| 国产欧美激情一区二区2| 久久久国产一区二区| 国产精品国产精品一区| 免费A级毛片无码A∨免费| 国产精品毛片v一区二区| 在线观看另类av| 日日噜噜夜夜狠狠久久丁香婷婷 | 色播欧美在线观看| 成人亚洲欧美视频在线观看| 中文字幕人妻一区二区免费| 国产精品一区二区av片香蕉| 久久久久久久久久久大尺度免| 亚洲成人日韩在线观看| 青春草免费观看视频大全下载| 十八禁免费在线网站| 成年视频人免费站| 福利国产一区二区三区| 精品久久久久久,| 国产一区二区在线观看了| 中国看毛片一区二区三区| 亚洲av片毛片成人观看你懂的| 综合激情丁香久久狠狠| 九九视频精品免费在线观看| 久久久久久久99精品观看| 国产看片在线| 欧美激情第一区| 中文字幕av在线免费m| 国产精品99久久99久久久不卡 | 在线观看不卡无码A片| 国产色版视频免费| 极品S乳私人玩物白丝自慰| 波多野结衣av高清一区二区| 玉蒲团之玉女心经| 国产亚洲成av片在线尤物| 亚洲一区二区三区在线观看电影| 亚洲精品乱码一区二区在线| 最近免费高清视频观看日本| 18禁动漫肉肉无遮挡无码| 在线观看成人一区二区三区四区| 日日噜噜夜夜狠狠久久丁香婷婷| 久久人人爽人人爽爽久久| 久久精品熟女亚洲av麻豆精品| 99精品只有久久精品免费| 欧美国产日韩资源精品| 老少配老妇老熟女中文普通话| 黄片免费在线免费播放| 欧洲一区二区三区网站| 欧美成人精品一区二区99| 成年人的黄色视频永久免费| 国产成人精品午夜视频免费| 人妻av乱片av出轨九九| 青春草在线伊人| 又大又黄又免费的视频| 亚洲国产一区二区免费在线观看| 精品国产欧美一区二区三区成人| 国产精品91av一区| 最近中文视频字幕大全| 极品粉嫩国产免费观看| 色播在线免费播放| 麻豆一区av| 成人深夜av亚洲一区二区 | 亚洲国产av黄色片| 国产精品久久电影中文字幕| 亚洲伊人论坛| 欧美三级日韩国产在线观看| a级毛片毛片免费观看永久| 青草久久国产| 午夜男人的天堂网| 久久99精品久久久久mm| 97香蕉久久国产在线观看| 一二三四社区在线视频观看| 在线看免费视频999| 久久99蜜桃| 天堂中文字幕av| 国产人妻高清国产拍精品| 在办公室被弄到了高潮视频| 免费日本黄色视频的网站| 高清欧美性猛交xxxx少妇| 网禁国产YOU女网站| 国产国产三级精品久久日本| 国内揄拍99国内精品人妻久久| 999精品视频在线免费看| 精品国产亚洲av日韩| 日韩a级福利视频| 91久久精品国产一区二区三区| 日韩在线看大片| 午夜精品成人在线观看视频| 国内精品人妻久久| 少妇人妻系列无码专区| 国产伦在线观看视频免费| 国产精品老妇久久精品老妇| 欧美狠狠久久免费看| 公交车被CAO得合不拢腿| 日韩专区亚洲| 女人18毛片免费看看| 午夜久久久久久久99精品婷婷| XXXXX无尽的毛茸茸撒尿| 懂色av天堂| 久久久精品欧美一区| 另类亚洲欧美激情| 久久综合国产亚洲精品| tobu8在线观看免费直播| 成人国产亚洲精品av天堂涩爱| 青草久久97| 亚洲国内久久精品| 国产成人av片在线观看| 亚洲精品美女一区二区| 国产精品国产三级国产av无| 欧美人与性动交q欧美精品| 久久免费精品视频99| 大陆一级a毛片| 91精品蜜桃| 亚洲成人av在线播放网站| 日韩国产精品久久久久久久蝴蝶| 天天躁夜夜躁狠狠久久av| 美女隐私羞羞视频在线观看| 国产人伦真实精品视频| 国产精品网站免费观看| 精品香蕉官网在线观看 | 91香蕉在线精品视频| 2亚洲精品一区二区三区| 欧美色逼逼视频| 欧美日韩一中文字幕| 欧美国产日韩一二三区| 国产高清黄色的视频在线播放| 久久人人妻人人爽人人澡| 91精品国产综合久久香蕉网| 极品美女写真| 天天狠天天透天干天天| 中文字幕aa在线| 中文字幕老妇偷乱视视频在线| av大片在线免费看| 所有免费永久岛国av网址| 亚洲精品自拍偷| 亚洲综合一区自偷自拍| 成人久久av一区二区三区| 日韩有码av电影在线免费观看 | 久久九九99综合一区二区| 国产精品原创巨作AV女教师| 国产福利不卡在线观看| 亚州一区二区三区av无玛| 久久婷婷国产综合精品| 亚洲免费观看av网址| 国产91色在线|日韩| 天堂网2021中文字幕公司| 国内少妇人妻精品视频三区四区| 一卡二卡国产免费| 亚洲av日韩av在线| 小美女毛多水多国产毛片| 日本在线视频www色| 久久这里只有精品777| 国产一区二区欧美日韩精品| 男女生羞羞视频网页| 国产一区二区三区午夜精品| 99久久99这里只有免费费精品视频| 欧美乱妇无乱码| 亚洲精品国产av成人精品| 精品国产日韩亚洲一区| 久久精品超碰国产酒店| 国产精品亚洲片在线观看不卡| 欧美性猛交aaaa乱大交| 亚洲专区av一区| 亚洲AV日韩AV永久无码| 亚洲AV日韩AV永久无码| 影音先锋啪啪AV资源网站APP| 国产偷窥熟女高潮精品视频| 学生和老师XXXXWWW| 影音先锋新男人AV资源站| a中文天堂在线官网| 狂干人妻少妇视频| 中文字幕日韩一区二区不卡| 99久久精品午夜| 国产一区二区人妻三区| 国产男女内射视频| 熟女少妇一区二区三区四区| 亚洲特黄视频免费观看| 亚洲精品午夜级久久久久| 欧美一二三四区精品中文字幕| 国产一区午夜福利在线观看| 亚洲免费观看国产| 亚洲av免费电影在线| 永久免费福利网站| 国产精品国产三级国产普通话三| 妺妺窝人体色WWW看美女| YOUJIZZCOM中国熟妇| 亚洲AV无码专区国产乱码| 一卡二卡三卡免费看| 一卡二卡三卡免费看| 美女露出尿口喷水无遮挡| 亚洲午夜精品av片| 日本熟妇午夜| 国产做a爰片久久毛| 久久综合色爱综合欧美| 久久国产99精品日本三级| 亚洲成人av在线播放网站| 天堂在线资源中文6| 综合色站国产| 亚洲熟妇av一区二区三区漫| 无码做爰欢H肉动漫网站在线看| 亚洲四区av| 亚洲码和欧洲码的区别三| 黑人巨大精品欧美欢迎你| 国产精品毛片av一区二区三区卜 | 欧美日韩一区二区在线观看国产 | 538在线观看一区二区三区| 怡红院亚洲免费| 欧美成人精品一区二区99| 日韩AV一区二区三区无码 | 色老99久久九九爱精品| 在线精品中文亚洲中文字幕| 日韩美女色视频网站全部免费| 日本人又色又爽的视频| 亚洲国产精品久久| y111111国产精品久久久| 久久久久久久婷婷免费视频| 午夜福利av线上网站| 99久久婷婷国产麻豆精| 人人澡人人妻人| 熟女人妻视频一区| 成年人黄色免费视频| 欧美性虐一区二区| 新搬来的四个爆乳邻居| 91av看片在线| 成人www网站视频在线看| 女人十八水多毛片免费看| 麻豆成人久久精品二区三区91| AV电影网站| 青草久久国产| 亚洲天堂高清在线观看| 国产亚洲精品久久一区二区| 午夜精品久久av蜜桃| 国产精品国产高清国产av| 久久久精品亚洲不卡av| 国产精品一区二区av麻豆蜜臀| 国产激情av女片自拍| 亚洲欧美日韩aⅴ一区二区三区| 老司机午夜十八禁免费视频| 人妻一区二区三区有码精品视频| 床震边吃奶边做视频| 国产内射999视频在线| 国产片特级美女逼逼视频| 欧美亚洲精品免费| 在线观看成人一区二区三区 | 蜜臀av伊在人亚洲香蕉精品区| 日韩精品中文字幕在线免费观看| 俄罗斯O|老太和小男| 久久精品亚洲精品国产色婷小说| 国产亚洲精品福利在线观看| 一级a爱做片观看免费网站| 又爽又黄a免费视频| 欧美黑人深夜精品免费| 国产高清网站在线观看91| 波多野在线AV一区二区荷花影视| 久久国产av人人做人人爽| 97人人模人人爽人人喊视频| 国产又粗又头又长又爽| 又粗又长又大又黄少妇毛片 | 亚洲av永久天码精品天堂dl| 999精品视频在线免费看| 欧美成人精品午夜视频播放| 午夜电影国产精品| 乱码午夜av噜噜噜噜| 中文字幕人妻有码av一区| 中文字幕av一道本| av天堂伊人网| 91亚洲精华国产精华精华乳| 麻豆国产av男优征选麻豆村| 国产日韩欧美精品另类一区二区| 中国浓毛少妇毛茸茸| 扒开女人内裤猛烈进入的视频| 婷婷在线中文字幕| 国产美女在线高潮精品| a级毛片免费观看大全| 苍井空高潮喷水在线观看| 熟女国产专区| 精品久久成人网| 人人妻人人爽av| 亚洲欧洲另类图片区| 好了AV四色综合无码久久| 亚洲国产综合欧美| 国产在线观看三级视频| 亚洲av成人精品一区久久| 特黄三级毛片在线观看| 善良迷人的女教师2中文| 熟女人妻波多野结衣av作品| 亚洲精选av在线观看| 亚洲国产404h网站| 2021在线精品自偷自拍无码| 欧美日韩成人伦理电影| A毛片免费全部播放无码| 美女被高潮视频网站| 欧美亚洲狠狠狠| 美女求操视频在线观看| 午夜黄色视频在线观看免费不卡| 黄网站视频大全影视频大| 九九久久精品影院| 免费乱理伦片在线观看2018| 96久久精品人人妻人人爽| 人妻在线中文字幕av| 三级黄色a大片| 国产毛片一卡二卡三卡四卡| 久久人成人国内精品| 欧美一区二区激情电影| 亚洲久久av在线| 国产高清黄色的视频在线播放 | 人妻有码一区二区三区| 日韩中文字幕欧美一区二区| 亚洲加勒比少妇无码AV| 国产成人精品一区二区三区福利| 放荡的女教师3在线观看| 成人国产精品福利| 日韩精品国产高清在线观看| 国产人妻熟女av久久| 午夜a级毛片| 人人妻人人做人人爽夜爽视频| 最大胆裸体人体牲交免费| 日韩亚洲av电影在线观看| 美女户外直播野战在线播放| 青春草在线观看国产| 亚洲成年人网站在线播放| 极品午夜福利1000在线| 人妻中文制服中文| 精品免费视频99| 欧美久久精品一级c片| 6一12泑女WWW雏| 亚洲国产精品亚洲一区二区三区 | 成人久久av爱乃娜美免费观看| 欧美日本一本线播放三区| 中文字幕亚洲无线码2020| 国产一区二区三区不卡免费视频| 美女露100%奶头18禁| 九九黄色视频免费看| 亚洲精品一区中文字幕在线| 女人在线免费视频| 国产成人一区二区三区西西视频 | 亚洲婷婷成人网| 久久国产xxx| 精品国产伦一区二区三区在| 久久久久人妻高清| 97人人模人人喊人人爽| 国产一级高清av| 免费在线观看永久免费av| 久久精品一区二区三区福利网| 美女久久亚洲150p| 精品一区久久三区二区| 亚洲加勒比少妇无码AV | 日本JAPANESE丰满同事| 日本高清色高清| 自线自在国产av| 色婷婷亚洲欧洲图片| 密臀久久精品久久久久酒店| 欧美日韩精品成人二区a∨| 人人妻人人澡人人偷拍| 永久中文免费在线观看| 日韩av一区二区三区| 热久久99热国产精品首页| 中文字幕精品一区二区三区我和| 久久av蜜桃一区二区三区| 国产视频福利免费看| 少妇毛片视频| 人妻毛片视频网站| 最近中文在线国语| 午夜免费福利不卡顿视频| a天堂tv香蕉| 5B肉蒲团之性战奶水国语| 免费一级做a爰片性色毛片| 免费1级a做爰片观看69网| 男女边摸边吃奶边做视频韩国| 天堂а√在线中文在线官网| 自拍视频一区二区三区| 国产黄色三级三级三级av| 最新91中文字幕在线观看| 久久精品麻豆国产| 男女视频福利网| 30位美人妻无套中出影片1 | 国产偷国产偷亚洲高清日韩| 日韩亚洲国产另类| 日韩成年人av| 午夜精品久久久久久久99桃| 天天干夜夜超| 男男啪啪激烈高潮漫画免费| 看a级黄片毛片| 男女视频免费观看| 久久国产乱子免费精品| 亚洲国产综合无码一区二区| 免费精品乱码一区二三区| av国产熟女丰满熟女| 国产avxxx| 久久九九有精品国产| 免费99精品视频| 综合久久亚洲精品| 把腿扒开做爽爽的事| 88av影院首页| 91人妻人人妻| 岛国视频在线观看免费一区二区| 小雪的13又嫩又紧又多水| 粉嫩极品国产在在线播| 日本在线视频2区| 91av永久免费视频| 亚洲av成人精品一二三区| 性老太大OLD| 人人妻人人狠人人爽欧美| 国产精品高清亚洲不卡| 亚洲在线观看国产精品| 在线视频 激情图片 三级| 国产精华免费| 国产精品久久电影中文字幕| 日本黄色视频网址在线| 午夜黄大色黄大片美女| 97热精品视频在线观看| 美国精品高潮呻吟久久av| 凸凹人妻人人澡人人添| 精品亚洲网站| 性在线观看视频| 国产丝袜美腿诱惑久久网| 国产熟女极品| 国产精品女人视频一区二区三区| 亚洲自偷自拍三级| 色婷婷国产精品久久包臀| 粗壮挺进邻居人妻| 幻女BBWXXXX国语| 成人毛片15女人毛片免费看| 骚妇视频在线免费观看| 亚洲中文字幕在线天堂| 最近2019在线中文字幕4| 国产精品毛片一区二区16| 久久久人人人人| 日日搞日日撸夜夜搞夜夜干夜夜插| 欧美最新免费一区二区三区| 麻豆果冻国产91在线极品| 亚洲加勒比少妇无码AV| 18禁美女裸身图无遮挡| 免费在线黄色av网站| 丝袜美腿免费av| 免费乱码欧美综合欧美日| 免费一区二区三区在在线视频| 亚洲黄色av日韩在线观看| av毛片免费看| 久久九九99综合一区二区| 一本精品久久a| 亚洲欧美日韩无卡精品| 搡女人真爽免费视频火全软件| 野花社区观看免费观看视频 | 手机a级毛片免费观看| 男人天堂国产精品| 亚洲成人a级电影| 欧美精品一区二区三区aⅴ在线| 一本色道久久综合婷婷| 视频一区图片区小说区| 制服诱惑二区| 色播欧美在线观看| 亚洲精品夜夜澡人人爽| 亚洲精品aa在线观看| 国产精品av久麻豆| ccc36免费在线视频| 国产日本欧美视频一区懂色| 在线视频网站中文字幕| 午夜精品成人在线观看视频| 亚洲高清视频在线观看免费| 亚洲另类国产日韩| 亚洲欧美一区二区三区| 中国看毛片一区二区三区| 黄色91影院| 插bb射bb福利免费午夜视频| 奇米影视7777久久精品人人 | 国产视频午夜福利| 疼死了大粗了放不进去视频| 午夜精品久久av蜜桃| 乱淫中文字幕视频| 亚洲精品,欧美精品| 日本妇女买淫| 中文字幕一区二区三区在线播放 | 亚洲欧美精品极品电影| 日日夜夜操撸| 免费看啪啪啪av| 亲子乱子伦XXXX| 从女儿九岁开始做| 亚洲永久免费视频最新网站| 午夜小视频一区二区三区| 自拍视频一区二区三区| 99热这里只有的精品99| 蜜桃在线视频一区| 亚洲欧美日韩久久精品狠狠| 台湾毛片在线播放| 午夜看片久久久久久| 亚洲成人av麻豆av| 国产日韩+欧美| 精品视频永久免费观看| 精品久久最新| 美女扒开双腿露出尿口无内裤图片| 精品va天堂亚洲国产麻豆高清| 国产在视频线在精品| 久久久国产午夜成年影院| 久久久久国产精品美女毛片| 黄频视频免费看| 一级做a爰片久久毛片50岁| 欧美日韩亚洲天堂久久| 亚洲成AV人的天堂在线观看| 四虎久久精品国产亚洲av| 大桥未久的av作品| 久久精品人人妻人人爽| 日韩人妻av一区二区三区| 日韩亚洲欧美三区中文字幕| 羞羞午夜福利在线| 国产又爽又黄的视频| 日本高清视频色视频免费| 国产成人精品综合久久久看| 产后漂亮奶水人妻无码| 日韩av情色在线| 女生毛片免费观看| 91中文字幕在线永久在线| 小h片在线播放免费毛片| 暖暖视频播放免费韩国| 国产精品酒店情侣视频| 欧洲极品诱惑另类视频| 亚洲av电影一区二三区| 国产高跟丝袜麻酥酥在线| 澳门一级a毛片| 91老司机视频在线| 日韩亚洲欧美精品综合| 天天躁夜夜躁狠狠综合2020| 欧美日韩一区二区三区色| 69一区二区三区在线观看| 中文字幕人妻精品乱又码| 不卡的熟女av| a√天堂中文资源在线| 男人的天堂久久激情| 国产超碰人人爽人人做| 91av永久免费视频| 亚洲欧美精品极品电影| 国产成人精品免费视频大全| 国内精品久久久久久久电影| 久久精品国产亚洲成人av| 老鸭窝在线视频免费观看网站| 开心五月激情综合婷婷| 扒开腿舔的视频| 三级精品国产| 波多野结衣久久中文字幕| 日本做受高潮好舒服视频| 精品亚洲一区| av在线不卡观看一本| 一本大道中文日本香蕉| 97人妻人人澡人人爽| 国产伦久视频免费观看+视频 | 国产精品午夜福利免费视频| 亚洲男同志GAY片可播放| 亚洲熟妇熟女久久| 亚洲欧美综合一区在线观看| 国内自拍偷国视频系列| 欧美亚洲久久综合精品| 久久视频伦理片| 轻轻挺进新婚少妇身体里| 日本一级a爰片| 精品亚洲999| 国产精品九九久久_久久国产| 97人人模人人爽人人喊免费从| 5566中文字幕在线观看| 亚洲中文字幕在线无码一区二区 | 久久久久久好爽爽久久| 黄色片一级av| 看黄色片子一级片| 午夜福利亚洲自偷自拍| 强d乱码中文字幕熟女| 亚洲日本不卡| 又黄又爽的视频免费99| 特黄特大黄色片| 国产精品免费福利久久| 亚洲一区二区三区最新网址 | 99午夜福利一区二区| 亚洲精品美女久久av网站| 最新欧美日韩国产在线观看| 99久久精品视香蕉蕉| 天堂av一区2区| 久久婷婷综合97色一本一本| 免费看小黄鸭av片成人| 国产精品久久男人的天堂| 中文字幕91专区欧美| 精品一区二区不卡视频| 漂亮人妻被同事疯狂玩弄| 一级做a爰片久久毛片免费看| 男女无遮挡羞羞视频| 亚洲欧美日韩在线播放一区二区| 女人十八水真多毛片| 毛茸茸厕所偷窥XXXX| 国产免费观看a大片的网站| 超碰人人人妻97| 成人久久亚洲欧美| 一区二区三区在线观看视频了| 一区二区日韩电影院| 国产亚洲午夜精品a| 国产在线一区二区不卡| 亚洲大色av| 国产av麻豆久久久久久久| 亚洲大色av| 国产97在线色/日韩| 久久久久久久久久成人精品| 国产97精品一区二区| 欧美九一精品久久久熟妇人妻| 最好看最新中文字幕2019| 国产精品一区二区hb在线观看| 漂亮人妻被夫部长强了| 日韩久久久黄色一区二区| 看免费黄色一级高清| 激情内射视频网站| 妓女影院永久免费| 国产午夜精品久久小说| 18禁动态无遮挡网站| 国产精品午夜视频一区| 高清视频免费观看一区二区| 春色成人av在线播放| 久久精品影院6| 一级做a爰片久久毛| 亚洲免费观看成人av| 中文字幕亚洲无线码欧美| 国产高潮美女av| 国产一区二区的免费在线播放| 久久久午夜人妻精品| 粉嫩精品国产尤物在线观看| 欧美狠狠久久免费看| 男人边吃奶边做好爽免费视频| 日本人人妻人人澡人人爽| 午夜黄大色黄大片美女| 最近免费中文字幕大全在线看| 舔舔嫩穴后中出| 国产三级av大全| 久久这里只有精品777| 亚洲av午夜精品久久看一区| av精品一区二区| 有码人妻中文字幕| 欧美日韩精品人人妻人人爽| 白饿罗斯毛片| 国产高清精品福利私拍国产写真| 手机亚洲av| 老司机91看片| 涨精装满肚子用塞子堵住| 麻豆成人观看| 精品欧美国产一区二区三| a级毛片1一级毛片免费看| 黑人巨大精品欧美黑寡妇| JAPANESE杂交人禽交| 91嫩草在线播放| 亚洲国产精品成人综合色| 欧美伦理在线第一页| 岛国av一区二区三区久久精品| av色黄色在线| 亚洲一区二区精品国产| 国产成人精品久久久久白嫩| a级毛片特级毛片免| 秋霞网站黄色一级国产毛片| 老司机精品视频免费入口| 精品国产专区| 免费吃奶摸下激烈视频青青网| a级毛片免费在线看| 久久中文字幕有码| 久久精品亚洲精品无码| 久久亚洲精品色一区| 99国产精品欲88av| 无码无遮拦午夜福利院| 亚洲精品午夜久久av蜜臀| 亚洲a级毛片免费播放| 一区二区三区久久电影院| 国产成人福利精品在线| 成人毛片18女人毛毛片免费看| 黄视频在线国产| 久久99就热视频精品| 亚洲天堂网资源在线| 老熟女xxxx| 欧美性群另类交| 看a级黄片毛片| 中文字幕熟女丝袜| 大肉大捧一进一出好爽免费视频| 男女爱爱好爽视频免费看| 天堂а√在线中文最新版官网| 亚洲国产综合自拍av| 久久久久久久久久精品尤物| 内射国产在线视频| 久久亚洲精品男人天堂| 亚洲毛深熟女射精精品hd| 久久精品免费视频17| 女人被躁到高潮嗷嗷叫免费| 美女被扒开内裤桶屁股眼| 国产精品欧美一区色| 欧美日韩国产另类久久| 欧美日韩综合久久久久久| 性夜影院爽黄A爽免费看| 美女高潮在线免费观看| 香蕉视频这里只有精品| 宅男噜噜99国产精品观看免费| av网址在线播放网站| 国语乱码一二| 午夜电影一区二区三区| 国产在线偷看观看免费| 成人av在线播放一区| 久久久久久88,久久se精品| 成年美女黄网站色奶头游戏| 国产美女精品免费视频 | 亚洲国产成人老色批| 国产精品久久久久精品| 搞黄视频大全在线| 亚洲欧洲日本精品一区二区三区| 久久午夜精品福利一区二区| 成人影院亚洲| 欧美www.一区二区三区| 69av精品久久久久久| 日本不卡一卡新区手机| 国产 国语对白 露脸| 国产一级精品午夜视频| 天堂av亚州av| 黄片视频无遮挡观看| 偷拍视频一区二区三区四区| 成人毛片a级毛片免费观看99| 久久99精品国语久久久| 日韩久久九九精品一本| 被领导添下面好爽| 美女极品粉嫩美鮑20P图| 国产福利精品免费在线观看| 人妻3p精品视频| 2021国产精品一区二区22| 国产成人亚洲综合色影视大| 人人妻人人澡人人爽人人| 熟妇人妻精品一区二区初频频| 国产毛片一区二区一区三区女| 亚洲av熟女| 国产成人亚洲无吗淙合青草| 日韩欧美亚洲久久| 东京热人妻无码人AV| 亚洲精品狼友访问紧急通知| 欧洲亚洲免费av| 免费乱码欧美综合欧美日| 亚洲欧美在线看看| 欧美交性又色又爽又黄麻豆| 午夜影院在线观看国产精品| 午夜黄色视频在线观看免费不卡| 国产精品99久久免费| 精品亚洲成人av| 青草永久福利| 啪啪啪亚洲综合色美利坚合众国| 97人人模人人爽人人喊谢绝| 精品黄色美女在线视频| 国产欧美亚洲在线| 国产午夜av毛片| 999欧美精品一区二区三区| 欧美精品一区二区三区四区五区高 | 日本免费一区二区三区视频日本| 精品99-区二区三区免费九一制片| 久久中文字幕有码| 国产91精品一区二区麻豆亚洲| 成人影院久久| 精品伊人久久大香线蕉| 永久久久国产精品电影| 女人18毛片a级18毛片水真mo| 久久久久久久久久精品尤物| 人人爽人人网| 人妻销魂中文字幕av| 极品少妇高潮喷水抽搐91| 国产一区二区三区在线观看免费| 日韩在线免费观看视频网站| 麻豆精品成人| 蜜桃一区二区三区精品av在线| 亚洲区日韩欧美国产| 国产精品zjzjzj在线观看| 国产熟女高清| 亚洲欧美熟女一区二区| 亚洲成人精品少妇| 深夜爽爽福利GIF动态图在线| 好男人看视频免费2019| 人妻电影久久| videos熟女内射| 久久久久久久久精品精品| 亚洲最大影院久久久久久| 一级片久久片| 欧美黑人a区| 18禁高潮娇喘喷水视频| 91麻豆天美精东密桃传媒| 亚洲欧美一区二区成人精| 免费成人激情在线视频| 最近中文在线国语| 嫩草影院av网址| 亚洲成人手机| 亚洲黄色av视频| 亚洲国产一二区在线观看| 国产麻豆成人专区| 欧美伦理在线第一页| 最近中文在线国语| 岛国视频午夜一区免费看| 日本禁漫免费观看| 日韩亚洲专区在线| 午夜免费爱爱视频| 美女被扒开内裤桶屁股眼| 欧美精品成人a多人在线观看| 国产精品久久人妻色戒| 国产一区二区三区不卡免费视频| 青草永久福利| 男女羞羞无遮挡免费| 午夜精品av国产| 北条麻妃人妻AV在线专区| 久久久久人妻精品一区果冻| 免费av大香蕉| 国产suv精品一区二区在线| 日本欧美v大码在线| 里番库全彩本子彩色本子H| 国产乱人偷精品视频| 久久国产av精品| 高清免费安全美女视频| 国产色综合视频在线| JAPONENSIS18一40| 精品欧美国产一区二区三| 日韩精品毛片久久久久久久| 有码一区二区三区四区| 浴室挺岳双腿之间| 亚洲视频在线观看视频在线观看| 欧美va亚洲va日韩∨a综合色| 欧美视频在线观看日本| 亚洲美女又黄又免费视频| 亚洲午夜精品国产电影院av| 久久人人爽人人爽人人爽vad给| 神马久久中文字幕| 在线影院中文字幕精品| 超薄丝袜脚在线| 国产免费的黄色视频完整版| 久久一区二区福利| 强睡年轻的女老板中文字幕| 久久精品亚洲av久久| 国产成人精品深夜福利视频| 人人做人人爱夜夜爽少妇| 国产高清一区二区老熟女| 久久久久久久高潮av| 人妻精品区一区二区三区| 成人免费av高清在线观看| 3D肉蒲团之极乐宝鉴| 欧美成人精品第一区二区三| 一区二区三区亚洲不卡| 最近中文字幕2018中文字幕 | 久久精品久久免费| 内射国产在线视频| 国产午夜一区二区三区影院| 五月天丁香色| 久久精品一区懂色| 亚洲av午夜福利精品| 午夜一区精品国产亚洲av| 亚洲精品一区中文字幕在线| 97久久草草超级碰碰碰| 国产亚洲人人爱| 99这里都是精品这里有精品| 国产日韩欧美综合一区| 美女双飞高潮直喷水在线播放| 亚洲精品成人av观看孕妇| 蜜臀av日韩精品一区二区| 日韩精品人妻系列一区二区三区| 国产av麻豆久久久久久久| 在线观看欧美国产成年人片| 丁香婷婷激情五月| 人妻激情视频一区二区三区| 国产麻豆精品无色码av| 亚洲国产精品毛片av| 欧美乱大交3| 女人被弄到高潮的视频| 久久国产精品成人| 岛国av毛片网站| 国产又色又爽又免费视频| av亚洲aⅴ久久精品| 黄色日韩一级| 国产av在精品二区| 亚洲高清av片在线观看| 亚洲av国产一级毛片久久| 亚洲va欧美va人人爽久| 日日噜噜夜夜狠狠视频欧美人| 不卡av在线精品| 国产欧美日韩精品A在线观看| 日本一级a爰片| 欧美日韩在线播放网址| 久久看电影久久久久国产| 精品久久人妻av中文字幕成人| 变态挤奶水AV大片| 国产夜视频在线观看| 欧美熟妇丰满XXXXX| 神马久久中文字幕| 熟女毛多熟妇人妻在线视频 | 99视频在线观看97| 99久久精品久久免费观看不卡 | 区二区三美女| 日本少妇高潮喷水免费观看| 一区二区亚洲乱码| 人妻av乱片av出轨九九| 国产日韩欧美一二三级| 国产精品成人在线免费视频| 成人免费av影片| 日日尻夜夜操| 精品国产极品美女在线观看| 污污汅18禁在线永久免费观看| 区二区三美女| 天天爱天天做天天做天天吃中文| 淫片久久久久久| 亚洲资源站无码AV网址| 欧美精品一区二区精油| 亚洲91av在线视频| 人妻少妇精品中文字幕91| 无码无套少妇毛多18P| 欧美一区二区三区六区| 精品国产一区二区久久| 夜夜叫夜夜爽夜夜高潮| 亚洲国产精品久久久久久久简单| 女人与公拘交的视频A片| 伊人 久久 中文字幕| 亚洲另类激情专区小说图片| 中文字幕国产懂色av| 玖草视频在线观看免费| 黄频视频免费看| 日韩人妻精品一区2区三区| 亚洲首页av网| 日本一区二区三区四区五区乱码| 午夜福利片1000无码免费| 亚洲影院一区二区三| 亚洲国产男同同性VIDEOS| 欧美亚洲另类一区二区| 国产精品一区二区免费欧美| 国产又色又爽又黄的免费| 久久国产精品国产精品99| 强睡年轻的女老板中文字幕| 成人高清av| 美国女孩儿成人毛片| 色婷婷精品视频| 夜夜撸日日撸| 久久电影免费福利网址一区| 色婷婷成人激情| 久久人人爽人人片浪潮a| 晚上开车又疼又叫的声音视频| 欧美亚洲另类一区二区| 日韩国产亚洲| 日本XXXX裸体XXXX自慰| 在线乱码av| 亚洲性夜夜综合久久| 日韩在线一区二区三区影视| 线在在线精品视频| 性欧美videosex18嫩| 天天狠天天透天干天天怕处| 最近免费中文字幕大全在线看| 伊人久久大香线蕉亚洲五| 怡红院亚洲色图| av国产熟女丰满熟女| 非洲黑人性xxxx精品又粗又长| 亚洲欧美日韩精美久久亚洲区| 啦啦啦免费观看在线视频| 久久资源中文字幕六区| 亚洲精品国产精品乱| 国产精品一区二区欧美日韩| 国产偷国产偷亚洲高清老太97 | 日韩亚洲区欧美精品| 国产欧美一区精品| 欧美色逼逼视频| 欧美亚洲美腿丝袜| 亚洲av网站在线观看免费| 欧美精品久久久久久久电影| 高级会所人妻互换| 欧美乱妇免费高清视频在线| 成人又黄又爽的免费观看| 午夜福利视频精品| 上别人丰满人妻| 亚洲一区av电影免费| 少妇高潮久久久久久久被弄| 亚洲中文日韩欧美视频| 新年贺岁档表兄妹麻豆| 青青草在线成人观看| 被舔到高潮喷水免费观看| 国产免费永久网址av| 美女MM131爽爽爽免费| 成人免费毛片在线播放视频| 18禁高潮娇喘喷水视频| 狠人色综合网亚洲| 国产亚洲最大av| 两根粗大一前一后好深| 91亚洲精品一区二区在线观看| 国产又粗又硬又爽又猛又黄视频| 久久伊人毛片av| 亚洲成av人影院一区二区| 2019亚洲欧美日韩在线| 香蕉182一区二区| 欧美日韩色一区二区三区| 99热精品这里只有精品| 亚洲欧美日韩成人综合网| av在线播放亚洲精品| 又黄又硬又粗又长又爽的视频| 大香蕉久久精品免费8| 亚洲熟女毛片儿| 97久热在线观看精品| 无码专区亚洲综合系列 | 亚洲av在线观看一区在线| 久久看片天天av免费观看| 欧美国产成人自拍| 99久久精品一区二区三区| 人妻系列无码专区久久五月天| 国产后入免费在线观看| 中文字幕亚洲无线码欧美| 亚洲av日韩av在线| 午夜小视频一区二区三区| 国产日韩欧美视频二区| 老司机午夜福利视频在线| 处膜破AV无码| 久久精品中文一区二区三区 | 老司机午夜视频十八福利| 天天躁久久躁夜夜狠狠躁| 亚洲中文字幕久久一区| 一级无遮挡毛片免费视频| 日韩成人免费电影网站| 中国精品18videosex性中国| 国产成人综合久久精品| 1级片视频免费看| 日本人69xXⅹ69护士| 99精品视频免费全部在线| 国外亚洲成av人片在线观看| 欧美交性又色又爽又黄麻豆| 亚洲欧洲日韩国码av系列天堂| 成人人妻精品一区二区三区| 成人性生交aaa大片| 亚洲欧洲国产综合考虑| 日韩高清一二三区| 丰满高跟丝袜老熟女HD| 人妻无码久久精品| 99久久精品这里只有精品| 欧美黑人巨大videos极品| 手机a级毛片免费观看| 可以免费在线观看的av网址| 欧美激情国产视频一区| 青春草视频在线观看一区| 成人又黄又爽的免费观看| 在线观看av天堂| 国产色/亚洲| 国产三级a午夜电影| 国产精品理论午夜福利91| 日韩精品欧美综合视频一区二区| 99re66久久在热青草| 三级毛片av免费| 亚洲精品国产色婷婷电影| 99久久精品久久免费观看不卡 | 天堂久久亚洲精品| 最好看的中文字幕久久| 一个月的宝宝边吃奶边哭还挣扎| 欧美日本一本线播放三区| 精品久久成人区二区| 国产精品影院99| 综合久久亚洲| 国产一级淫免费播放| 欧美多毛老妇人| 亚洲激情av色| 91一区二区三区在线观看| 999欧美精品一区二区三区| 九九热在线观看精品视频| 亚洲av蜜桃精品| 99亚洲激情| 裸体艳情视频在线观看网站| 亚洲成a人片在线观看视频| 青青草午夜影院| 美女黄a视频免费看| 黄片免费在线播| 新版中文资源天堂韩国在线| 在线日本视频一区| 久久精品国产99久久久 | 高清乱码一卡二卡| 精品三级视频播放| 444欧美日韩国产在线| 日韩亚洲欧美综合| 日韩欧美一区二区三区国产| 免费一级毛片全部免费播放视频 | 激情在线视频精品| 又黄又爽又色18禁网站| 免费成人精品久久| 黄色日韩一级| 亚洲精品无码鲁网午夜| 人妻av福利| 亚洲国内久久精品| 精品中文字幕久久久久久久| 亚洲免费中文字幕一区| 日少妇视频免费| 永久免费视频国产在线播放| 成人av哪个手机网站可以看| 午夜男女爽爽爽免费播放| 青草草免费视频网站| 又爽又黄a免费视频| 五十熟女人妻少妇精品视频| 国产三级精品国产三级人妇在线| 国语对白高潮视频在线免费观看| 日韩精品福利片午夜免费| 婷婷精品国产亚洲麻豆| 性欧美乱妇come| 精品久久国产一区二区三| 国产午夜福利大片免费看| 78m成人免费视频国产av| 亚洲国产色在线日韩| 在线观看av不卡| 口述一次疯狂刺激的交换经历| 插bb射bb福利免费午夜视频| 男人和女人三级黄色片| 成人av免费在线播放网站| 一级做a爰片久久毛| 成人黄色视频网站有哪些| CHINESE农村夫妇双飞| 亚洲av影院网站| 国产精品综合亚洲av久久久网站| 国产在线播放色| 久久精品国产亚洲av香:| 中文字幕人妻第一区| 免费日韩av网站| 亚洲精品嫩草研究院| 天堂а√在线地址中文在线| 亚洲欧美激情综合在线| 大量国产私密保健视频| 中文字幕视频免费看| 亚洲av人人爽人人爽人人爽| 老色鬼在线精品视频| 特级毛片免费视频播放| 精品中文字幕久久久久久久| 欧美激情国产视频一区| 情色av一区二区三区| 少妇的丰满3中文字幕| 国产午夜精品一区二区三区, | 人人妻人人爽人人添夜夜欢视频| 日韩欧美国产一成人av一区二区| 香蕉久久成人国产精品| 没有穿内衣的邻居爆乳| 国产成人精品久久二区二区免费 | 小辣椒导航福利入口| 亚洲欧洲精品一区二区精品久久久| 亚洲第一久久香蕉| 亚洲熟女电影一区二区三区| 中文字幕一区二区人妻| 亚洲欧美日韩快播| 另类天堂网av| 国产成人精品免费视频大全| 亚洲天天久久中文字幕| 99亚洲欧美国产另类在线| 精品免费观看国产一区二区视频| 理论片87福利理论电影| 在线激情亚洲欧美| 人人妻人人添人人爽日韩欧美| 欧美午夜视频在线观看| 午夜久久久久| 国产区精品福利区| 中文字幕精品综合第五页| 最新国语自产精品视频在| 欧美精品一级黄色| 99久久精品视香蕉蕉| 精华精华液有什么用| 男女午夜福利国产| 18禁无遮挡黄网| 午夜福利国产精品久久婷婷69| 在线观看99午夜福利| 亚洲老熟妇熟女| 东京热网站在线观看| 久久这里只有精品美女| 亚洲国产国产av| 国产av香港| 国产精品无码一区二区在线观看 | 久久精品天堂av| 久久自己只精产国品| 久久久精品欧美一区| 日本高清色视频免费观看| 国产色视频综合| 中文在线а√天堂官网| 极品少妇的粉嫩小泬视频| 查理与陈宝莲一级毛片| 午夜福利亚洲电影| 在线99视频在线| 久久精品一区二区三区在| a片太大太长太深好爽在线观看| 麻豆精品人妻一区二区三区| av网国产免费| 一边摸一边爽一边叫床全过程软件| 又黄又嫩又粗又硬视频| 黄色视频可看| 国产精品免费视频能看的| 2021最新国产精品网站 | 国产午夜福利久久精品app| 国产精品乱码一区二区三区四区| w波多野结衣人妻系列| 欧美色逼逼视频| av在线影片| 91伊人影院| 天天射日日干夜夜操| 18av天堂网| 99久久精品国产99| 伊人91影院| A片在线观看全免费| 一区二区毛片免费看| 疯狂伦交一女多男视频| 欧美靠逼视频免费| 三级黄a毛片| 国产精品毛片一区二区排 | JAPANESE55丰满熟妇| 熟女人妻免费视频一区二区三区 | 国产精品一区不卡在线| AV区无码字幕中文色| 欧美高清午夜视频| 高中男生自慰无遮挡网站| 日韩av成人一区二区三区在线| 亚洲区自拍视频在线观看| 亚洲老司机深夜私人影院| 久久精品免费看| 亚洲精品av一区二区三区在线| 亚洲最大AV资源网在线观看| 欧美日韩亚洲xx| 草蜢社区在线观看免费下载| 一区日韩二区人妻精品av| 国产精品激情综合| 亚洲男人手机天堂| 中文字幕久久精品第一页| 人人妻人人添人人爽日韩欧美| 久久精品国产成人综合婷婷免费| 久久精品四虎av| 午夜福利亚洲精品不卡| 邻居少妇很紧毛多水多| 最近亚洲av| 中文字幕91在线视频| 免费观看日本一区二区三区视频| 欧美日韩av天堂| av爱爱一区av| 一本色道久久综合婷婷| 深夜a级毛片| 天堂a中文在线资源库| e午夜精品久久久久久久| 亚洲国产日韩av网站在线观看| 婷婷99狠狠躁天天躁中文字幕| 女人与公拘交的视频A片| 欧美精品亚洲精品一区二区| 国产亚洲午夜高清国产拍精品| a级毛片免费进| 如何把女人舔到高潮| 无码小电影在线观看网站免费 | 中文亚洲爆乳无码专区| 亚洲无在线播放| 欧美精品亚洲一区二区| 岛国在线小视频| 久久大香蕉精品视频| 国产户外露出视频在线观看| 欧美日韩在线一区二区视频| 欧美巨大黑人精品videos| 国产老妇伦国产熟女老妇高清| 亚洲av免费毛片| 久久av电影| 亚洲国产日韩A在线乱码| 99这里只有精品久久| 亚洲国内久久精品| 女人爽得直叫免费视频| 亚洲av免费电影在线| 国产欧美久久久| 极品粉嫩学生国产在线| 久久电影国产精品一区二区三区| 久久aa大香蕉| 亚洲91av在线视频| 黑高跟丝袜一级毛片| 亚洲成人精品自拍| 国产精品综合久久久久久久免费| 国产99精品免费视频| 亚洲欧美日韩精品久久一区| 久久精品久久久久久噜| 2021在线精品自偷自拍无码| 六月丁香婷婷国产| 97碰碰人妻| no results were found for gay| 日本网址在线| 亚洲av不卡一| 99久久婷婷这里只有精品| 高清乱码一卡二卡| 欧美日韩国产在线一区二区三区| 亚洲精品一区二区三区r| 日本蜜桃免费观看mv| 男女亚洲视频在线观看| 骚妻内射免费视频| 69久久夜色精品国产| 人妻婬乱厨房波多野结| 国产精品亚洲精品777| 国产精品国产高清国产aⅴ| 久久最新地址中文字幕| 国产福利在线观看网站| 尤物av不卡在线观看| 18成人午夜网站在线观看| 久久不见久久见免费影院播放 | 成人黄色免费在线视频| 国产精品一区二区三区av.com| 大又大粗又爽又黄的a级毛片 | 少妇看A片自慰| 亚洲图片av一区| 高清国产天堂在线bt免费资源| 国产男女猛烈无遮挡91| 免费a级毛片:| 国产在线观看无码的免费网站| 久久成人一区二区精品| 青草久久97| 国产午夜精品一区二区电影| 黑人巨大精品一区二区在线| 朝鲜少妇漂亮毛茸茸| a级特级毛片在线播放| 国一区二区三区在线| 国产福利91精品一区| 中国国产av一级| 日韩综合一二三区| 欧洲毛片av片一级毛片| 经典人妻一区二区| 国产伦一区二区三区视频涩爱| 国产av电影区二区三| 久久久久久久久大av| 色哟哟在线永久网站| 性AV无码天堂VR专区| 要看免费黄色毛片| 色婷婷亚洲综合网| 午夜福利免费看黄片| 天堂va亚洲va欧美va| 欧美日韩视频高清一区二区三区二| 中国男男自慰GAY片免费观看| 精品一区二区三区乱码在线观看| 亚洲av男人的天堂日韩| 国产精品成人亚洲毛片| 99久久精品费精品蜜臀av| 精品一区久久三区二区| 精品人妻一区二区三区麻豆 | 精品人妻精品视频| 亚洲av国产精品午夜福利| 日本欧美国产在线播放| 精品亚洲成av人在观看| 被黑人猛烈30分钟视频| 成人做爰视频ww网站| 久久99精品国语久久久| av网站在线观看亚洲| 宝贝对着摄像头自己做| 国产精品亚洲综合久久系列| 正在播放国产美女福利vip| 中文字幕亚洲三区| av野外在线| 伊人精品成人一区二区三区四区| 蜜臀av日韩精品一区二区| 91性高潮久久久久久久久久久久| 亚洲午夜福利视频一区二区三区| 国产av自拍高清一区| 熟女电影av网| 高清不卡的av网站| 日本a视频在线| GAY男男自慰免费播放| 秋霞av午夜| 日本强伦姧护士MMM| 五月激激激综合网亚洲| 成 人3D动漫在线观看网站| 99午夜福利一区二区| 国产二区激情在线观看| 亚洲一区亚洲精品久久| 午夜鲁丝在线观看| 精品欧美一区二区三区四区五区| 女人高潮一级毛片| 久久视频伦理片| 欧美91激情亚洲| 国产精品欧美日韩三区四区| 国产亚洲精品网站在线视频| 国产精品午夜视频一区| 免费av尤物网址| 国产影片中文字幕久久| 青青草午夜影院| 可以免费在线观看的av网址| 免费在线观看不卡av网站| 国内精品欧美激情| 在线观看国产黄色片| 天堂а√在线地址中文在线| 人妻熟女av久久| 亚洲欧美日韩精品在| 黄爽无遮挡网站| 国产成人a亚洲精v品| 99热思思这里只有精品| 自拍偷区色综合| 国模生殖欣赏人体337| 新年贺岁档表兄妹麻豆| 亚洲综合色激情五月| 女人的超长巨茎人妖在线视频| 熟妇人妻中文AV无码| 国内精品美女久久久久久| 久久久久性生活片| 丰满人妻一区二区三区色91| 国产免费精品高清| 一本色道久久88综合亚州精品| www.av天堂色| 人人妻人人添人人爽欧美图片| 亚洲经典成人影院| 亚洲熟女老熟女自拍av| 黄片免费观看视频99| 2019亚洲欧美日韩在线| 亚洲最大AV资源网在线观看| 国产乱理伦片在线观看夜| GAY男同GV网站播放免费| 亚洲欧美综合区自拍另类| 444欧美日韩国产在线| 蜜臀av国产亚洲欧美xxx| 国产又爽又黄的激情精品视频| 久久99精品免费无毒不卡| 男女那个视频网站在线观看| 国产无遮挡又黄又爽久久| 欧美色播视频网| 国产一级高清av| 九九九中文无码AV在线播放| 成人在线观看免费一区二区| 成人二区视频| 国产精品乱码一区二区三区软件| 国产亚洲99精品| 欧美精品在线免费播放| 欧美综合在线一区二区三区| 各类成人av片在线观看| 亚洲深深色噜噜狠狠爱网站| 可以直接观看的一级毛片| 高清午夜精品一区二区三区| 亚洲av免费观看久久| 亚洲中文免费在线播放| 丝袜av网址| h动画网站在线观看| 91欧美人妻精品激情| 亚洲乱码中文字幕97精品| 亚洲成a人片8888一在线观看| 中文字幕在线91| 怡红院亚洲免费| 蜜芽亚洲AV无码精品色午夜| 色欲色香天天天综合网站| 最近中文字幕免费视频大全| 亚洲成人免费久久电影| 日韩高清av免费在线观看| 欧美三级日韩三级亚洲三级| 欧美国产成人精品| 亚洲avzai| 人妻少妇区一区二区三| 欧美亚洲狠狠狠| 国产精品亚洲第一区在线 | 日少妇视频免费| 久久—国产精品—爽—久久| 亚洲中文无码AV永久主页不卡| 日日噜夜夜操| 极品少妇高潮喷水抽搐91| 国产又爽又黄又免费观看视频| 在线播放av的网址| 最新国产高潮在线观看| 免费一级毛片正在播放| 精品国产亚洲一区观看| 免费av高潮喷水| 亚洲av少妇一区二区在线观看| 久久69国产精品久久69免费看| 久久成年精品视频| 国产亚洲一区二区三区av| 奈奈美被公侵犯到怀孕在线观看| 欧美多毛熟女视频| 香港三级电影国产精品| 欧美成人aa在线观看视频网站| 日本一卡二卡三卡四卡2021| 亚洲欧美洲成人1区二区| 91丝袜美腿诱惑| 99精品推荐免费视频| 亚洲AV无一区二区三区久久 | 久久99综合婷婷国产二区高清| 麻豆av激情在线观看| 永久免费av大片| 91精品国产国语对白视| 国产高清乱理伦片中文| 三级a毛片免费| 91po熟女| 久久天堂av在线| 国产日韩欧美亚洲二区| 成人精品99| 中文天堂在线一区删除| 女女同性AV片在线观看免费| 日韩欧美亚洲区| 可在线免费观看av网站| 国产亚洲精品av久久综合麻豆| 久久久久精品国产电影| 国产有码高清视频在线观看| 久久亚洲国产精品成人| aa级毛片免费网站| 久久一级毛片孕妇| 成人免费无码大片A毛片不卡| 波多野结衣av高清一区二区| 黑人操到高潮| 97在线观看免费视频| 久久久精品国产亚洲av一二三 | 波多野结衣四区| 国产美女被侮辱网站| 小14萝裸体洗澡视频免费网站| 日韩伦理电影一区二区三区| 好色妞色视频| 无码精品国产一区二区三区免费| 国产av零色中文网| 午夜精品av福利片| 中文字幕人妻一区二区免费| 精品亚洲999| 4399日本高清完整版在线观看 | 欧美国产亚洲浪潮| 88av影院首页| 亚洲码国产岛国毛片在线| CHINESE激烈高潮HD| 国产欧美一区色二区三区| 日韩在线看大片| 成人特级黄色片久久久久久久| a级毛片免费在线看| 晚上开车又疼又叫的声音视频| av天堂亚洲图| 亚洲精品国产有码| 久久黄色a级毛片| 欧美亚洲狠狠狠| 男男chinese男男资源| 亚洲精品一区国产二区三区四区| 先锋77xfplay色资源网站| 高潮直喷水呻吟| 手机中文字幕永久在线| 日本最新免费av| 欧美人与拘性视交免费| 99久久精品午夜| 一本久久综合精品| 美女露出尿口喷水无遮挡| 成人国产av精品久久电影| 美女国产精品理论片影院| 色婷婷色丁香| 精品欧美日韩国产一区| 毛片高清不卡| 久久国产男人的天堂| 国产黄a三级三级三级人| av中文字幕第一区| 国产黄色av片免费播放 | 粗暴蹂躏人妻av一区二区三区| 欧美精品成人a多人在线观看| 3344成人在线免费视频| 9420高清完整版在线观看| 免费黄色片一级毛片| 大号毛BBW厕所偷窥| 国产情侣色在线| 青柠在线观看免费高清完整版| 欧美一级日韩不卡播放免费| 第一次处破女18分钟| 丰满少妇猛男进入| 草蜢社区在线观看免费下载| 人人爽人人爽人人爽人人片av| 国产精品亚洲w码日韩中文| av网站网页在线观看| 张筱雨两腿玉门打开图| 久久日韩精品一区二区| 亚洲日本不卡专区| 亚洲av高清在线观看| 欧美日韩亚洲高清精品| 看a级黄片毛片| 亚洲欧美激情久久| 少妇久久精品一区二区夜夜嗨| 午夜小视频一区二区三区| 日韩欧美亚洲91| 精品国产女人久久久久| 再快点再深点我要高潮了视频 | 66美女精品网站| 娇小的AV|色学生| 成年大片40分钟免费视频播放| 日本网站在线免费看| 亚洲综合色另类| 被男子脱掉内裤的美女视频| 国产精品嫩草影院一二三区入口| 又黄又爽又色的视频| 午夜在线不卡精品国产| 国产成人秘史在线播放| 理论片87福利理论电影| 国产精品视频福利一区| 亚洲av国产一级毛片久久| 一区二区三区+在线播放| 风韵熟妇bbbbxxxx| 在线免费亚洲中文字幕| 在线视频精品999| 日韩久久人妻网| 美女下体隐私视频在线网站观看 | 公交车猛烈进出婷婷| 久久五月婷婷狠狠| 精品少妇一区二区三区视频日本电影 | 特级毛片A级毛片免费播放| 91在线电影免费观看| 一级毛片久久久久久久久女 | chinese国产HD中国情侣| 老司机免费午夜福利视频| 欧美一区二区三区老太婆性生活| 国产内射在线激情一区| 国产高清精品免费精品| 精品国产一区二区三区四区导航| 风骚老熟女视频| 99久久嫩草影院| 久久99热这里只有精品| 熟女人妻の波多野结衣在线| 国产乱久久亚洲国产精品| 天天躁日日躁狠狠躁欧美老妇 | 东方亚洲欧美色| 免费网站看v片在线毛| 精品亚洲在线播放| 99热这里只有的精品最新地址| 欧美日韩在线精品一区二区三区激情 | 欧美一级二级三级成人| A片试看120分钟做受视频大全| 日韩亚洲欧美三区中文字幕| 国产精品一区二区hb在线观看| 日韩高清网页| 国产福利免费在线视频观看| 女国产精品视频一区| 黑人欧美特级aaaaaa片| 国内亚洲精彩自拍视频| 午夜福利国产在线观看2| 亚洲精品黄片大全| 精品人妻一区二区四区| 亚洲精品一区二区91在线| 国产精品欧美1区| 久久婷婷人人爽人人干人人爱 | 亚洲婷婷成人网| 99久久精品久久久久清纯| 日韩久久九九精品一本| 人人做人人爱夜夜爽少妇| 日韩性色av一区二区三区| av毛片精品99久久av| 16女下面流水不遮视频| 国产精品成人欧美一区| 最新中文字幕永久视频播放| 伊人99热这里只有精品| 人妻少妇精品久久| 亚洲精品影院一区二区三区| 三上悠亚国产精品一区二区三区| 天天躁狠狠躁夜夜躁狠狠躁| 翁公和在厨房猛烈进出| 亚洲国产欧美乱码| 成人国产亚洲精品一区| 女同久久另类99精品国产91| 亚洲AV无码专区色爱天堂| 精品欧美日韩国产一区| 99这里只有精品久久| 久久久久久久高潮av|