RU2719277C1 - Method for determining time with small dispersion of arrival of impulse and device for determining time between sending and arrival of pulse, carrying out method - Google Patents
Method for determining time with small dispersion of arrival of impulse and device for determining time between sending and arrival of pulse, carrying out method Download PDFInfo
- Publication number
- RU2719277C1 RU2719277C1 RU2019123618A RU2019123618A RU2719277C1 RU 2719277 C1 RU2719277 C1 RU 2719277C1 RU 2019123618 A RU2019123618 A RU 2019123618A RU 2019123618 A RU2019123618 A RU 2019123618A RU 2719277 C1 RU2719277 C1 RU 2719277C1
- Authority
- RU
- Russia
- Prior art keywords
- signal
- sine
- arrival
- frequency
- analyzed
- Prior art date
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01H—MEASUREMENT OF MECHANICAL VIBRATIONS OR ULTRASONIC, SONIC OR INFRASONIC WAVES
- G01H5/00—Measuring propagation velocity of ultrasonic, sonic or infrasonic waves, e.g. of pressure waves
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N29/00—Investigating or analysing materials by the use of ultrasonic, sonic or infrasonic waves; Visualisation of the interior of objects by transmitting ultrasonic or sonic waves through the object
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N29/00—Investigating or analysing materials by the use of ultrasonic, sonic or infrasonic waves; Visualisation of the interior of objects by transmitting ultrasonic or sonic waves through the object
- G01N29/36—Detecting the response signal, e.g. electronic circuits specially adapted therefor
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)
Abstract
Description
ОБЛАСТЬ ТЕХНИКИFIELD OF TECHNOLOGY
Заявленное решение относится к измерительной технике, в частности к устройству и способу определения времени прихода импульса путем измерения параметров звуковых и прочих волн, и может быть использовано в области эхо-локации, УЗ позиционирования, в задачах изменения скоростей потока жидкостей или газов (УЗ расходомеров), УЗ дефектоскопии. По аналогии с УЗ, данный способ может применяться для определения времени прихода импульса ЭМ или прочей волны с наименьшей дисперсией, например, в задачах точного измерения времени между приемами (и/или посылками) порций сигнала (пучков импульсов или волн) или в задачах беспроводной синхронизации, тактирования, стробирования на дальних расстояниях и ЭМ позиционирования.The claimed solution relates to measuring equipment, in particular to a device and method for determining the time of arrival of an impulse by measuring the parameters of sound and other waves, and can be used in the field of echo location, ultrasonic positioning, in problems of changing flow rates of liquids or gases (ultrasonic flow meters) , Ultrasonic inspection. By analogy with ultrasound, this method can be used to determine the time of arrival of an EM pulse or other wave with the smallest dispersion, for example, in tasks of accurately measuring the time between receptions (and / or packages) of signal portions (pulse beams or waves) or in problems of wireless synchronization , timing, gating over long distances and EM positioning.
УРОВЕНЬ ТЕХНИКИBACKGROUND
Из уровня техники известно множество решений, направленных на определения времени прихода импульса. Наиболее близким аналогом к заявленном решению являются устройство компенсации погрешности измерения ультразвукового уровнемера, описанное в патенте RU 2470267 C1, опубл. 20.12.2012, выполненное с возможностью измерения временного интервала между излученным и принятым сигналами. Недостатками известного устройства является сложность конструкции, а также низкая точность измерения времени прихода сигнала и низкая помехозащищенность.The prior art many solutions are known aimed at determining the time of arrival of the pulse. The closest analogue to the claimed solution is a device for compensating the measurement error of an ultrasonic level gauge described in patent RU 2470267 C1, publ. 12/20/2012, made with the possibility of measuring the time interval between the emitted and received signals. The disadvantages of the known device is the design complexity, as well as low accuracy of measuring the time of arrival of the signal and low noise immunity.
раскрытие изобретенияdisclosure of invention
Технической задачей, на решение которой направленно заявленное решение, является создание устройства и способа, обеспечивающих более высокую точность определения времени прихода импульса и с наименьшей дисперсией, а также обеспечивающих более высокую повторяемость и помехоустойчивость к отраженным волнам за счет быстрого получения опорной точки на первом фронте огибающей импульсов не зависимо от фазы принятого сигнала до заметного влияния отраженных волн.The technical problem, the solution of which the stated solution is directed, is to create a device and method that provides higher accuracy in determining the time of arrival of the pulse and with the smallest dispersion, as well as providing higher repeatability and noise immunity to reflected waves due to the fast reception of the reference point on the first front of the envelope pulses regardless of the phase of the received signal to a noticeable effect of reflected waves.
Техническим результатом является повышение точности, повторяемости и скорости определения времени принимаемого сигнала не зависимо от его фазы, а также снижение требований к вычислительным ресурсам.The technical result is to increase the accuracy, repeatability and speed of determining the time of the received signal regardless of its phase, as well as reducing the requirements for computing resources.
Для достижения указанного технического результат разработан способ определения времени прихода сигнала, содержащий этапы, на которых: генерируют с помощью излучателя сигнал с заданным конечным количеством периодов волн (несколько импульсов); принимают с помощью приемника сгенерированный сигнал, оцифровывают и определяют огибающую амплитуды (амплитудную модулирующую функцию) сигнала в каждой точке семплирования (дискретизации); определяют анализируемый сигнал на основе огибающей амплитуды принятого сигнала, причем анализируемый сигнал пропорционален мощности принимаемого сигнала; на основе анализруемого сигнала определяют время прихода сигнала.To achieve the indicated technical result, a method has been developed for determining the time of arrival of the signal, comprising the steps of: generating a signal using a radiator with a given finite number of wave periods (several pulses); receive the generated signal using the receiver, digitize and determine the envelope of the amplitude (amplitude modulating function) of the signal at each sampling point (sampling); determine the analyzed signal based on the envelope of the amplitude of the received signal, and the analyzed signal is proportional to the power of the received signal; Based on the analyzed signal, the signal arrival time is determined.
Также разработано устройство определения времени прихода импульсов, содержащее: устройство обработки данных, включающее полосовой фильтр, синусный, косинусный фильтры, умножители (для возведения в квадрат сигналов с выходов синусного и косинусного фильров) и сумматор; приемник; АЦП; генератор тактового сигнала. Выход приемника соединен через АЦП с полосовым фильтром, к выходу которого параллельно подключены синусный и косинусный фильтры, выходы синусного и косинусного фильтров соединены с сумматором через умножители (для возведения в квадрат), при этом устройсво обработки данных выполнено с возможностью: определять огибающую амплитуды сигнала в каждой точке семплирования; определять анализируемый сигнал на основе огибающей амплитуды принятого сигнала, причем анализируемый сигнал пропорционален мощности принимаемого сигнала; и на основе анализруемого сигнала определять время прихода импульсов.Also developed is a device for determining the time of arrival of pulses, comprising: a data processing device including a band-pass filter, sine, cosine filters, multipliers (for squaring the signals from the outputs of the sine and cosine filters) and an adder; receiver; ADC; clock generator. The output of the receiver is connected through an ADC to a bandpass filter, to the output of which a sine and cosine filters are connected in parallel, the outputs of the sine and cosine filters are connected to the adder via multipliers (for squaring), while the data processing device is configured to: determine the envelope of the signal amplitude in each sampling point; determine the analyzed signal based on the envelope of the amplitude of the received signal, and the analyzed signal is proportional to the power of the received signal; and on the basis of the analyzed signal to determine the time of arrival of the pulses.
краткое описание чертежейbrief description of the drawings
Для лучшего представления и понимания сущности изобретения далее будет сделана ссылка, лишь в качестве примера, на прилагаемые чертежи, на которых изображены:For a better presentation and understanding of the invention, reference will hereinafter be made, by way of example only, to the accompanying drawings, in which:
фиг. 1а, 1б - схема устройства определения времени принимаемого сигнала или времени между посылкой и приемом сигналов;FIG. 1a, 1b is a diagram of a device for determining the time of a received signal or the time between sending and receiving signals;
фиг. 2 - форма сигнала напряжения, подаваемая на излучатель (преобразователь);FIG. 2 - voltage waveform supplied to the emitter (converter);
фиг. 3 - форма принимаемого приемником сигнала;FIG. 3 - the form of the signal received by the receiver;
фиг. 4 - форма анализируемого сигнала.FIG. 4 - a form of the analyzed signal.
Осуществление изобретенияThe implementation of the invention
Пример реализации заявленного решения далее подробно будет описан со ссылкой на фиг. 1а, на котором представлена упрощенная техническая схема устройства определения времени принимаемого сигнала, реализующего способ определения времени прихода сигнала, а также способ определения времеми между посылкой и приемом сигнала.An example implementation of the claimed solution will now be described in detail with reference to FIG. 1a, which shows a simplified technical diagram of a device for determining the time of a received signal that implements a method for determining the time of arrival of a signal, as well as a method for determining the time between sending and receiving a signal.
В соответствии с приведенным на фиг. 1а примером устройство определения времени прихода импульсов содержит приемную часть: устройство обработки данных 1, включающее полосовой фильтр 7, синусный 8 и косинусный 9 фильтры, умножители и сумматор 10; приемник 4; АЦП 5; генератор тактового сигнала 6; выход приемника 4 соединен через АЦП 5 с полосовым фильтром 7, к выходу которого параллельно подключены синусный 8 и косинусный 9 фильтры, выходы синусного 8 и косинусного 9 фильтров подключенны к сумматору 10 через умножители для возведения в квадрат, при этом устройсво обработки данных выполнено с возможностью: определять огибающую амплитуды сигнала; определять анализируемый сигнал на основе квадрата огибающей амплитуды принятого сигнала (представлен на выходе сумматора 10), причем анализируемый сигнал пропорционален мощности принимаемого сигнала; и на основе анализруемого сигнала определять время прихода сигнала и между приходами порций сигнала. In accordance with FIG. 1a, an example device for determining the time of arrival of pulses contains a receiving part: a
Для определения времени между посылкой и приходом (приемом) сигнала в частности для УЗ волны и других волн в случаях, когда измеряемый интервал времени больше интервала между отсчетами АЦП, устройство может содержать передающую часть: генератор напряжения 2, излучатель сигнала 3, причем выход устройства обработки данных 1 приемной части в таком случае соединен со входом генератора напряжения 2 передающей части, выход которого соединен с излучателем 3 сигнала, благодаря этой связи приемной и передающей части устройсво обработки данных дополнительно наделено возможностью определять или задавать начало передачи и определять время между посылкой и приходом сигнала с помощью таймера 11, которым наделено устройство обработки данных 1 или установленного между устройством обработки данных 1 и генератором напряжения 2. Генератор тактового сигнала 6, тактирующий АЦП 5 приемной части, может быть соединен с генератором напряжения 2 передающей части для более точного синхронного приема. Устройство обработки данных 1 реализовано на базе по меньшей мере одного процессора или микроконтроллера и в соответствии с заложенным в него программным алгоритмом или по команде, поступившей на вход устройства обработки данных 1, направляет сигналы управления на генератор напряжения 2, который в соответствии с полученными сигналами подает напряжение на излучатель сигнала 3 для генерирования сигнала с заданным количеством периодов волн. Пример формы сигнала напряжения, подаваемого на излучатель, приведен на фиг. 2. В данном случае предлагаемая характерная модель формы сигнала для простоты понимания данного способа изображена такой, чтобы ее огибающая была близка к форме синуса. To determine the time between sending and receiving (receiving) a signal, in particular for an ultrasonic wave and other waves in cases where the measured time interval is greater than the interval between the ADC readings, the device may include a transmitting part:
В качестве излучателя сигнала 3 могут быть использованы различные виды излучателей, например, ультразвуковые (УЗ) излучатели, электромагнитные (ЭМ) излучатели и т.д., соответственно в зависимости от используемого излучателя генерируемый им сигнал будет представлять пучок, состоящий из по меньшей мере одной УЗ, ЭМ или прочей волны.As the emitter of
Сгенерированный излучателем 3 сигнал поступает на приемник 4 и оцифровывается посредством АЦП 5. Оцифровываться принимаемый приемником 4 сигнал может синхронно с излучаемым сигналом с частотой, обеспечивающей восстановление сигнала с заданной (определенной) частотой несущей, например, равной четырехкратной частоте несущей. При таком условии коэффициенты синусного и косинусного фильтра принимают целые значения из ряда [-1, 0, 1] и значительно упрощаются операции вычисления, однако условия оцифровки могут быть выбраны другими. Пример формы принимаемого приемником сигнала для простоты понимания способа представлен на фиг. 3. Синхронность в работе генератора напряжения 2 и АЦП 5 обеспечивается посредством генератора тактового сигнала 6. Далее оцифрованный сигнал поступает на полосовой фильтр 7, который обеспечивает фильтрацию сигнала для увеличения соотношения сигнал/шум, после чего посредством устройства обработки данных 1 определяется анализируемый сигнал как пропорциональный мощности принимаемого сигнала в реальном времени и на основе анализруемого сигнала определяется время прихода сигнала.The signal generated by the
Таким образом, за счет того, что время прихода сигнала определяется на основе анализируемого сигнала, пропорционального мощности принимаемого сигнала, в заявленном решении не используется пиковый детектор, в отличие от решения, описанном в патенте RU 2470267 C1. Следовательно, исключаются помехи и наводки на принимаемый сигнал, возникающие вследствие работы пикового детектора, а также повышается быстродействие устройства. На основе данных анализируемого сигнала (пропорционального мощности) можно любым способом выбрать и определить характерную точку. Наиболее точный результат получается, если в качестве характерной точки выбрана точка перегиба на первом фронте анализируемого сигнала. Далее будет описан пример вычисления анализируемого сигнала посредством устройства обработки данных 1 с соответствующим программным обеспечением.Thus, due to the fact that the signal arrival time is determined based on the analyzed signal, which is proportional to the received signal power, the peak solution is not used in the claimed solution, unlike the solution described in patent RU 2470267 C1. Therefore, interference and interference to the received signal arising from the operation of the peak detector are eliminated, and the device’s speed is also increased. Based on the data of the analyzed signal (proportional power), you can select and determine a characteristic point in any way. The most accurate result is obtained if the inflection point on the first front of the analyzed signal is selected as a characteristic point. Next will be described an example of the calculation of the analyzed signal by the
Для вычисления анализируемого сигнала, пропорционального мощности принимаемого сигнала, выделяется пара квадратурных составляющих сигнала со сдвигом фаз на 90 град. Для этого последовательность отсчетов АЦП принятого сигнала (оцифрованный принятый сигнал) после фильтрации полосовым фильтром пропускают через синусный 8 и косинусный 9 КИХ фильтры (с конечной импульсной характеристкой), например, с коэффициентами Fsin = [0 1 0 -1], Fcos = [1 0 -1 0] для условия частоты дискретизации (семплирования) оцифрованного принятого сигнала, равной четырехкратной частоте несущей. Данные фильтры выполняют функцию квадратурного преообразователя и выделяют кваратурную (синусную) Xsin и синфазную (косинусную) Xcos составляющие сигнала. Эквивалентная схема квадратурного преобразователя посредством фильтров 8 и 9 показана на фиг. 1б. В общем случае коэффициенты синусного 8 и косинусного 9 фильтров представлют собой последовательноти, описывающие один или несколько полных периодов синуса и косинуса соответственно, и определяют (задают) частоту программного косинусно-синусного опорного генератора 12 (фиг. 1б), причем опорная частота этого генератора равна частоте несущей принимаемого сигнала. Процесс фильтрации синусным 8 и косинусным 9 КИХ фильтрами - это умножение принимаемого сигнала заданной частоты на сигналы с косинусно-синусного опорного генератора такой же частоты.To calculate the analyzed signal proportional to the power of the received signal, a pair of quadrature signal components with a phase shift of 90 degrees is extracted. For this, the sequence of ADC samples of the received signal (digitized received signal) after filtering with a band-pass filter is passed through a
Полученные сигналы Xsin(n) и Xcos(n) образуют идеально сопряженые по Гильберту последовательности (со сдвигом фаз ровно 90 градусов), что позволяет точно вычислить огибающую амплитуды принятого сигнала в любой точке независимо от фазы принятого сигнала. Фильтрованные сигналы Xsin(n) и Xcos(n) поступают на устройство обработки данных 1, которое с помощью умножителей и сумматора определяет исследуемый (анализируемый) сигнал по формуле: The received signals Xsin (n) and Xcos (n) form perfectly Hilbert-matched sequences (with a phase shift of exactly 90 degrees), which allows you to accurately calculate the envelope of the amplitude of the received signal at any point, regardless of the phase of the received signal. The filtered signals Xsin (n) and Xcos (n) are fed to a
P(n) = V(n)^2 = Xsin(n)^2 + Xcos(n)^2,P (n) = V (n) ^ 2 = Xsin (n) ^ 2 + Xcos (n) ^ 2,
где V(n) - огибающая сигнала, P(n) = V(n)^2 - исследуемый (анализируемый) сигнал, Xsin(n) - фильтрованный сигнал фильтром Fsin, Xcos(n) - фильтрованный сигнал фильтром Fcos.where V (n) is the envelope of the signal, P (n) = V (n) ^ 2 is the investigated (analyzed) signal, Xsin (n) is the filtered signal with the Fsin filter, Xcos (n) is the filtered signal with the Fcos filter.
Опорную координату устройство обработки данных 1 выделяет на участке с максимальной крутизной фронта исследуемого сигнала в точке перегиба (экстемума производной). Точка перегиба приходит раньше, чем локальный экстремум, поэтому посторонние отраженные волны меньше влияют на форму кривой. Также в точке перегиба производная максимальна, поэтому точность ее вычисления существенно выше, чем экстремума, где производная равна нулю. Физически мощность волны пропорциональна квадрату ее амплитуды, поэтому анализируемый сигнал в физическом смысле пропорционален мощности принятого сигнала. Такой способ дает минимальную погрешность выделения опорной координаты времени. Благодаря крутому фронту исследуемого сигнала (квадрата огибающей амплитуды сигнала), любая погрешность в оси координат приводит к значительно меньшей погрешности во временной оси.The reference coordinate, the
На фиг. 4 представлен пример формы анализируемого сигнала для определения точки максимальной крутизны фронта исследуемого сигнала, где n - моменты дискретизации, tоп - вычисленная опорная координата времени, P' max - значение исследуемого сигнала (мощности) в точке максимальной производной, вычисленной при интерполяции, t - ось времени, P - ось мощности сигнала.In FIG. Figure 4 shows an example of the shape of the analyzed signal to determine the point of maximum steepness of the front of the studied signal, where n are the sampling times, t op is the calculated reference time coordinate, P 'max is the value of the studied signal (power) at the point of maximum derivative calculated by interpolation, t - time axis, P - signal power axis.
Далее для определения точки перегиба (максимальной производной) на каждом дискретном участке исследуемого (анализируемого) сигнала строится кубическая интерполяция по четырем точкам n, n+1, n+2, n+3. Выбор кубической интерполяции обоснован простотой расчетов точки фронта огибающей пришедшей волны. Точки перегиба tоп (максимальной производной) фронта исследуемого сигнала считаем опорной координатой, принимаемого сигнала.Next, to determine the inflection point (maximum derivative) on each discrete section of the investigated (analyzed) signal, cubic interpolation is constructed at four points n, n + 1, n + 2, n + 3. The choice of cubic interpolation is justified by the simplicity of calculating the point of the front of the envelope of the incoming wave. The inflection points t op (maximum derivative) of the front of the investigated signal is considered the reference coordinate of the received signal.
Данный способ определения опорной координаты сигнала обладает относительно высокой повторяемостью, а также низким требованиям к вычислительным ресурсам и позволяет получить точность определения времени прихода импульса значительно больше интервала между моментами дискретизации.This method of determining the reference coordinate of the signal has a relatively high repeatability, as well as low requirements for computing resources and allows you to obtain the accuracy of determining the time of arrival of the pulse is significantly greater than the interval between sampling times.
Claims (36)
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
PCT/RU2018/000308 WO2019221629A1 (en) | 2018-05-16 | 2018-05-16 | Method and device for determining the arrival time of a pulse |
Publications (1)
Publication Number | Publication Date |
---|---|
RU2719277C1 true RU2719277C1 (en) | 2020-04-17 |
Family
ID=68540799
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
RU2019123618A RU2719277C1 (en) | 2018-05-16 | 2018-05-16 | Method for determining time with small dispersion of arrival of impulse and device for determining time between sending and arrival of pulse, carrying out method |
Country Status (2)
Country | Link |
---|---|
RU (1) | RU2719277C1 (en) |
WO (1) | WO2019221629A1 (en) |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
SU1384961A1 (en) * | 1986-10-20 | 1988-03-30 | Уральский политехнический институт им.С.М.Кирова | Device for measuring velocity of ultrasound waves |
SU1587347A1 (en) * | 1988-10-31 | 1990-08-23 | Уральский политехнический институт им.С.М.Кирова | Apparatus for measuring speed and spectral coefficient of damping of ultrasonic waves |
US5793704A (en) * | 1996-12-13 | 1998-08-11 | Solid Scientific Research And Development Ltd. | Method and device for ultrasonic ranging |
RU2470267C1 (en) * | 2011-07-29 | 2012-12-20 | Государственное образовательное учреждение высшего профессионального образования "Национальный исследовательский Томский политехнический университет" | Device to compensate for error in measurement by ultrasonic level metre |
US20140204700A1 (en) * | 2013-01-22 | 2014-07-24 | Schlumberger Technology Corporation | Automatic Processing of Ultrasonic Data |
-
2018
- 2018-05-16 RU RU2019123618A patent/RU2719277C1/en active
- 2018-05-16 WO PCT/RU2018/000308 patent/WO2019221629A1/en active Application Filing
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
SU1384961A1 (en) * | 1986-10-20 | 1988-03-30 | Уральский политехнический институт им.С.М.Кирова | Device for measuring velocity of ultrasound waves |
SU1587347A1 (en) * | 1988-10-31 | 1990-08-23 | Уральский политехнический институт им.С.М.Кирова | Apparatus for measuring speed and spectral coefficient of damping of ultrasonic waves |
US5793704A (en) * | 1996-12-13 | 1998-08-11 | Solid Scientific Research And Development Ltd. | Method and device for ultrasonic ranging |
RU2470267C1 (en) * | 2011-07-29 | 2012-12-20 | Государственное образовательное учреждение высшего профессионального образования "Национальный исследовательский Томский политехнический университет" | Device to compensate for error in measurement by ultrasonic level metre |
US20140204700A1 (en) * | 2013-01-22 | 2014-07-24 | Schlumberger Technology Corporation | Automatic Processing of Ultrasonic Data |
Also Published As
Publication number | Publication date |
---|---|
WO2019221629A1 (en) | 2019-11-21 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
RU2751088C2 (en) | Method and system for processing seismic data | |
Gueuning et al. | Accurate distance measurement by an autonomous ultrasonic system combining time-of-flight and phase-shift methods | |
Andria et al. | Digital signal processing techniques for accurate ultrasonic sensor measurement | |
WO1992007514A1 (en) | Using doppler modulation parameters for amplitude estimation | |
JP4930130B2 (en) | Active sonar device, received signal processing method for sonar, and signal processing program thereof | |
CN103308151A (en) | Heterodyne laser vibration measuring device and method | |
CN104062645A (en) | Method for measuring phase difference between parametric array difference-frequency wave and same-frequency small-amplitude wave signal | |
US4084148A (en) | Object recognition system | |
Maskell et al. | The estimation of subsample time delay of arrival in the discrete-time measurement of phase delay | |
RU2719277C1 (en) | Method for determining time with small dispersion of arrival of impulse and device for determining time between sending and arrival of pulse, carrying out method | |
Andria et al. | Interpolated smoothed pseudo Wigner-Ville distribution for accurate spectrum analysis | |
RU2293997C1 (en) | Method for correlation processing of signals, reflected from fast-moving targets | |
KR100739506B1 (en) | Ultrasonic distance measuring method using matched filter of reduced calculation | |
Andria et al. | Digital measuring techniques for high accuracy ultrasonic sensor application | |
CN109188378B (en) | Ultrasonic Doppler echo signal simulation method | |
RU2677835C1 (en) | Hydro-acoustic tone signals resolution method | |
D'Elia et al. | Software customization to provide digital oscilloscope with enhanced period-measurement features | |
Saulig et al. | Nonstationary signals information content estimation based on the local Rényi entropy in the time-frequency domain | |
CN112763023A (en) | Radar level gauge high-precision measurement output processing method based on optimized data model | |
Katakura et al. | Ultrasonic vector velocity measurement by projection computed velocimetry | |
Shumway et al. | Mixed signal processing for regional and teleseismic arrays | |
Hanh et al. | Source localization on solids utilizing logistic modeling of energy transition in vibration signals | |
Hayman et al. | Signal modelling techniques to extend the measurement frequency range downwards when echo-free time is limited | |
KR101550706B1 (en) | Thickness meter using ultrasonic wave and thickness measuring method of the same | |
Diaz et al. | Time delay measurement based on analytic signal for lightning interferometer |