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

Учебное пособие 1920

.pdf
Скачиваний:
14
Добавлен:
30.04.2022
Размер:
2.9 Mб
Скачать

нат будет в одном из его концов; длина провода будет l.

При этом сила тока i и напряжение v в любой точке провода есть функции координаты х и времени t. Ток i и напряжение v описываются дифференциальными уравнениями с частными производными. В этих уравнениях полагается, что емкость, сопротивление, самоиндукция и утечка вдоль провода определены непрерывно и равномерно и постоянные С, R, L и G рассчитаны на единицу длины провода. Пусть часть провода между двумя сечениями х=х1 и х=х2. По законуОма для этой части провода:

 

 

 

 

 

x2

x2

дi(x,t)

 

v(x1,t) v(x2,t) R i(x,t)dx L

dx,

 

 

 

 

 

 

 

x

x

 

дx

 

 

 

 

 

1

 

1

 

 

 

v(x ,t) v(x

2

,t) x2

дv(x,t)

dx,

 

1

 

 

 

 

x

дx

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

x2

 

дv

 

 

дi

 

 

 

 

 

 

(

L

Ri)dx 0.

 

 

 

 

 

 

 

 

x

 

дx

 

дt

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

В силу произвольности координат х1 и х2 имеем

дv L дi Ri 0.

дx дt

Количество электричества, протекающего по участку

x2 дi

(х12) провода за единицу времени i(x1t) i(x2,t) дx dx,

x1

равно сумме количества электричества, которое нужно для зарядки этого участка провода, и количества электричества, которое теряется в связи с несовершенством изоляции:

 

 

 

 

x

дv

 

x

 

 

 

 

 

 

 

 

 

C 2

 

dx G 2

vdx ,

 

 

 

 

 

 

дt

 

 

 

 

 

 

 

 

x1

 

x1

 

 

 

 

 

или

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x2

дi

 

дv

 

 

 

 

 

дi

 

дv

 

(

C

Gv)dx 0,

отсюда следует

C

Gv 0.

 

 

 

 

x

дx

дt

 

 

 

 

дx

 

дt

1

 

 

 

 

 

 

 

 

 

 

 

 

110

4.2. Телеграфное уравнение

Рассмотрим уравнение

д2v

 

д2v

 

дv

 

LC

 

 

(RC GL)

 

GRv,

дx2

дt

2

 

 

 

дt

а также уравнение для силы тока i:

д2i

 

д2i

 

дi

 

LC

 

(RC GL)

 

GRi.

дx2

дt2

 

 

 

дt

В результате получим, что напряжение v и сила тока i – это решение одного и того же типа. Уравнение называют телеграфным уравнением

 

 

 

 

 

д2w

a

 

д2w

2b

дw

 

c

w,

 

 

 

 

 

 

 

дx2

0 дt2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0 дt

0

 

 

 

 

 

где a0 LC, 2b0 RC GL, c0

GR.

 

 

 

 

 

 

 

Если рассмотреть новую функцию u(х,t), положив

 

 

 

 

 

 

 

 

 

 

 

 

 

b0

t

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

w e a0 u,

 

 

 

 

 

 

 

то уравнение будет иметь простую форму

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

д2u

 

2

д2u

 

 

 

где

1

 

 

 

 

b02 a0c0

 

 

 

 

a

 

 

 

bu,

 

a

 

 

 

, b

 

 

.

 

дt2

 

дx2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a0

 

 

 

 

a0

 

 

4.3. Интегрирование телеграфного уравнения по методу Римана

Метод Римана используется для решения телеграфного уравнения с начальными условиями

u

 

t 0

f (x),

дu

 

F(x).

 

 

дt

 

 

 

 

t 0

 

 

 

 

 

 

 

 

 

 

Приведем телеграфное уравнение к каноническому виду, и используя новые независимые переменные ξ и η:

111

ξb (x αt),η b (x αt). a a

Телеграфное уравнение примет вид

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

L(u)

 

 

 

 

д2u

 

 

1

 

u 0.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

д д 4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Прямая t=0 в новых переменных . Можно записать

 

 

 

 

 

 

 

 

 

 

связи x

a

 

 

 

 

,

t

1

 

 

,

тогда

 

 

 

 

 

 

 

b

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

b 2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

дu

 

дu

 

 

1

 

дu

,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

д

 

 

 

 

д

 

 

b

 

 

 

 

дt

 

 

 

 

 

 

 

 

 

 

 

 

и из начальных условий получим

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a

 

 

дu

 

дu

 

 

 

 

1

 

дu

 

 

 

 

1

F(x)

1

F(

a

),

 

а также u

 

 

f (

).

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

д

 

д

 

 

 

b дt

 

t 0

 

b

 

 

 

 

b

 

 

 

 

b

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

b

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Пусть в формуле Римана а=0, b=0, f= 0, тогда имеем

 

 

 

u( , )

(uv)P (uv)Q

 

 

 

1

 

 

 

v(

дu

дu

)d

1

 

 

u(

дv

 

дv

)d .

 

 

 

 

 

 

 

 

 

 

 

 

 

0

0

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

д д

2

 

д д

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

QP

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

QP

 

 

 

 

 

 

 

 

 

Функция Римана v( , ; 0, 0) удовлетворяет сопряжен-

ному уравнению

д2v

1

 

v 0

 

и равна 1 на характеристиках

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

д д

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

MP и MQ. Решение нашего уравнения имеет вид: v G(( 0)( 0)).

Введем λ= ( 0)( 0) . Функция v=G есть реше-

ние уравнения

G ( ) 1 G ( ) G( ) 0.

Частным решением будет функция Бесселя нулевого порядка:

112

J0

( ) 1

2

 

4

 

6

...

 

(2 4)2

(2 4 6)2

 

22

 

 

 

Если v=J0(λ), то получим решение уравнения, которое справедливо на характеристиках ξ=ξ0 , η=η0 значения 1 (здесь

λ=0) и v J0(( 0)( 0)).

Вычислим

дv

 

 

 

 

 

0.5( - 0 )/(

( 0 )( 0 )),

д

 

 

 

 

 

 

 

 

 

дv

 

 

 

 

 

dJ0

 

 

 

 

 

 

 

 

 

 

 

0.5( - )/(

 

 

 

 

 

 

 

)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(

0

)( )

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

д

 

 

 

 

 

 

 

 

d

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

и найдём выражение

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

дv

 

 

 

 

 

дv

 

 

 

 

 

 

0.5( 0

- 0 )/(

 

 

 

 

).

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

( 0 )( 0 )

 

 

 

 

 

д

 

 

 

д

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Получим

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

u(p) f (a 0 /b), u(Q) f (a 0 /b),

 

 

 

 

 

 

 

 

 

f(a /b) f(a /b)

1

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

J0(

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

u( 0, 0)

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

( 0)( 0))F(a /b)d

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2b

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

J/

(

 

 

 

 

 

 

 

 

 

 

 

 

)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

(

0

)( )

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

f (a /b)d .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(

0

)( )

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

В старых координатах имеем

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

f(x at) f(x at)

 

 

1

 

x at

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

u(x,t)

 

 

Ф(x,t,z)dz,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

2

 

x at

 

 

 

 

 

 

 

 

 

 

 

где

 

 

 

 

 

 

 

 

 

 

(

 

 

 

 

 

)

 

 

 

( )

 

 

( )

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

( )

(( / )

(

) ( )

).

Ф

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

( , ,

) =

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

+

 

 

 

 

 

 

 

( ) (

)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

113

ЭЛЕМЕНТЫ ВАРИАЦИОННОГО ИСЧИСЛЕНИЯ

1. ОСНОВЫ ВАРИАЦИОННОГО ИСЧИСЛЕНИЯ

1.1. Вариация и ее свойства

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

Величина z называется функцией х (z=f(x)), если любому х из области определения поставлено в соответствие значение z. Такое же определение имеется для функции нескольких переменных.

Величина v называется функционалом зависящим от у(х), что обозначается так v=v[у(х)], если функции у(х) соответствует значение v.

Приращением x функции f(x) называется разностьx =х х1. Так как х – независимое переменное, то дифференциал х равен приращению (dx= x).

Вариацией y аргумента у(х) функционала v[y(x)] назы-

вается y=у(х)–у1(х).

Функция f(x) называется непрерывной, если малому изменению х соответствует малое изменение f (х).

Функционал v[у(х)] будет непрерывным, если малому изменению y(x) соответствует малое изменение v[у(х)].

Функции у(х) и у1(х) будут близкими, если модуль разности у(х) –у1(х) мал для всех х. Можно считать близкими только те функции, для которых малы модули разностей

y(k) (x) y1(k) (x)..

Поэтому вводятся определения близости кривых у=у(х) и у1=у1(х), так кривые у и у1 имеют близость нулевого порядка, если мал модуль разности уy1. Кривые у=у(х) и у1=у1(х) имеют

114

близость k-го

порядка, если малы модули разно-

стей y(x) y1(x),

y (x) y1(x),..., y(k) (x) y1(k) (x).

Ниже изображены кривые (рис. 1.1), близкие в смысле близости нулевого порядка.

y

B

A

x

Рис. 1.1. Пример кривых с близостью нулевого порядка

y

B

A

x

Рис. 1.2. Пример кривых с близостью первого порядка

На рис. 1.2 кривые имеют близость первого порядка. Функция f(х) непрерывна при х=х0, если для любого >0

можно подобрать δ>0 такое, что

f (x) f (x0)

, при

x x0

.

Функционал v[y(x)] непрерывен при у=у0(х) (близость k- го порядка), если для любого >0 можно подобрать δ>0 такое, что

v[y(x)] v[y0(x)] .

115

Функционал v[у(х)] непрерывен при у= у0(х) (близость k- го порядка), если для любого >0 можно подобрать δ>0 такое, что

y(x) y0(x)

 

,...,

y(k)(x) y0(k)(x)

.

 

 

 

 

 

 

При этом функция у(х) берется из класса функций, на котором функционал v[у(х)] определен.

Если определено понятие расстояния (у1,у2) между кривыми y1=y1(x) и y2=y2(x) (x0 x x1),то близкие кривые это кривые, расстояние между которыми мало.

Если (у12)= max

y

(x) y

 

 

(x)

, то приходим к поня-

x x x

1

 

 

2

 

 

 

 

 

0

1

 

 

 

 

 

 

 

 

 

 

тию близости нулевого порядка.

 

 

 

 

 

 

 

Если считать

k

 

 

 

 

 

 

 

 

 

 

(у1,у2)=

 

 

 

 

 

y p

 

 

 

 

 

 

max

 

 

(x) y

p

(x)

 

 

 

x

x x

 

 

1

 

2

 

 

 

p 1

 

0

1

 

 

 

 

 

 

 

(у1 и у2 имеют непрерывные производные до порядка k включительно), то близость кривых будет k-го порядка.

Линейной функцией называется функция l(х): l(kx)=kl(х),

где k - произвольная постоянная и

l(х1 + х2)=l(х1)+l(х2).

Линейным функционалом называется функционал L[y(x)], если L[ky(x)]=kL[y(x)], где k – произвольная постоянная и L[y1(x) + y2(x)= L[y1(x)]+L[y2(x)].

Примером линейного функционала является

x1 x2

L[y(x)] p(x)y q(x)y ]dx.

x0 x1

Приращение функции

f f (x x) f (x)

можно представить в виде

f A(x) x (x, x) x,

где А(х) не зависит от x , (x, x) 0, при x 0.

116

Функция называется дифференцируемой, если существует ее производная. Линейная по отношению к x часть при-

ращения А(х) x называется дифференциалом функции и обо-

значается df. Можно получить А(х)=f’(х) и df = f'(x) x .

 

Приращение

функционала [y(x) y] [ y(x)]

можно представить в виде

 

L[y(x) y] (y(x), y)max

 

y

 

,

 

 

 

 

где L[y(x), y ]

- линейно по отношению к y. Пусть задан

функционал

mах| y | (максимальное значение

| y | и

(y(x), y) 0

при

max

 

y

 

0). Линейная величина по от-

 

 

ношению к значению | y |, т. е. L[y(x), y ] называется вариаци-

ей функционала ( v).

+αΔ )

при фиксированном х и

 

 

 

Рассмотрим

 

и .

Если α=1, то получим

, а при α=0 получим исходное

 

 

 

 

 

 

(

 

 

 

 

 

 

по

α

значение функции f(х). При этом производная

 

 

 

 

 

 

 

 

(

+Δ )

 

 

 

 

 

α

при α=0 равна дифференциалу функции f(х) в точке( +х, т. е).

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

f (x x)

 

f

 

 

 

 

 

 

 

 

 

(x x) x

0 f (x) x df (x).

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Для функции z=f(x1,x2,…,xn) можно получить

 

 

и при α=0

 

(

+

 

,…,

+

)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n

f

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

f (x1

x1,...,xn xn )

 

 

 

 

x df.

 

 

 

 

 

 

 

x

 

 

 

 

 

 

 

 

 

 

 

 

0

i 1

i

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Для функционалов вида v[y(x)], зависящих от нескольких неизвестных функций или от функций нескольких перемен-

ных, можно

ввести вариацию как производную от

[y(x) y]

по α при α=0, т. е.

v [y(x) y] v[y(x)] L(y, y) max y.

Производная [y(x) y] по α при α=0 будет

117

lim

v

lim

v

lim(L(y, y) [y(x) y]

 

 

 

max

 

y

 

)/

 

 

 

 

 

 

 

0

 

0

 

0

 

 

 

 

 

 

 

 

 

 

 

 

lim

L(y, y)

lim( [y(x) y]

 

 

 

max

 

y

 

)/ L(y, y)

 

 

 

 

 

 

 

 

0

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

в силу линейности

L(y,αδy)=αL(yy),

 

 

 

 

 

 

lim( [y(x) y]

 

 

max

y

)/ lim( [y(x) y]max

 

y

 

) 0,

0

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

т. к. (y(x), y) 0 при

 

 

 

 

 

 

 

 

 

Дифференциал

функции f(х):

 

α → 0.

 

 

 

 

 

 

 

f (x x)

,

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

а вариация функционала v[y(x)]:

v[y(x) y] .

0

Определение. Функционал v[y(x)] достигает на кривой у=у0(х) максимума, если функционал v[y(x)] на близкой к у=у0(х) кривой не больше v[y0(x)], то

lim [y0 (x)]v[y(x)] 0.

n

Теорема 1. Если функция f(х) (дифференцируемая) имеет максимум или минимум во внутренней точке x=x0, то в этой точке df=0.

Теорема 2. Если функционал v[у(х)] (имеющий вариацию) достигает максимума или минимума при y=y0(x), где у(x) - внутренняя точка функционала, то при y=y0(x), v 0.

Если функционал v[y(x)] достигает на кривой у=у0(х) максимума или минимума, для которых модуль разности у(х)–y0(х) мал, то максимум или минимум называется сильным.

Если же функционал v[y(x)] достигает на кривой у=у0(х) максимума или минимума по отношению к кривым у=у(х), близким к у=у0(х) (близость первого порядка), то максимум или минимум называется слабым.

118

1.2. Уравнение Эйлера

Рассмотрим функционал

x1

[y(x)] F(x, y(x), y/ (x))dx,

x0

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

у(х0)=у0 и у(х1)=у1 (рис. 1.3).

y В

А

y1

y0

x

0

x0

x1

Рис. 1.3. Функционал [y(x)] с закреплёнными концами

Функция F(x,у,у') будет трижды дифференцируема. Необходимым условием экстремума будет равенство ну-

лю вариации функционала.

Пусть экстремум достигается на кривой у=у(х). Возьмем близкую к у=у(х) кривую у= y (х) и рассмотрим

y(x, ) y(x) (y(x) y(x)).

При =0 получим кривую у=у(х), а при =1 имеем у= y (х) (рис. 1.4). Разность y(x) y(x) называется вариацией функции у(х) и обозначается y.

119