• <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)法
    超碰97精品在线观看| 十八禁人妻一区二区| 日本撒尿小便嘘嘘汇集6| 91大片在线观看| 美女国产高潮福利片在线看| 巨乳人妻的诱惑在线观看| 黄色视频,在线免费观看| 亚洲人成电影免费在线| 可以在线观看毛片的网站| 国产精品98久久久久久宅男小说| 在线观看免费视频网站a站| 精品国产乱码久久久久久男人| 女警被强在线播放| 91九色精品人成在线观看| 国产无遮挡羞羞视频在线观看| 国产欧美日韩一区二区三| 中文字幕人妻丝袜制服| 国产麻豆69| 啦啦啦在线免费观看视频4| 热re99久久精品国产66热6| 在线永久观看黄色视频| 18美女黄网站色大片免费观看| 精品日产1卡2卡| 欧美人与性动交α欧美软件| 美女大奶头视频| av在线天堂中文字幕 | 老司机在亚洲福利影院| 日韩欧美一区二区三区在线观看| 黄片大片在线免费观看| 人人妻人人添人人爽欧美一区卜| 天天躁夜夜躁狠狠躁躁| 亚洲狠狠婷婷综合久久图片| 女人精品久久久久毛片| 国产精品国产高清国产av| 美女大奶头视频| 久久人妻熟女aⅴ| 欧美丝袜亚洲另类 | 女人高潮潮喷娇喘18禁视频| 久久久久久久午夜电影 | 午夜日韩欧美国产| 精品久久久精品久久久| 国产免费男女视频| 搡老熟女国产l中国老女人| 成人精品一区二区免费| 亚洲在线自拍视频| 国产在线精品亚洲第一网站| 精品久久久久久久久久免费视频 | 亚洲 欧美 日韩 在线 免费| 色老头精品视频在线观看| 久久 成人 亚洲| 热re99久久国产66热| 欧美日韩亚洲高清精品| 国产亚洲欧美在线一区二区| 日本wwww免费看| 91精品国产国语对白视频| 少妇粗大呻吟视频| 国产精品自产拍在线观看55亚洲| 99久久综合精品五月天人人| 一二三四在线观看免费中文在| 搡老岳熟女国产| 18禁美女被吸乳视频| 久久精品国产99精品国产亚洲性色 | 男女高潮啪啪啪动态图| 亚洲专区国产一区二区| 亚洲成av片中文字幕在线观看| a级毛片在线看网站| 多毛熟女@视频| 九色亚洲精品在线播放| 男女床上黄色一级片免费看| 国产熟女午夜一区二区三区| 黑人操中国人逼视频| 亚洲精品在线观看二区| 少妇裸体淫交视频免费看高清 | www.999成人在线观看| 国产色视频综合| 久久热在线av| 色哟哟哟哟哟哟| 亚洲久久久国产精品| 国产欧美日韩一区二区精品| 搡老熟女国产l中国老女人| 亚洲自拍偷在线| 国产精品自产拍在线观看55亚洲| www.精华液| 亚洲精品一卡2卡三卡4卡5卡| 日韩欧美一区视频在线观看| 亚洲国产欧美一区二区综合| 国产aⅴ精品一区二区三区波| 一区二区三区激情视频| 女性生殖器流出的白浆| 极品人妻少妇av视频| 韩国av一区二区三区四区| 精品久久久久久,| 天堂中文最新版在线下载| 精品午夜福利视频在线观看一区| 老司机在亚洲福利影院| 男人舔女人的私密视频| ponron亚洲| 韩国精品一区二区三区| 在线观看一区二区三区| 日韩精品青青久久久久久| 亚洲 国产 在线| 涩涩av久久男人的天堂| 亚洲伊人色综图| 久热爱精品视频在线9| 黄网站色视频无遮挡免费观看| 丰满饥渴人妻一区二区三| 身体一侧抽搐| 久久人人97超碰香蕉20202| 久久精品影院6| 久久久精品国产亚洲av高清涩受| 欧美精品啪啪一区二区三区| 99精品欧美一区二区三区四区| 真人一进一出gif抽搐免费| 日日摸夜夜添夜夜添小说| 久久人人97超碰香蕉20202| 极品教师在线免费播放| 国产亚洲欧美精品永久| 丰满饥渴人妻一区二区三| 日韩三级视频一区二区三区| 久久草成人影院| 亚洲av五月六月丁香网| 天堂√8在线中文| 欧美日韩福利视频一区二区| 欧美 亚洲 国产 日韩一| 婷婷精品国产亚洲av在线| 国产免费现黄频在线看| 日本wwww免费看| 99精品在免费线老司机午夜| 巨乳人妻的诱惑在线观看| av网站免费在线观看视频| 正在播放国产对白刺激| 久久久久亚洲av毛片大全| 男女之事视频高清在线观看| 亚洲男人的天堂狠狠| 黄网站色视频无遮挡免费观看| 国产91精品成人一区二区三区| 91国产中文字幕| 亚洲在线自拍视频| 国产高清视频在线播放一区| 欧美在线一区亚洲| 免费一级毛片在线播放高清视频 | 99国产综合亚洲精品| 久久久久国产一级毛片高清牌| 水蜜桃什么品种好| 国产午夜精品久久久久久| 国产成人一区二区三区免费视频网站| 免费日韩欧美在线观看| 国产成年人精品一区二区 | 国产又色又爽无遮挡免费看| tocl精华| 一本综合久久免费| 精品乱码久久久久久99久播| 性欧美人与动物交配| 在线观看66精品国产| 啦啦啦 在线观看视频| 日本撒尿小便嘘嘘汇集6| 手机成人av网站| 中文字幕人妻丝袜一区二区| 亚洲欧美激情综合另类| 99国产精品一区二区三区| 国产xxxxx性猛交| 国产免费现黄频在线看| a在线观看视频网站| 国产精品秋霞免费鲁丝片| 美女扒开内裤让男人捅视频| 9191精品国产免费久久| 日韩中文字幕欧美一区二区| 日韩免费av在线播放| 国产精品久久久人人做人人爽| 国产精品爽爽va在线观看网站 | 中文字幕精品免费在线观看视频| 国产精品亚洲一级av第二区| 天堂俺去俺来也www色官网| 久久国产乱子伦精品免费另类| 亚洲激情在线av| 高清黄色对白视频在线免费看| 成人影院久久| 午夜精品国产一区二区电影| 一级,二级,三级黄色视频| 亚洲精品av麻豆狂野| 日韩成人在线观看一区二区三区| 嫁个100分男人电影在线观看| 国产成人精品久久二区二区免费| 亚洲片人在线观看| 国产麻豆69| av免费在线观看网站| 日本黄色日本黄色录像| 免费久久久久久久精品成人欧美视频| 成人av一区二区三区在线看| 麻豆av在线久日| 国产成人免费无遮挡视频| 免费少妇av软件| 国产精品成人在线| 欧美日韩av久久| 免费高清视频大片| www.精华液| 亚洲专区中文字幕在线| 99riav亚洲国产免费| 热re99久久国产66热| 亚洲片人在线观看| 97超级碰碰碰精品色视频在线观看| 侵犯人妻中文字幕一二三四区| 成人免费观看视频高清| 黑人欧美特级aaaaaa片| av欧美777| 午夜福利在线免费观看网站| 中文字幕最新亚洲高清| 美女大奶头视频| 后天国语完整版免费观看| av电影中文网址| 午夜免费观看网址| 中出人妻视频一区二区| 好男人电影高清在线观看| 啪啪无遮挡十八禁网站| netflix在线观看网站| 少妇被粗大的猛进出69影院| www国产在线视频色| 999精品在线视频| 免费女性裸体啪啪无遮挡网站| 久久国产精品人妻蜜桃| 国产又色又爽无遮挡免费看| 桃色一区二区三区在线观看| 侵犯人妻中文字幕一二三四区| 免费高清在线观看日韩| 99精品在免费线老司机午夜| 91麻豆精品激情在线观看国产 | 黄色女人牲交| 大码成人一级视频| 欧美亚洲日本最大视频资源| 黄色视频不卡| 99精品久久久久人妻精品| 91大片在线观看| 在线十欧美十亚洲十日本专区| 国产精品影院久久| 国产精品一区二区三区四区久久 | 一区二区三区国产精品乱码| 国产精品一区二区精品视频观看| 久久久久久人人人人人| 亚洲av美国av| 国产精品 国内视频| 最新在线观看一区二区三区| 可以在线观看毛片的网站| 国产乱人伦免费视频| 国产在线观看jvid| 香蕉国产在线看| 9191精品国产免费久久| 丝袜在线中文字幕| 国产精品一区二区在线不卡| 久久久久久人人人人人| 久久人妻熟女aⅴ| 亚洲在线自拍视频| 午夜免费成人在线视频| 久久99一区二区三区| 多毛熟女@视频| 欧美日韩亚洲高清精品| 精品免费久久久久久久清纯| 久久国产精品人妻蜜桃| 亚洲精品国产区一区二| 国产精品日韩av在线免费观看 | 亚洲七黄色美女视频| 国产成人一区二区三区免费视频网站| 在线观看舔阴道视频| 国产成人系列免费观看| 淫秽高清视频在线观看| 桃色一区二区三区在线观看| 又黄又爽又免费观看的视频| 日韩欧美三级三区| 国产极品粉嫩免费观看在线| 国产精品野战在线观看 | 在线观看免费日韩欧美大片| 成人影院久久| 国产亚洲av高清不卡| 午夜免费观看网址| 欧美日韩视频精品一区| 久久人人爽av亚洲精品天堂| 成人av一区二区三区在线看| 亚洲男人天堂网一区| 久久香蕉激情| 久久 成人 亚洲| 精品人妻1区二区| 国产有黄有色有爽视频| 亚洲avbb在线观看| ponron亚洲| 超色免费av| 色婷婷久久久亚洲欧美| 一夜夜www| 男男h啪啪无遮挡| 男女床上黄色一级片免费看| 丰满迷人的少妇在线观看| 一a级毛片在线观看| av片东京热男人的天堂| 婷婷精品国产亚洲av在线| 亚洲一卡2卡3卡4卡5卡精品中文| 久久久久久免费高清国产稀缺| 俄罗斯特黄特色一大片| 国产精品 欧美亚洲| 黑人巨大精品欧美一区二区mp4| 免费看a级黄色片| av有码第一页| 日韩有码中文字幕| 午夜91福利影院| 国产又色又爽无遮挡免费看| 91麻豆av在线| 巨乳人妻的诱惑在线观看| 精品一区二区三区四区五区乱码| www.自偷自拍.com| 国产熟女午夜一区二区三区| 每晚都被弄得嗷嗷叫到高潮| 久9热在线精品视频| 国产视频一区二区在线看| 狠狠狠狠99中文字幕| 精品国产国语对白av| 久久久水蜜桃国产精品网| 国产aⅴ精品一区二区三区波| 在线观看一区二区三区| 国产精品久久电影中文字幕| 国产精品香港三级国产av潘金莲| 婷婷丁香在线五月| 日韩 欧美 亚洲 中文字幕| 很黄的视频免费| 18禁观看日本| 午夜老司机福利片| 亚洲伊人色综图| 免费在线观看黄色视频的| 国产无遮挡羞羞视频在线观看| 亚洲一区二区三区欧美精品| 黑人巨大精品欧美一区二区mp4| av网站免费在线观看视频| 欧美国产精品va在线观看不卡| 大陆偷拍与自拍| 久久国产亚洲av麻豆专区| 99在线人妻在线中文字幕| 50天的宝宝边吃奶边哭怎么回事| 国产精品二区激情视频| 91老司机精品| 长腿黑丝高跟| 日韩大尺度精品在线看网址 | 十八禁人妻一区二区| av免费在线观看网站| 交换朋友夫妻互换小说| 国产精品1区2区在线观看.| 男女下面插进去视频免费观看| 国产精品久久久久成人av| 叶爱在线成人免费视频播放| 亚洲专区字幕在线| 日本a在线网址| 亚洲专区字幕在线| 日日干狠狠操夜夜爽| 国产精品久久久久成人av| 国内毛片毛片毛片毛片毛片| 高清欧美精品videossex| 叶爱在线成人免费视频播放| 中文字幕精品免费在线观看视频| 亚洲七黄色美女视频| 两性午夜刺激爽爽歪歪视频在线观看 | 亚洲精品美女久久久久99蜜臀| 美女扒开内裤让男人捅视频| 80岁老熟妇乱子伦牲交| 国产单亲对白刺激| 国产成人系列免费观看| 久久久久国产一级毛片高清牌| 新久久久久国产一级毛片| 久久久精品国产亚洲av高清涩受| 最好的美女福利视频网| 首页视频小说图片口味搜索| 国产在线观看jvid| 国产精品久久电影中文字幕| avwww免费| 成人三级做爰电影| 黄色女人牲交| 两个人看的免费小视频| 久久精品影院6| 欧美大码av| 宅男免费午夜| 免费在线观看日本一区| 久久久国产精品麻豆| 美女福利国产在线| 久久久精品欧美日韩精品| 国产激情欧美一区二区| 色综合站精品国产| 制服人妻中文乱码| 村上凉子中文字幕在线| 在线看a的网站| 在线播放国产精品三级| 老汉色∧v一级毛片| 久久精品人人爽人人爽视色| av网站在线播放免费| 女同久久另类99精品国产91| 久久中文字幕一级| 国产aⅴ精品一区二区三区波| 老汉色∧v一级毛片| 18禁美女被吸乳视频| 欧美一级毛片孕妇| 麻豆久久精品国产亚洲av | 韩国av一区二区三区四区| 宅男免费午夜| 少妇的丰满在线观看| 超碰97精品在线观看| 性色av乱码一区二区三区2| www.熟女人妻精品国产| 人人妻,人人澡人人爽秒播| 日本一区二区免费在线视频| 久久欧美精品欧美久久欧美| 久久草成人影院| av在线播放免费不卡| 美女扒开内裤让男人捅视频| 国产一卡二卡三卡精品| 99精品在免费线老司机午夜| 黑丝袜美女国产一区| 一本大道久久a久久精品| 亚洲国产精品一区二区三区在线| 一级黄色大片毛片| 久9热在线精品视频| 男人舔女人的私密视频| 性欧美人与动物交配| 午夜福利欧美成人| 亚洲精品国产色婷婷电影| 一进一出好大好爽视频| xxxhd国产人妻xxx| 成在线人永久免费视频| 亚洲午夜精品一区,二区,三区| 中亚洲国语对白在线视频| 成人特级黄色片久久久久久久| 如日韩欧美国产精品一区二区三区| 老熟妇仑乱视频hdxx| 精品久久久久久电影网| 日韩 欧美 亚洲 中文字幕| www.熟女人妻精品国产| 日本黄色日本黄色录像| 香蕉久久夜色| 久久久久久久久久久久大奶| 妹子高潮喷水视频| 精品无人区乱码1区二区| 国产野战对白在线观看| 多毛熟女@视频| 岛国视频午夜一区免费看| 一级毛片精品| 久久人妻福利社区极品人妻图片| 国产精品98久久久久久宅男小说| 男人操女人黄网站| 欧美日本亚洲视频在线播放| 夫妻午夜视频| 亚洲精品美女久久av网站| 9191精品国产免费久久| 一个人免费在线观看的高清视频| 国产亚洲精品一区二区www| 日日干狠狠操夜夜爽| 又紧又爽又黄一区二区| 在线观看免费午夜福利视频| 天天躁狠狠躁夜夜躁狠狠躁| 免费av中文字幕在线| 亚洲成av片中文字幕在线观看| 一区二区三区国产精品乱码| 黄频高清免费视频| 午夜两性在线视频| 99久久国产精品久久久| 天天影视国产精品| 欧美久久黑人一区二区| 日韩视频一区二区在线观看| 日韩精品中文字幕看吧| 久久性视频一级片| 亚洲精品一区av在线观看| 校园春色视频在线观看| 国产成年人精品一区二区 | 69av精品久久久久久| 亚洲熟女毛片儿| 欧美日韩福利视频一区二区| 丁香欧美五月| 黑人欧美特级aaaaaa片| 国产精品自产拍在线观看55亚洲| 免费日韩欧美在线观看| 免费观看人在逋| 国产亚洲精品综合一区在线观看 | 国产精品 国内视频| 美女福利国产在线| 免费女性裸体啪啪无遮挡网站| 一二三四在线观看免费中文在| 真人一进一出gif抽搐免费| 亚洲少妇的诱惑av| 欧洲精品卡2卡3卡4卡5卡区| 国产欧美日韩精品亚洲av| 精品久久久精品久久久| 久久精品成人免费网站| 大型黄色视频在线免费观看| 欧美日韩亚洲国产一区二区在线观看| 亚洲精品美女久久久久99蜜臀| 亚洲欧美一区二区三区久久| 色精品久久人妻99蜜桃| 男人舔女人的私密视频| 熟女少妇亚洲综合色aaa.| 日本精品一区二区三区蜜桃| 亚洲一区二区三区不卡视频| 丝袜人妻中文字幕| 日韩欧美国产一区二区入口| 熟女少妇亚洲综合色aaa.| 91精品三级在线观看| 国产伦一二天堂av在线观看| 国产成人影院久久av| 日韩精品青青久久久久久| 美女扒开内裤让男人捅视频| 啦啦啦 在线观看视频| 不卡av一区二区三区| 视频在线观看一区二区三区| 亚洲 国产 在线| 精品久久久久久成人av| 亚洲熟妇熟女久久| 欧美精品啪啪一区二区三区| 真人一进一出gif抽搐免费| 精品久久久久久,| 久久国产精品人妻蜜桃| 久久中文字幕一级| 一区二区三区激情视频| av片东京热男人的天堂| 宅男免费午夜| 好看av亚洲va欧美ⅴa在| 午夜福利在线观看吧| 精品电影一区二区在线| 欧美另类亚洲清纯唯美| 啦啦啦 在线观看视频| www日本在线高清视频| 制服诱惑二区| 国产又爽黄色视频| 一级毛片女人18水好多| 国产成人精品久久二区二区91| 老汉色∧v一级毛片| 精品久久蜜臀av无| 日本黄色日本黄色录像| 女性被躁到高潮视频| 一区在线观看完整版| 国产又爽黄色视频| 亚洲色图综合在线观看| cao死你这个sao货| 男人舔女人的私密视频| 成人18禁在线播放| 一个人免费在线观看的高清视频| 少妇 在线观看| 老司机午夜十八禁免费视频| 亚洲精品一二三| 亚洲欧美一区二区三区久久| 午夜免费激情av| av中文乱码字幕在线| 最新美女视频免费是黄的| 99精国产麻豆久久婷婷| 一区二区三区激情视频| 可以免费在线观看a视频的电影网站| 亚洲男人天堂网一区| 免费在线观看视频国产中文字幕亚洲| 男人操女人黄网站| 国产精品香港三级国产av潘金莲| 免费一级毛片在线播放高清视频 | www.熟女人妻精品国产| 黑丝袜美女国产一区| 黄网站色视频无遮挡免费观看| 亚洲一区高清亚洲精品| 亚洲精品美女久久久久99蜜臀| 热99国产精品久久久久久7| 久久香蕉激情| 精品人妻在线不人妻| 亚洲av成人一区二区三| 嫁个100分男人电影在线观看| 久久久国产欧美日韩av| 夜夜躁狠狠躁天天躁| 亚洲一区二区三区不卡视频| 久久九九热精品免费| 在线观看一区二区三区激情| 亚洲色图综合在线观看| 国产精品 国内视频| 国产av在哪里看| 91大片在线观看| 男女之事视频高清在线观看| 亚洲一区二区三区色噜噜 | 91成人精品电影| 我的亚洲天堂| 久久久久久久精品吃奶| 97碰自拍视频| 欧美久久黑人一区二区| 国产成人精品久久二区二区91| 91字幕亚洲| 级片在线观看| 精品免费久久久久久久清纯| 纯流量卡能插随身wifi吗| 高清在线国产一区| 少妇裸体淫交视频免费看高清 | 99riav亚洲国产免费| 久久国产乱子伦精品免费另类| 免费观看精品视频网站| 中文字幕人妻丝袜制服| 777久久人妻少妇嫩草av网站| 后天国语完整版免费观看| 黑人巨大精品欧美一区二区蜜桃| 欧美丝袜亚洲另类 | 一边摸一边抽搐一进一出视频| 国产极品粉嫩免费观看在线| 亚洲激情在线av| 欧美性长视频在线观看| 别揉我奶头~嗯~啊~动态视频| 日韩欧美一区二区三区在线观看| 亚洲一区二区三区不卡视频| 一级作爱视频免费观看| 他把我摸到了高潮在线观看| 欧美 亚洲 国产 日韩一| 女警被强在线播放| 精品乱码久久久久久99久播| 色老头精品视频在线观看| x7x7x7水蜜桃| 欧美黄色片欧美黄色片| 一边摸一边抽搐一进一出视频| aaaaa片日本免费| 长腿黑丝高跟| 亚洲欧美日韩另类电影网站| 可以免费在线观看a视频的电影网站| 国产精品九九99| 亚洲欧美激情在线| 精品日产1卡2卡| 久久久精品欧美日韩精品|