Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:
Кураев А.А., Попкова Т.Л., Синицын А.К. Электродинамика и распространение радиоволн.pdf
Скачиваний:
262
Добавлен:
02.05.2014
Размер:
4.84 Mб
Скачать

даемого негармоническими электрическими и магнитными токами источников.

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

сиального волновода не имеют изломов, т.е. bz и ϕb не имеют разрывов.

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

8.7. Уравнения возбуждения нерегулярных замедляющих систем

Поставим задачу следующим образом. Требуется найти решение уравнений Максвелла для гармонических процессов

 

r&

 

 

r&

 

r&

 

rotH = jωε

 

E

+δ

,

 

 

 

&

 

 

 

 

(8.65)

r

 

a

 

 

r

 

 

rotE&

= − jωµ&

a

H&

 

 

 

 

 

 

 

 

 

 

при граничных условиях импедансного типа

E&

= Z&(z).

(8.66)

&

τ

 

 

Hτ

Здесь координата z соответствует направлению оси замедляющей системы, поперечные координаты (в общем случае криволинейные) обозначим

q1 , q2.

Будем считать, что контур поперечного сечения замедляющей системы не зависит от z и её нерегулярность обусловлена только зависимостью от z импеданса стенок (или эквивалентных им боковых поверхностей, на которых задан Z&(z)).

В качестве

базисных функций

изберем «квазирегулярную» систему

r

r

 

 

 

функций ES , H S вида

 

 

r

r

(q1, q2 ,

S (z))ej hS dz ,

 

E&S = ES0

 

r

r

(q1,

q2 ,

S (z))ej hS dz ,

 

H&S = HS0

(8.67)

hS =

k 2 + 2S .

 

 

106

r r

Функции ES , H S удовлетворяют граничным условиям (8.66) в каждом сечении z' и являются решениями однородных уравнений (8.65) для регулярной системы с Z& = Z&(z) при всех z (соответственно и S (z)= Const в экви-

валентной регулярной системе). Система функций ErS ортогональна в ка-

H S

ждом сечении z, как и всякая система собственных волн регулярного волновода, т.е.

 

r

r

 

 

r

r

 

r

0,

 

p S

.

(8.68)

 

 

 

 

 

 

 

 

JS ,P = { ES ,

HP

EP ,

HS }z0dS =

 

 

S

 

 

 

 

 

 

 

 

NS

,

p = S

 

 

Поскольку, однако, S = S (z)Const , поля (8.67) не удовлетворяют однородным уравнениям Максвелла и система уравнений для них имеет вид

rotHrS = jωε&a ErS +ψrSe , rotErS = − jωµ&a HrS ψrSm ,

где

r

 

 

j

h dZ

 

 

 

r

 

 

r

 

 

ψ e

= −e

 

 

S

 

 

 

H

0

, z

 

 

,

 

 

 

 

 

 

 

S

 

 

 

 

 

 

z

r

S

 

 

0

 

r

 

j

 

h dZ

 

 

 

 

r

 

 

 

ψ m = e

 

 

S

 

 

E0

, z

 

 

.

 

 

 

 

 

 

S

 

 

 

 

 

z

S

 

 

0

 

 

 

Таким образом, ψrSe , ψrSm - чисто поперечные вектора, что существенно

в последующем выводе уравнений возбуждения.

Разделим все вектора на поперечные и продольные и запишем разложения для поперечных составляющих Et , Ht в виде

Ert = (C&S (z)ErSt + C&S (z)ESt ),

S

Hrt = (C&S (z)HrSt + C&S (z)H St ).

S

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

стве используется тот факт, что ψrSml =ψrSel = 0 ):

107

 

 

 

 

 

r

 

 

 

r

 

 

 

 

r

 

 

 

 

δ

l

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

E = (C&S (Z )ES + C&

S (Z )ES

)

 

 

 

 

,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

S

 

 

 

 

 

 

 

 

 

 

 

 

 

 

jωεα

 

 

 

 

 

 

 

 

(8.69)

 

 

 

 

 

H = (C&S (z)HS + C&S (z)HS ).

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

S

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Для определения коэффициентов разложения C&±S (z) воспользуемся

леммой Лоренца для бесконечно малого объема S dz

в волноводе, предпо-

лагая, что Z&(z) и соответственно S (z) - гладкие функции. В соответствии с

леммой Лоренца для dV = S dz можно записать

 

 

 

 

 

 

 

 

d

 

 

 

r r

 

 

r r

r

 

 

 

 

 

 

 

 

r

 

r

 

r r

r

r

r

r

 

)dS .

 

 

 

{[E , H

 

][E , H

 

]}Z dS

 

=

 

 

(δ e E

 

δ e E

δ m H

+δ m H

 

(8.70)

 

dz

 

 

 

 

 

 

 

 

1

2

 

 

2

 

1

 

 

0

 

 

 

 

 

 

1 2

2 1

 

1 2

 

2 1

 

 

 

 

 

 

S

 

 

 

 

 

 

 

 

S

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Полагая в качестве

E , H

 

 

 

поля (8.69) (δ e

= δ ,

δrm

= 0 ), а в качестве

 

 

r

 

r

 

 

 

 

 

r

 

 

 

 

r 1

 

1

 

 

 

δ m2

=ψrmmS )

1

 

1

 

 

 

 

 

E2

, H 2

поля

E±S , H ±S

(δ 2e =ψrmeS ,

и с учетом условия ортого-

нальности (8.68), из (8.70) получаем

 

 

 

 

 

 

 

 

 

 

 

 

 

d

(C&S NS )=

r r

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

δ ES dS + C&P

γ p, S

 

 

 

 

 

 

 

 

 

 

 

 

 

dz

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

S

 

 

 

 

 

 

P

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(8.71)

 

 

d

 

(C&S NS )=

r r

 

 

+ C&P

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

δ ES dS

 

γ p, S .

 

 

 

 

 

 

 

 

 

 

dz

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

S

 

 

 

 

 

 

P

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Здесь

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

γ p, ±S = γ mp, ±S γ ep, ±S = γ ±S, p , γ ep, ±S =

r

 

 

 

 

 

r

r

EPψ±eS dS , γ mp, ±S = H Pψ±mS dS .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

S

 

 

 

 

S

 

Система (8.71) представляет собой совершенно общую форму уравнений возбуждения для произвольной нерегулярной замедляющей системы.

Заметим, что (8.71) нетрудно видоизменить на случай, когда выделяет-

r

ρ

 

ся квазистатическая часть электрического поля Ecm = −gradΦ, 2Φ = −

.

 

 

εa

В этом случае в (8.71) δ необходимо заменить на δ ′ = δ jωεa gradΦ.

 

 

8.8. Нерегулярные волноводы с прямоугольным сечением.

108

Теория и приложения

Нерегулярные волноводы с прямоугольным сечением весьма широко используются как в различных устройствах СВЧ и антенной технике, так и в приборах СВЧ и КВЧ с пространственно - развитыми ленточными электронными потоками: ЛБВ, ЛОВО- и М-типов, оротронах, убитронах и других типах лазеров на свободных электронах.

В разделе 8.4 изложена теория возбуждения продольно - азимутально - нерегулярных волноводов. В ее основе лежит операция отображения внутренней поверхности нерегулярного волновода на регулярный круговой цилиндр единичного радиуса. Соответственно используются цилиндрические координаты r,ϕ, z . В случае нерегулярных волноводов с прямоугольным се-

чением этой теорией не только неудобно пользоваться из-за несоответствия естественных для прямоугольного волновода декартовых координат x, y, z с принятым в теории цилиндрическими, но и, строго говоря, невозможно, поскольку в угловых точках прямоугольного сечения имеет место разрыв производных b(z, ϕ)/ ϕ (b(z, ϕ) - радиус внутреннего контура волновода в

координатах r,ϕ, z).

Ниже развита строгая теория возбуждения произвольно – нерегулярных волноводов с прямоугольным поперечным сечением, основанная на операции отображения нерегулярной внутренней поверхности волновода на регулярный цилиндр с прямоугольным сечением и использовании прямоугольной системы координат. Приведены примеры расчета ЛБВ-О, заграждающего фильтра, а также дисперсионных характеристик периодического волновода на основе предложенной теории.

Преобразование координат. Рассмотрим (Рис. 8.1) произвольно (по х и

y

z

y

x

S1

S1 S2

x

S2

Рис.8.1.

у) нерегулярный прямоугольный волновод (x, y, z – компоненты исходной прямоугольной системы координат). На (Рис.8.2, а) изображено начальное поперечное сечение волновода z=0. Система координат выбрана таким обра-

109

 

 

 

 

y

 

 

 

 

y

 

 

Dyp(z)

 

v

 

 

 

 

 

 

 

Dy0

 

 

 

 

 

α

 

 

 

 

 

 

x

 

x

 

 

 

u

 

 

 

 

 

 

 

Dx0

 

-Dx0 Dxp(z)

 

 

Dxm(z)

1

 

-1

 

-Dy0

 

 

 

 

 

-α

 

 

 

 

 

Dym(z)

 

 

 

 

 

а)

 

б)

 

 

 

в)

Рис.8.2.

зом, что ее начало (точка x=y=z=0) соответствует центру поперечного сечения. В текущем поперечном сечении z нерегулярного волновода (Рис.8.2, б) это не так: в общем случае Dxp(z)-Dxm(z), Dyp(z)-Dym(z). Задача состоит в определении поля, возбуждаемого в волноводе источниками, заданными плотностью стороннего электрического тока

δr = xr0δx (x, y, z,t)+ yr0δ y (x, y, z,t)+ zr0δz (x, y, z,t)

и плотностью стороннего магнитного тока

δrм = xr δ м (x, y, z, t )+ yr δ м (x, y, z, t )+ zr δ м (x, y, z, t ).

0 x

0 y

0 z

Искомое поле должно удовлетворять граничным условиям на боковых стенках волновода S1,2 (потерями в стенках пренебрегаем, σ → ∞ )

r

r

 

S1,2

=0

(8.72)

 

n1,2

E

 

( nr1,2 - внешняя нормаль на вертикальных и горизонтальных стенках волно-

вода).

Для решения поставленной задачи воспользуемся методом преобразования координат, позволяющим преобразовать граничную задачу (8.72) (Рис. 8.1, Рис.8.2, б) к элементарной (Рис.8.2, в): контур поперечного сечения симметричен относительно z и не зависит от z. Введем следующие преобразования координат:

u =

x + b

,

v =

y + h

, z=z,

(прямое),

(8.73)

a

q

 

 

 

 

 

 

110

x = a u b, y = q v h

(обратное),

(8.74)

где a(z) = (Dxp (z) Dxm (z)) / 2,

b = (Dxp (z) + Dxm (z)) / 2,

 

q(z) = (Dyp (z) Dym (z)) / 2α, h = (Dyp (z) + Dym (z)) / 2, α = Dyo / Dxo.

При таком преобразовании в новых переменных внутренняя граничная поверхность волновода регулярна (Рис.8.2, в):

u = ±1

(S '

)

и ν = ±α

(S '

).

 

1

 

 

2

 

Учитывая обратное преобразование (8.74), для радиуса-вектора точки во внутренней области в новой системе координат u,v, z имеем

rr(u, v, z) = zr z + xr

 

(au + b)

+ yr

(qv + h).

 

 

 

 

 

 

 

0

0

 

 

 

 

0

 

 

 

 

 

 

Обозначим далее a

z

= da

, b

 

= db ,

q

z

= dq

,

h

= dh .

 

 

 

 

 

 

 

dz

z

 

dz

 

dz

 

z

dz

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Определим основную систему векторов косоугольной координатной

системы u,v, z :

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

r

 

rr

 

r

r

 

 

rr

 

 

r

 

 

 

 

 

 

 

a1

=

u = a x0 ,

a2

=

v = q y0

,

 

 

 

 

 

 

r

 

rr

r

 

 

 

 

r

 

 

 

 

 

r

 

 

 

 

a3

=

z = z0

+ (uaz + bz )x0

+ (vqz + hz )y0 .

 

 

 

Взаимная система векторов ar1, ar2 , ar3 находится следующим образом:

ar1 = ar2 V×a3 , ar2 = ar3V×a1 , ar3 = a1 V×a2 ,

V = a1(a2 ×a3 )= a2 (a3 ×a1 )= a3 (a1 ×a2 )= aq .

Производя указанные действия, имеем

ar1 = a1 (xr0 zr0 (uaz + bz )),

ar2 = q1 (yr0 zr0 (vqz + hz )),

111

 

ar3 = zr .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Найдем теперь элементы метрического тензора

 

 

g11 = ar1 ar1 = a2 (1 + (uaz + bz )2 ),

 

 

 

 

 

 

 

 

 

g22 = ar2 ar2 = q2 (1 + (vqz + hz )2 ),

 

 

 

 

 

 

 

 

 

g33 =1,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

g12 = ar1 ar2 = (ua

z

+ b

) (vq

z

+ h ) = g21

,

 

 

 

 

 

 

 

 

 

 

 

z

 

 

 

z

 

 

 

 

 

 

 

 

 

g13 = ar1 ar3 = −a1 (ua

z

+ b

) = g31 ,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

z

 

 

 

 

 

 

 

 

 

 

 

 

g23 = ar2 ar3 = −q1 (vq

z

+ h

).

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

z

 

 

 

 

 

 

 

 

 

 

 

r

Представим вектора электромагнитного поля и тока в виде (на примере

H ):

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

r

xr + H

yr + H

zr

 

= H

 

ar1 + H

 

ar2

+ H

 

ar3 ,

 

 

H = H

 

u

v

ζ

 

 

 

x 0

 

y 0

 

 

z

0

 

 

 

 

 

 

 

 

 

где Hu , Hv , Hζ

- ковариантные проекции H , которые связаны с исходны-

ми следующим образом:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Hx

= Hu / a, H y = Hv / q,

Hz

= Hζ

azu + bz Hu

qzv + hz

Hv .

(8.75)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a

 

 

 

 

 

q

 

Запишем теперь уравнения Максвелла в новой системе координат

1 Hζ

 

H

v

r

 

 

 

 

H

u

 

 

Hζ

 

r

 

H

v

H

u

r

 

 

 

 

 

 

 

 

 

 

 

 

a

 

+

 

 

 

 

 

 

 

 

a

+

 

 

 

 

 

 

a

 

 

=

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

3

 

V

 

v

 

 

z

 

 

 

 

 

z

 

 

 

u

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

u

 

v

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

E

r1

 

E

 

 

r2

 

 

 

Eζ

 

r3

 

 

 

 

 

r1

 

r2

 

r

r3

 

 

 

 

 

 

 

=ε0

u a

 

+

 

 

v

a

 

 

+

 

 

 

a

 

+δua

+δva

 

+δζ a

 

,

 

 

 

 

 

 

 

 

 

 

 

 

t

 

 

 

 

 

 

 

 

 

 

 

t

 

 

 

 

t

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1 Eζ

 

E

 

r

 

 

 

E

 

 

 

Eζ r

 

 

E

 

 

E

r

 

 

 

 

 

 

 

 

 

v

 

a

 

+

 

 

 

u

 

 

 

 

 

a

+

 

v

 

u

 

a

 

 

 

=

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

3

 

 

 

V

 

v

 

 

z

 

 

 

 

z

 

 

 

u

 

u

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

v

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(8.76)

(8.77)

112

 

H

u

r1

 

H

v

r2

 

Hζ

r3

 

м r1

м r2

м r3

 

µ0

 

a

+

 

a

+

 

a

 

δu a

δv a

δζ a

.

t

t

t

 

 

 

 

 

 

 

 

 

 

 

 

 

Используя далее свойство ортогональности основной и взаимной системы векторов ari ari = δij и умножая скалярно уравнения (8.76), (8.77) на ar1, ar2 , ar3 , получаем контравариантные проекции уравнений в форме

 

 

 

 

 

 

r

 

 

r1

 

 

 

 

 

 

 

 

 

 

 

 

r

1

=

ε0 gˆ

E1

 

 

,

 

 

 

 

 

 

 

 

 

rotH

 

 

+ gˆδ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

tr

 

 

 

 

 

 

 

 

 

 

 

 

 

(8.78)

 

r1

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

= −µ0 gˆ

H

 

 

r

 

 

 

 

 

 

 

 

rotE

t

gˆδ

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Здесь gˆ

=Vg,

 

 

 

 

r

 

 

r

r

r

,

r

r

r

r

(анало-

g = (gij ), r

= ux0

+ vy0

+ zz0

H 1

= Hu x0

+ Hv y0

+ Hζ z0

r

 

r

 

r

1m ).

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

гично E1

,δ 1,δ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Граничные условия (1) теперь принимают вид:

 

 

 

 

E1

 

vα

= 0, E1

 

u1

=0 .

 

 

 

 

 

 

(8.79)

 

 

 

 

 

 

 

 

u, z

 

 

v, z

 

 

 

 

 

 

 

 

 

 

Проекционные соотношения для амплитуд связанных волн нерегуляр-

ного волновода. Прежде чем переходить к решению задачи (8.78), (8.79), целесообразно для улучшения его сходимости выделить в (8.78) электростати-

ческую часть поля источников, содержащую разрыв первой производной E1 ,

и магнитостатическую, содержащую разрыв H 1 . При этом динамическая задача имеет вид

r

 

 

r

r

 

 

 

 

 

 

 

=ε0 gˆ E1

,

 

 

 

 

rotH1

+ gˆδ1

 

 

 

 

 

r

 

 

 

 

rt

r

 

 

 

 

 

 

(8.80)

 

 

 

 

 

 

 

 

 

 

rotE

= −µ

0

gˆ H1

gˆδ

 

м,

 

 

 

 

1

 

 

 

t

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

E1 u, z

 

vα =0,

E1 v, z

 

u1 =0 .

 

 

(8.81)

 

 

 

 

r

r

1 ε0 grad

(Фe

t ),

r

 

r

 

Здесь δ1 =δ

E1

= E1 + grad Фe ,

 

r

 

r

µ0 grad (Фм

t ),

r

 

r

м ,

δ1м =δ

H1

= H 1 + grad Ф

Фe, Фм – соответственно электрический и магнитный потенциалы источников.

113

Существенно, что Еr1 , Н1 - непрерывные на границе источников векторы и операция почленного дифференцирования представляющих их в решении рядов (rotE1,rotHr1 ) допустима, поскольку эти ряды сходятся равномер-

но.

Остановимся на решении задачи (8.80), (8.81), полагая режим установившимся (периодическим). Представим искомое решение в системе u,v, z в

виде разложения по собственным функциям регулярного волновода следующим образом:

 

r

 

 

 

 

 

 

r

 

 

 

 

r

 

 

 

 

 

 

 

 

 

 

 

 

 

r

 

 

 

 

 

 

 

 

 

 

 

r

 

 

 

 

r

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

E1 = Re(E&ts

 

+ E&ζ s )e jsωt , H1 = Re(H&ts + H&ζ s )e jsωt ,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

s

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

s

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

r&

 

 

M N

 

 

e

 

 

 

 

 

 

re

 

 

 

 

K L

 

м

 

 

rм

,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Ets = ∑∑Asmn (z)emn +

∑∑Askl(z)ekl

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

&

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

&

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

m=1 n=1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

k=0 l=0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

r

 

 

M

N

 

 

 

 

 

 

 

 

 

 

 

 

 

zr0 ,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

E&ζ s

= ∑∑C&smn (z)ψmne

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

m=1 n=1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

r

 

 

M

N

 

 

 

 

 

 

 

 

 

r

 

 

 

 

 

 

K

 

 

 

L

 

 

 

 

 

 

r

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

H&ts

= ∑∑B&smne (z)hmne

+ ∑∑B&skмl(z)hkмl

,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

m=1 n=1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

k=0 l=0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

r

 

 

K

L

 

 

 

 

 

 

 

 

 

 

м r

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Hζ s

 

= ∑∑Dskl(z)ψklz0 .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

&

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

k=0 l=0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Здесь ψmne

 

= sin

 

mπ (u +1)

sin

nπ (v +1)

,

ψkмl

 

= cos

kπ (u +1)

cos

lπ (v +1)

,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

2α

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

2α

 

 

 

 

 

 

 

 

 

 

 

 

r

 

 

 

 

π

 

 

 

 

 

mπ

(

u +1

)

 

 

 

 

nπ

(

v +α

)

 

r

 

 

 

 

π

 

 

 

 

mπ

(

u +1

 

 

 

nπ

(

v +α

)

 

r

 

 

ee

= ψe

 

 

= m

 

cos

 

 

 

 

 

 

 

 

 

sin

 

 

 

 

 

 

 

 

x

 

+

n

sin

 

 

 

 

 

 

 

 

)

cos

 

 

 

 

 

 

 

 

 

y

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

2α

 

 

 

2α

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

2α

 

 

 

 

 

mn

mn

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

r

r

 

 

 

 

 

lπ

 

 

 

 

kπ

(

u +1

 

 

 

 

 

lπ

(

v +α

)

 

r

 

 

π

 

 

 

kπ

(

u +1

)

 

 

 

 

lπ

(

v +α

)

 

r

e м

= × z

ψe =

 

 

 

 

cos

 

 

 

 

 

 

 

 

 

)

sin

 

 

 

 

 

 

 

 

x

k

sin

 

 

 

 

 

 

 

 

 

 

 

cos

 

 

 

 

 

 

 

 

 

y

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

kl

0

 

kl

 

 

2α

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

2α

 

 

 

 

 

0

 

 

2

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

2α

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

r

r

 

 

 

 

 

π

 

 

 

 

 

mπ u +1

 

 

 

nπ v +α

 

 

r

 

 

 

 

π

 

 

 

 

mπ u +1

 

 

 

 

 

nπ v +α

 

r

 

 

 

 

 

 

 

n

 

 

 

 

 

 

 

 

 

 

 

(

 

 

 

 

 

)

 

 

 

 

 

 

 

(

 

 

 

 

 

)

 

 

m

 

 

 

 

 

 

 

 

 

 

(

 

 

 

 

 

)

 

 

 

 

 

 

 

 

 

(

 

 

 

 

)

 

he

= × z

ψ м

=

 

 

 

sin

 

 

 

 

 

 

 

 

 

 

cos

 

 

 

 

 

 

 

 

x

 

 

cos

 

 

 

 

 

 

 

 

 

 

 

sin

 

 

 

 

 

 

 

 

 

y

2α

 

 

 

 

 

 

2

 

 

 

 

 

 

 

2α

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

2α

 

 

 

mn

 

0 mn

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

hrм

= − ψ

м = kπ sin

kπ (u +1)

cos

lπ (v +α)

xr

 

 

+

lπ

cos

kπ (u +1)

sin

lπ (v +α )

yr

 

 

 

 

 

 

 

kl

 

kl

2

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2α

 

 

 

 

0

 

 

 

2α

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2α

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

114

Применяя к решению (8.80), (8.81) проекционную процедуру, получим следующую систему проекционных соотношений для определения амплитуд:

1

α

 

 

 

 

 

r

 

 

r

 

 

 

 

 

 

 

 

 

 

 

 

 

r

r

 

 

 

 

{rot (H&ts + H&

ζ s )jsωε0 gˆ (E&ts + E&ζ s )}erpre dvdu =

(8.82)

1 α

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

2π 1 α

 

r re

 

jsωt

 

 

 

 

 

 

 

 

 

 

 

 

 

=

 

 

∫ ∫

gˆδ e

pr

e

 

 

 

 

 

dvdudωt ,

 

 

 

 

 

π

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

1 α

 

 

 

 

 

 

 

)

 

 

 

 

 

 

 

 

 

 

 

 

}

 

 

{

 

(

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

α

 

rot

 

r

 

 

r

 

 

 

 

jsωε

 

gˆ

(

r

r

 

erм dvdu =

(8.83)

 

 

 

 

 

 

H&

+ H&

 

 

 

0

E

+ E

 

1

α

 

 

 

 

 

ts

 

 

 

ζ s

 

 

 

 

 

 

 

 

 

ts

ζ s )

pr

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

2π 1

 

 

α

 

 

r rм jsωt

 

 

 

 

 

 

 

 

 

 

 

=

 

 

 

 

gˆδ e

pr

e

 

 

 

 

dvdudωt ,

 

 

 

 

π

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

1

 

α

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

α

 

 

 

 

r

 

 

r

 

 

 

 

 

 

 

 

 

 

 

 

 

r

r

 

 

 

 

{rot (H&ts

+ H&ζ s )jsωε0 gˆ (Ets

+ E&ζ s )}ψ epr zr0dvdu =

(8.84)

1

α

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

2π 1

 

 

α

 

rr

 

 

e jsωt

 

 

 

 

 

 

 

 

 

=

 

 

 

 

gˆδ z

ψ

pr

e

 

 

 

dvdudωt ,

 

 

 

π

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

1

 

α

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

α

 

 

 

 

 

 

r

 

 

r

 

 

 

 

 

 

 

 

 

 

r

r

 

r

 

 

{gˆ 1rot (E&ts

+ E&ζ s )

+ jsωµ0 (H&ts + H&

ζ s )}hpre dvdu =

(8.85)

1

α

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

2π 1

 

α

rм re jsωt

 

 

 

 

 

 

 

 

 

 

= −

 

 

 

 

δ1

hpre

 

 

 

 

dvdudωt ,

 

 

 

 

π

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

1

α

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

α

 

 

 

 

 

 

r

 

 

r

 

 

 

 

 

 

 

 

 

 

r

r

 

r

 

 

{gˆ 1rot (E&ts

+ E&ζ s )

+ jsωµ0 (H&ts + H&

ζ s )}hprмdvdu =

(8.86)

1

α

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

2π 1

 

α

rм rм jsωt

 

 

 

 

 

 

 

 

 

 

= −

 

 

 

 

δ1

hpre

 

 

 

 

dvdudωt ,

 

 

 

 

π

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

1

α

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

α

 

 

 

 

 

 

r

 

 

r

 

 

 

 

 

 

 

 

 

 

r

 

 

 

 

{gˆ 1rot (E&ts

+ E&ζ s )

+ jsωµ0

H&ts}ψ prм zr0dvdu =

(8.87)

1

α

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

2π 1

 

α

rм м r jsωt

 

 

 

 

 

 

 

 

= −

 

 

 

 

δ1

ψ pr zoe

 

 

 

 

dvdudωt .

 

 

 

π

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

1

α

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Уравнения (8.82) – (8.83) образуют полную систему ОДУ первого порядка, определяющую комплексные амплитуды связанных волн

115

A&smne , A&skмl, B&smne , B&skмl,C&smn , D&skl , т.е. они представляют собой систему уравнений возбуждения динамических полей нерегулярного волновода с прямоугольным сечением, возбуждаемых негармоническими электрическими и магнитными токами источников. Полученные уравнения возбуждения, строго говоря, справедливы только в случае, когда граничные поверхности нерегулярного волновода не имеют изломов, т.е. az , bz , qz , hz не имеют разрывов. Прак-

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

Самосогласованные уравнения возбуждения продольно-нерегулярного волновода электронным потоком. Взаимодействие электронного потока с возбуждаемым в волноводе электромагнитным полем будем описывать системой безразмерных уравнений Максвелла для комплексных амплитуд s- гармоники рабочей частоты

r

r

r

;

r

r

rr

 

= 0;

(8.88)

 

rotB&s = jsWE&s +δ&s

rotE&s = − jsWB&s ;

[nE]

 

 

 

 

 

 

 

 

 

S

 

 

 

 

 

 

 

 

 

 

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

пучка

 

r

 

 

 

 

r

 

 

 

 

 

r

 

 

 

 

 

 

 

 

 

 

 

 

 

 

r

 

r

 

 

r

 

r

 

 

 

 

 

 

 

dγi βi

 

 

 

 

 

] + S

}/

β

 

;

 

(8.89)

 

= −{(E

+[β

B])

+[β

F

E

zi

 

 

 

dz

 

 

 

 

 

 

i

 

 

 

 

 

i

 

 

q

qi

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

drr

 

βr

 

dθ

i

 

 

 

W

 

 

 

 

 

 

 

 

 

r2

 

 

 

i

 

=

i

;

 

 

=

 

 

 

;

 

 

γi =1/ 1βi

 

,

 

dz

 

 

 

 

 

 

 

 

 

 

 

 

βzi

 

dz

 

 

βzi

 

 

 

 

 

 

 

 

 

 

 

 

r

 

 

 

r

 

 

 

 

 

 

 

 

 

 

r

 

 

r

(rr i

, zi )e jsθi )

 

 

E = Re(E&s (rr i , zi )e jsθi )

 

B = Re(B&s

 

 

 

 

 

s

 

 

 

 

 

 

 

 

 

 

 

 

 

s

 

 

 

 

 

 

r

(0)

 

r

r

 

 

 

 

r

 

θ

(0)

= 2π(i 1) / Ne;

i =1...Ne .

β

= β0 ;

r

(0) = r 0

 

 

i

 

 

 

i

i

 

 

 

 

 

i

 

i

 

 

 

 

 

 

 

 

 

 

Здесь и далее приняты следующие безразмерные переменные (штрихом будем обозначать размерные величины, имеющие одинаковое написание с безразмерными):

′ ′ ′ ′ ′

′ ′ ′

 

ω0

 

, W =

ω

, θi =ωti ; ω

x, y, z,u,v,a,b,q,h = (x , y , z ,u ,v ,a ,b ,q ,h )

 

c

 

 

βi = vri

 

 

 

 

 

ω0

рабочая частота, ω0 – базовая,

/ c ,

 

vi

скорость крупных частиц, с -

скорость света,

116

r &e &м &

r

& e

&

м

&

 

r &e &м

r

e

 

м

 

 

(E, A , A ,C )

 

 

 

) Em ,

(B, B , B , D)= (H , B

 

, B

 

.,

= (E , A , A ,C

 

 

 

 

, D )cµ0 Em

r&

 

 

1

 

 

2π

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2eI0

 

 

 

 

 

 

 

J (z,ωt)e

jsωt

dωt, J

(z,t) - плотность тока в ЭП;

G0 =

 

 

δs =

 

 

 

 

 

 

 

;

ε ω

E

π

 

 

 

ε m c3

 

I0

0

0

 

m

0

 

Em = m0ω0 c e,

m0 ,e -

 

 

 

 

 

 

0 0

 

 

-

ток

пучка,

масса и

 

заряд электрона,

r

 

r

 

 

 

 

 

r

 

- индукция магнитостатического поля, Sq Erqi

 

 

 

 

F

= eB0 / m0ω0

, B0

 

- силовая со-

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

Используя проекционные соотношения (8.82-8.87) для уравнений (8.88), а также закон сохранения заряда для связи плотности тока электронного пучка с электронными траекториями, описываемыми уравнениями (8.89), получим следующую систему обыкновенных дифференциальных уравнений для безразмерных амплитуд возбуждаемых мод в нерегулярном волноводе:

 

&

 

 

 

 

 

&

e

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

π

 

 

 

 

 

&

e

 

 

 

az

 

1cs

 

 

bz

 

cs

 

 

 

 

 

 

 

 

Cspr jBspr

χpr

/ (sWaq)

+

 

 

 

Asmr m

 

 

Imp

+

 

 

 

Imp

+

 

 

 

 

 

2

 

a

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

m

 

 

 

 

 

 

 

 

 

 

 

a

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

&

 

 

 

 

 

qz

 

 

1cs

 

 

 

hz

 

 

 

 

cs

 

 

 

 

 

r

 

 

 

 

 

 

 

 

 

м

az

 

1cs

 

 

bz

 

cs

 

 

 

 

 

 

e

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

+ Aspn n

 

 

 

 

Inr

 

 

 

+

 

 

 

 

 

 

 

Inr

 

+

 

 

 

 

 

Askr

 

 

Ikp

 

+

 

 

 

 

 

Ikp

q

 

 

 

 

 

α q

 

α

 

 

 

a

 

 

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

k

 

 

 

 

 

a

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

&

 

 

 

 

 

q

 

 

 

1cs

 

 

 

h

 

 

 

 

 

 

cs

 

 

 

 

 

C

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

м

 

 

 

z

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

pAspl α

 

Ilr

 

 

 

+

 

z

 

Ilr

 

+ I pr ;

 

p = 1..M ; r = 1..N; .

 

 

 

q

 

 

 

 

q

 

 

 

 

 

 

 

l

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

jG

 

 

 

1

 

 

 

Ne

 

 

 

 

 

 

pπ

(u +1)

 

 

 

 

rπ

(v +α)

 

θ

 

 

 

 

 

 

 

 

 

I&Cpr

= −

 

 

 

 

 

 

0

 

 

 

 

 

 

 

sin

 

 

 

 

 

 

 

i

 

 

 

 

 

 

 

 

sin

 

 

 

 

 

i

 

 

 

 

e

js

i ;

 

 

 

 

 

 

sWaqα

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

2α

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Ne i=1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

&e

 

 

 

 

 

м

 

 

 

e

 

 

 

 

 

м

 

 

 

м

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dAspr

=

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

;

 

 

p

= 1..M ; r = 1..N;

 

 

 

 

 

 

 

 

 

 

dz

αχpr

f pr

+γ pr

f pr

 

χpr

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

&

м

 

 

 

 

e

 

 

e

 

 

 

 

 

 

e

 

 

 

м

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dAspr

=

 

 

 

 

 

 

 

 

 

 

 

 

 

 

;

 

 

p

= 0..K; r = 0..L;

 

 

 

 

 

 

 

 

 

 

 

 

dz

γ pr f pr

+αχpr

f pr

 

χpr

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

e

 

 

 

 

 

 

 

&

e

 

 

 

2

 

 

 

 

 

 

&

 

 

 

 

 

 

 

 

e

rπ

 

 

 

 

 

&

м

 

2

a

 

 

 

1cs

+

b

 

 

cs

 

+

 

f pr

= jsWBspr

χpsα + Cprαχpr

 

 

 

 

 

 

Askr χkr

 

z

Ikp

 

 

z

Ikp

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

k

 

 

 

 

 

 

 

q

 

 

 

 

 

 

 

q

 

 

 

 

 

 

 

 

pπ

 

 

 

 

 

 

&

м

 

 

 

 

2

 

qz

 

 

 

 

1cs

 

 

 

 

hz

 

 

 

cs

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

+

 

 

 

 

αAspl

χpl

 

 

 

 

 

 

αIlr

 

+

 

 

 

 

 

 

Ilr

;

 

p = 1..M ; r = 1..N;

 

 

 

 

 

2

 

 

 

a

 

 

a

 

 

 

 

 

 

 

 

 

 

 

 

 

 

l

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

f

м

 

 

 

 

 

 

 

 

&м

 

χ

2

 

 

 

 

 

 

&

 

 

 

 

γ

e

+

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

pr

= − jsWB

 

 

 

 

 

 

α C

spr

pr

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

spr

 

pr

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

α pπ

 

 

 

&

м

 

2

 

a

 

 

1cs

 

 

 

 

 

b

 

cs

 

 

 

rπ

 

 

 

 

&

м

2

q

 

 

 

cs

 

 

h

 

cs

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

z

 

 

 

 

 

 

 

 

 

 

z

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

z

 

 

 

 

 

 

 

 

z

 

 

+

 

 

 

Askr χkr

 

 

Ikp

 

 

+

 

 

 

Ikp

 

+

 

 

 

 

 

 

Aplχpl

 

 

αIlr

+

Ilr ;

 

2

 

q

 

 

 

 

q

 

 

 

 

2

 

a

a

 

 

 

k

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

l

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(8.90)

(8.91)

(8.92)

p=0..K;

117

r=0..L;

 

 

χ

2

 

=

pπ 2

+

 

rπ

2 ;

 

 

 

 

 

 

 

 

 

 

χ

 

= χe

 

χ

 

м

α2 γ e

γ м ;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

pr

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

2α

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

pr

 

 

 

 

 

 

 

pr

 

 

 

pr

 

 

 

 

 

 

 

 

 

 

 

pr

pr

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

e

 

prπ

 

 

 

 

 

a

 

q

 

 

 

если Epr

 

присутствует

 

 

 

 

 

 

 

 

e

 

 

 

 

 

 

 

 

q

 

pπ 2

 

 

a rπ 2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

γ pr =

4

 

 

 

 

 

 

q

 

 

 

a

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

;

 

 

 

 

 

χpr

=

 

 

 

 

 

 

 

 

 

 

+

 

 

 

 

 

 

 

;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a

 

2

 

 

 

 

2α

 

 

 

 

 

 

 

 

 

 

 

 

 

0 если Epr отсутствует

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

q

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

m

 

prπ

 

 

 

 

 

a

 

q

 

 

 

если H pr присутствует

 

 

 

 

 

 

 

м

 

 

 

 

 

 

a

pπ 2

 

 

 

q rπ 2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

γ pr =

4

 

 

 

 

 

 

q

 

 

 

a

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

;

 

 

 

 

χpr

 

=

 

 

 

 

 

 

 

 

 

+

 

 

 

 

 

 

 

;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

2α

 

 

 

 

 

 

 

 

 

 

 

 

 

0 если H pr

 

отсутствует

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

q

 

 

 

 

 

 

 

a

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

&

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

&м

 

 

2

 

(1

 

 

 

 

 

2

 

 

 

 

 

2

)/ aq +

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(8.93)

Dspr = − j

 

(sW ){Aspr

χpr

+ bz

 

+ hz

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

м

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

2cc

 

 

 

 

2azbz

 

 

 

1cc

 

 

 

π

 

 

 

 

 

&

м

 

 

 

 

 

az 1sc

 

bz

 

 

sc

 

 

 

 

 

 

+

 

 

&

 

 

 

 

 

2 az

 

 

 

 

 

 

 

 

 

 

 

 

dAskr

 

 

 

 

 

 

 

 

 

 

 

 

 

Askr χkr

aq

 

Ikp

 

+

 

 

aq

 

 

 

Ikp

 

 

 

 

 

2

 

 

 

 

dz

 

 

 

k

 

 

q

 

 

Ikp

+

q

 

Ikp

 

+

 

 

 

 

 

 

 

k

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

м

 

 

 

 

 

2

 

 

 

 

q

2

 

 

 

2cc

 

 

 

2q h

 

 

 

 

1cc

 

 

 

 

 

 

 

 

 

 

&

м

 

 

 

q

 

 

 

1sc

 

 

 

h

 

 

 

 

sc

 

 

 

 

 

 

 

 

 

 

 

&

 

 

 

 

 

 

 

 

 

 

α

 

 

 

 

 

 

 

 

 

 

 

 

π dAspl

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Ilr

 

+

 

 

 

 

 

 

 

 

 

Ilr

 

 

 

 

 

 

 

 

 

 

 

 

 

 

l

 

 

 

 

 

Ilr +

 

 

 

 

 

Ilr

 

 

 

+ Asklχplα

 

 

 

z

 

 

 

 

 

 

 

 

 

z

 

z

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

z

z

 

 

 

 

 

 

 

 

 

l

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

aq

 

 

 

 

 

 

 

 

 

 

 

aq

 

 

 

 

 

 

 

 

2

dz

a

 

 

 

 

 

 

 

 

aα

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

rπ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

&e

 

 

 

 

 

a

 

 

 

1sc

 

 

 

 

b

 

 

 

 

sc

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dA

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

C&smr

 

 

 

 

 

 

smr

 

 

 

 

 

 

 

z Imp +

 

 

 

z

 

Imp

 

 

+

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2α

 

 

 

 

 

 

 

 

 

 

 

 

 

 

q

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dz

 

 

 

 

 

 

 

q

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

m

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

pπ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

&e

 

 

 

 

q

 

α

 

 

 

 

1sc

 

 

 

 

 

h

 

 

 

sc

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dAspn

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

+

 

 

 

 

 

C&

spn

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

z

 

 

I

nr

 

 

+

 

 

 

z

I

nr

 

 

;

 

p

 

= 0..K; r = 0..L;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dz

 

 

 

 

 

a

 

 

 

 

 

 

 

 

 

 

a

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

&e

 

 

 

 

 

 

 

 

 

 

jsW π

 

2

 

 

 

e

 

q

 

 

 

 

 

 

 

 

2

 

 

 

 

2

 

 

 

 

a

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

r

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dBspr

 

 

 

 

 

 

 

 

&

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(1 + bz )p

 

 

 

 

 

 

 

 

 

(1 + hz

)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dz

 

 

 

 

 

=

 

χ

2

 

 

 

 

4

Aspr

 

a

 

 

+

 

q

 

α

 

2

+

 

 

 

 

 

 

 

 

 

 

 

 

 

(8.94)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

pr

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

&

м

 

 

 

 

 

pr q

 

 

 

 

 

 

 

 

 

2

 

 

 

 

a

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

&

 

e

 

 

 

 

 

 

 

 

q

 

 

 

 

 

 

 

 

 

 

 

 

 

1cc

 

 

 

2

 

 

 

2cc

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(

1 + bz )

 

 

 

 

(

 

 

 

 

 

 

 

 

 

) +

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(2azbz Imp

 

 

 

 

 

 

 

 

 

 

)+

 

 

+

 

Aspr

 

 

α

 

 

 

 

q

1 + hz

 

Asmr mp

a

+ az Imp

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

m

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

&

м

 

pr q

 

 

 

 

 

 

 

 

 

 

1cc

 

 

 

 

2

 

 

 

2cc

 

 

 

 

 

 

 

 

 

 

 

 

&

e

 

 

 

 

 

 

 

a

2qzhz

1cc

 

 

 

 

2

 

2cc

 

 

 

+ Askr

 

 

 

 

 

 

 

 

 

(2azbz Ikp

+ az Ikp

)+ Aspnrn

 

 

 

 

 

 

 

 

 

 

 

 

 

Inr

 

+ qz Inr

 

 

 

α a

 

 

 

 

α

 

 

 

 

 

 

 

 

 

 

k

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n

 

 

 

 

 

 

 

 

 

 

 

 

 

q

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

&

м

 

 

 

 

 

 

 

a

 

 

 

 

 

 

 

 

 

 

 

1cc

 

 

 

 

2

 

 

 

 

 

2cc

 

 

 

 

 

 

2

 

 

 

 

 

 

 

&

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

sc

 

 

 

 

 

 

1sc

 

 

 

 

 

 

 

Aspl pr

q

(2qzhz Ilr

 

 

 

+ qzαIlr

 

 

 

)

 

 

 

Csmr pq(bz Imp + az Imp

 

)

 

 

 

 

 

 

 

 

 

π

 

 

 

 

 

 

 

l

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

m

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

&

 

 

 

 

 

r

 

 

 

 

 

 

 

 

 

sc

 

 

 

 

 

 

 

 

 

1sc

 

 

 

 

 

 

 

 

 

 

 

&

e

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

&

м

 

 

2

 

 

 

 

&

Be

 

 

 

 

 

 

 

 

 

 

Cspn

 

 

 

 

 

a(hz Inr

+αqz Inr

 

 

)+ Asmnγmn + Asklγkl

 

 

+ I pr

;

 

 

 

 

 

 

π

 

 

α

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

mn

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

kl

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

118

1

mr

 

cs sc

+

np

sc cs

 

cs

1sc

 

sc

1cs

)+

γmn = bzhz

 

 

Imp Inr

α

Imp Inr

+ bzqz (mrImp Inr

+ npImp Inr

 

α

 

 

 

 

 

 

 

 

 

 

 

 

 

 

mr

1cs

 

sc

+

np

1sc

cs

 

1cs

1sc

 

1sc

1cs

);

 

+ azhz

Imp

Inr

α

Imp

Inr

+ azqz (mrImp

Inr

+ npImp

Inr

 

α

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

γ 2

= b h

lr

I cs I sc

kpI sc I cs

 

 

 

+ b q

 

 

 

lr

 

I cs I1sc

kpαI sc I 1cs +

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

z

 

 

 

 

 

 

 

 

 

kl

 

 

 

 

 

 

 

 

z z

α

 

 

kp

lr

 

 

 

 

 

 

 

 

 

 

kp

 

lr

 

 

 

 

 

 

 

z

 

 

 

 

α

 

kp

 

 

lr

 

 

 

 

 

 

 

 

 

 

 

kp lr

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

+ a h

 

 

lr

I 1cs I sc kpI1sc I cs

 

+ a

 

q

 

lr

 

I

1cs I

1sc

kpαI

1sc I 1cs

 

;

 

 

 

 

 

 

 

 

 

 

 

 

z

 

 

 

 

 

 

lr

 

 

 

 

 

 

 

 

 

 

z z

 

 

α

2

 

 

kp

 

 

 

 

lr

 

 

 

 

 

 

 

 

 

kp

 

 

 

 

lr

 

 

 

 

z

 

 

 

 

α

 

 

 

 

kp

 

 

 

 

 

 

 

 

 

 

 

 

kp lr

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Be

 

 

 

G 1

 

 

 

Ne

 

 

 

1

 

β

xi

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

pπ

 

 

 

 

 

 

 

 

 

 

pπ (ui +1)

 

 

 

 

 

 

rπ

(vi

+α )

 

 

I&pr

= −

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

uiaz

bz

 

 

 

 

 

 

 

 

 

cos

 

 

 

 

 

 

 

 

 

 

 

 

sin

 

 

 

 

 

 

+

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

βzi

 

 

 

2

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

2α

 

 

 

 

 

 

χprα Ne i=1

 

a

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

β

yi

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

π

 

 

 

 

 

 

 

pπ

 

(

u +1

)

 

 

 

 

 

 

rπ

(

v +α

 

)

 

 

 

 

 

 

 

 

 

 

 

+

 

 

 

 

 

 

 

 

viqz

 

hz

 

r

 

 

 

sin

 

 

 

 

 

 

 

 

i

 

 

 

 

 

 

cos

 

 

 

 

 

i

 

 

 

 

 

ejsθi

;

 

 

 

 

 

 

 

q

 

 

 

βzi

 

 

2α

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

2α

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

p =1..M ; r =1..N;

 

 

 

 

&м

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

jsW π

2

 

 

 

e

 

pr

 

 

q

 

 

 

 

 

 

 

 

 

 

2

 

 

 

a

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

dBspr

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

&

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(1 + bz

)

 

 

 

(1

 

 

 

 

 

 

 

) +

 

 

 

 

 

 

 

 

 

dz

 

 

 

 

= −Dspr

 

 

 

 

2

 

 

4

 

Aspr

 

α

 

 

 

 

 

q

+ hz

 

 

 

(8.95)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

χpr

 

 

 

 

 

 

 

 

 

 

 

 

a

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

&

м

 

 

 

q

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

r

2

 

 

 

 

 

 

a

 

 

 

 

 

 

 

 

2

 

 

 

 

2

 

 

 

 

 

 

 

 

 

&

e

 

 

q mr

 

 

 

 

1cc

 

 

 

 

2cc

2

 

 

 

 

+ Aspr

 

 

 

 

 

(

1

+ bz

)

α

2

 

+

 

q

(1 + hz )p + Aspr

 

a

α

 

(Imp 2azbz + Imp az )+

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

m

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

&

 

м

 

 

 

r2

 

 

q

 

 

 

 

1cc

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2cc

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

&

e

 

 

 

 

a

 

 

1cc

 

 

 

 

 

 

2cc

2

 

 

 

 

 

+ Askr

 

 

 

 

 

 

 

 

 

 

 

(Ikp

 

2azbz + Ikp

 

az

)Aspnnp

 

 

 

(Inr

 

 

2qzhz + Inr αqz

)+

 

 

 

 

 

α

2

 

a

 

 

q

 

 

 

 

 

 

 

 

k

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

&м

 

 

 

 

 

 

2

 

 

 

a

 

 

 

 

1cc

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2cc

 

 

 

2

 

 

 

 

 

 

 

 

2

 

 

 

 

&

 

 

 

 

qr

 

 

 

 

sc

 

 

1sc

 

 

 

 

 

 

+ Aspl p

α

q

(Ilr

 

 

2qzhz + Ils

 

 

αqz

)

 

 

 

Csmr

 

 

 

 

(Impbz + Imp az )+

 

 

 

 

 

 

 

 

 

π

α

 

 

 

 

 

 

 

 

 

 

l

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

m

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

+

 

2

 

 

C&spnap(Inrschz + Inr1scαqz )Asmne γmn3

 

Askмlγk4l

 

+ I≺

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

π

 

 

 

 

 

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

mn

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

kl

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

γ 3

 

= b h

mpI cs

 

I sc

nr

I sc

I cs

 

 

+ a

 

h

mpI

1cs I sc

nr

I

1sc I cs

 

+

 

 

 

 

 

 

 

 

 

z

 

 

 

 

 

 

mn

 

 

z

z

 

 

 

 

 

 

 

 

 

 

mp

 

 

nr

 

 

 

α2

 

 

mp

 

 

nr

 

 

 

 

 

 

 

 

z

 

 

 

 

 

 

 

 

 

 

mp

 

nr

 

 

 

α2

 

mp

 

nr

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

+αq

b

 

 

mpI cs I

1sc

nr

I sc

I

1cs

 

 

+

αa

q

 

 

 

mpI

1cs I

1sc

 

 

nr

I 1sc I 1cs

 

;

 

 

 

 

 

nr

 

nr

 

 

z

 

 

 

 

 

 

 

 

 

 

 

z

 

 

 

z

 

 

 

 

 

 

 

 

 

mp

 

 

 

 

 

α2

 

 

 

 

mp

 

 

 

 

 

 

 

 

 

 

z

 

 

 

 

 

 

 

 

mp

 

 

 

nr

 

 

 

 

 

α2

mp

nr

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

γk4l

= bzhz (lpIkpcs Ilscr

 

+ krIkpsc Ilcsr )

+

azhz

 

(lpIkp1cs Ilscr + krIkp1sc Ilcsr )+

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

α

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

α

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

+ qzbz (lpIkpcs Il1scr

+ krIkpsc Il1csr )+ azqz (lpIkp1cs Ilscr

+ krIkp1sc Ilcsr );

 

 

 

 

 

 

119

 

G 1

Ne

 

1

 

β

xi

rπ

 

pπ (ui +1)

 

rπ (vi +α )

 

I&pr

=

0

 

 

 

 

 

uiaz bz

 

cos

 

sin

 

2

 

 

 

βzi

2α

2

2α

 

 

χprα Ne i=1

a

 

 

 

 

 

 

1

βyi

pπ

 

pπ (ui +1)

 

rπ (vi +α )

jsθi

 

 

 

viqz hz

 

sin

 

cos

 

e

 

q

βzi

2

2

2α

 

 

 

 

 

 

 

 

p =0..K; r =0..L ;

здесь ui ,vi ,βxi ,βyi ,βzi - координаты и скорости частиц в сечении z,

 

1

iπ (x

+1)

 

 

jπ (x +1)

 

 

 

 

Iijcs = cos

sin

dx;

Iijsc = I csji ;

2

 

 

 

 

1

 

 

 

 

 

 

2

 

 

 

 

 

 

iπ (x +1)

 

 

jπ (x +1)

 

 

 

 

1

 

 

 

 

 

 

Iij1cs

= xcos

sin

dx;

Iij1sc = I 1csji ;

 

 

 

 

1

 

 

 

2

 

 

 

 

 

2

 

 

 

 

 

 

 

iπ (x

+1)

 

 

jπ (x +1)

 

 

 

1

 

 

 

 

 

 

Iij2cs

= x2 cos

sin

dx; Iij2sc = I 2csji ;

 

 

 

 

1

 

 

 

2

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

аналогично обозначаются интегралы Iij1cc , Iij2cc , Iij1ss , Iij2ss .

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

При задании граничных условий для уравнений (8.90)-(8.95) используем тот факт, что на регулярных участках волновода волновое поле представляет собой сумму прямой и обратной волн вида:

 

 

 

e&mne+

 

 

 

 

 

 

 

r

r

 

e

 

+ e&mne

 

 

 

 

r

r

 

e

 

 

 

 

 

ejK zmn (E&

+B&+ )mn

e+ jK zmn (E&

B&

)mn

+

(8.96)

 

 

 

mn

 

 

 

 

 

 

 

 

 

 

 

 

)

 

 

 

 

 

 

 

 

 

 

)

 

 

 

 

 

 

 

 

м+

e

jKzkl

 

r&+ r&

+

м

 

м

e

+ jKzkl

 

r&r&

м

 

 

 

 

 

+ekl

 

 

 

 

(E B

 

 

+ ekl

 

 

(E B

 

 

 

 

 

 

 

 

kl

&

 

 

 

 

 

 

 

 

 

 

 

 

 

kl

&

 

 

 

 

 

 

 

 

 

kl

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

&e±

&м±

- постоянные амплитуды,

r&± r&

±

e, м

- собственные функции Eij

где emn , ekl

(E B

 

)

 

и Hij

регулярного волновода, Kzij =

W 2 Ktij2

- продольное волновое число,

K

 

 

 

iπ

 

 

 

2

 

 

 

jπ

 

 

 

2

- поперечное волновое число.

 

 

=

 

 

 

 

 

 

 

 

+

 

 

 

 

 

 

 

 

 

D D

 

 

 

D

 

D

 

 

 

 

 

tij

 

 

 

 

 

 

 

yp

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

xp

 

 

ym

 

 

 

 

ym

 

 

 

 

 

 

 

 

 

 

 

 

 

Используя соотношения (8.75) а также предполагая, что при z 0 волновод регулярный, запишем начальные условия для амплитуд:

120

A&mne (0) = (e&mnee&mne+ ) Kzmn / Ktmn (0),

&e

 

ee+

(8.97)

Bmn (0) = (emn + emn ) W / Ktmn (0),

 

&

&

 

A&kмl(0) = (e&kмl+ + e&kмl) W / Ktkl (0),

B&kмl(0) = (e&kмl+ e&kмl) Kzkl / Ktkl (0),

Безразмерная мощность, переносимая волновым полем через поперечное сечение волновода в выбранных переменных имеет вид

 

 

2

&

e

&

e*

2

&

м

&

м*

 

(8.98)

P(z)= −α ∑ ∑χmnreal (Asmn Bsmn )+ χklreal (AsklBskl ) .

s

mn

 

 

 

 

 

kl

 

 

 

 

 

 

Исходя из представления (8.96) на регулярных участках, а также в точках волновода, где az=bz=qz=hz=0, мощности прямой и обратной волн в выбранных безразмерных переменных выражаются следующим образом:

 

 

 

 

 

 

 

 

 

 

j

 

&e

 

 

 

 

 

j

&e

*

P

 

 

 

 

 

 

 

 

 

dA

 

 

 

 

dB

 

+

±

= −α ∑ ∑χmn real Asmn ±

 

 

 

 

 

Bsmn ±

 

 

 

 

 

 

 

 

 

2

 

 

 

&e

 

 

 

 

smn

 

&e

 

 

 

 

smn

 

 

 

 

s

mn

 

 

 

 

Kzmn

 

 

dz

 

 

 

 

Kzmn

dz

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

j

 

&м

 

 

 

 

j

&м

*

 

 

 

 

 

 

2

 

&м

dA

&м

 

dB

 

 

 

 

 

 

 

 

 

 

 

 

 

skl

 

 

 

 

 

skl

 

 

 

 

 

 

 

+χkl real Askl ±

 

 

 

 

Bskl ±

 

 

 

 

 

 

 

 

 

 

(8.99)

Kzkl

 

dz

 

Kzkl

dz

 

 

 

 

kl

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Эффективность взаимодействия электронного потока с возбуждаемым электромагнитным полем оценивается величиной электронного КПД

 

1

Ne

γi (z)γi (0)

 

 

ηe =

;

(8.100)

 

1 γi (0)

 

Ne i=1

 

 

или волнового КПД

ηv =

P(z)P(0)

;

ηv± =

P± (z)P± (0)

;

G0 (γ0 1)

G0 (γ0 1)

 

 

 

 

Расчет эффективности по разным формулам позволяет контролировать погрешность вычислений, т.к. при отсутствии потерь в стенках волновода из-

121

за сохранения баланса энергии должно выполняться условие ηe =ηv , а также на регулярных участках ηv =ηv+ +ηv.

Ниже приведены примеры моделирования приборов и устройств СВЧ на основе развитой теории.

ЛБВ-О на нерегулярном гофрированном волноводе прямоугольного се-

чения. Схема ЛБВ приведена на Рис.8.3. Формируемый электронной пушкой

Рис. 8.3.

ленточный электронный поток, проходя в зазоре волновода с гофрированны-

 

y

 

 

Dyp

 

 

 

x

Dxm

Dym

Dxp

Рис.8.4

ми стенками (Рис. 8.4), взаимодействует с усиливаемой Е11-волной и другими связанными с ней волнами, после чего осаждается на коллекторе. Ввод сигнала и вывод СВЧ энергии осуществляется через согласованныеr переходы. В

области взаимодействия создано однородное магнитное полеF = F0 zr0 , обес-

печивающее минимизацию поперечных составляющих скорости пучка. Профиль гофра задавался в виде:

D

= −D ;

D

ym

= −D

yo

h

ym

(z)sin2 (n π z / L) ;

(8.101)

xm

xo

 

 

 

 

 

 

 

ym

 

D

= D ;

D

yp

= D

yo

+ h

yp

(z)sin2 (n π z / L) .

 

xp

xo

 

 

 

 

 

 

yp

 

В математической модели ЛБВ-О приняты трехмерныеr уравнения движения крупных частиц (8.89) с начальными условиями βi0 = β0 zr0 , rr0i = 0 . Начальные условия для амплитуд (8.87) задавались исходя из того, что на вход

122

подается сигнал в виде падающей Е11 - волны с амплитудой e11e+re = e0 , ампли-

туды всех остальных распространяющихся волн задавались равными нулю. Тем самым обеспечивается полное согласование на входном конце лампы. Исходя из условия селекции рабочей E11-волны на базовой частоте, были выбраны безразмерные (отнесенные к λ0 2π ) значения Dxo=10, Dyo=2 . Опти-

мизировались: глубина гофра hym(z)=hyp(z)=hy(z), длина гребенки L, количество периодов nym=nyp=ny, , ток пучка I0 амплитуда сигнала e0. Анализ показывает, что в волноводе при таких размерах распространяется 12 типов волн, однако в результате расчетов выяснилось, что основной "паразитной" волной в рассматриваемой ЛБВ-О является волна НЕ11, которая связана в нерегулярном волноводе с волной Е11 и оказывает существенное влияние на процесс усиления. Все другие распространяющиеся типы волн не возбуждаются.

Был произведен поиск оптимизированных по величине волнового КПД вариантов ЛБВ-О в диапазоне напряжений электронного пучка (0.5≤β00.85). Анализ полученных результатов показал необходимость учета НE11 - волны. Оказалось, что предварительные расчеты по упрощенной модели с учетом только Е11- волны дают значительную погрешность. При β0 >0.7 были най-

дены варианты с КПД 36% на регулярной гребенке и с КПД 57% на нерегулярной. Ниже приведены два варианта.

Вариант1(регулярная гребенка):

β0 =0.75,

I0 =100A,

L =71, ny = 40, hy =1.34,

e0 = 0.0025,

η = 32%,

Ku =13 Дб.

Вариант2: (нерегулярная гребенка)

β0 =0.85,

I0 = 386 A,

L = 88, ny = 40, η = 57%, Кu = 23 Дб.

В этом варианте глубина гофра изменяется в пределах 1.005 < hy <1.345.

За счет оптимального профилирования глубины гофра эффективность лампы может быть увеличена примерно в 1.5-2 раза.

Таким образом, полученная математическая модель позволяет производить расчеты и оптимизацию ЛБВ-О с нерегулярным профилем гофра замедляющей системы. Проведенные исследования показали высокую эффективность ЛБВ-О с замедляющей системой в виде гребенки с оптимально профилированным гофром. Причем, с увеличением напряжения электронного пучка КПД такой лампы возрастает.

Заграждающий фильтр на основе периодического гофра. Рассмотрим расчет заграждающего фильтра для волны Н10 - типа. Геометрия отрезка волновода, выполняющего функции фильтра представлена на Рис. 8.4. Профиль гофра задавался формулами (8.101).

123

При расчете фильтра граничные условия выбираются следующим образом. На входе задаются падающая Н10 волна с амплитудой e&01м+ =1 + j0 и от-

раженная e&01м= e0 e jϕ0 . Амплитуды всех других волн задавались равными

нулю. На выходном конце при z=L задается условие отсутствия встречной волны P(L)=0 и минимизируется величина проходящей мощности P+(L). В

уравнениях (19)-(24) полагается I&C = I&Be = I&=0 .

Задавались значения Dx0 , Dy0 и методом пристрелки подбирались зна-

чения e0, ϕ0, L, hy, ny из условия P(L)=0 и min P+(L).

Расчеты, в частности, показали, что при Dx0 < π, Dy0 < π, (т.е. пока не распространяются волны Н12, Е12), рассматриваемый периодический фильтр (hy=const) при ny=10 обеспечивает 25ДБ затухания. При этом волны Н01, Н11,

Е11 не возбуждаются.

Один из вариантов такого фильтра приведен ниже:

Dx0 = 1.75, Dy0 = 2.75, L =70.5 , nv = 10, hy=1.38, P+(L)/P+(0) =

0.0034.

На Рис. 8.5 приведено распределение мощностей прямой и встречной волн (8.99) P+(z), P(z). Полоса заграждения такого фильтра уменьшается с увеличением ny.

P+, P-

2

 

 

 

 

 

1,5

 

 

 

 

 

 

 

P+

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

0,5

 

 

 

 

 

 

0

 

 

 

 

 

 

-0,5

 

P-

 

 

 

 

-1

 

 

 

 

 

 

 

 

 

 

 

-1,5

 

 

 

 

 

 

-2

 

 

 

 

 

 

0

0,2

0,4

0,6

0,8

1

 

 

 

 

 

 

z/L

 

 

 

Рис 8.5.

 

 

 

Расчет дисперсионных характеристик периодического волновода. При расчете дисперсионных характеристик волн в периодических волноводах искомое решение для поля в волноводе подчиняется условию Флоке:

r

z + kd ),

r

r

r

 

(E&(x, y,

B&(x, y,

z + kd ))= (E&(x, y, z),

B&(x, y, z)) e jkϕ0 .

(8.102)

124

здесь d - период, 0 ϕ0 π - набег фазы нулевой пространственной гармоники поля. Постоянные распространения пространственных гармоник определяются через ϕ0 и d по формуле:

 

 

βφn =

 

d

 

 

 

; n =0,

±1, ±2, ...,

βφn = vφn / c

(8.103)

 

 

ϕ0 + 2nπ

vφn - фазовая скорость гармоники.

 

 

 

 

Hkl)

При расчете зависимости ϕ0(W) для волны определенного типа (Emn или

граничные

 

условия

при

z=0

задаются

следующим

образом:

&+

=1,

&

= eor . Методом пристрелки подбираются eor, ϕ0, для которых вы-

eij

eij

 

полняются условия при z=d , вытекающие из (8.102):

 

&

 

(d )

&

 

 

 

 

jϕ

0

;

&

 

 

&

jϕ0

 

 

 

Aij

= Aij (0) e

 

 

Bij (d )=

Bij (0) e .

 

 

 

 

 

Анализ решений этой задачи, в частности, показал, что для отрезка пе-

риодического

волновода нельзя

одновременно

удовлетворить

условию

(8.102)

и

 

 

условиям

полного

согласования

на входе и

выходе

&

 

 

&

 

(0)

&

 

 

&

(L)=0 .

Следовательно,

решение, удовлетворяю-

A

(0)= B

 

 

= A

(L)= B

 

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

В качестве примера расчета на Рис. 8.6. приведена зависимость ϕ0(W) для периодического гофрированного волновода, используемого в приведенном выше варианте 1 ЛБВ-0. Верхняя кривая соответствует расчету для Е11+НЕ11 волны, нижняя - для Е11 без учета НЕ11. Анализ показывает, что при

W& =1 нулевая пространственная гармоника Е11+НЕ11 имеет замедление βφ0 = 0.71, эта же гармоника Е11-волны имеет замедление βφ0 =0.755 . Таким

образом, усиление в этом варианте реализуется при синхронизме именно с пространственной гармоникой комбинированной Е11+НЕ11 волны.

Развитая строгая теория произвольно-нерегулярного волновода с прямоугольным сечением позволяет корректно решать широкий круг задач анализа и оптимизации электронных приборов и устройств СВЧ. Приведенные в статье примеры демонстрируют недопустимость использования в этих задачах упрощенных теорий, в которых не учитывается связь различных типов волн в нерегулярных периодических и непериодических волноводах.

125

3

 

 

 

 

 

 

ϕ0

 

 

 

 

 

 

2,5

 

E11+HE11

 

 

 

2

 

 

 

E11

 

 

 

 

 

 

 

 

1,5

 

 

 

 

 

 

1

0,85

0,9

0,95

1

1,05

1,1

0,8

 

 

 

 

 

 

W

 

 

 

Рис.8.6.

 

 

 

126