專利名稱：一種用聲學(xué)成像確定沖擊位置的方法
技術(shù)領(lǐng)域：
本發(fā)明涉及一種確定物體上沖擊位置的方法以及使用此方法的裝置。

背景技術(shù)：
一種現(xiàn)有的確定沖擊位置的方法是基于聲學(xué)成像。根據(jù)該方法，一個物體的表面被分成如同一個陣列的M個區(qū)域，另有N個聲學(xué)傳感器固定在這個表面上。每個聲學(xué)傳感器i用來接收由物體表面上的一個沖擊所生成的聲學(xué)信號，然后將一個感應(yīng)到的信號si(t)傳輸至處理器。
對于陣列中的每個區(qū)域k，該感應(yīng)信號si(t)被延遲一個時間延遲補(bǔ)償τik，然后被相加得到Pk。該時間延遲補(bǔ)償τik是根據(jù)區(qū)域k和傳感器i的各自位置決定的。如圖1a所示，如果該感應(yīng)信號si(t)被正確延遲，例如，如果沖擊發(fā)生于區(qū)域k中，延遲信號則為同向的，而Pk會得到一個脈沖波形。反之如圖1b所示，如果該感應(yīng)信號si(t)沒有被正確延遲，例如，如果沖擊未發(fā)生于區(qū)域k中，Pk的波形很寬且其最大值也低。這樣，對于陣列中的每個區(qū)域k，Pk的最大值都被計算，因此形成一個沖擊的聲波圖像。通過對Pk的比較就可以確定表面上沖擊的位置。
然而，在使用這種方法時，如果時間起始點并沒有被正確定義，就會導(dǎo)致了處理器的同步問題。繼而確定沖擊位置的產(chǎn)生就會的失敗。此外，這種方法的準(zhǔn)確性還依賴于沖擊信號的形狀，尤其是沖擊信號的持續(xù)時間。所以，沖擊信號持續(xù)時間越長，準(zhǔn)確性就越低。
本發(fā)明的目的是提供一種確定沖擊位置的方法，其準(zhǔn)確性不依賴于沖擊信號的波形。
本發(fā)明的另一個目的是提供一種確定沖擊位置的可靠的方法。


發(fā)明內(nèi)容
本發(fā)明提出一種確定表面上沖擊位置的方法，所述表面包括N個聲學(xué)傳感器和M個確定區(qū)域，N大于等于3，所述沖擊產(chǎn)生一個聲學(xué)信號，每個傳感器i接收所述聲學(xué)信號并將一個感應(yīng)信號si(t)傳送至處理器。所述方法包括以下步驟 (a)計算P個互相關(guān)乘積 其中， Si(ω)是傳感器i所感應(yīng)到的信號si(t)的傅立葉變換； Sj(ω)是傳感器j所感應(yīng)到的信號sj(t)的傅立葉變換； *是復(fù)共軛算子； (b)計算所述Pij(ω)的P個傅立葉反變換pij′(u)； (c)對于每個區(qū)域k，計算Pk(u)＝∑p′ij(u-τijk)以及Pk(u)的一個特征值，其中τijk是一個基于區(qū)域k及傳感器i和j各自位置的確定值。
(d)尋找k0，其中，當(dāng)k≠k0時，Pk0(u)的特征值大于相應(yīng)的Pk(u)的特征值。
在本發(fā)明的方法的不同實施例中，至少使用下列特征中的一個 -傅立葉反變化pij′(u)等于 其中Ω是一個波數(shù)系數(shù)；m(Ω)是對應(yīng)于波數(shù)系數(shù)Ω的頻率，取決于材料的色散關(guān)系； -u是距離，τijk是一個基于區(qū)域k和傳感器i之間距離，以及區(qū)域k和傳感器j之間距離的長度； -材料色散關(guān)系基本上等于其中α是取決于物體的一個系數(shù)； -材料色散關(guān)系基本上等于Ω＝αω，其中α是取決于物體的一個系數(shù)； -系數(shù)α取決于θik，其中θik是一個基于區(qū)域k和傳感器i之間位置的一個確定角度值； -系數(shù)α取決于物體的溫度α(T)，其中T是溫度值； -該方法還包括一種初始化模式，該初始化模式包括一個確定物體材料色散關(guān)系的步驟； -所述確定材料色散關(guān)系的步驟包括在所確定位置產(chǎn)生NSO沖擊的子步驟； -所述NSO沖擊產(chǎn)生于N個傳感器中的兩個傳感器之間的連線上的確定位置，并且均勻分布； -確定材料色散關(guān)系的步驟還包括以下子步驟 *根據(jù)α的測試值計算p′ij(u)； *對于每個α的測試值，對于每個沖擊m，計算 Pm(u)＝∑p′ij(u-τijm)，對于每個α的測試值，將所有Pm(u)相加； *選擇αopt以使∑Pm(u)達(dá)到最大值； -該方法還包括一種初始化模式，該初始化模式還包括在表面上確定N個傳感器位置的步驟； -互相關(guān)乘積Pij(ω)的個數(shù)等于 -Pk(u)的特征值是下列參數(shù)中的任一個 *Pk(u)的最大值； *Pk(u)的功率，或者最大幅值平方； *最大的峰-峰值； *Pk(u)的均方根；或者 *Pk(u)的能量值，其等同于 *Pk(u)的脈沖寬度參數(shù)，其等于其中REAL(x)是復(fù)數(shù)x的實部； -只有當(dāng)Pk0(u)的特征值大于預(yù)先設(shè)定的可靠閾值時沖擊才出現(xiàn)在區(qū)域k0； -只有在由比率f(Pk0(u))/MEAN(f(Pk，k≠k0(u)))定義的對比值大于預(yù)先設(shè)定的可靠閾值時，沖擊才出現(xiàn)在區(qū)域k0內(nèi)，其中MEAN是 一個平均值算子； -Pij(ω)被歸一化。
本發(fā)明的另一目的是提供一個用于確定物體表面沖擊位置的裝置，所述表面包括M個確定區(qū)域且所述沖擊產(chǎn)生一個聲學(xué)信號，所述裝置包括 -一個處理器； -N個位于所述表面的聲傳感器，其中N大于等于3，每個傳感 器i接收所述聲學(xué)信號并將感應(yīng)到的信號si(t)傳送到處理器，其中，所述處理器包括 (a)用于計算P個互相關(guān)乘積的裝置 其中， Si(ω)是傳感器i所感應(yīng)到的信號si(t)的傅立葉變換； Sj(ω)是傳感器j所感應(yīng)到的信號sj(t)的傅立葉變換； *是復(fù)共軛算子 (b)用于計算所述Pij(ω)的P個傅立葉反變換pij′(u)的裝置； (c)用于對每個區(qū)域k，計算Pk(u)＝∑p′ij(u-τijk)及Pk(u)的一個特征值的裝置，其中τijk是基于區(qū)域k及傳感器i和j各自位置的確定值； (d)用于確定k0的裝置，其中，當(dāng)k≠k0時，Pk0(u)的特征值大于相應(yīng)的Pk(u)的特征值。
因此，通過使用感知信號的互相關(guān)乘積，在處理信號時就不需要確定時間初始值。因此，實現(xiàn)本發(fā)明無需聲傳感器的準(zhǔn)確同步。此外，如將要解釋的，互相關(guān)乘積Pij(ω)的相位不依賴于沖擊的波形。因此，沖擊的波形不改變該方法的準(zhǔn)確度。



圖1a和1b，如上所述，是基于現(xiàn)有技術(shù)的描述聲學(xué)信號形成的計時圖； 圖2為本發(fā)明所述裝置的示意圖； 圖3為描述圖2所述裝置的基本操作的流程圖； 圖4為描述圖3中處理步驟的操作的流程圖； 圖5為描述初始化步驟的操作的流程圖。

具體實施例方式 如圖2所示，根據(jù)本發(fā)明的用于確定沖擊位置的裝置3包括N個聲傳感器，至少包括3個如圖所示的傳感器2a，2b，2c。所述裝置還包括一個連接到聲傳感器的處理器4和一個連接到處理器的存儲裝置，其中隨機(jī)存儲器存儲變量，只讀存儲器存儲程序和常量。此外，聲傳感器是固定的，例如，在一個物體表面1上。表面1可以被分為M個感應(yīng)區(qū)6。
如圖3所示，在操作前，一個初始化模式S100(INIT.)需要被執(zhí)行。這個初始化模式將在下文中被描述。在步驟S101中，一個信號測試(IMPACT？)用于檢測沖擊。例如，如果一個由傳感器感應(yīng)的信號超過一個確定閾值，由N個傳感器傳送的感應(yīng)信號si(t)會被采樣并在步驟S102(SAMPL.)中被存儲到處理器4中。之后，處理步驟S103(PROCES.)被執(zhí)行。
如圖4所示的更詳盡的描述，處理步驟從步驟S200的互相關(guān)乘積計算(INTERCORR.)開始。事實上，P個獨立的互相關(guān)乘積被計算。此處P等于1到N-1之間所有整數(shù)的和，例如，其中，N是傳感器的個數(shù)。每個互相關(guān)乘積Pij(ω)等于其中Si(ω)是傳感器i所感應(yīng)到的信號si(t)的傅立葉變換，Sj(ω)是傳感器j所感應(yīng)到的信號sj(t)的傅立葉變換，*是復(fù)共軛算子。所述傅立葉變換為，例如，只需要低運(yùn)算能力的快速傅立葉變化(FFT)。根據(jù)一個實施例，互相關(guān)乘積，例如在能量，被歸一化。
互相關(guān)乘積Pij(ω)不嚴(yán)格依賴于時間起始點和沖擊波形。事實上，如果考慮Si(ω)＝Ci(ω)·exp(-jkdi)·E(ω)，其中Ci是傳感器i的復(fù)頻響應(yīng)，di是沖擊位置和傳感器i之間的距離，E是沖擊波形的傅立葉變化。在本例中，波數(shù)k＝ω/C，其中C是聲傳播速率。那么，互相關(guān)乘積等于 或 當(dāng)E(ω)·E*(ω)是實數(shù)，時間互相關(guān)乘積pij(t)主要依賴于傳感器的響應(yīng)。一個大的傳感器響應(yīng)帶寬允許獲得非常短的脈沖信號，受限于沖擊激勵的頻率帶寬。在本例中，傳感器帶寬基本在[100Hz，7kHz]范圍內(nèi)。此外，如果互相關(guān)乘積是歸一化的，它就不依賴于感應(yīng)信號的幅度，和沖擊強(qiáng)度。
然后在步驟S201中，Pij(ω)的傅立葉反變換被計算。根據(jù)第一實施例，物體材料是不色散的。這說明聲學(xué)信號的傳播速率不依賴于這個信號的頻率。隨后，感應(yīng)信號在這種材料中的傳播波形基本等同于原始信號經(jīng)過一個延時，該延時是材料速率以及源和傳感器之間的距離的函數(shù)。在這種情況下，傅立葉反變換可能是一個常用的傅立葉反變換或者快速傅立葉反變換。
在步驟S202中，k被初始化為1，其目的是對表面上每個區(qū)域k執(zhí)行Pk(u)＝∑pij′(u-τijk)的計算(S203)，其中pij′(u)是Pij(ω)的傅立葉反變換。在現(xiàn)有的非色散材料的例子中，τijk是基于區(qū)域k至傳感器i的距離以及區(qū)域k至傳感器j的距離之間的差值計算出的延遲。這些值可以被存儲在設(shè)備1的存儲器5中。它們可以由操作人員在初始化模式S100中計算，或者在設(shè)計裝置時確定。
這樣，如果沖擊發(fā)生在區(qū)域k中，則所有的pij′(u-τijk)是同相的，Pk(u)＝∑pij′(u-τijk)實質(zhì)上是一個脈沖。相反，如果沖擊未發(fā)生于區(qū)域k中，則Pk(u)＝∑pij′(u-τijk)將會很寬而且幅度很小。這使得確定沖擊產(chǎn)生的位置更加容易。
對所有區(qū)域k的這種計算是通過步驟S204和S205的循環(huán)實現(xiàn)的。在步驟S206中，計算每個Pk(u)的特征值以及比較所有的特征值f(Pk(u))。這個特征值可以是 *Pk(u)的最大值； *Pk(u)的功率，或者最大幅值平方； *最大峰峰值； * Pk(u)的均方根； * Pk(u)的能量，其等于 *Pk(u)的脈沖寬度參數(shù)，其等于其中REAL(x)是復(fù)數(shù)x的實部； 那么，區(qū)域k0對應(yīng)于具有最大特征值f(Pk0(u))的區(qū)域。然后，在步驟S207中，最大值將與一個可靠閾值作比較，以保證檢測到的沖擊在區(qū)域k0內(nèi)。事實上，沖擊可能是一個干擾，由噪聲引起，或者是由在M區(qū)域以外的另一個區(qū)域中產(chǎn)生的。所以，這個最后的測試避免了裝置引起的非期望操作。
此外，這個測試可以被一個由比率f(Pk0(u))/MEAN(f(Pk，k≠k0(u)))定義的對比測試取代或合并，其中MEAN是除去k＝k0的所有f(Pk(u))值的平均值算子。如果對比值大于預(yù)設(shè)閾值時，結(jié)果被認(rèn)為是有效的。繼而，如果區(qū)域k0的特征值被認(rèn)為是有效的，則認(rèn)為沖擊發(fā)生在于這個區(qū)域中。裝置3就驅(qū)動一個操作，比如將這個裝置當(dāng)作接口使用。否則，裝置將忽略該沖擊。
此外，也可能用一個插值函數(shù)，例如來調(diào)整結(jié)果，允許對本方法的解決方案作改進(jìn)。
這個實施例提供了一種不需要精確同步來確定表面上沖擊位置的方法，所述方法尤其適用于非色散材料。
另外，本發(fā)明提供的第二個實施例針對于色散材料。本實施例與上第一述實施例的主要區(qū)別在于步驟S201中的傅立葉反變換是一個色散補(bǔ)償傅立葉變換。事實上，在色散材料中傳播時，聲學(xué)信號的傳播速率依賴于聲學(xué)信號的頻率。所以，感應(yīng)信號在這種材料中傳播時的波形是不同于原始信號波形的。在這種情況下，色散補(bǔ)償傅立葉變換，如下文所述，必須補(bǔ)償信號在傳播中的變形。
根據(jù)P.D.Wilcox&al在定量非破壞性評估的進(jìn)步回顧的學(xué)報上發(fā)表的標(biāo)題為“一種從導(dǎo)波信號中去除色散效果的信號處理技術(shù)”的文章，可以依靠以下公式的幫助來確定感應(yīng)到的聲學(xué)信號的信源距離 其中G(ω)是感應(yīng)到的聲學(xué)信號的傅立葉變換，k(ω)是取決于色散關(guān)系的波數(shù)。如果認(rèn)為k(ω)＝Ω且ω＝m(Ω)，可以從上述公式中導(dǎo)出 該公式是的傅立葉反變換。所獲信號的形狀被改變，但是延遲信息仍然保持相同，這是因為

通常是純實數(shù)。這樣，這種色散補(bǔ)償傅立葉反變換就可以用于改進(jìn)第一實施例，提供一種用于確定色散材料物體表面沖擊的方法。
所以，根據(jù)第二實施例的方法也基本相同，除了在步驟S201中需要計算的傅立葉反變換等于這樣，這個傅立葉變換補(bǔ)償了信號的色散，提供了可以被處理的信號。
在這個實施例中，步驟S203中，τijk值是基于區(qū)域k至傳感器i的距離以及區(qū)域k至傳感器j的距離之間的差值計算出的長度。本方法的其余步驟保持不變。
在大多數(shù)情況中，可以認(rèn)為一個基本平坦表面上的沖擊所產(chǎn)生的波形是蘭姆波。因為在本發(fā)明的本實施例中研究的波形是在低頻范圍內(nèi)的，可以認(rèn)為這些波形主要是零次蘭姆波。最后，當(dāng)表面上的沖擊發(fā)生時，所產(chǎn)生的波形為非對稱的零次蘭姆波。因此，在低頻范圍中，色散關(guān)系基本等于其中α是取決于材料的系數(shù)。
因此，就可以使用這種色散材料中的聲學(xué)圖像確定表面上沖擊的位置了。在滿足以下色散關(guān)系Ω＝αω時，第二實施例中的方法同樣也可以用在非色散物體表面。因此，這一實施例提供了一種可以適用于色散和非色散材料的通用方法。事實上，對于每個不同頻率范圍都可以很容易的改寫傳導(dǎo)波的色散關(guān)系。頻率范圍是由材料，應(yīng)用和傳感器的帶寬決定的。本實施例不僅限于當(dāng)前非對稱零次蘭姆波的例子。
此外，系數(shù)α可能依賴于材料溫度，這一參數(shù)也可能需要進(jìn)行考慮。例如，可以在物體上安裝一個溫度傳感器，系數(shù)α可能會根據(jù)這一參數(shù)有所變化。因此，可以提供一種考慮物體溫度的計算方法。同樣，如果材料是各相異性的，例如傳播速率依賴于信號的傳播方向，就可以被一個基于傳感器i和區(qū)域k之間的路徑的系數(shù)α所補(bǔ)償。例如，α依賴于θik，而θik是一個基于區(qū)域k和傳感器i的相對位置的確定的角度值。
根據(jù)圖3所示的流程圖，本方法由一個如圖5流程圖所示的初始化模式起始。在這一初始化模式中，可以通過步驟S300將傳感器的準(zhǔn)確位置輸入并存儲于裝置3中，但是這些傳感器的位置也可以在設(shè)計裝置的時候被輸入。
接下來，確定色散關(guān)系。通常，色散關(guān)系的算式是已知的。這個步驟的目的是為色散關(guān)系的差值計算提供正確的系數(shù)。在非對稱零次蘭姆波的情況下，我們需要確定系數(shù)α。但是，如果信號是在另一個頻率范圍內(nèi)被研究，可能需要一個色散關(guān)系的多項式插值運(yùn)算。那么所有的多項式系數(shù)都必須被確定。
此外，即使材料被認(rèn)為是一種非色散材料，這一步驟也可以被用來確定對應(yīng)于這一非色散材料的正確的系數(shù)α。為了完成對系數(shù)α的確定，可以在材料表面上生成NSO沖擊(步驟S301至S305)并且計算對應(yīng)于每個沖擊m的延遲τijm并將其保存在存儲器中。在這一實施例中，沖擊產(chǎn)生于兩個傳感器的連線上，且均勻分布。這樣很容易找到系數(shù)α的最優(yōu)值。
然后，α值將通過步驟S306至步驟S314在確定區(qū)間[αmin∶αmax]被測試。測試值α以增量Δα累加(步驟S313)。對于每個測試值，根據(jù)測試值α計算與沖擊m以及傳感器i和j有關(guān)的互相關(guān)乘積pijm′(u)(步驟S308)，并將其累加得到(步驟S309)。這樣，所有的沖擊Pm(u)被相加，這個和的最大值通過步驟S312被保存。因此，對于每個測試值α，有一個最大值F(α)被保存。
然后，在步驟S315中，最大值max(F(α))被認(rèn)為是最佳值，這樣相應(yīng)的α被選擇為最優(yōu)系數(shù)αopi并將被用于裝置的操作中。此外，單獨的二分法或者是其與上述方法的結(jié)合也可以被用于尋找最佳系數(shù)。
因此，可以提供一種通過充分材料色散補(bǔ)償和基于聲學(xué)圖像確定表面上沖擊位置的方法。
根據(jù)另一實施例，初始化模式包括一個確定M區(qū)域位置和相關(guān)延時的步驟。例如，一個操作員需要在每個敏感區(qū)生成沖擊，裝置將相關(guān)延時保存在它的只讀存儲器5中。
權(quán)利要求
1.一種用于確定物體表面(1)沖擊位置的方法，所述表面(1)包括N個聲學(xué)傳感器(2a，2b，2c)和M個確定的區(qū)域，并且N大于等于3，所述沖擊產(chǎn)生一個聲學(xué)信號，每個傳感器i接收所述聲學(xué)信號并且將感應(yīng)到的信號si(t)傳送到一個處理器(4)，所述方法包括以下步驟
(a)計算P個互相關(guān)乘積
其中，
Si(ω)是傳感器i所感應(yīng)到的信號si(t)的傅立葉變換；
Sj(ω)是傳感器j所感應(yīng)到的信號sj(t)的傅立葉變換；
*是復(fù)共軛算子；
(b)計算所述Pij(ω)的P個傅立葉反變換pij′(u)；
(c)對每個區(qū)域k，計算Pk(u)＝∑p′ij(u-τijk)及Pk(u)的一個特征值，其中τijk是一個基于區(qū)域k及傳感器i和j各自位置的確定值；
(d)確定k0，，其中，當(dāng)k≠k0時，Pk0(u)的特征值大于相應(yīng)的Pk(u)的特征值。
2.如權(quán)利要求1所述的方法，其特征在于，傅立葉反變換pij′(u)等于
其中Ω是波數(shù)系數(shù)；
m(Ω)是波數(shù)系數(shù)Ω的頻率，決定于材料的色散關(guān)系。
3.如權(quán)利要求2所述的方法，其特征在于，u是距離，τijk是一個基于區(qū)域k與傳感器i之間的距離以及區(qū)域k與傳感器j之間的距離的長度。
4.如權(quán)利要求3所述的方法，其特征在于，該材料色散關(guān)系等于其中α是取決于物體的一個系數(shù)。
5.如權(quán)利要求3所述的方法，其特征在于，該材料色散關(guān)系等于Ω＝αω，其中α是取決于物體的一個系數(shù)。
6.如權(quán)利要求4或5所述的方法，其特征在于，該系數(shù)α取決于θik，θik是基于區(qū)域k和傳感器i的位置的確定角度值。
7.如權(quán)利要求4至6中任一項所述的方法，其特征在于，該系數(shù)α取決于物體的溫度a(T)，其中，T是溫度值。
8.如權(quán)利要求2至7中任一項所述的方法，其特征在于，還包括一個初始化模式，所述初始化模式包括一個確定物體材料色散關(guān)系的步驟。
9.如權(quán)利要求8所述的方法，其特征在于，所述確定材料色散關(guān)系步驟包括在所確定位置產(chǎn)生NSO沖擊的子步驟。
10.如權(quán)利要求9所述的方法，其特征在于，所述NSO沖擊產(chǎn)生于N個傳感器中的兩個傳感器之間的連線上的確定位置，并且均勻分布。
11.如權(quán)利要求9或10所述方法，其特征在于，所述確定材料色散關(guān)系的步驟還包括以下子步驟
*用α測試值計算pij′(u)；
*對于每個α測試值，對于每個沖擊m，計算Pm(u)＝∑p′ij(u-τijm)，對于每個α測試值，將所有Pm(u)相加；
*選擇αopt以使∑Pm(u)達(dá)到最大值。
12.如上述權(quán)利要求中任一項所述的方法，其特征在于，還包括初始化模式，所述初始化模塊包括在表面上(1)確定N個傳感器(2a，2b，2c)位置的步驟。
13.如上述權(quán)利要求中任一項所述的方法，其特征在于，互相關(guān)乘積Pij(ω)的個數(shù)等于
14.如上述權(quán)利要求中任一項所述的方法，其特征在于，Pk(u)的特征值是下列參數(shù)中的任一個
-Pk(u)的最大值；
-Pk(u)的功率，或者最大幅值平方；
-最大峰峰值；
-Pk(u)的均方根；
-Pk(u)的能量，其等于
-Pk(u)的脈沖寬度參數(shù)，其等于，其中
REAL(x)是復(fù)數(shù)x的實部。
15.如上述權(quán)利要求中任一項所述的方法，其特征在于，只有在Pk0(u)的特征值大于預(yù)先設(shè)定的可靠閾值時，沖擊才出現(xiàn)在區(qū)域k0內(nèi)。
16.如上述權(quán)利要求中任一項所述的方法，其特征在于，只有在由比率f(Pk0(u))/MEAN(f(Pk，k≠k0(u)))定義的對比值大于預(yù)先設(shè)定的可靠閾值時，沖擊才出現(xiàn)在區(qū)域k0內(nèi)，其中MEAN是一個平均值算子。
17.如上述權(quán)利要求中任一項所述的方法，其特征在于，Pij(ω)被歸一化。
18.一種用于確定物體表面(1)上沖擊位置的裝置，所述表面(1)包括M個確定區(qū)域(6)，所述沖擊產(chǎn)生聲學(xué)信號，所述裝置包括
-一個處理器(4)；
-N個位于所述表面(1)的聲傳感器(2a，2b，2c)，其中N大于等于3，每個傳感器i接收所述聲學(xué)信號并將感應(yīng)到的信號si(t)傳送到處理器(4)，其中，所述處理器(4)包括
(a)用于計算P個互相關(guān)乘積的裝置
其中，
Si(ω)是傳感器i所感應(yīng)到的信號si(t)的傅立葉變換；
Sj(ω)是傳感器j所感應(yīng)到的信號sj(t)的傅立葉變換；
*是復(fù)共軛算子；
(b)用于計算所述Pij(ω)的P個傅立葉反變換pij′(u)的裝置；
(c)用于對每個區(qū)域k，計算Pk(u)＝∑p′ij(u-τijk)及Pk(u)的一個特征值的裝置，其中τijk是基于區(qū)域k及傳感器i和j各自位置的確定值；
(d)用于確定k0的裝置，其中，當(dāng)k≠k0時，Pk0(u)的特征值大于相應(yīng)的Pk(u)的特征值。
全文摘要
本發(fā)明提供了一種確定表面(1)上沖擊位置的方法，表面(1)上包括N個聲學(xué)傳感器(2a，2b，2c)，其用于發(fā)送一個感應(yīng)到的信號si(t)至處理器(4)。該方法包括以下步驟(a)計算P個互相關(guān)乘積；(b)計算P個傅立葉反變換值pij′(u)；(c)對于每個區(qū)域k，計算Pk(u)＝∑p′ij(u-τijk)；(d)尋找k0，其中，Pk0(u)的特征值大于一個給定的閾值。
文檔編號G01N29/14GK101198861SQ200580050098
公開日2008年6月11日 申請日期2005年4月13日 優(yōu)先權(quán)日2005年4月13日
發(fā)明者羅斯·基里·英 申請人:傳感器公司