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

隨機共振現(xiàn)象與微弱信號接收

2015-09-23 08:31:26石盛超李廣俠李志強朱文明張衛(wèi)同
數(shù)字通信世界 2015年3期
關(guān)鍵詞:信號系統(tǒng)

石盛超,李廣俠,李志強,朱文明,張衛(wèi)同

(解放軍理工大學通信工程學院天基信息教研中心,南京 210007)

隨機共振現(xiàn)象與微弱信號接收

石盛超,李廣俠,李志強,朱文明,張衛(wèi)同

(解放軍理工大學通信工程學院天基信息教研中心,南京 210007)

本文從隨機共振(Stochastic Resonance)可增強信號信噪比的現(xiàn)象出發(fā),介紹了隨機共振非線性雙穩(wěn)勢阱模型和參數(shù)調(diào)節(jié)隨機共振理論。在此基礎(chǔ)上,分析了隨機共振現(xiàn)象在微弱信號接收中兩種可能的應(yīng)用途徑,通過仿真驗證了應(yīng)用的可行性。最后,分析提出了將隨機共振理論應(yīng)用于弱信號接收需要突破的三大關(guān)鍵技術(shù),即:大參數(shù)信號變尺度與還原技術(shù)、參數(shù)自適應(yīng)調(diào)節(jié)技術(shù)和調(diào)制信號隨機共振技術(shù),還提出了初步解決的思路。為拓寬隨機共振現(xiàn)象的應(yīng)用領(lǐng)域和探究解決低信噪比條件下的信號接收提供新的方法和途徑。

隨機共振;微弱信號接收;低載噪比;參數(shù)調(diào)節(jié)

1 引言

意大利學者R.Benzi在研究古氣象冰川演化問題時發(fā)現(xiàn):在過去的70萬年中,地球的冰川期和暖氣候期大概以10萬年為周期交替出現(xiàn)。對這一時期地球氣候的研究結(jié)果表明,地球繞太陽轉(zhuǎn)動偏心率的變化周期大約為10萬年,這意味著太陽對地球施加了周期信號。然而,這一周期信號本身很小,不足以使地球氣候產(chǎn)生從冰川期到暖氣候期的大幅度變化。為此,1981年R.Benzi[1]-[3]等人首次提出了隨機共振的概念,對這種現(xiàn)象進行了很好的解釋,他們認為這一周期信號與地球本身能夠產(chǎn)生冷態(tài)和暖態(tài)的非線性條件,以及地球所受到的隨機力(如太陽的各種無規(guī)則變化)三者之間達到了協(xié)同,通過隨機共振引起了地球古氣象的大幅度變動。

隨機共振理論一經(jīng)提出,不僅得到了物理、化學、生物醫(yī)學、控制、通信和傳感等自然科學界的關(guān)注,甚至被用于解決經(jīng)濟和社會科學界面臨的一些關(guān)鍵問題,成為近三十年來持續(xù)的研究熱點。

隨機共振理論[4]-[7]指出,當淹沒在強噪聲背景中的微弱信號通過一個非線性系統(tǒng)時,當非線性系統(tǒng)、信號和噪聲達到匹配時,噪聲的能量會增強信號,使得處理后的信號信噪比增加。這為豐富信號檢測方法提供了新的思路。與傳統(tǒng)經(jīng)典的消除或者抑制噪聲的處理方法不同,一般的隨機共振用于弱信號檢測的方法[8]-[13]是:設(shè)計基于隨機共振的非線性檢測器,通過調(diào)整系統(tǒng)參數(shù)或改變進入檢測器的噪聲,使得非線性系統(tǒng)、信號和噪聲達到隨機共振狀態(tài),以隨機共振輸出信號為依據(jù),可以實現(xiàn)微弱信號檢測的目的。

2 隨機共振理論

隨機共振現(xiàn)象發(fā)生的三個基本要素是非線性系統(tǒng)、輸入信號和噪聲??紤]最簡單的非線性雙穩(wěn)態(tài)系統(tǒng),它可以由朗之萬(langevin)微分方程[8]-[13]描述,如式(1)所示:

式中,x(t)為隨機共振系統(tǒng)輸出信號;s(t)為輸入信號;a和b為正值的系統(tǒng)參數(shù);n(t)為加性高斯白噪聲,且滿足E[n(t)]=0,D[n(t)]=σ2。

定義V(x)為雙穩(wěn)態(tài)系統(tǒng)的勢函數(shù),當不存在周期信號s(t)時,其表達式為:

,一個勢壘點xb=0,如圖1所示。但是當系統(tǒng)輸入周期信號時,勢函數(shù)會隨著周期信號的變化而發(fā)生變化,如圖2所示。

從力學角度出發(fā),式(1)可用一個布朗粒子在雙穩(wěn)態(tài)勢阱中的運動來形象地解釋。當周期信號不存在時,勢函數(shù)左右對稱,噪聲驅(qū)動粒子在其中一個勢阱內(nèi)隨機波動,或者能夠躍過勢壘以克萊默斯(Kramers)速率在兩個勢阱之間來回穿梭。當加入周期信號后,勢函數(shù)的變化導致勢阱具有了非對稱性,粒子從高勢阱躍遷到低勢阱的概率變大。在適量噪聲能量的幫助下,當粒子來回躍遷周期與信號周期一致時,就出現(xiàn)了噪聲、微弱信號和粒子運動協(xié)同的現(xiàn)象,這就是隨機共振。

圖1 雙穩(wěn)系統(tǒng)勢函數(shù)圖

圖2 雙穩(wěn)系統(tǒng)勢函數(shù)變化圖

在研究的早期階段,噪聲調(diào)節(jié)隨機共振得到了較多關(guān)注,即通過改變噪聲強度來實現(xiàn)非線性系統(tǒng)的隨機共振。而在實際應(yīng)用中,噪聲調(diào)節(jié)隨機共振方法可操作性較差,并且當噪聲的強度已經(jīng)超過發(fā)生隨機共振的閾值條件時,進一步增加噪聲強度只會使系統(tǒng)性能惡化,而減小噪聲強度在實際應(yīng)用中又難以實現(xiàn)。這樣,通過調(diào)節(jié)參數(shù)來實現(xiàn)隨機共振的理論[14]應(yīng)運而生。

Xu[14]等于2002年提出了參數(shù)調(diào)節(jié)隨機共振理論。簡單地說,參數(shù)調(diào)節(jié)隨機共振就是對于確定的輸入信號和加性噪聲,通過調(diào)節(jié)非線性系統(tǒng)的參數(shù),逐步改變最小勢高的值,促使發(fā)生隨機共振。Duan等[14]的研究證明,增加噪聲產(chǎn)生的隨機共振現(xiàn)象和通過參數(shù)調(diào)節(jié)隨機共振的現(xiàn)象具有一致性,調(diào)節(jié)參數(shù)相應(yīng)地調(diào)節(jié)噪聲的大小。

仍然考察雙穩(wěn)態(tài)系統(tǒng),對系統(tǒng)作如式(3)變換

則方程(1)變?yōu)?/p>

式(4)所示的雙穩(wěn)態(tài)系統(tǒng)的勢函數(shù)為:

勢壘的高度為

式(6)說明,增大噪聲強度σ和增大系統(tǒng)參數(shù)b均可降低系統(tǒng)勢壘的高度,二者具有相同的調(diào)節(jié)效果。參數(shù)調(diào)節(jié)可以更靈活地改變最小勢高,進而加速粒子實現(xiàn)成功躍遷的概率,隨機共振現(xiàn)象就得以發(fā)生。參數(shù)調(diào)節(jié)隨機共振既可以加強非線性系統(tǒng)的作用,又可以減弱非線性系統(tǒng)的作用。

3 兩種基于隨機共振的微弱信號接收方法

強噪聲背景下微弱信號的接收一直是通信中亟需解決的問題。常規(guī)的微弱信號檢測方法主要是基于時域和頻域兩種,有時域的相關(guān)方法、取樣積分方法和頻域的譜分析方法等。然而,這些方法都有一定的局限性,主要表現(xiàn)在所能檢測到的微弱信號的信噪比門限值較高。因此,人們一直在尋找新的微弱信號接收方法以彌補傳統(tǒng)方法的不足。隨機共振現(xiàn)象中噪聲增強信號的特性給研究人員提供了新的微弱信號接收的思路,對于噪聲由傳統(tǒng)的抑制噪聲變?yōu)殡S機共振處理中利用噪聲。將噪聲N(t)和信號S(t)同時通入非線性處理器,調(diào)節(jié)系統(tǒng)參數(shù)使非線性系統(tǒng)、信號和噪聲產(chǎn)生隨機共振,繼而從系統(tǒng)輸出信號y(t)中提取信號。

基于現(xiàn)有的相關(guān)研究工作[10][20],將隨機共振理論應(yīng)用于微弱信號接收有兩種基本的方法和途徑:單頻微弱信號隨機共振接收;基帶微弱信號隨機共振接收。兩種方法分別描述如下:

根據(jù)單頻信號的頻譜特點,在對單頻信號進行接收中,往往對隨機共振系統(tǒng)的輸出信號 進行頻譜分析,并以此來恢復單頻信號得到 。具體原理框圖如圖3所示。

圖3 基于隨機共振的單頻信號接收原理框圖

將噪聲N(t)和信號S(t)同時通入隨機共振模塊,通過調(diào)節(jié)系統(tǒng)參數(shù)a和b使系統(tǒng)達到隨機共振,然后對隨機共振模塊的輸出信號進行頻譜分析得到單頻微弱信號的頻率特性進而恢復單頻信號,隨機共振模塊的原理框圖如圖4。圖中,(A)子模塊是加法器,能夠?qū)斎胄盘柤拜敵鲂盘柕姆答伅B加;(B)子模塊是積分器,起到對疊加信號的積分作用;(C)子模塊是一個運算器,對輸出信號進行立方運算并且取反;(D)(E)子模塊均是倍乘器??梢?,輸入信號與輸出信號的a次倍以及輸出信號取反后立方值的b次倍共同經(jīng)過積分,最終得到經(jīng)隨機共振作用的輸出信號。這與朗之萬方程是吻合的。

圖4 參數(shù)調(diào)節(jié)隨機共振模塊原理框圖

為了驗證方法的有效性,本文選取頻率f=20Hz的單頻信號,在信噪比SNR=-18dB條件下進行了仿真,時域仿真結(jié)果如圖5。

圖5 單頻信號隨機共振處理時域效果圖

仿真結(jié)果顯示,通過隨機共振處理之后的加噪單頻信號,部分噪聲能量轉(zhuǎn)化為信號能量,淹沒在強噪聲下的信號特征恢復效果明顯。

(2)基帶微弱信號隨機共振接收

早期隨機共振處理的信號僅僅局限于周期信號,20世紀90年代初,Collins[19]提出了非周期隨機共振的理論,對描述系統(tǒng)特性的朗之萬方程進行了修改如式(7),取得了很好的效果。

非周期隨機共振理論的提出使得隨機共振與信息理論結(jié)合起來,人們將隨機共振用于強噪聲背景的微弱基帶信號接收,取得了很好的效果,基本原理框圖如圖6。

印度消費的石油有多達80%為進口石油,這使得它面對價格波動時比其他進口國更加脆弱。最近幾個月,由于盧比貶值導致石油進口價格膨脹,更加劇了其脆弱性。八月,印度政府估算,如果石油價格居高不下,2018/2019財年的進口石油支出將上升至高達260億美元。

圖6 基于隨機共振的基帶信號接收原理圖

將噪聲N(t)和信號S(t)同時通入非周期隨機共振模塊,通過調(diào)節(jié)系統(tǒng)參數(shù)A使系統(tǒng)達到隨機共振,然后對隨機共振輸出信號進行積分判決恢復信息序列。為了適應(yīng)非周期信號特性,對隨機共振模塊進行了相應(yīng)的修改,具體結(jié)構(gòu)如圖7,取消了立方取反模塊之后的倍乘器,在輸入端增加(D)相乘取反模塊,這與修改后的朗之萬方程相對應(yīng)。

圖7 非周期隨機共振模塊原理框圖

圖8 基于隨機共振的PAM信號接收

本文對比特率r=1000b/s,噪聲方差σ2=16的PAM信號進行了仿真,具體結(jié)果如圖8。仿真結(jié)果表明經(jīng)過隨機共振處理之后PAM信號得到了很好的恢復。

4 需突破的關(guān)鍵技術(shù)與初步思路

雖然隨機共振現(xiàn)象應(yīng)用于微弱信號檢測可取得優(yōu)于傳統(tǒng)微弱信號檢測的效果,但是發(fā)生隨機共振現(xiàn)象有著苛刻的前提條件:信號必須是慢變信號并且振幅不能過大。因此,為了將隨機共振理論應(yīng)用于實際通信、導航或其他信息傳輸系統(tǒng),必須突破以下三個關(guān)鍵技術(shù):為了突破傳統(tǒng)隨機共振小參數(shù)條件的限制,必須研究新的變尺度隨機共振算法;為了應(yīng)對信號未知且多變的實際情況,實現(xiàn)對信號的快速捕獲與連續(xù)跟蹤,必須研究適用于信號傳輸?shù)淖赃m應(yīng)隨機共振算法;針對通信系統(tǒng)中常用的調(diào)制方式(如PSK,F(xiàn)SK,MSK等),必須研究相應(yīng)的調(diào)制信號隨機共振處理方法。

(1)大參數(shù)信號變尺度隨機共振算法

針對通信信號是大參數(shù)的特點,研究將大參數(shù)信號變換到小參數(shù)信號從而觸發(fā)隨機共振的變尺度算法。在研究二次采樣算法[17][18]的基礎(chǔ)上,提出一種基于互質(zhì)欠采樣的變尺度隨機共振算法,使該算法適用于通信系統(tǒng)隨機共振傳輸系統(tǒng)。

對加噪信號進行采樣,采樣頻率為:

采用互質(zhì)欠采樣的方式對第一次采樣得到的序列進行欠采樣,欠采樣頻率為:

N,M互質(zhì),且滿足:

采用互質(zhì)欠采樣得到的類小信號頻率為:

選擇合適的采樣系數(shù)N和M,可以使類小參數(shù)信號滿足隨機共振條件。對小參數(shù)信號進行隨機共振處理得到頻率信息后,即可再按公式還原大參數(shù)信號。

(2)參數(shù)自適應(yīng)隨機共振技術(shù)

針對通信中信號和噪聲的不確定性問題,研究提出新的自適應(yīng)隨機共振算法。根據(jù)隨機共振的輸出,以輸出信噪比、波形與原始波形相關(guān)性、互信息量等為性能準則,自適應(yīng)地調(diào)節(jié)系統(tǒng)參數(shù)。在可實現(xiàn)計算復雜度的基礎(chǔ)上,能夠完成信號的快速捕獲和動態(tài)跟蹤,滿足通信系統(tǒng)信息可靠傳輸?shù)男枨蟆?/p>

針對已有的自適應(yīng)隨機共振算法冗余度高、計算量大等問題,擬在隨機共振處理前增加自適應(yīng)預(yù)調(diào)節(jié)模塊,根據(jù)對噪聲強度和信號頻率的初步估計調(diào)整系統(tǒng)參數(shù),在隨機共振處理后的自適應(yīng)再調(diào)節(jié)模塊根據(jù)輸出的信號進行再調(diào)節(jié),最終實現(xiàn)隨機共振。具體結(jié)構(gòu)如圖9所示。

圖9 雙自適應(yīng)參數(shù)調(diào)節(jié)隨機共振示意圖

(3)通信系統(tǒng)調(diào)制信號隨機共振算法

針對通信系統(tǒng)中典型的調(diào)制方式(PSK,F(xiàn)SK,MSK等),研究與隨機共振結(jié)合進行信息傳輸?shù)目尚行?,分析可獲得的性能增益,通過比較,得出適合于通信系統(tǒng)調(diào)制信號隨機共振的調(diào)制方式。

寬帶調(diào)制信號信息速率高,頻譜帶寬較寬,這使得隨機共振的增強作用并不明顯。據(jù)此,在發(fā)端加入弱導頻或殘留載波輔助,在接收端對導頻信號或殘留載波進行隨機共振處理,以此進行載波恢復,用于調(diào)制信號解調(diào),具體結(jié)構(gòu)如圖10所示。

圖10 有輔助導頻PSK調(diào)制信號隨機共振

窄帶調(diào)制信號瞬時頻譜為窄帶,可直接對窄帶調(diào)制信號進行隨機共振處理,具體結(jié)構(gòu)如圖11所示。

圖11 無輔助導頻FSK調(diào)制信號隨機共振

5 總結(jié)與展望

隨機共振理論在非線性信號處理中發(fā)揮著重要的作用。通過建立參數(shù)調(diào)節(jié)雙穩(wěn)隨機共振模型,文章對現(xiàn)有的基于隨機共振的微弱信號接收的兩種常見途徑分別進行了計算機仿真,仿真結(jié)果表明了基于隨機共振的微弱信號接收方法的有效性和在此應(yīng)用領(lǐng)域的巨大應(yīng)用前景。

經(jīng)典隨機共振理論只能夠適用于慢變?nèi)跣盘枺@大大限制了隨機共振的應(yīng)用范圍。為此,本文針對提出的需要突破大參數(shù)信號變尺度隨機共振、參數(shù)自適應(yīng)隨機共振和調(diào)制信號隨機共振三大關(guān)鍵技術(shù),在現(xiàn)有研究成果的基礎(chǔ)上,分別提出了初步的解決方案,為下一步問題的解決提供了參考。

隨機共振理論和關(guān)鍵技術(shù)的不斷完善勢必將大大拓寬隨機共振應(yīng)用范圍,為探究解決低信噪比條件下的信號處理提供新的思路和途徑。

[1]R.Benzi,A.Sutera,and A.Vulpiani,“The mechanism of stochastic resonance ”. J.Phys. A:Math. Gen. 14(1981) L453-L457

[2]R.Benzi, “Stochastic resonance in climatic change”,Tellus,IEEE,vol.34,pp.10-16,1982

[3]B.McNamara and K.Wiesenfeld,“Theory of stochastic resonance”Physical Review A, vol.39, no.9, pp.4854-4869,May 1989

[4]S.Mitaim and B.Kosko,”Adaptive stochastic resonance”,Proceedings of the IEEE, vol.86, no.11, pp.2152-2183,Nov.1988

[5]F.C.Blondeau,”Input-output gains for signal in noise in stochastic resonance”,Elsevier Physics Lettersm232(1997):41-48

[6]K.Loerincz,Z.Gingl,L.B.Kiss,”A stochastic resonator is able to greatly improve signal-to-noise ratio”Elsevier Physical Letters A,224(1996):63-67

[7]Z.Ginal,L.B.Kiss,and F.Moss,”Non-dynamical stochastic resonance:theory and experiments with white and arbitrarily coloured noise”,Europhys.Lett.29(3),pp.191-196,1995

[8]楊祥龍,汪樂宇.隨機共振技術(shù)在弱信號檢測中的應(yīng)用.電路與系統(tǒng)學報,2001,6(2):94-97

[9]段江海,宋愛國,王一清.隨機共振理論在微弱信號檢測中的應(yīng)用研究.信號處理,2003,19(6) :569-572

[10]梁軍利,楊樹元,唐志峰.基于隨機共振的微弱信號檢測.電子與信息學報,2006,28(6):1068-1072

[11]何大海,趙文禮,梅曉俊.基于隨機共振原理的微弱信號檢測與應(yīng)用.機電工程,2008,25(4 ) :71-74

[12]H.Chen.P.K.Varshney,S.Kay,and J.H.Michels, ”Theory of the stochastic resonance effect in signal detection:part I-fixed detectors” IEEE Trans. Signal process.,vol.55,no.7,pp.3172-3184,July,2007

[13]H.Chen.P.K.Varshney, ”Theory of the stochastic resonance effect in signal detection:part I-variable detectors”, IEEE Trans .Signalprocess.,vol.56,no.10,pp.5031-5041,Oct,2008

[14]Xu Bohou, Duan Fabing, Bao Ronghao, et al. Stochastic resonance with tuning system parameters:the application of bistable systems in signal processing[J]. Chaos, Solitons&Fractals, 2002, 13(4):633-644

[15]F.Duan,D. Abbott,”Binary modulated signal detection in a bistable receiver with stochastic resonance”,Physica A:Statistical Mechanics and its Applications,2007;376:176-190

[16]F.Duan,B.Xu,”Parameter-induced stochastic resonance and baseband binary PAM signals transmission over an AWGN channel”,International Journal of Bifurcation and Chaos,vol.13,no.2(2007)411-425

[17]Leng Y.G.,Leng Y.S.,Wang T.Y.,et al,“Numerical analysis and engineering application of large parameter stochastic resonance”,Journal of Sound and Vibration,2006,292(3-5):788-801

[18]Leng Y.G.,Wang T.Y.,Guo Y.,et a1.,“Engineering signal processing based on bistable stochastic resonance”,Mechanical System and Signal processing,2007,21:138-150

[19]J.J.Collins,Carson C. Chow,and Thomas T. Imhoff. “Aperiod stochastic resonance in excitable systems, ”Phys. Rev. E,1995,52(4),:R3321-3324

[20]段法兵.參數(shù)調(diào)節(jié)隨機共振在數(shù)字信號傳輸中的應(yīng)用.博士學位論文.浙江大學,2002

Stochastic Resonance Phenomenon and Weak Signal Receiving

Shi Shengchao, Li Guangxia, Li Zhiqiang, Zhu Wenming, Zhang Weitong
(College of Communications Engineering PLA University of Science and Technology, Nanjing, 210007)

On the basis of the phenomenon that Stochastic Resonance (SR) can enhance the signal noise ratio (SNR), the Stochastic Resonance nonlinear Bistable-well model and parameter-induced Stochastic Resonance theory are introduced in this paper. Moreover, two weak signal receiving methods based on Stochastic Resonance are analyzed and simulation results are given. In addition, three key techniques including the scale-transforming SR algorithm, the adaptive parameter-adjusting SR algorithm and the SR algorithm for modulated signal are proposed. The elementary solving ideas are also put forward. The related researches will broaden the application area of Stochastic Resonance and provide novel ways for weak signal receiving under very low signal noise ratio circumstance.

stochastic resonance; weak signal receiving ; low carrier to noise ratio (CNR) ; parameter-induced

10.3969/j.issn.1672-7274.2015.03.003

TN75,TN85

A

1672-7274(2015)03-0013-06

國家自然科學基金資助項目(61032004,91338201);國家高技術(shù)研究發(fā)展計劃(“863”計劃)資助項目(2012AA121605,2012AA01A 503,2012AA01A510)。

猜你喜歡
信號系統(tǒng)
Smartflower POP 一體式光伏系統(tǒng)
信號
鴨綠江(2021年35期)2021-04-19 12:24:18
WJ-700無人機系統(tǒng)
ZC系列無人機遙感系統(tǒng)
北京測繪(2020年12期)2020-12-29 01:33:58
完形填空二則
基于PowerPC+FPGA顯示系統(tǒng)
半沸制皂系統(tǒng)(下)
孩子停止長個的信號
連通與提升系統(tǒng)的最后一塊拼圖 Audiolab 傲立 M-DAC mini
基于LabVIEW的力加載信號采集與PID控制
主站蜘蛛池模板: 国产丝袜啪啪| 精品91在线| 欧美中文字幕一区| 天堂av综合网| 老色鬼欧美精品| 免费不卡视频| 日日拍夜夜操| JIZZ亚洲国产| 国产激情影院| 谁有在线观看日韩亚洲最新视频| 欧美精品亚洲精品日韩专| 国产毛片片精品天天看视频| 丰满的少妇人妻无码区| 国产福利一区在线| 久久精品人人做人人爽| 免费jjzz在在线播放国产| 午夜日b视频| 欧美日韩亚洲综合在线观看| 免费观看精品视频999| 亚洲无码高清一区| 香港一级毛片免费看| 亚洲第一综合天堂另类专| 久久国产乱子| 久久国产拍爱| 色成人综合| 影音先锋丝袜制服| 亚洲一区毛片| 国产91蝌蚪窝| 欧美人与牲动交a欧美精品 | 一本大道AV人久久综合| 制服无码网站| 亚洲综合第一区| 国产SUV精品一区二区| 高清免费毛片| 亚洲美女一级毛片| 日韩毛片基地| 伊人91在线| 色噜噜综合网| 亚洲午夜福利在线| 中文字幕伦视频| 成人免费黄色小视频| 国产一级无码不卡视频| 波多野结衣中文字幕一区| 在线观看91香蕉国产免费| 国产精品林美惠子在线观看| 精品福利视频网| 天天躁夜夜躁狠狠躁图片| 久久www视频| 亚洲日韩第九十九页| 热久久国产| 国产美女主播一级成人毛片| 高潮毛片免费观看| 国产第八页| 亚洲人成高清| 欧美日韩一区二区在线免费观看| 国产免费网址| 无码AV日韩一二三区| 国产精品一区在线麻豆| 国产精品毛片一区视频播| 亚洲国产午夜精华无码福利| 亚洲日韩图片专区第1页| 亚洲第一色网站| 久久亚洲高清国产| 91青草视频| 国产真实乱子伦视频播放| 亚洲第一黄色网| 国产精品精品视频| 在线欧美国产| AV老司机AV天堂| 日本午夜影院| 一本色道久久88亚洲综合| 亚洲综合久久成人AV| 综合色在线| 久久毛片网| 国模沟沟一区二区三区| 欧美a在线看| 精品国产一二三区| 深爱婷婷激情网| 国产97视频在线| 免费毛片全部不收费的| 波多野结衣的av一区二区三区| 亚洲福利视频一区二区|