На правах рукописи
КОЛЯДОВ Дмитрий Валерьевич
РАДИОПОЛЯРИМЕТРИЯ ПРИРОДНЫХ СРЕД И ТЕХНОГЕННЫХ ОБЪЕКТОВ ПРИ ИХ МОНИТОРИНГЕ
Специальность - 01.04.03 УРадиофизикаФ
АВТОРЕФЕРАТ
диссертации на соискание ученой степени доктора физико-математических наук
Москва - 2009 г.
Работа выполнена на кафедре УАвиационных радиоэлектронных системФ Федерального государственного образовательного учреждения высшего профессионального образования УМосковский государственный технический университет гражданской авиацииФ МГТУ ГА.
Научный консультант: Заслуженный деятель науки и техники РФ, зав. кафедрой АРЭС, профессор, доктор физико-математических наук Козлов Анатолий Иванович
Официальные оппоненты: Лауреат Государственной премии СССР и РФ, Заслуженный деятель науки и техники РФ, зав.
абораторией ИРЭ РАН, профессор, доктор физико-математических наук Кутуза Борис Георгиевич Лауреат Государственной премии РФ, зам. директора по НР ФГУП УГидрометпоставкаФ профессор, доктор физико-математических наук Горелик Андрей Габриэлович Лауреат Государственной премии СССР, декан факультета информационных систем и компьютерных технологий РосНОУ, профессор, доктор физико-математических наук Крюковский Андрей Сергеевич
Ведущая организация: ОАО УНПО УРадар-ММСФФ
Защита диссертации состоится У 24 У июня 2009 г. на заседании диссертационного совета Д 212.156.06 в Московском физико-техническом институте по адресу: 117393, г. Москва, ул. Профсоюзная, д. 84/32, корп. В-2.
Отзывы направлять по адресу: 141700 г. Долгопрудный, Московская обл., Институтский переулок, д.9, МФТИ, Диссертационный совет Д 212.156.06.
С диссертацией можно ознакомиться в библиотеке МФТИ.
Автореферат разослан У Ф 2009 г.
Ученый секретарь диссертационного совета Д. 212.156.06:
доцент, кандидат технических наук Чубинский Н.П.
Общая характеристика работы
Актуальность работы. В настоящее время в различных областях жизнедеятельности человека интенсивно развиваются технологии дистанционного зондирования, направленные на извлечение и использование некоординатной информации о техногенных объектах и различных природных средах при их мониторинге. Под некоординатной информацией имеются в виду сведения о типе и свойствах исследуемого объекта, конфигурации и его линейных размерах (физические, электрофизические и геометрические параметры).
Развитие указанных технологий не в последнюю очередь объясняется негативным влиянием человеческой деятельности на экологическую обстановку и необходимостью проведения различных геофизических и гидрофизических исследований.
Основными методами решения различных задач мониторинга техногенных объектов и земных покровов являются разнообразные методы дистанционного зондирования, среди которых можно выделить оптические, активные и пассивные радиолокационные методы. Однако при их использовании возникают проблемы обеспечения высокой достоверности извлекаемой некоординатной информации и увеличения числа анализируемых информационных параметров.
Использование оптических методов серьезно ограничивается условиями, в которых производится мониторинг исследуемых объектов.
Существующие радиотехнические методы дистанционного зондирования протяженных объектов и природных сред в большей мере основаны на анализе энергетических характеристик электромагнитной волны, отраженной от анализируемых объектов.
Методы дистанционного зондирования, основанные на анализе только энергетических параметров отраженной электромагнитной волны, в большей степени исчерпали свои возможности. По этой причине использование информации, заложенной в поляризационных свойствах отраженной электромагнитной волны и в статистических параметрах ее поляризационных характеристик, дает дополнительные возможности для улучшения процедуры извлечения как координатной, так и некоординатной информации и повысить эффективность дистанционных исследований, конечной целью которых является определение физических, химических, геометрических и электрофизических характеристик объектов мониторинга.
Как показывают теоретические разработки и подтверждают экспериментальные исследования, при использовании радиополяриметрии открываются новые возможности для определения физических (соленость, влажность, состав грунта и т.п.), электрофизических (диэлектрическая проницаемость и проводимость), геометрических (структура, неровность, шероховатость и т.п.) характеристик объектов мониторинга, получение которых при помощи известных методов дистанционного зондирования крайне затруднено. В этой связи, диссертационная работа, содержащая решение научной проблемы разработки и использования методов радиополяриметрии для мониторинга природных сред и техногенных объектов, является актуальной.
Цели и задачи исследования. Целью работы является разработка методов определения физических, химических, электрических и геометрических характеристик природных сред и техногенных объектов на основе анализа рассеянного и собственного радиоизлучения объектов мониторинга.
Для достижения поставленной цели необходимо было решить следующие задачи:
1.Провести анализ традиционных методов дистанционного определения физических, химических, электрических и геометрических характеристик природных сред и техногенных объектов мониторинга.
2.Оценить возможность и целесообразность применения радиополяриметрии при дистанционном определении электрофизических и геометрических характеристик природных сред и техногенных объектов мониторинга.
3.Оценить возможность применения радиополяриметрии для определения электрофизических характеристик объектов мониторинга, обладающих свойствами нелинейного рассеивателя.
4.Определить взаимосвязи электродинамических и поляризационных характеристик отраженных объектами мониторинга электромагнитных волн.
5.Обосновать возможность применения методов радиополяриметрии для определения электрофизических характеристик объектов мониторинга при их дистанционном зондировании.
6. Обосновать возможность применения методов радиополяриметрии для улучшения обнаружения техногенных объектов на фоне природных сред при их мониторинге.
7.Обосновать возможность применения методов радиополяриметрии для визуализации, классификации и идентификации природных средств и техногенных объектов.
Научная новизна работы заключается в том, что в ней:
1.Установлена прямая функциональная связь между электродинамическими (комплексная диэлектрическая проницаемость) и поляризационными характеристиками (элементы матрицы рассеяния) объектов мониторинга.
2.Предложен метод визуализации объектов мониторинга при помощи KLL-сферы, позволяющий реализовать их идентификацию и классификацию.
3.Разработан метод определения комплексной диэлектрической проницаемости широкого класса природных объектов на основе анализа поляризационного состояния отраженных этими объектами электромагнитных волн.
5.Экспериментально определены статистические характеристики элементов матрицы рассеяния некоторых классов подстилающих покровов.
6.Разработан метод оценки геометрических характеристик техногенных объектов, находящихся на фоне природных сред, путем управления поляризационными характеристиками электромагнитных волн.
7.Разработаны математические модели неоднородных техногенных объектов, обладающих нелинейными свойствами, при рассеянии электромагнитных волн на таких объектах, дающие возможность проводить их анализ.
Практическая значимость работы состоит в том, что полученные в ней результаты позволяют:
1. Определять при дистанционном зондировании окружающей среды физические характеристики подстилающих покровов (соленость, влажность, плотность, температуру и т.п.) 2. Оценивать при дистанционном зондировании окружающей среды параметры геометрических неоднородностей.
3.Увеличивать степень различения техногенных объектов мониторинга, находящихся на поверхности природных сред.
4.Оценивать геометрические размеры объектов мониторинга.
На защиту выносятся:
1. Методы определения физических (соленость, влажность, состав грунта и т.п.), электрофизических (диэлектрическая проницаемость и проводимость), геометрических (структура, неровность, шероховатость и т.п.) характеристик природных покровов при их дистанционном зондировании путем управления поляризационным состоянием электромагнитных волн.
2. Методы определения комплексной диэлектрической проницаемости при наличии полной и неполной информации о поляризационном состоянии рассеянных объектами мониторинга электромагнитных волн.
3. Методы визуализации, классификации и идентификации объектов мониторинга.
4. Методы оценки геометрических параметров поверхности объектов мониторинга при когерентном и некогерентном рассеянии электромагнитных волн.
5. Методы оценки геометрических параметров поверхности объектов мониторинга путем управления поляризационными характеристиками приемных антенн.
6. Методы определения поляризационных характеристик объектов мониторинга, обладающих нелинейными свойствами.
Внедрение результатов. Основные результаты работы нашли применение в разработках предприятий МКБ УКомпасФ, УИнтеррадиоФ и ЦНИИ УРадиосвязьФ, ОАО УНПП УРадар-ММСФФ, о чем имеются соответствующие акты о внедрении.
Апробация работы. Материалы диссертации докладывались на LVI научной сессии, посвященной дню радио (Москва, 2004), Международной научно-технической конференции MIKON 2006 (Microwave Week, Warsawa, Poland, May 2006), Международной научно-технической конференции EuRAD (EuMW) 2006 (European Microwave Week, European Radar, United Kingdom, Manchester, 2006), 2-й Международной конференции УАкустооптические и радиолокационные методы измерений и обработки информацииФ, (Суздаль, 2007), XX научной конференция МФТИ УСовременные проблемы фундаментальных и прикладных наукФ (Москва, 2007), Международной научно-технической конференции УСовременные научнотехнические проблемы гражданской авиацииФ (Москва 2002, 2004, 2006, 2008 гг.), на научном семинаре Института радиоэлектроники РАН (Москва, 2008), а также на научно-технических семинарах кафедры физикоматематических проблем МФТИ (Москва 2008 г.), кафедр УАвиационных радиоэлектронных системФ и УТехнической эксплуатации радиотехнического оборудования и связиФ МГТУ ГА (Москва 2000-2008 гг.).
По материалам диссертации опубликовано 45 работ. Из них статьей в научных журналах и изданиях, рекомендованных ВАК Минобразования России для публикации основных научных результатов диссертаций на соискание ученой степени доктора наук, 30 статей в иных журналах и изданиях, и в виде 4 параграфов объемом 13 стр. в монографии УПоляризация радиоволн. Поляризационная структура радиолокационных сигналовФ (Москва, УРадиотехникаФ, 2005), 3 параграфов объемом 37 стр. в монографии УПоляризация радиоволн. Радиолокационная поляриметрияФ (Москва, УРадиотехникаФ, 2007), 2 параграфов объемом 33 стр. в монографии УПоляризация радиоволн. Радиополяримерия сложных по структуре сигналовФ (Москва, УРадиотехникаФ, 2008).
Работа состоит из введения, пяти глав, заключения, списка цитируемой литературы (150 наименований). Общий объем диссертации 260 листов сквозной нумерации, основной текст изложен на 250 листах. Диссертация содержит 98 рисунков и 4 таблицы.
Содержание работы Во введении проведено обоснование актуальности темы работы, определена главная цель исследования, рассматриваются известные подходы к решению поставленных задач, определяется научная новизна и практическая ценность полученных результатов, выделены основные научные положения, выносимые на защиту и дана общая характеристика работы.
В первой главе диссертационной работы рассматриваются природные и техногенные объекты как объекты радиополяриметрического мониторинга, определяются функциональные связи между электродинамическими и электрофизическими характеристиками объектов мониторинга, с одной стороны, и между электродинамическими и поляризационными характеристиками с другой стороны, а также зависимость между геометрическими и поляризационными характеристиками объектов мониторинга.
Задачи рассеяния электромагнитной волны являются классическими в электродинамике. Однако использование соответствующих решений уравнений Максвелла в большинстве случаев не дает возможности непосредственного определения информации об объектах мониторинга. Это, в первую очередь, связано со значительным числом факторов, влияющих на электродинамические характеристики объектов мониторинга. Во вторую очередь, с необходимостью учета случайного характера рельефа поверхностей и изменения электродинамических характеристик как вдоль, так и внутрь объектов мониторинга.
На практике является достаточным наличие оценочной информации об электродинамических характеристиках объектов, полученной на основе обработки статистической информации об объектах радиолокационного мониторинга.
Как известно, в качестве основной электродинамической характеристики объектов мониторинга выступает комплексная диэлектрическая про ницаемость = - i, которая определяется большим числом физических характеристик: температура, влажность, соленость, плотность, пористость и т.д. При этом выделить зависимость комплексной диэлектрической проницаемости от одного какого-либо параметра достаточно проблематично.
Большинство природных объектов мониторинга можно рассматривать как многокомпонентные среды (смеси). Например, снег является смесью льда и воздуха, горные породы включают кроме породообразных минералов еще газ или нефть и воду. Диэлектрическая проницаемость смеси зависит от ее строения.
В работе рассматриваются следующие виды природных объектов мониторинга: пресная и соленая вода, пресный и морской лед, сухой и влажный снег, минералы, горные породы, сухая и влажная почва, растительные покровы и лесные породы. В работе проводится анализ известных аналитических зависимостей между комплексной диэлектрической проницаемостью и физическими характеристиками таких объектов мониторинга.
При проведении измерений, направленных на определение физической структуры и геометрических характеристик объектов мониторинга, основное значение имеет измерение элементов матрицы рассеяния волны таких объектов. По этой причине в работе решается задача установления прямой связи между комплексной диэлектрической проницаемостью и поляризационными характеристиками объектов мониторинга, под которыми понимаются элементы матрицы рассеяния, степень поляризационной анизотропии и др.
Для решения задач мониторинга важно знать поляризационные характеристики отраженных от объектов электромагнитных волн. В первую очередь, имеется в виду знание элементов матрицы рассеяния облучаемых площадок на поверхности. При изменении физических свойств таких элементов поверхности будет происходить изменение основной электродинамической характеристики поверхности - ее комплексной диэлектрической проницаемости, изменение которой приводит к изменению отражательных характеристик (элементов матрицы рассеяния).
В работе рассматривается наиболее используемый случай дистанционного зондирования объекта мониторинга. Проводится наблюдение относительно гладкого участка поверхности под углом . Размеры участка определяются шириной диаграммы направленности антенны.
В работе показывается, что в рассматриваемом случае, матрица рассеяния объекта мониторинга может быть представлена в виде произведения двух множителей: первый из них зависит только от комплексной диэлектрической проницаемости, а второй определяется только геометрией задачи и характеристиками антенны:
RHP 01 1S = = (1) ( ) ( )1+ 0 RVP 0 1.
где R(H,V)P - коэффициенты отражения Френеля на горизонтальной и вертикальной поляризациях; - некоторый геометрический множитель.
( ) Множитель в соотношении (1) связан с геометрической кон( ) фигурацией облучаемого участка поверхности, а диагональные элементы матрицы рассеяния представляют собой коэффициенты Френеля.
В работе получено соотношение, непосредственно связывающее между собой коэффициенты отражения Френеля на различных поляризациях:
RVP cos2 - RHP f = =, (2) RHP 1- RHP cos2 где комплексная величина f носит название поляризационного отношения.
Выражение (2) имеет значение для определения различных характеристик объектов мониторинга при анализе поляризационного состояния отраженных электромагнитных волн и связывает коэффициенты отражения Френеля на различных поляризациях и угол, под которым наблюдается поверхность. Таким образом, необходимость в проведении измерений на вертикальной поляризации нет, поскольку элементы матрицы рассеяния на любой поляризации могут быть определены путем прямого расчета.
Кроме элементов матрицы рассеяния, среди поляризационных характеристик можно отметить степень поляризационной анизотропии и разность аргументов между элементами матрицы рассеяния :
sin2 при Br, cos - sin2 = (3) cos - sin2 при >Br, sin2 RHP sinsin2 2 tg =, (4) RHP 1+ cos2 cos - RHP () (1+ )cos2 где Br - угол Брюстера; = arg RHP.
Проведенный в первой главе работы анализ показывает, что для широкого класса объектов мониторинга оказывается справедливым выражение (1) для матрицы рассеяния в виде произведения двух сомножителей. Как показывает выражение (2) для определения комплексной диэлектрической проницаемости достаточно проводить только относительные измерения коэффициентов отражения Френеля (элементов матрицы рассеяния) на различных поляризациях.
Полученные в первой главе результаты служат основой для применения радиополяриметрии при дистанционном определении электрофизических и геометрических характеристик объектов мониторинга.
Во второй главе проанализированы возможности радиополяриметрии для визуализации, классификации и идентификации объектов мониторинга, а также предложены новые методы, позволяющие определять электрофизические характеристики природных сред и техногенных объектов путем анализа рассеянного и собственного излучения этих объектов.
При использовании радиополяриметрии для дистанционного зондирования объектов мониторинга необходимо решить две фундаментальные задачи: прямую задачу, т. е. определить характеристики отраженной электромагнитной волны в зависимости от вида, типа и состояния объектов мониторинга; обратную задачу, когда по характеристикам отраженной электромагнитной волны определяется тип, вид и состояние объектов мониторинга. При мониторинге природных сред и техногенных объектов наибольший интерес представляет решение обратной задачи, которая относится к классу так называемых некорректных задач, при этом конечной целью мониторинга является определение электрофизических (комплексная диэлектрическая проницаемость и проводимость) и геометрических (форма и размеры) параметров объектов.
Определение электрофизических характеристик связано с анализом характеристик поляризационного состояния рассеянного и собственного излучения объектов мониторинга, которые непосредственно связаны комплексной диэлектрической проницаемостью исследуемых объектов. Результатом этого анализа является установленная взаимосвязь между поляризационными характеристиками отраженных электромагнитных волн и комплексной диэлектрической проницаемостью объектов мониторинга.
Поляризационное отношение, представленное в выражении (2), дает возможность выразить коэффициенты отражения RHP, RVP, т.е. диаго( ) нальные элементы матрицы рассеяния S11 и S22, через параметр f :
cos2- f S11 = RHP =, 1- f cos2 (5) cos2- f S22 = RVP = f.
1- f cos2 В этом случае степень анизотропии:
1- f 1+ f при Br, = (6) 1+ f при >Br, 1- f Анализ выражений (5) показывает, что путем измерений ортогональных составляющих матрицы рассеяния можно найти коэффициенты отражения Френеля, а, следовательно, и комплексную диэлектрическую проницаемость объектов мониторинга. Для этого необходимо вычислить отношения амплитуд и разность фаз таких составляющих.
В работе показано, что путем проведения относительных измерений возможно определять значение комплексной диэлектрической проницаемости:
1+ 4 f sin2 tg2.
= (7) 1- f ( ) Для оценки возможностей практического использования формулы (7) в работе проведена оценка относительной погрешности определения комплексной диэлектрической проницаемости, которая показала принципиальную перспективность предлагаемого подхода.
Соотношение (7) является также основой для решения задачи визуализации, классификации и идентификации объектов мониторинга. Для этого выразим комплексную диэлектрическую проницаемость непосредственно через параметры эллипса поляризации (коэффициент эллиптичности - r, угол наклона эллипса поляризации - ) отраженной электромагнитной волны:
1- r =1+ 1+ e4itg2sin2 r (8) Re= 1- r cos4tg2sin2 1+ 1+ r При использовании выражений (7) и (8) необходимо обеспечить последовательное излучение электромагнитных волн одинаковой мощности на горизонтальной и вертикальной поляризациях, при этом каждый раз измерять значение соответствующей ортогональной составляющей, после чего, необходимо определить параметр f.
С другой стороны, если осуществить излучение линейно поляризованной электромагнитной волны, то можно наблюдать аналогичную ситуацию.
Сигналы ортогональных составляющих могут быть использованы для определения поляризационного отношения f, который в конечном случае полностью описывает эллипс поляризации отраженной электромагнитной волны.
Следовательно, каждому значению f будет соответствовать некоторая поляризация, которая, как известно, взаимно однозначно отображается соответствующей точкой на сфере Пуанкаре. С другой стороны, параметр f однозначно определяет комплексную диэлектрическую проницаемость, т.е. вид объекта мониторинга. Следовательно, каждому виду объектов мониторинга можно поставить во взаимно-однозначное соответствие некоторую точку на сфере, каждая точка будет определять тот или иной вид объекта. Эта сфера носит название KLL-сферы (в соответствии с первыми буквами фамилий ее авторов - А.И. Козлова, А.И. Логвина, Л.П. Лихарда) (рис.1).
Рис.1. KLL-сфера С помощью KLL-сферы открываются новые возможности для классификации объектов мониторинга, при этом основным классификационным признаком будет выступать вид поляризации при облучении объекта электромагнитной волной, имеющей линейную поляризацию.
В работе получены количественные соотношения, позволяющие связывать параметр f с координатами точки на поверхности KLL-сферы (долгота и широта).
Особую роль при использовании предлагаемого подхода имеет угол наблюдения . Если угол наблюдения = 0, то в этом случае горизонтальная и вертикальная поляризации для любого вида объекта мониторинга не различимы между собой, при этом f =1. Данному значению f соответствует линейная поляризация с углом 45. При этом все объекты отобразятся в одну единственную точку (точка E). Поверхности в этом случае будут не различимы. Во втором предельном случае = 2, для любого ( ) вида объекта параметр f =-1. В этом случае все объекты также отобразятся в одну точку (точка F). Точка F диаметрально противоположна точке E. Она соответствует линейной поляризации с углом 135. При изменении угла наблюдения точки, соответствующие различным видам объектов, начнут описывать собственные траектории, каждая из которых начинается в точке E и заканчивается в точке F, при этом объекты мониторинга, представляющие собой диэлектрики отображаются на экваторе KLL-сферы.
Для иллюстрации возможностей применения KLL-сферы для визуализации, классификации и идентификации объектов мониторинга в работе представлены зависимости географических координат точек на поверхности KLL-сферы, отображающих различные типы земных подстилающих покровов. Пример таких зависимостей представлен на рис.2. Как видно из рис.2, для каждого типа подстилающих покровов на KLL-сфере будут вычерчиваться некоторые кривые линии, исходящие из точки, соответствующей координатам 2 = 0, 2 = 90, и заканчивающиеся в точке с координатами 2 = 0, 2 = 90.
Для иллюстрации выбраны поверхности, для которых отношение действительной части комплексной диэлектрической проницаемости к ее мнимой части было бы существенно меньше 1, существенно больше 1 или составляло бы величину одного порядка.
Существующие дистанционные методы определения связаны с непосредственными измерениями коэффициентов отражения и сопряжены с рядом принципиальных трудностей, вызванных необходимостью проведения абсолютных измерений. В результате коэффициент отражения и, следовательно, комплексная диэлектрическая проницаемость, могут быть определены с существенной ошибкой.
В работе предложен новый метод - амплитудно-фазовый, лишенный указанных недостатков, использующийся для определения комплексной диэлектрической проницаемости при наличии полной информации о поляризационном состоянии электромагнитных волн (информации об амплитудах и разностях фаз ортогональных составляющих).
Рис.2. Зависимость географических координат 2, 2 точек на KLLЦсфере для раз( ) личных типов подстилающих поверхностей от угла наблюдения:
Х - увлажненный песок; - морская вода; - кукуруза Относительные измерения ортогональных составляющих позволяют определять как абсолютные значения коэффициентов отражения, так и действительную и мнимую части комплексной диэлектрической проницаемости:
(1 + f )cos - 2 f sin tg2, Re = + 4 f (1 - 2 f cos + f ) (9) (1 - f )sin sin tg2, Im = 4 f (1 - 2 f cos + f ) где f и - модуль и аргумент поляризационного отношения.
В работе в качестве примера проводится расчет поляризационного отношения для зеленой травы высотой 80 см в X-диапазоне при различных условиях наблюдения. Полученные данные показали, что при дистанционном зондировании исследуемой поверхности под различными углами можно варьировать точностью и применимостью использованной модели и измерительной системы; в указанном диапазоне частот затруднительно получить точные результаты для значений разности фаз ортогональных составляющих электромагнитной волны.
Измерения разности фаз ортогональных составляющих представляют собой самостоятельную, технически весьма сложную задачу, в то время как измерения отношения амплитуд этих составляющих могут быть выполнены с достаточно высокой степенью точности. С другой стороны не представляет особых трудностей проведение амплитудных измерений на разных углах зондирования.
В работе показано, что и в условиях отсутствия необходимой информации о поляризационном состоянии электромагнитных волн возможно определять комплексную диэлектрическую проницаемость объектов мониторинга. Для этого используется модификация предлагаемого в работе амплитудно-фазового метода, только с учетом информации об отношении амплитуд ортогональных составляющих электромагнитной волны.
Сущность амплитудного метода определения комплексной диэлектрической проницаемости заключается в следующем. Любому значению аргумента (согласно (9)) будет соответствовать точка на плоскости с координатами Im,Re. При изменении аргумента от нуля до 360 получа( ) ется замкнутая кривая, каждая точка которой определяет значения действительной и мнимой частей комплексной диэлектрической проницаемости при определенном значении аргумента . Аналогичные кривые могут быть получены для всех значений угла наблюдения. Точка пересечения таких кривых соответствует реальному значению комплексной диэлектрической проницаемости. Координаты этой точки могут быть также получены в результате решения системы соответствующих уравнений.
Для подтверждения возможностей амплитудного метода в работе проведены расчеты действительной и мнимой частей комплексной проницаемости путем обработки экспериментальных данных для следующих типов подстилающих поверхностей: зеленая трава, сухой и мокрый снег, поле без растительности, деревья. Примером такого расчета может служить номограмма для Im и Re зеленой травы, представленная на рис.3.
Рис.3. Номограмма для определения Im и Re зеленой травы Кроме рассеянного радиоизлучения, информация о физических, электрофизических и геометрических параметрах объектов мониторинга содержится в собственном радиоизлучении этих объектов. Прием сигналов радиотеплового излучения, позволяет определять температуру, влажность, оценить уровень шероховатости, измерения обнаружить неровности и возвышенности на поверхности объектов мониторинга. Учет поляризационных параметров радиотеплового излучения повышает информативность проводимых дистанционных исследований.
Технологии оценки собственного и рассеянного излучения объектов мониторинга связаны между собой. В работе приведены основные соотношения, связывающие между собой излучательную способность (радиояркостную температуру) природных сред с их физическими характеристиками и комплексной диэлектрической проницаемостью. Для определения комплексной диэлектрической проницаемости объектов мониторинга при микроволновой радиополяриметрии, как и при рассеянном радиоизлучении, возможно использование методов решения прямых и обратных задач.
Для описания поляризационного состояния собственного радиоизлучения в работе использована матрица потерь. Элементы матрицы потерь косвенно несут информацию о свойствах объектов мониторинга и представляют собой, с одной стороны, излучательные способности на различных поляризациях, а с другой стороны - координаты некоторых точек на поверхности сферы Пуанкаре.
Матрица потерь имеет следующий вид:
1+ cos2 - e2i sin 2 L = П , (10) - e-2i sin 2 1- cos2 где - полная излучательная способность; и - сферические координаты на сфере Пуанкаре.
Принципиальное значение представления матрицы потерь в виде (10) состоит в том, что оно содержит непосредственно параметры исследуемого объекта. В работе показан физический смысл элементов матрицы потерь.
Верхний диагональный элемент определяет излучательную способность на поляризации 2 и 2, нижний равен излучательной способности на ортогональной поляризации. Максимальное значение излучательной способности достигается на собственной поляризации (2 = 0 ):
1 max = П 1+ , (11) ( ) а минимальное - на ортогональной:
2 min = П 1- . (12) ( ) При изменении поляризации на сфере Пуанкаре так, что соответствующая ей точка перемещается по параллели = const, излучательная способность не изменяется, при перемещении по меридиану от Северного к Южному полюсу ( = const ) она монотонно убывает от max до min.
Экстремальные значения величины , как видно из (11), являются собственными значениями матрицы потерь. Из равенств (11) и (12) можно выразить степень поляризационной анизотропности через экстремальные значения излучательной способности:
= 1 -2 1 +2. (13) ( ) ( ) В работе показано, что для построения матрицы потерь, т. е. проведения полного поляризационного анализа, достаточно измерений на любых четырех поляризациях, УобразыФ которых не лежат на одном большом круге сферы Пуанкаре. Такими поляризациями могут быть горизонтальная, вертикальная, линейная с углом наклона 45 по отношению к горизонтальной и круговая.
Полностью определить матрицу потерь можно путем определения четырех радиояркостных температур на указанных выше поляризациях.
Этот этап соответствует решению прямой задачи. Обратная задача подразумевает определение электрофизических и геометрических характеристик объектов мониторинга по имеющимся значениям радиояркостных температур. Под неизвестными геометрическими параметрами понимают угол наклона поверхности объекта мониторинга и угол скрутки.
В работе показано, что при известном угле визирования диэлектрическая проницаемость определяется по формуле = 1+ 1-Г -2 cos2 + sin2. (14) ( ) Г Угол визирования и диэлектрическую проницаемость находят следующим образом. Берется какое-то значение 0, и из формул Френеля определяются соответствующие значения cos для вертикальной и горизонтальной поляризаций. Если 0 задана с погрешностью, то значения углов визирования, вычисленные по этим формулам, будут различаться; при этом разность углов в первом приближении пропорциональна отклонению 0 от истинного значения диэлектрической проницаемости и (рис.4). На рис.5 0 соответствует истинному значению угла.
Таким образом, знание матрицы потерь, определяющей радиотепловое излучение поверхности объекта мониторинга с известной температурой, позволяет найти ориентацию исследуемой поверхности и диэлектрическую проницаемость, которая предполагается вещественной. Приведенные соотношения относились к случаю поляризационного измерения собственного излучения с одного направления. В работе получены аналогичные результаты для случая поляризационных измерений с нескольких направлений. Такие измерения позволяют получить более подробную информацию об объекте мониторинга, например, определить мнимую часть комплексной диэлектрической проницаемости.
Рис.4. Зависимость излучательных способностей на вертикальной и горизонтальной поляризациях от угла :
1 - истинное значение и ; 2 - первоначальное приближение и В третьей главе проанализированы возможности радиополяриметрии при дистанционном определении геометрических характеристик. Проводится анализ объектов мониторинга с геометрическими неоднородностями и особенности определения геометрических характеристик пространственно-протяженных объектов дистанционного зондирования.
Источником информации об геометрических параметрах объектов мониторинга служит как собственное, так и рассеянное радиоизлучение объектов мониторинга. В качестве геометрических характеристик объектов мониторинга в работе рассматриваются степень шероховатости (при анализе рассеянного радиоизлучения) и угол визирования, угол скрутки (при анализе собственного радиоизлучения). В общем случае провести разделение факторов, связанных с поляризационным состоянием электромагнитных волн, на зависящие от электродинамических характеристик и от геометрической конфигурации объектов не представляется возможным.
Однако, как показано в работе, для достаточно широкого класса подстилающих поверхностей амплитуду отраженной электромагнитной волны можно представить в виде произведения двух множителей, один из которых определяется только геометрическими параметрами, а второй - коэффициентами отражения Френеля.
Как частный случай, в работе рассматриваются возможности радиополяриметрии при определении геометрической конфигурации объектов мониторинга, имеющих протяженные линейные размеры с точки зрения дистанционного наблюдения. Показано, что путем управления поляризационным состоянием электромагнитных волн возможно оценивать линейные размеры таких объектов.
В общем случае поверхность реальных объектов мониторинга представляет собой статистически неровную двумерно шероховатую поверхность z = x, y с диэлектрической проницаемостью x, y. Как и в слу( ) ( ) чае гладкой однородной поверхности, необходимо определить элементы матрицы рассеяния таких объектов, т.е. найти решение прямой задачи дистанционного зондирования.
В работе проводится анализ четырех наиболее распространенных моделей неровных (шероховатых) поверхностей (рис.5):
1.крупномасштабная, плавная, пологая, с произвольной шероховатостью, в среднем плоская поверхность;
2.мелкомасштабная, пологая, слабошероховатая, в среднем плоская поверхность;
3.сложная поверхность, покрытая мелкомасштабными и крупномасштабными неоднородностями;
4.модель в виде некогерентных независимых рассеивателей.
Рис.5. Модели неровных поверхностей:
а - модель 1; б - модель 2; в - модель 3; г - модель Параметрами неоднородностей поверхности выступают радиусы корреляции и кривизны поверхности и среднеквадратичная высота неровностей. С помощью перечисленных моделей могут быть описаны реальные поверхности природных сред и техногенных объектов.
При помощи решения задачи дифракции электромагнитной волны на поверхности в приближении Кирхгофа и использовании теоремы Грина, в работе получены выражения для матриц рассеяния трех перечисленных выше моделей. Матрицу рассеяния для четвертой модели из-за общего характера структур, описываемых этой моделью, представить в явном виде затруднительно. Матрицы рассеяния рассматриваемых моделей неровных поверхностей имеют вид:
a2nz a1a2nz () () 1- a2nz 2RVP - a1 RHP 1- a2nz RHP + RVP S1 =, (15) aa2nz a2nz RHP + RVP a1 RVP - 2RHP () () 1- a2nz 1- a2nz 1+ RVP L , () ( ) S2 =, (16) 01+ RVP 1- RHP sin2 + L , () () ( ) () 2q2nz 1+ a2nz q - a2nz 2a1qnz -2 a2nz + q a2nz + q nz a2nz - q () 2a2nz,(17) S3 = a2nz + q () 2a1qnz 2qa1 nz 1+ a2nz q - a2nz - 2 a2nz + q a2nz + q nz a2nz - q2 () где a1,a2,nz,,q - параметры.
Кроме этого, для каждой из трех моделей определены координаты (широта и долгота) собственного поляризационного базиса на поверхности сферы Пуанкаре. Например, для модели 1 координаты собственного поляризационного базиса определяются только характеристиками шероховатости поверхности объекта. Наличие мелкомасштабных неровностей на поверхности объекта (модель 2) приводит к тому, что координаты собственного поляризационного базиса равны нулю. Совокупность мелкомасштабных и крупномасштабных неоднородностей на поверхности объекта (модель 3) определяет однозначную зависимость координат собственного поляризационного базиса от характеристик шероховатости.
Значительная часть техногенных объектов мониторинга, исследуемых в работе, с точки зрения дистанционного наблюдения могут считаться протяженными объектами. Такие объекты мониторинга удобно представить в виде набора пространственно разнесенных блестящих точек, которые являются достаточно чувствительными к изменению поляризационного состояния облучаемой электромагнитной волны. Это означает, что такие участки поверхности объектов мониторинга меняют свои отражательные характеристики при изменении вида поляризации облучаемой электромагнитной волны. При изменении поляризации источник отраженной электромагнитной волны (фазовый центр) будет перемещаться в пространстве по некоторой кривой. Конфигурация и размеры такой пространственной кривой связана с пространственным расположением блестящих точек.
В работе введено понятие пространственной кривой в виде пространственного поляризационного контура и предложена процедура его построения - поляризационное оконтуривание.
В работе показано, что линейные размеры пространственного поляризационного контура связаны с геометрическими характеристиками объектов мониторинга, что дает возможность оценивать эти характеристики.
В связи с предлагаемой процедурой поляризационного оконтуривания, в работе рассматриваются два случая: некогерентное и когерентное рассеяние электромагнитных волн от поверхности объектов мониторинга.
При описании процедуры поляризационного оконтуривания в работе рассматривается модель объекта мониторинга в виде набора из N блестящих точек с пространственными координатами xi, yi. Матрица рассеяния ( ) i-й блестящей точки имеет вид: Si = Smni.
При некогерентном рассеянии электромагнитных волн суммарная мощность всех блестящих точек представляет собой сумму мощностей электромагнитной волны от каждой блестящей точки. В свою очередь мощность для каждой блестящей точки пропорциональна значению s11i.
Для решения задачи оценки геометрических характеристик в работе используется следующий алгоритм:
1. Определяются пространственные координаты фазового центра отраженной электромагнитной волны:
NN N N Xфц = Pi i i s11i, (18) x P = x s11i i i=1 i=1 i=1 i=NN N N Yфц = yiPi i yi s11i s11i. (19) P = i=1 i=1 i=1 i=2. При изменении поляризационного состояния электромагнитной волны матрица рассеяния Si подвергается конгруэнтному преобразованию с помощью унитарной матрицы. По этой причине вместо s11i будет иметь место:
q11 = s11 e-2i cos2 + s22 e2i sin2 + s12 sin 2, (20) где , - параметры, определяющие конкретный вид поляризации, на которой производится наблюдение объекта мониторинга.
3. Строится кривая, соответствующая перемещению положения фазового центра в пространстве (пространственный поляризационный контур).
В работе проведен анализ изменения формы пространственного поляризационного контура объекта мониторинга в зависимости от числа блестящих точек и их отражательных характеристик.
Примеры пространственного поляризационного контура приведены на рис.6.
а) б) в) Рис.6. Пространственный поляризационный контур:
а - для случая 4 блестящих точек; б, в - трансформация при изменении параметра Изменение числа блестящих точек и их отражательных характеристик (параметров и ) приводит к деформации исходного контура при сохранении геометрических размеров контура.
Для случая когерентного рассеяния электромагнитных волн объект мониторинга представляет собой совокупность большого числа статистически независимых между собой отражателей. Такая система отражателей вполне пригодна для описания большей части природных сред и техногенных объектов.
По аналогии с формулами (18), (19) и с учетом рассматриваемой модели, выражения для средних значений координат фазового центра отраженной электромагнитной волны имеет вид:
NNN N N N п Xфц == x Pi P xiPi P =x Pi P = 0, (21) i i j i j i=1 i=1 i=1 j=1 i=1j=NN N N N N п Yфц = yiPi i = yi Pi j P yiPi P = P = 0. (22) j i=1 i=1 i=1 j=1 i=1j= В работе для рассматриваемой модели когерентного отражения получены соотношения для среднего квадратичного отклонения положения фазового центра от своего среднего значения (величины X фц и Y фц ). Показано, что это отклонение от геометрического центра системы независимых отражателей будет:
X фц = Y фц 0,4L. (23) Для проверки правильности полученных теоретических соотношений в работе было проведено моделирование. Объект мониторинга представлялся в виде системы из 100 отражателей, случайным образом распределенных на участке поверхности объекта размером 2L 2L. Каждый отражатель характеризовался излучательной способностью, которая является случайной величиной с известным законом распределения. С помощью формул (18), (19), (21), (22), (23) определялось положение фазового центра и средние квадратические отклонения его положения. По результатам моделирования производились оценки соответствия полученных размеров пространственного поляризационного контура и реальных геометрических характеристик объектов мониторинга.
Для получение целостной картины при оценки геометрических характеристик объектов мониторинга путем управления поляризационным состоянием электромагнитной волны в работе рассмотрены результаты поляризационного оконтуривание различных моделей пространственно протяженных объектов дистанционного наблюдения. Были рассмотрены модели: комбинированная (случай некогерентного и когерентного рассеяния), двухточечная, трехточечная и многоточечная.
В качестве иллюстрации проведенных исследований на рис.7 изображены положения фазового центра при моделировании, а на рис.8 - влияние вида поляризации на положение среднего значения фазового центра и вид поляризационного контура.
Рис.7. Положение фазового центра при моделировании (сплошная линия - пространственный поляризационный контур; маркер - положения блестящих точек) а) б) Рис.8. Поляризационное оконтуривание (маркер - положения блестящих точек):
а - влияние поляризации на положение фазового центра;
б - определение среднего положения фазового центра Задача определения параметров крупных геометрических неоднородностей на поверхности объектов мониторинга может быть успешно решена путем анализа поляризационного состояния собственного радиоизлучения этих объектов. Для анализа возможностей микроволновой радиополяриметрии для оценки геометрических размеров неоднородностей, в диссертационной работе используются ранее введенные модели неровных (шероховатых) поверхностей.
Основными геометрическими характеристиками при анализе собственного радиоизлучения являются угол скрутки 0 и угол визирования .
Угол скрутки представляет собой угол между плоскостью условной вертикальной поляризации и плоскостью, перпендикулярной к поверхности и проходящей через направление наблюдения. Угол - угол между направлением наблюдения и нормалью к поверхности.
При дистанционном наблюдении объектов мониторинга с использованием анализа собственного радиоизлучения используют измеренные значения излучательной способности на трех поляризациях:
В = 0,5 1 +2 + 0,5 1 -2 cos 2cos20, ( ) ( ) Г = 0,5 1 +2 - 0,5 1 -2 cos2cos20, (24) () () Н = 0,5 1 +2 - 0,5 1 -2 cos2sin 20.
() () где 0 - угол скрутки.
В работе показано, что наличие неоднородностей на поверхности объекта мониторинга приводит к изменению радиояркостных температур на вертикальной и горизонтальной поляризации по сравнению с теми же параметрами для гладкой поверхности. Таким образом, в работе удобно использовать известный параметр шероховатости S, который определяется только разностями указанных температур для неровной и гладкой поверхностей, и носит название степени шероховатости. Использование степени шероховатости, с одной стороны, позволяет отказаться от использования моделей, определяющих тип поверхности, а с другой стороны, проводить их определенную классификацию. Классификация поверхностей связаны с различными значениями угла визирования. Кроме этого, параметр S связан с различными характеристиками геометрических неоднородностей.
В работе показано, что степень шероховатости S =2 2Г - 2 B -Г 1-Г cos2 -B 1-Г cos2 2 -B , (25) ( ) ( ) Г где Г и B - излучательные способности шероховатой поверхности на горизонтальной и вертикальной поляризациях, является количественным критерием неровности поверхности.
В работе также проведен анализ изменения величины S при изменении излучательных способностей на вертикальной и горизонтальной поляризации.
Этот анализ показал, что существует минимальное значение степени шероховатости, которое соответствует условию B =Г =0. Кроме этого, на значение параметра S оказывает влияние величина угла визирования.
Таким образом, проведенные в третьей главе исследования подтверждают возможности использования анализа поляризационного состояния рассеянного и собственного радиоизлучения объектов мониторинга для определения геометрических характеристик таких объектов. Для обоих видов радиоизлучения при определении геометрических характеристик используются адекватные модели неровных поверхностей объектов мониторинга. Геометрические характеристики протяженных объектов также могут быть оценены с высокой степенью точности путем реализации процедуры полного поляризационного сканирования. В результате этой процедуры строится пространственный поляризационный контур, геометрическая конфигурация которого связана с геометрическими характеристиками самих объектов. Проведение поляризационных измерений собственного излучения объектов наблюдения на горизонтальной и вертикальной поляризациях дает возможность определять геометрические параметры неоднородностей на поверхности этих объектов.
В четвертой главе рассмотрены особенности радиополяриметрии при нелинейном рассеянии поляризованных электромагнитных волн объектами мониторинга, определены поляризационные характеристики объектов мониторинга, обладающих нелинейными свойствами рассеяния.
Нелинейные эффекты при рассеянии поляризованных электромагнитных волн позволяют оценить различные параметры достаточно распространенного класса природных или техногенных объектов на гармониках или комбинационных составляющих излучаемой электромагнитной волн. Примером нелинейных эффектов могут служить процессы распространения электромагнитных волн в ионосфере и рассеяния электромагнитных волн объектами, содержащими контакты металл-диэлектрик-металл с нелинейными вольтамперными характеристиками и полупроводниковые элементы.
Уравнение движения простейшего негармонического осциллятора имеет следующий вид:
d x t dx t ( ) ( ) me 2 (26) + + me0x t +x2 t = qeE cos t -kr, ( ) ( ) ( ) dt2 dt где me, qe - масса и заряд осциллятора соответственно; - коэффициент, определяющий нелинейность; - коэффициент затухания; x t - смещение ( ) заряда. Если атом расположен в центре симметрии, то появляется нелинейность более высокого порядка из-за присутствия в последнем дифференциальном уравнении не только члена ax2 t, но и bx3 t, а также более ( ) ( ) высоких мультипольных членов. Все это приводит к установлению наведенных дипольных моментов на частотах 2, 3 и выше.
В работе приведена классификация объектов мониторинга, обладающих нелинейными свойствами. К таким объектам, в большей степени, относятся объекты техногенной природы.
Для оценки возможностей радиополяриметрии при определении различных характеристик объектов мониторинга в работе были рассмотрены возможные модели объектов, обладающих нелинейными свойствами. Основные результаты такого анализы были получены для двухвибраторной модели объекта мониторинга - нелинейного рассеивателя.
Рассматриваемая модель представляется линейно-поляризованной антенной. Напряжение, индуцируемое падающей электромагнитной волной, определяется его амплитудной характеристикой:
n sc U = an Ein, (27) где an - постоянный коэффициент; n - номер гармоники.
Электромагнитная волна, рассеянная нелинейным объектом, имеет следующий вид:
in usc(t) = f ( E )cos(ntin + nin +n+ nkinr +n), (28) где n - величина, принимающая значение 0 или 1 в зависимости от знака соответствующего коэффициента разложения вольт-амперной характеристики нелинейного объекта; n фазовая добавка, связанная с инерционными свойствами нелинейного элемента; nkinr и n - фазовые слагаемые.
Поляризационные характеристики объекта мониторинга, обладающего нелинейными свойствами, рассматриваются в предположении изменении поляризационного состояния облучаемой электромагнитной волны. Рассеянный электромагнитная волна будет определяться компонентами sс uГ (y) = U cos(nint + nin + Г), Г (29) sс uВ (y) = U cos(nint + nin + В).
В Изменения одного из коэффициентов Г или приводит к тому, что В поляризационное состояние рассеянной электромагнитной волны изменится на ортогональное.
Анализ двухвибраторной модели нелинейного рассеивателя показывает, что можно получить необходимую поляризацию рассеянной электромагнитной волны при заданной поляризации облучаемой. При этом проявляются следующие особенности, связанные с нелинейным характером рассеяния:
при смене поляризации облучаемой электромагнитной волны не сохраняется величина интенсивности рассеянной волны;
на поляризацию рассеянной волны оказывает существенное влияние вид амплитудной характеристики модели;
для четных гармоник фаза рассеянной волны зависит от ориентации модели в пространстве, поэтому на поляризацию рассеянной волны существенное влияние оказывает ориентация образующих модель вибраторов;
исчезает возможность однозначного определения поляризационного состояния облучаемой электромагнитной волны по поляризации рассеянной, и существует ограничение на некоторые виды поляризации рассеянной электромагнитной волны при заданной поляризации облучаемой электромагнитной волны и фиксированных параметрах модели нелинейного объекта.
По аналогии с классическим (линейным) случаем радиополяриметрии при дистанционном наблюдении объектов мониторинга в работе определяется основная поляризационная характеристика - матрица рассеяния.
При этом квадраты модулей для элементов сформированной матрицы рассеяния будут выступать аналогом ЭПР на частоте рабочего нелинейного продукта, а аргументы этих элементов должны отсчитываться от фазы того же нелинейного продукта, но полученного путем преобразования в специальном (эталонном) устройстве. Сформированная таким образом нелинейная ЭПР становится функцией от величины плотности потока мощности излучаемой электромагнитной волны, падающей на нелинейный рассеиватель, и не имеет физического смысла без указания того значения мощности, которое в настоящий момент определяет поддержание процесса нелинейного рассеяния.
В работе проведен анализ существующих методов описания объектов мониторинга, обладающих нелинейными свойствами: электродинамические методы, метод модулирующих функций и функциональные методы.
Среди предлагаемых методов наиболее целесообразно использовать функциональные методы, использующих ряды ВольтерраЦВинера:
y(t) = S0 t + S1 t, u d + S2 t,1,2 u 1 u 2 d1d2 +...+ ( ) ( ) ( ) ( ) ( ) ( ) - - - + Sn t,1,2,...,n u 1 u 2 u n d1d2 dn, (30) () ( ) ( ) ( ) - - где u(t) - сигнал, действующий на входе исследуемого объекта; у(t) - отклик объекта; S1(.), S2(.), Е, Sn(.), - ядра функционального ряда. С помощью функциональных рядов Вольтерра-Винера решается задача разделения видов нелинейности для последующей классификации и устанавливается относительно простая связь используемых на практике критериев нелинейности с характеристиками ядер S2(.), Е, Sn(.), Е ряда (30).
Для получения соотношений для элементов матрицы рассеяния объектов мониторинга, обладающих нелинейными свойствами, в работе были рассмотрены процессы рассеяния плоских волн на нелинейном рассеивателе и стационарное нелинейное рассеяния электромагнитных волн. Для этого облучающая объект мониторинга поляризованная электромагнитная волна может быть представлена в виде совокупности задержанной на различное время копий одной и той же электромагнитной волны:
in Esc r,t = S r,t E t - ,r d + ( ) ( ) ( ) - in + S r,1,2 E t - 1,r Ein t - 2,r d1d2 +...+ () () () - - in + S... r,1,2,...,n E t -1,r Ein t -n,r d1d2 dn. (31) () () () - - В предположении отсутствия временной дисперсии S(r,t) = S(r)(t), S(r,t1,t2) = S (r)(t1)(t2) и т.д., ряд (31) становится алгебраическим:
in in Еsc(r,t) = S(r)E (r,t) + S (r)E (r,t)Ein (r,t). (32) В случае применяемой в теории рассеяния нормировке (падающая волна - плоская, размерность напряженности - Увольт на метрФ, рассеянная, будучи диаграммой рассеяния имеет размерность - вольт), размерность элементов матрицы S - УметрФ, S - Уметр квадратный на вольтФ, и вообще, для n-мерного ядра - Уметр в степени (n-1) на вольт в степени (n2)Ф. Такие размерности обеспечивают независимость введенных характеристик рассеивающей способности от поляризационного состояния падающей волны, в том смысле, что эти ядра, как и линейная матрица рассеяния S, в совокупности определяют теперь уже нелинейный оператор рассеяния (функтор).
В работе рассматривается случай стационарного нелинейного рассеивания электромагнитных волн, поэтому для облучающей электромагнитной волны разложение имеет вид:
Ein (r,t) = Еin (,k)еi(t-kr)ddk - и аналогично для рассеянной волны справедливо выражение Esc(r,t) = Еsc(,k)еi(t-kr)ddk.
- Комплексный вектор плоской волны для рассеянного поля, определяемый в (6) вкладом ядра S (.), получается равным:
in Esc(,k) = S (, )Е ( - ,k-k )Еin (,k )d dk, (33) - где S ( ) = S (1,2)еi( + ) d1d2.
, - Из (33) следует, что в рассеянном электромагнитном поле появляется составляющая суммарной частоты, распространяющаяся в направлении, определяемом суммой волновых векторов составляющих плоских волн падающего электромагнитного поля. В работе показано, что и для составляющей нелинейности более высокого порядка соотношения между параметрами, воздействующими на рассеиватель плоских волн и параметрами рассеянного поля, имеют такой же характер.
На основе анализа реакции объекта мониторинга, обладающего нелинейными свойствами, на облучение его эллиптически поляризованной электромагнитной волной, в работе определяется выражение для матрицы рассеяния такого объекта наблюдения. Такое выражение может быть получено путем анализа преобразования поляризационного базиса рассеивателем вида:
N m Esc = Am Ein. (34) ( ) m=В работе показано, что для описания рассеивающих свойств нелинейного рассеивателя в общем случае требуется 8 комплексных чисел, т.е.
16 действительных. Эти числа образуют прямоугольную матрицу размером 4х2.
Выражение для элементов матрицы рассеяния объектов мониторинга, обладающих нелинейными свойствами, получено путем применения матрицы преобразования вида:
Q(2) bC,bH = e2i cos2 0,5e2i(+) sin 2 0,5e2i(+) sin 2 e2i(+2) sin2 -2i sin 2 cos2 -sin2 0,5e2i sin 2 = e2i -0,5e.
-0,5e-2i sin 2 -sin2 cos2 0,5e2i sin 2 e-2i(+2) sin2 -0,5e-2i(+) sin 2 -0,5e-2i(+) sin 2 e-2i cos2 В пятой главе диссертационной работы проведено описание экспериментальных исследований по определению электрофизических характеристик и геометрической конфигурации некоторых объектов мониторинга.
Рассматриваются результаты экспериментальных исследований с точки зрения их соответствия предлагаемых в работе методов.
Одна часть проведенных экспериментальных исследований касалась определения поляризационных характеристик широкого класса подстилающих покровов. Эксперимент проводился на трассовом радиолокаторе УСкалаФ (рис.9), работающий в дециметровом диапазоне волн. При проведении эксперимента поляризация излучаемой электромагнитной волны последовательно изменялась с горизонтальной до вертикальной, принимая значения 0, 30, 45, 60 и 90. Электромагнитная волна, отраженная от подстилающей поверхности, поступала на приемное устройство. Последующая обработка отраженных электромагнитных волн позволила определить элементы матрицы рассеяния подстилающих поверхностей и местных предметов.
Для определения поляризационных характеристик в процессе эксперимента были получены соответствующие реализации отраженной электромагнитной волны. Выведены соотношения, связывающие между собой амплитуду отраженной электромагнитной волны и параметры матрица рассеяния.
Элементы матрицы рассеяния в любом поляризационном базисе по сравнению с исходным базисом (базис Угоризонтальная - вертикальная поляризацииФ) определяются основным уравнением для элементов матрицы рассеяния:
e-2iP11 = e2i(+) S11e-2i cos2 + S22e2i sin2 + S21sin 2, (35) ( ) где Pij,Sij i =1,2 - элементы матрицы рассеяния в произвольном и ис( ) ходном базисах соответственно; и - параметры поляризационного базиса. Например, для значения = 45 будем иметь:
22 2 2 P11 45 = S11 + S22 + S12 + S11 S22 cos 11 - 22 + () 44 2 (36) + S11 S12 cos 11 - 12 + S22 S12 cos 22 - 11 ;
() () Рис.9. Радиолокационная станция УСкалаФ В результате первичной обработки экспериментальных данных были получены реализации для следующих параметров электромагнитной волны: элементы матрицы рассеяния S11, S12, S22 ; фазы этих элементов 11, 12, 22 и их разности; собственные значения матрицы рассеяния 1 и 2.
На следующем этапе были определены законы распределения этих случайных величин. Каждый параметр вычислялся трижды по трем значениям объема выборки на азимутальных углах 0 - 360, 0 -100 и 280- 360, 100- 280.Объемы выборки для рассматриваемых параметров составляли в среднем 1440.
На основе методики, предложенной в работе, были рассчитаны числовые характеристики статистических законов перечисленных выше параметров (в общей сложности 14 параметров): математическое ожидание, дисперсия, среднее квадратическое отклонение, асимметрия и эксцесс. В работе проведен анализ числовых характеристик различных поляризационных параметров. Этот анализ показал, что наиболее устойчивыми к объему выборки являются не амплитудные, а именно фазовые параметры и степень анизотропии. Наиболее устойчивыми к изменению азимутального угла являются координаты собственного поляризационного базиса ( 2 и 2 ) и степень анизотропии.
При аппроксимации статистического закона распределения рассматриваемых случайных величин следует выбирать такой закон, который учитывает большинство изменений статистических характеристик и параметров распределения. Это соответствует применению физикостатистического подхода к выбору типа аппроксимирующего распределения, который предполагает выбор предельного распределения поляризационных характеристик на основе физических условий их формирования и электродинамических характеристик подстилающих покровов. Изменение окружающих условий вызывает изменение параметров статистической модели, не изменяя закона распределения, основанного на физикостатистическом подходе. При обработке результатов эксперимента в качестве общего вида распределения в работе выбрано семейство кривых Пирсона, определяемое дифференциальным уравнением:
dy x + b = dx.
y C0 + C1x + C2xВ работе были построены гистограммы распределения элементов S11, S12, S22. Путем использования метода наименьших квадратов были построены плотности и функции распределения вероятностей этих параметров, аппроксимирующие эти гистограммы.
Анализ гистограмм и соответствующих им статистических распределений параметров матрицы рассеяния подстилающих покровов показывает, что измеренные поляризационные характеристики достаточно компактно группируются вокруг соответствующих средних значений. Гистограммы и аппроксимирующие их законы распределения показывают, что для элементов S11 и S22 характерно наличие двумодального, а иногда и трехмодального вида распределения. Резко отличается вид закона распределения для кроссового элемента S12 ; для него характерна унимодальность и мала ширина кривой плотности распределения. Аналогичные гистограммы и аппроксимирующие их распределения получены для соотношений фаз элементов матрицы рассеяния и углов, характеризующих положение собственного поляризационного базиса на сфере Пуанкаре, собственных элементов матрицы рассеяния. На рис.10 и рис. 11 изображены гистограмма и плотность (функция) распределения элемента S11, которые аппроксимируют эту гистограмму.
Одним из возможных подходов к описанию распределений на базе системы математических формул является способ, основанный на отыскании семейства кривых Пирсона, при помощи которого можно удовлетворительно представить встречающиеся распределения. В работе на основе известного алгоритма проведен расчет значения параметра согласования для всех 14 параметров элементов матрицы рассеяния, соответствующий тип кривых Пирсона, а также соответствующий им вид распределения.
При этом наиболее УраспространеннымФ является третий тип кривых Пирсона и соответствующее - распределение параметров. Для оценки соответствия эмпирического распределения теоретическому был использован критерий Колмогорова.
Рис.10. Гистограмма элемента SРис.11. Плотность и функция распределения элемента SВторая часть экспериментальных исследований проводилась с целью подтверждения возможности оценки геометрических характеристик объектов мониторинга, имеющих протяженные линейные размеры, с помощью процедуры поляризационного оконтуривания.
Экспериментальная установка (рис.12) включала в себя передающую и приемную части. В качестве передающих и приемных антенн использовались рупорные антенны (три передающих и одна приемная). Передающие рупоры, моделирующие блестящие точки реального объекта мониторинга, располагались в вершинах равнобедренного треугольника. Имитация изменения матриц рассеяния блестящих точек при изменении поляризации зондирующего сигнала производилась вращением передающих антенн вокруг своей оси на произвольный угол. На приемной стороне производилась настройка приемной антенны на максимум принимаемого сигнала по показаниям измерительного прибора с фиксацией угловых положений измерительной антенны.
а) б) Рис.12. Установка для проведения экспериментальных исследований возможностей поляризационного оконтуривания В результате проведенного эксперимента было получено свыше 1реализаций. Далее с помощью известных соотношений были рассчитаны координаты фазового центра принимаемой электромагнитной волны и построен пространственный поляризационный контур. Описанная процедура проводилась для случаев когерентного и некогерентного рассеяния электромагнитных волн.
На рис.13 и рис.14 представлены пространственные поляризационные контуры, построенные по результатам обработки экспериментальных данных. Первый рисунок соответствует случаю когерентного рассеяния электромагнитных волн, а второй - некогерентного рассеяния.
Рис.13. Поляризационный контур, полученный в результате обработки экспериментальных данных (случай когерентного рассеяния):
маркер - расположение антенн;
сплошная линия - пространственный поляризационный контур Рис.14. Поляризационный контур, полученный в результате обработки экспериментальных данных (случай некогерентного рассеяния):
маркер - расположение антенн;
сплошная линия - пространственный поляризационный контур Анализ полученных экспериментальных результат позволяет сделать вывод об адекватности предлагаемого подхода по оценки геометрической конфигурации объектов мониторинга путем управления поляризационным состоянием электромагнитных волн.
Форма и расположение пространственного поляризационного контура совпадают с аналогичными, полученными в результате теоретических расчетов.
Заключение В диссертационной работе на основе проведенного анализа:
- традиционных методов дистанционного определения физических, химических, электрических и геометрических характеристик природных сред и техногенных объектов мониторинга, - возможности и целесообразности применения радиополяриметрии при дистанционном определении электрофизических и геометрических характеристик природных сред и техногенных объектов мониторинга, - возможности применения радиополяриметрии для определения электрофизических характеристик объектов мониторинга, обладающих свойствами нелинейного рассеивателя, а также путем проведенного в работе - определения взаимосвязи электродинамических и поляризационных характеристик отраженных объектами мониторинга электромагнитных волн;
- обоснования возможности применения методов радиополяриметрии для определения электрофизических характеристик объектов мониторинга при их дистанционном зондировании;
- обоснования возможности применения методов радиополяриметрии для улучшения обнаружения техногенных объектов на фоне природных сред при их мониторинге;
- обоснования возможности применения методов радиополяриметрии для визуализации, классификации и идентификации природных средств и техногенных объектов установлена прямая функциональная связь между электродинамическими и поляризационными характеристиками объектов мониторинга, разработаны методы определения электрофизических и оценки геометрических характеристик этих объектов и разработаны математические модели неоднородных техногенных объектов, обладающих нелинейными свойствами, что соответствует поставленной цели диссертационной работы.
В ходе выполнения работы были получены следующие новые научные результаты:
- установлена прямая функциональная связь между электродинамическими (комплексная диэлектрическая проницаемость) и поляризационными характеристиками (элементы матрицы рассеяния) объектов мониторинга;
- предложен метод визуализации объектов мониторинга при помощи KLL-сферы, позволяющий реализовать их идентификацию и классификацию;
- разработан метод определения комплексной диэлектрической проницаемости широкого класса природных объектов на основе анализа поляризационного состояния отраженных этими объектами электромагнитных волн;
- экспериментально определены статистические характеристики элементов матрицы рассеяния некоторых классов подстилающих покровов;
- разработан метод оценки геометрических характеристик техногенных объектов, находящихся на фоне природных сред, путем управления поляризационными характеристиками электромагнитных волн;
- разработаны математические модели неоднородных техногенных объектов, обладающих нелинейными свойствами, при их взаимодействии с поляризованными волнами.
Полученные результаты позволяют:
- определять при дистанционном зондировании окружающей среды физические характеристики подстилающих покровов (соленость, влажность, плотность, температуру и т.п.);
- оценивать при дистанционном зондировании окружающей среды параметры геометрических неоднородностей;
- увеличивать степень различения техногенных объектов мониторинга, находящихся на поверхности природных сред;
- оценивать геометрические размеры объектов мониторинга.
Публикации по теме диссертационной работы Результаты исследований, представленные в диссертационной работе, опубликованы в следующих основных печатных работах автора.
Статьи в научных журналах, рекомендованных ВАК Минобразования России для публикации основных научных результатов диссертаций на соискание ученой степени доктора наук 1. Колядов Д.В. Анализ взаимосвязи между электрофизическими и электродинамическими характеристиками подстилающих покровов для решения задач дистанционного зондирования. - Научный вестник МГТУ ГА. Сер. Радиофизика и радиотехника, 2000, №24.
2. Колядов Д.В. Анализ влияния поляризационных характеристик целей на их различимость. - Научный вестник МГТУ ГА. Сер. Радиофизика и радиотехника, 2001, №36.
3. Колядов Д.В. Некоторые принципы классификации радиолокационных целей. - Научный вестник МГТУ ГА. Сер. Радиофизика и радиотехника, 2001, №36.
4. Колядов Д.В. Поляризационная обработка радиолокационных сигналов, отраженных от протяженных объектов. - Научный вестник МГТУ ГА. Сер. Радиофизика и радиотехника, 2002, №54.
5. Колядов Д.В. Связь между электродинамическими и поляризационными характеристиками подстилающих покровов. - Научный вестник МГТУ ГА. Сер.
Радиофизика и радиотехника, 2000, №24.
6. Колядов Д.В., Поляризационная обработка радиолокационных сигналов, отраженных от протяженных объектов. - Научный вестник МГТУ ГА. Сер. Радиофизика и радиотехника, 2002, №54.
7. Колядов Д.В. Амплитудный метод определения комплексной диэлектрической проницаемости подстилающих поверхностей при их дистанционном зондировании. - Научный вестник МГТУ ГА. Сер. Радиофизика и радиотехника, 2002, №54.
8. Козлов А.И., Колядов Д.В. Матрица рассеяния нелинейных отражателей.
Ц Научный вестник МГТУ ГА. Сер. Радиофизика и радиотехника, 2004, №79.
9. Козлов А.И., Колядов Д.В. Основное уравнение нелинейной радиолокации. - Научный вестник МГТУ ГА. Сер. Радиофизика и радиотехника, 2004, №79.
10. Козлов А.И., Колядов Д.В. Уравнение дальности при нелинейной радиолокации. - Научный вестник МГТУ ГА. Сер. Радиофизика и радиотехника, 2004, №79.
11. Козлов А.И., Колядов Д.В. Эффективная площадь рассеяния нелинейных отражателей. - Научный вестник МГТУ ГА. Сер. Радиофизика и радиотехника, 2004, №79.
12. Козлов А.И., Логвин А.И., Качалкин М.В., Колядов Д.В. Оконтуривание радиолокационных целей и возможность оценки их геометрических размеров путем управления поляризационными характеристиками излучаемых электромагнитных волн (случай когерентного рассеяния). - Научный вестник МГТУ ГА. Сер. Радиофизика и радиотехника, 2005, №96.
13. Козлов А.И., Логвин А.И., Качалкин М.В., Колядов Д.В. Оконтуривание радиолокационных целей и возможность оценки их геометрических размеров путем управления поляризационными характеристиками излучаемых электромагнитных волн (случай некогерентного рассеяния). - Научный вестник МГТУ ГА. Сер. Радиофизика и радиотехника, 2005, №93.
14. Козлов А.И., Логвин А.И., Колядов Д.В., Качалкин М.В. Оконтуривание радиолокационных целей и возможность оценки их геометрических размеров путем управления поляризационными характеристиками излучаемых электромагнитных волн (случай смешанного рассеяния). - Научный Вестник МГТУ ГА. Сер. Радиофизика и радиотехника, 2005, №96.
15. Козлов А.И., Логвин А.И., Колядов Д.В., Качалкин М.В. Оконтуривание радиолокационных целей и возможность оценки их геометрических размеров путем управления поляризационными характеристиками излучаемых электромагнитных волн (случай случайного рассеяния). - Научный Вестник МГТУ ГА. Сер. Радиофизика и радиотехника, 2005, №93.
Книги и монографии 1. Колядов Д.В. Графические представления поляризации радиоволн. В кн. А.И.
Козлов, А.И. Логвин, В.А. Сарычев УПоляризация радиоволн. Поляризационная структура радиолокационных сигналовФ - М.: Радиотехника, 2005.
2. Колядов Д.В. Матрица рассеяния нелинейных безынерционных рассеивателей.
В кн. А.И. Козлов, А.И. Логвин, В.А. Сарычев УПоляризация радиоволн. Радиолокационная поляриметрияФ - М.: Радиотехника, 2007.
3. Колядов Д.В. Применение поляризационного анализа для улучшения различения радиолокационных целей методами поляризационной адаптации. В кн. А.И. Козлов, А.И. Логвин, В.А. Сарычев УПоляризация радиоволн. Радиополяриметрия сложных по структуре сигналовФ - М.: Радиотехника, 2008.
Статьи и работы в иных отечественных журналах и изданиях 1. Козлов А.И., Колядов Д.В. УО возможности определения электрофизических и геометрических характеристик объектов радиолокационного мониторинга при их дистанционном зондировании методами радиополяриметрииФ. Труды 2-й Международной конференции УАкустооптические и радиолокационные методы измерений и обработки информацииФ, Суздаль, 2007, с. 143-146.
2. Колядов Д.В. О возможности определения электрофизических и геометрических характеристик радиолокационных объектов методами радиполяриметрии. Труды XX научной конференции МФТИ УСовременные проблемы фундаментальных и прикладных наукФ, Москва, 2007, с.98-100.
3. Козлов А.И., Логвин А.И., Колядов Д.В. Поляризационные эффекты при нелинейной радиолокации. Доклады Сибирского поляризационного семинара УСибПол-2004Ф.
Ц Томск - Сургут, 2004.
Статьи и работы в иных иностранных журналах и изданиях 1. Kozlov A.I., Ligthart L.P., Logvin A.I. Modeling and Verification of Earth-Based Radar Objects. Vol.2, УElectrodynamic characteristics derived from physical characteristicsФ, Moscow, MSTUCA, Delft, TUD, 1998 (Kolyadov D.V. - pp.13-17).
2. Kozlov A.I., Ligthart L.P., Logvin A.I. Modeling and Verification of Earth-Based Radar Objects. Vol.2, УElectrodynamic characteristics derived from physical characteristicsФ, Moscow, MSTUCA, Delft, TUD, 1998 (Kolyadov D.V. - pp.33-35).
3. Kozlov A.I., Ligthart L.P., Logvin A.I. Modeling and Verification of Earth-Based Radar Objects. Vol.4, УDeterministic and Stochastic modeling of objectsФ, Moscow, MSTUCA, Delft, TUD, 1998 (Kolyadov D.V. - pp.7-14).
4. Kozlov A.I., Ligthart L.P., Logvin A.I. Modeling and Verification of Earth-Based Radar Objects. Vol.9, УMethod to increase the radar contrastФ, Moscow, MSTUCA, Delft, TUD, 1998 (Kolyadov D.V. - pp.15-18).
5. Kozlov A.I., Ligthart L.P., Logvin A.I. Modeling and Verification of Earth-Based Radar Objects. Vol.12, УCriteria for testing the radar functionsФ, Moscow, MSTUCA, Delft, TUD, 1999 (Kolyadov D.V. - pp.23-25).
6. Kozlov A.I., Ligthart L.P., Logvin A.I. Modeling and Verification of Earth-Based Radar Objects. Vol.12, УMethods of parameter evaluationФ, Moscow, MSTUCA, Delft, TUD, 1999 (Kolyadov D.V. - pp.16-19).
7. Kozlov A.I., Ligthart L.P., Logvin A.I. Modeling and Verification of Earth-Based Radar Objects. Vol.14, УMeasurement campaigns using an 1,8 cm and 3,2 cm coherent radar with controlled polarization capabilitiesФ, Moscow, MSTUCA, Delft, TUD, 2000 (Kolyadov D.V. - pp.8-10).
8. Kozlov A.I., Ligthart L.P., Logvin A.I. Modeling and Verification of Earth-Based Radar Objects. Vol.15, УData processing and data analysis of experimentsФ, Moscow, MSTUCA, Delft, TUD, 2000 (Kolyadov D.V. - pp.21-24).
9. Kozlov A.I., Ligthart L.P., Logvin A.I. Modeling and Verification of Earth-Based Radar Objects. Vol.16, УComparisons between theory and experimentФ, Moscow, MSTUCA, Delft, TUD, 2000 (Kolyadov D.V. - pp.21-23).
10. Kozlov A.I., Ligthart L.P., Logvin A.I. Modeling and Verification of Earth-Based Radar Objects. Vol.17, УRefinement of theory and experimentФ, Moscow, MSTUCA, Delft, TUD, 2001 (Kolyadov D.V. - pp.10-13).
11. Kozlov A.I., Ligthart L.P., Logvin A.I. Modeling and Verification of Earth-Based Radar Objects. Vol.18, УMeasurement campaigns using an 1,8 cm and 3,2 cm coherent radar with controlled polarization capabilities (continuation of 14th vol.)Ф, Moscow, MSTUCA, Delft, TUD, 2001 (Kolyadov D.V. - pp.11-13).
12. Kozlov A.I., Ligthart L.P., Logvin A.I. Modeling and Verification of Earth-Based Radar Objects. Vol.19-20, УConclusion, recommendation, new areas for applicationФ, Moscow, MSTUCA, Delft, TUD, 2001 (Kolyadov D.V. - pp.25-30).
13. Kozlov A.I., Ligthart L.P., Logvin A.I. Remote Sensing of the Earth Surface based on Radar Polarimetry. Vol.1, УRelations between the polarization characteristics of the electromagnetic waves reflected from the geophysical objects and the characteristics of these geophysical objectsФ, Moscow, MSTUCA, Delft, TUD, 2003 (Kolyadov D.V. - pp.22-25).
14. Kozlov A.I., Ligthart L.P., Logvin A.I. Remote Sensing of the Earth Surface based on Radar Polarimetry. Vol.2, УSensitivity analysis on the measurement accuracy of the various polarization parameters to distinguish geophysical objectsФ, Moscow, MSTUCA, Delft, TUD, 2003 (Kolyadov D.V. - pp.15-20).
15. Kozlov A.I., Ligthart L.P., Logvin A.I. Remote Sensing of the Earth Surface based on Radar Polarimetry. Vol.3, УDerivation of accurate estimates of the polarization parameters of reflected radio wavesФ, Moscow, MSTUCA, Delft, TUD, 2003 (Kolyadov D.V. - pp.7-10).
16. Kozlov A.I., Ligthart L.P., Logvin A.I. Remote Sensing of the Earth Surface based on Radar Polarimetry. Vol.4, УCriteria for evaluation of methods to distinguish targets and geophysical objects on the basis for various polarization parametersФ, Moscow, MSTUCA, Delft, TUD, 2003 (Kolyadov D.V. - pp.17-24).
17. Kozlov A.I., Ligthart L.P., Logvin A.I. Remote Sensing of the Earth Surface based on Radar Polarimetry. Vol.5, УDevelopment of statistical models describing the process of reflected radiowaves from various geophysical objects; verification of the models by experiments and determination of statistical characteristics the target and the reflected radiowave parametersФ, Moscow, MSTUCA, Delft, TUD, 2004 (Kolyadov D.V. - pp.25-28).
18. Kozlov A.I., Ligthart L.P., Logvin A.I. Remote Sensing of the Earth Surface based on Radar Polarimetry. Vol.6, УApplication of the KLL-sphere and evaluation of efficiency of this application for distinguishing geophysical objectsФ, Moscow, MSTUCA, Delft, TUD, 2004 (Kolyadov D.V. - pp.18-21).
19. Kozlov A.I., Ligthart L.P., Logvin A.I. Remote Sensing of the Earth Surface based on Radar Polarimetry. Vol.7, УEffects of antenna polarization characteristics of distinguishing geophysical objectsФ, Moscow, MSTUCA, Delft, TUD, 2005 (Kolyadov D.V. - pp.14-19).
20. Kozlov A.I., Ligthart L.P., Logvin A.I. Remote Sensing of the Earth Surface based on Radar Polarimetry. Vol.8, УDevelopment of adaptive algorithms using polarization parameters for optimization to distinguishing geophysical objectsФ, Moscow, MSTUCA, Delft, TUD, 2005 (Kolyadov D.V. - pp.16-25).
21. Kozlov A.I., Ligthart L.P., Logvin A.I. Modeling and Algorithm Developments of Surface and Sub-Surface Sensing based Polarimetric Scatterometry. Vol.1, УMethods to distinguish geophysical objects based on their special physical properties (non-linear reflection)Ф, Moscow, MSTUCA, Delft, TUD, 2006 (Kolyadov D.V. - pp.26-29).
22. Kozlov A.I., Ligthart L.P., Logvin A.I. Modeling and Algorithm Developments of Surface and Sub-Surface Sensing based Polarimetric Scatterometry. Vol.2, УAnalysis of polarimetric radar possibilities to improve the>
23. Kolyadov D.V., Ligthart L.P., Kozlov A.I.. Amplitude-based measurement technique in polarimetric radar remote sensing for determining the dielectric permittivity of earth media. Book of abstracts, MIKON 2006, Warsaw, p.157.
24. Kolyadov D.V., Ligthart L.P., Kozlov A.I.. Amplitude-phase method allowing the determination of the complex dielectric permittivity of underlying surfaces using polarimetric radar remote sensing. Book of abstracts, EuRAD 2006, Manchester, p.154.
Колядов Дмитрий Валерьевич автореферат диссертации Радиополяриметрия природных сред и техногенных объектов при их мониторинге Подписано в печать 29.01.09 г.
Печать офсетная Формат 60х84/16 2,08 уч.-изд. л.
2,2 усл.печ.л. Заказ № 744/ Тираж 100 экз.
Московский государственный технический университет ГА 125993 Москва, Кронштадтский бульвар, д. Редакционно-издательский отдел 125493 Москва, ул. Пулковская, д.6а й Московский государственный технический университет ГА, 20 Авторефераты по всем темам >> Авторефераты по физике