CC BY
40
Техніка та пристроїНВЧдіапазону. Антенна техніка
УДК 621.372
КОЭФФИЦИЕНТ ПЕРЕДАЧИ СИСТЕМЫ ЦИЛИНДРИЧЕСКИХ ДИЭЛЕКТРИЧЕСКИХ РЕЗОНАТОРОВ, СВЯЗАННЫХ С НЕСИММЕТРИЧНОЙ ПОЛОСКОВОЙ ЛИНИЕЙ
Трубаров И. В., аспирант,
Национальный технический университет Украины "Киевский политехнический институт ", г. Киев, Украина
Введение. Постановка задачи
При проектировании антенных решеток, в которых излучающими элементами являются диэлектрические резонаторы (ДР), возникает необходимость расчета АЧХ антенны. Данная задача сводится к расчету частотной зависимости коэффициента передачи системы ДР, связанных с питающей их линией. Целью настоящей работы является расчет зависимости коэффициента передачи системы ДР, связанных с несимметричной полосковой линией (НПЛ). При этом ДР размещены в ортогональной ориентации относительно НПЛ (рис. 1). Ось симметрии ДР лежит в плоскости металлической подложки линии. В самой линии, как видно из рис. 1, удалена часть заземленной металлической поверхности и диэлектрика. Будем далее называть подобную линию краевой НПЛ (КНПЛ).
Рассматривается система из множества цилиндрических ДР, координаты центров которых х и у в системе координат, изображенной на рис. 1 (начало координат по оси х соответствует середине полоска), одинаковы для
всех ДР. Резонаторы, таким образом, разнесены лишь вдоль оси z и находятся на равном удалении от линии.
Методика аналитического расчета коэффициента передачи
Для расчета коэффициента передачи системы ДР воспользуемся подходом, изложенным в [1]. Поле системы ДР можно однозначно определить, вычислив коэффициенты связи ДР между собой, открытым пространством
Рис. 1. Система из двух ДР, связанных с КНПЛ
100
Вісник Національного технічного університету України "КПІ" Серія — Радіотехніка. Радіоапаратобудування.-2011.-№45
Техніка та пристроїНВЧдіапазону. Антенна техніка
и линией передачи. Пусть последние образуют матрицу коэффициентов связи К размерности NxN, где N - количество ДР в исследуемой системе. Пусть X - вектор собственных значений матрицы К, а B=(bij) - матрица, столбцы которой являются собственными векторами К, размещенными в порядке соответствия собственным числам в X. Тогда введем N матриц Cn = (cf), n = 1, N, определенных следующим образом:
Jj — , ------
С — Л с=1
n S=1
N
£ ЬшКг, если j = n;
с=1
bj, если j * n,
где ksi - коэффициент связи s-того и i-того ДР по распространяющимся
волнам линии.
Тогда выражение для коэффициента передачи системы ДР Tf) как отношения выходной амплитуды напряженности электрического поля ко входной принимает вид:
г (f)=і - SL £ id,
|B | n=1 Qn (f)
Qnn (f) = QD Im К + f + 2iQD f0
f 1
-2--1 — Re К
fo 2 n
(1)
где Q - добротность материала диэлектрика, из которого изготовлен ДР (предполагается, что все ДР идентичны); f - частота; f0 - собственная частота основной моды ДР.
Элементы матрицы коэффициентов связи К=(кj определяются следующим образом:
kj = <
Г • ikt + i(1 - rj)kOs, если i = j;
Щ
(ikft irAz‘J + kj ) + (1 -ц)кц , если i * j
(2)
где n=0.5 - коэффициент (щє[0;1]), характеризующий часть переизлучае-мой резонатором мощности, отдаваемой в линию (в данном случае 50% запасаемой ДР энергии излучается в линию, 50% - в открытое пространство); к - коэффициент связи ДР с КНПЛ, вычисляемый в предположении, что вся запасаемая мощность отдается в линию; Г - постоянная распространения КНПЛ; AziJ=\zi-Zj\ - расстояние между центрами i-того и j-го ДР вдоль оси z; k0os - коэффициент связи ДР с открытым пространством [1]; kj - коэффициент связи между ДР по не распространяющимся волнам линии; kjos - коэффициент взаимной связи между ДР в свободном пространстве [2].
Значение коэффициента kl рассчитывается из выражения, которое может быть записано в виде:
kl = Af 2(0)Hx(% zoX
где А - величина, зависящая от параметров резонатора и поля основной
Вісник Національного технічного університету України "КПІ" Серія — Радіотехніка. Радіоапаратобудування.-2011.-№45
101
Техніка та пристроїНВЧдіапазону. Антенна техніка
моды его собственных колебаний; f(^) - функция связи, определенная в [3]; Hx(x0, у0, z0) - соответствующая компонента (в системе координат, указанной на рис. 1) напряженности магнитного поля НПЛ в точке размещения центра цилиндрического ДР.
В выражении (2) диагональные элементы ku характеризуют связь l-того ДР с линией и свободным пространством. Недиагональные элементы описывают взаимную связь между резонаторами. Можно также видеть, что матрица К является взвешенной суммой (с весами п и 1-ц) двух матриц коэффициентов связи, одна из которых рассчитывается в предположении, что вся запасаемая ДР энергия отдается в линию (берется с коэффициентом п), а другая - в предположении, что ДР размещены в открытом пространстве в отсутствие линии передачи (берется с коэффициентом 1-п).
При объединении нескольких идентичных ДР в систему в ее АЧХ наблюдается расщепление частот относительно резонансной частоты одного ДР f0 (см. ниже). При сравнении экспериментально измеренных и рассчитанных теоретически кривых коэффициента передачи системы нескольких ДР было замечено, что расщепление частот, предсказываемое моделью, меньше имеющего место в действительности. Данное несоответствие исчезает при введении в выражение для недиагональных элементов в (2) параметра kj (действительная величина), который может быть истолкован как коэффициент связи двух ДР по не распространяющимся волнам КНПЛ, поскольку в регулярной линии резонаторы связываются между собой как по распространяющимся, так и по не распространяющимся волнам [1], а слагаемое ikle~lT^lJ соответствует связи только по распространяющимся волнам.
В отсутствие аналитического выражения для kJ для количественной оценки этой величины можно использовать аппроксимирующую кривую kiJ0(Azij), построенную по нескольким экспериментально измеренным точкам. По построенной кривой затем можно оценивать данную величину для расчета системы любого количества ДР, размещенных на любом расстоянии друг от друга вдоль оси z.
Нахождение аппроксимирующей кривой для зависимости kij0(Azij) по экспериментально измеренным значениям
Для экспериментального измерения коэффициента kj использовалась система из двух резонаторов (ДР1 и ДР2). Таким образом, для нескольких значений расстояния между центрами ДР1 и ДР2 Az12 необходимо вычислить значения k12 , и затем аппроксимировать полученные точки подходящей кривой.
АЧХ системы из двух ДР будет иметь два локальных минимума, расположенных на частотах f и f2. Наиболее простой способ определения значения k12 состоит в измерении разности частот локальных минимумов АЧХ системы Af= f1 - f2. Тогда можно в явном виде записать выражение
102
Вісник Національного технічного університету України "КПІ" Серія — Радіотехніка. Радіоапаратобудування.-2011.-№45
Техніка та пристроїНВЧдіапазону. Антенна техніка
для действительной части собственных чисел матрицы К, зависящее от искомого параметра k120 , причем Re X(kj02) = - Re X(k12). Тогда из соотношения
f - f
[1] Re X = 2——— непосредственно следует f0
|Re \2(k°)|= = f
f0 f0
Подставляя в последнее уравнение явное выражение для Re X и решая
его относительно k12 , получаем:
k = — Л12
п
1 (а/
----nkl sinГАг-(1 -пПRek
f0
12
(3)
Таким образом, для нескольких значений расстояния между ДР Az12 по АЧХ измеряется А/, а затем из (3) находится соответствующее этому расстоянию значение k12 .
Однако данный способ имеет существенный недостаток: при удалении ДР друг от друга один из локальных минимумов АЧХ системы становится все менее выраженным, и при превышении некоторого значения расстояния между ДР А/не может быть измерено из АЧХ (см. ниже). Таким образом, k120 может быть измерен только на незначительных удалениях ДР друг от друга, что недостаточно.
Для измерения k120 для любых расстояний между ДР применялся второй, более сложный, подход, состоящий в использовании метода наименьших квадратов (МНК). Ниже изложим суть данного подхода и полученные с его помощью результаты.
Пусть в выражении (2) величина k120=k210=k0 - дополнительная переменная. Тогда определяемая из (1) функция коэффициента передачи является функцией двух переменных: частоты/и коэффициента k0: T(f,k0). Рассмотрим измерение величины k0 для одного значения Az12. Для данного расстояния между двумя ДР измеряется АЧХ системы в M точках. Тогда имеем набор из M пар (fm;Lm), m = 1,M, где/m - частота в m-той точке измерения, Lm - уровень затухания сигнала на данной частоте, дБ. Перейдем от
Lm
логарифмических единиц к натуральным: Tmx = 10 20 , Tmxp є[0; 1]. Тогда имеем M пар (fm; Tmexp). Дальнейшая задача состоит в отыскании такого значения k0=k0min, при котором функция T(f,k°min) минимально отличалась бы от экспериментально измеренных значений. Поскольку очевидно, что k0 входит в Tf,k0) нелинейно и к линейной зависимости нельзя также прийти логарифмированием (1), то реализация МНК в виде решения системы линейных уравнений оказывается невозможной. В таком случае решение будем искать путем минимизации функции суммы квадратов разностей теоретических и экспериментальных значений коэффициента передачи в точках измерения. Данная функция имеет вид:
Вісник Національного технічного університету України "КПІ" Серія — Радіотехніка. Радіоапаратобудування.-2011.-№45
103
Техніка та пристроїНВЧдіапазону. Антенна техніка
M 2
s(к”) = £(гUm, k0) - г;г) .
m=1
На рис. 2 приведены примеры нахождения величины k0min для двух случаев, когда расстояние между ДР невелико и на АЧХ видны локальные минимумы (рис. 2, а), и когда расстояние значительно и четко различим лишь один минимум (рис. 2, в), а также теоретические и экспериментальные кривые АЧХ при использовании в (2) определенных таким образом значений к0 (рис. 2, б, г). На рис. 2 величина d=Az12-2r0 (где r0 - радиус цилиндрического ДР) - расстояние между ближайшими точками резонаторов (легко измеримая величина).
1.5
0.5
& т
0.02
0.04
-5
- 10
- 15
/. -
'А уОЧг- .-■у
і L ■ 1 /
Theory -3dB Experiment 1
&)d=1.3 мм; k0min=0.016
4xl0S 4.1х109 4.2x10s 4.3x10s £ Hz
б) АЧХ системы при d=1.3 мм
1.5
0.5
S(k)
0.02
0.04
в) d=4.2 мм; k0min=0.004
4.3x10s £ Hz
г) АЧХ системы при d=4.2 мм
Рис. 2. Определение величины k0min при нахождении минимума функции S(k ). а, в - кривые функции S(k0) для двух различных значений расстояния между ДР d; б, г - соответствующие им АЧХ исследуемой системы двух резонаторов
Таким образом для каждого из нескольких значений расстояния между ДР d находится величина k0min. Оказывается, что при увеличении d величина k уменьшается экспоненциально, т. е. полученная зависимость k0min(d)аппроксимируется экспоненциальной функцией. На рис. 3 изображены измеренные значения k min и проведена аппроксимация их экспонен-
104
Вісник Національного технічного університету України "КПІ" Серія — Радіотехніка. Радіоапаратобудування.-2011.-№45
Техніка та пристроїНВЧдіапазону. Антенна техніка
циальной функцией.
\ 1 і ♦♦♦ Experiment
\ \ \ 2 Approximation
\ ♦х ч Ч
^ ___
0 2 4 fid, mtn S
Рис. 3. Аппроксимация измеренных значений k0min экспоненциальной функцией
Сравнение теоретических и экспериментальных результатов
Для проверки теоретической модели измерялась АЧХ системы из 3-х цилиндрических ДР (ДР1, ДР2, ДР3), размещенных друг от друга на расстояниях Azi2= Az23=0.5Az13=11.8 мм (что соответствует d=1.7 мм). Параметры ДР: высота L=2.5мм, диаметр 2г0=10.1мм, диэлектрическая проницаемость материала ДР є=81. Параметры КНПЛ: : высота диэлектрической подложки h=1.2мм, шириной полоска w=2мм, диэлектрической проницаемостью подложки є=2 (фторопласт). ДР были размещены вплотную к КНПЛ, так что x1=x2=x3=w/2 (система координат на рис. 1). Для линии и ДР с этими же параметрами проводились расчеты, описываемые выше.
На рис. 3 аппроксимирующая функция имеет вид k0(d) = 0.03е-045d (d в мм). Тогда при расчете системы из 3-х резонаторов в (2) подставлялись значения kj = k 0(Azy. - 2r0).
На рис. 4 изображены теоретическая и экспериментальная АЧХ данной системы.
Рис. 4. Сравнение теоретической и экспериментальной кривых частотной зависимости коэффициента передачи системы 3-х ДР, связанных с КНПЛ
Вісник Національного технічного університету України "КПІ" Серія — Радіотехніка. Радіоапаратобудування.-2011.-№45
105
Техніка та пристроїНВЧдіапазону. Антенна техніка
Выводы
В работе представлена модель расчета АЧХ системы N цилиндрических ДР, возбуждаемых КНПЛ. Основной трудностью при расчете является отсутствие аналитического выражения для вычисления коэффициента связи ДР по не распространяющимся волнам линии. Для того, чтобы иметь возможность рассчитывать АЧХ системы любого количества ДР, необходимо сделать несколько измерений АЧХ системы из двух ДР. При этом было выяснено, что искомая зависимость имеет экспоненциальный характер, поэтому достаточно произвести лишь два измерения (при двух различных расстояниях между ДР), чего будет достаточно для определения неизвестных коэффициентов в уравнении экспоненциальной функции, аппроксимирующей искомую зависимость коэффициента от расстояния между ДР.
Литература
1. М. Е. Ильченко, А. А. Трубин. Электродинамика диэлектрических резонаторов. -Киев: Наукова думка - 2004.
2. Трубін О. О., Шмиглюк Г. С. Моделювання параметрів антенної решітки на циліндричних діелектричних резонаторах// Вісник Національного технічного університету України "КПІ". Серія - Радіотехніка. Радіоапаратуробудування. - 2006. - №33. - С. 101-108.
3. Трубин А. А. Расчет связи цилиндрического диэлектрического резонатора со щелевой линией передачи // Радиоэлектроника, №4, 1996, с. 61-66.
Трубаров І. В.Коефіцієнт передачі системи циліндричних діелектричних резонаторів, зв’язаних з несиметричною смужковою лінією. Викладено методику аналітичного розрахунку коефіцієнту передачі системи циліндричних діелектричних резонаторів у випадку їх ортогональної орієнтації відносно крайової несиметричної смуж-ковоїлінії. Наведено два способи експериментального визначення величини коефіцієнту зв ’язку резонаторів по хвилям, що не розповсюджуються, смужкової лінії передачі: шляхом виміру величини розщепленнячастот в А ЧХ системи та шляхом використання метода найменших квадратів. Показано, що залежність коефіцієнту зв ’язку резонаторів по хвилях, що не розповсюджуються, від відстані між ними може бути апрок-симована експоненціальною функцією. Проведено порівняння експериментальних і теоретичних даних.
Ключові слова: діелектричний резонатор, коефіцієнт передачі, НВЧ, АЧХ, коефіцієнт зв’язку, несиметрична смужкова лінія, антена, зв’язані коливання, розщеплення частот.
Трубаров И. В. Коэффициент передачи системы цилиндрических диэлектрических резонаторов, связанных с несимметричной полосковой линией. Изложена методика аналитического расчета коэффициента передачи системы цилиндрических диэлектрических резонаторов в случае их ортогональной ориентации относительно краевой несимметричной полосковой линии. Приведены два способа экспериментального определения величины коэффициента связи резонаторов по не распространяющимся волнам полосковой линии передачи: путем измерения величины расщепления частот в А ЧХ системы и путем использования метода наименьших квадратов. Показано, что
106
Вісник Національного технічного університету України "КПІ" Серія — Радіотехніка. Радіоапаратобудування.-2011.-№45
Техніка та пристроїНВЧдіапазону. Антенна техніка
зависимость коэффициента связи резонаторов по не распространяющимся волнам от расстояния между ними может быть аппроксимирована экспоненциальной функцией. Проведено сравнение экспериментальных и теоретических данных.
Ключевые слова: диэлектрических резонатор, коэффициент передачи, СВЧ, АЧХ, коэффициент связи, несимметричная полосковая линия, антенна, связанные колебания, расщепление частот.
Trubarov I. V. Transmission gain of a system of cylindrical dielectric resonators coupled with a microstrip line. The method of analytical calculation of gain of a system of cylindrical dielectric resonators in case of their orthogonal orientation relative to an edge microstrip line is stated. Two ways of experimental determination of value of coupling coefficient between the resonators by evanescent waves of microstrip line are stated: measuring frequency splitting value and by force of least squares method. It is disclosed, that the dependence of coupling coefficient between resonators by evanescent waves versus a distance between them can be approximated by exponential function. The comparison between experimental and theoretical data is provided.
Keywords: dielectric resonator, transmission gain, microwave frequencies, amplitude-frequency characteristic, coupling coefficient, microstrip line, antenna, coupled oscillation, frequency splitting.
Вісник Національного технічного університету України "КПІ" Серія — Радіотехніка. Радіоапаратобудування.-2011.-№45
107
