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

    基于自適應(yīng)EEMD和盲辨識(shí)算法的橋梁結(jié)構(gòu)模態(tài)參數(shù)識(shí)別

    2016-05-22 02:33:23
    關(guān)鍵詞:盲源幅值預(yù)處理

    陳 永 高

    (浙江工業(yè)職業(yè)技術(shù)學(xué)院,浙江 紹興 312000)

    基于自適應(yīng)EEMD和盲辨識(shí)算法的橋梁結(jié)構(gòu)模態(tài)參數(shù)識(shí)別

    陳 永 高

    (浙江工業(yè)職業(yè)技術(shù)學(xué)院,浙江 紹興 312000)

    基于集合經(jīng)驗(yàn)?zāi)B(tài)分解(EEMD)在信號(hào)的預(yù)處理上的不足之處,提出了一種基于自適應(yīng)EEMD分解的盲源分離算法,即:先根據(jù)原始信號(hào)自身的特點(diǎn)確定加入白噪聲的幅值標(biāo)準(zhǔn)差和集成次數(shù),再進(jìn)行EEMD分解,對(duì)所得IMF分量進(jìn)行模糊綜合評(píng)價(jià)以選出有效的IMF分量,構(gòu)造IMF分量矩陣,最后利用盲源分離算法對(duì)其進(jìn)行盲辨識(shí),完成對(duì)信號(hào)的分解與重構(gòu)。分別通過(guò)模擬信號(hào)和橋梁實(shí)測(cè)振動(dòng)信號(hào)對(duì)該算法進(jìn)行驗(yàn)證。結(jié)果表明:所提算法具有可行性,且能運(yùn)用于實(shí)際結(jié)構(gòu)信號(hào)的預(yù)處理。

    橋梁工程;自適應(yīng)EEMD;盲辨識(shí);模糊綜合評(píng)價(jià)法;參數(shù)識(shí)別

    0 引 言

    橋梁結(jié)構(gòu)作為我國(guó)基礎(chǔ)設(shè)施的重要組成部分之一,隨著經(jīng)濟(jì)的飛速發(fā)展,橋梁建設(shè)得到了長(zhǎng)足的進(jìn)步。近年來(lái),隨著自然災(zāi)害發(fā)生次數(shù)的不斷增加,橋梁結(jié)構(gòu)的安全性[1]、耐久性與正常使用日漸成為了人們關(guān)注的問(wèn)題之一?,F(xiàn)階段,國(guó)內(nèi)外學(xué)術(shù)界以及工程界已經(jīng)將橋梁結(jié)構(gòu)的健康監(jiān)測(cè)和智能控制作為了重點(diǎn)研究對(duì)象。工程中,可以通過(guò)對(duì)橋梁結(jié)構(gòu)進(jìn)行模態(tài)參數(shù)識(shí)別[2]來(lái)達(dá)到對(duì)橋梁結(jié)構(gòu)的健康監(jiān)測(cè)以及損傷診斷,而在進(jìn)行參數(shù)識(shí)別之前,則需要對(duì)信號(hào)進(jìn)行預(yù)處理。對(duì)于橋梁結(jié)構(gòu)而言,結(jié)構(gòu)的振動(dòng)響應(yīng)是由結(jié)構(gòu)上傳感器記錄得到的,由于傳感器是在環(huán)境激勵(lì)下工作的,以至所測(cè)振動(dòng)響應(yīng)具有幅值小、隨機(jī)性強(qiáng)、容易受噪聲影響和數(shù)據(jù)量巨大的特點(diǎn)?;诖?,筆者提出了基于自適應(yīng)EEMD[3]和盲源分離算法[4]的信號(hào)預(yù)處理方法,并通過(guò)仿真信號(hào)檢驗(yàn)該算法的準(zhǔn)確性與可行性。為了檢驗(yàn)該方法能否運(yùn)用于實(shí)際工程中,筆者最后利用協(xié)方差驅(qū)動(dòng)隨機(jī)子空間算法[5]對(duì)信號(hào)進(jìn)行參數(shù)識(shí)別,并通過(guò)對(duì)比Hilbert-Huang譜和穩(wěn)定圖檢驗(yàn)算法的可行性。

    1 自適應(yīng)集合經(jīng)驗(yàn)?zāi)B(tài)分解

    集合經(jīng)驗(yàn)?zāi)B(tài)分解[3](EEMD),即:一種噪聲輔助數(shù)據(jù)分析方法。具體的算法步驟可見(jiàn)文獻(xiàn)[3]。EEMD是針對(duì)經(jīng)驗(yàn)?zāi)B(tài)分解(EMD)的不足而被提出的一種分解算法。通過(guò)不斷的研究發(fā)現(xiàn),其依然存在以下幾點(diǎn)不足:

    1)需要自定擬白噪聲的幅值標(biāo)準(zhǔn)差;

    2)分解的次數(shù),即EEMD的集成次數(shù)沒(méi)有統(tǒng)一的標(biāo)準(zhǔn),需要人為事先確定;

    3)不能事先確定輸入信號(hào)與本征模態(tài)函數(shù)(IMFs)之間的相對(duì)誤差;

    4)IMF的有效選取需要人為確定,且不能實(shí)現(xiàn)信號(hào)的自我重構(gòu)。

    針對(duì)EEMD存在的不足,筆者提出了自適應(yīng)EEMD分解算法,具體實(shí)現(xiàn)如下。

    1.1 白噪聲的自適應(yīng)選擇

    為了實(shí)現(xiàn)白噪聲的自適應(yīng)選擇,則需要在分解信號(hào)之前分析信號(hào)自身存在的特點(diǎn),以便確定最佳的白噪聲幅值標(biāo)準(zhǔn)差。研究表明,加入的白噪聲應(yīng)滿足兩個(gè)條件[6]:①不能影響原始振動(dòng)信號(hào)中高頻成分極值點(diǎn)的具體分布;②能減小原始振動(dòng)信號(hào)中低頻成分的極值點(diǎn)間隔,并使間隔的分布更為均勻。主要目的在于減小在使用三次樣條函數(shù)進(jìn)行擬合包絡(luò)時(shí),局部均值的誤差。如何根據(jù)原始信號(hào)自適應(yīng)確定加入白噪聲的幅值標(biāo)準(zhǔn)差(σn),具體分析如下:

    Step1: 計(jì)算振動(dòng)信號(hào)的幅值標(biāo)準(zhǔn)差,即σ0。

    Step2: 對(duì)振動(dòng)信號(hào)進(jìn)行高通濾波分解,并計(jì)算高頻分量的幅值標(biāo)準(zhǔn)差,即σh。

    1.2 集成次數(shù)的確定

    分解過(guò)程中,白噪聲的幅值比值系數(shù)與集成次數(shù)之間存在如下的關(guān)系[6]:

    (1)

    式中:e為輸入信號(hào)與IMFs的相對(duì)誤差;σn為加入的白噪聲幅值標(biāo)準(zhǔn)差;σ0為原始信號(hào)幅值標(biāo)準(zhǔn)差;M為在EEMD方法中集成的次數(shù)。

    當(dāng)加入白噪聲的幅值標(biāo)準(zhǔn)差和輸入信號(hào)與IMFs的相對(duì)誤差e被確定時(shí),便能根據(jù)式(1) 計(jì)算出具體的集成次數(shù)。

    1.3 有效IMF的選取

    利用模糊綜合評(píng)價(jià)法[7]對(duì)所有的IMFs進(jìn)行判別分析。對(duì)于每個(gè)IMF都能計(jì)算出系統(tǒng)的頻率f,阻尼比ξ和振型m,所以選擇這3個(gè)參數(shù)作為評(píng)價(jià)因子。IMF與原始信號(hào)之間的模糊相似系數(shù)ri可由式(2)計(jì)算:

    (2)

    式中: ri為第i個(gè)IMF與原始信號(hào)之間的模糊相似系數(shù);wf,wξ,wm分別為頻率f,阻尼比ξ和振型m的權(quán)重[7]; f,ξ,m分別為原始信號(hào)的頻率、阻尼比和振型;fi,ξi,mi分別為第i個(gè)IMF的頻率、阻尼比和振型。

    考慮到頻率和振型會(huì)隨著計(jì)算階次的升高而逐漸穩(wěn)定,而阻尼一般會(huì)發(fā)生較大波動(dòng),所以文中阻尼比的權(quán)重相對(duì)較小,式(2)中wf=0.5,wξ=0.2,wm=0.3。ri越接近1則表示該第i個(gè)IMF分量與原始信號(hào)的相識(shí)程度越高,當(dāng)ri>0.3時(shí)[7],則認(rèn)為該IMF為有效IMF。

    2 基于自適應(yīng)EEMD分解的盲源分離算法

    信號(hào)分析中的“盲源分離”,即:在沒(méi)有混合系統(tǒng)和原始信號(hào)的先驗(yàn)知識(shí)前提下,從一組未知混合信號(hào)中分離出原始信號(hào)。筆者以廣義特征值盲源分離算法為研究對(duì)象,其具體實(shí)現(xiàn)步驟見(jiàn)文獻(xiàn)[8]。該分離算法雖然算法簡(jiǎn)單且分離速度快,但也有局限性。當(dāng)信號(hào)為非平穩(wěn)混合信號(hào)時(shí),由于源信號(hào)自身的頻譜具有不確定性,且信號(hào)之間存在頻譜疊加現(xiàn)象,這便增加了信號(hào)分離的難度,以至分離結(jié)果不具可靠性。

    針對(duì)這一問(wèn)題,提出了基于自適應(yīng)EEMD分解的盲源分離算法。利用自適應(yīng)EEMD算法處理非平穩(wěn)信號(hào)的優(yōu)勢(shì),先對(duì)混合信號(hào)進(jìn)行分解,得到一系列IMF分量,每個(gè)分量自身的頻率段都不同,且能夠以從高到低的順序進(jìn)行排列,這便很好地避免了“頻譜混疊”的現(xiàn)象。基于自適應(yīng)EEMD分解的盲源分離算法,其實(shí)質(zhì)是先利用自適應(yīng)EEMD在頻域內(nèi)對(duì)混合信號(hào)進(jìn)行一次分離,再利用廣義特征值盲源分離法對(duì)有效的IMF進(jìn)行第二次分離,并對(duì)信號(hào)進(jìn)行重構(gòu)。

    具體的算法步驟如下:

    1)假設(shè)X為混合信號(hào)組成的混合矩陣,通過(guò)自適應(yīng)EEMD分解該信號(hào),得到一系列的有效IMF分量,并構(gòu)造矩陣Cx:

    Cx=[IMFa]

    (3)

    式中:a為有效IMF的具體編號(hào)。

    2)對(duì)矩陣Cx進(jìn)行線性變化,即:

    Cy=CxHT

    (4)

    3)分別計(jì)算Cx和Cy的自相關(guān)矩陣RCx和RCy,并構(gòu)造矩陣束RC:

    (5)

    4)對(duì)矩陣束RC進(jìn)行廣義特征分解,即:

    RC yEC=RCxECDC

    (6)

    式中:EC為特征向量矩陣;DC為特征值矩陣。

    (7)

    上述算法的具體的流程見(jiàn)圖1。

    圖1 自適應(yīng)EEMD分解的盲源分離算法Fig.1 Blind source separation algorithm of adaptive EEMD

    3 協(xié)方差驅(qū)動(dòng)隨機(jī)子空間算法

    隨機(jī)子空間法(SSI)是一種時(shí)域模態(tài)參數(shù)識(shí)別算法[9]。主要被運(yùn)用于線性系統(tǒng),優(yōu)點(diǎn)在于不僅能有效地識(shí)別環(huán)境激勵(lì)下的結(jié)構(gòu)響應(yīng)數(shù)據(jù),同時(shí)其識(shí)別精度高,且不需要事先由響應(yīng)數(shù)據(jù)得到自由衰減曲線。鑒于此,該方法越來(lái)越得到人們的重視。隨機(jī)子空間法主要有兩種,即協(xié)方差驅(qū)動(dòng)隨機(jī)子空間算法和數(shù)據(jù)驅(qū)動(dòng)隨機(jī)子空間算法。筆者以協(xié)方差驅(qū)動(dòng)隨機(jī)子空間為研究對(duì)象,就其算法步驟進(jìn)行簡(jiǎn)單分析,詳細(xì)分析可見(jiàn)文獻(xiàn)[9]。

    1)利用結(jié)構(gòu)響應(yīng)數(shù)據(jù)構(gòu)造Toeplitz矩陣:

    (8)

    2)對(duì)T1|i進(jìn)行奇異值分解(SVD):

    (9)

    3)利用U1,S1求得系統(tǒng)擴(kuò)展可觀矩陣Ti和擴(kuò)展可控矩陣Δi:

    (10)

    4)計(jì)算系統(tǒng)的狀態(tài)轉(zhuǎn)換矩陣A和輸出矩陣C:

    (11)

    5)利用系統(tǒng)的狀態(tài)轉(zhuǎn)換矩陣A和輸出矩陣C計(jì)算出系統(tǒng)的模態(tài)參數(shù)。

    基于筆者所提算法的協(xié)方差驅(qū)動(dòng)隨機(jī)子空間模態(tài)參數(shù)識(shí)別算法的具體流程見(jiàn)圖2。

    圖2 模態(tài)參數(shù)識(shí)別算法流程Fig.2 Flowchart of modal parameters identification algorithm

    4 測(cè)試與驗(yàn)證

    利用自適應(yīng)EEMD分解橋梁結(jié)構(gòu)測(cè)試信號(hào)后,再利用模糊綜合評(píng)價(jià)法選出有效的IMF分量,最后對(duì)有效的IMF分量矩陣進(jìn)行盲辨識(shí),實(shí)現(xiàn)對(duì)橋梁結(jié)構(gòu)動(dòng)力測(cè)試信號(hào)的處理。分別用模擬信號(hào)和實(shí)際橋梁健康監(jiān)測(cè)系統(tǒng)的動(dòng)力測(cè)試信號(hào)對(duì)筆者所提方法進(jìn)行驗(yàn)證,以說(shuō)明其可行性。

    4.1 模擬信號(hào)

    模擬信號(hào)由1,3,5 Hz的3個(gè)正弦信號(hào)疊加噪聲水平約為10%的隨機(jī)噪聲組成:

    首先利用自適應(yīng)EEMD分解算法對(duì)原始信號(hào)進(jìn)行分解,結(jié)果見(jiàn)圖3。圖中第1行代表原始信號(hào)的時(shí)域,圖3中顯示共有10個(gè)IMF分量,可知這其中一定存在虛假模態(tài)分量。為了剔除虛假的模態(tài)分量,利用模糊綜合評(píng)價(jià)法計(jì)算出每個(gè)IMF分量與原始信號(hào)之間的模糊相似系數(shù)。同時(shí)為了進(jìn)一步驗(yàn)證模糊相似系數(shù)能很好地選出有效的IMF分量,也計(jì)算了每個(gè)IMF分量與3個(gè)仿真信號(hào)之間的相似系數(shù),具體數(shù)據(jù)結(jié)果見(jiàn)表1。

    表1 IMF分量與仿真信號(hào)以及原始信號(hào)之間的模糊相似系數(shù)

    圖3 自適應(yīng)EEMD分解原始信號(hào)Fig.3 Original signal of adaptive EEMD

    表1中IMF2,IMF3以及IMF5是有效IMF分量(系數(shù)大于0.3)。根據(jù)IMF分量與每個(gè)仿真信號(hào)之間的相似系數(shù)可知,與5 Hz正弦信號(hào)最為相似的是IMF2,與3 Hz正弦信號(hào)最為相似的是IMF3,與1 Hz正弦信號(hào)最為相似的是IMF5??梢?jiàn)通過(guò)相似系數(shù)和模糊相似系數(shù)都能分辨出有效IMF分量。

    將選取出來(lái)的有效IMF分量矩陣作為盲源分離算法的輸入并計(jì)算,可得圖4、圖5所示結(jié)果。

    圖4 盲源分離結(jié)果Fig.4 Results of blind source separation

    由圖4可知,有效IMF分量經(jīng)過(guò)盲源分離之后,每個(gè)有效IMF分量與其對(duì)應(yīng)的原仿真信號(hào)之間很接近。

    圖5 分解結(jié)果對(duì)比Fig.5 Comparison diagram of decomposition results

    由圖5可知,重構(gòu)信號(hào)與原始仿真信號(hào)很接近。利用相關(guān)性分析[10],計(jì)算得這兩者之間的相似系數(shù)為0.976。相比原始信號(hào)(含噪聲),重構(gòu)信號(hào)中只含很少部分的噪聲。

    通過(guò)以上分析可知:可以將IMF分量與原始信號(hào)之間的模糊相似系數(shù)作為選取有效IMF分量的依據(jù);該算法能對(duì)原始信號(hào)進(jìn)行有效地降噪;對(duì)選取出來(lái)的有效IMF分量矩陣進(jìn)行盲源分離可以實(shí)現(xiàn)仿真信號(hào)的精確提取,提取得到的單獨(dú)信號(hào)與模擬信號(hào)中加入的正弦信號(hào)能一一對(duì)應(yīng),且相似系數(shù)都在0.95以上。

    4.2 斜拉橋?qū)崪y(cè)動(dòng)力信號(hào)

    筆者以長(zhǎng)江上某雙塔雙索面斜拉橋?yàn)樽R(shí)別對(duì)象,該斜拉橋全長(zhǎng)2 088 m,其主跨為1 088 m,具體孔跨布置見(jiàn)圖6。橋上共布置豎向加速度傳感器14個(gè),位于主梁主跨1/6截面和次邊跨1/2截面的上游和下游處,具體位置見(jiàn)圖6。加速度信號(hào)采樣頻率為20 Hz,測(cè)試時(shí)間為48 h??紤]到實(shí)際結(jié)構(gòu)都處于環(huán)境激勵(lì)下,所以傳感器采集得到的信號(hào)中會(huì)含有一定的噪聲。為了更好地識(shí)別結(jié)構(gòu)的模態(tài)參數(shù)則需要事先對(duì)信號(hào)進(jìn)行預(yù)處理。利用本文算法對(duì)原始信號(hào)進(jìn)行預(yù)處理,得到重構(gòu)信號(hào),如圖7。為了檢驗(yàn)筆者所提方法的可行性,分別將經(jīng)過(guò)筆者提出方法處理的信號(hào)和未經(jīng)處理的信號(hào)作為協(xié)方差驅(qū)動(dòng)隨機(jī)子空間的輸入,通過(guò)對(duì)比兩種不同輸入情況下的Hilbert-Huang譜和穩(wěn)定圖來(lái)判別筆者提出方法是否可行。

    圖6 傳感器的布置Fig.6 Sensor layout in the SHM

    圖7 實(shí)測(cè)信號(hào)與重構(gòu)信號(hào)對(duì)比Fig.7 Comparison diagram of the measured signal and the reconstructed signal

    由圖7可知,筆者所提方法能夠有效地去除實(shí)測(cè)信號(hào)中的噪聲,并能對(duì)其進(jìn)行重構(gòu)。為了檢驗(yàn)筆者所提算法的可行性,分別得到了未經(jīng)預(yù)處理和經(jīng)過(guò)預(yù)處理信號(hào)的Hilbert-Huang譜[11],如圖8。

    圖8 Hilbert-Huang譜Fig.8 Hilbert-Huang spectrums

    由圖8可知,經(jīng)過(guò)預(yù)處理之后的信號(hào)得到的Hilbert-Huang譜中的瞬時(shí)頻率更為連續(xù)和清晰。為了更進(jìn)一步驗(yàn)證所提算法能很好地消除原始信號(hào)中的噪聲,并保留結(jié)構(gòu)自身的信息,利用協(xié)方差驅(qū)動(dòng)隨機(jī)子空間算法對(duì)信號(hào)進(jìn)行識(shí)別,得到的穩(wěn)定圖和前3階振型,如圖9、圖10。

    圖9 信號(hào)的穩(wěn)定圖Fig.9 Signal stability diagram

    圖10 該斜拉橋前3階振型Fig.10 First three orders vibration chart of the cable-stayed bridge

    由圖9可見(jiàn),用筆者所提算法處理后的信號(hào)能識(shí)別到更多的頻率值;預(yù)處理之后能得到更為穩(wěn)定的頻率值,且穩(wěn)定軸更為清晰;預(yù)處理之后的信號(hào)含有的結(jié)構(gòu)信號(hào)更為豐富。

    圖10是協(xié)方差驅(qū)動(dòng)隨機(jī)子空間識(shí)別預(yù)處理過(guò)的實(shí)測(cè)信號(hào)而得到的該斜拉橋前3階模態(tài)振型。由圖10可知前3階模態(tài)振型圖與實(shí)際振型圖很相似,相似度在95%左右,進(jìn)一步驗(yàn)證筆者所提算法的可行性。

    5 結(jié) 論

    經(jīng)模擬信號(hào)與實(shí)橋數(shù)據(jù)的驗(yàn)證,可得如下結(jié)論:

    1)筆者所提的自適應(yīng)EEMD分解能根據(jù)信號(hào)自身的特點(diǎn),確定加入白噪聲的幅值標(biāo)準(zhǔn)差、EEMD集成次數(shù)等,以實(shí)現(xiàn)對(duì)信號(hào)的自適應(yīng)分解。

    2)通過(guò)對(duì)比計(jì)算每個(gè)IMF與仿真信號(hào)的相似系數(shù)以及用模糊綜合評(píng)價(jià)法計(jì)算得到的結(jié)果,可知能利用模糊綜合評(píng)價(jià)法實(shí)現(xiàn)對(duì)有效IMF的自動(dòng)選取。

    3)利用盲源分離算法能實(shí)現(xiàn)對(duì)有效IMF分量的二次分離,并能實(shí)現(xiàn)對(duì)信號(hào)的重構(gòu)。

    4)通過(guò)對(duì)比Hilbert-Huang譜和穩(wěn)定圖可知,筆者所提算法能對(duì)橋梁結(jié)構(gòu)的動(dòng)力測(cè)試信號(hào)進(jìn)行有效的分解和降噪,且提取的結(jié)構(gòu)信息更為豐富準(zhǔn)確。

    5)將筆者所提算法運(yùn)用于實(shí)橋測(cè)試數(shù)據(jù),計(jì)算結(jié)果表明筆者所提的基于自適應(yīng)EEMD的盲源分離算法能用于實(shí)際橋梁的動(dòng)力測(cè)試分析中。

    [1] 漆景星,侯艷紅. 橋梁可靠性評(píng)估綜述[J]. 浙江交通職業(yè)技術(shù)學(xué)院學(xué)報(bào),2011,12(1):20-22. QI Jingxing, HOU Yanhong. Bridge reliability evaluation review [J].JournalofZhejiangInstituteofCommunications,2011,12(1):20-22.

    [2] LIN Ping, ZHANG Nanxiong, NI Bin. On-line modal parameter monitoring of bridges exploiting multi-core capacity by recursive stochastic subspace identification method[C]∥2008AmericanControlConference. Washington,D.C., USA:[s.n.],2008.

    [3] 鄭近德,程軍圣,楊宇.改進(jìn)的EEMD算法及其應(yīng)用研究[J].振動(dòng)與沖擊,2013,32(21):21-26. ZHENG Jingde, CHENG Junsheng, YANG Yu. Modified EEMD algorithm and its applications[J].JournalofVibrationandShock,2013,32(21):21-26.

    [4] 高穎.盲源分離理論及其在地球物理勘探中的應(yīng)用[D].長(zhǎng)春:吉林大學(xué),2008. GAO Ying.BlindSourceSeparationTheoryandItsApplicationinGeophysicalExploration[D]. Changchun: Jilin University,2008.

    [5] LARDIES J , MINH-NGI T. Modal parameter identification of stay cables from output-only measurements[J].MechanicalSystemsandSignalProcessing,2011,25(1):133-150.

    [6] 蔡艷平,李艾華,徐斌,等.集成經(jīng)驗(yàn)?zāi)B(tài)分解中加入白噪聲的自適應(yīng)準(zhǔn)則[J].振動(dòng)、測(cè)試與診斷,2011,31(6):709-714. CAI Yanping, LI Aihua, XU Bin, et al. The adaptive rule of ensemble empirical mode decomposition with added white noise [J].JournalofVibration,Measurement&Diagnosis,2011,31(6):709-714.

    [7] 許順國(guó),牟瑞芳,張雪梅. 模糊數(shù)學(xué)綜合評(píng)判法在水質(zhì)評(píng)價(jià)中的應(yīng)用——以成都市府河為例[J]. 唐山師范學(xué)院學(xué)報(bào),2007,29(2):68-70. XU Shunguo, MOU Ruifang, ZHANG Xuemei. The application of fuzzy mathmatical comprehensive judgment to the evaluation of water quality: about the Fu river in Chengdu city[J].JournalofTangshanTeachersCollege,2007,29(2):68-70.

    [8] 黃青華.基于源信號(hào)模型的盲分離技術(shù)研究及應(yīng)用[D].上海:上海交通大學(xué),2007. HUANG Qinghua.ResearchandApplicationofBlindSeparationTechnologybasedonSourceSignalModel[D]. Shanghai: Shanghai Jiaotong University,2007.

    [9] PEETERS B, ROECK G D. Reference-based stochastic subspace identification for output-only modal analysis[J].MechanicalSystemsandSignalProcessing,1999,13(6):855-878.

    [10] 張渝.基于DEA的典型相關(guān)分析在商業(yè)銀行信用風(fēng)險(xiǎn)評(píng)估中的應(yīng)用[J].數(shù)學(xué)的實(shí)踐與認(rèn)識(shí),2007,37(24):33-39. ZHANG Yu. Use of CCA in the DEA context for credit risk assessment in commercial banks [J].MathematicsinPracticeandTheory,2007, 37(24): 33-39.

    [11] HANG Jing,YUAN Haiqing, ZHAO Yi, et al. Applying Hilbert-Huang transform to identifying structural modal parameters[C]//EducationTechnologyandTraining&GeoscienceandRemoteSensing.[S.l.]:IEEE,2008:617-621.

    Modal Parameter Identification of Bridge Structure Based on Adaptive EEMD and Blind Identification Algorithm

    CHEN Yonggao

    (Zhejiang Industry Polytechnic College, Shaoxing 312000, Zhejiang, P.R.China)

    Due to the defects of ensemble empirical mode decomposition (EEMD) in signal pretreatment, a blind source separation algorithm based on adaptive EEMD decomposition was proposed, that was: firstly, the amplitude standard deviation and integration times of the added white noise were confirmed according to the characteristics of the original signal, and then a EEMD decomposition was carried out; secondly, a fuzzy comprehensive evaluation on the obtained IMF components was carried out to select out the effective components of the IMF, and then IMF component matrix was established; finally, blind source separation algorithm was used for blind identification, and the signal decomposition and reconstruction was completed. The proposed algorithm was verified by analog signal and measuring vibration signal of bridges respectively. The results show that the proposed algorithm is feasible and can be applied to the signal preprocessing of actual structure.

    bridge engineering; adaptive EEMD; blind identification; fuzzy comprehensive evaluation; parameter identification

    10.3969/j.issn.1674-0696.2016.03.03

    2015-06-17;

    2015-10-09

    浙江省教育廳科研項(xiàng)目(Y201432555);浙江省住建廳科研項(xiàng)目(2014ZI26);紹興市科技計(jì)劃項(xiàng)目(2014B70003)

    陳永高(1984—),男,江蘇鹽城人,工程師,主要從事土木工程建造與管理方面的研究。E-mail:higaoge@163.com。

    U446.3

    A

    1674-0696(2016)03-011-06

    猜你喜歡
    盲源幅值預(yù)處理
    基于干擾重構(gòu)和盲源分離的混合極化抗SMSP干擾
    基于預(yù)處理MUSIC算法的分布式陣列DOA估計(jì)
    基于S變換的交流電網(wǎng)幅值檢測(cè)系統(tǒng)計(jì)算機(jī)仿真研究
    電子制作(2017年7期)2017-06-05 09:36:13
    淺談PLC在預(yù)處理生產(chǎn)線自動(dòng)化改造中的應(yīng)用
    正序電壓幅值檢測(cè)及諧波抑制的改進(jìn)
    絡(luò)合萃取法預(yù)處理H酸廢水
    低壓電力線信道脈沖噪聲的幅值與寬度特征
    改進(jìn)的互信息最小化非線性盲源分離算法
    基于零序電壓幅值增量的消弧線圈調(diào)諧新方法
    盲源分離技術(shù)在化探數(shù)據(jù)處理中的思考與探索
    久久精品91无色码中文字幕| 香蕉久久夜色| av电影中文网址| www日本在线高清视频| 欧美日韩精品网址| 亚洲自拍偷在线| 丁香欧美五月| 久久久国产精品麻豆| 大型av网站在线播放| 真人一进一出gif抽搐免费| 欧美中文日本在线观看视频| 中文字幕高清在线视频| 中文亚洲av片在线观看爽| 亚洲精品一区av在线观看| 午夜a级毛片| 久久久久久国产a免费观看| 欧美成人性av电影在线观看| 女警被强在线播放| 法律面前人人平等表现在哪些方面| 免费在线观看日本一区| 中文亚洲av片在线观看爽| 亚洲精品一区av在线观看| 天堂√8在线中文| 国产精华一区二区三区| 久久久久精品国产欧美久久久| 一边摸一边抽搐一进一出视频| 无遮挡黄片免费观看| 黑丝袜美女国产一区| 久久性视频一级片| 两个人看的免费小视频| 身体一侧抽搐| 成人免费观看视频高清| 午夜精品久久久久久毛片777| 精品第一国产精品| 亚洲三区欧美一区| 97人妻天天添夜夜摸| 最近最新中文字幕大全电影3 | 亚洲三区欧美一区| 午夜a级毛片| 成人亚洲精品一区在线观看| 18禁裸乳无遮挡免费网站照片 | 在线视频色国产色| 看免费av毛片| 成熟少妇高潮喷水视频| 最近最新中文字幕大全电影3 | 大型av网站在线播放| 在线观看日韩欧美| 国产精品一区二区免费欧美| 搞女人的毛片| 国产精品影院久久| 久久久久久久午夜电影| 午夜精品在线福利| 天堂影院成人在线观看| 9热在线视频观看99| 一边摸一边抽搐一进一出视频| 成人国产一区最新在线观看| 欧美日本中文国产一区发布| 麻豆国产av国片精品| 国产精品 欧美亚洲| 亚洲国产欧美网| 夜夜看夜夜爽夜夜摸| 在线视频色国产色| 亚洲九九香蕉| x7x7x7水蜜桃| 91精品国产国语对白视频| 97人妻精品一区二区三区麻豆 | 真人做人爱边吃奶动态| 色哟哟哟哟哟哟| 亚洲片人在线观看| 亚洲av熟女| 亚洲欧美激情综合另类| 麻豆一二三区av精品| 国产亚洲精品久久久久久毛片| 国产亚洲精品综合一区在线观看 | 女警被强在线播放| 欧美一级a爱片免费观看看 | 老司机午夜十八禁免费视频| 窝窝影院91人妻| 一边摸一边抽搐一进一出视频| 日韩欧美国产一区二区入口| 美女高潮喷水抽搐中文字幕| 久久亚洲真实| 色av中文字幕| 天天一区二区日本电影三级 | 大型av网站在线播放| 亚洲人成伊人成综合网2020| 丰满人妻熟妇乱又伦精品不卡| 黄片播放在线免费| 不卡一级毛片| 欧美日本中文国产一区发布| 亚洲一卡2卡3卡4卡5卡精品中文| 久久青草综合色| 97人妻天天添夜夜摸| 91成人精品电影| 19禁男女啪啪无遮挡网站| 中国美女看黄片| 成人三级做爰电影| 日韩国内少妇激情av| 99国产极品粉嫩在线观看| 天堂√8在线中文| 久久久久九九精品影院| 欧美中文综合在线视频| 亚洲欧美精品综合一区二区三区| 国产一区在线观看成人免费| 久久精品国产亚洲av高清一级| 国产成人av教育| 高潮久久久久久久久久久不卡| 精品福利观看| 老司机靠b影院| 久久精品国产清高在天天线| 夜夜爽天天搞| 性少妇av在线| 亚洲av片天天在线观看| 国产亚洲精品第一综合不卡| 99国产综合亚洲精品| 国产av又大| 色尼玛亚洲综合影院| 美国免费a级毛片| 嫩草影视91久久| 亚洲精品美女久久久久99蜜臀| a级毛片在线看网站| 久久草成人影院| 香蕉久久夜色| av在线播放免费不卡| 少妇裸体淫交视频免费看高清 | 中文亚洲av片在线观看爽| 在线观看免费视频网站a站| 亚洲av熟女| 久久久久国产一级毛片高清牌| 中文字幕高清在线视频| 亚洲国产毛片av蜜桃av| 午夜免费激情av| 高潮久久久久久久久久久不卡| 不卡一级毛片| 欧美日韩中文字幕国产精品一区二区三区 | 亚洲成人精品中文字幕电影| 久久久精品欧美日韩精品| 日韩欧美国产在线观看| 看片在线看免费视频| 两性午夜刺激爽爽歪歪视频在线观看 | 日韩一卡2卡3卡4卡2021年| 亚洲欧美激情综合另类| 中出人妻视频一区二区| 欧洲精品卡2卡3卡4卡5卡区| 少妇的丰满在线观看| 国内久久婷婷六月综合欲色啪| 久久国产精品人妻蜜桃| 窝窝影院91人妻| 欧美一级a爱片免费观看看 | 久久天躁狠狠躁夜夜2o2o| 一个人观看的视频www高清免费观看 | 女同久久另类99精品国产91| 丝袜在线中文字幕| 亚洲精品一区av在线观看| 国产av一区二区精品久久| 久久天堂一区二区三区四区| 十八禁网站免费在线| 国产精品自产拍在线观看55亚洲| 国内毛片毛片毛片毛片毛片| 大香蕉久久成人网| 成人三级做爰电影| 国产亚洲精品综合一区在线观看 | 一进一出抽搐动态| 午夜免费鲁丝| 免费在线观看影片大全网站| 啪啪无遮挡十八禁网站| 人妻丰满熟妇av一区二区三区| 黑人巨大精品欧美一区二区蜜桃| 黄色视频不卡| 视频在线观看一区二区三区| 亚洲一卡2卡3卡4卡5卡精品中文| 最新在线观看一区二区三区| 法律面前人人平等表现在哪些方面| 国内久久婷婷六月综合欲色啪| 中文字幕av电影在线播放| 亚洲全国av大片| 亚洲美女黄片视频| 啦啦啦 在线观看视频| 中亚洲国语对白在线视频| 国产单亲对白刺激| 黄色毛片三级朝国网站| 亚洲成人免费电影在线观看| 国产成年人精品一区二区| 国产精品美女特级片免费视频播放器 | 伦理电影免费视频| 村上凉子中文字幕在线| 成年人黄色毛片网站| 99久久99久久久精品蜜桃| 久久久久亚洲av毛片大全| 国产精品久久久久久人妻精品电影| 怎么达到女性高潮| 亚洲国产精品999在线| 国产成人影院久久av| 国产成人精品久久二区二区免费| 亚洲男人天堂网一区| 97超级碰碰碰精品色视频在线观看| av电影中文网址| 久久人人97超碰香蕉20202| 国产一卡二卡三卡精品| 88av欧美| 亚洲精品中文字幕一二三四区| 久久久久九九精品影院| 成年女人毛片免费观看观看9| 日日夜夜操网爽| 一级毛片高清免费大全| 曰老女人黄片| 国产精品久久久久久亚洲av鲁大| av天堂在线播放| 亚洲,欧美精品.| 亚洲黑人精品在线| 最新在线观看一区二区三区| 中文字幕av电影在线播放| 亚洲精品国产区一区二| 日韩欧美三级三区| cao死你这个sao货| 欧美一级毛片孕妇| 国产激情久久老熟女| 久久青草综合色| 久久国产亚洲av麻豆专区| 亚洲国产欧美网| 一边摸一边做爽爽视频免费| 国产午夜福利久久久久久| 97人妻精品一区二区三区麻豆 | 欧美日本亚洲视频在线播放| 在线天堂中文资源库| 69精品国产乱码久久久| 最近最新中文字幕大全电影3 | 韩国av一区二区三区四区| 波多野结衣高清无吗| 欧美日本中文国产一区发布| 久久九九热精品免费| 久久午夜综合久久蜜桃| 国产精品亚洲一级av第二区| 亚洲激情在线av| 亚洲aⅴ乱码一区二区在线播放 | 欧美+亚洲+日韩+国产| 精品电影一区二区在线| av天堂在线播放| 亚洲av美国av| 又大又爽又粗| 欧美国产日韩亚洲一区| 美女大奶头视频| 青草久久国产| 丝袜美足系列| 日本 av在线| 免费在线观看亚洲国产| 欧美亚洲日本最大视频资源| www.www免费av| 中文字幕精品免费在线观看视频| 成人亚洲精品一区在线观看| 啦啦啦韩国在线观看视频| 久久婷婷人人爽人人干人人爱 | 国产精品香港三级国产av潘金莲| x7x7x7水蜜桃| 大型黄色视频在线免费观看| 视频在线观看一区二区三区| 成年版毛片免费区| 少妇裸体淫交视频免费看高清 | 无人区码免费观看不卡| 久久婷婷成人综合色麻豆| 老鸭窝网址在线观看| 在线国产一区二区在线| 日本精品一区二区三区蜜桃| 国产精品一区二区在线不卡| 国产av一区二区精品久久| 国产97色在线日韩免费| 97碰自拍视频| 久久九九热精品免费| 亚洲va日本ⅴa欧美va伊人久久| 一二三四社区在线视频社区8| 亚洲av电影不卡..在线观看| 在线视频色国产色| 欧美色欧美亚洲另类二区 | e午夜精品久久久久久久| 久久久久久免费高清国产稀缺| 日日干狠狠操夜夜爽| 看黄色毛片网站| 一区福利在线观看| av网站免费在线观看视频| 极品人妻少妇av视频| 两个人看的免费小视频| 午夜福利影视在线免费观看| 男女午夜视频在线观看| 中文字幕另类日韩欧美亚洲嫩草| 成年女人毛片免费观看观看9| 日韩精品中文字幕看吧| 精品人妻在线不人妻| 黄片大片在线免费观看| 九色国产91popny在线| 精品卡一卡二卡四卡免费| 色哟哟哟哟哟哟| 国产精品一区二区三区四区久久 | 国产精品免费一区二区三区在线| 在线观看午夜福利视频| 亚洲欧美日韩高清在线视频| 精品久久久久久,| 啪啪无遮挡十八禁网站| 午夜久久久久精精品| 亚洲国产欧美日韩在线播放| 女生性感内裤真人,穿戴方法视频| 男人操女人黄网站| 精品一区二区三区视频在线观看免费| 一区福利在线观看| 国产欧美日韩一区二区精品| 又黄又粗又硬又大视频| 在线观看免费视频网站a站| 亚洲五月婷婷丁香| 免费一级毛片在线播放高清视频 | 老鸭窝网址在线观看| 99香蕉大伊视频| 九色国产91popny在线| 又紧又爽又黄一区二区| 国产精品一区二区三区四区久久 | 午夜福利一区二区在线看| 国产成人精品久久二区二区91| 狂野欧美激情性xxxx| 午夜激情av网站| 99精品欧美一区二区三区四区| 18禁黄网站禁片午夜丰满| 又黄又粗又硬又大视频| 午夜福利一区二区在线看| 久9热在线精品视频| 精品国产国语对白av| 在线视频色国产色| 麻豆av在线久日| 国产精品野战在线观看| 97人妻精品一区二区三区麻豆 | 免费高清在线观看日韩| 中文字幕色久视频| av免费在线观看网站| 精品一品国产午夜福利视频| 一个人观看的视频www高清免费观看 | 黄片小视频在线播放| 一级作爱视频免费观看| 亚洲自偷自拍图片 自拍| 91麻豆精品激情在线观看国产| 欧美日韩一级在线毛片| 亚洲最大成人中文| 欧美黄色片欧美黄色片| 三级毛片av免费| 777久久人妻少妇嫩草av网站| 免费高清视频大片| 欧美一级毛片孕妇| 麻豆久久精品国产亚洲av| 久久人妻熟女aⅴ| 午夜免费鲁丝| 国产又色又爽无遮挡免费看| 别揉我奶头~嗯~啊~动态视频| 亚洲自偷自拍图片 自拍| 亚洲性夜色夜夜综合| 99精品在免费线老司机午夜| 女生性感内裤真人,穿戴方法视频| 日韩大尺度精品在线看网址 | 亚洲午夜理论影院| 久久人妻福利社区极品人妻图片| 亚洲午夜理论影院| 国产伦人伦偷精品视频| 国产精品国产高清国产av| 日本免费a在线| 国产亚洲精品一区二区www| 黄色视频,在线免费观看| 亚洲av片天天在线观看| 国内毛片毛片毛片毛片毛片| 在线十欧美十亚洲十日本专区| 18禁美女被吸乳视频| 欧美性长视频在线观看| av网站免费在线观看视频| 成年版毛片免费区| 日韩大码丰满熟妇| 妹子高潮喷水视频| 免费在线观看黄色视频的| 亚洲专区中文字幕在线| 免费不卡黄色视频| 搡老岳熟女国产| 欧美成人免费av一区二区三区| 欧美在线黄色| 国产av一区在线观看免费| 一进一出好大好爽视频| 50天的宝宝边吃奶边哭怎么回事| 国产精品九九99| 精品久久久久久久人妻蜜臀av | 日韩欧美国产在线观看| 久久久国产欧美日韩av| a在线观看视频网站| 国产免费男女视频| 十八禁网站免费在线| 欧美另类亚洲清纯唯美| 欧美亚洲日本最大视频资源| 51午夜福利影视在线观看| 精品第一国产精品| 丁香欧美五月| 在线观看日韩欧美| 麻豆一二三区av精品| 国产精品久久久久久精品电影 | 91麻豆av在线| 亚洲成人免费电影在线观看| 午夜亚洲福利在线播放| 丁香六月欧美| 丰满的人妻完整版| 久久精品国产99精品国产亚洲性色 | 999久久久国产精品视频| 中国美女看黄片| www日本在线高清视频| 91成年电影在线观看| 夜夜躁狠狠躁天天躁| 一区二区日韩欧美中文字幕| 99香蕉大伊视频| 亚洲人成网站在线播放欧美日韩| 99国产综合亚洲精品| 亚洲第一电影网av| 伊人久久大香线蕉亚洲五| av网站免费在线观看视频| 久久青草综合色| 国产xxxxx性猛交| 久久久精品国产亚洲av高清涩受| 亚洲人成网站在线播放欧美日韩| 欧美日韩乱码在线| 欧美最黄视频在线播放免费| 丰满的人妻完整版| 欧美一级毛片孕妇| av天堂在线播放| 一级作爱视频免费观看| 国产精品亚洲一级av第二区| 久久精品国产亚洲av高清一级| 亚洲av片天天在线观看| 十八禁网站免费在线| 丰满人妻熟妇乱又伦精品不卡| 成人三级黄色视频| 两人在一起打扑克的视频| 悠悠久久av| av网站免费在线观看视频| 在线观看66精品国产| 午夜精品在线福利| 国产av在哪里看| 欧美成人免费av一区二区三区| av视频免费观看在线观看| 黄片播放在线免费| 欧美日韩中文字幕国产精品一区二区三区 | 一区二区三区高清视频在线| 亚洲人成伊人成综合网2020| av在线播放免费不卡| 国产精品一区二区在线不卡| av有码第一页| 少妇熟女aⅴ在线视频| 亚洲国产欧美一区二区综合| 亚洲欧美激情在线| 国产精品 国内视频| 一区二区三区国产精品乱码| 国产私拍福利视频在线观看| 欧美日韩亚洲国产一区二区在线观看| 日韩免费av在线播放| 禁无遮挡网站| 国产精品1区2区在线观看.| 亚洲七黄色美女视频| 男男h啪啪无遮挡| 久久久水蜜桃国产精品网| 亚洲男人的天堂狠狠| 国产精品久久久av美女十八| 老司机午夜十八禁免费视频| 少妇粗大呻吟视频| 男人的好看免费观看在线视频 | 日本撒尿小便嘘嘘汇集6| 日本免费一区二区三区高清不卡 | 身体一侧抽搐| 国产精华一区二区三区| 国产亚洲精品久久久久久毛片| a级毛片在线看网站| 欧美精品亚洲一区二区| 中亚洲国语对白在线视频| 欧美亚洲日本最大视频资源| 首页视频小说图片口味搜索| 午夜免费观看网址| 色播在线永久视频| 啪啪无遮挡十八禁网站| 丝袜在线中文字幕| 可以免费在线观看a视频的电影网站| 国产蜜桃级精品一区二区三区| 脱女人内裤的视频| 黄频高清免费视频| 亚洲中文字幕日韩| 国内精品久久久久久久电影| 精品人妻在线不人妻| 免费在线观看日本一区| 精品免费久久久久久久清纯| 激情视频va一区二区三区| 非洲黑人性xxxx精品又粗又长| 亚洲国产精品999在线| 中文字幕最新亚洲高清| 亚洲av电影在线进入| 欧美 亚洲 国产 日韩一| 18禁裸乳无遮挡免费网站照片 | 在线观看一区二区三区| 免费观看精品视频网站| 免费女性裸体啪啪无遮挡网站| 国产视频一区二区在线看| 国产欧美日韩综合在线一区二区| 国产一区在线观看成人免费| 中文字幕久久专区| 久久精品国产亚洲av高清一级| 亚洲 欧美 日韩 在线 免费| 夜夜夜夜夜久久久久| 精品国产美女av久久久久小说| 久久 成人 亚洲| 亚洲一码二码三码区别大吗| 老司机靠b影院| 久久草成人影院| 美女高潮喷水抽搐中文字幕| 成人永久免费在线观看视频| 亚洲人成电影观看| 女人高潮潮喷娇喘18禁视频| 制服人妻中文乱码| 亚洲av电影在线进入| 一二三四在线观看免费中文在| 可以在线观看的亚洲视频| 美女免费视频网站| 啦啦啦观看免费观看视频高清 | 黑人巨大精品欧美一区二区mp4| 国产麻豆成人av免费视频| 精品一品国产午夜福利视频| 日日爽夜夜爽网站| 日韩av在线大香蕉| 99久久久亚洲精品蜜臀av| 老司机午夜福利在线观看视频| 88av欧美| 在线十欧美十亚洲十日本专区| 制服诱惑二区| 中文字幕人成人乱码亚洲影| 啦啦啦韩国在线观看视频| av视频在线观看入口| 色哟哟哟哟哟哟| 婷婷丁香在线五月| 欧美日本中文国产一区发布| 亚洲午夜精品一区,二区,三区| av福利片在线| 精品人妻在线不人妻| 亚洲中文日韩欧美视频| 18美女黄网站色大片免费观看| 嫩草影视91久久| 国产av一区在线观看免费| or卡值多少钱| 亚洲中文字幕日韩| 可以免费在线观看a视频的电影网站| 国产精品二区激情视频| www日本在线高清视频| 国产熟女午夜一区二区三区| 国产又爽黄色视频| 中文字幕av电影在线播放| 久久久久久久久久久久大奶| 日韩欧美国产一区二区入口| 午夜久久久久精精品| 久久中文字幕人妻熟女| 亚洲 欧美 日韩 在线 免费| 亚洲第一欧美日韩一区二区三区| 天堂√8在线中文| 老熟妇乱子伦视频在线观看| 国产成人av教育| 久久天堂一区二区三区四区| 久久精品亚洲精品国产色婷小说| 久久午夜亚洲精品久久| 人人妻人人澡欧美一区二区 | 黄色视频,在线免费观看| 日韩中文字幕欧美一区二区| 亚洲第一欧美日韩一区二区三区| 变态另类丝袜制服| 九色亚洲精品在线播放| 亚洲成a人片在线一区二区| 深夜精品福利| 成人国产综合亚洲| 精品一区二区三区四区五区乱码| 亚洲第一青青草原| 欧美国产精品va在线观看不卡| 精品熟女少妇八av免费久了| 国产蜜桃级精品一区二区三区| 咕卡用的链子| 久久九九热精品免费| 精品国产一区二区三区四区第35| 成人精品一区二区免费| 一本久久中文字幕| 涩涩av久久男人的天堂| 国产在线精品亚洲第一网站| 亚洲一码二码三码区别大吗| 日韩精品青青久久久久久| 久久国产精品人妻蜜桃| 亚洲av熟女| 波多野结衣高清无吗| 中文字幕最新亚洲高清| 香蕉丝袜av| 一进一出抽搐gif免费好疼| 国产激情久久老熟女| 日韩视频一区二区在线观看| 激情视频va一区二区三区| 欧美日韩亚洲国产一区二区在线观看| 国产成人精品在线电影| 国产午夜福利久久久久久| 成年女人毛片免费观看观看9| 国产欧美日韩精品亚洲av| 波多野结衣巨乳人妻| 这个男人来自地球电影免费观看| 亚洲国产精品sss在线观看| 可以在线观看毛片的网站| 日韩高清综合在线| 国产真人三级小视频在线观看| 热re99久久国产66热| 高清毛片免费观看视频网站| 97人妻精品一区二区三区麻豆 | 老熟妇乱子伦视频在线观看| aaaaa片日本免费| 老汉色av国产亚洲站长工具| 变态另类成人亚洲欧美熟女 | 88av欧美| 97人妻天天添夜夜摸| 在线av久久热| 久久久久久久午夜电影| 在线观看免费视频日本深夜| 日本vs欧美在线观看视频|