функцию. Конкретно, записывая в следующем виде
[pic]
(19)
заметим, что [pic] регулярна в верхней полуплоскости, а [pic] регулярна
в нижней полуплоскости. Функции в соотношениях (16) обладают свойствами
регулярности, о которых говорилось здесь, поэтому коэффициенты
разложений по формулам (14) не являются произвольными, а их нужно
определять таким образом, чтобы исключить полюса в каждой полуплоскости.
После выполнения необходимых преобразований коэффициенты могут быть
заданы в следующем виде:
[pic] (20)
Если допустить иные разложения, чем задаваемые формулами (14), то
сохранение описанных здесь свойств регулярности становится невозможным.
Таким образом ясно, что способ разложения по формулам (14) оказывается
достаточным для рассматриваемой задачи.
УРАВНЕНИЯ ВИННЕРА-ХОПФА
В предыдущем разделе было установлено, что используя только Фурье-
компоненты рассеянной волны (конкретно, [pic]) на граничной плоскости |
y |=b, можно представить Фурье-компоненты рассеянной волны в каждой из
областей ([pic]) таким образом, чтобы удовлетворялись граничные
условия (В1), (В2), (В3). Таким образом, если в конечном счете удастся
отыскать эти неизвестные Фурье-компоненты так, чтобы удовлетворить
граничному условию (В4), то тем самым поставленная задача будет решена.
Как следует из уравнений (3), граничное условие (В4) можно свести к
непрерывности производной [pic] при | y |=b. Если записать это
требование, используя формулы (9) и (16), то окончательный результат
после выполнения необходимых преобразований, учитывающих свойства
четности Фурье-компонент, может быть записан в следующем виде (q= c, s):
[pic] (21)
Здесь кернфункции (ядра)задаются следующими формулами:
[pic] (21)
[pic]
(22)
Уравнения (21) по существу представляют собой систему уравнений
Винера-Хопфа. Этот факт становится еще более очевидным, если применить
преобразования по приводимым ниже формулам. А именно, если умножить эти
уравнения на [pic] или [pic], выполнить преобразования с использованием
соотношений (19), то можно получить два следующих соотношения:
[pic] (24)
При этом имеют место следующие соотношения:
[pic] (25)
[pic] (26)
Тот факт, что правые части в формулах (24), а в конечном счете,
правые части в формулах (26), являются регулярными в верхней
полуплоскости (области U ) или в нижней полуплоскости (области L ),
можно установить следующим образом. Ясно, что с первого взгляда можно
заключить, что особенностью функции [pic] являются только полюса в L.
Однако, эти полюса исключаются в силу соотношений, представляемых
формулой (20), так что эта функция оказывается регулярной в нижней
полуплоскости. Аналогично, функция [pic] является регулярной в верхней
полуплоскости.
Что касается системы уравнений Винера-Хопфа, представленной
формулой (24), то ее решение можно найти, выполняя интегрирование вдоль
разрезов от точек ветвления в разложениях керн-функций на множители. При
выполнении расчетов возникают определенные затруднения, однако вывод
решений проводится автоматически. Сначала выполняем разложение керн-
функций на частное (произведение) функций, регулярных в верхней
полуплоскости и в нижней полуплоскости, и, сверх того, не имеющих нулей
(j=0, 1, q= c, s ):
[pic] (27)
Затем, делим (умножаем) правые части и левые части уравнений (24)
на эти функции разложения, исключаем их полюса и выполняем
интегрирование вдоль разрезов от точек ветвления. Если выполнить
описанные действия, то обе части уравнений (24) можно разложить таким
образом, что они окажутся регулярными соответственно в верхней
полуплоскости (область U) или в нижней полуплоскости (область L). Более
того, полученные соотношения окажутся справедливыми в общей области Д.
Следовательно, по теореме Лиувилля обе части вместе оказываются
функциями, регулярными во всей плоскости, т.е. постоянными. Однако, эти
постоянные в силу граничного условия (В5) концевой точки, оказываются
равными нулю, так что решение уравнений (24) определяется единственным
образом. Здесь в качестве граничного условия концевой точки принимается
условие [pic]-const при | x | (a, | y |(b. Окончательный результат
представляется в следующем виде:
[pic] (28) Здесь [pic] выражаются через интегралы вдоль разрезов
от точек ветвления (рис. 3, 4). А именно, если считать, что функция
регулярна в точках пути интегрирования [pic] то тогда [pic] определяются
по формулам:
[pic] (29)
Наконец, [pic] обозначают нули керн-функций, а Res - вычет.
[pic]
Рис.3. Путь интегрирования [pic]
[pic]
Рис.4. Путь интегрирования [pic]
Далее
[pic] (30)
N представляет собой число нулей в нижней полуплоскости. Эти нули
соответствуют собственным значениям плоской волны вдоль пластины из
диэлектрика с потерями толщины 2b, если интерпретировать их с физической
точки зрения.
Формулы (28) представляют формальное решение на основании точного
исследования. Если в этом решении перейти к пределу [pic]то [pic] а
[pic] будет совпадать с решением для идеального проводника (подробности
доказательства опускаются).
ПРИБЛИЖЕННЫЕ РЕШЕНИЯ
Рассмотрим приближенные решения для случая, когда в среде имеются
потери, а ширина прямоугольного цилиндра велика по сравнению с длиной
волны. Если считать, что в диэлектрике имеются потери, а ширина
прямоугольного цилиндра 2а мала, то можем положить:
[pic]
(31)
Если выполнять расчеты для бетона , используя полученные к настоящему
времени данные на основании экспериментальных результатов, то даже в
диапазоне ультравысоких частот при ширине в одну длину волны абсолютные
величины в (31) ниже [pic]. Отсюда следует, что можно пренебречь
экспоненциальными членами в соотношениях (20) и членами, связанными с
плоской волной, в соотношениях (28), если ширина 2а больше длины волны.
Это означает, что почти можно пренебречь взаимодействием рассеянных волн
на обеих торцевых плоскостях (x= a и x= -a), обусловленных
электромагнитной волной, распространяющейся в среде с потерями, и
плоской волной, возбужденной в пластине из диэлектрика с потерями,
поскольку затухание вследствие потерь велико. Кроме того, функциональные
величины [pic] необходимо определять алгебраически посредством
подстановки [pic] в формулы (28), но и в этом случае можно пренебречь
всеми экспоненциальными членами. Следовательно, выполнение расчетов в
среде с потерями, в качестве которой предполагается бетон, значительно
упрощается.
Рассмотрим далее интегралы вдоль разрезов от точек ветвления, а
именно, оценим интегралы (29). Интеграл [pic] представляет собой
величину, характерную для задачи с конечной шириной; если ширина 2а
мала, то путем выбора пути интегрирования так, чтобы он был параллелен
мнимой оси, оценку можно получить по методу наибыстрейшего спуска. С
другой стороны, интеграл [pic] является величиной характерной для
полубесконечной пластины из диэлектрика. В этом случае интеграл можно
найти аналитически, если выбрать путь интегрирования параллельным
действительной оси и предположить, что подынтегральная функция такова,
что к ее регулярной части применимо приближение ломаными линиями.
Выполняя расчеты численно с использованием описанных выше методов
приближения, можно сравнительно просто отыскать неизвестные Фурье-
компоненты. Переходя в полученных результатах к обратному преобразованию
Фурье, можем найти рассеянную волну. В частности, ограничиваясь
дальнейшей рассеянной волной и используя преобразование переменных x= r
cos (, y= r sin ( и метод перевала, можно легко выполнить необходимые
расчеты.
ПРИМЕРЫ РАСЧЕТОВ И ПРИМЕРЫ ЭКСПИРИМЕНТОВ
На рис.5,6 приведены примеры, в которых найдены дальние рассеянные
волны при помощи использования приближенных методов, описанных в
предыдущем разделе. На приведенных рисунках значения представлены
нормированными на максимальное значение в дБ. Результаты, изображенные
сплошными линиями относятся к диэлектрику с потерями, а изображенные
пунктирными линиями - к идеальному проводнику. Значение[pic]
представляет один из примеров измеренного значения комплексной
относительной диэлектрической проницаемости бетона при частоте 100 МГц.
Размеры рассеивающих тел, результаты для которых приведены на рис.5,
порядка длины волны ( 2( х 2(), а на рис.6 - сравнительно большие по
отношению к длине волны (приблизительно 13( х 12(), однако, общим для
обоих рисунков является то, что на теневой стороне нет различий между
средой с большими потерями и идеальным проводником. Ясно, что отраженная
волна подавляется средой с потерями.
[pic]
[pic]
Рис.5. Дальняя рассеянная волна [pic] (дБ):
1 диэлектрик с большими потерями; 2
идеальный проводник
[pic]
[pic]
Рис.6. Дальняя рассеянная волна [pic] (дБ)
1 диэлектрик с большими потерями; 2
идеальный проводник.
Для подтверждения результатов численных расчетов были выполнены
модельные эксперименты с использованием микроволн частоты 15 ГГц.
Упрощенное изображение экспериментального оборудования, установленного в
помещении, изолированном от влияния электрических волн, показано на
рис.7. Излучающая пирамидоидальная рупорная антенна, центр рассеивающего
тела (конечной длины) и приемная пирамидоидальная рупорная антенна
устанавливались в одной и той же горизонтальной плоскости. Излучающая
антенна находилась в фиксированном положении, приемная антенна могла
передвигаться по рельсам, проложенным по окружности с центром на
центральной оси рассеивающего тела. При этом расстояние между
центральной осью рассеивающего тела и излучающей антенной составляет 3
м. (150 ().
[pic]
Рис.7. Оборудование для измерения рассеянной волны:
1 - передатчик (излучающая антенна); 2 - приемник (приемная
антенна)
То обстоятельство, что при выполнении сравнения экспериментальных
результатов и результатов численных расчетов следует соблюдать
определенную осторожность, связано с тем, что в экспериментах происходит
интерференция падающей и рассеянной волн позади рассеивающего тела и
измерение только рассеянной волны с технической точки зрения сопряжено
со значительными трудностями. Одним из методов является выполнение
сравнения только с отраженной волной. Однако, в данной работе в качестве
одной из попыток решено провести сравнение с волной, которая получается
в результате умножения падающей плоской волны на весовую функцию:
[pic]
(33)
При этом функция f(() представляет собой функцию, которая зависит от
приемопередающих характеристик измерительной системы, а именно от угла,
под которым происходит прием в электрическом поле принимаемых сигналов
приемной антенной при отсутствии рассеивающего тела. В данной работе
используется аппроксимация этой функции тригонометрическими функциями
так, чтобы при отклонении от точки на одной прямой с передающим рупором
больше, чем ( 30( происходило ослабление на - 20 дБ.
На рис.8 представлены измеренные и рассчитанные значения для
рассеянной волны в том случае, когда параллелепипед из бетона с
поперечным сечением в виде квадрата (25,5 см. х 25,5 см.) облучается
электрической волной при угле падения 60(. При этом максимум в
рассчитанных значениях равен 0 дБ, а измеренные значения представляют
собой значения, которые сопоставляются рассчитанным значениям через
максимальный уровень отраженной волны. Использованное в расчетах
значение комплексной относительной диэлектрической проницаемости [pic]
представляет собой значение, найденное по методу кратковременного
открытия (short(open) с заполнением микроволнового волновода на 15 ГГц
бетоном. Это значение, будучи сравнено с результатами измерений,
выполненных другими исследователями, представляется правильным.
Из рис.8 видно, что измеренные значения и рассчитанные значения
для отраженной волны хорошо согласуются. С другой стороны, в теневой
области (30(-90() обнаруживаются и чрезвычайно отчетливые различия в
этих значениях. В качестве первой причины этих различий можно указать на
то обстоятельство, что падающая волна не является плоской волной в
реальных экспериментах, а близка к сферической волне. Заключение об этом
можно с делать так же и по тому факту, что, вследствие распространения
фронта волны, в измеренных значениях более всего проявляется теневая
темная часть. В качестве второй
[pic]
[pic]
Рис.8. рассеянная волна на диэлектрике с большими потерями: а) -
результаты измерений (бетон); b) - результаты расчетов (при
[pic]); 1 - ((град); 2 - (дБ).
В качестве второй причины можно считать то обстоятельство, что
описанная экспериментальная система не является вполне двухмерной
моделью. В экспериментах в качестве параллелепипеда используется тело
конечной длины (1 м.), установленное на подставке; это приводит к тому,
что нельзя пренебрегать влиянием волны, отраженной от подставки. Эти
влияния проявляются в заметной интерференции измеренных значений при 10(
- 110(. В качестве других причин можно отметить, что рассеянная волна в
формуле (33) представляет собой величину, полученную применением просто
метода перевала, вряд ли являющуюся хорошим приближением.
Для изучения различий между бетоном и проводником на рис.9
приводится пример результатов для случая , когда рассеивающее тело
заменено на проводник с теми же параметрами, что и на рис.8. Измеренные
значения относятся к случаю алюминиевой пластины толщиной 1 мм.,
изготовленной для параллелепипеда, а рассчитанные значения относятся к
случаю идеального проводника. И в этом примере обнаруживается, что
измеренные значения для отраженной волны и рассчитанные значения хорошо
согласуются. Кроме того очевидно, что предположение о том, идеальным
проводником является даже алюминиевая пластина, оказывается достаточно
правильным. Наконец, сравнивая рис.8, 9, можем заключить, что подавление
отраженной волны наблюдается в среде с большими потерями. Это заключение
совпадает с заключением, сделанном на основании рис.5, 6.
На рис.10 представлены результаты, которые относятся к случаю,
когда ширина равна ширине рассеивающего тела, приведенной на рис.8, а
толщина в два раза меньше; угол падения выбран равным 45(. Очевидно, что
в той мере, насколько мала толщина, отраженная волна в направлении (=
135( слабее волны, отраженной под углом (= -45(. И в этом примере
измеренные значения для отраженной волны хорошо согласуются с
рассчитанными значениями.
[pic]
[pic]
Рис.9. Рассеянная на проводнике волна: а) - результаты измерений
(алюминий);
b) - результаты расчетов (идеальный проводник); 1 - ( (град.); 2 -
(дБ).
[pic]
[pic]
Рис.10. Рассеянная волна диэлектриком с большими потерями: а) -
результаты измерений (бетон); b) - результаты расчетов (при [pic]); 1 -
( (град.); 2 - (дБ)
7. ЗАКЛЮЧЕНИЕ
Приведена точная формулировка задачи рассеяния плоской
электромагнитной волны (Е-волны) на параллелепипеде из диэлектрика с
большими потерями, в которой используется преобразование Фурье. В
терминах преобразований Фурье приведено решение задачи. В том случае,
когда ширина рассеивающего тела сравнительно велика по отношению к длине
волны, а в среде этого тела имеются большие потери порядка потерь в
бетоне, как показывает исследование, расчеты можно значительно
упростить. Обсуждены различия по сравнению со случаем идеального
проводника.
-----------------------
P
r
y
-a
b
(
(
0
a
x
-b
[pic]
[pic]
