книги из ГПНТБ / Амитей Н. Теория и анализ фазированных антенных решеток
.pdfОсновные формулировки граничной задачи |
81 |
Другой важный вид преобразований уравнения |
(44) связан |
с переводом этого интегрального уравнения в дискретную форму в смысле гильбертова пространства [6]. Например, неизвестную функцию Е( — точное решение уравнения (44) — можно опреде лить как континуум значений, заданных на совокупности точек г
вобласти А'. Если известна дискретная, но бесконечная система значений vt и A t -f- R t в выражениях (26) и (28), то функция Ег полностью определена. Полезность интегрального уравнения состоит даже не в том, что оно записано в координатах г, а скорее
втом, что оно содержит в компактной форме все граничные условия п может служить отправной точкой для получения различных эквивалентных представлений. Эквивалентность можно понимать,
вчастности, в том смысле, что если г — непрерывная координата
и |
Е, (г) — искомая |
функция, то |
Фг — дискретная |
координата, |
а |
Е ( (Фг) = Vi (i = |
1, 2, . . ., |
оо) — дискретный |
эквивалент |
Е,(г).
Преобразование уравнений (44), (55) и (93) в эквивалентную дискретную форму оказывается в особенности полезным, когда используется конечная дискретная аппроксимация по методу Ритца или по методу Галеркпна [21] (см. также гл. 3).
2.6. Интегральное уравнение для случая возбуждения одного излучателя
Напомним, что если удалось определить путем расчета, [выра жение (59) и рис. 2.5], либо экспериментально коэффициенты вза имной связи элементов антенной решетки с одним возбужденным элементом, то решение задачи для ФАР, в которой возбуждены все излучатели, можно найти методом суперпозиции. Подробно этот вопрос освещен в гл. 8. В данном разделе нас интересует то, что в задаче с единственным возбужденным излучателем можно составить интегральное уравнение относительно тангенциальной компоненты электрического поля. Более того, для расчета коэффи циентов взаимной связи можно применить вариационный принцип.
Вывод интегрального уравнения и вариационного выражения при возбуждении в антенной решетке только одного элемента подобен выводу для случая возбуждения всех элементов. Для простоты рассмотрим линейную антенную решетку из плоско параллельных волноводов с толстыми стенками (рис. 2.9) [24]. Результаты можно обобщить на плоские решетки с элементами других типов [24, 25]. Отметим, что и в этом случае, так же как и при выводе вариационного выражения (85), должны быть выпол нены условия симметрии элементов антенной решетки.
Рассмотрим случай, когда только один элемент с индексом О антенной решетки (рис. 2.9) возбуждается падающей волной, амплитуда которой равна 1. Поле внутри каждого волновода
0 - 0 1 6 8
82 |
Глава 2 |
можно разложить по обычным волноводным гармоникам, а поле излучения в свободном пространстве можно связать с полем в раскрывах волноводов с помощью преобразования Фурье. Пусть функция
ф n(x — pb), |
I X— pb I < - | , |
Фрп (х) = |
(97) |
О, |
— pb\, |
представляет собой п-ю ортонормпрованную гармонику в р-ом волноводе, а Ф„ (х) — п-ю гармонику в нулевом волноводе. Пред полагается, что волноводы имеют ширину а и их центры находятся
Рпс. 2.9. Схема линейной антенной решетки из плоскопараллельпых волноводов с толстыми стенками.
на расстоянип Ъ один от другого. Пусть Ст обозначает коэффи циент овязн по току на п-й гармонике возбуждаемого волповода и волновода с индексом р. Напряженность магнитного поля в раскрыве решетки (—оо < х < оо) можно представить в виде
ОО
Я я (*)= 2 [ (^рс — Ср1) ФР 1 (х) +
ОО
со |
|
pb+a/2 |
+ 2 |
ФРп(а:) j Фрп(х') нх(х') dx'^ (98) |
|
п = |
2 |
р Ь -а /2 |
при z = 0" и
Н х (х) = 2Г J dx' J d k J k^ - x,)H x {x') |
(99) |
|
Основные формулировки |
граничной |
задачи |
83 |
при z — 0+. |
Тогда электрическое |
поле в |
раскрыве |
решетки |
(— go z |
оо) имеет вид |
|
|
|
оо
Еу(х) — 2 [( — ро— Ср1) ZjФр1(х)-(-
Р — — ОО
ооp b - j- a /2
+ 2 2„ФР„ (*) |
J |
фрп (*') нх (х') dx'] |
(100) |
|
n = 2 |
|
p b — a / 2 |
|
|
при z = 0~ и |
|
|
|
|
°0 |
OO |
jfo |
|
|
- H r J a*' J ‘“>(Уг1Гя“м |
(101) |
|||
— |
OO |
— OO |
'X |
|
при z = 0+, где волновые |
сопротивления Zn определяются по |
|||
формуле |
|
|
|
|
2_______ top
У/ с 2 — ( я я /a )2
Условия непрерывности электромагнитного поля в раскрыве
[Нх {х, z = 0~) = Hx (x, z = 0+) и Еу (х, z = 0~) = Еу (х, z = 0+)]
приводят к интегральным уравнениям Фредгольма первого и вто рого рода относительно тангенциальной составляющей магнит ного поля в раскрыве
ООоо
2Z^t>0i (х) = 2 |
2 2 пФрп (*) 1рп + |
Р = —ОО 7 1 = 1
|
ОО |
|
+ ^ |
[ H f { k \ x - x ' \ ) H x {x')dx' |
(102) |
|
— оо |
|
и
ОООО
НХ(Х)= 2Ф01( * ) - 2 |
2 ^ |
^ |
фрп(х)1рп- |
|
р = — ОО |
7 1 = 2 |
|
|
|
|
|
оо |
|
|
- |
^ |
j Н ^( |
(к \ х - х' I) Нх (х') dx', |
(103) |
где |
рЪ+а/2 |
|
|
|
|
|
|
|
|
Ipn= |
j |
Фрп(х') Hx (x')dx', |
(103а) |
|
рЬ—а/ 2
аЯо2) (ж) — функция Гаикеля второго рода нулевого порядка. Уравнение (103) можно преобразовать к уравнению для собст
6*
84 |
Глава 2 |
венных значений:
а /2
J Ф01{х')Нх (х')йх' = —af2
|
со |
|
|
а/2 |
|
= |
2 |
ZiCDpi (я) |
^ |
Фр1 (я ) Нх {%) dxl -)- |
|
|
Р — — оо |
|
— а / 2 |
|
|
|
со |
со |
|
а/2 |
|
+ |
2 |
2 z "clv |
(*) |
J ф „п(*') я* (ж') d®' + |
|
р ~ —СО 7 1 = 2 |
|
-а /2 |
|||
|
|
|
|
СО |
|
|
|
+ |
i | i |
{ |
п Ч \ к \ х - х ' \ ) Н х {х')йх', (104) |
где 2 ? означает, |
что опущено |
слагаемое с индексом р — 0, а |
|||
Zx (1 + С01)/(1 — |
С01) представляет собой нормированное входное |
||||
сопротивление возбужденного элемента решетки. Поскольку ядро уравнения (104) является комплексно-симметричным, для входного сопротивления, а следовательно, и для С01 можно написать ва риационное выражение
|
|
2 2 |
Z A + |
|
р = —оо |
р = —оо п =2 |
|
+ |
оо |
оо |
(к \ хJ- х ' I) IIх (х') dx' dx, (105) |
|
Нх (®) IIfJ |
||
в котором 1рп задано формулой (103а). |
|||
Чтобы получить вариационное |
выражение для Cpi при р Ф 0, |
||
рассмотрим возбуждение двух элементов антенной решетки — эле
ментов с индексами 0 и q. Пусть Еу (х’) и II^ (х') обозначают распределения электрического и магнитного полей в раекрыве
решетки. |
Тогда |
|
|
|
|
|
ОО |
|
|
|
|
Н х ( х) = |
2 |
— R p) Ф Р1(X) + |
|
|
|
р = —ОО |
|
|
|
|
|
|
|
ОО |
ОО |
оо |
|
|
+ |
2 2 |
ф* п И |
J ф Рп(*,) я ^ (® ')^ ' |
(106) |
р г = — ОО 7 1 = 2 |
— ОО |
Основные формулировки граничной задачи |
85 |
п
СО
К ( * ) = 2 |
( - 6 р о - б л 9 - Л р ) Ф р 1 ( * ) + |
|||
Р = — СО |
|
|
|
|
|
СО |
СО |
|
оо |
+ |
2 |
2 ^пФрп(х) |
j <bPn { x ' ) H l { x ' ) d x ' = |
|
|
р——оо п = 2 |
|
—оо |
|
|
|
|
оо |
|
|
|
= |
j |
Я (о2) ( k \ x - x ' \ ) H i (я') dx', (107) |
где R p= Ср + Cp_q. [Если в обозначении Срп опущеи второй индекс, это означает, что коэффициент относится к волне основного типа, т. е. Ср‘п = Ср[.] Комбинируя выражения (106) и (107), получаем уравнение
p b - fa / 2
C T ^ A ,i ( * ) |
j |
Ф01 (* ')# £ (*')<**'+ |
|
|
||||
|
p b — a / 2 |
|
|
|
|
|
|
|
|
1 +Rq |
|
p b + a / 2 |
|
|
|
||
|
ZiOql(x) |
j |
Oqi(x') Hx (x’) dx' = |
|
||||
|
i - R q |
|
p b — a / 2 |
|
|
|
||
|
|
|
|
|
|
|
||
|
° o ( 0 , q) |
|
p b + a / 2 |
|
|
|
||
= |
2 |
z $>pl{x) |
j |
Фpl{x')Hl{x')dx' + |
|
|||
|
p = - o o |
|
p b — a/2 |
|
|
|
||
|
OO |
OO |
|
pb-j-a /2 |
|
|
|
|
+ |
2 |
2 |
2 »ф*»(*) |
J |
ф р п ( * ') ^ ( я ') ^ ' + |
|
||
p = — oo n = 2 |
p b —a / 2 |
|
|
|
||||
|
|
|
|
OO |
|
|
|
|
|
|
|
+ - ^ |
j |
H ^ ^ x - x ' ^ H l i x ^ d x ' , |
(108) |
||
oo (0, 3) |
|
|
|
|
|
|
|
|
в котором 2 |
|
означает исключение из суммы слагаемых с индек- |
||||||
р = —со |
|
|
|
|
|
|
|
|
сами 0 и q. |
|
|
R g = С0 + Cg, левую часть уравнения (108) |
|||||
Поскольку i?0 = |
||||||||
можно записать в виде произведения суммы ср01/^ + |
[где |
|||||||
коэффициент |
1^п определяется как 1рп, |
если вместо Н х исполь |
||||||
зовать Н % в |
выражении |
(103а)] |
на |
выражение Zx (1 + |
С0 + |
|||
+ Cq)/(1 — С0— Cq), представляющее собой входное сопротивление одного из возбуждаемых элементов (нулевого илид-го). Вариацион ное выражение получить легко, так как ядро уравнения (108)
86 |
|
|
Глава 2 |
|
комплексно-симметричное. Имеем |
|
|||
Z i \ t c l t c l ( |
^ |
+ ^ |
2)== 2 2 1[/J1]2 + |
|
|
|
|
р — — СО |
|
|
оо |
со |
оо |
оо |
+ |
2 |
2 |
г п и 1 п ? + ^ j |
J л хт(х)х |
|
Р=—со п = 2 |
—оо —оо |
||
X Яо2) (Л I ® —я' I) //J (х') dx' dx. (109)
Выражение (109) можно упростить, если использовать симметрию антенной решеткн относительно излучателя, который располо
жен |
посередине между |
элементами с индексами |
р — 0 и р = д. |
||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
Л И Т Е Р А Т У Р А |
|
|
|
|
|
|
|
|
|
||||
1. |
Collin R. Е. F ield Theory of G uided W a v es, M cG raw -H ill, N ew Y ork, 1960. |
||||||||||||||||||||||
2. Amitay N., |
Galindo |
Г. The A n alysis of Circular W aveguide Phased |
Arrays, |
||||||||||||||||||||
|
«B ell System |
T ech. |
J.», |
1968, |
v . 47, |
N o. |
9. |
|
|
|
|
|
|
|
|
||||||||
3. |
Marcuvitz N. (ed.). |
|
W aveguide |
H andbook, |
M IT |
R ad iation |
Laboratory |
||||||||||||||||
|
Series, v . 10, |
M cG raw -H ill, N ew |
Y ork, 1951. |
|
|
|
|
|
|
|
|
||||||||||||
4 . |
Titchmarsh E. C. The Theory |
of |
F unctions, |
O xford |
U n iv ersity |
Press, Lon |
|||||||||||||||||
|
don, Paragraphs 1.7, 13.5, and 13.53 |
(p. 419 — 421), |
1952. |
|
|
|
|
||||||||||||||||
5. |
Sommerfeld A. J. O ptics, |
A cadem ic Press, |
N ew Y ork, 1954. |
|
|
|
|||||||||||||||||
6. |
Friedman B. |
P rinciples |
and |
T echniques |
of |
A pplied |
M athem atics, |
John |
|||||||||||||||
|
W iley and |
Sons, N ew |
Y ork, |
1960. |
|
|
|
|
|
|
|
|
|
|
|
||||||||
7. DuFort |
E. |
C. |
F in ite |
Scattering |
M atrix |
for |
an |
In fin ite |
A ntenna |
Array, |
|||||||||||||
|
«R adio Science», 1967, v . 2, N o. 1, p. 19— 27. |
M. |
|
|
|
|
|
|
|
||||||||||||||
8. |
Montgomery |
C. |
G., |
Dicke |
R. |
II., |
Purcell |
E. |
P rinciples of |
M icrowave |
|||||||||||||
|
C ircuits, |
R ad iation |
Laboratory Series, v . 8, M cG raw -H ill, N ew Y ork, 1948. |
||||||||||||||||||||
9. |
Amitay N., |
Galindo |
V. A |
N ote on |
the |
R ad iation |
C haracteristics and |
Forced |
|||||||||||||||
|
Surface |
W ave |
Phenom ena |
in |
T riangular Grid |
Circular |
W aveguide |
Phased |
|||||||||||||||
|
Arrays, «IE E E |
Trans. A ntennas and Propagation», 1968, v . A P -16, p. 760— |
|||||||||||||||||||||
|
762. |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
10. |
Hannan |
P. |
W., Meier P. J., |
Balfour M. A. S im u lation |
of Phased |
Array |
|||||||||||||||||
|
A ntenna |
Im pedance in W aveguide, «IE E E |
|
Trans. A ntennas and |
Propagati |
||||||||||||||||||
|
on», 1963, v . |
A P -11, |
N o. |
6. |
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
11. |
Balfour |
M. |
A. |
Phased |
Array |
Sim ulators |
in |
W aveguide |
for a |
Triangular |
|||||||||||||
|
A rrangem ent of E lem ents, «IE E E |
Trans. A ntennas and Propagation», 1965, |
|||||||||||||||||||||
|
v . A P -13, N o. 3, p. |
475 — 476. |
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||
12. |
Hannan |
P. |
W., Balfour M. A. S im u lation |
of |
a |
Phased-A rray |
A ntenna in |
||||||||||||||||
|
W aveguide, |
«IE E E |
Trans. A ntennas |
and |
Propagation», |
1965, |
v . |
A P -13, |
|||||||||||||||
|
N o. 3, p. 342— 353. |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
13. |
Balfour |
M. |
A. |
A ctive |
Im pedance of a |
Phased |
Array |
A ntenna |
E lem ent |
||||||||||||||
|
Sim ulated b y a Sin gle |
E lem ent in W aveguide, «IE E E Trans. A ntennas and |
|||||||||||||||||||||
|
Propagation», |
1967, |
v . A P -15, |
N o. 2, |
p. 313 — 314. |
|
|
|
|
|
|
||||||||||||
14. |
Hildebrand |
F. B. M ethods of |
A p plied |
M athem atics, |
P ren tice-H all, |
E ngle |
|||||||||||||||||
|
wood C liffs, |
N .J ., 1952. |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
15. |
Amitay |
N., |
Galindo |
|
F ., |
|
Wu |
С. |
P. |
On |
Certain |
E igen valu e P roblem s in |
|||||||||||
|
E .M . W ave |
Scattering, |
private |
m em oranda, |
ava ila b le |
|
on request. |
|
|||||||||||||||
16. |
Stratton |
J. |
A. |
E lectrom agnetic |
T heory, |
M cG raw -H ill, |
N ew York, |
1941, |
|||||||||||||||
|
p. |
486 — 487. |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
Основные формулировки, граничной |
задачи |
|
|
|
|
87 |
|||||||||||||||||
17. Galindo V., Wu С. Р. |
Integral |
E quations |
and V ariation al |
E xpressions for |
|||||||||||||||||||||||
|
A rbitrary |
Scanning |
of |
R egular |
In fin ite |
Arrays, |
«IE E E |
Trans. |
A ntennas |
||||||||||||||||||
|
and Propagation», |
1966, |
v. A P -14, |
N o. 3, p. |
392 — 394. |
|
|
|
|
|
|||||||||||||||||
18. |
Collin R. E. A |
N ote on |
W aveguide |
Im ago |
T echniques, |
Case |
In stitu te of |
||||||||||||||||||||
|
T echnology |
C leveland, |
O hio, Sci. R ept. 19, |
A D -250931, |
1960. |
|
|
|
|||||||||||||||||||
19. |
Saxon D. |
S. |
N otes |
on |
Lectures |
b y |
Julian |
Schw inger — |
D isco n tin u ities |
||||||||||||||||||
|
in W aveguides, to |
be published |
shortly; |
see |
|
also |
R ef. |
3. |
|
|
|
|
|
||||||||||||||
20. |
Galindo V. On |
the |
S tation ary |
Properties |
of |
|
the |
Integral |
E quations |
of Infi |
|||||||||||||||||
|
nite |
Phased |
Arrays, |
«IE E E |
Trans. |
[Antennas |
and |
Propagation», |
1969, |
||||||||||||||||||
|
v. A P -17, |
N o. |
3. |
Krylov V. I. A pproxim ate |
|
|
|
|
|
|
|
|
|
|
|||||||||||||
21. |
Kantorovich L. |
V., |
M ethods of H igher |
A n alysis, |
|||||||||||||||||||||||
|
Interscience |
P ublishers, |
|
N ew |
Y ork, 1964. |
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||
22. |
Jones D. S. A Critique |
of the |
V ariation al |
M ethod |
|
in |
Scattering |
Problem s, |
|||||||||||||||||||
|
«IR E |
Trans. A ntennas |
and Propagation», |
1956, |
v . A P -4, p. |
297— 301. |
|||||||||||||||||||||
23 . Amitay N., Galindo |
V. |
A p p lication |
of |
a |
N ew |
M ethod |
for |
A pproxim ate |
|||||||||||||||||||
|
Solu tion s |
and |
Error |
E stim ates |
to |
W aveguide |
D isco n tin u ity |
and |
Phased |
||||||||||||||||||
|
Array Problem s, |
«Radio |
Science», |
1968, |
v . 3, |
N o. |
8. p. |
830— 843. |
|
||||||||||||||||||
24 . Galindo V., Wu |
С. |
P. |
A V ariation al |
E xpression |
for |
the D om inant M ode |
|||||||||||||||||||||
|
C oupling |
C oefficients B etw een the |
E lem ents |
in |
an |
In fin ite |
Array, |
«IE E E |
|||||||||||||||||||
|
Trans. A ntennas |
and |
Propagation», |
1966, v . A P -14, |
N o. 5, |
p. |
637— 639. |
||||||||||||||||||||
25. |
Galindo V., Wu С. P. |
N um erical |
S olu tion |
|
for |
an |
|
In fin ite Phased |
Array |
||||||||||||||||||
|
of R ectangular W aveguides w ith |
T hick W alls, |
«IE E E |
Trans. A ntennas and |
|||||||||||||||||||||||
|
Propagation», |
1966, |
v . A P -14, p. 149— 158. |
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||
26. |
Wu С. P. N ote on |
Integral E quations and V ariation al E xpression for A rbi |
|||||||||||||||||||||||||
|
trary |
Scanning |
of |
R egular In fin ite Arrays, |
«IE E E |
Trans. |
A ntennas and |
||||||||||||||||||||
|
Propagation», |
1968, |
v. A P -16, |
p. 136— 138. |
|
|
|
|
|
|
|
|
|
|
|
||||||||||||
3. Методы решения
1.ВВЕДЕНИЕ
Вгл. 3 выведены интегральные уравнения для граничных за дач, относящихся к бесконечным волноводным решеткам. Особен ностью этих уравнений является то, что ядра имеют внд бесконеч ных сумм, представляющих собой вклады типов волн в двух суще ственно различающихся «волноводах». Из-за сложности задачи такие интегральные уравнения не могут быть решены аналити
чески, за исключением одного-двух частных случаев. Поэтому в этой главе основное внимание уделено численным методам решения.
Для применения численных методов необходимо сначала с помощью метода моментов перейти от интегрального уравненпя к системе линейных алгебраических уравнений. После этого мож но воспользоваться известными методами решения такой системы, среди которых наиболее распространенным является метод обра щения матрицы. Хотя обращение матрицы легко выполняется на электронных вычислительных машинах, эффективность решения зависит от ряда факторов, таких, как простота вывода системы линейных алгебраических уравнений н вычисления элементов матрицы. Известен ряд хороших функциональных представлений для искомого решения. В частности, во многих случаях базис волноводных типов волн дает ряд преимуществ по сравнению с другими базисамп.
Систему линейных алгебраических уравнений можно получать также путем применения процедуры Релея — Ритца к вариацион ному выраяшнию для входного импеданса (или проводимости) решетки. Условия справедливости вариационного принципа об суждены в гл. 2. В данной главе показано, что применение метода моментов при решении интегрального уравнения п процедура Релея — Ритца приводят к одной и той же системе линейных алгебраических уравнений. Предполагается, что решения, полу ченные методом моментов, автоматически удовлетворяют соотно шению взаимности независимо от их точности. Это имеет ваяшое значение для проверки верности решений.
Другой подход к решению интегральных уравнений основан на приближенной замене ядра вырожденным ядром. Из теории интег ральных уравнений хорошо известно, что интегральное уравнение
Методы решения |
89 |
с вырожденным ядром можно свести к конечной системе линейных алгебраических уравнений. Можно показать, что при соответ ствующем выборе базиса приближенного решения системы линей ных алгебраических уравнений, полученные этим методом и мето дом моментов, совпадают. Следовательно, замена ядра вырожден ным ядром эквивалентна применению метода моментов.
Важным аспектом применения численных методов является вопрос о корректности и точности решения. Некоторые способы проверки верности решения подробно описаны в этой главе.
Для частного впда решетки из тонких параллельных пластин можно найти точное решение уравнений методом Винера — Хопфа или методом вычетов. Это решение, полезное во многих отношени ях, подробно рассмотрено в разд. 6. Поскольку уравнения, исполь зуемые в методе вычетов, выводятся из интегральных уравнений с помощью соответствующего выбора базиса, метод вычетов рас смотрен после метода моментов.
Вместо непосредственного обращения матрицы можно искать решение системы линейных алгебраических уравнений построе нием ряда Неймана, который быстро сходится во многих случаях. Этот метод также полезен для определения более точных границ погрешности решения.
2. ПРЕДСТАВЛЕНИЕ РЕШЕНИЯ С ПОМОЩЬЮ ДИСКРЕТНОГО БАЗИСА
Граничные задачи, возникающие в теории антенных решеток, можно сформулировать в впде интегральных уравнений Фредголь ма первого пли второго рода. Поскольку интегральные уравнения первого рода более удобны для решения этих задач, мы рассмот рим главным образом эти интегральные уравнения. Интегральное уравнение Фредгольма первого рода в операторной форме имеет вид
у = Кх, |
(1) |
где К — линейный оператор, определяемый только геометрией решетки, у — известная функция, определяемая полем возбуждающего сигнала, х — неизвестная функция. Обычно к и у — векторные величины, а К — тензор *).
Область определения оператора, т. е. систему функций, на кото рую данный оператор действует, будем называть пространством
х) При постановке математической задачи возникают два вопроса: суще ствует ли решенпе и, если существует, то является ли оно единственным. Доказательство существования решения уравнения (1) довольно сложно, и оно выходит за рамки данной книги. Однако, имея дело с физической зада чей, можно считать, что решение существует, если задача верно поставлена. Что же касается единственности .решения, то следует помнить, что исходными уравнениями для таких задач являются уравнения Максвелла, для которых доказана единственность решения.
90 |
Глава 3 |
S D. (Так как мы имеем дело с интегральными операторами, то граничные условия учитываются автоматически; в случае диффе ренциальных операторов при задании области определения необ ходимо учесть граничные условия.) Если оператор действует иа функцию, он порождает новую функцию. Набор всех функций, порождаемых оператором при его действии иа любой элемент в области определения оператора, будем называть простран ством S L. Свойство лнпейности оператора означает, что еслп даны две функции хг и .г\, ц два произвольных числа сц и а 2, то выпол няется соотношение
A (CL^X1-|- 0^2^) == |
-{- OL^K‘-3'2' |
Предположим теперь, что мы нашли каким-то образом систему линейно независимых функций, перекрывающую *) пространство S Dн другую систему, перекрывающую пространство S ь. Запишем эти две системы функций в виде {fn: /„ 6 S в } и {gn: gn £ FL) соответственно. Можно считать, что функции этих двух систем ортонормпровапы. Это не является ограничением, так как любую линейно независимую систему функций можно ортонормировать с помощью метода Грама — Шмидта.
Так как неизвестная функция х является элементом простран ства S д, то ее можно представить как линейную комбинацию функ ций /„:
я = X ап1п. |
(2) |
где ап — неизвестные коэффициенты. Такое представление обычно содержит бесконечное число членов, если требуется найти точное решение уравнения (1). В этом случае системы функций {/„} л {gn} должны быть полными в соответствующих пространствах. Для нахождения приближенных решений суммы (2) обычно берут ся конечными. Система функций {/„} необязательно должна быть полной, хотя для получения верных результатов важно включать определенные элементы пространства. Это условие более подробно рассмотрено ниже. Сначала рассмотрим случай точного решения.
Подставляя выражение (2) и уравнение (1) найдем
у = 2 «»*/»• |
(3) |
Разность между левой и правой частями выражения (3) является нулевым элементом в пространстве S L. Чтобы убедиться в этом,
потребуем, чтобы элемент (у — У anKfn) был ортогонален ко всем элементам системы, перекрывающей S L:
_________ |
{gm, y — '2<ZnKfn) = 0 для всех т, |
(4а) |
г) Говорят, что система функций «перекрывает» пространство, еслп лю бую функцию из этого пространства можно, согласно соответствующей мере пространства, представить как линейную комбинацию функций системы.
