RU2541682C1 - Method of analysis of acoustic field of complex source - Google Patents
Method of analysis of acoustic field of complex source Download PDFInfo
- Publication number
- RU2541682C1 RU2541682C1 RU2013141773/28A RU2013141773A RU2541682C1 RU 2541682 C1 RU2541682 C1 RU 2541682C1 RU 2013141773/28 A RU2013141773/28 A RU 2013141773/28A RU 2013141773 A RU2013141773 A RU 2013141773A RU 2541682 C1 RU2541682 C1 RU 2541682C1
- Authority
- RU
- Russia
- Prior art keywords
- signals
- field
- received signals
- acoustic
- antenna
- Prior art date
Links
Images
Landscapes
- Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
- Monitoring And Testing Of Transmission In General (AREA)
Abstract
Description
Изобретение относится к измерительной технике и может быть использовано для поиска зон повышенного акустического излучения по длине транспортных средств - на автомобильном или железнодорожном транспорте, а также на судах различного назначения при их диагностическом обследовании.The invention relates to measuring equipment and can be used to search for areas of increased acoustic radiation along the length of vehicles - in automobile or railway transport, as well as on ships for various purposes during their diagnostic examination.
Известен способ аддитивного направленного приема, основанный на использовании интерференции сигналов, снимаемых с выходов отдельных приемников при падении на антенну волны от акустического источника (См.: Новиков А.К. Статистические измерения в судовой акустике. Л.: Судостроение, 1985, стр.263-266).A known method of additive directional reception, based on the use of interference of signals from the outputs of individual receivers when a wave is incident on an antenna from an acoustic source (See: A. Novikov, Statistical Measurements in Ship Acoustics. L .: Sudostroenie, 1985, p. 263 -266).
Сущность способа поясним на примере испытаний движущихся транспортных средств, когда транспортное средство проходит на близком расстоянии от приемной системы измерительных преобразователей. В основе способа - определение звукового давления, обусловленного источниками акустического поля - различными виброактивными механизмами, распределенными по длине испытуемого транспортного средства. При этом положение источника звука (шум машин и механизмов), а следовательно, и положение источника энергии гидроакустического поля определяется путем сопоставления максимумов суммарного звукового давления с теми частями транспортного средства, которые в момент возникновения этих максимумов оказываются в створе диаграммы направленности системы отдельных приемников.The essence of the method is illustrated by the example of tests of moving vehicles when the vehicle passes at a close distance from the receiving system of the measuring transducers. The basis of the method is the determination of sound pressure due to sources of acoustic field - various vibroactive mechanisms distributed along the length of the test vehicle. In this case, the position of the sound source (noise of machines and mechanisms), and therefore the position of the energy source of the hydroacoustic field, is determined by comparing the maximums of the total sound pressure with those parts of the vehicle that, at the time of the occurrence of these maximums, are in the alignment of the radiation pattern of the individual receiver system.
Недостатком способа является то, что расстояния между приемными элементами направленной системы должны быть меньше половины длины волны принимаемого сигнала. При поиске зон повышенного акустического излучения протяженного объекта измерений разрешающая способность способа определяется апертурой антенны.The disadvantage of this method is that the distance between the receiving elements of the directional system must be less than half the wavelength of the received signal. When searching for areas of increased acoustic radiation of an extended measurement object, the resolution of the method is determined by the aperture of the antenna.
Наиболее близким по технической сущности к предлагаемому решению является способ направленного приема, основанный на использовании интерференции сигналов, снимаемых с выходов отдельных приемников при падении на антенну акустической волны от источника с последующим умножением частот принятых сигналов, суммированием исходных сигналов, суммированием сигналов с умноженными частотами и перемножением сумм сигналов, полученных в результате детектирования (см. патент Российской Федерации №2478982, МПК G01S 15/00, 2011 г., Гарин В.Ю., Стефанский В.М.). Этот способ определения звукового давления принят за прототип. При поиске зон доминирующего акустического излучения объекта измерений необходимо относительное перемещение системы приемников и протяженного источника акустического поля.The closest in technical essence to the proposed solution is the method of directional reception, based on the use of interference of signals taken from the outputs of individual receivers when an acoustic wave is incident from the source upon an antenna, followed by multiplying the frequencies of the received signals, summing the original signals, summing the signals with multiplied frequencies and multiplying the sum of the signals obtained as a result of detection (see the patent of the Russian Federation No. 2478982, IPC G01S 15/00, 2011, Garin V.Yu., Stefansky V.M.). This method of determining sound pressure is taken as a prototype. When searching for zones of dominant acoustic radiation from the measurement object, relative movement of the receiver system and the extended source of the acoustic field is necessary.
Сущность способа-прототипа для аддитивной антенны в режиме приема сводится к следующим операциям:The essence of the prototype method for an additive antenna in receive mode is reduced to the following operations:
1) производится многоканальный прием сигналов излучения;1) multichannel reception of radiation signals;
2) производится фильтрация принятых сигналов на частотах повышенного излучения;2) filtering the received signals at high radiation frequencies;
3) производится умножение частот принятых сигналов в число раз, соответствующее количеству приемников в системе;3) the frequencies of the received signals are multiplied by the number of times corresponding to the number of receivers in the system;
4) производится суммирование исходных сигналов;4) the summation of the source signals;
5) производится суммирование сигналов с умноженными частотами;5) the summation of signals with multiplied frequencies;
6) выделяется огибающая суммы сигналов с исходными частотами путем амплитудного детектирования;6) the envelope of the sum of the signals with the initial frequencies is highlighted by amplitude detection;
7) выделяется огибающая суммы сигналов с умноженными частотами путем амплитудного детектирования;7) the envelope of the sum of signals with multiplied frequencies is extracted by amplitude detection;
8) производится перемножение огибающих сумм сигналов;8) the envelopes of the sums of the signals are multiplied;
9) осуществляется регистрация результата перемножения сумм сигналов.9) the result of multiplying the sums of the signals is recorded.
В последнем случае поиск максимума результата перемножения сумм сигналов с исходными частотами и сумм сигналов с умноженными частотами, соответствующего положению источника звука по длине транспортного средства, производится в результате относительного движения транспортного средства и неподвижной приемной системы.In the latter case, the search for the maximum result of multiplying the sum of the signals with the original frequencies and the sum of the signals with the multiplied frequencies corresponding to the position of the sound source along the length of the vehicle is performed as a result of the relative motion of the vehicle and the stationary receiving system.
Недостаток способа-прототипа состоит в необходимости относительного движения объекта измерений и приемной системы. Для условий измерений в ограниченных акваториях, например в гидроакустическом бассейне, этот режим нереализуем.The disadvantage of the prototype method is the need for relative movement of the measurement object and the receiving system. For measurement conditions in limited waters, for example, in a hydroacoustic basin, this mode is not feasible.
Задачей изобретения является обеспечение возможности поиска доминирующих зон акустического поля без относительного перемещения сложного источника и приемной системы, что упрощает процесс измерения.The objective of the invention is to enable the search for dominant zones of the acoustic field without the relative movement of a complex source and receiving system, which simplifies the measurement process.
Это достигается тем, что после полосовой фильтрации запоминают принятые сигналы, а затем многократно воспроизводят с ними операции способа, причем в каждом цикле воспроизведения принятые сигналы и сигналы с умноженными частотами задерживают в каждом канале на время распространения акустической волны от контролируемой точки поля до соответствующего приемника дискретной антенны. После чего запоминают результат перемножения огибающих сигналов, полученных после детектирования сигналов для контролируемых точек поля. Изменение координат контролируемых точек поля в очередном цикле воспроизведения и регистрация совокупности результатов измерений позволяют получить оценку интенсивности акустического поля, обусловленного излучением сложного источника.This is achieved by the fact that after bandpass filtering, the received signals are stored and then the method operations are repeatedly reproduced with them, and in each playback cycle, the received signals and signals with multiplied frequencies are delayed in each channel for the duration of the propagation of the acoustic wave from the controlled field point to the corresponding discrete receiver antennas. After that, the result of multiplying the envelopes of the signals obtained after detecting the signals for the controlled points of the field is stored. Changing the coordinates of the controlled field points in the next playback cycle and recording the totality of the measurement results allow us to estimate the intensity of the acoustic field due to the radiation of a complex source.
Сущность предложенного технического решения поясняется фигурами, где:The essence of the proposed technical solution is illustrated by figures, where:
на фиг.1 представлена структурная схема устройства, реализующего предлагаемый способ;figure 1 presents the structural diagram of a device that implements the proposed method;
на фиг.2 представлена планарная диаграмма направленности антенны с лепестком, направленным в центральную точку поля, формируемого сложным источником;figure 2 presents a planar radiation pattern of the antenna with a petal directed to the center point of the field formed by a complex source;
на фиг.3 представлена планарная диаграмма направленности антенны с лепестком, направленным в центральную точку поля, формируемого сложным источником, полученная в результате перемножения двух диаграмм: диаграммы направленности, построенной по принятым сигналам, и диаграммы направленности, построенной по сигналам с умноженными мгновенными частотами.figure 3 presents a planar radiation pattern of an antenna with a petal directed to the center point of the field formed by a complex source, obtained by multiplying two diagrams: a radiation pattern constructed from received signals and a radiation pattern constructed from signals with multiplied instantaneous frequencies.
на фиг.4 представлена планарная диаграмма направленности антенны, сфокусированной в заданную точку анализируемого поля, полученная до умножения частоты сигналов в каждом приемнике;figure 4 presents a planar radiation pattern of the antenna focused at a given point in the analyzed field, obtained before multiplying the frequency of the signals in each receiver;
на фиг.5 представлена планарная диаграмма направленности антенны, сфокусированной в заданную точку анализируемого поля после умножения частоты сигналов в каждом приемнике;figure 5 presents a planar radiation pattern of the antenna focused at a given point in the analyzed field after multiplying the frequency of the signals in each receiver;
на фиг.6 представлено трехмерное изображение результирующей планарной диаграммы направленности антенны, полученное по предлагаемому способу;figure 6 presents a three-dimensional image of the resulting planar radiation pattern of the antenna obtained by the proposed method;
на фиг.7 представлена трехмерная характеристика, полученная в результате взаимодействия направленных свойств антенны, сформированная по предлагаемому способу, и акустического поля, обусловленного излучением сложного источника.figure 7 presents a three-dimensional characteristic obtained by the interaction of the directed properties of the antenna, formed by the proposed method, and the acoustic field due to radiation from a complex source.
Устройство, реализующее предлагаемый способ (см. фиг.1), содержит блок приемников 1, соединенный через блок полосовых фильтров 2 с запоминающим устройством 3, выход которого подключен к входу первого блока задержек 4 и блока умножения частоты 5, соединенного с входом второго блока задержек 6, а выходы блоков задержек 4 и 6 подключены к входам сумматоров 7 и 8, выходы которых соединены с входами блока детектирования и перемножения 9, выход которого подключен к входу регистратора 10, причем к управляющим входам блоков задержек 4 и 6 подключен выход блока управления 11.A device that implements the proposed method (see Fig. 1) contains a
С использованием описанного устройства предложенный способ реализуется следующим образом.Using the described device, the proposed method is implemented as follows.
Входной сигнал принимается блоком приемников 1, в качестве которого рассматривается дискретная антенна, состоящая из дискретных измерительных преобразователей - микрофонов или гидрофонов. После полосовой фильтрации сигналов, снимаемых с выходов измерительных преобразователей блока приемников 1, в блоке фильтров 2, сигналы запоминают в запоминающем устройстве 3 и при воспроизведении умножают по частоте в блоке умножения частоты 5. После чего осуществляется их задержка в блоках 4 и 6 на время распространения акустической волны от контролируемой точки поля до соответствующего приемника дискретной антенны и их суммирование в сумматоре 7. Эта операция позволяет сфокусировать диаграмму направленности на контролируемую точку акустического поля, формируемого сложным источником. Расчет требуемых задержек для блоков задержек 4 и 6 осуществляется при помощи устройства управления 11. После умножения частоты в блоке умножения частоты 5 и суммирования сигналов с умноженными частотами в сумматоре 8 формируется диаграмма направленности, характеризуемая появлением боковых лепестков, равных по величине основному. Сигналы с сумматоров 7 и 8 детектируются и перемножаются в блоке детектирования и перемножения 9. В результате этой операции подавляется влияние боковых лепестков, а результирующая диаграмма направленности обеспечивает повышенное пространственное разрешение. В результате формируется распределение амплитуд результирующих сигналов для совокупности контролируемых точек поля сложного источника.The input signal is received by the
Рассмотрим пример компьютерного моделирования работы устройства, реализующего предлагаемый способ.Consider an example of computer simulation of a device that implements the proposed method.
В результате моделирования рассчитывались сигналы, наблюдаемые на выходе антенны, состоящей из девяти эквидистантных приемных элементов, от точечного монохроматического излучателя, перемещающегося относительно приемников антенны на различных удалениях. В результате сканирования формируется традиционная планарная диаграмма направленности, пример которой показан на фиг.2. На фиг.3 представлена планарная диаграмма направленности антенны с лепестком, направленным в центральную точку поля, формируемого сложным источником, полученная в результате перемножения двух диаграмм: диаграммы направленности, построенной по принятым сигналам, и диаграммы направленности, построенной по сигналам с умноженными мгновенными частотами. Как и в способе-прототипе, умножение сигнала по частоте осуществляется путем выделения мгновенной амплитуды и мгновенной фазы сигнала с ее умножением в заданное число раз. При этом осуществляют умножение мгновенной частоты. С этой целью используется преобразование Гильберта узкополосного сигнала и модуляция гармонического сигнала на средней частоте полосового фильтра по амплитуде и по фазе после ее умножения. Применение способа-прототипа позволяет увеличить эффективность направленного приема для различных частот без изменения апертуры антенны, но не обеспечивает возможность анализа его пространственной структуры. Введение требуемых задержек в устройство по способу-прототипу и в устройство по предлагаемому способу (фиг.4) позволяет сфокусировать дискретную антенну в произвольную точку поля. На фиг.5 показан пример планарной диаграммы направленности с повышенным пространственным разрешением, основной лепесток которой направлен в заданную точку поля. На фиг.6 в качестве примера приведено трехмерное изображение планарной диаграммы направленности, полученной по предлагаемому способу.As a result of the simulation, the signals observed at the output of the antenna, consisting of nine equidistant receiving elements, from a point monochromatic emitter moving relative to the antenna receivers at different distances were calculated. As a result of the scan, a traditional planar radiation pattern is formed, an example of which is shown in FIG. Figure 3 presents a planar radiation pattern of an antenna with a lobe directed to the center point of the field formed by a complex source, obtained by multiplying two diagrams: a radiation pattern constructed from the received signals and a radiation pattern constructed from signals with multiplied instantaneous frequencies. As in the prototype method, the frequency is multiplied by highlighting the instantaneous amplitude and instantaneous phase of the signal with its multiplication by a specified number of times. In this case, the instantaneous frequency is multiplied. For this purpose, the Hilbert transform of a narrow-band signal and the modulation of the harmonic signal at the middle frequency of the band-pass filter in amplitude and phase after its multiplication are used. The use of the prototype method allows to increase the efficiency of directional reception for various frequencies without changing the aperture of the antenna, but does not provide the ability to analyze its spatial structure. The introduction of the required delays in the device according to the prototype method and in the device according to the proposed method (figure 4) allows you to focus the discrete antenna at an arbitrary point on the field. Figure 5 shows an example of a planar radiation pattern with increased spatial resolution, the main lobe of which is directed to a given point in the field. Figure 6 as an example shows a three-dimensional image of a planar radiation pattern obtained by the proposed method.
Эффективность предлагаемого способа подтверждается результатами следующего модельного эксперимента.The effectiveness of the proposed method is confirmed by the results of the following model experiment.
Источник тонального сигнала на частоте 1000 Гц перемещался в пределах плоскости, ограниченной размерами Х=±10 м и Y=±10 м. Исследуемый сигнал принимался антенной из 9 приемных элементов, расположенных на окружности. Расстояние между приемными элементами составляло 0,75 м. Дистанция между осью антенны и плоскостью, в которой перемещался источник, равнялась 8 м. Оценка амплитуды сигнала на выходе антенны производилась через каждые 0,2 м. В результате совместного представления измерительных данных строилась планарная диаграмма направленности. На приемные элементы антенны воздействовали сигналы, обусловленные работой 5 излучателей, расположенных в пределах квадрата 2 м × 2 м.The source of the tonal signal at a frequency of 1000 Hz moved within a plane limited by the dimensions X = ± 10 m and Y = ± 10 m. The signal under study was received by an antenna of 9 receiving elements located on a circle. The distance between the receiving elements was 0.75 m. The distance between the axis of the antenna and the plane in which the source moved was 8 m. The signal amplitude at the antenna output was estimated every 0.2 m. As a result of the joint presentation of the measurement data, a planar radiation pattern was constructed . The receiving elements of the antenna were affected by signals due to the work of 5 emitters located within a square of 2 m × 2 m.
На фиг.7 показан результат взаимодействия направленных свойств антенны, сформированный по предлагаемому способу, и акустического поля, обусловленного излучением сложного источника.7 shows the result of the interaction of the directional properties of the antenna, formed by the proposed method, and the acoustic field due to radiation from a complex source.
Проведенное моделирование показало, что использование новых операций в соответствии с предлагаемым способом анализа акустического поля позволило определять положение доминирующих зон излучения в пределах заданной плоскости.The simulation showed that the use of new operations in accordance with the proposed method for the analysis of the acoustic field made it possible to determine the position of the dominant radiation zones within a given plane.
Основное преимущество предлагаемого способа определения звукового давления перед способом-прототипом состоит в том, что достигается изменение чувствительности приемной системы по пространству путем введения компенсационных задержек в тракты каналов отдельных приемников. В этом случае не требуется относительного перемещения сложного источника относительно приемной системы.The main advantage of the proposed method for determining sound pressure over the prototype method is that a change in the sensitivity of the receiving system in space is achieved by introducing compensation delays in the channel paths of individual receivers. In this case, the relative movement of the complex source relative to the receiving system is not required.
Claims (1)
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
RU2013141773/28A RU2541682C1 (en) | 2013-09-11 | 2013-09-11 | Method of analysis of acoustic field of complex source |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
RU2013141773/28A RU2541682C1 (en) | 2013-09-11 | 2013-09-11 | Method of analysis of acoustic field of complex source |
Publications (2)
Publication Number | Publication Date |
---|---|
RU2541682C1 true RU2541682C1 (en) | 2015-02-20 |
RU2013141773A RU2013141773A (en) | 2015-03-20 |
Family
ID=53285475
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
RU2013141773/28A RU2541682C1 (en) | 2013-09-11 | 2013-09-11 | Method of analysis of acoustic field of complex source |
Country Status (1)
Country | Link |
---|---|
RU (1) | RU2541682C1 (en) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
RU2625617C2 (en) * | 2015-10-22 | 2017-07-17 | Российская Федерация, От Имени Которой Выступает Министерство Промышленности И Торговли Российской Федерации | Method of measuring reflection coefficient of sound-absorbing structure |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
RU2116705C1 (en) * | 1997-04-17 | 1998-07-27 | Государственное предприятие "Всероссийский научно-исследовательский институт физико-технических и радиотехнических измерений" | Method for modeling hydroacoustic fields |
WO2001098796A2 (en) * | 2000-06-22 | 2001-12-27 | Lockheed Martin Corporation | Acoustic surveillance system with signature and bearing data display |
RU2246132C2 (en) * | 2003-01-09 | 2005-02-10 | Военно-морской институт радиоэлектроники | Method and device for fast computing of discontinuous wavelet-conversion of signal with random discretization step of scale coefficients |
RU2305297C2 (en) * | 2005-07-05 | 2007-08-27 | Тихоокеанский военно-морской институт им. С.О. Макарова | Device with correlation shaper of directivity characteristics for detection of signals and determination of direction to their source |
RU2346294C2 (en) * | 2007-04-09 | 2009-02-10 | ОАО "Концерн "Океанприбор" | Method for correcting amplitude-phase distribution of multichannel hydroacoustic antenna excitation |
RU2367970C2 (en) * | 2007-11-28 | 2009-09-20 | Государственное образовательное учреждение высшего профессионального образования Военно-морская академия имени Адмирала Флота Советского Союза Н.Г. КУЗНЕЦОВА | Device for detection of narrow-band noise hydroacoustic signals based on calculation of integral wavelet-spectrum |
-
2013
- 2013-09-11 RU RU2013141773/28A patent/RU2541682C1/en active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
RU2116705C1 (en) * | 1997-04-17 | 1998-07-27 | Государственное предприятие "Всероссийский научно-исследовательский институт физико-технических и радиотехнических измерений" | Method for modeling hydroacoustic fields |
WO2001098796A2 (en) * | 2000-06-22 | 2001-12-27 | Lockheed Martin Corporation | Acoustic surveillance system with signature and bearing data display |
RU2246132C2 (en) * | 2003-01-09 | 2005-02-10 | Военно-морской институт радиоэлектроники | Method and device for fast computing of discontinuous wavelet-conversion of signal with random discretization step of scale coefficients |
RU2305297C2 (en) * | 2005-07-05 | 2007-08-27 | Тихоокеанский военно-морской институт им. С.О. Макарова | Device with correlation shaper of directivity characteristics for detection of signals and determination of direction to their source |
RU2346294C2 (en) * | 2007-04-09 | 2009-02-10 | ОАО "Концерн "Океанприбор" | Method for correcting amplitude-phase distribution of multichannel hydroacoustic antenna excitation |
RU2367970C2 (en) * | 2007-11-28 | 2009-09-20 | Государственное образовательное учреждение высшего профессионального образования Военно-морская академия имени Адмирала Флота Советского Союза Н.Г. КУЗНЕЦОВА | Device for detection of narrow-band noise hydroacoustic signals based on calculation of integral wavelet-spectrum |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
RU2625617C2 (en) * | 2015-10-22 | 2017-07-17 | Российская Федерация, От Имени Которой Выступает Министерство Промышленности И Торговли Российской Федерации | Method of measuring reflection coefficient of sound-absorbing structure |
Also Published As
Publication number | Publication date |
---|---|
RU2013141773A (en) | 2015-03-20 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106680825A (en) | Acoustic array imaging system and method thereof | |
RU2014133018A (en) | METHOD FOR SIMULTANEOUS INCLUSION OF MARINE VIBRATORS | |
Anderson et al. | Time reversal techniques | |
CN103149274B (en) | A kind of concrete defect detection method | |
CN109061654B (en) | Single-circular-ring-array active three-dimensional positioning method in deep sea environment | |
CN104655728A (en) | Acoustic phased array imaging method | |
US10111644B2 (en) | Method of coherent flow imaging using synthetic transmit focusing and acoustic reciprocity | |
RU2541682C1 (en) | Method of analysis of acoustic field of complex source | |
CN110269643A (en) | Self-adapting clutter filtering in the ultrasonic imaging based on acoustic radiation power | |
CA2915303C (en) | Saft analysis of defects close to the surface | |
Ma et al. | Spatiotemporal two-dimensional deconvolution beam imaging technology | |
Touzé et al. | Double-Capon and double-MUSICAL for arrival separation and observable estimation in an acoustic waveguide | |
RU69646U1 (en) | PARAMETRIC ECHO-PULSE LOCATOR | |
Matveev et al. | Forward scattering observation with partially coherent spatial processing of vertical array signals in shallow water | |
Bianchi et al. | A plenacoustic approach to acoustic signal extraction | |
Lobréau et al. | Hemispherical double-layer time reversal imaging in reverberant and noisy environments at audible frequencies | |
CN112834623B (en) | Acoustic covering layer echo reduction measurement method based on compressed sensing technology | |
Li et al. | Research on the imaging of concrete defect based on the pulse compression technique | |
Prada | Detection and imaging in complex media with the DORT method | |
RU2478982C2 (en) | Method for determining acoustic pressure of moving extended source of acoustic field | |
RU2525701C1 (en) | Method to suppress reverberation interference | |
Lu et al. | Imaging a crack in a thin plate by reverse time migration with dispersive flexural waves | |
To et al. | Power estimation of sound sources on low-speed electric trains using a deconvolution approach | |
Tao et al. | The lateral resolution of three-dimensional underground imaging by using amplitude correlation synthesis processing method | |
JP2544535B2 (en) | Measuring device and measuring method for moving sound source |