本發(fā)明涉及，特別涉及一種脈沖響應(yīng)信號(hào)聲強(qiáng)測(cè)量的方法。
背景技術(shù)：
：聲強(qiáng)表征聲波能量的強(qiáng)度。聲強(qiáng)測(cè)量方法廣泛應(yīng)用于聲功率級(jí)計(jì)算、振動(dòng)源識(shí)別等領(lǐng)域。工程師通過(guò)觀(guān)察聲強(qiáng)頻譜上各頻率點(diǎn)對(duì)應(yīng)的聲強(qiáng)值，判斷各振動(dòng)部件實(shí)際貢獻(xiàn)的聲能量，進(jìn)而采取針對(duì)性的降噪措施。聲強(qiáng)測(cè)量技術(shù)一直是人們關(guān)注的問(wèn)題。目前，最常用的聲強(qiáng)測(cè)量方法是雙傳聲器互譜法，傅里葉變換是互譜法測(cè)聲強(qiáng)的核心步驟，傅里葉變換將雙傳聲器a、b通道測(cè)得的時(shí)間域聲壓信號(hào)pa(t)、pb(t)變換成頻域聲壓信號(hào)pa(ω)、pb(ω)，計(jì)算其互譜后代入頻域聲強(qiáng)計(jì)算公式求出聲強(qiáng)。對(duì)于穩(wěn)態(tài)正、余弦噪聲波，采用互譜法測(cè)量計(jì)算其聲強(qiáng)，獲取的頻譜能夠有效并合理地反映出聲波的特征。但是在實(shí)際工程中，聲波存在形式眾多，脈沖沖擊是一種常見(jiàn)的激勵(lì)形式，沖擊產(chǎn)生脈沖響應(yīng)。采用互譜法計(jì)算脈沖響應(yīng)波的聲強(qiáng)時(shí)，傅里葉變換將沖擊響應(yīng)信號(hào)分解為多個(gè)不同頻率的正、余弦信號(hào)，導(dǎo)致頻譜上的聲強(qiáng)值分散在多個(gè)不同的頻率點(diǎn)上，甚至可能在部分頻率點(diǎn)上出現(xiàn)負(fù)聲強(qiáng)，測(cè)量計(jì)算結(jié)果反映脈沖響應(yīng)噪聲聲強(qiáng)不夠簡(jiǎn)便，物理特征不明顯。本發(fā)明提出一種新的聲強(qiáng)測(cè)量計(jì)算方法。該方法不以傅里葉變換為核心，而是按照脈沖響應(yīng)信號(hào)其本身固有特征來(lái)構(gòu)造函數(shù)，采用內(nèi)積方法計(jì)算出脈沖響應(yīng)信號(hào)的各項(xiàng)特征參數(shù)，根據(jù)獲取的特征參數(shù)重構(gòu)脈沖響應(yīng)信號(hào)，最后將重構(gòu)信號(hào)直接代入時(shí)域聲強(qiáng)積分公式計(jì)算，求解出脈沖響應(yīng)聲波的聲強(qiáng)值。本發(fā)明方法可以有效地測(cè)量識(shí)別出脈沖響應(yīng)聲波的特征，物理意義明確。技術(shù)實(shí)現(xiàn)要素：針對(duì)現(xiàn)有技術(shù)中存在不足，本發(fā)明提供了一種脈沖響應(yīng)信號(hào)聲強(qiáng)測(cè)量的方法，可以將每個(gè)脈沖響應(yīng)信號(hào)從原始信號(hào)中提取出來(lái)，直接計(jì)算其聲強(qiáng)，這樣清晰地測(cè)量識(shí)別出脈沖響應(yīng)信號(hào)的聲強(qiáng)值，有助于準(zhǔn)確地判斷脈沖響應(yīng)聲波的特性。本發(fā)明是通過(guò)以下技術(shù)手段實(shí)現(xiàn)上述技術(shù)目的的。一種脈沖響應(yīng)信號(hào)聲強(qiáng)測(cè)量的方法，包括如下步驟：s01：用雙傳聲器探頭a、b采集脈沖響應(yīng)信號(hào)聲壓數(shù)據(jù)，將測(cè)得的聲壓數(shù)據(jù)導(dǎo)入計(jì)算機(jī)，分別記為p'a(t)、pb'(t)；s02：記錄信號(hào)p'a(t)和p'b(t)的時(shí)間長(zhǎng)度為t，對(duì)信號(hào)p'a(t)和p'b(t)做離散傅里葉變換，通過(guò)頻譜的峰值個(gè)數(shù)判定通道中包含的脈沖響應(yīng)信號(hào)數(shù)量為m，其中探頭a第i個(gè)脈沖響應(yīng)信號(hào)設(shè)為pai(t)，探頭b第i個(gè)脈沖響應(yīng)信號(hào)設(shè)為pbi(t)，pa(t)和pb(t)分別記為：其中0<t<t；s03：確定脈沖響應(yīng)信號(hào)pai(t)的固有頻率值和脈沖響應(yīng)信號(hào)pbi(t)的固有頻率值構(gòu)建信號(hào)：當(dāng)初始值i＝1時(shí)，la0(t)＝pa(t)、lb0(t)＝pb(t)；構(gòu)建信號(hào)lai-1(t)做傅里葉變換，確定最大幅值及對(duì)應(yīng)頻率點(diǎn)，確定構(gòu)建信號(hào)lai-1(t)的頻率估計(jì)值確定頻率搜索區(qū)間在頻率搜索區(qū)間內(nèi)獲取h個(gè)均分的細(xì)分頻率點(diǎn)其中，w＝0,1,…,h；構(gòu)造函數(shù)將構(gòu)造函數(shù)與構(gòu)建信號(hào)lai-1(t)做內(nèi)積運(yùn)算，尋找內(nèi)積模極值點(diǎn)，根據(jù)極值點(diǎn)處的頻率，獲取響應(yīng)信號(hào)pai(t)的固有頻率值構(gòu)建信號(hào)lbi-1(t)做傅里葉變換，確定最大幅值及對(duì)應(yīng)頻率點(diǎn)，確定構(gòu)建函數(shù)lbi-1(t)的頻率估計(jì)值確定頻率搜索區(qū)間在頻率搜索區(qū)間內(nèi)獲取h個(gè)均分的細(xì)分頻率點(diǎn)其中，w＝0,1,…,h；構(gòu)造函數(shù)將構(gòu)造函數(shù)與構(gòu)建信號(hào)lbi-1(t)做內(nèi)積運(yùn)算，尋找內(nèi)積模極值點(diǎn)，根據(jù)極值點(diǎn)處的頻率，獲取響應(yīng)信號(hào)pbi(t)的固有頻率值s04：根據(jù)脈沖響應(yīng)信號(hào)pai(t)的固有頻率值構(gòu)造在時(shí)間軸上移動(dòng)的函數(shù)：根據(jù)脈沖響應(yīng)信號(hào)pbi(t)的固有頻率值構(gòu)造在時(shí)間軸上移動(dòng)的函數(shù)：其中：δt為函數(shù)的時(shí)間長(zhǎng)度，t為變量且t∈(0,t-δt)；構(gòu)建信號(hào)lai-1(t)與構(gòu)造函數(shù)cai(t)做內(nèi)積運(yùn)算，構(gòu)建函數(shù)lbi-1(t)與構(gòu)造函數(shù)cbi(t)做內(nèi)積運(yùn)算，確定內(nèi)積模極大值，根據(jù)確定的內(nèi)積模極大值，從而確定脈沖響應(yīng)信號(hào)pai(t)的起始時(shí)間信號(hào)pbi(t)的起始時(shí)間s05：截取構(gòu)建信號(hào)lai-1(t)在時(shí)間段上的聲壓信號(hào)，記為截取構(gòu)建信號(hào)lbi-1(t)在時(shí)間段上的聲壓信號(hào)，記為構(gòu)造函數(shù)：其中，nai和nbi為變量；對(duì)構(gòu)造函數(shù)歸一化處理得到：yai(t)＝uai(t)/||uai(t)||ybi(t)＝ubi(t)/||ubi(t)將函數(shù)yai(t)與截取信號(hào)進(jìn)行內(nèi)積運(yùn)算，將函數(shù)ybi(t)與截取信號(hào)進(jìn)行內(nèi)積運(yùn)算，得到xai、xbi的表達(dá)式：采用逐級(jí)細(xì)化方法尋找內(nèi)積模||xai||和||xbi||的最大值，根據(jù)xai、xbi的表達(dá)式計(jì)算出衰減系數(shù)和s06：根據(jù)內(nèi)積值xmaxai、xmaxbi和構(gòu)造函數(shù)uai(t)、ubi(t)計(jì)算出數(shù)值：其中：xmaxai、xmaxbi分別為||xai||、||xbi||取得最大值時(shí)對(duì)應(yīng)的數(shù)值；根據(jù)計(jì)算得出cai和cbi的相位和幅值：其中：arg(cai)表示求cai的輻角，arg(cbi)表示求cbi輻角；|cai|表示求復(fù)數(shù)cai的絕對(duì)值，|cbi|表示求復(fù)數(shù)cbi的絕對(duì)值；s07：根據(jù)脈沖響應(yīng)信號(hào)pai(t)的起始時(shí)間固有頻率衰減系數(shù)幅值初始相位脈沖響應(yīng)信號(hào)pbi(t)的起始時(shí)間固有頻率衰減系數(shù)幅值初始相位參數(shù)重構(gòu)脈沖響應(yīng)信號(hào)，分別記為其中：s08：將重構(gòu)脈沖響應(yīng)信號(hào)代入瞬時(shí)聲強(qiáng)計(jì)算公式：其中，ρ：空氣密度；d：探頭a和b的間距；取信號(hào)瞬時(shí)聲強(qiáng)的時(shí)間平均值，得到此脈沖響應(yīng)信號(hào)的聲強(qiáng)：s09：若i<m時(shí)，則i＝i+1，對(duì)lai-1(t)和lbi-1(t)執(zhí)行第s03-s08，計(jì)算出剩下m-1個(gè)脈沖響應(yīng)信號(hào)的聲強(qiáng)；當(dāng)i＝m時(shí)，循環(huán)結(jié)束。進(jìn)一步，所述s03步驟中本發(fā)明的有益效果：1.本發(fā)明所述的脈沖響應(yīng)信號(hào)聲強(qiáng)測(cè)量的方法，可以不使用傅里葉變換分解脈沖響應(yīng)信號(hào)，而是將每個(gè)脈沖響應(yīng)信號(hào)從原始信號(hào)中提取出來(lái)，直接計(jì)算其聲強(qiáng)。2.本發(fā)明所述的脈沖響應(yīng)信號(hào)聲強(qiáng)測(cè)量的方法，可以清晰地測(cè)量識(shí)別出脈沖響應(yīng)信號(hào)的聲強(qiáng)值，有助于準(zhǔn)確地判斷脈沖響應(yīng)聲波的特性。3.本發(fā)明所述的脈沖響應(yīng)信號(hào)聲強(qiáng)測(cè)量的方法，可以避免出現(xiàn)負(fù)聲強(qiáng)。附圖說(shuō)明圖1為本發(fā)明所述的脈沖響應(yīng)信號(hào)聲強(qiáng)測(cè)量的方法的流程圖。圖2為本發(fā)明所述的實(shí)施例的脈沖響應(yīng)信號(hào)。圖3為本發(fā)明所述的使用互譜法測(cè)量實(shí)施例得到的聲強(qiáng)級(jí)頻譜圖。具體實(shí)施方式下面結(jié)合附圖以及具體實(shí)施例對(duì)本發(fā)明作進(jìn)一步的說(shuō)明，但本發(fā)明的保護(hù)范圍并不限于此。如圖1所示，一種脈沖響應(yīng)信號(hào)聲強(qiáng)測(cè)量的方法，包括如下步驟：s01：用雙傳聲器探頭a、b采集脈沖響應(yīng)信號(hào)聲壓數(shù)據(jù)，將測(cè)得的聲壓數(shù)據(jù)導(dǎo)入計(jì)算機(jī)，分別記為p'a(t)、pb'(t)；s02：記錄信號(hào)pa(t)和pb(t)的時(shí)間長(zhǎng)度為t，對(duì)信號(hào)pa(t)和pb(t)做離散傅里葉變換，通過(guò)頻譜的峰值個(gè)數(shù)判定通道中包含的脈沖響應(yīng)信號(hào)數(shù)量為m，其中探頭a第i個(gè)脈沖響應(yīng)信號(hào)設(shè)為pai(t)，探頭b第i個(gè)脈沖響應(yīng)信號(hào)設(shè)為pbi(t)，pa(t)和pb(t)分別記為：其中0<t<t；s03：確定脈沖響應(yīng)信號(hào)pai(t)的固有頻率值和脈沖響應(yīng)信號(hào)pbi(t)的固有頻率值構(gòu)建信號(hào)：當(dāng)初始值i＝1時(shí)，la0(t)＝pa(t)、lb0(t)＝pb(t)；構(gòu)建信號(hào)lai-1(t)做傅里葉變換，確定最大幅值及對(duì)應(yīng)頻率點(diǎn)，確定構(gòu)建信號(hào)lai-1(t)的頻率估計(jì)值確定頻率搜索區(qū)間在頻率搜索區(qū)間內(nèi)獲取h個(gè)均分的細(xì)分頻率點(diǎn)其中，w＝0,1,…,h；構(gòu)造函數(shù)將構(gòu)造函數(shù)與構(gòu)建信號(hào)lai-1(t)做內(nèi)積運(yùn)算，尋找內(nèi)積模極值點(diǎn)，根據(jù)極值點(diǎn)處的頻率，獲取響應(yīng)信號(hào)pai(t)的固有頻率值構(gòu)建信號(hào)lbi-1(t)做傅里葉變換，確定最大幅值及對(duì)應(yīng)頻率點(diǎn)，確定構(gòu)建函數(shù)lbi-1(t)的頻率估計(jì)值確定頻率搜索區(qū)間在頻率搜索區(qū)間內(nèi)獲取h個(gè)均分的細(xì)分頻率點(diǎn)其中，w＝0,1,…,h；構(gòu)造函數(shù)將構(gòu)造函數(shù)與構(gòu)建信號(hào)lbi-1(t)做內(nèi)積運(yùn)算，尋找內(nèi)積模極值點(diǎn)，根據(jù)極值點(diǎn)處的頻率，獲取響應(yīng)信號(hào)pbi(t)的固有頻率值s04：根據(jù)脈沖響應(yīng)信號(hào)pai(t)的固有頻率值構(gòu)造在時(shí)間軸上移動(dòng)的函數(shù)：根據(jù)脈沖響應(yīng)信號(hào)pbi(t)的固有頻率值構(gòu)造在時(shí)間軸上移動(dòng)的函數(shù)：其中：δt為函數(shù)的時(shí)間長(zhǎng)度，t為變量且t∈(0,t-δt)；構(gòu)建信號(hào)lai-1(t)與構(gòu)造函數(shù)cai(t)做內(nèi)積運(yùn)算，構(gòu)建函數(shù)lbi-1(t)與構(gòu)造函數(shù)cbi(t)做內(nèi)積運(yùn)算，確定內(nèi)積模極大值，根據(jù)確定的內(nèi)積模極大值，從而確定脈沖響應(yīng)信號(hào)pai(t)的起始時(shí)間信號(hào)pbi(t)的起始時(shí)間s05：截取構(gòu)建信號(hào)lai-1(t)在時(shí)間段上的聲壓信號(hào)，記為截取構(gòu)建信號(hào)lbi-1(t)在時(shí)間段上的聲壓信號(hào)，記為構(gòu)造函數(shù)：其中，nai和nbi為變量；對(duì)構(gòu)造函數(shù)歸一化處理得到：yai(t)＝uai(t)/||uai(t)||ybi(t)＝ubi(t)/||ubi(t)||將函數(shù)yai(t)與截取信號(hào)進(jìn)行內(nèi)積運(yùn)算，將函數(shù)ybi(t)與截取信號(hào)進(jìn)行內(nèi)積運(yùn)算，得到xai、xbi的表達(dá)式：采用逐級(jí)細(xì)化方法尋找內(nèi)積模||xai||和||xbi||的最大值，根據(jù)xai、xbi的表達(dá)式計(jì)算出衰減系數(shù)和s06：根據(jù)內(nèi)積值xmaxai、xmaxbi和構(gòu)造函數(shù)uai(t)、ubi(t)計(jì)算出數(shù)值：其中：xmaxai、xmaxbi分別為||xai||、||xbi||取得最大值時(shí)對(duì)應(yīng)的數(shù)值；根據(jù)計(jì)算得出cai和cbi的相位和幅值：其中：arg(cai)表示求cai的輻角，arg(cbi)表示求cbi輻角；|cai|表示求復(fù)數(shù)cai的絕對(duì)值，|cbi|表示求復(fù)數(shù)cbi的絕對(duì)值；s07：根據(jù)脈沖響應(yīng)信號(hào)pai(t)的起始時(shí)間固有頻率衰減系數(shù)幅值初始相位脈沖響應(yīng)信號(hào)pbi(t)的起始時(shí)間固有頻率衰減系數(shù)幅值初始相位參數(shù)重構(gòu)脈沖響應(yīng)信號(hào)，分別記為其中：s08：將重構(gòu)脈沖響應(yīng)信號(hào)代入瞬時(shí)聲強(qiáng)計(jì)算公式：其中，ρ：空氣密度；d：探頭a和b的間距；取信號(hào)瞬時(shí)聲強(qiáng)的時(shí)間平均值，得到此脈沖響應(yīng)信號(hào)的聲強(qiáng)：s09：若i<m時(shí)，則i＝i+1；對(duì)lai-1(t)和lbi-1(t)執(zhí)行第s03-s08，計(jì)算出剩下m-1個(gè)脈沖響應(yīng)信號(hào)的聲強(qiáng)；當(dāng)i＝m時(shí)，循環(huán)結(jié)束。所述s03步驟中時(shí)，計(jì)算出的固有頻率值較準(zhǔn)確。實(shí)施例：如圖2所示，假設(shè)存在一由3個(gè)脈沖響應(yīng)成分組成的平面聲波：其中，ai、ni、ti、fi、分別為第i(i＝1,2,3)個(gè)脈沖響應(yīng)信號(hào)的幅值、衰減系數(shù)、起始時(shí)間、固有頻率、初相位。具體參數(shù)值見(jiàn)表1。響應(yīng)信號(hào)2、3同時(shí)發(fā)生，響應(yīng)信號(hào)1在后面的時(shí)間上發(fā)生。響應(yīng)信號(hào)固有頻率(hz)起始時(shí)間(s)幅值(pa)衰減系數(shù)相位信號(hào)1450.211.3204.20.2信號(hào)2250.820.3204.60.2信號(hào)3100.150.3205.30.5表1：各脈沖響應(yīng)聲波信號(hào)的特征參數(shù)s01：用雙傳聲器探頭a、b采集脈沖響應(yīng)信號(hào)聲壓數(shù)據(jù)，將測(cè)得的聲壓數(shù)據(jù)導(dǎo)入計(jì)算機(jī)，探頭a在前端，探頭b在后端，分別記為p'a(t)、p'b(t)；信號(hào)pb(t)比信號(hào)pa(t)滯后一微小時(shí)間段其中d為a、b探頭的間距(單位：m)，c為空氣中聲音傳播速度(單位：m/s)。s02：記錄信號(hào)p'a(t)和p'b(t)的時(shí)間長(zhǎng)度為t為3s，對(duì)信號(hào)做離散傅里葉變換，通過(guò)頻譜的峰值個(gè)數(shù)判定通道中包含3個(gè)沖擊響應(yīng)信號(hào)，其中探頭a第i個(gè)脈沖響應(yīng)信號(hào)設(shè)為pai(t)，探頭b第i個(gè)脈沖響應(yīng)信號(hào)設(shè)為pbi(t)，pa(t)和pb(t)分別記為：其中0<t<3；s03：確定脈沖響應(yīng)信號(hào)pai(t)的固有頻率值和脈沖響應(yīng)信號(hào)pbi(t)的固有頻率值構(gòu)建信號(hào)：當(dāng)初始值i＝1時(shí)，la0(t)＝pa(t)、lb0(t)＝pb(t)；構(gòu)建信號(hào)lai-1(t)做傅里葉變換，確定最大幅值及對(duì)應(yīng)頻率點(diǎn)，確定構(gòu)建信號(hào)lai-1(t)的頻率估計(jì)值確定頻率搜索區(qū)間其中，在頻率搜索區(qū)間內(nèi)獲取h個(gè)均分的細(xì)分頻率點(diǎn)，當(dāng)h＝100時(shí)，則其中，w＝0,1,…,100；構(gòu)造函數(shù)將構(gòu)造函數(shù)與構(gòu)建信號(hào)lai-1(t)做內(nèi)積運(yùn)算，尋找內(nèi)積模極值點(diǎn)，根據(jù)極值點(diǎn)處的頻率，獲取響應(yīng)信號(hào)pai(t)的固有頻率值構(gòu)建信號(hào)lbi-1(t)做傅里葉變換，確定最大幅值及對(duì)應(yīng)頻率點(diǎn)，確定構(gòu)建函數(shù)lbi-1(t)的頻率估計(jì)值確定頻率搜索區(qū)間其中，在頻率搜索區(qū)間內(nèi)獲取100個(gè)均分的細(xì)分頻率點(diǎn)其中，w＝0,1,…,100；構(gòu)造函數(shù)將構(gòu)造函數(shù)與構(gòu)建信號(hào)lbi-1(t)做內(nèi)積運(yùn)算，尋找內(nèi)積模極值點(diǎn)，根據(jù)極值點(diǎn)處的頻率，獲取響應(yīng)信號(hào)pbi(t)的固有頻率值s04：根據(jù)脈沖響應(yīng)信號(hào)pai(t)的固有頻率值構(gòu)造在時(shí)間軸上移動(dòng)的函數(shù)：根據(jù)脈沖響應(yīng)信號(hào)pbi(t)的固有頻率值構(gòu)造在時(shí)間軸上移動(dòng)的函數(shù)：其中：δt為函數(shù)的時(shí)間長(zhǎng)度，取值0.02s。t為變量且t∈(0,t-δt)；構(gòu)建信號(hào)lai-1(t)與構(gòu)造函數(shù)cai(t)做內(nèi)積運(yùn)算，構(gòu)建函數(shù)lbi-1(t)與構(gòu)造函數(shù)cbi(t)做內(nèi)積運(yùn)算，確定內(nèi)積模極大值，根據(jù)確定的內(nèi)積模極大值，從而確定脈沖響應(yīng)信號(hào)pai(t)的起始時(shí)間信號(hào)pbi(t)的起始時(shí)間s05：截取構(gòu)建信號(hào)lai-1(t)在時(shí)間段上的聲壓信號(hào)，記為截取構(gòu)建信號(hào)lbi-1(t)在時(shí)間段上的聲壓信號(hào)，記為構(gòu)造函數(shù)：其中，nai和nbi為變量；對(duì)構(gòu)造函數(shù)歸一化處理得到：yai(t)＝uai(t)/||uai(t)||ybi(t)＝ubi(t)/||ubi(t)||將函數(shù)yai(t)與截取信號(hào)進(jìn)行內(nèi)積運(yùn)算，將函數(shù)ybi(t)與截取信號(hào)進(jìn)行內(nèi)積運(yùn)算，得到xai、xbi的表達(dá)式：采用逐級(jí)細(xì)化方法尋找內(nèi)積模||xai||和||xbi||的最大值，根據(jù)xai、xbi的表達(dá)式計(jì)算出衰減系數(shù)和s06：根據(jù)內(nèi)積值xmaxai、xmaxbi和構(gòu)造函數(shù)uai(t)、ubi(t)計(jì)算出數(shù)值：其中：xmaxai、xmaxbi分別為||xai||、||xbi||取得最大值時(shí)對(duì)應(yīng)的數(shù)值；根據(jù)計(jì)算得出cai和cbi的相位和幅值：其中：arg(cai)表示求cai的輻角，arg(cbi)表示求cbi輻角；|cai|表示求復(fù)數(shù)cai的絕對(duì)值，|cbi|表示求復(fù)數(shù)cbi的絕對(duì)值；s07：根據(jù)脈沖響應(yīng)信號(hào)pai(t)的起始時(shí)間固有頻率衰減系數(shù)幅值初始相位脈沖響應(yīng)信號(hào)pbi(t)的起始時(shí)間固有頻率衰減系數(shù)幅值初始相位參數(shù)重構(gòu)脈沖響應(yīng)信號(hào)，識(shí)別結(jié)果見(jiàn)表2，分別記為其中：重構(gòu)信號(hào)如圖2所示；表2：本發(fā)明方法提取的各脈沖響應(yīng)信號(hào)特征參數(shù)s08：將重構(gòu)脈沖響應(yīng)信號(hào)代入瞬時(shí)聲強(qiáng)計(jì)算公式：其中，ρ＝1.293：空氣密度(單位：kg/m3)；d＝0.01：探頭a和b的間距(單位：m)；取信號(hào)瞬時(shí)聲強(qiáng)的時(shí)間平均值，得到此脈沖響應(yīng)信號(hào)的聲強(qiáng)：s09：若i<m時(shí)，則i＝i+1，對(duì)lai-1(t)和lbi-1(t)執(zhí)行步驟3-步驟8，計(jì)算出剩下m-1個(gè)脈沖響應(yīng)信號(hào)的聲強(qiáng)，計(jì)算結(jié)果見(jiàn)表3；響應(yīng)信號(hào)理論聲強(qiáng)(w/m2)本發(fā)明計(jì)算結(jié)果(w/m2)相對(duì)誤差信號(hào)10.018131090.018067880.3486％信號(hào)20.016550580.016636220.5174％信號(hào)30.014275260.014339910.4528％表3：本發(fā)明方法聲強(qiáng)計(jì)算結(jié)果當(dāng)i＝m時(shí)，循環(huán)結(jié)束。對(duì)比采用雙傳聲器互譜法計(jì)算沖擊響應(yīng)信號(hào)的聲強(qiáng)，互譜法使用了很多項(xiàng)描述脈沖響應(yīng)信號(hào)的聲強(qiáng)，無(wú)法在表格中表示。采用1/3倍頻程方法表示。如圖3所示，可以看到聲強(qiáng)分散在不同的頻率點(diǎn)上，甚至出現(xiàn)了負(fù)聲強(qiáng)。而本方法計(jì)算出的聲強(qiáng)值結(jié)果準(zhǔn)確，與各個(gè)脈沖響應(yīng)信號(hào)相互對(duì)應(yīng)，物理意義明確。所述實(shí)施例為本發(fā)明的優(yōu)選的實(shí)施方式，但本發(fā)明并不限于上述實(shí)施方式，在不背離本發(fā)明的實(shí)質(zhì)內(nèi)容的情況下，本領(lǐng)域技術(shù)人員能夠做出的任何顯而易見(jiàn)的改進(jìn)、替換或變型均屬于本發(fā)明的保護(hù)范圍。當(dāng)前第1頁(yè)12