一種強(qiáng)背景噪聲環(huán)境下的模態(tài)參數(shù)識別方法
【專利摘要】一種強(qiáng)背景噪聲環(huán)境下的模態(tài)參數(shù)識別方法，先通過力錘敲擊實(shí)驗(yàn)測得一級脈沖響應(yīng)信號，然后通過譜減法對一級脈沖響應(yīng)信號進(jìn)行初步降噪，得到二級脈沖響應(yīng)信號，再使用最小均方差短時(shí)譜估計(jì)方法對二級脈沖響應(yīng)信號進(jìn)行二次降噪，得到理想的脈沖響應(yīng)信號，最后采用模態(tài)參數(shù)識別算法對理想的脈沖響應(yīng)信進(jìn)行模態(tài)參數(shù)識別，本發(fā)明具有自適應(yīng)、計(jì)算速度快，強(qiáng)噪聲環(huán)境下模態(tài)參數(shù)識別準(zhǔn)確等優(yōu)點(diǎn)。
【專利說明】
一種強(qiáng)背景噪聲環(huán)境下的模態(tài)參數(shù)識別方法
技術(shù)領(lǐng)域
[0001] 本發(fā)明涉及模態(tài)測試領(lǐng)域，特別涉及一種強(qiáng)背景噪聲環(huán)境下的模態(tài)參數(shù)識別方 法。
【背景技術(shù)】
[0002] 模態(tài)參數(shù)識別技術(shù)在故障診斷、動力學(xué)響應(yīng)分析以及模型修正等領(lǐng)域發(fā)揮著重要 的作用，識別參數(shù)包括固有頻率，阻尼比和振型。常用的模態(tài)參數(shù)識別方法分為頻域和時(shí)域 兩大類。頻域和時(shí)域識別方法都以測試得到的脈沖響應(yīng)信號作為基礎(chǔ)數(shù)據(jù)進(jìn)行模態(tài)參數(shù)識 另IJ，因此獲得準(zhǔn)確可靠的脈沖響應(yīng)信號至關(guān)重要。由于信號在測量、傳輸過程中不可避免地 受到不同程度的噪聲污染，特別是在大型結(jié)構(gòu)的模態(tài)測試中，周圍電機(jī)、栗、空調(diào)等部件產(chǎn) 生的強(qiáng)噪聲對脈沖響應(yīng)信號影響較大，模態(tài)參數(shù)識別的準(zhǔn)確性受到極大影響。
[0003] 針對上述問題，目前主要采用統(tǒng)計(jì)平均、小波降噪和奇異值分解等方式來消除噪 聲對脈沖響應(yīng)信號的干擾，提高模態(tài)參數(shù)識別準(zhǔn)確性。其中，統(tǒng)計(jì)平均必須保證在測量過程 中參數(shù)不發(fā)生變化，否則采用統(tǒng)計(jì)平均沒有意義。而小波降噪和奇異值分解技術(shù)都需要根 據(jù)噪聲強(qiáng)弱設(shè)置不同閥值參數(shù)，否測降噪效果受影響。與此同時(shí)，目前沒有研究結(jié)果顯示， 當(dāng)強(qiáng)噪聲中的周期成分和固有頻率接近或重疊時(shí)，現(xiàn)有的方法能夠準(zhǔn)確識別模態(tài)參數(shù)。因 此，目前的模態(tài)參數(shù)識別方法存在局限性，還有改進(jìn)空間。

【發(fā)明內(nèi)容】

[0004] 為了克服上述現(xiàn)有技術(shù)的缺點(diǎn)，本發(fā)明的目的在于提供一種強(qiáng)背景噪聲環(huán)境下的 模態(tài)參數(shù)識別方法，提高模態(tài)參數(shù)識別精度，具有自適應(yīng)強(qiáng)，計(jì)算速度快，強(qiáng)噪聲環(huán)境下固 有頻率識別準(zhǔn)確等優(yōu)點(diǎn)。
[0005] 為了達(dá)到上述目的，本發(fā)明采取的技術(shù)方案為：
[0006] -種強(qiáng)背景噪聲環(huán)境下的模態(tài)參數(shù)識別方法，包括以下步驟：
[0007] 步驟一，通過力錘敲擊實(shí)驗(yàn)，測得一級脈沖響應(yīng)信號y(t)，測得的一級脈沖響應(yīng)信 號y(t)包含理想的脈沖響應(yīng)信號ip)和背景噪聲信號n(t)，t代表時(shí)間；
[0008] 步驟二，采用最優(yōu)平滑算法、最小統(tǒng)計(jì)的噪聲估計(jì)算法，或者最小受控遞歸平均法 中任意一種獲取一級脈沖響應(yīng)信號y(t)的背景噪聲估計(jì)值
[0009] 步驟三，對測得的一級脈沖響應(yīng)信號y(t)和一級脈沖響應(yīng)信號y(t)的背景噪聲估 計(jì)仉進(jìn)行分幀，加漢寧窗和FFT處理，得到y(tǒng)(t)的幅度譜Y(p，k)和相位譜Φ(ρ，1〇,以及 ??的幅度譜々&,句，1(?， 1〇表示7(仂的第?幀信號的第1^根譜分量幅值，<1)(?， 1〇表示7 (t)的第ρ幀信號的第k根譜分量相位表示丨的第ρ幀信號的第k根譜分量幅值， 分幀時(shí)，每幀長度為0.02倍的采樣頻率，前后幀之間采用50%的重疊度；
[0010] 步驟四，通過譜減法對測得一級脈沖響應(yīng)信號y(t)的幅度譜Y(p，k)進(jìn)行降噪處理，得 到初步降噪后幅度譜A ο，計(jì)算公式為
其中，譜減系數(shù)α = 9，β = 0.05;
[0011] 步驟五，將步驟三得到的相位譜Φ(ρ，1〇和步驟四中得到的初步降噪后幅度譜 $ 結(jié)合，通過傅立葉逆變換和重疊相加，得到二級脈沖響應(yīng)信號yi(t);
[0012] 步驟六，采用最優(yōu)平滑算法、最小統(tǒng)計(jì)的噪聲估計(jì)算法，或者最小受控遞歸平均法 中任意一種獲取二級脈沖響應(yīng)信號yi(t)的背景噪聲估計(jì)值
[0013] 步驟七，對二級脈沖響應(yīng)信號71(〇和二級脈沖響應(yīng)信號yi(t)的背景噪聲估計(jì)值 進(jìn)行分幀，加漢寧窗和FFT處理，得到二級脈沖響應(yīng)信號 yi(t)的幅度譜心化少)和相位 譜巾1(?，1〇，以及爲(wèi)(《)的幅度譜次（/^)，￥1(?，1〇表示 71(〇的第?幀信號的第奸艮譜分量幅 值，Φ 1 (p，k)表不y 1⑴的第p幀信號的第k根譜分量相位，來:(表不?1⑴的第p幀信號的 第k根譜分量幅值，分幀時(shí)，每幀長度為0.2倍的采樣頻率，前后幀之間采用50%的重疊度；
[0014] 步驟八，計(jì)算二級脈沖響應(yīng)信號yi(t)的每一幀后驗(yàn)信噪比估計(jì)值，計(jì)算公式為：
,其中，代表汶(於太)的方差，來如，^)與Yi(P，k)已在步 驟五中計(jì)算得到；
[0015] 步驟九，通過最小均方差短時(shí)譜估計(jì)方法對二級脈沖響應(yīng)信號yi(t)的幅度譜￥: (P，k)進(jìn)行二次降噪處理，得到二次降噪后幅度譜該計(jì)算過程采用循環(huán)方式進(jìn) 行，循環(huán)從第一幀開始到最后一幀結(jié)束，依次進(jìn)行下面①~⑤步計(jì)算，步驟如下：
[0016] ①使用經(jīng)典DD算法，計(jì)算二級脈沖響應(yīng)信號yi(t)的第p幀先驗(yàn)信噪比估計(jì)值，計(jì) 算公式為：
[0017]
，其中，衰減系數(shù)y = 〇.98;4(p,0已在步驟八中計(jì)算得到；的初始值設(shè)置為響量，以后的每一幀 使用第⑤步計(jì)算結(jié)果；（/; + U)代表第p+1幀信號的第k根譜分量的后驗(yàn)信噪比估計(jì) 值，該值已在步驟八中計(jì)算得到；
[0018] ②計(jì)算二級脈沖響應(yīng)信號yi(t)的第p幀譜分量的噪聲抑制因子，計(jì)算公式為： ，其中使用第①步計(jì)算結(jié)果；
[0019] ③使用先驗(yàn)信噪比估算方法，計(jì)算二級脈沖響應(yīng)信號yi(t)的第p幀先驗(yàn)信噪比估 計(jì)值，計(jì)算公式為：
[0020]
其中，衰 減系數(shù)f = ο. 98;Υι(p，k)已在步驟七中計(jì)算得到；j;;! (/·?.々）已在步驟八中計(jì)算得到;Gdd(p， k)使用第②步計(jì)算結(jié)果;,+ U ；)已在步驟八中計(jì)算得到；
[0021]④計(jì)算二級脈沖響應(yīng)信號yi(t)的第p幀譜分量的噪聲抑制因子，計(jì)算公式為：
其中使用第③步計(jì)算結(jié)果；
[0022]⑤通過最小均方差短時(shí)譜估計(jì)方法對心化氺)的第p幀譜分量進(jìn)行二次降噪處理， 得到二次降噪后幅度譜i(/U)，計(jì)算公式為：義(#) = (;(/〇)<(#)，其中6(?，10使用 第④步計(jì)算結(jié)果;Yi(p，k)已在步驟七中計(jì)算得到；
[0023] 步驟十，將步驟七中得到的相位譜ΦΚρ,?Ο和采用最小均方差短時(shí)譜估計(jì)方法后 得到的二次降噪后幅度譜結(jié)合，通過傅立葉逆變換和重疊相加，得到理想的脈沖響 應(yīng)信號:
[0024] 步驟十一，通過模態(tài)參數(shù)識別方法對理想的脈沖響應(yīng)信號進(jìn)行模態(tài)參數(shù)識 別。
[0025] 本發(fā)明的有益效果:研究發(fā)現(xiàn)，該方法可有效消除測試脈沖響應(yīng)信號中的平穩(wěn)信 號部分和非平穩(wěn)信號中的連續(xù)部分，保留非平穩(wěn)信號中的瞬態(tài)部分，即理想的脈沖響應(yīng)信 號，從而提高模態(tài)參數(shù)識別精度。
[0026] 相比于現(xiàn)有方法，該方法具有自適應(yīng)強(qiáng)，計(jì)算速度快，強(qiáng)噪聲環(huán)境下固有頻率識別 準(zhǔn)確等優(yōu)點(diǎn)。該方法即使在測試噪聲強(qiáng)烈，測試噪聲中周期性成分和結(jié)構(gòu)固有頻率極為接 近情況下，也可以準(zhǔn)確識別結(jié)構(gòu)固有頻率。
【附圖說明】
[0027] 圖1為用于模態(tài)參數(shù)識別的汽輪發(fā)電機(jī)殼體示意圖。
[0028] 圖2為振動測點(diǎn)的一級脈沖響應(yīng)信號y(t)時(shí)域圖。
[0029]圖3為振動測點(diǎn)的一級脈沖響應(yīng)信號y(t)頻域圖。
[0030] 圖4為振動測點(diǎn)的一級脈沖響應(yīng)信號y(t)經(jīng)本發(fā)明處理后的時(shí)域圖。
[0031] 圖5為振動測點(diǎn)的一級脈沖響應(yīng)信號y(t)經(jīng)本發(fā)明處理后的頻域圖。
【具體實(shí)施方式】
[0032]下面結(jié)合附圖和實(shí)施案例對本發(fā)明進(jìn)行詳細(xì)描述。
[0033] 某核電廠的發(fā)電機(jī)組在運(yùn)行過程中存在發(fā)電機(jī)殼體結(jié)構(gòu)共振問題，為了解決共振 問題，需求準(zhǔn)確識別發(fā)電機(jī)殼體的結(jié)構(gòu)固有頻率，由于現(xiàn)場環(huán)境噪聲強(qiáng)烈，并且測試噪聲中 周期性成分和發(fā)電機(jī)殼體固有頻率極為接近，現(xiàn)有的方法不能準(zhǔn)確識別發(fā)電機(jī)殼體固有頻 率。下面通過本發(fā)明解決該問題。
[0034] -種強(qiáng)背景噪聲環(huán)境下的模態(tài)參數(shù)識別方法，包括以下步驟：
[0035] 步驟一，通過力錘敲擊試驗(yàn)，測得發(fā)電機(jī)殼體上振動測點(diǎn)的一級脈沖響應(yīng)信號y (t)，采樣頻率為500Hz，參照圖1，圖1顯示的是汽輪發(fā)電機(jī)殼體示意圖，發(fā)電機(jī)3的機(jī)殼僅依 靠重力放置在安裝基礎(chǔ)4上，發(fā)電機(jī)3的轉(zhuǎn)子5通過滑動軸承連接在發(fā)電機(jī)端蓋2上，發(fā)電機(jī)3 的機(jī)殼上部連接有2個(gè)發(fā)電機(jī)氫冷器1，振動測點(diǎn)位置設(shè)置在1個(gè)發(fā)電機(jī)氫冷器1上；圖2和圖 3分別是通過該步驟后獲得的發(fā)電機(jī)殼體上振動測點(diǎn)的一級脈沖響應(yīng)信號時(shí)域波形和頻域 波形圖，由于汽輪發(fā)電機(jī)周圍的環(huán)境噪聲較大，圖2中顯示的一級脈沖響應(yīng)信號已經(jīng)被噪聲 信號污染，圖3中顯示的一個(gè)突出頻率(49.82Hz)為周期性噪聲；
[0036] 步驟二，采用最優(yōu)平滑算法獲取一級脈沖響應(yīng)信號y(t)的背景噪聲估計(jì)值
[0037] 步驟三，對測得的一級脈沖響應(yīng)信號y(t)和一級脈沖響應(yīng)信號y(t)的背景噪聲估 計(jì)值Μ"進(jìn)行分幀，加漢寧窗和FFT處理，得到y(tǒng)(t)的幅度譜Y(p，k)和相位譜Φ(ρ，1〇,以及 /ψ)的幅度譜ihM，分幀時(shí)，每幀長度為〇.02倍的采樣頻率，前后幀之間采用50%的重 疊度；
[0038]步驟四，通過譜減法對一級脈沖響應(yīng)信號y(t)的幅度譜Y(p，k)進(jìn)行降噪處理，得 到初步降噪后幅度譜$〇?，幻，計(jì)算公式為：
[0040] 其中，譜減系數(shù)α = 9，β = 〇·〇5;
[0041] 步驟五，將步驟三得到的相位譜Φ(ρ，1〇和步驟四中得到的初步降噪后幅度譜 f if"結(jié)合，通過傅立葉逆變換和重疊相加，得到二級脈沖響應(yīng)信號yi(t);
[0042] 步驟六，采用最優(yōu)平滑算法獲取二級脈沖響應(yīng)信號71(〇的背景噪聲估計(jì)值^
[0043] 步驟七，對二級脈沖響應(yīng)信號71(〇和二級脈沖響應(yīng)信號yi(t)的背景噪聲估計(jì)值 %(0進(jìn)行分幀，加漢寧窗和1^處理，得到二級脈沖響應(yīng)信號 71(〇的幅度譜￥1(?，1〇和相位 譜Φ i(P，k)，以及木（?)的幅度譜汶（/^)，分幀時(shí)，每幀長度為〇. 2倍的采樣頻率，前后幀之 間采用50%的重疊度；
[0044] 步驟八，計(jì)算二級脈沖響應(yīng)信號yi(t)的每一幀后驗(yàn)信噪比估計(jì)值，計(jì)算公式為：
其中，代表的方差與Υι(Ρ，1〇已在步驟 五中計(jì)算得到；
[0045] 步驟九，通過最小均方差短時(shí)譜估計(jì)方法對二級脈沖響應(yīng)信號yi(t)的幅度譜￥: (P，k)進(jìn)行二次降噪處理，得到二次降噪后的幅度譜該計(jì)算過程采用循環(huán)方式進(jìn) 行，循環(huán)從第一幀開始到最后一幀結(jié)束，依次進(jìn)行下面①~⑤步計(jì)算，步驟如下：
[0046] ①使用經(jīng)典DD算法，計(jì)算二級脈沖響應(yīng)信號yi(t)的第ρ幀先驗(yàn)信噪比估計(jì)值，計(jì) 算公式為：
[0047]
其中，衰減系數(shù)f = 0.98;^〇M〇已在步驟八中計(jì)算得到；的初始值設(shè)置為響量，以后的每一幀 使用第⑤步計(jì)算結(jié)果；(/; + 1,/Γ)代表第ρ+l幀信號的第k根譜分量的后驗(yàn)信噪比估計(jì) 值，該值已在步驟八中計(jì)算得到；
[0048] ②計(jì)算二級脈沖響應(yīng)信號yi(t)的第p幀譜分量的噪聲抑制因子，計(jì)算公式為：
，其中使用第①步計(jì)算結(jié)果；
[0049] ③使用改進(jìn)的先驗(yàn)信噪比估算方法，計(jì)算二級脈沖響應(yīng)信號yi(t)的第ρ幀先驗(yàn)信 噪比估計(jì)值，計(jì)算公式為：
[0050]
其中， 衰減系數(shù)f =〇.98;￥1(口，1〇已在步驟七中計(jì)算得到；(^1(/;>,0已在步驟八中計(jì)算得到 ;6〇〇 (P，k)使用第②步計(jì)算結(jié)果; STkp, (/M- U)已在步驟八中計(jì)算得到；
[0051] ④計(jì)算二級脈沖響應(yīng)信號yi(t)的第p幀譜分量的噪聲抑制因子，計(jì)算公式為：
，其中(菸芍使用第③步計(jì)算結(jié)果；
[0052]⑤通過最小均方差短時(shí)譜估計(jì)方法對心化氺)的第p幀譜分量進(jìn)行二次降噪處理， 得到二次降噪后幅度譜，計(jì)算公式為：.(_(/以）=6'(/，人丨4(/以），其中6化，10使用 第④步計(jì)算結(jié)果;Yi(p，k)已在步驟七中計(jì)算得到；
[0053] 步驟十，將步驟七中得到的相位譜ΦΚρ,ΙΟ和采用最小均方差短時(shí)譜估計(jì)方法后 得到的二次降噪后幅度譜1(八/〇結(jié)合，通過傅立葉逆變換和重疊相加，得到理想的脈沖響 應(yīng)信號，參照圖4,圖4是理想的脈沖響應(yīng)信號交的的時(shí)域波形圖；
[0054] 步驟十一，對理想的脈沖響應(yīng)信號進(jìn)行傅立葉變換，從對應(yīng)的頻域圖中識別 發(fā)電機(jī)殼體固有頻率，圖4顯示一級脈沖響應(yīng)信號y(t)中的環(huán)境噪聲已被去除；圖5是理想 的脈沖響應(yīng)信號的頻域圖，從圖5中可以清晰的觀察到發(fā)電機(jī)殼體存在48.99Hz的固有 頻率，進(jìn)而可判斷發(fā)電機(jī)殼體結(jié)構(gòu)共振問題是由于發(fā)電機(jī)殼體的固有頻率(48.99Hz)和工 作轉(zhuǎn)頻(50Hz)過于接近所致。
[0055] 以上內(nèi)容是結(jié)合具體的優(yōu)選實(shí)施方式對本發(fā)明所作的進(jìn)一步詳細(xì)說明，不能認(rèn)定 本發(fā)明的【具體實(shí)施方式】僅限于此，對于本發(fā)明所屬技術(shù)領(lǐng)域的普通技術(shù)人員來說，在不脫 離本發(fā)明構(gòu)思的前提下，還可以做出若干簡單的推演或替換，都應(yīng)當(dāng)視為屬于本發(fā)明由所 提交的權(quán)利要求書確定專利保護(hù)范圍。
【主權(quán)項(xiàng)】
1. 一種強(qiáng)背景噪聲環(huán)境下的模態(tài)參數(shù)識別方法，其特征在于，包括W下步驟： 步驟一，通過力鍵敲擊實(shí)驗(yàn)，測得一級脈沖響應(yīng)信號y(t)，測得的一級脈沖響應(yīng)信號y (t)包含理想的脈沖響應(yīng)信號和背景噪聲信號n(t)，t代表時(shí)間； 步驟二，采用最優(yōu)平滑算法、最小統(tǒng)計(jì)的噪聲估計(jì)算法，或者最小受控遞歸平均法中任 意一種獲取一級脈沖響應(yīng)信號y(t)的背景噪聲估計(jì)值Α??); 步驟Ξ，對測得的一級脈沖響應(yīng)信號y(t)和一級脈沖響應(yīng)信號y(t)的背景噪聲估計(jì)值 進(jìn)行分帖，加漢寧窗和FFT處理，得到y(tǒng)(t)的幅度譜Y(p，k)和相位譜(1)(9，1〇，^及&(。 的幅度譜·?，Y(P，k)表示y(t)的第P帖信號的第k根譜分量幅值，Φ (p，k)表示y(t)的 第P帖信號的第k根譜分量相位，表示的的第P帖信號的第k根譜分量幅值，分帖 時(shí)，每帖長度為0.02倍的采樣頻率，前后帖之間采用50 %的重疊度； 步驟四，通過譜減法對測得一級脈沖響應(yīng)信號y(t)的幅度譜Y(p，k)進(jìn)行降噪處理，得 到初步降噪后幅度譜，計(jì)算公式為：其中，譜減系數(shù)α = 9，β = 0.05; 步驟五，將步驟Ξ得到的相位譜Φ (p，k)和步驟四中得到的初步降噪后幅度譜《（戶，叫 結(jié)合，通過傅立葉逆變換和重疊相加，得到二級脈沖響應(yīng)信號yi(t); 步驟六，采用最優(yōu)平滑算法、最小統(tǒng)計(jì)的噪聲估計(jì)算法，或者最小受控遞歸平均法中任 意一種獲取二級脈沖響應(yīng)信號yi(t)的背景噪聲估計(jì)值為(句; 步驟屯，對二級脈沖響應(yīng)信號yi(t)和二級脈沖響應(yīng)信號yi(t)的背景噪聲估計(jì)值卻（?) 進(jìn)行分帖，加漢寧窗和FFT處理，得到二級脈沖響應(yīng)信號yi(t)的幅度譜Yi(p，k)和相位譜Φι (口，1〇，^及/1(〇的幅度譜々1(護(hù)*)訊(口，1〇表示71(〇的第口帖信號的第村良譜分量幅值， Φι(Ρ，1〇表示yi(t)的第Ρ帖信號的第k根譜分量相位，瑪.(扔fc):表示的第Ρ帖信號的第 k根譜分量幅值，分帖時(shí)，每帖長度為0.2倍的采樣頻率，前后帖之間采用50%的重疊度； 步驟八，計(jì)算二級脈沖響應(yīng)信號yi(t)的每一帖后驗(yàn)信噪比估計(jì)值，計(jì)算公式為：其中，妨(護(hù)A)代表成(扔巧的方差，瑪(巧&)與Yi(p，k)已在步驟 五中計(jì)算得到； 步驟九，通過最小均方差短時(shí)譜估計(jì)方法對二級脈沖響應(yīng)信號yi(t)的幅度譜Yi(p，k) 進(jìn)行二次降噪處理，得到二次降噪后幅度譜乂(、,"，叫，該計(jì)算過程采用循環(huán)方式進(jìn)行，循 環(huán)從第一帖開始到最后一帖結(jié)束，依次進(jìn)行下面①~⑤步計(jì)算，步驟如下： ①使用經(jīng)典DD算法，計(jì)算二級脈沖響應(yīng)信號yi(t)的第P帖先驗(yàn)信噪比估計(jì)值，計(jì)算公式 為：其中，衰減系數(shù)β/ = 0.98;巧(礦A)已在步驟八中計(jì)算得到；l(f-：U)的初始值設(shè)置為響量，W后的每一帖使 用第⑤步計(jì)算結(jié)果;+ 代表第P+1帖信號的第k根譜分量的后驗(yàn)信噪比估計(jì)值， 該值已在步驟八中計(jì)算得到； ② 計(jì)算二級脈沖響應(yīng)信號yi(t)的第P帖譜分量的噪聲抑制因子，計(jì)算公式為：其中沒?Κ篇〇α)使用第①步計(jì)算結(jié)果； ③ 使用先驗(yàn)信噪比估算方法，計(jì)算二級脈沖響應(yīng)信號yi(t)的第Ρ帖先驗(yàn)信噪比估計(jì)值， 計(jì)算公式為：其中，畏減系 數(shù)e'=〇.98;Yi(p，k)已在步驟屯中計(jì)算得到;巧批句已在步驟八中計(jì)算得到;GDD(p，k)使 用第②步計(jì)算結(jié)果;S,i/j十U ：)已在步驟八中計(jì)算得到； ④ 計(jì)算二級脈沖響應(yīng)信號yi(t)的第P帖譜分量的噪聲抑制因子，計(jì)算公式為：其中巧咸黑A yt)使用第③步計(jì)算結(jié)果. ⑤ 通過最小均方差短時(shí)譜估計(jì)方法對Yi(P，k)的第P帖譜分量進(jìn)行二次降噪處理，得到 二次降噪后幅度譜樂?/α)，計(jì)算公式為：去(、/^)二6'(、.化/:).}^/^)，其中6(9，1〇使用第 ④步計(jì)算結(jié)果;Yi(P，k)已在步驟屯中計(jì)算得到； 步驟十，將步驟屯中得到的相位譜Φι(Ρ，1〇和采用最小均方差短時(shí)譜估計(jì)方法后得到 的二次降噪后幅度譜^(/，，/〇結(jié)合，通過傅立葉逆變換和重疊相加，得到理想的脈沖響應(yīng)信 號《(〇; 步驟十一，通過模態(tài)參數(shù)識別方法對理想的脈沖響應(yīng)信號.^0進(jìn)行模態(tài)參數(shù)識別。
【文檔編號】G06K9/00GK106096530SQ201610397420
【公開日】2016年11月9日
【申請日】2016年6月7日 公開號201610397420.7, CN 106096530 A, CN 106096530A, CN 201610397420, CN-A-106096530, CN106096530 A, CN106096530A, CN201610397420, CN201610397420.7
【發(fā)明人】王琇峰, 和丹, 曾志豪, 林京, 雷亞國
【申請人】西安交通大學(xué)