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

        ?

        基于雙波長(zhǎng)時(shí)域合成技術(shù)的微波光子波形產(chǎn)生*

        2019-09-21 06:35:20麻艷娜王文睿宋開臣于晉龍馬闖張華芳
        物理學(xué)報(bào) 2019年17期
        關(guān)鍵詞:信號(hào)

        麻艷娜 王文睿 宋開臣 于晉龍 馬闖 張華芳

        1) (浙江大學(xué)生物醫(yī)學(xué)工程與儀器科學(xué)學(xué)院, 杭州 310027)

        2) (浙江大學(xué)航空航天學(xué)院, 杭州 310027)

        3) (天津大學(xué)電氣自動(dòng)化與信息工程學(xué)院, 天津 300072)

        1 引 言

        近年來(lái), 多種微波波形, 比如正弦波、方波、三角波和鋸齒波, 被廣泛應(yīng)用于醫(yī)學(xué)成像、超寬帶有線和無(wú)線通信、光學(xué)雷達(dá)和光學(xué)傳感等系統(tǒng)中[1-4],其波形質(zhì)量很大程度上決定了系統(tǒng)的性能.由于采樣率的限制, 采用電學(xué)方式生成的微波信號(hào)具有頻率低和帶寬窄的缺點(diǎn).對(duì)于要求高頻率大帶寬的應(yīng)用場(chǎng)合, 則需要借助光學(xué)器件高速、大帶寬和抗電磁干擾的優(yōu)勢(shì).

        目前, 利用微波光子技術(shù)產(chǎn)生高速微波波形是各國(guó)人員的研究熱點(diǎn)[5-9].其中, 傅里葉光脈沖整形技術(shù)[10-13]是常用的方案, 它通過(guò)頻譜整形技術(shù)與頻率-時(shí)間映射技術(shù)的結(jié)合來(lái)實(shí)現(xiàn), 因此也被稱為譜線逐行脈沖整形.而利用Sagnac環(huán)的偏振調(diào)制也可以產(chǎn)生任意波形微波信號(hào).Liu和Yao[14]利用Sagnac環(huán)路內(nèi)的偏振調(diào)制器結(jié)構(gòu)生成三角波、鋸齒波和方波波形.Gao等[15]通過(guò)環(huán)路相位調(diào)制器控制信號(hào)的強(qiáng)度和相位, 并引入色散補(bǔ)償, 實(shí)現(xiàn)重復(fù)頻率可調(diào)諧的三角波和方波信號(hào).

        此外, 通過(guò)控制光電調(diào)制的調(diào)制參數(shù), 也可以實(shí)現(xiàn)正弦微波信號(hào)到方波和三角波等信號(hào)的轉(zhuǎn)變[16-18].張華芳等[19]利用偏振延時(shí)干涉技術(shù), 通過(guò)調(diào)節(jié)調(diào)制器調(diào)制參數(shù)和和差分延時(shí)線的延時(shí)量,得到了5 GHz的方波和三角波輸出.Jiang等[20,21]采用時(shí)域疊加的方式, 通過(guò)調(diào)節(jié)微波波形譜線間相位關(guān)系, 實(shí)現(xiàn)了方波、三角波和鋸齒波信號(hào)輸出.Chen等[22]通過(guò)集成的偏振復(fù)用雙驅(qū)動(dòng)馬赫-曾德爾調(diào)制方案, 避免了大調(diào)制系數(shù)和光濾波結(jié)構(gòu), 實(shí)現(xiàn)了5 GHz重復(fù)頻率的全占空比的三角波輸出.Li等[23]采用正弦信號(hào)驅(qū)動(dòng)的雙平行馬赫-曾德爾調(diào)制器(MZM)和帶阻濾波器, 不需要特定的調(diào)制系數(shù), 即可實(shí)現(xiàn)全占空比的倍頻三角波輸出.Li等[24]通過(guò)雙偏振正交相移鍵控調(diào)制器方案, 通過(guò)改變X軸相移鍵控的偏置和調(diào)節(jié)偏振控制器, 產(chǎn)生了3和6 GHz的三角波和方波波形.

        本文提出了一種可調(diào)諧的微波多波形產(chǎn)生方法, 通過(guò)雙波長(zhǎng)時(shí)域合成技術(shù), 得到與調(diào)制正弦波信號(hào)同頻的方波、三角波和鋸齒波信號(hào).本文在理論分析微波多波形產(chǎn)生原理的基礎(chǔ)上, 實(shí)驗(yàn)生成了2.5和5.0 GHz重復(fù)頻率的方波、鋸齒波和三角波波形.在調(diào)制器工作帶寬內(nèi), 可以方便地實(shí)現(xiàn)輸出波形重復(fù)頻率的可調(diào)諧, 為任意波形的微波信號(hào)產(chǎn)生提供了一種有效的解決方案.因此, 該方案具有良好的應(yīng)用前景.

        2 理論分析

        微波波形發(fā)生器的原理圖如圖1.兩個(gè)單驅(qū)動(dòng)MZM分別偏置在線性工作點(diǎn)和最低功率點(diǎn)位置,當(dāng)不同波長(zhǎng)的兩束光波在時(shí)域疊加合成時(shí), 可以產(chǎn)生不同的微波波形輸出.因此通過(guò)調(diào)節(jié)調(diào)制器的偏置電壓以及調(diào)制深度, 可以實(shí)現(xiàn)調(diào)制光信號(hào)相位和幅度的控制, 最終實(shí)現(xiàn)光子微波波形產(chǎn)生.

        圖1 基于雙波長(zhǎng)時(shí)域合成技術(shù)的微波波形發(fā)生器原理圖, 圖中LD為激光器, WDM為波分復(fù)用器, PC為偏振控制器, OC為3 dB光耦合器, MZM為馬赫-曾德爾調(diào)制器,ODL為光延時(shí)線, AMP為微波放大器Fig.1.Schematic diagram of the proposed microwave waveform generator based on dual-wavelength time domain synthesis technology.LD, laser diode; WDM, wavelength division multiplexer; PC, polarization controller; OC, 3 dB optical coupler; MZM, Mach-Zehnder modulator; ODL, optical delay line; AMP, amplifier.

        假設(shè)MZM的微波調(diào)制輸入為[17-19]

        其中Vm為調(diào)制幅度,ωm為微波信號(hào)的角頻率.若光載波的角頻率為ω0, 則MZM的輸出光場(chǎng)可被近似描述為

        其中E0為輸出光場(chǎng)幅度,φ=πVbias/Vπ為直流偏置電壓Vbias決定的相位偏移,Vπ為調(diào)制器半波電壓.

        將(2)式進(jìn)行Jacobi Anger擴(kuò)展, 可得到,

        其中Jn為第一類n階的貝塞爾函數(shù);β為MZM的調(diào)制系數(shù), 定義式為β=πVm/Vπ.根據(jù)貝塞爾函數(shù)的特點(diǎn), 高階分量可以忽略不計(jì), (3)式可以只用前三階頻率分量近似表述為

        當(dāng)MZM的電壓偏置固定于線性工作點(diǎn), 即φ=π/2時(shí), MZM的輸出光場(chǎng)可被簡(jiǎn)化為

        該調(diào)制光信號(hào)輸入到光電探測(cè)器后, 由(2)式可得,檢測(cè)到的微波信號(hào)為

        通過(guò)Jacobi Anger擴(kuò)展, 可將MZM的輸出經(jīng)過(guò)PD進(jìn)行光電轉(zhuǎn)換后的微波信號(hào)表述為

        而當(dāng)MZM的電壓偏置為最低功率點(diǎn)時(shí), 將φ=0代入(4)式, 同樣可以得到MZM的輸出經(jīng)過(guò)PD進(jìn)行光電轉(zhuǎn)換后的微波信號(hào):

        因此, 當(dāng)MZM工作在線性工作點(diǎn)時(shí), 產(chǎn)生的微波波形只具有奇數(shù)階頻率, 而工作在最低功率點(diǎn)時(shí),則只具有偶數(shù)階頻率.

        2.1 方波波形產(chǎn)生

        通常情況下, 歸一化方波的傅里葉級(jí)數(shù)展開可表示為

        對(duì)比(9)式和(7)式的前三階頻率分量, 得出實(shí)現(xiàn)方波輸出需要滿足以下條件:

        因此調(diào)制系數(shù)β被設(shè)定為2.30, 此時(shí)產(chǎn)生的微波信號(hào)的前三階分量可被描述為

        根據(jù)(9)式和(11)式, 可以得到圖2所示的仿真結(jié)果, 其中虛線表示DC=0時(shí)理想的方波波形輸出, 實(shí)線表示β=2.30 時(shí)僅包含前三階分量的方波波形, 較好地近似了方波波形輸出.

        圖2 β=2.30 時(shí)生成的三階方波(實(shí)線)與理想方波(虛線)的波形圖Fig.2.Comparison of ideal square waveform (dotted line)and three-order square waveform with β of 2.3 (solid line).

        2.2 鋸齒波產(chǎn)生

        對(duì)于鋸齒波形, 其傅里葉級(jí)數(shù)展開可表示為

        因此, 鋸齒波波形的一階和三階頻率分量與上述方波波形相同, 可以通過(guò)同樣的方式獲取.而波形的二階頻率分量則通過(guò)偏置在最低功率點(diǎn)的MZM2產(chǎn)生.因此, 結(jié)合兩束光信號(hào)即可獲得需要的鋸齒波波形.MZM1的調(diào)制深度仍然設(shè)置為2.30, 假設(shè)MZM2的調(diào)制深度為β2, 忽略3階以上的頻率分量, 則產(chǎn)生的微波信號(hào)的交流分量可被描述為

        其中θ為兩束光信號(hào)之間的相對(duì)傳輸延時(shí)導(dǎo)致的相位誤差, 通過(guò)光可調(diào)延時(shí)線進(jìn)行調(diào)節(jié).為了獲取鋸齒波波形的三階近似描述, 要求θ為 π /2 , 且其中θ可通過(guò)光可調(diào)延時(shí)線進(jìn)行調(diào)節(jié), 二階分量的幅值可通過(guò)兩束光載波之間的光功率比或者M(jìn)ZM2的調(diào)制深度進(jìn)行調(diào)節(jié).而反向?qū)ΨQ的鋸齒波波形產(chǎn)生則通過(guò)調(diào)節(jié)θ至 3 π/2 , 或者設(shè)置MZM2的偏置點(diǎn)至最大值來(lái)實(shí)現(xiàn).

        根據(jù)(12)式和(13)式, 可以得到圖3所示的仿真結(jié)果, 其中虛線表示DC=0時(shí)理想的鋸齒波波形輸出, 實(shí)線表示包含僅前三階分量的鋸齒波波形, 較好地近似了鋸齒波波形輸出.

        2.3 三角波波形產(chǎn)生

        歸一化三角波的傅里葉級(jí)數(shù)展開可被表示為

        同樣將(14)式和(7)式的前三階頻率分量進(jìn)行對(duì)比, 為了滿足諧波幅值之間的關(guān)系, 要求:

        圖3 三階鋸齒波(實(shí)線)與理想鋸齒波(虛線)的波形圖Fig.3.Comparison of ideal sawtooth waveform (dotted line and three-order sawtooth waveform (solid line).

        因此調(diào)制系數(shù)β被設(shè)置為1.51, 則(7)式中的交流分量可以被改寫為

        然而一階和三階分量之間還存在附加的π相移.為了滿足該相位關(guān)系, 三角波的產(chǎn)生需要將雙波長(zhǎng)的光載波與1/4周期時(shí)間延遲相結(jié)合.假設(shè)兩個(gè)信號(hào)具有同樣的調(diào)制深度和相同的光場(chǎng)幅度, 則結(jié)合的信號(hào)可被描述為

        根據(jù)(14)式和(17)式, 可以得到圖4所示的仿真結(jié)果, 其中虛線表示DC=0時(shí)理想的三角波波形輸出, 實(shí)線為β=1.51 時(shí)僅包含前三階分量的三角波波形, 較好地近似了三角波波形輸出.

        圖1所示的基于雙波長(zhǎng)時(shí)域合成技術(shù)的微波光子波形產(chǎn)生方案中生成的方波、鋸齒波和三角波波形與理想波形的前三階分量一致, 僅在高階分量的系數(shù)上存在差異.若僅考慮波形前三階, 如圖2-圖4所示, 可實(shí)現(xiàn)波形的近似輸出.而該方案產(chǎn)生的波形除了含有前三階分量, 還存在高階分量, 因此生成的波形相比于上述仿真更接近理想值.波形的各階分量與理想值的誤差通過(guò)后續(xù)的實(shí)驗(yàn)測(cè)量進(jìn)行分析和討論.

        圖4 β=1.51 時(shí)生成的三角波(實(shí)線)與理想三角波(虛線)的波形圖Fig.4.Comparison of ideal triangle waveform (dotted line)and three-order triangle waveform with β of 1.511 (solid line).

        3 實(shí)驗(yàn)與結(jié)果分析

        圖1為基于雙波長(zhǎng)時(shí)域合成技術(shù)的微波波形發(fā)生器實(shí)驗(yàn)裝置圖.兩個(gè)分布反饋激光器(LD1和LD2)分別產(chǎn)生中心波長(zhǎng)為 1551.2 nm (λ1)和1549.6 nm (λ2)的直流光信號(hào).兩路光信號(hào)的頻率間隔約為200 GHz, 遠(yuǎn)高于產(chǎn)生的微波波形的頻率和光電探測(cè)器的帶寬, 從而保證激光器的拍噪聲不會(huì)對(duì)最終的波形造成影響.其中, LD2的輸出被3 dB光耦合器2均分為2路, 一路與LD1的輸出經(jīng)由3 dB光耦合器1合路后進(jìn)入馬赫-曾德爾型的 LiNbO3光強(qiáng)度調(diào)制器 MZM1(Optilab, IM-1550-12), 另一路則進(jìn)入另一個(gè)馬赫-曾德爾型的LiNbO3光強(qiáng)度調(diào)制器MZM2(Optilab, IM-1550-12).其中, MZM1的偏壓固定在線性工作點(diǎn)上, 而MZM2的偏壓固定在最低功率點(diǎn)上.兩個(gè)強(qiáng)度調(diào)制器由同一個(gè)微波信號(hào)源產(chǎn)生的正弦信號(hào)驅(qū)動(dòng), 驅(qū)動(dòng)功率通過(guò)調(diào)節(jié)可調(diào)微波放大器(JDSU, H301)的增益進(jìn)行調(diào)節(jié).偏振控制器PC1, PC2, PC3通過(guò)調(diào)節(jié)MZM輸入光信號(hào)的偏振態(tài), 調(diào)節(jié)調(diào)制器對(duì)相應(yīng)波長(zhǎng)信號(hào)的插損, 最終實(shí)現(xiàn)調(diào)整各路調(diào)制光信號(hào)之間的相對(duì)強(qiáng)度的作用.在MZM1的輸出端, 通過(guò)波分復(fù)用器將λ1和λ2兩個(gè)波長(zhǎng)的光分開, 可調(diào)光延時(shí)線ODL (OZ, ODL-100)則用于調(diào)節(jié)二者的相對(duì)時(shí)延.

        通過(guò)調(diào)節(jié)調(diào)制器的驅(qū)動(dòng)電壓, 可以得到需要的方波輸出.而通過(guò)將兩個(gè)波長(zhǎng)的光通過(guò)3 dB光耦合器耦合, 并調(diào)節(jié)ODL引入的相對(duì)延時(shí)量, 則可以得到三角波輸出.而將λ1波長(zhǎng)的光與MZM2輸出的λ2波長(zhǎng)的光經(jīng)過(guò)光耦合器4耦合, 并調(diào)節(jié)驅(qū)動(dòng)電壓和相對(duì)延時(shí)量, 則可以得到鋸齒波輸出.系統(tǒng)的輸出端由光電示波器(Agilent, 86100A)觀測(cè)輸出波形.同時(shí)由一個(gè)3 dB帶寬為20 GHz的光電探測(cè)器(Optilab, PD-20)將該光信號(hào)轉(zhuǎn)換至電域, 由電譜儀(Agilent, 8564E)觀測(cè)產(chǎn)生的微波信號(hào)頻譜.

        在實(shí)驗(yàn)過(guò)程中, 首先將微波源輸出頻率設(shè)置為2.5 GHz.調(diào)節(jié)MZM1的驅(qū)動(dòng)功率, 當(dāng)驅(qū)動(dòng)功率為19 dBm時(shí), 在示波器可以觀測(cè)到如圖5(a)和圖5(b)所示的方波波形.圖5(c)為該方波對(duì)應(yīng)的電譜圖, 可以看出電譜中1階分量(2.5 GHz)與3階分量之間的功率差為9.7 dB, 與理想方波波形的頻譜特性基本符合, 5階分量與理想值相差12.7 dB.此外, 電譜中2階分量比1階分量低36.8 dB, 4階分量接近電譜底噪, 偶數(shù)階分量得到了很好的抑制.

        而后保持MZM1的參數(shù)不變, 將MZM2的驅(qū)動(dòng)電壓功率設(shè)置為16.8 dBm, 此時(shí)MZM2的輸出端可以得到如圖6(a)所示的倍頻信號(hào), 該倍頻信號(hào)的電譜圖如圖6(b)所示.由于MZM2工作在調(diào)制曲線的最低點(diǎn), 奇數(shù)階頻率分量被有效抑制, 偶數(shù)階分量, 尤其是5 GHz倍頻分量得到增強(qiáng).

        進(jìn)一步將MZM1的λ1波長(zhǎng)信號(hào)和MZM2的λ2信號(hào)耦合, 并調(diào)節(jié)MZM2調(diào)制深度以調(diào)整兩路信號(hào)的相對(duì)強(qiáng)度關(guān)系, 可以得到如圖7所示的鋸齒波波形, 圖7(b)為對(duì)應(yīng)的鋸齒波電譜圖.電譜中1階分量與2階分量之間的功率差為6.0 dB, 與3階分量之間的功率差為9.7 dB, 與理想鋸齒波波形的頻譜特性基本符合.此外, 4階分量與理想值相差3.8 dB, 5階段分量與理想值相差8.4 dB.因此生成的鋸齒波信號(hào)與理想值符合較好.

        此時(shí), 進(jìn)一步調(diào)節(jié)光延時(shí)線, 經(jīng)過(guò)大約200 ps(π相位)后可以得到如圖8所示的反向?qū)ΨQ鋸齒波波形, 各階電譜分量與圖7基本一致.

        然后降低MZM1的驅(qū)動(dòng)功率至15.5 dBm, 此時(shí)調(diào)節(jié)延時(shí)線的延時(shí)量, 并將兩路波長(zhǎng)強(qiáng)度調(diào)整到近似一致.在耦合器3輸出得到如圖9所示的三角波波形圖.圖9(b)為該三角波對(duì)應(yīng)的電譜圖, 可以看出電譜中1階分量與3階分量之間的功率差為19.3 dB, 與理想三角波波形的頻率特性基本符合,而5階分量比理想值相差5.0 dB.此外, 電譜中2階分量比1階分量低35.7 dB, 4階分量接近電譜底噪, 偶數(shù)階分量得到了很好的抑制.

        圖5 2.5 GHz重復(fù)頻率的方波信號(hào) (a), (b)時(shí)域波形圖; (c)頻譜圖Fig.5.(a), (b) Generated square waveform with a repetition rate of 2.5 GHz, and (c) its spectrum.

        圖6 MZM2的輸出端的倍頻信號(hào) (a) 時(shí)域波形圖; (b) 頻譜圖Fig.6.(a) Generated double-frequency waveform and (b) its spectrum in output port of MZM2.

        圖7 2.5 GHz重復(fù)頻率的鋸齒波信號(hào) (a) 時(shí)域波形圖; (b) 頻譜圖Fig.7.(a) Generated sawtooth waveform with a repetition rate of 2.5 GHz, and (b) its spectrum.

        圖8 2.5 GHz重復(fù)頻率的反向?qū)ΨQ鋸齒波信號(hào) (a)時(shí)域波形圖; (b)頻譜圖Fig.8.(a) Generated reverse sawtooth waveform with a repetition rate of 2.5 GHz, and (b) its spectrum.

        圖9 2.5 GHz重復(fù)頻率的三角波信號(hào) (a)時(shí)域波形圖; (b)頻譜圖Fig.9.(a) Generated triangular waveform with a repetition rate of 2.5 GHz, and (b) its spectrum.

        圖10 5 GHz重復(fù)頻率的(a)方波信號(hào), (b)三角波信號(hào), (c)三角波信號(hào)時(shí)域波形圖Fig.10.Generated (a) square waveform, (b) sawtooth waveform, (c) triangular waveform with a repetition rate of 5 GHz.

        此外, 該微波光子波形產(chǎn)生方案還具有可調(diào)諧性, 即可以通過(guò)改變微波源的頻率靈活地控制產(chǎn)生微波波形的頻率.在實(shí)驗(yàn)中, 將微波源的頻率設(shè)置為5 GHz, 采用上述方式, 通過(guò)調(diào)節(jié)調(diào)制器的驅(qū)動(dòng)電壓、兩路信號(hào)的相對(duì)強(qiáng)度關(guān)系和延時(shí)線的相對(duì)延時(shí)量, 得到重復(fù)頻率為5 GHz的方波、鋸齒波和三角波波形 (見圖10).

        上述實(shí)驗(yàn)中波形的產(chǎn)生主要通過(guò)信號(hào)的時(shí)域疊加, 不需要復(fù)雜的倍頻或?yàn)V波技術(shù), 即可直觀且方便地實(shí)現(xiàn)高頻方波、鋸齒波和三角波輸出, 具有良好的實(shí)用性和重要性.通過(guò)改變微波源的調(diào)制頻率即可以實(shí)現(xiàn)上述微波波形的重復(fù)頻率可調(diào)諧.在實(shí)際測(cè)試中, 該方案的操作帶寬僅由所使用的MZM調(diào)制器和光電探測(cè)器的帶寬決定.

        4 結(jié) 論

        本文提出了一種基于雙波長(zhǎng)時(shí)域合成技術(shù)的微波光子波形產(chǎn)生方案, 首先分析了該方案的理論模型, 而后搭建實(shí)驗(yàn), 通過(guò)對(duì)系統(tǒng)中MZM調(diào)制參數(shù)和延時(shí)線相對(duì)延時(shí)量的調(diào)節(jié), 實(shí)現(xiàn)了2.5 GHz和5 GHz的方波、鋸齒波和三角波波形輸出.方案中波形的產(chǎn)生均通過(guò)信號(hào)的時(shí)域合成直接完成, 不需要考慮諧波分量之間的相位和幅度關(guān)系, 也不需要復(fù)雜的倍頻過(guò)程, 因此可以更直觀、更方便地生成所需的微波波形.而產(chǎn)生的微波波形的重復(fù)頻率主要受到調(diào)制器和光電探測(cè)器帶寬的限制, 因此采用更高帶寬的器件可以產(chǎn)生更高重復(fù)頻率的多波形信號(hào).基于雙波長(zhǎng)時(shí)域合成技術(shù)的微波光子波形產(chǎn)生方案具有結(jié)構(gòu)簡(jiǎn)單、成本較低、操作帶寬大和頻率可調(diào)諧的優(yōu)點(diǎn).

        猜你喜歡
        信號(hào)
        信號(hào)
        鴨綠江(2021年35期)2021-04-19 12:24:18
        完形填空二則
        7個(gè)信號(hào),警惕寶寶要感冒
        媽媽寶寶(2019年10期)2019-10-26 02:45:34
        孩子停止長(zhǎng)個(gè)的信號(hào)
        《鐵道通信信號(hào)》訂閱單
        基于FPGA的多功能信號(hào)發(fā)生器的設(shè)計(jì)
        電子制作(2018年11期)2018-08-04 03:25:42
        基于Arduino的聯(lián)鎖信號(hào)控制接口研究
        《鐵道通信信號(hào)》訂閱單
        基于LabVIEW的力加載信號(hào)采集與PID控制
        Kisspeptin/GPR54信號(hào)通路促使性早熟形成的作用觀察
        成人av毛片免费大全| 国产免费专区| 中文字幕日本特黄aa毛片| 亚洲欧美日韩国产精品专区| 中文字幕亚洲无线码| 伊人精品在线观看| 蜜桃一区二区三区在线看| 日本免费三片在线播放| 老女老肥熟女一区二区| 亚洲av无码成人专区片在线观看| 无码人妻品一区二区三区精99| 538在线啪在线观看| 少妇一级aa一区二区三区片| 亚洲大片一区二区三区四区| av免费观看网站大全| 人成综合视频在线播放| 国产成a人亚洲精品无码樱花 | 亚洲精彩av大片在线观看| 亚洲香蕉av一区二区三区| 中文字幕日韩精品有码视频| 白又丰满大屁股bbbbb| 亚洲精品国产国语| 国产免费精品一品二区三| 久久免费看的少妇一级特黄片| 久久国产成人精品av| 亚洲综合色区另类av| 天码av无码一区二区三区四区 | 北条麻妃国产九九九精品视频| 东方aⅴ免费观看久久av| 试看男女炮交视频一区二区三区| 精品人妻久久av中文字幕| 蜜乳一区二区三区亚洲国产| 国产成人小视频| 国产精品jizz视频| 2017天天爽夜夜爽精品视频| 国产午夜福利av在线麻豆| 色综合天天综合网国产成人网| 欧美成人一区二区三区| 久久综合给合久久狠狠狠9| 亚洲天堂av在线免费看| 成年人干逼视频水好多|