Method for determining the position of impacts
View Patent ↗A method for determining the position of impacts on an object comprising two acoustic sensors, and N active areas of said object, comprises the steps of: (a) receiving two acoustic signals S 1 (t) and S 2 (t); (b) calculating a sample signature function SIG S (ω)=S 1 (ω)−S 2 (ω)*, where S 1 (ω) and S 2 (ω) are the respective Fourier transforms of S 1 (t) and S2 (t), (c) comparing SIG S (ω) with N predetermined reference signature functions SIG R, (ω) corresponding to the predetermined area j for j from 1 to N; (d) determining the active area in which the impact occurred, on the basis of the comparison of step (c).
1. A method for determining the position of impacts on an object, said object comprising:
two acoustic sensor;
N predetermined active areas, whereby N is an integer at least equal to 1;
said method comprising the steps of:
(a) receiving two acoustic signals S 1 (t) and S 2 (t) originating respectively from said acoustic sensors and generated by an impact received on said object;
(b) calculating, using at least one processing unit, a sample signature function:
SIG S (ω)= S 1 (ω)· S 2 (ω)*,
where S 1 (ω) and S 2 (ω) are the respective Fourier transforms of S 1 (t) and S 2 (t), and where * is the complex conjugate operator;
(c) comparing SIG S (ω) with N predetermined reference signature functions SIG Rj (ω) each corresponding to a predetermined active area j, for j from 1 to N;
(d) determining an active area in which the impact occurred, on the basis of the comparison of step (c).
2. The method as claimed in claim 1 , wherein each reference signature function equals:
SIG Rj (ω)= R 1j (ω)· R 2j (ω)*
where R 1j (ω) and R 2j (ω) are Fourier transforms of acoustic signals r 1j (t) and r 2j (t) received by each of the respective acoustic sensors when an impact occurs on the predetermined area j.
3. The method as claimed in claim 1 , wherein step (c) comprises the calculation of a similarity estimator α j representing a function of the phase φ(COR j (ω)) of COR j (ω)=SIG R j (ω)·SIG S (ω)*.
4. The method as claimed in claim 3 , wherein step (c) comprises the calculation of a function δ j (ω) for j from 1 to N, wherein
δ j (ω)=ε k , if φ(COR j (ω)) belongs to I k ,
where ε k is a predetermined value and I k a corresponding angular interval for k from 1 to n, where n is an integer greater than 1.
5. The method as claimed in claim 4 , wherein the values ε i are not greater than 1.
6. The method as claimed in claim 4 or claim 5 , wherein:
if |φ(COR j (ω))|≦a1, then δ j (ω)=ε 1 ,
. . .
if n is greater than 2 and a k-1 <|φ(COR j (ω))|≦a k , then δ j (ω)=ε k , for k=2 . . . n−1,
. . .
if |φ(COR j (ω))|>a n−1 , then δ j (ω)=ε n ,
wherein a k is increasing with k and ε k is decreasing with k.
7. The method as claimed in claim 6 , wherein:
if Re(COR j (ω))≧|Im(COR j (ω))|/tan(a 1 ), then δ j (ω)=ε 1 ,
. . .
if n is greater than 2 and |Im(COR j (ω))|/tan(a k-1 )>Re(COR j (ω))≧Im(COR j (ω))|/tan(a k ), then δ j (ω)=ε k , for k=2 . . . n−1,
. . .
else, δ j (ω)=ε n ,
and wherein Re(COR j (ω)) is the real part of COR j (ω) and Im(COR j (ω)) is the imaginary part of COR j (ω).
8. The method according to claim 6 , wherein:
if γ j (ω)≧β j (ω)/tan(a 1 ), then δ j (ω)=ε 1 =1;
. . .
if n is greater than 2 and β j (ω)/tan(a k-1 )>γ j (ω)≧β j (ω)/tan(a k ), then δ j (ω)=ε k , for k=2, . . . n−1;
. . .
else, δ j (ω)=ε n ,
and wherein:
β j (ω)=|1+{A j (ω)/B j (ω)}·{(D(ω)/C(ω)}|,
γ j (ω)=sign B j (ω)·sign C(ω)·[{A j (ω)/B j (ω)}−{D(ω)/C(ω)}],
sign B j (ω)=1 if B j (ω) is positive and =−1 if B j (ω) is negative,
sign C(ω)=1 if C(ω) is positive and =−1 if C(ω) is negative,
A j (ω) and B j (ω) are respectively the real part and the imaginary part of each reference signature function SIG Rj (ω),
C(ω) and D(ω) are respectively the real part and the imaginary part of the complex conjugate sample signature function SIG S (ω)*.
9. The method as claimed in claim 6 , wherein ε 1 =1 and ε n =0.
10. The method as claimed in claim 9 , wherein n is greater than 2 and ε k =cos(a k-1 ), for k=2, . . . n−1.
11. The method as claimed in claim 4 , wherein
α j =K·∫ B δ j (ω)· dω,
where B is a frequency interval and K a constant.
12. The method as claimed in claim 11 , wherein B=[ω min ,ω max ] and K is proportional to 1/(ω max −ω min ).
13. The method as claimed in claim 3 , wherein the active area j 0 in which the impact occurred is determined such that α j 0 is the greatest similarity estimator among the N calculated similarity estimators α j .
14. The method as claimed in claim 13 , wherein it is concluded that the impact occurred in the area j 0 only if α j 0 is considered as valid in step (d).
15. The method of claim 14 , wherein α j 0 is considered as valid if α j 0 is greater than a predetermined threshold of confidence.
16. The method as claimed in claim 14 , wherein it is concluded that there is no impact if α j 0 is not considered as valid in step (d).
17. The method as claimed in claim 1 , wherein step (a) begins if the acoustic signals s 1 (t) and S 2 (t) are above a predetermined trigger threshold.
18. The method as claimed in claim 1 , wherein the acoustic signals s 1 (t) and S 2 (t) in step (a) are delayed relative to real audio signals.
19. The method as claimed in claim 1 , wherein the Fourier transform is a Fast Fourier transform.
20. A device for determining the position of impacts on an object, comprising:
two acoustic sensors adapted to be borne by said object for receiving acoustic signals S 1 (t) and s 2 (t) generated by an impact on said object;
memory means comprising N reference signature functions, corresponding to N predetermined active areas of said object, whereby N is an integer at least equal to 1;
calculation means for calculating a sample signature function
SIG S (ω)= S 1 (ω)· S 2 (ω)*,
where S 1 (ω) and S 2 (ω) are the respective Fourier transforms of S 1 (t) and S 2 (t), and where * is the complex conjugate operator;
comparison means for comparing SIG S (ω) with N predetermined reference signature functions SIG R j (ω) for j from 1 to N;
processing means for determining an active area in which the impact occurred, on the basis of results given by said the comparison means.
21. The device as claimed in claim 20 , wherein each reference signature function equals:
SIG R j (ω)= R 1j (ω)· R 2j (ω)*
where R 1j (ω) and R 2j (ω) are Fourier transforms of acoustic signals r 1j (t) and r 2j (t) received by each of the respective acoustic sensors when an impact occurs on the predetermined area j.
22. The device as claimed in claim 20 , wherein said comparison means are adapted to calculate a similarity estimator α j representing a function of the phase φ(COR j (ω)) of COR j (ω)=SIG R j (ω)·SIG S (ω)*.
23. The device as claimed in claim 22 , wherein said comparison means are adapted to calculate a function δ j (ω) for j from 1 to N, wherein
δ j (ω)=ε k , if φ(COR j (ω)) belongs to I k , where ε k is a predetermined value and I k a corresponding angular interval for k from 1 to n, where n is an integer greater than 1.
24. The device as claimed in claim 23 , wherein the values S 1 are not greater than 1.
25. The device as claimed in claim 23 , wherein said comparison means are adapted so that:
if |φ(COR j (ω))|≦a1, then δ j (ω)=ε 1 ,
. . .
if n is greater than 2 and a k-1 <|φ(COR j (ω))|≦a k , then δ j (ω)=ε k , for k=2 . . . n−1,
. . .
if |φ(COR j (ω))|>a n−1 , then δ j (ω)=ε n ,
wherein a k is increasing with k and ε k is decreasing with k.
26. The device as claimed in claim 25 , wherein said comparison means are adapted so that:
if Re(COR j (ω))≧|Im(COR j (ω))|/tan(a 1 ), then δ j (ω)=ε 1 ,
. . .
if n is greater than 2 and |Im(COR j (ω))|/tan(a k-1 )>Re(COR j (ω))≧|Im(COR j (ω))|/tan(a k ), then δ 1 (ω)=ε k , for k=2 . . . n−1,
. . .
else, δ j (ω)=ε n ,
and wherein Re(COR j (ω)) is the real part of COR j (ω) and Im(COR j (ω)) is the imaginary part of COR j (ω).
27. The device according to claim 25 , wherein said comparison means are adapted so that:
if γ j (ω)≧β j (ω)/tan(a 1 ), then δ j (ω)=ε 1 =1;
. . .
if n is greater than 2 and β j (ω)/tan(a k-1 )>γ j (ω)≧β j (ω)/tan(a k ), then δ j (ω)=ε k , for k=2, . . . n−1;
. . .
else, δ j (ω)=ε n ,
and wherein:
β j (ω)=|1+{A j (ω)/B j (ω)}·{(D(ω)/C(ω)}|,
γ j (ω)=sign B j (ω)·sign C(ω)·[{A j (ω)/B j (ω)}−{D(ω)/C(ω)}],
sign B j (ω)=1 if B j (ω) is positive and =−1 if B j (ω) is negative,
sign C(ω)=1 if C(ω) is positive and =−1 if C(ω) is negative,
A j (ω) and B j (ω) are respectively the real part and the imaginary part of each reference signature function SIG Rj (ω),
C(ω) and D(ω) are respectively the real part and the imaginary part of the complex conjugate sample signature function SIG S (ω)*.
28. The device as claimed in claim 25 , wherein ε 1 =1 and ε n =0.
29. The device as claimed in claim 28 , wherein n is greater than 2 and ε k =cos(a k-1 ), for k=2, . . . n−1.
30. The device as claimed in claim 23 , wherein said comparison means are adapted so that:
α j =K·∫ B δ j (ω)· dω,
where B is a frequency interval and K a constant.
31. The device as claimed in claim 22 , wherein the processing means are adapted to determine the active area j 0 in which the impact occurred such that α j 0 is the greatest similarity estimator among the N calculated similarity estimators α j .
32. The device as claimed in claim 31 , wherein the processing means are adapted to determine that the impact occurred in the area j 0 only if α j 0 is considered as valid.
33. The device of claim 32 , wherein the processing means are adapted to determine that α j 0 is valid if α j 0 is greater than a predetermined threshold of confidence.
34. The device as claimed in claim 32 , wherein the processing means are adapted to determine that there is no impact if α j 0 is not considered as valid.