Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

Кураев А.А., Попкова Т.Л., Синицын А.К. Электродинамика и распространение радиоволн

.pdf
Скачиваний:
342
Добавлен:
02.05.2014
Размер:
5 Мб
Скачать
☆

Систему функций {Φmβ }будем использовать для разложения искомого

поля ϕ. Фактическое построение Φmβ , в частности вычисления Rnm(α) для мно-

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

Для первой задачи искомое поле в канале с номером α обозначим символом ϕα и представим разложением по функциям (10.32):

 

 

∞

N

 

 

 

 

∞

N

 

 

 

 

 

 

ϕα = ∑ ∑Cm(β )Φαβm

= ∑ ∑Cm(β ) ×

 

(10.33)

 

 

m=1 β =1

 

 

 

 

m=1 β =1

 

 

 

 

 

 

∞

 

 

jΓ

Z

 

 

 

 

− jΓ

Z

 

×

 

∑Rmα ϕ

 

e

+δαδϕ

 

e

 

 

α

n(α ) α

α

m(α ) α

 

 

n( )

n( )

 

 

 

 

 

m( )

 

 

 

 

n=1

 

 

 

 

 

 

 

 

 

 

 

Для второй задачи разложения (10.33) модифицируется следующим образом:

∞N

ϕ= Φ0 + ∑ ∑Cm(β )Φmβ .

m=1 β =1

где Φ0 - частное решение уравнения (3.1), удовлетворяющее граничному условию (10.31) на поверхности S0.

Вычисление неизвестных коэффициентов Cm(β) производится при помощи граничных условий на поверхностях Sα. Будем использовать для этого метод моментов (для определенности рассмотрим его на примере задачи (1), заключающийся в переходе от условий (10.30) при zα=lα к цепочке соотношений ортогональности

∫ϕαϕn(α )dS = 0, α =1,2,...N;

n =1,2,...

(10.34)

Sα

 

 

 

 

 

Подставляя в эти равенства (3.5) и интегрируя по Sα, получаем однородную

систему линейных алгебраических уравнений для Cm(β)

]= 0 .

 

∞ N

 

 

 

∑ ∑Cm(β )[Rnm(α)e2 jΓn(α )lα +δαβδmn

(10.35)

m=1 β

 

 

 

 

 

Спектральный параметр p входит в (3.7) через Rm

 

и Γ

 

. Характеристиче-

n(α)

n(α)

 

 

ские числа p отыскиваются из условий существования нетривиального решения системы (10.35).

Справедливо утверждение (теорема эквивалентности), что собственные значения задачи 1 и характеристические числа системы (10.35) совпадают, а собственные функции задачи 1 и нетривиальные решения системы (10.35)

связаны соотношением (10.33).

Следует отметить, что как изложенная выше схема решения, так и теорема эквивалентности остаются в силе и для других типов граничных условий: Неймана, 3-го рода и т.д. Однако в последнем случае, если коэффициент в граничном условии комплексный, рассматриваемая краевая задача стано-

171

вится несамосопряженной. Это может потребовать добавления в (10.33) присоединенных волн. Аналогичные схемы можно построить и в случае краевых задач для системы уравнений Максвелла.

Для обоснования различных вычислительных схем на основе функций (10.32) существенны вопросы о полноте и линейной независимости таких функциональных систем. Рассмотрим эти свойства.

Пусть все lα таковы, что некоторое действительное число p не является собственным значением задачи 1 в области V. Тогда система функций {Φmβ (p)}линейно независима на S' и полна в L2(S').

Свойство полноты будет доказано, если из условий

 

N

 

∫ϕΦmβ *dS = ∑ ∫Φmβα* dS = 0,

β =1,2,...N; m =1,2,...∞. (10.36)

S ′

α =1Sα

 

вытекает ϕ=0 на S'. Подставим в (10.36) разложение (10.33) и получим, что набор равенств (10.36) эквивалентен системе линейных алгебраических уравнений (символ” ” означает комплексное сопряжение).

N ∞

−2 jΓn(α )lα

 

 

∑∑dn(α)Rnm(α* )e

+ dn(β )δαβδmn = 0 .

(10.37)

α =1n=1

относительно

dn(α) = ∫ϕϕn(α)dS

Sα

Видно, что в условиях сделанного утверждения система (10.37) имеет только тривиальное решение: все dn(α)=0. При этом следует учесть, что матрицы систем (10.37) и (10.35) являются эрмитово-сопряженными, а область V

нерезонансна. Но из равенства dn(α)=0 следует ϕ=0 на S', т.к. системы {ϕn(α)}, n=1, 2, ... , для каждого α полны в соответствующих пространствах L2(Sα).

Для доказательства линейной независимости приравняем к нулю некоторое разложение по функциям Φmβ на S':

N ∞ ~

= 0 на S'.

(10.38)

∑ ∑Cm(β )Φmβ

β =1m=1

Далее прейдем от (10.38) к эквивалентной цепочке соотношений ортогональ-

~

ности вида (10.34). Тогда для Cm(β ) получим вновь систему (10.35), которая,

как показано выше, имеет лишь тривиальное решение.

Аналогично можно рассмотреть также уравнение (10.29) со смешанными граничными условиями на S (условие Дирихле на S0, условие Неймана на S'). В предложении, что p не является собственным значением данной

 

m

 

 

 

 

∂Φβ

(p)

 

краевой задачи, можно утверждать, что система функций

 

 

 

линейно

∂n

 

 

 

 

 

 

 

 

независима и полна в L2(S').

172

Рассмотрим теперь область ~ , ограниченную 0 и гладкими, но не

V

V

S

плоскими поверхностями Σα (рис. 3.1а). Пусть выполнены условия теорем о линейной независимости и полноте, доказанных выше, и Σα Sα= . Справедливо следствие, состоящее в том, что система функций {Φmβ (p)}полна в

L2(Σ'),

 

 

 

N

 

 

 

 

 

 

Σ′ = UΣα .

 

 

 

 

 

 

α =1

~

(Μ),

~

 

 

 

 

Это следствие доказывается следующим образом: пусть ϕ

Μ V , - ре-

шение уравнения (3.1), удовлетворяющее условию Дирихле на S0 и условию

~

(Μ)

 

Σ′ = f (Μ0 ), а f (Μ0 ) - заданная функция, Μ0 Σ′.

 

 

 

 

 

ϕ

 

 

 

 

 

 

~

f (Μ0 ) L2

(Σ′). Для лю-

 

 

 

Функция ϕ(Μ) существует для произвольной

бого набора Σα имеется S' такая, что соответствующий ей объем V нерезонансен. Аналогично по теореме Купрадзе можно показать, что существует

F(Q) L2(S'), для которой

 

 

 

~

′

′

,

(10.39)

ϕ

(Μ)∫F (Q)G (Μ, Q)dSQ ,

Q S

S ′

где G (M, Q) - функция Грина в области V, подчиненная указанным выше смешанным граничным условиям. Используя (3.11) и неравенство Буняков- ского-Шварца, получаем, что

f (Μ0 )−

где ˆm(β )

C

n0 N

ˆ

m

 

 

 

n0 N

ˆ

∂Φmβ

 

 

 

 

 

 

∑ ∑Cm(β )Φβ

(p)

 

≤ const

F(Q)− ∑ ∑ Cm(β )

 

 

 

 

(10.40)

 

∂n

 

 

n=1β =1

 

 

 

L2 (Σ′)

 

n=1β =1

 

 

L2

(Σ′)

 

 

 

 

 

 

 

 

 

 

- произвольные коэффициенты. Выбором n0 и подбором этих ко-

эффициентов правая часть (3.12) может быть сделана сколь угодной малой, откуда и вытекает утверждение следствия. Совершенно аналогично можно

доказать полноту {Φmβ (p)}в L2, на произвольной замкнутой поверхности внутри V.

Доказанное выше утверждение о полноте {Φmβ }позволяют использо-

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

173

10.6.2 Расчет физических параметров полосковых линий передачи

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

В квази -Т- приближении задача сводится к вычислению потенциала ϕ (x, y),

удовлетворяющего

уравнению

Лапласа

2 ϕ = 0

в

области

 

 

 

xy

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Рис. 10.12.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

 

< a , − b

≤ z ≤ b

 

 

Граничные условия имеют вид

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

2

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

± a

, z

 

 

 

 

 

 

 

(x,− b2 )=ϕ(x,b1 )= 0,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ϕ

 

 

 

 

=ϕ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ϕ(x,−(∆ + ∆

 

))

 

 

 

 

 

 

d 2 = 0, ϕ(x,0)

 

 

 

 

 

 

d1 =1.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

2

 

 

 

x

 

 

 

 

 

 

 

 

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

≤ 2

 

 

 

 

≤ 2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Кроме того, должны быть выполнены условия непрерывности Z-

d

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

компонента вектора электрического смещения D = −ε ϕ

при z=0

 

x

≥

1

,

 

2

 

 

 

 

 

 

 

 

 

 

 

d2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

z=–∆1, z=–(∆1+∆2)

 

 

x

 

≥

 

(при этом ε=1 при z>0,

z ≤ −(∆1 + ∆2 ), ε=ε1, при

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0 ≥ z ≥ −∆1, ε=ε2, при − ∆ ≥ z ≥ −(∆1 + ∆2 ).

 

 

 

 

 

 

 

 

 

 

 

 

 

Наибольший прикладной интерес представляет погонная емкость линии, определяемая соотношением

174

 

 

 

d1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

C =ε

 

 

2

 

 

 

 

ε

∂ϕ

 

 

 

 

 

−

∂ϕ

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

∫d

 

 

 

 

1 ∂z

 

 

 

Z =−0

∂z

 

Z

=+0

dx ,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

−

1

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где ε0 - диэлектрическая проницаемость вакуума.

 

Через нее выражается ε

 

=

 

C

 

, где C =

lim C

 

 

 

C

 

 

 

эфф

 

 

 

 

 

0

 

ε1, 2 →1

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

и волновое сопротивление Z

 

 

=120πε−1

(C

 

C )−

 

В

 

2 .

 

 

 

 

 

 

 

 

 

 

0

 

 

0

 

 

 

 

Конкретизируем для данной задачи способ построения базисных функ-

ций из п.10.6.1. Рассмотрим полосу

 

x

 

≤

α , − ∞ ≤ Z ≤ ∞ (рис. 10.12б) и по-

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

строим семейство решений краевых задач для уравнения Лапласа Um(X, Z), m=1, 2,... , при граничных условиях

 

a

,Z

 

 

= 0, Um (X ,0)= 0,

 

X

 

≤

d

,

 

 

 

 

 

 

 

Um ±

2

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(10.41)

 

∂Um

 

 

∂Um

 

 

 

 

 

 

 

 

 

 

 

Z =+0 =ε

 

Z =−0.

 

 

 

 

 

 

 

 

 

 

 

Потребуем, чтобы

 

∂Z

 

 

 

∂Z

 

 

 

 

 

 

 

 

 

 

lim [Um − cos(γ m x)eγ m Z ]= 0 ,

 

lim Um =0 и

 

 

Z →−∞

 

 

 

 

Z →−∞

 

 

 

 

 

 

 

 

 

 

 

где γ m = (2m +1)π . Функции Um(X, Z) можно представить в виде разложе-

a

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ний:

 

 

 

 

 

∞

 

 

 

 

Um (X ,Z )= cos(γ m x)eγ m Z

 

 

±γ p Z .

 

+ ∑C(pm) cos(γ p x)e

(10.42)

 

 

 

 

 

 

p =0

 

 

 

 

(знаки “+” и “-” соответствуют областям Z >0 и Z< 0), а для коэффициентов

C(pm) получить из граничных условий (10.41) системы парных сумматорных уравнений вида

 

 

 

 

 

∞

 

 

(pm)cos(γ p x)= −cos(γ m x),

x ≤ d ,

 

 

 

 

 

 

 

∑ C

 

 

 

 

 

∞

p =0

 

 

 

 

 

 

 

2

 

(10.43)

 

 

 

 

 

 

 

 

1 − ε cos(γ m x), x

≥ d .

 

 

∑ C(pm)γ p cos(γ p x)=γ m

 

 

p =0

 

 

 

 

1 + ε

 

 

 

2

 

Уравнения (10.43) могут быть точно решены методом задачи Римана-

Гильберта; в результате имеем

 

 

 

 

 

 

 

 

 

 

1

 

 

 

p

 

 

 

 

 

 

 

 

 

2 m +

 

 

(−1)

 

 

 

R

 

 

 

 

 

2

 

 

p

 

 

 

C(pm) =

 

 

 

 

 

 

Vpm −Vp−(m+1) −

 

(Vσm

−Vσ−(m+1))

−δ pm (−1)p , (10.44)

 

 

 

 

1

 

R

 

 

 

 

 

 

 

 

 

 

 

 

 

(ε +1) p +

 

 

 

 

σ

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

175

где V m , V m , R

 

, R

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

π(d − a)

p

- специальные функции аргумента U = cos

,

p

σ

 

σ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a

выражаемые через полиномы Лежандра:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(u)P (u)], p ≠ m;

 

 

 

 

V m =

 

p +1

 

 

 

[P (u)P

(u)− P

 

 

 

 

 

 

2

(p − m)

 

 

 

 

 

 

p

 

 

 

 

p

 

m+1

 

 

 

p +1

 

 

 

m

 

 

 

 

 

 

 

 

Vmm (u)=

1

m+1

 

 

 

(u)Pk −m−1(u), m ≥ 0,

 

 

 

 

 

 

 

 

∑ρm+1−k

 

где ρ0 (u)=1,

ρ1(u)= −u,

 

 

 

 

 

2 k =0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ρn (u)= Pn (u)− 2uPn−1(u)+ Pn−2 (u),

n ≥ 2 ;

 

 

 

 

 

 

 

R

p

= 1 P(u), R

 

= π

P

 

(u);

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

σ

2

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

π

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

V m = −

 

 

 

 

 

 

 

 

P

1

(u)P

 

(u)

− P

(u)P (u) .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

σ

 

 

 

 

 

1

−

 

m+1

 

 

 

 

 

1

 

m

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4 m

+

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Кроме того, введем семейство функций Vm(X, Z) удовлетворяющих уравнению Лапласа, граничным условиям (10.41), но с условиями на бесконечности вида

lim Vm = 0 и

lim [Vm − cos(γ m x)e−γ m Z ]= 0 .

 

Z →∞

Z →∞

 

 

Функции Vm(X, Z) даются разложениями

 

 

 

∞

(pm)cos(γ p z)emγ p Z

 

Vm (X ,Z )= cos(γ m z)e−γ m Z + ∑B

(10.45)

p =0

а коэффициенты B(pm) вычисляются аналогично и даются выражением

B(pm) =ε C(pm) −δmp (1 −ε). Функции (10.42) и (10.44) будут использованы для

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

Для построения алгоритма решения разобьем исследуемую область на две частичные подобласти:

область 1:

 

 

 

x

 

≤ a

, − ∆ ≤ z ≤ b ;

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

область 2:

 

x

 

≤ a ,

− b ≤ z ≤ −∆

1

.

 

 

 

 

 

 

 

 

2

2

 

 

 

 

 

 

 

 

Будем искать решение в виде (область 1)

 

 

ϕ(X ,Z )=V (X , Z )+ ∑∞ [αm(1)Um(1)(X , Z )+ βm(1)Vm(1)(X , Z )]

m=0

 

 

 

и (область 2)

 

 

 

ϕ(X ,Z )= ∑∞ [αm(2)Um(2)(a − x,−(z + ∆1 + ∆2 ))+

m=0

+ βm(2)Vm(2)(a − x,−(z + ∆1 + ∆2 ))];

176

αm(j ), βm(j ) - неизвестные коэффициенты; Um(j ), Vm(j ) - функции (10.42), (10.44)

при d=d1, ε=ε1 для j=1 и d=d2, ε=ε2 для j=2; V(X, Z) - решение уравнения Лапласа, убывающее при z → ± ∞, удовлетворяющее условию

V (X ,0)=1,

 

X

 

≤

d1

, а в остальном - граничным условиям (10.41). Оно имеет

 

 

 

вид

 

 

 

2

V (X ,Z )= ∑∞ Ap cos (γ p x)e±γ p Z =

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

p =0

 

 

 

 

 

 

∞

4(−1)p R

 

 

 

 

 

 

= ∑

p

cos(γ p x)emγ p Z .

 

 

 

 

 

 

(2 p +1)R

 

 

 

 

 

 

p =0

σ

Для нахождения коэффициентов αm(1), βm(1) следует наложить условия

непрерывности ϕ и z-компоненты D при z= –∆1, а также граничные условия для ϕ при z=b1, z= –b2. Тогда получим систему линейных алгебраических уравнений вида

 

ξ11

ξ12

0

0

 

α(1)

 

ν1

 

 

 

 

 

 

 

 

ξ 23

 

 

 

 

 

0

0

ξ 24

 

β (1)

 

0

 

 

 

−ξ33

 

=

.

(10.46)

 

ξ31

ξ32

ξ34

α(2)

ν2

 

ξ31 ξ32 −

2I ξ33 ε2 ξ 44

 

β (2)

 

ν2

 

 

 

 

 

ε1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где α(j), β(j) - бесконечномерные векторы коэффициентов αn(j),

βn(j), a ξij -

матрицы следующей структуры:

 

 

 

 

 

 

 

ξmn11 =δmn + Cn(m) e−2γ nb1 ,

ξmn12 = (δmn + Bn(m) )e−2γ nb1 ,

ξmn31 = (δmn + Cn(m) )e−2γ n ∆1 ,

ξmn32 = δmn + Bn(m) e−2γ n ∆1 ,

ξ33 = −δ + C(nm) e−γ n (∆1 +∆2 ),

mn mn

 

 

 

 

 

 

 

 

 

ξmn34 = {−δmn e−γ n (∆2 −∆1 ) −

 

 

n(m) e−γ n (∆1 +∆2 )},

 

 

 

 

 

 

 

B

 

 

 

 

 

 

 

44

 

 

ε

2

 

−γ

 

(∆

 

−∆

 

)

 

 

 

 

(m)

e

−γ

 

(∆

 

+∆

 

)

,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ξmn =

 

−δmn e

 

n

 

2

 

1

 

+ Bn

 

n

 

1

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

ε1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Imn =δmn , ξ 23 , ξ 24

получаются из ξ11, ξ12 заменами

 

 

 

 

 

 

 

(m)

,

 

 

(m) →

 

 

(nm),

 

(nm)

и b1→ b1 - ∆1 - ∆2 (одна либо две черты над

 

 

 

 

 

 

 

 

 

 

 

C

B

C

 

B

 

n

 

 

 

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

177

Cn(m), Bn(m) означают их принадлежность, соответственно, первой и второй

частичным областям). Коэффициенты ν1,2 в правой части (10.46) даются равенствами

ν1=–An e–2γnb1, ν2=–An e–2γn∆1,

Система (10.46) является Фредгольмовой, что следует из асимптотических оценок Cn(m), Bn(m), An. Для нее строго обоснован метод редукции.

В качестве простейшего примера рассмотрим экранированную микрополосковую линию (рис. 10.13а).

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

например: ε1→ ε2, d2→ a, либо ε1→ ε2, ∆1+∆2→ b2, либо ∆2→ 0, d2→ a (рис. 10.13б, в, г).

На рис.10.14 приведены результаты расчета погонной емкости микро-

Рис. 10.13.

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

линия 1:

b2

= 0.25;

b1 + b2

= 0.5 ; ε=10.5;

a

a

 

 

 

178

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Рис. 10.14.

линия 2:

 

b2

= 0.1;

 

b1 + b2

= 0.9

;

ε=10.5;

 

a

a

 

 

 

 

 

 

 

 

линия 3:

b2

 

= 0.25;

 

b1 + b2

 

= 0.5

;

ε=1.0;

 

 

a

 

 

a

 

 

 

 

 

 

линия 4:

b2

 

= 0.1;

b1 + b2

 

= 0.

ε=1.0.

 

 

 

 

a

 

 

 

a

 

 

 

В табл. 10.1 приведены примеры характеристик полосковых линий, часто применяемых на практике

179

 

 

 

 

 

 

 

 

Табл.10.1

 

 

 

 

 

 

 

 

 

 

ε

а,

d1,

b1,

h,

N

I

εэфф

ZВ,

 

 

мм

мм

мм

мм

4

 

 

Ом

 

9.8

1

0.02

20

1

2

300

5.403383

136.3948

 

9.8

1.5

0.05

20

1

3

300

5.430295

118.5018

 

9.8

2

0.1

20

1

5

300

5.492849

103.9079

 

9.8

8

0.6

20

1

5

300

6.229931

64.07808

 

9.8

10

1

20

1

5

300

6.484238

50.78413

 

9.8

15

3

20

1

10

300

7.28545

26.13688

 

9.8

20

8

20

1

10

300

8.107649

12.14053

 

3.78

2

0.1

20

1

5

300

2.419478

156.5620

 

3.78

8

0.6

20

1

5

300

2.662139

98.02476

 

3.78

10

1

20

1

5

300

2.744443

78.06032

 

3.78

18

5

20

1

10

300

3.136074

28.01408

 

1

20

8.87

0.5

0.5

10

500

1

10.18

 

1

25

7.31

0.5

0.5

10

500

1

12.27

 

1

25

6.2

0.5

0.5

10

500

1

14.33

 

1

23

5.37

0.5

0.5

10

500

1

16.40

 

1

20

4.72

0.5

0.5

10

500

1

18.46

 

1

15

3.13

0.5

0.5

10

500

1

26.73

 

1

5

1.21

0.5

0.5

10

500

1

57.98

 

1

3

0.19

0.5

0.5

10

500

1

162.02

 

180