IPC분류정보
국가/구분 |
United States(US) Patent
등록
|
국제특허분류(IPC7판) |
|
출원번호 |
US-0092252
(2011-04-22)
|
등록번호 |
US-8861311
(2014-10-14)
|
발명자
/ 주소 |
- Ledeczi, Akos
- Volgyesi, Peter
- Sallai, Janos
|
출원인 / 주소 |
|
대리인 / 주소 |
Novak Druce Connolly Bove + Quigg LLP
|
인용정보 |
피인용 횟수 :
0 인용 특허 :
6 |
초록
▼
Systems and methods for estimating projectile trajectory and projectile source location are provided. A method for estimating location information associated with a supersonic projectile propelled from a source includes recording sound at a first location using a single microphone during travel of t
Systems and methods for estimating projectile trajectory and projectile source location are provided. A method for estimating location information associated with a supersonic projectile propelled from a source includes recording sound at a first location using a single microphone during travel of the supersonic projectile to produce an acoustic recording. The method further includes estimating a miss distance between the first location and a trajectory of the projectile based on the shockwave length. Locating a projectile source includes concurrently recording sound at multiple locations and generating data sets associated with the locations, each of the plurality of data sets including a miss distance, a range, a time of arrival of a muzzle blast from the source, and a time of arrival of a shockwave produced by the projectile. Additionally, the method includes calculating an approximate location of the source at each of the locations based on the data sets.
대표청구항
▼
1. A method for estimating location information associated with a supersonic projectile propelled from a source, the method comprising: recording sound at a first location using a single microphone during travel of the supersonic projectile to produce an acoustic recording; andestimating a miss dist
1. A method for estimating location information associated with a supersonic projectile propelled from a source, the method comprising: recording sound at a first location using a single microphone during travel of the supersonic projectile to produce an acoustic recording; andestimating a miss distance (b) between the first location and a trajectory of the supersonic projectile based on a length of a shockwave generated by the supersonic projectile (shockwave length),wherein the estimating comprises identifying the shockwave length based on the acoustic recording and determining the miss distance based on the shockwave length and a plurality of projectile characteristics, andwherein the determining of the miss distance comprises performing one of: selecting the miss distance based on a table of predefined miss distance values as a function of a period of the shockwave based on the acoustic recording, orevaluating: b=l(TcdM(1.82))4(M2-1)3/2where l is the projectile length, d is the projectile diameter, T is the shockwave length, c is the speed of sound, and M is the Mach number associated with the projectile. 2. The method of claim 1, wherein said recording comprises providing mechanical damping of the sound at the first location. 3. The method of claim 1, wherein said recording further comprises sampling the sound at a sampling rate of at least 48 kHz. 4. The method of claim 1, further comprising: calculating a time difference of arrival (TDOA) between a muzzle blast at the source and the shockwave generated by the supersonic projectile based on the acoustic recording;estimating a range between the first location and the source based on the miss distance and the TDOA. 5. The method of claim 4, wherein said estimating said range further comprises: computing an angle between the trajectory and a surface of the shockwave; andcalculating the range based on the TDOA, the miss distance, and the angle. 6. The method of claim 5, wherein said calculating said range (dS,M) further comprises evaluating: dS,M=12(c4-v4)(P-2Q)where P and Q are defined as P=-2v3dQ,Mv2+c2-2(tmb-tsw)c3v2+2c2dQ,Mvv2+c2-2(tmb-tsw)cv4Q=-2c4v4dQ,M2+2(tmb-tsw)2c6v4+(tmb-tsw)2c4v6-2c7dQ,M(tmb-tsw)vv2+c2+c8(tmb-tsw)2v2+2c8dQ,M2+2v5dQ,Mv2+c2(tmb-tsw)c3,and where tmb−tsw is the TDOA, dQ,M is the miss distance, v is the velocity of the projectile, and c is the speed of sound. 7. A method for locating a source of a projectile, comprising: concurrently recording sound at three or more locations during travel of the supersonic projectile to produce acoustic recordings at each of said locations;generating a plurality of data sets associated with said locations, each of said plurality of data sets comprising a miss distance, a range, a time of arrival of a muzzle blast from said source, and a time of arrival of a shockwave produced by said projectile;sharing said plurality of data sets among said plurality of locations; andindependently calculating an approximate location of said source at each of said locations. 8. The method of claim 7, wherein said each of said plurality of data sets further comprises time and location information. 9. The method of claim 7, wherein said recording comprises providing mechanical damping of the sound at said locations. 10. The method of claim 7, wherein said recording further comprises sampling the sound at a sampling rate of at least 48 kHz. 11. The method of claim 7, wherein said generating further comprises: identifying a length of the shockwave (shockwave length) based on the acoustic recording;determining the miss distance based on the shockwave length and a plurality of projectile characteristics. 12. The method of claim 11, wherein determining the miss distance (b) comprises evaluating: b=l(TcdM(1.82))4(M2-1)3/2where l is the projectile length, d is the projectile diameter, T is the shockwave length, c is the speed of sound, and M is the Mach number associated with the projectile. 13. The method of claim 11, wherein determining the miss distance comprises selecting the miss distance based on a table of predefined miss distance values as a function of a period of the shockwave. 14. The method of claim 7, wherein said generating further comprises: computing an angle between the trajectory and a surface of the shockwave; andcalculating the range based on a time difference of arrival between the muzzle blast and the shockwave (TDOA), the miss distance, and the angle. 15. The method of claim 14, wherein said calculating said range (dS,M) further comprises evaluating: dS,M=12(c4-v4)(P-2Q)where P and Q are defined as P=-2v3dQ,Mv2+c2-2(tmb-tsw)c3v2+2c2dQ,Mvv2+c2-2(tmb-tsw)cv4Q=-2c4v4dQ,M2+2(tmb-tsw)2c6v4+(tmb-tsw)2c4v6-2c7dQ,M(tmb-tsw)vv2+c2+c8(tmb-tsw)2v2+2c8dQ,M2+2v5dQ,Mv2+c2(tmb-tsw)c3,and where tmb−tsw is the TDOA, dQ,M is the miss distance, v is the velocity of the projectile, and c is the speed of sound. 16. The method of claim 7, wherein said independently calculating at each of said locations further comprises performing a sensor fusion process based on a time difference of arrival between the muzzle blast and the shockwave (TDOA), the miss distance, and the range computed for each of said locations. 17. A system for determining location information associated with a supersonic projectile, comprising at least one sensor device, the sensor device comprising: a transducer system configured for sampling sound at a first location during travel of the supersonic projectile to produce a single channel acoustic recording of said sound;a storage element for storing the acoustic recording; anda processing element configured for calculating a time difference of arrival (TDOA) between a muzzle blast at the source and a shockwave generated by the supersonic projectile based on the acoustic recording, estimating a miss distance (b) between the first location and estimating a range between the first location and the source based on the miss distance and the TDOA,wherein the processing element is configured for estimating the miss distance by identifying the shockwave length based on the acoustic recording and determining the miss distance based on the shockwave length and a plurality of projectile characteristics, andwherein the determining of the miss distance comprises performing one of: selecting the miss distance based on a table of predefined miss distance values as a function of a period of the shockwave based on the acoustic recording, orevaluating: b=l(TcdM(1.82))4(M4-1)3/2.where l is the projectile length, d is the projectile diameter, T is the shockwave length, c is the speed of sound, and M is the Mach number associated with the projectile. 18. The system of claim 17, wherein said transducer system comprises a microphone and a mechanical damper disposed over said microphone. 19. The system of claim 17, wherein said transducer system is configured for sampling the sound at a sampling rate of at least 48 kHz. 20. The system of claim 17, wherein said processing element is further configured for said estimating of said range by: computing an angle between the trajectory and a surface of the shockwave; andcalculating the range based on the TDOA, the miss distance, and the angle. 21. The method of claim 17, wherein processing element is further configured for said calculating said range (dS,M) by evaluating: dS,M=12(c4-v4)(P-2Q)where P and Q are defined as P=-2v3dQ,Mv2+c2-2(tmb-tsw)c3v2+2c2dQ,Mvv2+c2-2(tmb-tsw)cv4Q=-2c4v4dQ,M2+2(tmb-tsw)2c6v4+(tmb-tsw)2c4v6-2c7dQ,M(tmb-tsw)vv2+c2+c8(tmb-tsw)2v2+2c8dQ,M2+2v5dQ,Mv2+c2(tmb-tsw)c3,and where tmb−tsw is the TDOA, dQ,M is the miss distance, v is the velocity of the projectile, and c is the speed of sound. 22. The system of claim 17, wherein said at least one sensor device comprises three or more the sensor devices at a plurality of sensor locations, and wherein each of said sensor devices further comprises at least one transceiver, wherein said transceiver is configured for sharing a time of arrival (TOA) of the shockwave, a TOA of the muzzle blast, said miss distance, and said range among said sensor devices, and wherein said processing element is further configured for estimating a location of said source and a trajectory of said projectile based on said TOA of said shockwave, said TOA of the muzzle blast, said miss distance, and said range from said three or more sensor devices. 23. The system of claim 22, wherein the processing element is further configured for estimating said location by: determining a consistency function for a plurality of grid locations based on said TOA of said shockwave, said TOA of the muzzle blast, said miss distance, and said range from said three or more sensor devices;selecting a one of said plurality of grid locations associated with a maximum value for said consistency function. 24. The system of claim 23, wherein determining the consistency function comprises generating miss distance circles and range circles for each of said plurality of sensor locations based on said miss distance and said range associated with said three or more sensors;determining a number of miss distance circles whose tangent is within a pre-defined miss distance error angle; anddetermining a number of range circles intersecting each of said plurality of grids. 25. The system of claim 24, wherein determining the consistency function further comprises: identifying a number of said TOA of said muzzle blast range from said three or more sensor devices for each of said plurality of grids; andidentifying a number of said TOA of said shockwave from said three or more sensor devices for each of said plurality of grids. 26. A method for locating a source of a projectile, comprising: obtaining, for a plurality of sensor locations, a length of a shockwave associated with the projectile (shockwave lengths), times of arrival for the shockwave (shockwave TOAs), and time of arrival for a muzzle blast associated with the projectile (muzzle blast TOAs);estimating a trajectory for the projectile based at least one the plurality of sensor locations, the shockwave lengths, and the shockwave TOAs, the estimating comprising providing an error function relating miss distance to shockwave length and calculating the trajectory using a set of miss distance values for the plurality of locations that minimize the error function;determining a final source location for the source based on the trajectory and the plurality of sensor locations. 27. The method of claim 26, wherein the error function comprises: ∑wi(bi-kTi4)2,i=1…n,where n is the number of the plurality of sensor locations, bi is the distance of ith one of the plurality of locations from the trajectory, and wiε{0, 1) are binary weights assigned such that only the N smallest error values are considered. 28. The method of claim 26, wherein the error function comprises: ∑wi(bi-kTi4)2var({αi,j}),i,j=1…nwhere n is the number of the plurality of sensor locations, bi is the distance of ith one of the plurality of locations from the trajectory, αi,j is a shockwave angle computed from a time delay of arrival (TDOA) between the muzzle blast and the shockwave between the ith one of the plurality of locations and a jth one of the plurality of locations, and wiε{0, 1) are binary weights assigned such that only the N smallest error values are considered. 29. The method of claim 26, further comprising: computing a scaling factor based on the set of miss distances that minimize the error function;computing a velocity of the projectile based on the shockwave angle computed from time delay of arrival (TDOA) between the muzzle blast and the shockwave between a first of the plurality of locations and a second of the plurality of locations; andclassifying the projectile using a projectile coefficient derived from the scaling factor and the velocity of the projectile. 30. The method of claim 29, further comprising: computing a muzzle velocity for a source of the projectile based on the velocity of the projectile; andclassifying the weapon using the muzzle velocity. 31. The method of claim 30, wherein computing the muzzle velocity (vmuzzle) further comprises evaluating for at least one of the plurality of locations: vmuzzle=vprojectile2-2arwhere vprojectile is the velocity of the projectile, a is a deceleration factor, and r is an estimated range from the one of the plurality of locations to the source of the projectile. 32. The method of claim 26, wherein the step of determining further comprises: providing an error function relating a Euclidean distance from a one of the plurality of sensor locations to a candidate location for the source of the projectile and a range from the one of the plurality of sensor locations to the source of the projectile; andcalculating the final source location using the set of Euclidean distance values for the plurality of sensor location that minimize the error function. 33. The method of claim 32, wherein the error function comprises: ∑wi(ri-dS,Mi)2,i=1…nwhere n is the number of the plurality of locations, ri is the range from an ith one of the plurality of locations to the candidate location, and dS,Mi is the Euclidean distance for between the ith one of the plurality of locations and the candidate location, and wiε{0, 1) are binary weights assigned such that only the N smallest error values are considered. 34. The method of claim 33, wherein said calculating the range if the ith one of the plurality of locations further comprises evaluating: r=12(c4-υ4)(P-2Q)where P and Q are defined as P=-2v3dQ,Mv2+c2-2(tmb-tsw)c3v2+2c2dQ,Mvv2+c2-2(tmb-tsw)cv4Q=-2c4v4dQ,M2+2(tmb-tsw)2c6v4+(tmb-tsw)2c4v6-2c7dQ,M(tmb-tsw)vv2+c2+c8(tmb-tsw)2v2+2c8dQ,M2+2v5dQ,Mv2+c2(tmb-tsw)c3,and where tmb−tsw is a time difference of arrival (TDOA) between the muzzle blast and the shockwave for the ith one of the plurality of locations, dQ,M is the miss distance for the ith one of the plurality of locations, v is the velocity of the projectile based on the shockwave angle, and c is the speed of sound. 35. The method of claim 26, further comprising: computing an initial source location for the source based on multilateration using the muzzle blast TOAs; andperforming the step of estimating based on the initial source location. 36. A system for locating a source of a projectile, comprising: a processor;a first module to configure the processor to obtain, for plurality of sensor locations, a length of a shockwave associated with the projectile (shockwave lengths), times of arrival (TOAs) for the shockwave (shockwave TOAs), and TOAs for a muzzle blast associated with the projectile (muzzle blast TOAs);a second module to configure the processor to estimate a trajectory for the projectile based at least the plurality of sensor locations, the shockwave lengths, and the shockwave TOAs;a third module to configure the processor to determine a final source location for the source based on the trajectory and the plurality of sensor locations. 37. The system of claim 36, wherein the first module is configured for: subdividing a search space into a plurality of points;for each of the plurality of points, calculating shot times for the plurality of sensor locations based on the muzzle blast TOAs; andselecting as the initial source location a one of the plurality of points associated with comparable the shot times among the plurality of sensor locations. 38. The system of claim 37, wherein the second module is configured for: obtaining an error function relating miss distance to shockwave length; andcalculating the trajectory using a set of miss distance values for the plurality of locations that minimize the error function. 39. The system of claim 38, wherein the error function comprises: ∑wi(bi-kTi4)2,i=1…n,where n is the number of the plurality of sensor locations, bi is the distance of ith one of the plurality of locations from the trajectory, and wiε{0, 1) are binary weights assigned such that only the N smallest error values are considered. 40. The system of claim 38, wherein the error function comprises: ∑wi(bi-kTi4)2var({αi,j}),i,j=1…nwhere n is the number of the plurality of sensor locations, bi is the distance of ith one of the plurality of locations from the trajectory, αi,j is a shockwave angle computed from the shockwave time delay of arrival (TDOA) between the ith one of the plurality of locations and a jth one of the plurality of locations, and wiε{0, 1) are binary weights assigned such that only the N smallest error values are considered. 41. The system of claim 38, further comprising a fourth module to configure the processor to: compute a scaling factor based on the set of miss distances that minimize the error function;compute a velocity of the projectile based on the shockwave angle computed from a time delay of arrival (TDOA) between the muzzle blast and the shockwave between a first of the plurality of locations and a second of the plurality of locations; andclassify the projectile using a projectile coefficient derived from the scaling factor and the velocity of the projectile. 42. The system of claim 41, further comprising a fifth module to configure the processor to: computing a muzzle velocity for a source of the projectile based on the velocity of the projectile; andclassifying the weapon using the muzzle velocity. 43. The system of claim 42, wherein the computing the muzzle velocity (vmuzzle) further comprises evaluating for at least one of the plurality of locations: vmuzzle=vprojectile2-2arwhere vprojectile is the velocity of the projectile, a is a deceleration factor, and r is an estimated range from the one of the plurality of locations to the source of the projectile. 44. The system of claim 36, wherein the step of determining further comprises: providing an error function relating a Euclidean distance from a one of the plurality of sensor locations to a candidate location for the source of the projectile and a range from the one of the plurality of sensor locations to the source of the projectile; andcalculating the final source location using the set of Euclidean distance values for the plurality of sensor location that minimize the error function. 45. The system of claim 44, wherein the error function comprises: ∑wi(ri-dS,Mi)2,i=1…nwhere n is the number of the plurality of locations, ri is the range from an ith one of the plurality of locations to the candidate location, and dS,Mi is the Euclidean distance for between the ith one of the plurality of locations and the candidate location, and wiε{0, 1) are binary weights assigned such that only the N smallest error values are considered. 46. The system of claim 45, wherein said calculating the range if the ith one of the plurality of locations further comprises evaluating: r=12(c4-v4)(P-2Q)where P and Q are defined as P=-2v3dQ,Mv2+c2-2(tmb-tsw)c3v2+2c2dQ,Mvv2+c2-2(tmb-tsw)cv4Q=-2c4v4dQ,M2+2(tmb-tsw)2c6v4+(tmb-tsw)2c4v6-2c7dQ,M(tmb-tsw)vv2+c2+c8(tmb-tsw)2v2+2c8dQ,M2+2v5dQ,Mv2+c2(tmb-tsw)c3,and where tmb−tsw is a time difference of arrival (TDOA) between the muzzle blast and the shockwave for the ith one of the plurality of locations, dQ,M is the miss distance for the ith one of the plurality of locations, v is the velocity of the projectile based on the shockwave angle, and c is the speed of sound. 47. The system of claim 36, further comprising: computing an initial source location for the source based on multilateration using the muzzle blast TOAs; andperforming the step of estimating based on the initial source location. 48. A method for locating a source of a projectile, comprising: obtaining, for a plurality of sensor locations, a length of a shockwave associated with the projectile (shockwave lengths), times of arrival for the shockwave (shockwave TOAs), and muzzle blast TOAs;estimating a trajectory for the projectile based at least one the plurality of sensor locations, the shockwave lengths, and the shockwave TOAs;determining a final source location for the source based on the trajectory and the plurality of sensor locations,wherein the estimating comprises: subdividing a search space into a plurality of points,for each of the plurality of points, calculating shot times for the plurality of sensor locations based on the muzzle blast TOAs, andselecting as the initial source location a one of the plurality of points associated with comparable the shot times among the plurality of sensor locations.
※ AI-Helper는 부적절한 답변을 할 수 있습니다.