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

        ?

        高斯脈沖調(diào)制Chirp信號(hào)分?jǐn)?shù)階Fourier變換的解析表示

        2010-05-13 09:17:24瑜,陳
        現(xiàn)代電子技術(shù) 2009年19期
        關(guān)鍵詞:參數(shù)估計(jì)調(diào)頻高斯

        范 瑜,陳 蘭

        摘 要:分?jǐn)?shù)階Fourier變換是Fourier變換的一種廣義形式,揭示了信號(hào)從時(shí)域到頻域變化過程中所呈現(xiàn)的特征,運(yùn)用FrFT分析Chirp信號(hào)具有極佳的能量積聚效果。介紹分?jǐn)?shù)階Fourier變換的基本定義,給出高斯脈沖Chirp信號(hào)的FrFT以及最優(yōu)化高斯脈沖Chirp信號(hào)分析角度,對(duì)于該類信號(hào)的低噪聲檢測與參數(shù)估計(jì)具有良好的應(yīng)用前景。

        關(guān)鍵詞:Fourier變換;分?jǐn)?shù)階Fourier變換;Chirp信號(hào);高斯脈沖調(diào)制Chirp信號(hào)

        中圖分類號(hào):TN911文獻(xiàn)標(biāo)識(shí)碼:A

        文章編號(hào):1004-373X(2009)17-069-04

        Fractional Fourier Transform of Gauss Pulse Modulated Chirp Signal and

        Its Analytical Representation

        FAN Yu1,CHEN Lan2

        (1.Changshu Institute of Technology,Changshu,215500,China;2.Software Institute,Southeast University,Suzhou,215123,China)

        Abstract:Fractional Fourier transform is a version of generalized Fourier transform,which reveals the features of signal from time and frequency simultaneously.The analysis of Chirp signal using FrFT can bring excellent energy accumulation effect.The definition of fractional Fourier transform is presented.Fractional Fourier transform of Gauss pulse Chirp signal is analyzed and the design of optimal angle is also studied in Gauss pulse modulated Chirp signal analysis,which have a good future in signal detection and parameter estimation when in low signal noise ratio situation.

        Keywords:Fourier transform;fractional Fourier transform;Chirp signal;Gauss pulse modulated Chirp signal

        0 引 言

        Fourier變換(Fourier Transform,FT)作為最主要的信號(hào)分析工具主要用于處理頻率不隨時(shí)間變化的平穩(wěn)信號(hào),在時(shí)頻平面時(shí)間軸與頻率軸相互垂直,即Fourier變換是從時(shí)間域旋轉(zhuǎn)π/2到頻率域。但是Fourier變換通常無法表述信號(hào)的時(shí)頻域性質(zhì),不能表示某種頻率分量發(fā)生在哪個(gè)時(shí)間,而這種性質(zhì)恰恰是非平穩(wěn)信號(hào)最關(guān)鍵的性質(zhì),這對(duì)非平穩(wěn)信號(hào)十分重要。分?jǐn)?shù)階Fourier變換(Fractional Fourier Transform,FrFT)是Fourier變換的一種推廣形式,揭示了信號(hào)從時(shí)間域到頻率域變化過程中所呈現(xiàn)的特征,即從時(shí)間域和頻率域同時(shí)表示信號(hào)旋轉(zhuǎn)π/2分?jǐn)?shù)倍時(shí)的特征,從而克服了傳統(tǒng)Fourier變換不能反映非平穩(wěn)信號(hào)的統(tǒng)計(jì)量隨時(shí)間變化的缺陷[1-4]。

        分?jǐn)?shù)階Fourier變換是一種線性變換,與經(jīng)典的Fourier變換有著天然的聯(lián)系,又提供了Fourier變換所不具備的某些特點(diǎn),而且它與小波變換、Wigner-Ville分布(WVD)都有密切的關(guān)系[3],因此近年來分?jǐn)?shù)階Fourier變換受到了研究者的廣泛關(guān)注,相應(yīng)的離散變換算法也相繼提出[4]。此外,有關(guān)分?jǐn)?shù)階Fourier變換應(yīng)用的研究也越來越多,如[5-10]等文獻(xiàn)就介紹了分?jǐn)?shù)階Fourier變換在信號(hào)處理和通信系統(tǒng)方面的應(yīng)用。

        線性調(diào)頻(Chirp)信號(hào)是一種典型的非平穩(wěn)信號(hào),它的瞬時(shí)頻率隨時(shí)間呈線性變化,常見于雷達(dá)、聲納和移動(dòng)通信等系統(tǒng)中。對(duì)Chirp信號(hào)的研究,是非平穩(wěn)隨機(jī)信號(hào)處理理論及方法的基礎(chǔ)。研究證明,分?jǐn)?shù)階Fourier變換就是一種很適合處理Chirp信號(hào)的變換。近年來,將分?jǐn)?shù)階Fourier變換用于線性調(diào)頻信號(hào)(包括時(shí)不變、時(shí)變幅度線性調(diào)頻信號(hào))的檢測和參數(shù)估計(jì)引起了越來越多的關(guān)注[11,12]。本文介紹了FrFT的定義、性質(zhì)和簡單應(yīng)用,分析了高斯脈沖Chirp信號(hào)的FrFT以及最優(yōu)化分析角度的選取問題。

        1 分?jǐn)?shù)階Fourier變換的基礎(chǔ)

        分析和處理平穩(wěn)信號(hào)最常用和最主要的方法是Fourier變換(FT)。Fourier變換建立了信號(hào)整個(gè)時(shí)域與整個(gè)頻域的對(duì)應(yīng)關(guān)系,其中:

        X(ω)=1/2π∫∞-∞x(t)e-jωtdt(1)

        x(t)=1/2π∫∞-∞X(ω)ejωtdt(2)

        時(shí)域和頻域構(gòu)成了分析一個(gè)信號(hào)的兩種表達(dá)方式,Fourier變換在整體時(shí)域上將信號(hào)分解為不同的頻率分量,但它沒有將時(shí)域和頻域組合成一個(gè)域,不能提供時(shí)間和頻率的聯(lián)合信息。譜X(ω)只是顯示任一頻率ω包含在信號(hào)x(t)內(nèi)的總強(qiáng)度,無法表述信號(hào)的時(shí)頻域性質(zhì),不能表示某種頻率分量發(fā)生在哪個(gè)時(shí)間,而這種性質(zhì)恰恰是非平穩(wěn)信號(hào)中的最根本和最關(guān)鍵的性質(zhì)。如果x(t)是由幾個(gè)非平穩(wěn)分量組成的,那么時(shí)間上的任何變化都會(huì)改變X(ω)。此時(shí),傳統(tǒng)的Fourier變換就不能滿足信號(hào)分析的要求了。

        時(shí)頻分析的基本思想就是設(shè)計(jì)時(shí)間和頻率的聯(lián)合函數(shù),用它描述信號(hào)在不同時(shí)間和頻率上的能量密度或強(qiáng)度,從而克服傳統(tǒng)Fourier變換不能反映非平穩(wěn)信號(hào)的統(tǒng)計(jì)量隨時(shí)間變化的缺陷。分?jǐn)?shù)階Fourier變換(FrFT)是借用時(shí)頻面的概念,以時(shí)間和頻率分別為橫軸和縱軸,旋轉(zhuǎn)一定的角度進(jìn)行的一種線性變換。

        傳統(tǒng)的Fourier變換X(ω)就是x(t)旋轉(zhuǎn)π/2,即x(t)由時(shí)間軸t變到頻率軸ω的表示形式。令:

        α=pπ/2(3)

        并且定義線性算子:

        Rα=Rpπ/2(4)

        記作Fp。F2相當(dāng)于t軸連續(xù)兩次逆時(shí)針旋轉(zhuǎn)π/2,得到x(-t);F3相當(dāng)于t軸連續(xù)三次逆時(shí)針旋轉(zhuǎn)π/2,得到指向-ω軸的函數(shù);F4表示t軸連續(xù)四次逆時(shí)針旋轉(zhuǎn)π/2,得到原函數(shù)。因此線性算子Rα有以下數(shù)學(xué)性質(zhì):

        (1) 零旋轉(zhuǎn):

        R0=x(t)(5)

        (2) 與Fourier變換等價(jià):

        Rπ/2=F[x(t)](6)

        (3) 旋轉(zhuǎn)相加性:

        RαRβ=Rα+β(7)

        (4) 2π旋轉(zhuǎn)(恒等變換):

        R2π=x(t)(8)

        如果角度α以π/2的非整數(shù)倍進(jìn)行旋轉(zhuǎn),則得到函數(shù)x(t)的廣義Fourier變換,記作:

        {Rαx}(u)={Fpx}(u)=Xp(u)(9)

        這就是x(t)的分?jǐn)?shù)階Fourier變換:

        Xp(u)={Fpx}(u)=∫∞-∞x(t)Kp(t,u)dt(10)

        式中:變換核Kp(t,u)定義為:

        Kp(t,u)=1-jcot α2πejt2+u22cot α-jutcsc α, α≠nπ

        δ(t-u), α=2nπ

        δ(t+u),α=(2n+1)π(11)

        式中:n為整數(shù)。

        從以上的定義可以看出,分?jǐn)?shù)階Fourier變換是經(jīng)典Fourier變換的廣義形式,它包含了信號(hào)的時(shí)間域和頻率域表示。旋轉(zhuǎn)角度為π/2時(shí),即階數(shù)為1的分?jǐn)?shù)階Fourier變換就是傳統(tǒng)的Fourier變換;不旋轉(zhuǎn)或旋轉(zhuǎn)的角度為2π的整數(shù)倍則為信號(hào)本身;當(dāng)旋轉(zhuǎn)角度不在以上兩個(gè)位置即p為分?jǐn)?shù)時(shí),它同時(shí)從時(shí)間域和頻率域給出了信號(hào)的特征。在討論分?jǐn)?shù)階Fourier變換時(shí),由于角度以2π為模,所以只需要考慮0≤p≤2的旋轉(zhuǎn)階數(shù)。α角度的分?jǐn)?shù)階Fourier反變換對(duì)應(yīng)著-α角度的分?jǐn)?shù)階Fourier變換,即:

        x(t)=∫∞-∞Xp(u)K-p(t,u)du(12)

        Fourier變換在線性系統(tǒng)分析、光學(xué)系統(tǒng)、信息處理系統(tǒng)等方面起著核心作用,并應(yīng)用于眾多的工程技術(shù)領(lǐng)域。作為廣義Fourier變換,FrFT具有比Fourier變換更普遍的特性和更廣的應(yīng)用場合,尤其是普通Fourier變換技術(shù)不能解決問題的場合,FrFT更顯其優(yōu)越性。在FrFT研究的啟示下,許多學(xué)者推廣了分?jǐn)?shù)階的概念,得到分?jǐn)?shù)階卷積和分?jǐn)?shù)階相關(guān),還把分?jǐn)?shù)階的概念應(yīng)用到Hadamard變換、Hartley變換等,得到相應(yīng)的分?jǐn)?shù)階變換。而且,分?jǐn)?shù)階Fourier變換已應(yīng)用到圖象處理的優(yōu)化圖象恢復(fù)方面。此外,將分?jǐn)?shù)階Fourier變換算子的分?jǐn)?shù)冪推廣至復(fù)數(shù)冪也成為目前研究的一個(gè)熱點(diǎn)。

        在20世紀(jì)90年代中期,分?jǐn)?shù)階Fourier變換被引入了信號(hào)處理領(lǐng)域。在信號(hào)處理中,FrFT有很多應(yīng)用,其中兩個(gè)典型的應(yīng)用是信號(hào)濾波和信號(hào)分離[5]。實(shí)驗(yàn)已經(jīng)證實(shí)分?jǐn)?shù)階Fourier變換濾波的效果明顯優(yōu)于Fourier變換;特別是對(duì)于線性調(diào)頻Chirp信號(hào),分?jǐn)?shù)階Fourier變換能夠獲得最佳的能量積聚效果,是分?jǐn)?shù)階Fourier變換最合適的應(yīng)用領(lǐng)域之一。

        2 高斯脈沖Chirp信號(hào)的分?jǐn)?shù)階Fourier變換

        線性調(diào)頻Chirp信號(hào)是一種特殊的非平穩(wěn)信號(hào),它的瞬時(shí)頻率隨時(shí)間呈線性變化,廣泛地出現(xiàn)在通信、雷達(dá)、聲納和地震勘探等系統(tǒng)中。在工程實(shí)踐中,高斯調(diào)制Chirp信號(hào)有著廣泛的應(yīng)用,而FrFT對(duì)于Chirp信號(hào)良好的檢測和分析效果自然讓我們想到能否將此分析用具引入到高斯調(diào)制的Chirp信號(hào)分析中。本文對(duì)此進(jìn)行了深入研究,給出了相關(guān)解析結(jié)果。

        高斯信號(hào)的標(biāo)準(zhǔn)形式為:

        f(t)=12πδ2e-(t-t0)22δ2(13)

        這是一個(gè)服從正態(tài)分布(t0,δ)的函數(shù),它的Fourier變換和分?jǐn)?shù)階Fourier變換都具有高斯信號(hào)的形式[12]。利用已有的這些性質(zhì),可以研究高斯脈沖Chirp信號(hào)的分?jǐn)?shù)階Fourier變換。

        令高斯脈沖Chirp信號(hào)為:

        x(t)=e-(t-t0)22δ2ej(at2+bt+c)(14)

        將高斯脈沖Chirp信號(hào)寫成幅度和相位的函數(shù),得到:

        x(t)=x(t)ejθ(ω)(15)

        其中:

        x(t)=abs(e-(t-t0)22δ2ej(at2+bt+c))=e-(t-t0)22δ2

        θ(ω)=at2+bt+c(16)

        式(14)所表示的高斯脈沖Chirp信號(hào)的分?jǐn)?shù)階Fourier變換為:

        Xp(u)=1-jcot α2πeju22cot α∫e-(t-t0)22δ2ej(at2+bt+c)ejt22cot αe-jutcsc αdt

        =1-jcot α2πejceju22cot αe-t202δ2∫e-[(12δ2-ja-jcot α2)t2-(t0δ2+jb-jucsc α)t]dt

        =P∫e-(M2t2-2MNt)dt(17)

        其中:

        P=1-jcot α2πejceju22cot αe-t202δ2

        M=12δ2-ja-jcot α2

        N=t0δ2+jb-jucsc α12δ2-ja-jcot α2(18)

        所以式(17)化為:

        Xp(u)=P∫e-(M2t2-2MNt+N2-N2)dt

        =PeN2∫e-(Mt-N)2dt(19)

        令Mt-N=Q,則t=Q+NM,dt=dQM,因此得到:

        Xp(u)=PeN2M∫e-Q2dQ(20)

        然后取積分限為(-∞,∞),由于∫∞-∞e-Q2dQ =π,所以可得:

        Xp(u)=πPeN2M(21)

        再將P,M,N的值代入式(21)得:

        Xp(u)=1-jcot α2ejceju22cot αe-t202δ2?

        112δ2-ja-jcot α2et0δ2+jb-jucsc α22δ2-j4a-j2cot α(22)

        將(t0/δ2+jb-jucsc α)22/δ2-j4a-j2cot α的分子分母同時(shí)乘以2/δ2+j4a+j2cot α,則式(22)化為:

        Xp(u)=1-jcot α2ejceju22cot αe-t20 2δ2112δ2-ja-jcot α2?

        et0δ2+jb-jucsc α22δ2+j4a+j2cot α4δ4+(4a + 2cot α)2(23)

        若t0=0,可以得到:

        Xp(u)=1-jcot α2112δ2-ja-jcot α2?

        ejceju22cot αe-j(4a+2cot α)(ucsc α-b)24δ4+(4a+2cot α)2e-2δ2(ucsc α-b)24δ4+(4a+2cot α)2(24)

        將Chirp信號(hào)的分?jǐn)?shù)階Fourier變換寫成幅度和相位的函數(shù),得到:

        Xp(u)=Xp(u)ejρ(ω)(25)

        因?yàn)閑jA=1,在計(jì)算Xp(u)的幅值時(shí)可以忽略所有模為1的部分,這樣就能得到:

        Xp(u)=

        abs1-jcotα2112δ2-ja-jcot α2e-(u-bsin α)22sin2α[1δ2+δ2(2a+cot α)2](26)

        令: σ2=sin2α[1δ2+δ2(2a+cot α)2](27)

        得:

        Xp(u)=abs1-jcot α2

        112δ2-ja-jcot α2e-(u-bsin α)22σ2(28)

        由式(24)和式(25)得:

        ejρ(ω)=ejc+u22cot α-(4a+2cot α)(ucsc α-b)24δ4+(4a+2cot α)2(29)

        因此Xp(u)的相位:

        ρ(ω)=c+u22cot α-(4a+2cot α)(ucsc α-b)24δ4+(4a+2cot α)2(30)

        由式(28)可得Xp(u)具有高斯函數(shù)的形式,服從正態(tài)分布(bsin α,σ2)。當(dāng)方差σ2最小時(shí),高斯脈沖Chirp信號(hào)的分?jǐn)?shù)階Fourier變換Xp(u)的能量最大,因此這時(shí)的α就是最優(yōu)化角度αopt。從式(27)可知,σ2最小時(shí),有2a+cot αopt=0,即:

        αopt=-arctan12a(31)

        同樣,由式(3)可知,此時(shí)Xp(u)的最優(yōu)化階數(shù)popt為:

        popt=-2πarctan12a(32)

        由式(31)和式(32)對(duì)比文獻(xiàn)[12]可得,Chirp信號(hào)和高斯脈沖調(diào)制Chirp信號(hào)的αopt,popt是一樣的,這說明加上高斯窗的Chirp信號(hào)(即高斯脈沖Chirp信號(hào))的最優(yōu)化參數(shù)并沒有改變。這就為高斯Chirp信號(hào)的FrFT最優(yōu)化分析奠定了堅(jiān)實(shí)的理論基礎(chǔ),可以充分利用FrFT對(duì)于LFM信號(hào)最優(yōu)的能量積聚性質(zhì),將其應(yīng)用于高斯Chirp信號(hào)的分析中,在低信噪比的環(huán)境中這一點(diǎn)更具有實(shí)際意義。

        3 結(jié) 語

        分?jǐn)?shù)階Fourier變換是Fourier變換的廣義形式。作為一種新的時(shí)頻分析工具,分?jǐn)?shù)階Fourier變換既與經(jīng)典的Fourier變換有著天然的聯(lián)系,又提供了Fourier變換所不具備的某些特點(diǎn),而且它與小波變換、Wigner-Ville分布(WVD)都有密切的關(guān)系[12,13],是一種很適合處理線性調(diào)頻Chirp信號(hào)等非平穩(wěn)信號(hào)的變換。

        在Chirp信號(hào)的檢測和參數(shù)估計(jì)方面,分?jǐn)?shù)階Fourier變換具有重要的應(yīng)用價(jià)值。尤其是它所具有的線性特性,可以方便地用于多分量線性調(diào)頻信號(hào)的檢測和參數(shù)估計(jì)而不受交叉項(xiàng)的干擾,大量被用于恒定幅度的線性調(diào)頻信號(hào),并取得了良好的效果。事實(shí)上,時(shí)變幅度線性調(diào)頻信號(hào)在工程實(shí)際中有著更為廣泛的應(yīng)用背景,在工程實(shí)踐中,會(huì)遇到很多經(jīng)過幅度調(diào)制的線性調(diào)頻信號(hào),其中的高斯窗調(diào)制就是典型的一種。

        從分析結(jié)果看,經(jīng)過高斯窗調(diào)制的Chirp信號(hào)保持了原本的特性,在FrFT這個(gè)分析工具下能夠獲得優(yōu)異的性能。

        參考文獻(xiàn)

        [1]于鳳芹,姚旭輝,曹家麟.分?jǐn)?shù)階傅立葉變換的若干問題[J].江南大學(xué)學(xué)報(bào):自然科學(xué)版,2002,1(4):349-353.

        [2]Ahmed I Zayed.On the Relationship between the Fourier and Fractional Fourier Transforms[J].IEEE Signal Processing Letters,1996,12(3):310-311.

        [3]Jian Jiun Ding,Soo Chang Pei.Fractional Fourier Transforms and Wigner Distribution Functions for Stationary and Non-Stationary Random Process[A].ICASSP 2006[C].2006,(3):428-431.

        [4]平先軍,陶然,周思永.一種新的分?jǐn)?shù)階傅立葉變換快速算法[J].電子學(xué)報(bào),2001,3(3):406-408.

        [5]姚紅玉.分?jǐn)?shù)傅立葉變換在信號(hào)處理中的應(yīng)用研究[J].數(shù)值計(jì)算與計(jì)算機(jī)應(yīng)用,2007,28(3):198-214.

        [6]陳恩慶,陶然,張衛(wèi)強(qiáng).一種基于分?jǐn)?shù)階傅立葉變換的時(shí)變信道參數(shù)估計(jì)方法[J].電子學(xué)報(bào),2005(12):2 101-2 104.

        [7]薛文虎,張明敏,唐勁松,等.Chirp信號(hào)參數(shù)估計(jì)算法性能比較[J].海軍工程大學(xué)學(xué)報(bào),2007,19(2):1-5.

        [8]Wei Yu,Huang Chang,Sun De bao,et al.The Detection of the Linear Frequency modulated Signal by Fractional Fourier Transform[A].Proceedings of the Second International Conference on Information Technology for Application[C].2004:322-325.

        [9]楊秀麗,馮小平.基于分?jǐn)?shù)階傅里葉變換的Chirp信號(hào)參數(shù)估計(jì)及恢復(fù)[J].空間電子技術(shù),2007(3):78-81.

        [10]牛虹,齊林,宋家友.基于FrFT的時(shí)變幅度Chirp信號(hào)的參數(shù)估計(jì)[J].現(xiàn)代雷達(dá),2007,29(11):37-43.

        [11]Olcay Akay G,Faye Boudreaux-Bartels.Fractional Autocorrelation and Its Application To Detection and Estimation of Linear FM Signals[A].Time-Frequency and Time-Scale Analysis,1998.Proceedings of the IEEE-SP International Symposium[C].1998:213-216.

        [12]Capus C,Brown K.Fractional Fourier Transform of the Gaussian and Fractional Domain Signal Support[J].IEEE Proc.Vis.Image Signal Process,2003,150(2):99-106.

        [13]Luis Almeida B.The Fractional Fourier Transform and Time-Frequency Representations[J].IEEE Trans.on Signal Processing,1994,42(11):3 084-3 091.

        猜你喜歡
        參數(shù)估計(jì)調(diào)頻高斯
        小高斯的大發(fā)現(xiàn)
        基于新型DFrFT的LFM信號(hào)參數(shù)估計(jì)算法
        考慮頻率二次跌落抑制的風(fēng)火聯(lián)合一次調(diào)頻控制
        能源工程(2021年5期)2021-11-20 05:50:42
        天才數(shù)學(xué)家——高斯
        Logistic回歸模型的幾乎無偏兩參數(shù)估計(jì)
        基于向前方程的平穩(wěn)分布參數(shù)估計(jì)
        基于競爭失效數(shù)據(jù)的Lindley分布參數(shù)估計(jì)
        調(diào)頻發(fā)射機(jī)技術(shù)改造
        調(diào)頻激勵(lì)器干擾的排除方法
        有限域上高斯正規(guī)基的一個(gè)注記
        少妇夜夜春夜夜爽试看视频| 中国免费av网| 亚洲av伊人久久综合性色| 日本免费三片在线视频| 亚洲av免费不卡在线观看| 肉色欧美久久久久久久免费看| 亚洲av成人无码网站大全| 国内精品久久久久久久久齐齐| 爱v天堂在线观看| 日本av一区二区三区四区| 免费亚洲老熟熟女熟女熟女| 久久久久无码国产精品一区| 欧美疯狂做受xxxx高潮小说| 人妻少妇一区二区三区| 精品国产夫妻自拍av| 日本一区二区视频在线| 中文字幕一精品亚洲无线一区| 久久人人玩人妻潮喷内射人人| 毛片av在线播放亚洲av网站| 福利视频自拍偷拍视频| 经典三级免费看片天堂| 日本少妇浓毛bbwbbwbbw| 五月婷婷俺也去开心| 国产v精品成人免费视频400条| 蜜桃成人精品一区二区三区| 东北女人一级内射黄片| 国产啪亚洲国产精品无码 | 手机在线中文字幕av| 少妇人妻中文久久综合| 丰满人妻一区二区三区视频53| 少妇人妻偷人精品无码视频| 日本人妻av在线观看| 久久精品久99精品免费| 青青青爽在线视频观看| 国产精品成人av在线观看| 亚洲精品无人区一区二区三区| 老女人下面毛茸茸的视频| 一二区成人影院电影网| 精品国产精品久久一区免费式| 国产在线白浆一区二区三区在线| 国产极品大奶在线视频|