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

    基于MATLAB實(shí)現(xiàn)的AR模型功率譜估計(jì)

    2017-09-07 06:37:46劉明曉王旭光
    電子設(shè)計(jì)工程 2017年17期
    關(guān)鍵詞:譜估計(jì)階次協(xié)方差

    劉明曉,王旭光

    (1.西安鐵路職業(yè)技術(shù)學(xué)院 陜西 西安 710014;2.西安航天自動(dòng)化股份有限公司 陜西 西安 710014)

    基于MATLAB實(shí)現(xiàn)的AR模型功率譜估計(jì)

    劉明曉1,王旭光2

    (1.西安鐵路職業(yè)技術(shù)學(xué)院 陜西 西安 710014;2.西安航天自動(dòng)化股份有限公司 陜西 西安 710014)

    現(xiàn)代功率譜估計(jì)和經(jīng)典功率譜估計(jì)是兩種常用的功率譜估計(jì),同時(shí)也是分析隨機(jī)信號的常用方法。本文詳細(xì)介紹了現(xiàn)代功率譜估計(jì)中有關(guān)AR模型參數(shù)的功率譜估計(jì),具體包括自相關(guān)算法、Burg算法、協(xié)方差算法以及改進(jìn)的協(xié)方差算法,在此基礎(chǔ)上又對四種不同功率譜估計(jì)方法的性能指標(biāo)進(jìn)行了比較分析。在MATLAB仿真軟件平臺上對AR模型參數(shù)的四種不同功率譜估計(jì)算法進(jìn)行了仿真,同時(shí)對功率譜估計(jì)結(jié)果進(jìn)行了分析比較并得到了預(yù)期的譜估計(jì)效果。最后從實(shí)際應(yīng)用角度出發(fā)討論了AR模型參數(shù)不同功率譜估計(jì)算法的特點(diǎn),以便在應(yīng)用中能夠選擇合適的功率譜估計(jì)算法。

    功率譜估計(jì);AR模型;算法;MATLAB

    在頻譜分析中功率譜估計(jì)是常被采用的一種重要方法,功率譜估計(jì)廣泛應(yīng)用于人們的日常生活、聲納、生物醫(yī)學(xué)、雷達(dá)以及農(nóng)業(yè)活動(dòng)等不同領(lǐng)域的信號處理過程中[1-2]。如何從頻帶較寬的信號中檢測出頻帶較窄的信號也是功率譜估計(jì)在信號處理中的重要應(yīng)用[3]。此外功率譜估計(jì)是用有限長的數(shù)據(jù)來對信號進(jìn)行功率譜估計(jì),是分析隨機(jī)信號時(shí)常被采用的一種重要方法。

    為了能夠良好地解決經(jīng)典功率譜估計(jì) (包括自相關(guān)法和周期圖法)中方差性能不好以及分辨率較低等問題提出了現(xiàn)代功率譜估計(jì)[4],其涉及的領(lǐng)域以及學(xué)科相當(dāng)廣泛且內(nèi)容比較豐富[5]?,F(xiàn)代功率譜估計(jì)包括非參數(shù)模型估計(jì)和參數(shù)模型估計(jì)兩種,前者有最小方差法,多分量的MUSIC方法等;后者有AR模型、MA模型、ARMA模型、PRONY指數(shù)模型等[6-7]。文中對現(xiàn)代功率譜估計(jì)中有關(guān)AR模型參數(shù)的幾種典型求解算法進(jìn)行分析比較,同時(shí)在MATLAB仿真軟件平臺上對AR模型參數(shù)中不同算法的功率譜估計(jì)進(jìn)行了仿真,最后從實(shí)際應(yīng)用的角度出發(fā)討論了不同功率估計(jì)算法的特點(diǎn),以便在實(shí)際應(yīng)用中能夠選擇最合適的功率譜估計(jì)算法[8]。

    1 現(xiàn)代功率譜估計(jì)

    現(xiàn)代功率譜估計(jì)又稱參數(shù)譜估計(jì)[9],是先對數(shù)據(jù)進(jìn)行觀察后再進(jìn)一步對參數(shù)模型進(jìn)行估計(jì),然后根據(jù)得到參數(shù)模型輸出的功率來對信號進(jìn)行譜估計(jì)。采用現(xiàn)代功率譜估計(jì)能夠較好地改善經(jīng)典譜估計(jì)中方差性能差以及分辨率低等問題[10]?,F(xiàn)代功率譜估計(jì)中AR模型是應(yīng)用較多的一種模型。

    1.1 AR模型的建立

    AR模型是全極點(diǎn)模型同時(shí)又稱為自回歸模型[11]??捎霉剑?)所示的差分方程來描述:

    其中u(n)是一個(gè)白噪聲序列,其均值為0、方差為 σ2,p 是 AR 模型的階數(shù),ap(i),i=1,2,…,p。AR 模型系統(tǒng)的轉(zhuǎn)移函數(shù)H(z)為:

    由公式(1)和公式(2)可得到AR模型參數(shù)的功率譜估計(jì)計(jì)算公式:

    AR模型的參數(shù)和x(n)自相關(guān)函數(shù)有如下的關(guān)系:

    將上式寫成矩陣的形式:

    式(5)為AR模型的正則方程,又稱尤拉-沃克(Yule-Walker)方程[12]。

    1.2 AR模型參數(shù)估計(jì)的典型算法

    1.2.1 自相關(guān)法

    自相關(guān)法是AR模型參數(shù)估計(jì)中較為簡單的一種功率譜估計(jì)方法。按照模型階數(shù)由小到大的順序進(jìn)行計(jì)算,先計(jì)算階次 m=1 時(shí)的系數(shù) am(k)=a1(1)和,再計(jì)算 m=2 時(shí)的 a2(1),a2(2)和,按此方法依次計(jì)算到階次 m=p 時(shí)的 ap(1),ap(2),…,ap(p),及,當(dāng)?shù)木葷M足要求時(shí)即可停止。遞推公式為:

    1.2.2 Burg算法

    Burg算法進(jìn)行功率譜估計(jì)時(shí)其前后項(xiàng)序列誤差的前后都不加窗,使用Levinson-Durbin遞推可快速求解AR模型參數(shù)的系數(shù)。

    Burg算法與自相關(guān)算法不同之處在于使序列x(n)的前后項(xiàng)功率預(yù)測誤差之和:

    再根據(jù)Levinson-Durbin進(jìn)行遞推得出AR模型系數(shù)如下:

    Burg算法不需先計(jì)算自相關(guān)函數(shù)進(jìn)而計(jì)算速度有所提高同時(shí)計(jì)算也相對簡單。與自相關(guān)算法相比分辨率也有所提高。但在處理白噪聲加正弦信號時(shí)功率譜線會(huì)出現(xiàn)分裂現(xiàn)象。

    1.2.3 協(xié)方差算法

    協(xié)方差算法與其他算法的不同之處在于預(yù)測誤差功率求和的上下限取值范圍不同。在協(xié)方差算法區(qū)間[0,N-1]外 x(n)樣本不為 0,故表達(dá)式中的 x(n-k)總是落在區(qū)間[0,N-1]中,因此預(yù)測誤差功率求和的上下限必然會(huì)落在[p,N-1]之間。得到的矩陣是半正定矩陣,且不具有Toeplitz性質(zhì),故不能采用Levinson-Durbin遞歸算法進(jìn)行參數(shù)求解,同時(shí)得到的AR模型也不穩(wěn)定。

    1.2.4 改進(jìn)協(xié)方差算法

    采用改進(jìn)的協(xié)方差算法進(jìn)行功率譜估計(jì)時(shí),為了使前后項(xiàng)功率誤差之和最小其前后項(xiàng)都不加窗,得到的協(xié)方差矩陣不是Toeplitz矩陣,故不能用Levinson遞推算法對模型參數(shù)進(jìn)行求解。Marple在1980年提出了改進(jìn)協(xié)方差方程求解的快速算法,該算法能夠大大提高功率譜估計(jì)的性能。

    1.3 用于功率譜估計(jì)的數(shù)據(jù)

    文中采用由復(fù)數(shù)噪聲加上4個(gè)復(fù)正弦組成的數(shù)據(jù)對功率譜進(jìn)行估計(jì)。其中歸一化頻率分別是f1=0.15,f2=0.16,f3=0.252,f4=-0.16。 該數(shù)據(jù)的 128 點(diǎn)復(fù)序列如圖1所示。

    圖1 功率譜估計(jì)數(shù)據(jù)

    2 AR模型功率譜估計(jì)的MATLAB仿真

    借助MATLAB對AR模型功率譜進(jìn)行估計(jì)時(shí),上述4種算法分別采用pyulear函數(shù)、pburg函數(shù)、pcov函數(shù)和pmcov函數(shù)來完成功率譜估計(jì)。分別用這4個(gè)函數(shù)對圖1中的隨機(jī)序列進(jìn)行功率譜估計(jì)便可得到相應(yīng)的功率譜線。圖2~4是自相關(guān)算法得出的階次p分別為10,20和30時(shí)的AR模型功率譜曲線。

    圖2 階次p=10時(shí)自相關(guān)算法功率譜估計(jì)

    從圖2~4的功率譜估計(jì)圖形可以得出,在階次p較低時(shí) (如圖2所示)功率譜估計(jì)的分辨率不理想。當(dāng)階次p提高到30時(shí),在頻率f1和f2處兩個(gè)正弦信號剛好被分離開,在頻率f3和f4處也能夠分離出兩個(gè)正弦信號。在自相關(guān)算法中由于對序列的前后項(xiàng)加窗降低了分離效果,數(shù)據(jù)越短分離效果越不好。

    圖3 階次p=20時(shí)自相關(guān)算法功率譜估計(jì)

    圖4 階次p=30時(shí)自相關(guān)算法功率譜估計(jì)

    圖5~圖6是Burg算法對隨機(jī)數(shù)列信號進(jìn)行功率譜估計(jì)時(shí)的功率譜曲線。

    圖5 階次p=10時(shí)Burg算法功率譜估計(jì)

    圖6 階次p=15時(shí)Burg算法功率譜估計(jì)

    比較分析圖5和圖6可以得出在階次p=15時(shí)得到了較好的譜估計(jì)。另外對比圖4和圖5可以看到Burg算法的分辨率明顯優(yōu)于自相關(guān)算法,在階次較低時(shí)(如圖5所示)也能較好地分辨出間隔較小的頻率點(diǎn)。

    圖7為階數(shù)p=10時(shí)協(xié)方差算法對隨機(jī)數(shù)列信號進(jìn)行估計(jì)得到的功率譜曲線。

    圖7 階次p=10時(shí)協(xié)方差算法功率譜估計(jì)

    從圖7可得出功率譜估計(jì)在信號源頻率f1=0.15,f2=0.16,f3=0.252,f4=-0.16 處功率譜線相對較窄,在其他頻率處功率譜線起伏劇烈。采用協(xié)方差算法對功率譜進(jìn)行估計(jì)能夠較真實(shí)地反映信號實(shí)際模型。

    圖8~圖9是階數(shù)p=10和p=15時(shí)改進(jìn)的協(xié)方差算法得出的AR模型功率譜估計(jì)曲線。

    圖8 階次p=10時(shí)改進(jìn)協(xié)方差算法功率譜估計(jì)

    圖9 階次p=15時(shí)改進(jìn)協(xié)方差算法功率譜估計(jì)

    比較圖7~圖9可以得出兩種協(xié)方差算法的功率譜估計(jì)大致相同,但改進(jìn)的協(xié)方差算法在信號頻率處的峰值更加尖銳、更加突出且更有利于辨識。分析圖8和圖9可以得出在階次較高時(shí)能得到非常滿意的功率譜估計(jì)曲線。

    3 結(jié) 論

    文中介紹了現(xiàn)代功率譜估計(jì)中AR模型參數(shù)的幾種功率譜估計(jì)算法。其中自相關(guān)算法需要對前項(xiàng)預(yù)測的誤差序列前后加窗,加窗降低了自相關(guān)算法的分辨率,數(shù)據(jù)越短分辨率越低。與自相關(guān)算法相比Burg算法有較好的分辨率,但有時(shí)會(huì)出現(xiàn)譜估計(jì)曲線分裂現(xiàn)象。改進(jìn)的協(xié)方差算法性能最好,但計(jì)算與編程都相對繁瑣。在對功率譜進(jìn)行估計(jì)時(shí),需要依據(jù)實(shí)際數(shù)據(jù)序列選擇更合適的算法進(jìn)行功率譜估計(jì)。

    [1]鄧澤懷,劉波波,李彥良.常見的功率譜估計(jì)方法及其Matlab仿真[J].電子科技,2014(2):50-52.

    [2]李軍,王凱.三種現(xiàn)代功率譜估計(jì)方法性能研究[J].科技信息,2010(30):554-555.

    [3]李文超.現(xiàn)代功率譜估計(jì)在設(shè)備故障診斷應(yīng)用[J].傳感器世界,2012(2):25-29.

    [4]王春興.基于Matlab實(shí)現(xiàn)現(xiàn)代功率譜估計(jì)[J].現(xiàn)代電子技術(shù),2011,34(16):65-67.

    [5]范澤華,白鐵成.基于改進(jìn)差分演化的高效圖像分割算法[J].現(xiàn)代電子技術(shù),2016,39(14):87-91.

    [6]楊蓉.多普勒雷達(dá)測速系統(tǒng)設(shè)計(jì)及信號處理方法研究[D].武漢:華中科技大學(xué),2012.

    [7]汪勝波.功率譜估計(jì)在寬帶ADCP信號檢測中的研究與應(yīng)用[D].青島:中國海洋大學(xué),2011.

    [8]閆慶華,程兆剛,段云龍.AR模型功率譜估計(jì)及Matlab實(shí)現(xiàn) [J].計(jì)算機(jī)與數(shù)字工程,2010,38(4):154-156.

    [9]郭會(huì)娟.基于AR參數(shù)模型和小波分析的電機(jī)故障診斷[D].鄭州:華北水利水電大學(xué),2013.

    [10]崔曉榮.電力系統(tǒng)間諧波的檢測算法研究與實(shí)現(xiàn)[D].成都:西華大學(xué),2013.

    [11]張峰,伊麗,石現(xiàn)峰.振動(dòng)信號的最小方差譜估計(jì)算法[J].計(jì)算機(jī)技術(shù)與發(fā)展,2015(10):76-79.

    [12]袁華,趙越,藺永政,等.基于最小方差的功率譜估計(jì)方法的研究[J].?dāng)?shù)字通信世界,2015(9):254-255.

    The realization of AR model power spectrum estimation based on MATLAB

    LIU Ming-xiao1,WANG Xu-guang2
    (1.Xi'an Railway Vocational&Technical Institute,Xi'an 710014,China; 2.Xi'an Aerospace Automation CO,.LTD.,Xi'an 710014,China)

    Power spectrum estimation can be divided into modern spectral estimation and classical spectral estimation,it is an method for analyzing random signal.It describes the autocorrelation algorithm,Burg algorithm,covariance algorithm and improved covariance algorithm of AR model parameters in Modern Spectral Estimation,and their perform indicators are analyzed in this paper.These methods for AR model parameter algorithm are implemented by MATLAB,the results of simulation are analyzed and get a better spectrum estimation.Moreover,the advantages and disadvantages of these methods are discussed from the experimental view so as to make reasonable selection in practical work.

    PSD estimation; AR model; algorithm;MATLAB

    TN015

    A

    1674-6236(2017)17-0129-04

    2016-07-23稿件編號:201607167

    劉明曉(1987—),女,河南焦作人,碩士,助教。研究方向:電力系統(tǒng)及其自動(dòng)化。

    猜你喜歡
    譜估計(jì)階次協(xié)方差
    階次分析在驅(qū)動(dòng)橋異響中的應(yīng)用
    基于Vold-Kalman濾波的階次分析系統(tǒng)設(shè)計(jì)與實(shí)現(xiàn)*
    基于齒輪階次密度優(yōu)化的變速器降噪研究
    不確定系統(tǒng)改進(jìn)的魯棒協(xié)方差交叉融合穩(wěn)態(tài)Kalman預(yù)報(bào)器
    一種基于廣義協(xié)方差矩陣的欠定盲辨識方法
    經(jīng)典功率譜估計(jì)方法的研究
    高維隨機(jī)信號THREE功率譜估計(jì)及其仿真
    Welch譜估計(jì)的隨機(jī)誤差與置信度
    縱向數(shù)據(jù)分析中使用滑動(dòng)平均Cholesky分解對回歸均值和協(xié)方差矩陣進(jìn)行同時(shí)半?yún)?shù)建模
    關(guān)于協(xié)方差的U統(tǒng)計(jì)量檢驗(yàn)法
    男女下面进入的视频免费午夜| 91久久精品国产一区二区成人 | 亚洲精品久久国产高清桃花| 亚洲一区高清亚洲精品| 真人做人爱边吃奶动态| 两个人的视频大全免费| 国产精品一及| 深爱激情五月婷婷| 亚洲欧美日韩高清专用| 久久久久精品国产欧美久久久| 欧美+日韩+精品| 国产免费一级a男人的天堂| 草草在线视频免费看| 欧美性感艳星| 欧美日韩福利视频一区二区| av片东京热男人的天堂| 亚洲av成人精品一区久久| 成年版毛片免费区| 男女下面进入的视频免费午夜| 欧美一区二区亚洲| 久久久国产成人免费| www.熟女人妻精品国产| 黄色日韩在线| 久久精品国产自在天天线| 免费看美女性在线毛片视频| 精品无人区乱码1区二区| 精品一区二区三区av网在线观看| 在线观看免费视频日本深夜| 国产精品av视频在线免费观看| 最近视频中文字幕2019在线8| 日韩亚洲欧美综合| 怎么达到女性高潮| 女同久久另类99精品国产91| 中文字幕人成人乱码亚洲影| 亚洲av电影在线进入| 99在线视频只有这里精品首页| www.熟女人妻精品国产| 搡老熟女国产l中国老女人| 中文字幕av成人在线电影| 久久精品国产清高在天天线| 欧美xxxx黑人xx丫x性爽| 欧美丝袜亚洲另类 | 亚洲不卡免费看| 国产精品亚洲美女久久久| 国产高清视频在线播放一区| 亚洲精品美女久久久久99蜜臀| 久久天躁狠狠躁夜夜2o2o| 婷婷丁香在线五月| 少妇高潮的动态图| 精品久久久久久,| 一个人免费在线观看的高清视频| 午夜激情福利司机影院| 好看av亚洲va欧美ⅴa在| 亚洲电影在线观看av| 免费在线观看影片大全网站| 成人av一区二区三区在线看| 亚洲中文字幕一区二区三区有码在线看| 亚洲国产欧美网| 天天一区二区日本电影三级| 成人av一区二区三区在线看| 3wmmmm亚洲av在线观看| 久久99热这里只有精品18| 九色成人免费人妻av| 91久久精品国产一区二区成人 | 夜夜躁狠狠躁天天躁| 欧美日韩国产亚洲二区| 国产精品亚洲一级av第二区| 亚洲欧美日韩无卡精品| 欧美日韩中文字幕国产精品一区二区三区| 免费观看人在逋| 久久久色成人| 国产v大片淫在线免费观看| 国产成人啪精品午夜网站| 午夜免费男女啪啪视频观看 | av欧美777| 久久这里只有精品中国| 国产精品av视频在线免费观看| 亚洲狠狠婷婷综合久久图片| 国产亚洲av嫩草精品影院| 在线看三级毛片| 国产中年淑女户外野战色| 亚洲欧美日韩高清专用| 日本撒尿小便嘘嘘汇集6| 麻豆成人午夜福利视频| 欧美日韩福利视频一区二区| 精品一区二区三区视频在线 | 欧美最新免费一区二区三区 | 看片在线看免费视频| 天堂av国产一区二区熟女人妻| 99精品在免费线老司机午夜| 国产精品99久久久久久久久| 一进一出抽搐gif免费好疼| 在线观看一区二区三区| 国产一区二区激情短视频| 国语自产精品视频在线第100页| 99久国产av精品| 好男人在线观看高清免费视频| x7x7x7水蜜桃| 两人在一起打扑克的视频| 国产av一区在线观看免费| xxx96com| 亚洲中文日韩欧美视频| 国产免费男女视频| 国产免费男女视频| 国产黄片美女视频| 国产老妇女一区| 99国产极品粉嫩在线观看| 日本 av在线| 国内揄拍国产精品人妻在线| 少妇人妻一区二区三区视频| 国产高清激情床上av| 久久久久久久久大av| 日日干狠狠操夜夜爽| 少妇人妻精品综合一区二区 | 老熟妇乱子伦视频在线观看| 日本成人三级电影网站| 变态另类成人亚洲欧美熟女| 在线观看66精品国产| 九色国产91popny在线| 美女高潮的动态| 熟女人妻精品中文字幕| 桃红色精品国产亚洲av| 日韩欧美免费精品| av视频在线观看入口| 色精品久久人妻99蜜桃| 好男人在线观看高清免费视频| 国产单亲对白刺激| 欧美中文日本在线观看视频| 亚洲av不卡在线观看| 中出人妻视频一区二区| 国产v大片淫在线免费观看| 国产高清videossex| 精品人妻一区二区三区麻豆 | 国产成人av教育| 女警被强在线播放| 免费在线观看亚洲国产| 国产成人aa在线观看| 国产中年淑女户外野战色| 桃色一区二区三区在线观看| 午夜福利成人在线免费观看| 成熟少妇高潮喷水视频| 黑人欧美特级aaaaaa片| 久久人妻av系列| 一级毛片女人18水好多| 韩国av一区二区三区四区| 国产精品98久久久久久宅男小说| 午夜福利高清视频| 3wmmmm亚洲av在线观看| 男人的好看免费观看在线视频| 天堂网av新在线| 91久久精品电影网| 最近最新中文字幕大全免费视频| 天堂网av新在线| av中文乱码字幕在线| 男女之事视频高清在线观看| 校园春色视频在线观看| 国产极品精品免费视频能看的| 亚洲精品色激情综合| 别揉我奶头~嗯~啊~动态视频| 一级毛片高清免费大全| 亚洲国产精品999在线| 狠狠狠狠99中文字幕| 亚洲久久久久久中文字幕| 免费av不卡在线播放| 久久中文看片网| 国产欧美日韩一区二区精品| 国产成+人综合+亚洲专区| eeuss影院久久| 一本久久中文字幕| 国产伦在线观看视频一区| 搡女人真爽免费视频火全软件 | 国产亚洲精品av在线| 丁香欧美五月| 给我免费播放毛片高清在线观看| 久久精品亚洲精品国产色婷小说| 一进一出好大好爽视频| 热99re8久久精品国产| 少妇熟女aⅴ在线视频| 国产精品久久久久久精品电影| 欧美成狂野欧美在线观看| 国产伦精品一区二区三区四那| 免费看a级黄色片| 夜夜爽天天搞| 日韩大尺度精品在线看网址| 搡老岳熟女国产| 啪啪无遮挡十八禁网站| 成人精品一区二区免费| 男女做爰动态图高潮gif福利片| 真人一进一出gif抽搐免费| 久久久久久久久中文| 久久精品国产自在天天线| a在线观看视频网站| 亚洲人成网站高清观看| 亚洲国产精品sss在线观看| 久久久久久久午夜电影| 俺也久久电影网| 别揉我奶头~嗯~啊~动态视频| 美女 人体艺术 gogo| 国产真实乱freesex| 一级黄片播放器| 亚洲不卡免费看| 免费观看精品视频网站| 国产高清三级在线| 一本综合久久免费| 日本黄色视频三级网站网址| 亚洲真实伦在线观看| 男人舔女人下体高潮全视频| 国产蜜桃级精品一区二区三区| 少妇的逼好多水| 少妇裸体淫交视频免费看高清| 18禁裸乳无遮挡免费网站照片| 中文字幕精品亚洲无线码一区| 天堂影院成人在线观看| 日本黄色视频三级网站网址| 久久久精品欧美日韩精品| 成年免费大片在线观看| 欧美日韩综合久久久久久 | 午夜免费激情av| 免费在线观看日本一区| 男女午夜视频在线观看| 性色av乱码一区二区三区2| 成熟少妇高潮喷水视频| 变态另类成人亚洲欧美熟女| 亚洲精品色激情综合| 精品一区二区三区视频在线 | 哪里可以看免费的av片| 少妇的丰满在线观看| 国产精品1区2区在线观看.| 草草在线视频免费看| 最近最新免费中文字幕在线| 尤物成人国产欧美一区二区三区| 国产亚洲精品一区二区www| 国产伦一二天堂av在线观看| 国语自产精品视频在线第100页| 国产精品自产拍在线观看55亚洲| 有码 亚洲区| 日韩欧美一区二区三区在线观看| 久久久国产精品麻豆| 久久6这里有精品| 午夜激情福利司机影院| 一本精品99久久精品77| 两性午夜刺激爽爽歪歪视频在线观看| 国产一级毛片七仙女欲春2| 99久久精品热视频| 亚洲av一区综合| 在线免费观看不下载黄p国产 | 一个人看视频在线观看www免费 | 国产中年淑女户外野战色| 亚洲av五月六月丁香网| 床上黄色一级片| 国产97色在线日韩免费| 1000部很黄的大片| 亚洲aⅴ乱码一区二区在线播放| 精品一区二区三区av网在线观看| 国产欧美日韩一区二区三| 99精品久久久久人妻精品| 少妇高潮的动态图| 51午夜福利影视在线观看| 怎么达到女性高潮| 欧美黄色片欧美黄色片| 亚洲不卡免费看| 男女午夜视频在线观看| 在线免费观看不下载黄p国产 | 日韩国内少妇激情av| 波多野结衣高清作品| 精品久久久久久成人av| 搡老岳熟女国产| 51国产日韩欧美| 变态另类成人亚洲欧美熟女| 亚洲国产精品成人综合色| 在线播放无遮挡| 97超视频在线观看视频| 色尼玛亚洲综合影院| 亚洲真实伦在线观看| 午夜精品在线福利| 国内精品一区二区在线观看| 久久久久久国产a免费观看| 亚洲天堂国产精品一区在线| ponron亚洲| 天天添夜夜摸| 欧美最黄视频在线播放免费| 国产男靠女视频免费网站| 亚洲国产精品成人综合色| 91九色精品人成在线观看| 五月玫瑰六月丁香| 真实男女啪啪啪动态图| 国产免费av片在线观看野外av| e午夜精品久久久久久久| 国内精品美女久久久久久| 久久天躁狠狠躁夜夜2o2o| 亚洲中文日韩欧美视频| 亚洲国产欧洲综合997久久,| 法律面前人人平等表现在哪些方面| 国产精品久久久久久亚洲av鲁大| 国产淫片久久久久久久久 | 亚洲中文字幕一区二区三区有码在线看| 午夜亚洲福利在线播放| 真人一进一出gif抽搐免费| 搞女人的毛片| 九九热线精品视视频播放| 精品日产1卡2卡| 免费看a级黄色片| 夜夜爽天天搞| 手机成人av网站| av欧美777| 欧美xxxx黑人xx丫x性爽| 亚洲在线观看片| 一进一出好大好爽视频| 日韩有码中文字幕| 999久久久精品免费观看国产| 偷拍熟女少妇极品色| 可以在线观看毛片的网站| 嫩草影视91久久| 久久午夜亚洲精品久久| 久久精品亚洲精品国产色婷小说| 欧美性感艳星| 国产成人啪精品午夜网站| 亚洲精品日韩av片在线观看 | 亚洲人成网站高清观看| 免费看a级黄色片| 黄色片一级片一级黄色片| 成人精品一区二区免费| 综合色av麻豆| 久久精品国产99精品国产亚洲性色| h日本视频在线播放| 国产一区在线观看成人免费| 成年免费大片在线观看| 国产午夜福利久久久久久| 国模一区二区三区四区视频| 欧美乱码精品一区二区三区| 亚洲精品国产精品久久久不卡| av片东京热男人的天堂| 小说图片视频综合网站| 99riav亚洲国产免费| 亚洲人成网站高清观看| 久久中文看片网| 麻豆国产av国片精品| 九九在线视频观看精品| 日本精品一区二区三区蜜桃| 久久久久久久亚洲中文字幕 | 少妇的逼好多水| 一进一出抽搐动态| 999久久久精品免费观看国产| 国产成人影院久久av| 欧美高清成人免费视频www| 亚洲内射少妇av| 性欧美人与动物交配| 日韩精品青青久久久久久| 岛国在线免费视频观看| 少妇熟女aⅴ在线视频| 免费电影在线观看免费观看| 亚洲av成人精品一区久久| 精品国产美女av久久久久小说| 美女 人体艺术 gogo| 亚洲 国产 在线| 欧美在线一区亚洲| 精品福利观看| 一个人免费在线观看的高清视频| 在线观看舔阴道视频| 国产精品av视频在线免费观看| 青草久久国产| 岛国在线观看网站| 国产欧美日韩精品亚洲av| 亚洲人与动物交配视频| 成人av一区二区三区在线看| 啪啪无遮挡十八禁网站| 最新美女视频免费是黄的| 男女午夜视频在线观看| 色在线成人网| 日韩欧美国产一区二区入口| 亚洲精品影视一区二区三区av| 久久久久国产精品人妻aⅴ院| 十八禁网站免费在线| 国产成人av教育| 欧美一区二区精品小视频在线| 狠狠狠狠99中文字幕| 欧美最黄视频在线播放免费| 国产精品三级大全| 亚洲色图av天堂| 日日干狠狠操夜夜爽| 国产成人a区在线观看| 美女黄网站色视频| 两人在一起打扑克的视频| 日本黄色视频三级网站网址| 99精品久久久久人妻精品| 亚洲在线自拍视频| 757午夜福利合集在线观看| 久久久久久久午夜电影| 中文在线观看免费www的网站| 久久精品国产亚洲av香蕉五月| 美女 人体艺术 gogo| 窝窝影院91人妻| 亚洲国产欧洲综合997久久,| 国产精品久久电影中文字幕| 999久久久精品免费观看国产| 波多野结衣高清作品| 午夜福利在线观看吧| 黄色丝袜av网址大全| 国内精品美女久久久久久| 有码 亚洲区| 一级作爱视频免费观看| 午夜久久久久精精品| 欧美绝顶高潮抽搐喷水| 人妻夜夜爽99麻豆av| 中文字幕高清在线视频| 嫩草影视91久久| 十八禁网站免费在线| 18禁裸乳无遮挡免费网站照片| 亚洲人成网站高清观看| 欧美av亚洲av综合av国产av| 国内毛片毛片毛片毛片毛片| 长腿黑丝高跟| 熟女人妻精品中文字幕| 欧美日本视频| 日本 av在线| 亚洲18禁久久av| 99精品久久久久人妻精品| 国产高清三级在线| 久久婷婷人人爽人人干人人爱| 99久国产av精品| 午夜免费成人在线视频| 亚洲av电影不卡..在线观看| 婷婷精品国产亚洲av在线| 亚洲成a人片在线一区二区| 欧美日本亚洲视频在线播放| 黄色片一级片一级黄色片| 亚洲精品一区av在线观看| 国产97色在线日韩免费| 欧美+亚洲+日韩+国产| 国产97色在线日韩免费| 美女大奶头视频| 久99久视频精品免费| 日韩欧美国产一区二区入口| 成年免费大片在线观看| 欧美一区二区精品小视频在线| 很黄的视频免费| 国产精品乱码一区二三区的特点| 51午夜福利影视在线观看| 国产精品国产高清国产av| 日本在线视频免费播放| 日本免费一区二区三区高清不卡| 成年免费大片在线观看| 亚洲av成人不卡在线观看播放网| 成年免费大片在线观看| 老熟妇仑乱视频hdxx| 99热6这里只有精品| 国产精品98久久久久久宅男小说| www.www免费av| 国产精品亚洲美女久久久| 最近最新中文字幕大全免费视频| 免费av观看视频| 亚洲av第一区精品v没综合| 怎么达到女性高潮| 亚洲狠狠婷婷综合久久图片| 国产色爽女视频免费观看| 一级a爱片免费观看的视频| 法律面前人人平等表现在哪些方面| 此物有八面人人有两片| 人妻久久中文字幕网| 麻豆久久精品国产亚洲av| 国产麻豆成人av免费视频| 日韩欧美 国产精品| 午夜福利视频1000在线观看| 在线观看av片永久免费下载| 午夜激情欧美在线| 中文字幕av成人在线电影| 久久久国产成人精品二区| 啦啦啦免费观看视频1| 亚洲av不卡在线观看| 91av网一区二区| 成人18禁在线播放| 一级a爱片免费观看的视频| 级片在线观看| 黄色女人牲交| 可以在线观看毛片的网站| 黄片小视频在线播放| 国产av麻豆久久久久久久| 一进一出好大好爽视频| 国产av麻豆久久久久久久| 97超视频在线观看视频| 久久精品国产综合久久久| 美女高潮喷水抽搐中文字幕| 夜夜躁狠狠躁天天躁| 久久久久久大精品| 男女那种视频在线观看| 国产欧美日韩一区二区三| 成熟少妇高潮喷水视频| 久久人人精品亚洲av| 丰满人妻熟妇乱又伦精品不卡| 五月玫瑰六月丁香| 又黄又粗又硬又大视频| 听说在线观看完整版免费高清| 欧美激情久久久久久爽电影| 青草久久国产| 日本成人三级电影网站| 非洲黑人性xxxx精品又粗又长| 男人舔女人下体高潮全视频| 熟女少妇亚洲综合色aaa.| 88av欧美| 久久久久免费精品人妻一区二区| 天天添夜夜摸| 女同久久另类99精品国产91| 女警被强在线播放| 中文字幕av成人在线电影| 午夜福利成人在线免费观看| www.999成人在线观看| 校园春色视频在线观看| 99热只有精品国产| 亚洲一区二区三区不卡视频| 超碰av人人做人人爽久久 | 国产蜜桃级精品一区二区三区| 亚洲不卡免费看| 女生性感内裤真人,穿戴方法视频| 日本熟妇午夜| 性色av乱码一区二区三区2| 成人特级黄色片久久久久久久| 一区二区三区免费毛片| 男人舔奶头视频| 欧洲精品卡2卡3卡4卡5卡区| 九九久久精品国产亚洲av麻豆| 国产伦精品一区二区三区四那| 69人妻影院| 亚洲欧美日韩东京热| 男插女下体视频免费在线播放| 久久国产精品影院| 老司机午夜十八禁免费视频| 身体一侧抽搐| 国产美女午夜福利| 男女下面进入的视频免费午夜| 欧美一级毛片孕妇| 国产高清视频在线观看网站| 久久精品国产亚洲av涩爱 | 中文字幕高清在线视频| 俄罗斯特黄特色一大片| 一级黄色大片毛片| 亚洲最大成人手机在线| 色在线成人网| 亚洲人成伊人成综合网2020| 久久精品国产清高在天天线| 少妇裸体淫交视频免费看高清| 日本黄色片子视频| 麻豆成人av在线观看| 在线免费观看不下载黄p国产 | 欧美丝袜亚洲另类 | 久久6这里有精品| 色哟哟哟哟哟哟| 欧美性感艳星| 国产探花极品一区二区| 欧美性感艳星| 欧美日韩黄片免| 免费观看的影片在线观看| 人人妻人人澡欧美一区二区| 97超级碰碰碰精品色视频在线观看| 午夜免费观看网址| 淫妇啪啪啪对白视频| 18禁美女被吸乳视频| 亚洲va日本ⅴa欧美va伊人久久| 亚洲欧美日韩高清专用| 在线观看66精品国产| 给我免费播放毛片高清在线观看| netflix在线观看网站| 国产真实伦视频高清在线观看 | 一夜夜www| 尤物成人国产欧美一区二区三区| 伊人久久大香线蕉亚洲五| 色播亚洲综合网| 国产激情偷乱视频一区二区| 五月伊人婷婷丁香| 嫁个100分男人电影在线观看| 熟女少妇亚洲综合色aaa.| 日韩欧美精品v在线| 麻豆成人午夜福利视频| 国产乱人伦免费视频| 舔av片在线| 国产亚洲精品综合一区在线观看| 欧美在线黄色| 老司机午夜十八禁免费视频| 免费av不卡在线播放| av在线天堂中文字幕| 久久久久久久亚洲中文字幕 | 国产精品日韩av在线免费观看| 日日夜夜操网爽| 久久亚洲精品不卡| 久久精品国产99精品国产亚洲性色| 久久精品国产自在天天线| 午夜激情欧美在线| 悠悠久久av| 99热6这里只有精品| 亚洲电影在线观看av| 免费人成视频x8x8入口观看| 国产精品98久久久久久宅男小说| 色综合站精品国产| 国产在视频线在精品| 国产乱人视频| 日本黄色片子视频| 一区二区三区高清视频在线| 在线观看午夜福利视频| 午夜福利在线观看免费完整高清在 | 亚洲精品乱码久久久v下载方式 | 欧美一区二区国产精品久久精品| 午夜激情福利司机影院| 久久香蕉精品热| 嫩草影院入口| 中文亚洲av片在线观看爽| 51午夜福利影视在线观看| 久久久久久久久中文| 18禁在线播放成人免费| 国内精品久久久久精免费| 久久亚洲精品不卡| 色综合欧美亚洲国产小说| 国产主播在线观看一区二区| 亚洲国产色片| 99视频精品全部免费 在线| 女人十人毛片免费观看3o分钟| 老熟妇乱子伦视频在线观看|