Модель входного устройства радиометрического приемника для измерения профиля глубинных температур биологических сред методом неинвазивной радиотермометрии

УДК 621.37 А.В. Филатов

Модель входного устройства радиометрического приемника для измерения профиля глубинных температур биологических сред методом неинвазивной радиотермометрии

Рассмотрена модель входного блока микроволнового радиометрического приемника для дистанционных измерений температуры в глубине биологической среды неинвазивным способом, в которой применен ответвитель на связанных линиях передачи. Ответвление части высокочастотного сигнала генератора шума происходит из основного канала в оба рабочих плеча вторичного канала с различными значениями коэффициентов переходного ослабления, разность которых определяет размах диапазона измерения. Это позволяет на стадии моделирования параметров ответвителя создавать радиотермометры с заданными диапазонами измерений. Использование в основе функционирования модели модифицированного метода нулевого приема позволяет снизить влияние на измерения коэффициента отражения в месте приложения аппликаторной антенны к телу и повысить динамические характеристики при картировании биологического объекта.

Ключевые слова: неинвазивный метод радиотермометрии, микроволновая радиотермометрия, измерение глубинных температур биологических сред. ао1: 10.21293/1818-0442-2018-21-4-21-27

В медицинской диагностике при дистанционных измерениях глубинной температуры биологических объектов неинвазивным способом применение микроволновых радиометров имеет неоспоримые преимущества перед инвазивными методами, заключающимися в непосредственном введении в тело объекта миниатюрных термодатчиков, например термопар. К достоинствам радиоволнового метода можно отнести бескровность, раннюю диагностику заболеваний; контроль за ходом лечения (метод безвреден для пациентов, измерения могут повторяться неоднократно), уникальную способность обнаруживать быстрорастущие опухоли (удельное тепловыделение в опухоли прямо пропорционально скорости ее роста).

В настоящее время применяемые в медицинских измерениях радиотермометры [1-7], в основе функционирования которых используется нулевой метод, имеют основной недостаток - низкое быстродействие (измерение в одной точке составляет около 20 с), что значительно снижает возможности радиометрического метода, например по сканированию объекта в ходе его радиотеплового картирования. В таких радиотермометрах состояние нулевого приема достигается в ходе нагрева-охлаждения опорной согласованной нагрузки, шумовая температура которой выравнивается с шумовой температурой антенны. Тепловая инерция снижает динамические свойства приемной радиометрической системы. Во входных блоках известных радиотермометров установлен циркулятор, работающий в режиме вентиля. Ферритовые вентили относятся к дорогостоящим изделиям сверхвысокочастотной (СВЧ) техники. Кроме того, имеют достаточно узкую полосу частот, что для радиометрических измерений является крайне нежелательным, а возникающие в них потери при прохождении сигнала дополнительно ухудшают точность.

В данной работе рассмотрена модель входного блока радиотермометра для дистанционного определения температуры в глубине биологического объекта радиометрическим методом без использования в нем вентильного устройства, содержащая канал формирования подшумливающего объект сигнала для реализации нулевого метода измерений. Под-шумливание объекта осуществляется импульсами, длительность которых зависит от величины сигнала от объекта. Для этого во входном устройстве осуществляется дополнительно к амплитудно-импульсной модуляции широтно-импульсная. Достоинством такой схемы перед известными является более высокое быстродействие. Регулировка нулевого баланса осуществляется по модифицированному методу нулевого приема [8, 9] и заключается не в изменении опорного шумового сигнала до совпадения с входным сигналом, а в изменении времени поступления этого сигнала в измерительный тракт. Электромагнитный сигнал от объекта, непосредственно несущий информацию о глубинной температуре, определяется косвенным образом через длительность ши-ротно-импульсного сигнала, и на этот сигнал не влияет коэффициент отражения в месте приложения аппликаторной антенны к объекту.

В составе модели входного блока используется ответвитель, представляющий собой две линии передачи с распределенной электромагнитной связью и с различными значениями переходного ослабления при ответвлении сигнала из основного (первичного) канала во вспомогательный (вторичный) канал в прямом и противоположном направлениях, т.е., ко -гда возбуждаются оба выходных плеча вторичного канала ответвителя.

На рис. 1 приведена структурная схема модели входного блока, включающего антенну аппликатор-ного типа А, сверхвысокочастотный ключ Кл, ответ-витель, в основной канал которого через переключатель ПК поступает сигнал генератора шума ГШ для

подшумливания объекта исследования. В канале формирования подшумливающего сигнала также установлены две согласованные нагрузки СН и

Граница «объект - антенна»

СН2. Сигнал с выхода ответвителя поступает на радиометрический приемник, собственная шумовая температура которого равна Тш.

Рис. 1. Структурная схема входного блока

РТгш

в1Тгш

Антенна

в1Тгш _

Генератор шума, Тгш

РТ^ —

Приемник

Т

Антенна

Т(1 - Я)(1 - р - РО -

1

Генератор шума

б

Приемник

Рис. 2. Схемы направлений распространения сигналов во вспомогательном канале ответвителя: а - при поступлении сигнала генератора шума на различные входы основного канала; б - при ответвлении части полезного сигнала антенны в основной канал

Возникающее в толще объекта электромагнитное излучение характеризуется эффективной шумовой температурой Т, которая непосредственно определяет его интегральную по глубине термодинамическую температуру. В месте приложения аппликаторной антенны часть этого излучения, равная ТЯ, отражается обратно в объект, где Я - коэффициент отражения по мощности на границе антенны с объектом. Другая часть, Т(1 - Я), принимается антенной и является для радиометрического приемника входным сигналом.

Во входном блоке осуществляются два вида импульсных модуляций: амплитудная и широтная. Амплитудно-импульсная модуляция сигнала антенны выполняется в СВЧ-ключе по меандровому сигналу. Для этого поступающие на вход управления ключом импульсы длительностью tАИМ следуют со скважностью 2. Широтно-импульсной моду-

ляции подвергается сигнал генератора в канале подшумливания. Модуляция осуществляется по управляющему сигналу длительностью tШИМ и заключается в переключении входных плеч основного канала ответвителя. На рис. 2, а показана схема направлений распространения сигналов во вспомогательном канале ответвителя при поступлении высокочастотного сигнала генератора шума с эффективной температурой шумов Тгш на различные входы основного канала. При питании входного плеча 1 основного канала, поступающая во вспомогательный канал часть сигнала генератора шума, равная рТгш, распространяется в сторону антенны, другая часть, р1Тгш, - в сторону входа приемника, где в и р1 - значения переходного ослабления связанных линий, для которых во всей полосе принимаемых частот выполняется условие в1 < в < 1. Аналогичным образом, если сигнал генератора шума поступает во входное плечо 2 основного канала ответвителя, то его часть, равная вТгш, распространяется в сторону радиометрического приемника, а другая часть, в1Тгш, - в сторону антенны.

Принцип работы заключается в следующем. Во входном блоке по управляющему сигналу tАИМ замыкается высокочастотный ключ и сигнал антенны поступает на вход приемника (рис. 3). При замкнутом ключе синхронно с амплитудно-импульсной модуляцией осуществляется широтно-импульсная по сигналу tшиМ. В первом полупериоде амплитудно-импульсной модуляции, в начале импульса длительностью tАИМ, до появления сигнала управления широтной модуляцией, в переключателе выход генератора шума скоммутирован на первый вход основного канала ответвителя. В результате большая часть этого сигнала (вТгш) через СВЧ-ключ поступает в антенну, и отражаясь на границе с объектом, возвращается на вход приемника. Меньшая часть сигнала генератора шума (в1 Тгш), поступая во вторичный канал ответвителя, непосредственно распространяется в сторону при-

1

2

а

в

в

2

емника. При появлении импульса /ШИМ происходит переключение входов основного канала ответвите-ля, они меняются местами. От генератора запиты-вается вход 2, соответственно, большая часть сигнала генератора шума поступает на вход приемника, меньшая - передается в сторону антенны.

I I

Рис. 3. Временные диаграммы управляющих сигналов амплитудно-импульсной и широтно-импульсной модуляциями

Во втором полупериоде амплитудно-импульсной модуляции, когда импульсный сигнал /АИМ отсутствует, ключ во входном блоке размыкается. В этом полупериоде импульсный сигнал управления широтной модуляцией не вырабатывается и сигнал генератора шума в течение всего полупериода поступает во входное плечо 1 основного канала ответвителя.

Таким образом, в зависимости от комбинации управляющих модуляциями сигналов на входе линейного приемника имеют место три сигнала, пропорциональных соответствующим шумовым температурам:

- сигнал А:

Т

- сигнал В:

Т

- сигнал С:

ТгшР + ТгшР: + Тш + Тш = Тгш

Входящие в формулы (1)-(3) эффективные температуры принадлежат различным шумовым источникам, являются некоррелированными и их мощности складываются (шумовая температура является энергетическим параметром поступающего в антенну электромагнитного излучения). В формулах не представлен интерференционный член от взаимодействия собственных шумов приемника и этих же шумов, прошедших через измерительный тракт входного блока и отраженных от антенны. Условием его отсутствия является выполнение требования для минимальной длины волноведущих линий, которая находится из следующего равенства: sмИН = с/(2Д/^е), где с - скорость света в вакууме, е - диэлектрическая проницаемость подложки, на которой размещаются пассивные элементы входного узла, Д/ - полоса принимаемых частот.

В основе функционирования радиотермометра используется модифицированный нулевой метод [9]. Как следует из [10, 11], для данной модификации математическая модель, через которую косвенно определяется измеряемый сигнал антенны через длительность импульса, управляющего широтно-импульсной модуляцией, имеет вид

С - В

А - В

-/а

где А, В, С - значения сигналов из (1)-(3). После подстановки этих значений в (4) и простых преобразований получим

/ш

Тш + Тгшв - Т

тгш (в - Р:;

/ А

Из (5) выразим определяемый сигнал объекта исследований

Т = Т ш + Т гшР - Тгш (Р - М

шим

АИМ

Для крайних значений длительности /ШИМ определяем границы диапазона измерений

Тмакс = Т ш + Т гшР

Тмин = Т ш + Т гшР1

Размах диапазона измерений соответственно составляет величину

<Т = Тмакс - ТмИН = Тгш (р - р1). (9) Из (7)-(9) следует, что необходимый диапазон измерений можно устанавливать не соответствующим изменением эффективных температур опорных источников шума, в частности, согласованной нагрузки вентиля, как это выполнено в [6], а изменением параметров ответвителя, изменяя значения переходного ослабления при передаче сигнала генератора шума из основного канала в выходные плечи вторичного канала. То есть, задача изменения границ диапазона измерений является параметрической и заключается в изменении геометрических размеров ответвителя, обеспечивающих в заданной полосе частот необходимые значения переходного ослабления. На рис. 4 приведена диаграмма диапазона измерения, из которой наглядно следует возможность настройки диапазона на выбранные границы изменением значений коэффициентов переходного ослабления Р и Р! в процессе проектирования ответвителя. При увеличении Р1 диапазон измерений сужается и при Р = Р1, когда равные мощности сигнала генератора шума ответвляются в различных направлениях, сольется в точку.

Р1Тгш

РТП

0 Тш

Т

Т

м

Рис. 4. Диаграмма диапазона измерения

Изменения значений переходного ослабления ß и ßj ответвителя можно выполнить варьированием длины области связи линий передачи и изменением зазора между ними. На рис. 5 приведен график изменения переходного ослабления ß в связанных полосковых линиях передачи, расположенных на подложке с относительной диэлектрической проницаемостью 10, при изменении зазора между ними от 1 до 3 мм. График получен для частоты 3 ГГц в результате моделирования ответвителя в программной среде Microwave Office. На этом графике значения ß были получены для таких длин связанных линий, при которых поступающие сигналы в оба плеча вспомогательного канала одинаковы, т.е. когда выполняется условие равенства значений переходного ослабления ß = ßj.

-15 TR

-40

h, мм

0,5 J 1,5 2 2,5 3 3,5 Рис. 5. Изменение переходного ослабления в связанных линиях передачи от зазора между ними

При исследовании биологических сред наибольший интерес представляют температуры, находящиеся, как правило, в границах ~ 20-50 °С, что в термодинамических температурах составляет ~ 295-325 К [12-15]. С позиции радиометрических измерительных систем данный диапазон можно классифицировать как узкий. С учетом этого рассмотрим последовательность определения геометрических характеристик ответвителя для заданных минимальной Тмин и максимальной Тмакс границ измерительного диапазона.

В качестве исходной информации для моделирования вначале выбирается значение переходного ослабления ß* ответвителя в диапазоне от -20 до -35 дБ. Предпочтение следует отдавать более слабой связи между линиями, чтобы минимизировать ответвление сигнала из вспомогательного канала в основной и тем самым минимизировать потери полезного антенного сигнала при его прохождении через ответвитель (влияние величины переходного ослабления на точность измерений будет рассмотрено ниже). Для выбранного ß* определяется зазор между линиями. Например, для частоты 3 ГГц можно воспользоваться графиком на рис. 5 (аналогичные графики можно построить и для других частот).

После выбора зазора между связанными линиями, дальнейшее моделирование параметров

ответвителя производится изменением размеров области связи линий передачи с использованием метода подбора. Для узкого диапазона измерений параметры ß и ßJ имеют близкие значения, поэтому определение размеров начинается от длины, которая соответствует выбранному значению ß*, когда выполняется равенство ß = ßJ (в этой точке осуществляется смена режима работы ответвите-ля, когда он переходит из противонаправленного режима в сонаправленный). Уменьшение длины связанных линий выполняется с шагом 0,0j мм. На каждом шаге с использованием программной среды моделирования определяется пара значений переходного ослабления ß и ßJ. На этом же шаге для известной эффективной температуры Тш собственных шумов приемника в ходе расчета находится значение ßJ, расч по формуле, полученной из совместного решения уравнений (7) и (8):

ß

J,расч

Т - Т

мин ш

Т - Т

макс ш

ß.

Подстановкой в (J0) полученного из программы значения ß определяется ßJ, расч и сравнивается с полученным в программе значением ßJ. Если ßJ > ßJ, расч, переходят к следующему шагу изменения длины связанных линий, и выше приведенные операции повторяются. Если ßJ ~ ßJ, расч, на этом процесс моделирования ответвителя завершается.

По данному алгоритму был смоделирован от-ветвитель на частоту 3 ГГц для согласованного основного канала. Ответвитель проектировался для диапазона измеряемых температур 295-325 К и собственных шумов с эффективной температурой 70 К. Для полученных значений переходного ослабления ß = -26,05 дБ, ßJ = -26,68 дБ, длина связанных линий составила 13,42 мм при зазоре между ними 2,04 мм.

Ответвитель относится к взаимному реактивному устройству, поэтому в нем происходит не только передача электромагнитной энергии из основного канала во вспомогательный, но и отбор некоторой части распространяющегося во вспомогательном канале высокочастотного сигнала в основной. На рис. 2, б приведена схема направлений ответвления сигнала антенны из вспомогательного канала в основной. На входе приемника измеряемый сигнал уменьшается на величину, пропорциональную J - ß - ßJ. То же самое имеет место и для других сигналов, распространяющихся в обоих направлениях. В этом случае уровни сформированных во входном блоке сигналов будут отличаться от (1)-(3) и соответственно равны:

- сигналА:

Т (J - R)a - ß - ßj) + Trmß + ТГШМ(1 - ß - ßj) + +Тш R(j - ß - ßj)2 + Тш ,

- сигнал В:

T (1 - R)(l - ß - ßj) + TmßR(l - ß - ßj) + Т mßj +

+Тш R(j - ß - ßj)2 + Тш ,

- сигнал С:

ТгшР(1 -р-р1) + ТгшР1 +ТУВЧ (1 -Р-Р1)2 +Тш . Проведя аналогичные выше приведенным преобразования, получим выражения для определения длительности широтно-импульсного сигнала и входного сигнала:

t - (1 - R)V Т+ Т™ß ~ Т t

'ШИМ -

Т = ТшУ + ТгшР - Тгш

(1 - Я)Т 'АИМ

где у = 1 - Р - Р1.

Максимальная Тмакс и минимальная Тмин границы диапазона измерений соответственно определяются из (12) для /шим = 0 и /шим = /аим:

T - Т Y + Т ß

макс ш / гш F > Тмин - Тш Y + Т™ß ~ Тгш

1 - Ry

(1 - '

В соотношениях (13), (14) конечное значение ответвления части сигнала из вторичного канала в первичный непосредственно связано с величиной 1 - Р - Рь Если данным ответвлением сигналов пренебречь, т.е. 1 - Р - Р1 ~ 1, то в этом случае получим (7) и (8).

Рассмотрим влияние вторичного канала от-ветвителя на первичный и возникающие при этом погрешности. Для этого вычитанием (12) из (6) определим величину абсолютной погрешности:

st - тш(ß + ß1; + тгш(ß -ß1)-

Р + Р1 ^шим (15) (1 - Я)У ^АИМ '

Из полученной формулы (15) следует, что абсолютная погрешность зависит от длительности /шим и при /шим = /АИМ принимает максимальное значение, соответствующее нижней границе диапазона измерений:

sТмакс - (ß + ß1)(Tш +

Тгш

(1 - r)y

).

По формуле (16) можно оценить погрешность для известных значений шумовых температур радиометрического приемника, генератора шума, учитывая коэффициенты переходного ослабления сигналов в ответвителе и отражения на границе антенны с объектом. С другой стороны, из этой формулы можно найти такие значения переходного ослабления между связанными линиями, для которых погрешность не выходила бы за заданный предел, определяемый ЛТмакс. Для этого решим равенство (16) относительно коэффициентов переходного ослабления, принимая во внимание, что различия между их значениями имеют небольшую величину и можно считать, что Р = Р1 = Р*. Учитывая выражение (9), получим, что значение переходного ослабления должно удовлетворять следующему неравенству:

ß* <

ST м

Т - Т

+ т + макс мин •

2(sT м

В (17) входит коэффициент отражения Я на границе антенны с объектом. При измерениях значение коэффициента зависит от многих факторов, таких как диэлектрическая проницаемость биологического объекта в месте приложения апплика-торной антенны, качество контакта и т.д., т.е. относится к трудно учитываемым параметрам. Его величина может варьироваться в широких пределах и при каждом измерении может принимать различные значения. В частности, при исследовании биологических сред Я, как правило, изменяется в диапазоне от 0,1 до 0,5. Из (17) следует, что чем больше коэффициент отражения, тем меньшее значение Р* необходимо выбирать, чтобы обеспечить заданную точность. После подстановки в (17) Я = 0,5 получим окончательное выражение для определения переходного ослабления Р*, которое задается в качестве исходной информации в начале моделирования ответвителя:

с 'т7

Р* <-—макс-.

2

Если Р* будет удовлетворять неравенству (18), тогда влиянием погрешности, возникающей в результате передачи сигналов из вспомогательного канала в основной, можно пренебречь.

При моделировании данного блока номинальная выходная мощность генератора шума, выраженная в значениях эффективной шумовой температуры, определяется по формуле, полученной из совместного решения выражений (7), (8):

Т - Т

Т — макс мин

ß - ß1

1 - R

После подстановки в (19) полученных выше значений переходного ослабления ответвителя выходной сигнал генератора шума составил величину 105 К.

Таким образом, если в радиотермометре с установленным на входе вентильным устройством границы диапазона измерений определяет термодинамическая температура согласованной нагрузки и изменения границ можно осуществить только ее нагревом или охлаждением, что требует времени и определенных энергетических затрат, в радиотермометре с данным входным блоком этот недостаток отсутствует, так как изменение нижней и верхней границ можно выполнить соответствующим изменением характеристик ответвителя.

Работа выполнена при поддержке РФФИ

Литература

1. Barrett A. Detection of breast cancer by microwave radiometre / A. Barrett, P.C. Myers, N.L. Sadowsky // Ra-dioSci. - 1977. - Vol. 12, № 68. - P. 167-171.

2. О возможности диагностики рака молочной железы путем измерения собственного электромагнитного излучения тканей (радиотермометрии) / Л.М. Бурдина, А.В. Вайсблат, С.Г. Веснин, Н.Н. Тихомирова // Маммология. - 1997. - № 2. - С. 17-22.

3. Klemetsen Q. Design of medical radiometer frontend for improved performance / Q. Klemetsen, Y. Birkelund, S.K. Jacobsen, P.F. Maccarini, P.R. Stauffer // Progress In Electromagnetics Research B. - 2011. - Vol. 27. - P. 289-306.

4. Asimakis N.P. Non-invasive microwave radiometric system for intracranial applications. A study using the con-formal L-notch microstrip patch antenna / N.P Asimakis, I.S. Karanasiou, N.K. Uzunoglu // Progress In Electromagnetics Research. - 2011. - Vol. 117. - P. 83-101.

5. Годик Э.Э. Человек «глазами радиофизики» / Э.Э. Годик, Ю.В. Гуляев // Радиотехника. - 1991. - № 8. -С. 51-62.

6. Вайсблат А.В. Медицинский радиотермометр РТМ-01-РЭС // Биомедицинские технологии и радиоэлектроника. - 2001. - № 8. - С. 3-9.

7. Маречек С.В. Влияние структуры биоткани на результаты СВЧ-термометрических измерений / С.В. Маречек, В.М. Поляков // Успехи современной радиоэлектроники. - 2001. - № 11. - С. 21-30.

8. Filatov A.V. Application Concept of Zero Method Measurement in Microwave Radiometers // Modern Instrumentation. - 2015. - Vol. 4, No. 3. - P. 19-31. - doi: 10.4236/mi.2015.43003

9. Филатов А.В. Нулевой метод в радиометрических измерениях. - Томск: ТУСУР, 2010. - 206 с.

10. Филатов А.В. Новый подход к построению радиотехнических СВЧ-устройств пассивной локации на принципе синхронного совмещения двух видов импульсной модуляции // Доклады ТУСУР. - 2011. -№ 2(24), ч. 3. - С. 20-26.

11. Filatov A.V. A Microwave Radio Thermometer for Measuring the Depth Temperatures of Biological Objects by a Noninvasive Method / A.V. Filatov, A.G. Lekhanov // Measurement Techniques. - 2015. - Vol. 58, No. 4. -P. 438-445. - doi: 10.1007/s11018-015-0731-1

12. Филатов А.В. Приемный блок нулевого модифицированного микроволнового радиометра для исследования объектов в ближней зоне антенны // Приборы и техника эксперимента. - 2015. - № 1. - С. 82-86.

13. Filatov A.V. A New Principle for Construction of Microwave Multireceiver Radiometers Using a Modified Method of Zero Measurement // Radiophysics and Quantum Electronics. - 2016. - Vol. 59, No. 5. - P. 382-392. -DOI: 10.1007/s11141-016-9707-5

14. Filatov A.V. A Microwave Radiometer for Deep Noninvasive Diagnostics of Thermal Fields Inside a Biological Object / A.V. Filatov, B.V. Utkin, S.E. Tarasov, S.R. Gazitov, M.N. Anishin // Instruments and Experimental Techniques. - 2017. - Vol. 60, No. 3. - P. 367-371. -doi: 10.1134/S002044121703006X

15. Филатов А. В. Новый принцип построения высокочувствительных микроволновых радиометров по многоприемниковой схеме с применением нулевого метода измерений // Физические основы приборостроения. - 2017. - Т. 6, № 2 (24). - С. 67-75.

Филатов Александр Владимирович

Д-р техн. наук, профессор каф. телекоммуникаций и основ радиотехники (ТОР) ТУСУРа Ленина пр-т, д. 40, г. Томск, Россия, 634050 Тел.: +7-952-899-96-41 Эл. почта: filsash@mail.ru

Filatov A.V.

Model of the input device of the radiometric receiver for measuring the profile of the deep temperatures of biological media by non-invasive radio thermometry

A model of the input unit of a microwave radiometric receiver is considered for remote temperature measurements in the depth of the biological environment in a non-invasive way, in which a coupler is used on connected transmission lines. The branch of a part of the high-frequency signal of the noise generator comes from the main channel to both working arms of the secondary channel with different values of the transient attenuation coefficients, the difference of which determines the span of the measurement range. This allows to create radio thermometers with specified measurement ranges at the stage of modeling the parameters of the coupler. Using the model of the modified zero-reception method at the basis of the functioning allows reducing the influence on the measurements of the reflection coefficient at the place of application of the applicator antenna to the body and increasing the dynamic characteristics when mapping a biological object. Keywords: non-invasive method of radio thermometry, microwave radiothermometry, measurement of deep temperatures of biological media. Doi: 10.21293/1818-0442-2018-21-4-21-27

References

1. Barrett A., Myers P.C., Sadowsky N.L. Detection of breast cancer by microwave radiometre. RadioSci, 1977, vol. 12, no. 68, pp. 167-171.

2. Byrdina L.M., Vaisblat A.V., Vesnin S.G., Tixomi-rova N.N. On the possibility of diagnosing breast cancer by measuring the tissue's own electromagnetic radiation (radio thermometry). Mammalogy, 1997, no. 2, pp. 17-22 (in Russ.).

3. Klemetsen Q., Birkelund Y., Jacobsen S.K., Maccarini P.F., Stauffer P.R. Design of medical radiometer frontend for improved performance. Progress In Electromagnetics Research B, 2011, vol. 27, pp. 289-306.

4. Asimakis N.P., Karanasiou I.S., Uzunoglu N.K. Non-invasive microwave radiometric system for intracranial applications. A study using the conformal L-notch microstrip patch antenna. Progress In Electromagnetics Research, 2011, vol. 117, pp. 83-101.

5. Godik J.J., Gyliyev Y.V. Man "through the eyes of radiophysics". Radioengineering, 1991, no. 8, pp. 51-62 (in Russ.).

6. Vaisblat A.V. Medical radiometer RTM-01-RES. Biomedical technology and electronics, 2001, no. 8, pp. 3-9 (in Russ.).

7. Marechek S.V. Poliykov V.M. The influence of the structure of biological tissue on the results of microwave thermometric measurements. Successes of modern radio electronics, 2001, no. 11, pp. 21-30 (in Russ.).

8. Filatov A.V. Application Concept of Zero Method Measurement in Microwave Radiometers. Modern Instrumentation, 2015, vol. 4, no. 3, pp. 19-31.

9. Filatov A.V. Zero method in radiometric measurements [Zero method in radiometric measurements]. Tomsk, TUSUR Publ., 2010. 206 p.

10. Filatov A.V. A new approach to the construction of radio microwave devices, passive location on the principle of synchronous combination of two types of pulse modulation. TUSUR reports, 2011, no. 2(24), ch. 3, pp. 2026 (in Russ.).

11. Filatov A.V., Lekhanov A.G. A Microwave Radio Thermometer for Measuring the Depth Temperatures of Biological Objects by a Noninvasive Method. Measurement Techniques, 2015, vol. 58, no. 4, pp. 438-445.

12. Filatov A.V. Receiving unit of zero modified microwave radiometer for examining objects in the near zone of the antenna. Instruments and Experimental Technique, 2015, no. 1, pp. 82-86 (in Russ.).

13. Filatov A.V. A New Principle for Construction of Microwave Multireceiver Radiometers Using a Modified Method of Zero Measurement. Radiophysics and Quantum Electronics, 2016, vol. 59, no. 5, pp. 382-392.

14. Filatov A.V., Utkin B.V., Tarasov S.E., Gazitov S.R., Anishin M.N. A Microwave Radiometer for Deep Noninvasive Diagnostics of Thermal Fields Inside a Biological Object // Instruments and Experimental Techniques, 2017, vol. 60, no. 3, pp. 367-371.

15. Filatov A.V. A new principle of building highly sensitive microwave radiometers using a multi-receiver circuit using the zero measurement method. Physical Foundations of Instrument Engineering, 2017, no. 2 (24), pp. 67-75.

Alexander V. Filatov

Doctor of Technical Sciences, Professor, Department of Telecommunications and Basic Principles of Radio Engineering, Tomsk State University of Control Systems and Radioelectronics 40, Lenin pr., Tomsk, Russia, 634050 Phone: +7 (382-2) 51-05-30 Email: filsash@mail.ru