RU2108061C1 - Remote diagnosis device for diagnosing respiration and heart beat processes - Google Patents
Remote diagnosis device for diagnosing respiration and heart beat processes Download PDFInfo
- Publication number
- RU2108061C1 RU2108061C1 RU92015438A RU92015438A RU2108061C1 RU 2108061 C1 RU2108061 C1 RU 2108061C1 RU 92015438 A RU92015438 A RU 92015438A RU 92015438 A RU92015438 A RU 92015438A RU 2108061 C1 RU2108061 C1 RU 2108061C1
- Authority
- RU
- Russia
- Prior art keywords
- signal
- antenna
- generator
- output
- conversion unit
- Prior art date
Links
Images
Landscapes
- Measurement Of The Respiration, Hearing Ability, Form, And Blood Characteristics Of Living Organisms (AREA)
Abstract
Description
Изобретение относится к областям медицины и ветеринарии, а именно к технике дистанционной диагностики параметров дыхания и сердцебиения, а также к прикладным областям, связанным с контролем состояния человека или животного в ответственных и экстремальных условиях. The invention relates to the fields of medicine and veterinary medicine, in particular to the technique of remote diagnostics of respiration and heartbeat parameters, as well as to applied areas related to monitoring the condition of a person or animal in critical and extreme conditions.
Известно устройство для контроля за движением (заявка Великобритании N 21664773, 1986, кл. A 61 B 8/00). A device for controlling movement is known (UK application N 21664773, 1986, class A 61
Известное устройство является устройством для дистанционной диагностики, в основном дыхательной функции объекта наблюдения. Устройство построено на принципе доплеровской эхолокации подвижного участка поверхности тела объекта с выделением синхронным детектором сигнала разностной частоты отраженного и излученного сигналов. Целью изобретения является повышение точности и достоверности определения параметров дыхания и сердцебиения за счет обеспечения линейности выходного эффекта измерителя величине смещения поверхности тела под влиянием диагностируемого процесса. The known device is a device for remote diagnostics, mainly respiratory function of the object of observation. The device is based on the principle of Doppler echolocation of a moving portion of an object’s body surface with the synchronous detector detecting the difference frequency signal of the reflected and emitted signals. The aim of the invention is to increase the accuracy and reliability of determining the parameters of respiration and heartbeat by ensuring the linearity of the output effect of the meter to the magnitude of the displacement of the body surface under the influence of the diagnosed process.
Предложенные устройства, как и прибор-прототип, построены на принципе анализа информации, заложенной в принимаемом ультразвуковом сигнале и обусловленной модуляцией этого сигнала на подвижной от дыхания и сердцебиения поверхности тела объекта. Принципиальным является то, что предложенные устройства осуществляют слежение за текущей фазой принимаемого сигнала в широком диапазоне ее отклонений и тем самым обеспечивают измерение с последующим спектральным анализом непосредственно исследуемой функции смещения поверхности тела объекта Δt вместо "свернутой" функции вида cos[γΔ(t)] без потери чувствительности устройства за счет увеличения (для избежания нелинейности) рабочей длины волны λo . Предложено 3 варианта конфигурации устройства, соответствующие используемому методу относительной фазовой дальнометрии.The proposed devices, as well as the prototype device, are based on the principle of analyzing the information embedded in the received ultrasonic signal and due to the modulation of this signal on the surface of the body that is movable from breathing and heartbeat. It is important that the proposed devices monitor the current phase of the received signal in a wide range of its deviations and thereby provide a measurement with subsequent spectral analysis of the directly investigated function of the displacement of the body surface of the object Δt instead of a “convoluted” function of the form cos [γΔ (t)] without loss of sensitivity of the device due to the increase (to avoid non-linearity) of the working wavelength λ o . 3 variants of the device configuration are proposed, corresponding to the used method of relative phase ranging.
Устройство по п. 1 формулы, как и устройство-прототип, работает по сигналу, отраженному от поверхности тела объекта. Состав устройства включает излучающую антенну и принимающую антенну, неподвижные одна относительно другой (установленные в едином корпусе или разнесенные в пространстве), узел когерентного преобразования, узел спектральной обработки и опорный генератор, обеспечивающий формирование излучаемого и опорного сигналов. Узел когерентного преобразования выполнен в виде линейного измерителя фазы (ЛИФ) принимаемого сигнала. The device according to
Устройство по п. 2 формулы реализует однопросветный вариант измерения смещения поверхности тела. В дополнение к узлам устройства по п. 1 в состав устройства введен автономный генератор, формирующий излучаемый сигнал и связанный с излучающей антенной. Этот генератор и излучающая антенна выполнены в виде автономного конструктивного элемента - датчика. Датчик выполнен в миниатюрном исполнении и размещен на объекте таким образом, что его излучающая антенна, подвержена смещению под влиянием исследуемого(ых) процесса(ов), например, в области солнечного сплетения при одновременной диагностике параметров дыхания и сердцебиения. Прием сигнала, модулированного смещением излучающей антенны, осуществляется неподвижной принимающей антенной с распространением сигнала только в одном направлении. Требуемая когерентность излучаемого и опорного сигналов обеспечивается внутренней стабильностью соответствующих независимых генераторов. The device according to
Третий вариант устройства (устройство по п. 3 формулы) в двухпросветной схеме измерения, подобной схеме первого варианта, предусматривает использование для дальнометрии вместо сигнала, отраженного от тела объекта, ответного сигнала, формируемого на подвижном участке поверхности тела объекта специальным датчиком - активным когерентным ультразвуковым приемоответчиком. Когерентность датчика предполагает сохранение в ответном сигнале структуры фазовой модуляции, заложенной во входном сигнале, в общем случае с умножением полной фазы входного сигнала на коэффициент трансформации, равный или отличный от 1. The third version of the device (the device according to
В состав устройства входят все узлы устройства по п. 1. Вводимый в схему датчик содержит последовательно соединенные антенну приема, когерентный преобразователь (преобразователь частоты или просто усилитель при обеспечении пространственной развязки между антеннами) и антенну излучения. Как и в устройстве по п. 2, датчик размещается на объекте таким образом, что его антенны подвержены смещению под влиянием диагностируемого процесса. The device includes all nodes of the device according to
Во всех предложенных устройствах излучающей антенны излучается непрерывный немодулированный сигнал. Длина волны сигнала соизмерима с амплитудой смещения поверхности тела объекта под влиянием исследуемого процесса. In all proposed devices of the radiating antenna, a continuous unmodulated signal is emitted. The wavelength of the signal is comparable with the amplitude of the displacement of the surface of the body of the object under the influence of the investigated process.
ЛИФ в устройствах по п.п. 1-3 охарактеризован на уровне функционального обобщения как устройство, обеспечивающее измерение изменяющейся фазы принимаемого сигнала в широком (более ±2π ) диапазоне отклонений относительно фазы опорного сигнала с линейной зависимостью выходного эффекта от величины отклонения фазы. LIF in devices according to 1-3 is characterized at the level of functional generalization as a device that provides measurement of the changing phase of the received signal in a wide (more than ± 2π) range of deviations relative to the phase of the reference signal with a linear dependence of the output effect on the magnitude of the phase deviation.
ЛИФ может быть построен на различных схемно-технических принципах, например на основе замкнутой следящей системы, на принципе коррекции значений фазы, измеряемых по модулю 2π , при достижении углов, кратных 2π , с использованием схем с фазовой автоподстройкой частоты, на основе цифрового измерителя фазы интегрирующего типа с квадратурными каналами приема. Квадратурные каналы также могут формироваться различными способами, например, с использованием аналоговых квадратурных смесителей, со взятием мгновенных отсчетов принимаемого узкополосного сигнала непосредственно на рабочей частоте устройства и другими. В 4 формулы представлен частный случай аналого-цифрового исполнения ЛИФ устройств по п.п. 1 или 2, или 3 с аналоговым формированием квадратурных сигналов и последующей цифровой оценкой текущей фазы по отсчетам значений квадратурных сигналов. The LIF can be built on various circuit-technical principles, for example, on the basis of a closed-loop tracking system, on the principle of correcting phase values measured modulo 2π, when angles are multiples of 2π, using phase-locked loop, based on a digital phase meter integrating type with quadrature reception channels. Quadrature channels can also be formed in various ways, for example, using analog quadrature mixers, taking instantaneous samples of the received narrow-band signal directly at the operating frequency of the device and others. In 4 formulas presents a special case of analog-to-digital execution of the LIF devices according to paragraphs. 1 or 2, or 3 with analog generation of quadrature signals and subsequent digital estimation of the current phase from the samples of values of quadrature signals.
ЛИФ по п. 4 формулы содержит два квадратурных канала приема сигнала с выхода антенны. Каждый из каналов включает последовательно соединенные смеситель, фильтр, усилитель и схему выборки и хранения. В состав ЛИФ входят также фазовращатель на π/2 опорного сигнала одного из смесителей, мультиплексор выборок квадратурных сигналов обоих каналов, аналого-цифровой преобразователь этих выборок и вычислитль текущего значения фазы принимаемого сигнала. The LIF according to
Технический результат от исполнения изобретения заключается в следующем:
для устройства по п. 1 - обеспечено адекватное представление исследуемой функции Δt на входе узла спектральной обработки. В результате исключается искажение временной и спектральной структур исследуемого процесса, повышается уровень исследуемых составляющих спектра по отношению к уровню флюктуационных шумов и помех от побочных движений, исключается подавление малого сигнала, обусловленного сердцебиением, при совместном наблюдении процессов дыхания и сердцебиения, исключается зависимость в результатах измерений от глубины дыхания и расстояния от прибора до объекта;
для устройств по п.п. 2 и 3 - то же, что для устройства по п. 1; получен энергетический выигрыш в отношении сигнал / шум за счет использования активного излучающего датчика и (для устройства по п. 2) одностороннего распространения сигнала; практически полностью снята помеха в виде отражения от побочных участков тела объекта и окружающих предметов (использование датчиков эквивалентного случаю отражения от "блестящей точки" на теле объекта) при снятии необходимости в формировании узких диаграмм направленности.The technical result of the invention is as follows:
for the device according to
for devices according to 2 and 3 - the same as for the device according to
Предложенные технические решения имеют одинаковое назначение, направлены на достижение одного и того же технического результата (повышение точности измерений), решаются одним и тем же путем и не могут быть охарактеризованы в одном независимом пункте. Это послужило основанием для их представления как группы изобретений в качестве вариантов выполнения одного устройства. The proposed technical solutions have the same purpose, are aimed at achieving the same technical result (improving the accuracy of measurements), are solved in the same way and cannot be characterized in one independent point. This served as the basis for their presentation as a group of inventions as embodiments of a single device.
На фиг. 1 изображена структурная схема устройства по п. 1; на фиг. 2 - структурная схема устройства по п. 2; на фиг. 3 - структурная схема устройства по п. 3; на фиг. 4 - структурная схема устройства по п. 1 с аналого-цифровым ЛИФ по п. 4 формулы. In FIG. 1 shows a block diagram of a device according to
Устройство по п. 1 (фиг. 1) содержит опорный генератор 1, излучающую антенну 2 (излучатель), принимающую антенну 3 (микрофон), линейный измеритель фазы (ЛИФ) 4, узел спектральной обработки 5 и узел вывода 6. Антенны 2 и 3 неподвижны относительно друг друга. The device according to claim 1 (Fig. 1) comprises a
Устройство работает по сигналу, излучаемому антенной 2 и отражаемому от тела объекта. На излучающую антенну 2 подается непрерывный немодулированный сигнал рабочей частоты от опорного генератора опорного. Этот же сигнал поступает на ЛИФ 4 в качестве опорного. Диаграмма излучения-приема устройства (антенн 2 и 3) направляется на исследуемый объект в область, подверженную смещению под влиянием дыхания и/или сердцебиения, например в область солнечного сплетения. При этом отраженный от объекта сигнал оказывается промодулированным по фазе в соответствии с изменением расстояния от прибора до отражающего участка поверхности тела объекта под воздействием исследуемого процесса. Отраженный сигнал принимается принимающей антенной 3. Выход антенны 3 связан с ЛИФ 4, который обеспечивает измерение с устранением неоднозначности сдвига фазы принимаемого сигнала относительно опорного. Выходной сигнал ЛИФ подается на узел спектральной обработки 5, где производится частотная фильтрация сигнальной (полезной) составляющей на фоне шумов и помех от побочных движений, взаимное частотное разделение сигнальных составляющих, обусловленных дыханием и сердцебиением, их тонкий спектральный анализ, высокочастотная коррекция временных выходных сигналов и пр. Узел вывода 6 обеспечивает оценку числовых параметров исследуемого процесса и отображение выходной информации. Вывод информации может осуществляться на дисплей с отображением временных разверток и спектрограмм процессов дыхания и/или сердцебиения и их параметров, исполнительное устройство, устройство аудио- или визуализации наличия динамической составляющей и выходном сигнале устройства (для устройства, предназначенного для поиска людей в экстремальных ситуациях) и другие приборы в зависимости от конкретного назначения устройства. The device operates on a signal emitted by
Устройство по п. 2 (фиг. 2) содержит датчик 1, включающий в свой состав генератор 2 сигнала рабочей частоты и излучающую антенну 3, принимающую антенну 4. ЛИФ 5, узел спектральной обработки 6, узел вывода 7 и опорный генератор 8. Датчик 1 выполнен в миниатюрном виде и размещен на теле объекта так, что его антенна 3 подвержена смещению под влиянием дыхания и/или сердцебиения. Устройство работает по сигналу, излучаемому датчиком 1. В датчике от генератора 2 на антенну 3 подается непрерывный немодулированный сигнал. Излучаемый антенной 3 ультразвуковой сигнал промодулирован в соответствии с ее смещениями, в том числе под влиянием дыхания и сердцебиения. Этот сигнал принимается принимающей антенной 4. Выход антенны 4 связан с ЛИФ 5, который обеспечивает линейное измерение сдвига фазы принимаемого сигнала относительно опорного. Опорный сигнал поступает на ЛИФ от опорного генератора. Выход ЛИФ связан с узлом спектральной обработки 6 и далее с узлом вывода 7. Функции узлов 6 и 7 аналогичны функциям рассмотренных узлов 5 и 6 фиг. 1. The device according to claim 2 (Fig. 2) contains a
Устройство по п. 3 (фиг. 3) содержит опорный генератор 1, излучающую антенну 2, датчик 3, включающий в свой состав антенну приема 4, когерентный преобразователь 5 и антенну излучения 6, принимающую антенну 7. ЛИФ 8, узел спектральной обработки 9 и узел вывода 10. Подобно датчику устройства по п. 2, датчик 3 выполнен в миниатюрном виде и размещен на теле объекта так, что его антенны 4 и 6 подвержены смещению под влиянием исследуемого процесса. The device according to claim 3 (Fig. 3) contains a
Устройство работает по сигналу от излучающей антенны 2, когерентно переизлученному датчиком 4. На излучающую антенну 2 непрерывный немодулированный сигнал рабочей частоты излучения подается от опорного генератора 1. От этого же генератора на ЛИФ 8 подается опорный сигнал, когерентный сигналу рабочей частоты. Датчик 3 представляет собой активный когерентный приемоответчик. Его антенна приема принимает сигнал от излучающей антенны 2. Сигнал с выхода антенны 4 усиливается и когерентно преобразуется по частоте с сохранением структуры фазовой модуляции (в частном случае при обеспечении пространственной развязки между антеннами - просто усиливается) в преобразователе когерентном 5 и переизлучается антенной излучения 6. Ответный сигнал датчика 3 на входе принимающих антенны 7 промодулирован по фазе в соответствии со смещениями антенн 4 и 6 датчика 3 под воздействием дыхания и сердцебиения объекта. В ЛИФ 8 обеспечивается линейное измерение сдвига фазы принимаемого сигнала относительно опорного. Выходной сигнал ЛИФ подается на узел спектральной обработки 9, где, как и в схеме фиг. 1, производится частотная фильтрация и анализ компонентов исследуемого процесса. В узле вывода 10 проводится оценка числовых параметров исследуемого процесса и отображения выходной информации. The device operates on a signal from a
Аналого-цифровое построение ЛИФ по п. 4 формулы требует использования вычислителя, реализующего расчет текущей фазы сигнала из значений выборок сигналов квадратурных каналов. В схеме устройства по п. 1 с ЛИФ по п. 4 (фиг. 4) указанный вычислитель реализован программным путем в микроЭВМ, позволившей решить ряд других задач устройства в целом. Схема фиг. 4 включает в себя ранее рассмотренные опорный генератор 1, излучающую антенну 2 и принимающую антенну 3. Схема содержит два квадратурных канала приема 4 с и 4s, каждый из которых включает последовательно соединенные смеситель 5, фильтр 6, усилитель 7 и схему выборки и хранения (СВХ) 8. Опорный сигнал на смеситель канала 4s подается через фазовращатель (ФЭВР) 9 на π/2 . Схема содержит также мультиплексор 10, аналого-цифровой преобразователь (АЦП) 11, микроЭВМ 12 и аппаратуру вывода 13. В микроЭВМ 12 реализованы следующие программные блоки: блоки оценки текущей фазы сигнала (блок ОТФС) 14: блок спектральной обработки (блок CO) 15; блок подготовки данных к выводу (блок ПДКВ) 16. The analog-digital construction of the LIF according to
Блок ПДКВ совместно с аппарурой вывода 13 образует узел вывода 17. Block PDKV together with the
Опорный генератор 1 выполнен в виде единого генератора рабочего сигнала и тактовых импульсных сигналов управления СВХ 8, мультиплексора 10, АЦП 11, микроЭВМ 12 и аппаратуры вывода 13. Объединенные информационные входы смесителей квадратурных каналов и входы этих смесителей по опорным сигналам (для канала 4s), с включением Ф3ВР 9) образуют соответственно информационный и опорный входы ЛИФ, а выход программного блока ОТФС - его выход. The
Снимаемый принимающей антенной 3 входной информационный сигнал преобразуется смесителями 5 и фильтрами 6 в два низкочастотных сигнала Uc(t) и Us(t) - аналоговых по выходу усилителей 7 и дискретизированных с частотой дискретизации Fg по выходу СВХ 8. При фазовом сдвиге в опорных сигналах смесителей 5, равном π/2 , квадратурные сигналы Uc(t) и Us(t) представляют собой огибающие косинусной и синусной составляющих принимаемого высокочастотного сигнала и адекватно отражают структуру его модуляции. Мультиплексор 10 за один такт дискретизации поочередно коммутирует выборки квадратурных сигналов на групповое АЦП 11, откуда их квантованные значения поступают на программный блок ОТФС 14.The input information signal recorded by the receiving
В блоке ОТФС использован алгоритм оценки фазы интегрирующего типа:
,
где
φ(tn) - оценка текущей фазы принимаемого сигнала на момент времени tn = n/ng, в начале цикла измерения φ(to) устанавливается равным 0;
n, k - номер такта дискретизации на цикле измерения;
Δφ(tn) - приращение фазы на n-м такте дискретизации.In the OTFS block, an integrating type phase estimation algorithm is used:
,
Where
φ (t n ) - assessment of the current phase of the received signal at time t n = n / n g , at the beginning of the measurement cycle φ (t o ) is set to 0;
n, k is the number of the sampling clock on the measurement cycle;
Δφ (t n ) is the phase increment at the nth sampling cycle.
Величина Δφ(tn) на каждом такте дискретизации рассчитывается по следующей формуле:
где
Uc,s(tn) - значения косинусной и синусной выборок сигналов на входе блока ОТФС на n-м такте дискретизации,
Δ Uc,s(tn)= Uc,s(tn)-Uc,s(tn-i) - приращения соответствующих выборок на n-м такте дискретизации.The value Δφ (t n ) at each sampling cycle is calculated by the following formula:
Where
U c, s (t n ) are the values of the cosine and sine samples of the signals at the input of the OTFS block at the nth sampling cycle,
Δ U c, s (t n ) = U c, s (t n ) -U c, s (t ni ) are the increments of the corresponding samples at the nth sampling cycle.
Выходная информация ЛИФ в виде оценки текущей фазы сигнала поступает в программный блок спектральной обработки 15. Функции блока CO 15 и узла вывода 17 рассмотрены ранее при описании соответствующих узлов схемы фиг. 1. The output of the LIF in the form of an estimate of the current phase of the signal enters the spectral
Включение ЛИФ по п. 4 в схемы устройств по п.п. 2 и 3 полностью идентично приведенному на фиг. 4.2 The inclusion of the LIF according to
Claims (4)
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
RU92015438A RU2108061C1 (en) | 1992-12-30 | 1992-12-30 | Remote diagnosis device for diagnosing respiration and heart beat processes |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
RU92015438A RU2108061C1 (en) | 1992-12-30 | 1992-12-30 | Remote diagnosis device for diagnosing respiration and heart beat processes |
Publications (2)
Publication Number | Publication Date |
---|---|
RU92015438A RU92015438A (en) | 1995-04-20 |
RU2108061C1 true RU2108061C1 (en) | 1998-04-10 |
Family
ID=20134780
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
RU92015438A RU2108061C1 (en) | 1992-12-30 | 1992-12-30 | Remote diagnosis device for diagnosing respiration and heart beat processes |
Country Status (1)
Country | Link |
---|---|
RU (1) | RU2108061C1 (en) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
RU2496410C1 (en) * | 2012-04-20 | 2013-10-27 | Федеральное государственное бюджетное образовательное учреждение высшего профессионального образования "Рязанский государственный радиотехнический университет" | Device for distance registration of processes of patient's heart beats and respiration |
RU2648602C2 (en) * | 2013-02-28 | 2018-03-26 | Конинклейке Филипс Н.В. | Device and method for detection of subjects on basis of vital signs |
-
1992
- 1992-12-30 RU RU92015438A patent/RU2108061C1/en active
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
RU2496410C1 (en) * | 2012-04-20 | 2013-10-27 | Федеральное государственное бюджетное образовательное учреждение высшего профессионального образования "Рязанский государственный радиотехнический университет" | Device for distance registration of processes of patient's heart beats and respiration |
RU2648602C2 (en) * | 2013-02-28 | 2018-03-26 | Конинклейке Филипс Н.В. | Device and method for detection of subjects on basis of vital signs |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Li et al. | A review on recent progress of portable short-range noncontact microwave radar systems | |
US7753849B2 (en) | Doppler radar cardiopulmonary sensor and signal processing system and method for use therewith | |
Nosrati et al. | High-accuracy heart rate variability monitoring using Doppler radar based on Gaussian pulse train modeling and FTPR algorithm | |
Li et al. | Robust overnight monitoring of human vital signs by a non-contact respiration and heartbeat detector | |
JP3877783B2 (en) | A method for finding the position of a living organism and a microwave probe using the | |
Girbau et al. | Remote sensing of vital signs using a Doppler radar and diversity to overcome null detection | |
Bakhtiari et al. | A real-time heart rate analysis for a remote millimeter wave IQ sensor | |
US9993175B2 (en) | Pulsed ultra-wideband sensor and the method thereof | |
Will et al. | Local pulse wave detection using continuous wave radar systems | |
Gao et al. | Radius correction technique for Doppler radar noncontact periodic displacement measurement | |
RU2313108C2 (en) | Mode of detection living objects and an arrangement for its execution | |
US20080275337A1 (en) | Breathing detection apparatus and method | |
Prat et al. | Collimated beam FMCW radar for vital sign patient monitoring | |
Rice et al. | A wireless multifunctional radar-based displacement sensor for structural health monitoring | |
JP2019525195A (en) | Method for processing signals originating from coherent riders and associated rider systems | |
Ambrosanio et al. | A multi-channel ultrasound system for non-contact heart rate monitoring | |
CN113608206B (en) | Method and apparatus for detecting movement of a subject | |
Park et al. | Center tracking quadrature demodulation for a Doppler radar motion detector | |
Michler et al. | Micrometer sensing with microwaves: Precise radar systems for innovative measurement applications | |
Das et al. | Antenna evaluation of a non-contact vital signs sensor for continuous heart and respiration rate monitoring | |
Michler et al. | Pulse wave velocity detection using a 24 GHz six-port based Doppler radar | |
Vinci et al. | 24 GHz six-port medical radar for contactless respiration detection and heartbeat monitoring | |
RU2108061C1 (en) | Remote diagnosis device for diagnosing respiration and heart beat processes | |
Kontou et al. | Heartbeat and respiration detection using a low complexity cw radar system | |
Hati et al. | $ W $-band vibrometer for noncontact thermoacoustic imaging |