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

    基于MSLDSMC窗的諧波與間諧波參數(shù)估計(jì)算法*

    2023-12-21 10:41:14陳乾明張志威李明皓汪旭明雷可君
    關(guān)鍵詞:旁瓣基頻參數(shù)估計(jì)

    陳乾明,張志威,李明皓,汪旭明,雷可君

    (吉首大學(xué)通信與電子工程學(xué)院,湖南 吉首 416000)

    隨著我國新能源及智能電網(wǎng)技術(shù)的不斷發(fā)展,電力諧波對(duì)電力系統(tǒng)的影響日益擴(kuò)大.為了減小設(shè)備損耗、提高電能質(zhì)量,需要采取有效措施來應(yīng)對(duì)諧波問題[1-3],而準(zhǔn)確檢測(cè)電力諧波的各項(xiàng)參數(shù)是解決諧波問題中的重要一環(huán).目前,應(yīng)用較多的電力諧波檢測(cè)方法有瞬時(shí)無功功率理論[4]、小波變換[5]、人工神經(jīng)網(wǎng)絡(luò)[6]及快速傅里葉變換(Fast Fourier Transform,FFT)[7-8]等.其中,FFT算法具有估計(jì)精度高、計(jì)算簡單及運(yùn)算速度快等優(yōu)點(diǎn),廣泛應(yīng)用于諧波檢測(cè)中.然而,在非同步采樣過程中,該算法不可避免地會(huì)出現(xiàn)柵欄效應(yīng)和諧波泄露,降低諧波估計(jì)的精度.針對(duì)FFT算法的不足,通常利用窗函數(shù)來抑制頻譜泄漏,同時(shí)結(jié)合多種插值算法以減少柵欄效應(yīng)造成的誤差[9-10].常用的窗函數(shù)有Hanning窗、Blackman窗、Nuttall窗、Rife-Vincent窗及Kaiser窗等[11].特性優(yōu)良的窗函數(shù)具有主瓣寬度窄、旁瓣峰值電平低和旁瓣衰減速率快等特點(diǎn),而旁瓣特性決定了窗函數(shù)抑制頻譜泄漏的能力.為了進(jìn)一步提升窗函數(shù)的旁瓣性能,學(xué)者相繼構(gòu)建了自乘窗[12]及自卷窗[13].這2類窗函數(shù)的旁瓣性能大幅提升,抑制頻譜泄漏的能力得到進(jìn)一步改善,但同時(shí)也存在一些問題:自乘窗的主瓣寬度會(huì)隨著自乘計(jì)算次數(shù)的增加而增加,導(dǎo)致頻率分辨率降低[14];自卷窗的序列長度會(huì)隨著卷積計(jì)算次數(shù)的增加而增加,導(dǎo)致計(jì)算量增大[15].此外,當(dāng)所采用的母窗旁瓣性能不夠理想時(shí),即便經(jīng)過乘積或卷積計(jì)算后,窗函數(shù)的旁瓣性能得到了一定的改善,但估計(jì)精度的提升仍然是有限的[16-17].應(yīng)用較多的插值算法有雙譜線[18]、三譜線[19]、四譜線[20]和六譜線[21]插值算法.它們都是利用目標(biāo)譜線附近多根譜線獲取信息,不過雙譜線插值算法對(duì)譜線信息的利用不夠充分,而四譜線插值與六譜線插值算法盡管利用了更多的譜線信息,但是計(jì)算復(fù)雜度也隨之提高[22].相較而言,三譜線插值算法能兼顧估計(jì)精度和計(jì)算復(fù)雜度的要求[23].

    綜上考慮,筆者擬利用具有優(yōu)越旁瓣性能的6項(xiàng)最大旁瓣衰減(Maximum Side Lobe Decay,MSLD)窗[24]并綜合自乘窗與自卷窗的優(yōu)勢(shì),構(gòu)建MSLD自乘-卷積(MLSD Self-Multiplication-Convolution,MSLDSMC)窗,再在此基礎(chǔ)上,與三譜線插值算法相結(jié)合,設(shè)計(jì)一種新的高精度三譜線插值電力諧波檢測(cè)算法,推導(dǎo)信號(hào)幅值、頻率及相位的估計(jì)公式,以期進(jìn)一步提高諧波與間諧波各項(xiàng)參數(shù)檢測(cè)的精度.

    1 最大旁瓣衰減自乘-卷積窗

    1.1 最大旁瓣衰減窗

    這里:a0=0.246 093 750,a1=0.410 156 250,a2=0.234 375 000,a3=0.087 890 625,a4=0.019 531 250,a5=0.001 953 125.將MSLD窗與幾種經(jīng)典余弦組合窗函數(shù)進(jìn)行比較,它們的頻域特性參數(shù)見表1,幅頻特性曲線如圖1所示.

    表1 MSLD窗與幾種經(jīng)典余弦組合窗函數(shù)的頻域特性對(duì)比Table 1 Frequency Domain Characteristics Comparison of MSLD Window Function and Several Classic Cosine Combination Window Functions

    圖1 MSLD窗與幾種經(jīng)典余弦組合窗函數(shù)的幅頻特性對(duì)比Fig. 1 Amplitude-Requency Characteristics Comparison Between MSLD Window Function and Several Classic Cosine Combination Window Functions

    由表1和圖1可知:窗函數(shù)主瓣寬度會(huì)隨著窗函數(shù)項(xiàng)數(shù)的增加而增加,頻率分辨率會(huì)相應(yīng)降低;Hanning窗和Blackman窗有較小的主瓣寬度,但旁瓣峰值電平太高,4項(xiàng)1階Nuttall窗旁瓣峰值電平最小,但旁瓣衰減速率太慢,因此它們的旁瓣特性都不夠理想;與4項(xiàng)3階Nuttall窗、4項(xiàng)5階Nuttall窗和5項(xiàng)1階Rife-Vincent窗相比,MSLD窗的旁瓣峰值電平更低,旁瓣衰減速率更快,能更有效地減少頻譜泄漏造成的計(jì)算誤差.因而,本研究采用MSLD窗作為母窗,構(gòu)建MSLDSMC窗對(duì)信號(hào)進(jìn)行加權(quán)處理.

    1.2 MSLDSMC窗的構(gòu)建

    p個(gè)長度為M的MSLD窗,在時(shí)域經(jīng)p-1次乘積計(jì)算后得到p階MSLD自乘窗.考慮到隨著乘積運(yùn)算階數(shù)的增加,窗函數(shù)主瓣寬度也會(huì)增加,因此為了確保頻率分辨率,本研究只進(jìn)行1次乘積運(yùn)算.乘積運(yùn)算不改變窗函數(shù)長度,故2階MSLD自乘(MLSD Self-Multiplication,MSLDSM)窗的時(shí)域表達(dá)式為

    (1)

    其中bi為窗函數(shù)系數(shù),具體取值見表2.

    表2 MSLDSM窗函數(shù)系數(shù)Table 2 Coefficients of MSLDSM Window Function

    為了兼顧算法復(fù)雜度和估計(jì)精度要求,利用wMSLDSM(m)進(jìn)行1次卷積計(jì)算后得到MSLDSMC窗,其時(shí)域表達(dá)式為wMSLDSMC(n)=wMSLDSM(m)*wMSLDSM(m),其中n=0,1,…,N-1,N為MSLDSMC窗函數(shù)的長度.為了便于FFT實(shí)現(xiàn),經(jīng)過卷積計(jì)算后在序列末尾進(jìn)行補(bǔ)0操作,使得離散MSLDSMC窗的長度N=2M.

    1.3 MSLDSMC窗頻域特性

    由(1)式可得MSLDSM窗的頻域表達(dá)式

    其中:ω為角頻率;WR(ω)為矩形窗的離散頻譜函數(shù),

    根據(jù)卷積定理,函數(shù)在時(shí)域進(jìn)行卷積等效于在頻域進(jìn)行乘積.利用離散傅里葉變換性質(zhì),MSLDSMC窗的頻域表達(dá)式為

    WMSLDSMC(ω)=(WMSLDSM(ω))2.

    (2)

    令ω=2πk/N,k為離散頻譜線的序號(hào),k=0,1,…,N-1.由于在實(shí)際應(yīng)用中N通常遠(yuǎn)遠(yuǎn)大于1,因此(2)式可近似表示為

    其中c=0,±1,±2,….當(dāng)c=±1時(shí)可得,MSLDSMC窗的主瓣寬度BMSLDSMC=44π/M=88π/N.

    選擇長度M=64的MSLD窗構(gòu)建MSLDSM窗和MSLDSMC窗,幅頻特性如圖2所示.

    圖2 MSLD,MSLDSM和MSLDSMC窗幅頻特性Fig. 2 Amplitude-Frequency Characteristics of MSLD.MSLDSM and MSLDSMC Windows

    由圖2可知:MSLDSMC窗與MSLDSM窗的主瓣寬度相等,且相較于MSLD窗的主瓣寬度更大;MSLDSMC窗的旁瓣峰值(-301.83 dB)遠(yuǎn)低于MSLD窗的旁瓣峰值(-87.94 dB)和MSLDSM窗的旁瓣峰值(-152.26 dB),且MSLDSMC窗具有更快的旁瓣衰減速度.由此可知,MSLDSMC窗相比MSLD窗和MSLDSM窗的旁瓣特性改善很大,能更有效地抑制頻譜泄漏.

    2 基于最大旁瓣衰減自乘-卷積窗的三譜線插值算法

    設(shè)離散化的單頻信號(hào)為

    (3)

    其中A0,f0,θ0分別為信號(hào)的幅值、頻率和相位,fs為采樣頻率.由于復(fù)譜序列具有對(duì)稱性,因此一般忽略負(fù)頻點(diǎn)峰值的影響.對(duì)(3)式疊加MSLDSMC窗后進(jìn)行傅里葉變換,可得

    (4)

    其中ω0=2πf0/fs.將(4)式離散化,令ω=2πk/N,ω0=2πk0/N,可得

    (5)

    (6)

    可以看出γ是關(guān)于ε的函數(shù),記為γ=g(ε).理論上,通過求解(6)式可得ε=g-1(γ),但這種操作使得計(jì)算復(fù)雜度高,不利于實(shí)時(shí)估計(jì),而采用最小二乘法得到ε的多項(xiàng)式估計(jì)式的計(jì)算復(fù)雜度大大降低.ε的表達(dá)式為

    ε=g-1(γ)=5.505 949 177 228 290 0γ-0.360 510 245 886 963 0γ3+

    0.047 234 928 916 999 2γ5-0.007 707 827 943 080 4γ7,

    從而可確定信號(hào)頻率的修正公式為f0=k0Δf=(k1+ε)fs/N,信號(hào)相位的修正公式為

    θ0=arg(X(k1))-π(ε-0.5).

    利用y1,y2,y3作幅值修正,因?yàn)閗1最靠近真實(shí)峰值頻譜,對(duì)參數(shù)估計(jì)的影響最大,所以對(duì)這根譜線賦予權(quán)值2,而對(duì)譜線k2,k3分別賦予權(quán)值1.結(jié)合(5)式,可得信號(hào)幅值的修正公式

    (7)

    考慮到(7)式是關(guān)于ε的復(fù)雜函數(shù),不利于實(shí)時(shí)計(jì)算,故仍利用最小二乘法求得簡化后的多項(xiàng)擬合式A0=N-2(y2+2y1+y3)h(ε),其中

    h(ε)=65.954 764 271 573 910+2.995 806 576 073 435ε2+0.069 064 953 701 636 3ε4+

    0.001 083 987 688 066 51ε6.

    3 仿真分析

    3.1 整次諧波信號(hào)參數(shù)估計(jì)實(shí)驗(yàn)

    為了驗(yàn)證新算法的有效性,對(duì)含有14次諧波成分的弱幅值信號(hào)進(jìn)行參數(shù)估計(jì).假設(shè)該信號(hào)表達(dá)式為

    (8)

    其中:f為基波頻率,f=50.1 Hz;fs為采樣頻率,fs=2 500 Hz;N為采樣長度,N=2 048;h為諧波次數(shù);Ah,θh分別為諧波信號(hào)的幅值和相位.具體參數(shù)設(shè)置見表3.

    表3 諧波信號(hào)參數(shù)設(shè)置Table 3 Harmonic Signal Parameter Setting

    將MSLDSMC窗與采樣長度相同的MSLD窗和MSLDSM窗分別對(duì)(8)式所示的信號(hào)進(jìn)行仿真分析,并采用三譜線插值算法進(jìn)行諧波參數(shù)估計(jì),估計(jì)所得的基波各次諧波信號(hào)的幅值、頻率和相位的相對(duì)誤差如圖3所示.

    圖3 復(fù)雜諧波信號(hào)幅值、頻率和相位的相對(duì)誤差Fig. 3 Relative Errors of Amplitude,Frequency and Phase of Complex Harmonic Signals

    由圖3可知:在相同條件下,采用新算法對(duì)諧波進(jìn)行分析,估計(jì)所得的幅值、頻率和相位的相對(duì)誤差絕對(duì)值范圍分別為10-13~10-11(圖3(a)),10-16~10-14(圖3(b)),10-15~10-10(圖3(c)),與MSLD窗和MSLDSM窗的三譜線插值算法相比具有更高的估計(jì)精度;圖3(b)中,新算法估計(jì)所得第3次諧波的頻率相對(duì)誤差為0;第11,13,14次諧波的幅值都偏小,但是新算法估計(jì)所得的幅值相對(duì)誤差絕對(duì)值分別為1.21×10-11,1.18×10-11,3.55×10-12,說明新算法估計(jì)整次弱幅值信號(hào)時(shí)精度較高.

    3.2 間諧波信號(hào)參數(shù)估計(jì)實(shí)驗(yàn)

    電網(wǎng)在實(shí)際運(yùn)行中產(chǎn)生的間諧波不僅會(huì)造成電壓、電流波動(dòng),還會(huì)干擾用電設(shè)備的正常運(yùn)行,危及電力系統(tǒng)的安全運(yùn)行,故進(jìn)行間諧波參數(shù)的高精度估計(jì)非常有必要.為了驗(yàn)證新算法估計(jì)間諧波參數(shù)的有效性,仿真過程選用6次間諧波組成的信號(hào),信號(hào)表達(dá)式為

    (9)

    其中:p為諧波次數(shù);Ap為諧波幅值;fp為諧波頻率;θp為諧波相位;fs為采樣頻率,fs=2 500 Hz.間諧波信號(hào)幅值、頻率和相位的具體參數(shù)設(shè)置見表4.

    表4 間諧波信號(hào)參數(shù)設(shè)置Table 4 Inter Harmonic Signal Parameter Setting

    利用新算法和基于MSLD窗三譜線插值算法分別對(duì)(9)式描述的信號(hào)進(jìn)行參數(shù)估計(jì),仿真結(jié)果見表5和表6.

    表5 新算法估計(jì)間諧波信號(hào)各參數(shù)的相對(duì)誤差Table 5 New Algorithm for Estimating the Relative Error of Parameters of Inter-Harmonic Signals

    表6 基于MSLD窗三譜線插值算法估計(jì)間諧波信號(hào)各參數(shù)的相對(duì)誤差Table 6 Three Spectral Line Interpolation Algorithm Based on MSLD Window for Estimating the Relative Error of Parameters of Inter-Harmonic Signals

    對(duì)比表5和表6可知:與基于MSLD窗三譜線插值算法相比,新算法估計(jì)所得的幅值相對(duì)誤差降低了2個(gè)數(shù)量級(jí),頻率相對(duì)誤差降低了2~3個(gè)數(shù)量級(jí),相位相對(duì)誤差數(shù)量級(jí)降低了2~4個(gè)數(shù)量級(jí);頻率相對(duì)誤差結(jié)果中,新算法估計(jì)所得的第4次間諧波信號(hào)頻率的相對(duì)誤差為0.由此可見,新算法估計(jì)間諧波信號(hào)參數(shù)時(shí)精度較高.

    3.3 基頻變動(dòng)實(shí)驗(yàn)

    信號(hào)基波頻率(基頻)的變動(dòng)會(huì)導(dǎo)致非同步采樣發(fā)生,進(jìn)而影響信號(hào)諧波估計(jì)精度.因此,驗(yàn)證新算法在基頻變動(dòng)情況下參數(shù)估計(jì)的穩(wěn)定性是十分重要的.電力系統(tǒng)基頻偏差的最大范圍為0.5 Hz,故設(shè)頻率fi是f在范圍49.5~50.5 Hz內(nèi)以步長為0.1 Hz變動(dòng).對(duì)(8)式描述的諧波信號(hào)進(jìn)行參數(shù)估計(jì),各次諧波信號(hào)參數(shù)的相對(duì)誤差如圖4所示.

    圖4 基頻變動(dòng)時(shí)幅值、頻率和相位的相對(duì)誤差Fig. 4 Relative Errors of Amplitude,Frequency and Phase with Fundamental Frequency Changes

    由圖4可知,在基頻變動(dòng)的情況下,各次諧波信號(hào)幅值、頻率和相位的相對(duì)誤差絕對(duì)值分別不超過1.20×10-11(圖4(a))、9.13×10-15(圖4(b))和2.89×10-10(圖4(c)).可見,新算法在基頻變動(dòng)時(shí)同樣估計(jì)精度較高.

    5 結(jié)語

    鑒于MSLD窗比經(jīng)典的余弦窗具有更優(yōu)異的旁瓣性能,筆者利用乘積和卷積運(yùn)算構(gòu)建了MSLDSMC窗,并設(shè)計(jì)了一種基于MSLDSMC窗三譜線插值的諧波與間諧波參數(shù)估計(jì)算法.與MSLD窗函數(shù)相比,MSLDSMC窗具有更低的旁瓣峰值和更快的旁瓣衰減速度.整次諧波信號(hào)參數(shù)估計(jì)實(shí)驗(yàn)結(jié)果表明,與基于MSLD窗和MSLDSM窗的三譜線插值諧波與間諧波參數(shù)估計(jì)算法相比,新算法對(duì)整次諧波信號(hào)參數(shù)的估計(jì)精度明顯提高.間諧波信號(hào)參數(shù)估計(jì)實(shí)驗(yàn)結(jié)果表明,與基于MSLD的三譜線插值諧波參數(shù)估計(jì)算法相比,新算法對(duì)間諧波信號(hào)參數(shù)的估計(jì)精度更高.基頻變動(dòng)實(shí)驗(yàn)結(jié)果表明,新算法對(duì)間諧波信號(hào)參數(shù)和在基頻變動(dòng)環(huán)境下對(duì)整次諧波信號(hào)參數(shù)的估計(jì)精度都較高.

    猜你喜歡
    旁瓣基頻參數(shù)估計(jì)
    語音同一認(rèn)定中音段長度對(duì)基頻分析的影響
    基于圓柱陣通信系統(tǒng)的廣義旁瓣對(duì)消算法
    基于新型DFrFT的LFM信號(hào)參數(shù)估計(jì)算法
    基于時(shí)域的基頻感知語音分離方法?
    橋面鋪裝層對(duì)中小跨徑橋梁基頻影響分析
    一種基于線性規(guī)劃的頻率編碼旁瓣抑制方法
    基于加權(quán)積分旁瓣最小化的隨機(jī)多相碼設(shè)計(jì)
    45000kHz基頻晶體濾波器
    電子制作(2017年20期)2017-04-26 06:57:35
    Logistic回歸模型的幾乎無偏兩參數(shù)估計(jì)
    基于向前方程的平穩(wěn)分布參數(shù)估計(jì)
    国产精品一区二区精品视频观看| 丝瓜视频免费看黄片| 在线观看免费视频网站a站| 国产精品乱码一区二三区的特点 | 国产高清国产精品国产三级| 亚洲熟妇熟女久久| 热99国产精品久久久久久7| 人人澡人人妻人| 夜夜夜夜夜久久久久| 欧美 日韩 精品 国产| 亚洲午夜理论影院| 大码成人一级视频| 男人的好看免费观看在线视频 | 欧美日韩一级在线毛片| 99热只有精品国产| 777久久人妻少妇嫩草av网站| 一级毛片女人18水好多| 亚洲精品国产色婷婷电影| 久久国产精品人妻蜜桃| 国产极品粉嫩免费观看在线| 另类亚洲欧美激情| 亚洲欧美精品综合一区二区三区| 亚洲色图综合在线观看| 日本vs欧美在线观看视频| 亚洲国产欧美一区二区综合| 国产精品免费一区二区三区在线 | 热re99久久精品国产66热6| www.熟女人妻精品国产| 老鸭窝网址在线观看| 国产精品二区激情视频| 免费黄频网站在线观看国产| 又紧又爽又黄一区二区| 伦理电影免费视频| 波多野结衣一区麻豆| 极品教师在线免费播放| 亚洲成av片中文字幕在线观看| 国产在线一区二区三区精| 岛国在线观看网站| 国内久久婷婷六月综合欲色啪| 亚洲精华国产精华精| 亚洲性夜色夜夜综合| 久久人人97超碰香蕉20202| 国产成人啪精品午夜网站| 精品国产一区二区三区久久久樱花| 丰满迷人的少妇在线观看| 国产精品一区二区在线观看99| 国产成+人综合+亚洲专区| 高清毛片免费观看视频网站 | 欧美激情 高清一区二区三区| 日日爽夜夜爽网站| 精品乱码久久久久久99久播| av免费在线观看网站| 制服人妻中文乱码| 国产97色在线日韩免费| 一级片免费观看大全| avwww免费| 又紧又爽又黄一区二区| 99久久精品国产亚洲精品| 国产成人精品无人区| 狠狠狠狠99中文字幕| 国产亚洲欧美在线一区二区| 亚洲国产看品久久| av天堂在线播放| 天天影视国产精品| 新久久久久国产一级毛片| 久久婷婷成人综合色麻豆| 中文字幕精品免费在线观看视频| 一级毛片女人18水好多| 超碰成人久久| 一边摸一边抽搐一进一出视频| 午夜激情av网站| 亚洲国产欧美日韩在线播放| 国产xxxxx性猛交| 亚洲男人天堂网一区| 国产伦人伦偷精品视频| 十八禁人妻一区二区| 99久久精品国产亚洲精品| 丰满人妻熟妇乱又伦精品不卡| 在线观看免费高清a一片| 国产精品免费视频内射| 大型黄色视频在线免费观看| 一本一本久久a久久精品综合妖精| 欧美日韩中文字幕国产精品一区二区三区 | 国产精品一区二区在线不卡| 国产激情欧美一区二区| 免费在线观看视频国产中文字幕亚洲| 午夜激情av网站| 999久久久精品免费观看国产| 久久精品国产a三级三级三级| av视频免费观看在线观看| 1024香蕉在线观看| 大陆偷拍与自拍| 国产伦人伦偷精品视频| 一边摸一边抽搐一进一小说 | 亚洲欧洲精品一区二区精品久久久| 久久午夜综合久久蜜桃| 成人18禁在线播放| 日韩一卡2卡3卡4卡2021年| 国产成人av教育| 国产精华一区二区三区| 在线观看一区二区三区激情| 纯流量卡能插随身wifi吗| 国产精品国产高清国产av | 男人舔女人的私密视频| 亚洲精品中文字幕在线视频| 黄色成人免费大全| 欧美不卡视频在线免费观看 | 国产精品二区激情视频| 亚洲国产毛片av蜜桃av| 美女福利国产在线| xxxhd国产人妻xxx| 啦啦啦在线免费观看视频4| 久久ye,这里只有精品| 校园春色视频在线观看| 国产97色在线日韩免费| bbb黄色大片| 性少妇av在线| 十分钟在线观看高清视频www| 国产成人免费观看mmmm| 亚洲精品国产一区二区精华液| aaaaa片日本免费| 亚洲人成电影免费在线| 欧美日韩福利视频一区二区| 免费观看a级毛片全部| 国产成人精品在线电影| 久久精品国产亚洲av高清一级| 国产欧美日韩一区二区三| 亚洲精品成人av观看孕妇| 韩国精品一区二区三区| 亚洲成人免费电影在线观看| 又黄又爽又免费观看的视频| 高清欧美精品videossex| 免费不卡黄色视频| 日韩一卡2卡3卡4卡2021年| 三级毛片av免费| 免费在线观看亚洲国产| 又大又爽又粗| 大陆偷拍与自拍| 亚洲成av片中文字幕在线观看| 久热这里只有精品99| 精品一品国产午夜福利视频| 国产在线观看jvid| 成年动漫av网址| 亚洲欧美激情在线| 下体分泌物呈黄色| 亚洲午夜理论影院| 9色porny在线观看| 丁香欧美五月| ponron亚洲| 精品少妇一区二区三区视频日本电影| 亚洲午夜理论影院| 欧美乱色亚洲激情| 在线观看免费视频日本深夜| 久久久精品区二区三区| 亚洲精品中文字幕在线视频| 精品电影一区二区在线| 国产乱人伦免费视频| 国产成+人综合+亚洲专区| 久久久久精品国产欧美久久久| av一本久久久久| 男男h啪啪无遮挡| 精品国产乱子伦一区二区三区| 国内久久婷婷六月综合欲色啪| 久久人妻av系列| 成人精品一区二区免费| 亚洲av片天天在线观看| 欧美精品啪啪一区二区三区| 午夜激情av网站| 在线观看免费视频日本深夜| 成人特级黄色片久久久久久久| 999精品在线视频| 18禁观看日本| 人人妻,人人澡人人爽秒播| 啦啦啦免费观看视频1| 在线观看一区二区三区激情| 国产xxxxx性猛交| 美女视频免费永久观看网站| 午夜视频精品福利| 一a级毛片在线观看| 一区二区日韩欧美中文字幕| 欧美黄色片欧美黄色片| 久久久水蜜桃国产精品网| 亚洲,欧美精品.| 欧美激情高清一区二区三区| 90打野战视频偷拍视频| 在线看a的网站| 成年人黄色毛片网站| 电影成人av| 桃红色精品国产亚洲av| 极品人妻少妇av视频| 中出人妻视频一区二区| 亚洲avbb在线观看| 精品国内亚洲2022精品成人 | 亚洲欧美一区二区三区久久| 韩国av一区二区三区四区| 日韩免费高清中文字幕av| 亚洲国产欧美网| 不卡av一区二区三区| 91精品国产国语对白视频| 色综合婷婷激情| 97人妻天天添夜夜摸| 亚洲国产欧美日韩在线播放| 极品少妇高潮喷水抽搐| av视频免费观看在线观看| 久久久久久久午夜电影 | 久热这里只有精品99| 国产欧美日韩精品亚洲av| 成在线人永久免费视频| 黄色视频,在线免费观看| 老鸭窝网址在线观看| 18禁黄网站禁片午夜丰满| 中文字幕色久视频| 亚洲综合色网址| 黄色a级毛片大全视频| 最近最新中文字幕大全电影3 | 亚洲av第一区精品v没综合| 国产精品久久久久久人妻精品电影| 亚洲av第一区精品v没综合| 嫁个100分男人电影在线观看| 亚洲色图 男人天堂 中文字幕| x7x7x7水蜜桃| 老司机福利观看| 正在播放国产对白刺激| 亚洲精品国产色婷婷电影| 精品高清国产在线一区| 亚洲自偷自拍图片 自拍| 色老头精品视频在线观看| 国产精品亚洲一级av第二区| 精品人妻1区二区| av欧美777| 人人妻人人澡人人爽人人夜夜| 性色av乱码一区二区三区2| 亚洲av第一区精品v没综合| 中亚洲国语对白在线视频| 青草久久国产| 美女 人体艺术 gogo| 国产亚洲欧美精品永久| 免费一级毛片在线播放高清视频 | 亚洲国产精品sss在线观看 | 免费观看精品视频网站| 国产精品99久久99久久久不卡| 欧美精品亚洲一区二区| 精品国产超薄肉色丝袜足j| 国产在视频线精品| 欧美久久黑人一区二区| 啦啦啦视频在线资源免费观看| 成年人午夜在线观看视频| 每晚都被弄得嗷嗷叫到高潮| 国产野战对白在线观看| 免费少妇av软件| 久久中文看片网| 精品国产乱码久久久久久男人| 后天国语完整版免费观看| 中出人妻视频一区二区| 久久午夜综合久久蜜桃| 女人爽到高潮嗷嗷叫在线视频| 一级黄色大片毛片| 久久狼人影院| 久久香蕉国产精品| a在线观看视频网站| 免费不卡黄色视频| 久久ye,这里只有精品| 精品高清国产在线一区| 老司机亚洲免费影院| 国产高清视频在线播放一区| 亚洲片人在线观看| 51午夜福利影视在线观看| 国产国语露脸激情在线看| 午夜久久久在线观看| 国产激情欧美一区二区| 免费在线观看视频国产中文字幕亚洲| 久久香蕉精品热| 国产亚洲精品第一综合不卡| 久久久精品免费免费高清| 少妇被粗大的猛进出69影院| 天天添夜夜摸| 久9热在线精品视频| 国产激情久久老熟女| 人妻丰满熟妇av一区二区三区 | 欧美精品亚洲一区二区| 午夜福利视频在线观看免费| 天天添夜夜摸| 99香蕉大伊视频| 国产激情久久老熟女| 美女福利国产在线| 女同久久另类99精品国产91| 波多野结衣av一区二区av| 久久中文看片网| 免费女性裸体啪啪无遮挡网站| 在线看a的网站| 午夜福利乱码中文字幕| 欧美激情高清一区二区三区| 另类亚洲欧美激情| 亚洲人成77777在线视频| 免费观看人在逋| 一夜夜www| 母亲3免费完整高清在线观看| 成人手机av| 操出白浆在线播放| 多毛熟女@视频| 在线国产一区二区在线| 999精品在线视频| 人人妻人人爽人人添夜夜欢视频| 一级片免费观看大全| 久久久久视频综合| 岛国在线观看网站| 日韩免费av在线播放| 国产精品99久久99久久久不卡| 婷婷精品国产亚洲av在线 | 久久国产精品人妻蜜桃| 亚洲午夜精品一区,二区,三区| 国产成人系列免费观看| 久久中文字幕人妻熟女| 中出人妻视频一区二区| 亚洲欧美激情在线| 免费高清在线观看日韩| 自拍欧美九色日韩亚洲蝌蚪91| 久久国产精品影院| 亚洲av成人av| e午夜精品久久久久久久| av电影中文网址| 国产淫语在线视频| 黄片大片在线免费观看| 91大片在线观看| 亚洲成人免费电影在线观看| 90打野战视频偷拍视频| 国产亚洲欧美98| 欧美不卡视频在线免费观看 | 亚洲精品久久午夜乱码| 欧美性长视频在线观看| 十八禁高潮呻吟视频| 中文字幕另类日韩欧美亚洲嫩草| 十分钟在线观看高清视频www| 精品高清国产在线一区| ponron亚洲| 免费在线观看完整版高清| 18在线观看网站| 人人澡人人妻人| 亚洲精品一卡2卡三卡4卡5卡| 欧美在线黄色| 51午夜福利影视在线观看| 国产精品成人在线| 成年女人毛片免费观看观看9 | 麻豆成人av在线观看| 久久国产精品男人的天堂亚洲| 男女之事视频高清在线观看| 久久久国产一区二区| 色播在线永久视频| av线在线观看网站| 亚洲精品美女久久久久99蜜臀| 交换朋友夫妻互换小说| 黑人猛操日本美女一级片| 成人av一区二区三区在线看| videosex国产| 亚洲黑人精品在线| 成人国语在线视频| 欧美日韩国产mv在线观看视频| 黑人巨大精品欧美一区二区蜜桃| 高清在线国产一区| 国产精品 欧美亚洲| 丁香六月欧美| 欧美激情久久久久久爽电影 | 亚洲色图综合在线观看| 91在线观看av| 欧美日韩瑟瑟在线播放| 国产视频一区二区在线看| 黄色 视频免费看| 国产在线观看jvid| 国产精品二区激情视频| 99热只有精品国产| 国产精品一区二区在线观看99| 国产精品一区二区在线观看99| 色精品久久人妻99蜜桃| 国产单亲对白刺激| 久久人妻av系列| 夜夜躁狠狠躁天天躁| 两性午夜刺激爽爽歪歪视频在线观看 | 又大又爽又粗| 国产高清videossex| 亚洲欧美激情在线| 免费在线观看完整版高清| 大香蕉久久成人网| 脱女人内裤的视频| 十八禁高潮呻吟视频| 国产精品av久久久久免费| 女人被狂操c到高潮| 欧美日韩精品网址| 国产一卡二卡三卡精品| tocl精华| 在线观看免费高清a一片| 欧美亚洲日本最大视频资源| 免费在线观看完整版高清| 欧美精品一区二区免费开放| 在线国产一区二区在线| 精品亚洲成a人片在线观看| 俄罗斯特黄特色一大片| 大片电影免费在线观看免费| 色老头精品视频在线观看| 日韩成人在线观看一区二区三区| 精品免费久久久久久久清纯 | 午夜精品国产一区二区电影| 国产精品久久视频播放| 老熟妇仑乱视频hdxx| 1024香蕉在线观看| 国产精品二区激情视频| 丝袜人妻中文字幕| 亚洲精品国产一区二区精华液| 精品卡一卡二卡四卡免费| 日韩免费av在线播放| 在线免费观看的www视频| 麻豆av在线久日| 欧洲精品卡2卡3卡4卡5卡区| 午夜福利影视在线免费观看| 久久婷婷成人综合色麻豆| 国产精品乱码一区二三区的特点 | 最新在线观看一区二区三区| 日韩欧美在线二视频 | 欧美另类亚洲清纯唯美| 91精品国产国语对白视频| 一进一出抽搐gif免费好疼 | 久久中文看片网| 两人在一起打扑克的视频| 水蜜桃什么品种好| 精品人妻熟女毛片av久久网站| tocl精华| 免费一级毛片在线播放高清视频 | 国产色视频综合| 精品少妇久久久久久888优播| 伊人久久大香线蕉亚洲五| 大陆偷拍与自拍| 欧美大码av| 日韩中文字幕欧美一区二区| 午夜福利一区二区在线看| 亚洲人成电影观看| 亚洲欧美日韩高清在线视频| 久久精品国产a三级三级三级| 亚洲aⅴ乱码一区二区在线播放 | 欧美日韩亚洲高清精品| 亚洲精品美女久久久久99蜜臀| 精品一区二区三区四区五区乱码| 天天躁狠狠躁夜夜躁狠狠躁| 国产午夜精品久久久久久| 欧美亚洲 丝袜 人妻 在线| 十八禁网站免费在线| 天天操日日干夜夜撸| 91老司机精品| 一进一出抽搐gif免费好疼 | 国产高清激情床上av| 欧美精品人与动牲交sv欧美| 成人18禁高潮啪啪吃奶动态图| 欧美日韩亚洲高清精品| 成在线人永久免费视频| 一级片'在线观看视频| 日本黄色视频三级网站网址 | 成人18禁高潮啪啪吃奶动态图| 亚洲全国av大片| 大型av网站在线播放| 欧美久久黑人一区二区| 一边摸一边做爽爽视频免费| 日本黄色日本黄色录像| 大码成人一级视频| 变态另类成人亚洲欧美熟女 | 欧美另类亚洲清纯唯美| 老汉色∧v一级毛片| 天堂中文最新版在线下载| 精品第一国产精品| 少妇粗大呻吟视频| 精品免费久久久久久久清纯 | 丝瓜视频免费看黄片| 最新美女视频免费是黄的| 飞空精品影院首页| 午夜福利,免费看| 后天国语完整版免费观看| 免费一级毛片在线播放高清视频 | 亚洲欧美一区二区三区黑人| 国产亚洲欧美精品永久| 国产精品亚洲一级av第二区| 黑丝袜美女国产一区| 成人18禁在线播放| 嫩草影视91久久| 久久天躁狠狠躁夜夜2o2o| 日日摸夜夜添夜夜添小说| 日韩精品免费视频一区二区三区| 中文字幕人妻熟女乱码| 黄色毛片三级朝国网站| 一级a爱片免费观看的视频| 在线观看免费日韩欧美大片| 99久久99久久久精品蜜桃| 99精品欧美一区二区三区四区| 久久这里只有精品19| 一级黄色大片毛片| 90打野战视频偷拍视频| √禁漫天堂资源中文www| 亚洲国产欧美日韩在线播放| 俄罗斯特黄特色一大片| 欧美日韩成人在线一区二区| 精品少妇久久久久久888优播| 免费黄频网站在线观看国产| 一边摸一边抽搐一进一小说 | tube8黄色片| 夜夜爽天天搞| 亚洲情色 制服丝袜| 午夜影院日韩av| 桃红色精品国产亚洲av| 少妇猛男粗大的猛烈进出视频| 日韩 欧美 亚洲 中文字幕| 国产成人影院久久av| 国产在线精品亚洲第一网站| 午夜福利视频在线观看免费| 久久精品国产亚洲av高清一级| 久久九九热精品免费| 国产真实乱freesex| 国产aⅴ精品一区二区三区波| 成人永久免费在线观看视频| 久久精品人妻少妇| 欧美午夜高清在线| 欧美黑人巨大hd| tocl精华| 亚洲在线自拍视频| 欧美日韩瑟瑟在线播放| 此物有八面人人有两片| 少妇丰满av| 丰满乱子伦码专区| 亚洲精品影视一区二区三区av| 性欧美人与动物交配| 欧美最黄视频在线播放免费| 一卡2卡三卡四卡精品乱码亚洲| 亚洲av成人精品一区久久| 免费看光身美女| 欧美一区二区精品小视频在线| 91av网一区二区| 日本免费a在线| 99久久精品国产亚洲精品| 亚洲成人免费电影在线观看| 久久精品91无色码中文字幕| 9191精品国产免费久久| 色综合亚洲欧美另类图片| 婷婷丁香在线五月| 国内少妇人妻偷人精品xxx网站| 国产欧美日韩一区二区精品| 精品福利观看| 久久亚洲真实| 国产精品野战在线观看| 国产伦精品一区二区三区视频9 | 亚洲国产色片| 久久精品亚洲精品国产色婷小说| 99国产极品粉嫩在线观看| 欧美xxxx黑人xx丫x性爽| 日韩欧美国产在线观看| 中文亚洲av片在线观看爽| 极品教师在线免费播放| 成年女人毛片免费观看观看9| 国产乱人伦免费视频| 午夜福利免费观看在线| 婷婷丁香在线五月| 日韩中文字幕欧美一区二区| 波多野结衣高清作品| 亚洲不卡免费看| 亚洲国产精品sss在线观看| 欧美日韩福利视频一区二区| 国产老妇女一区| 性色av乱码一区二区三区2| 可以在线观看毛片的网站| 日本精品一区二区三区蜜桃| 一本精品99久久精品77| 久久久成人免费电影| 国产精品久久电影中文字幕| 精品一区二区三区人妻视频| 黄色视频,在线免费观看| 看片在线看免费视频| 欧美在线黄色| 亚洲av不卡在线观看| 亚洲一区二区三区色噜噜| 天天躁日日操中文字幕| 天堂影院成人在线观看| ponron亚洲| 午夜精品在线福利| 亚洲国产日韩欧美精品在线观看 | 亚洲av日韩精品久久久久久密| 国产成人啪精品午夜网站| 欧美三级亚洲精品| 岛国在线免费视频观看| 窝窝影院91人妻| av欧美777| 少妇的丰满在线观看| 在线观看av片永久免费下载| 色精品久久人妻99蜜桃| 国产精品国产高清国产av| 变态另类丝袜制服| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 91av网一区二区| 欧美日韩福利视频一区二区| 国产高潮美女av| 免费搜索国产男女视频| 国产欧美日韩精品一区二区| 99热只有精品国产| 又爽又黄无遮挡网站| 精品不卡国产一区二区三区| 亚洲黑人精品在线| 午夜精品一区二区三区免费看| 老熟妇仑乱视频hdxx| 在线观看美女被高潮喷水网站 | www.色视频.com| 久久久久九九精品影院| 18+在线观看网站| 国产真人三级小视频在线观看| 国产美女午夜福利| 真人做人爱边吃奶动态| 久久久久久国产a免费观看| 99热6这里只有精品| 丁香六月欧美| 午夜福利欧美成人| 超碰av人人做人人爽久久 | 动漫黄色视频在线观看| 在线a可以看的网站|