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

киб_3_семестр_матан_лекции

.pdf
Скачиваний:
20
Добавлен:
10.05.2015
Размер:
4.03 Mб
Скачать

Если к правой части равенства (4) добавить f(x,t), получится неоднородное уравнение; – некоторая константа.

Струны с положением равновесия вдоль оси Ох:

Концы струны закреплены в т. x=0; x=l. Надо найти U(x,t) – положение точки струны с абсциссой х в момент времени t.

Постановка задачи Решить уравнение (4)

(уравнение колебания струны) При граничных условиях Граничные условия : U(0,t)=U(ℓ,t)=0 -концы струны

закреплены; (5)

Начальные условия :

 

U(x,0)=f(x)

– форма струны в начальный момент времени;

 

 

 

(6)

 

U (x,0) = Y (x)

- скорость, которую сообщили струне в начальный

 

t

 

 

 

 

 

0≤ x ≤ ℓ

момент времени

Решение задачи методом Фурье:

Ищем решение в виде:

 

1) ,

 

(7)

2) подставим (7)формулу в (4), получим уравнение с разделяющимися

переменными

 

 

T"(t)

 

=

X "(x)

= -λ

(8)

 

a 2 ×T (t)

 

 

 

X (x)

 

T"(t) + λ × a 2 ×T (t) = 0

(9)

X "(x) + λ × X (x) = 0

(10)

Решим (10), используя формулу (5) с учетом (7), получим

X (0) = X (l) = 0

(11)

Задача Штурма-Лиувилля :

Определение1(задачи Штурма-Лиувилля). Найти нетривиальные решения (10) с условием (11) называется задачей Штурма-Лиувилля.

Определение2: (собственных значения и собственных функций)

λ удовлетворяющее нетривиальному решению задачи Штурма-Лиувилля

называются собственными значениями, а X(x), полученные при этих λ ,

носят название собственных функций.

Решение уравнения колебаний струны методом Фурье

2U

= a 2 2U

(1)

t 2

x 2

 

Граничные условия:

 

U(0,t)=U(ℓ,t)=0(2)

(2)

Начальные условия:

U (x,0) = f (x)

U (x,0) = ϕ (L) ∂t

0 ≤ x L

Метод Фурье:

U (x, t) = X (x) ×T (t)

T "(t) = X "(x) = −λ a 2T (t) X (x)

T "(t) + λa 2T (t) = 0

Задача Штурма-Лиувилля:

X "(x) + λX (x) = 0

X (0) = X (L) = 0

(3)

(4)

(5)

(6)

(7)

(8)

Нужно найти нетривиальное решение задачи(7), удовлетворяющее условию

(8) :

1) λ < 0, λ = −μ 2 , μ > 0

В этом случае (7) имеет решение: X (x) = C1e− μx + C2eμx

Из условия (8) вытекает система: X (0) = C1 + C2 = 0

X (l) = C1e−μl + C2eμl = 0

Мы получили однородную систему уравнений относительно С1 иС2 . составляем определитель этой системы относительно С1 иС2 .

D =

 

1

1

¹ 0

С1=0

X (x) ≡ 0

e

−μl

e

μl

С2 =0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2)второй случай, когда λ = 0

X (x) = C1 + C2 X

X (0) = C1 = 0 , X (l) = C1 + C2l = 0

C2 = 0

X ( x) ≡ 0

3)третий случай, когда λ > 0

λ= μ 2 , μ = λ

X (x) = C1 cos μx + C2 sin μx = C1 cos λ x + C2 sin λ x

X (0) = C1 = 0

X (l) = C1 cos λl + C2 sin λl = 0 → C2 sin λl = 0

Поскольку С2 не должно роняться нулю, то sin λl = 0 . Надо найти такие λ , при которых sin λl = 0

 

 

π 2k 2

λl = πk ; λk =

l

2

 

 

 

Равенство (9) задает нам собственные значения задачи Ш.-Л.

Теперь нужно найти собственные функции X k (x) . Для дальнейших

выкладок положим C2 = 1 , тогда X k (x) = sin πk lx

Дальше ищем решения задачи (6)при собственных значениях

(9)

(10)

T (t) = A cos

kπa

t + B sin

kπa

t

 

 

 

 

,

(11)

k

k

l

k

l

 

 

 

 

 

Где Ak и Bk неопределенные коэффициенты, которые еще нужно найти.

Подставляем (11) и (10) в уравнение (4)

U

 

(x,t) = ( A cos

kπa

t + B sin

kπa

t) sin

πk x

 

k

 

 

(12)

 

k

l

k

l

l

 

 

 

 

 

Равенство (12) уже является решение уравнения колебания струны, удовлетворяющее граничным условиям (2).

Поскольку (1)-линейное и однородное, то

 

kπa

 

 

kπa

t) sin πk x (13)

U (x, t) = U k (x, t) = ( Ak

cos

t + Bk

sin

 

l

k =1

k =1

 

l

 

l

(при условии, что ряд сходится равномерно, т.к. тогда его можно дважды дифференцировать при подстановке в (1) )

3) найдём , используя

начальные условия (3)

 

, 0

sin

 

 

 

 

 

 

 

 

 

 

 

A

=

2

l f (x) sin πk xdx

 

 

 

 

l

 

 

 

 

где k

 

l

 

 

 

(15)

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

Используя второе начальное условие, продифференцируем ряд

 

∂U

 

 

πka

πk

 

 

 

(x,0)

= ϕ (x) =

Bk

sin

x

 

(16)

∂t

 

 

k =1

l

l

 

 

 

- этот ряд (16) так же является р.Фурье

 

 

2

%

! ∙ sin

 

2

%

 

17

 

&

"

 

! ∙ sin

"

 

 

 

 

 

&

 

 

Вывод: решение уравнения (1) с граничными условиями (2) и начальными условиями (3) задается формулой (13), где коэффициенты Ak , Bk ищутся по формулам соответственно (15) и (17), как коэф.Фурье.

Замечание 1: если начальная скорость ϕ ( x) = 0 , то Bk = 0

 

kπa

 

 

 

U (x, t) = Ak

cos

t × sin

x =

(18)

l

 

k =1

 

 

l

 

1

 

kπ

 

kπ

 

 

1

[ f (x at) + f (x + at)]

 

= Ak

 

sin

 

(x at) + sin

 

(x + at)

=

 

(19)

2

l

l

2

k =1

 

 

 

 

 

 

Где под f(x) понимается продолжение начального условия f(x) нечетно на (-ℓ;0), а затем периодически повторяется на всю прямую и с физической точки зренияэто представляет собой полусумму двух разбегающихся волн. На этом результате (19) основан метод Даламбераэто метод бегущих волн.

Замечание2: ( о физическом смысле решения задачи методом Фурье) Из(13) следует (в более компактной форме), что

 

 

 

 

 

 

 

sin πk x ×sin(πka t + ϕk ) ,

 

 

U (x, t) = Dk

 

(20)

 

 

 

 

 

 

k =1

 

l

 

l

 

 

 

 

 

 

 

 

 

 

 

=

 

 

 

= arctg

 

Ak

 

 

 

 

 

где

D

k

 

A2 + B 2

k

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

k

k

 

 

Bk

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Каждое слагаемое суммы (20) – есть гармоника;

 

C

k

= D sin πk x

- амплитуда гармоники;

 

 

 

 

 

k

 

 

l

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

πka

 

 

 

 

 

 

 

 

 

 

 

 

 

ω

k

 

=

- собственная частота гармоники

 

 

 

 

 

 

 

L

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1) найдем точки нулевой амплитуды

 

 

 

 

 

 

 

 

πk x = nπ

 

Xn =

n

l - узловые точки (или узлы)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

l

 

 

 

 

 

 

 

k

%

 

%

 

 

 

 

 

 

 

 

(n=0,1,2…k); & 0,

; …

;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Т.е. струна длины l разбивается на k участков, колеблющихся по закону (12) 2) точки max Сk

πkx = π (2n + 1) l 2

Xnmax = 2n + 1 l ((n=0)….(k-1)) – точки пучности (середины отмечены

2k

выше участков) колеблются с наибольшей амплитудой Ck.

3) на соседних участках отклонения струны различны по знаку.

Определение 3: (стоячие волны)

Явление, для которого выполняются все три условия, называется

стоячими волнами.

Гармоническое колебание с ω1 получается при k=1)

ω = πa

 

1

l

- основной тон струны, а все остальные (к=2,…)- обертоны.

 

U1 (x, t)

Основной тон определяется составляющими U1(x,t). Суммарное действие отдельных тонов дает тембр звука струны, т.е. суммарной гармоники, накладываясь друг на друга, определяют тембр струны.

k=1

k=2

k=3

k=4

Лекция 19.

Понятие гильбертова пространства. Ортогональные и

ортонормированные системы функций. Обобщенный

ряд Фурье по полной ортогональной системе

функций.

Рассмотрим примеры линейных бесконечномерных функциональных пространств.

L : 1) линейное пространство ( бесконечно мерное) функций L=С[a,b] непрерывных на [a,b];

2) L = [a,b] –линейное пространство функций, кусочнонепрерывных на [a,b];

Со скалярным произведением в следующем виде:

b

 

( f , g) = f (x)g(x)dx

(1)

a

Заданием на дом:

Проверить, что для (1) выполняются все свойства скалярного произведения. Замечание 1. В правой части (1) стоит интеграл Римана для функций с интегрируемым квадратом, т.е. выполняется

 

b

b

f 2 (x)dx < +∞

g 2 (x)dx < +∞

 

a

a

В пространство L введем норму (длину элемента пространства)

Определение 4: (нормы)

 

 

 

 

 

=

 

 

 

 

 

f

 

( f , f )

(2)

Свойства нормы:

1) f > 0, f = 0 , тогда и только тогда, когда f=0.

2)λf = λ f , λ - вещественное число.

3) f + g < f + g (неравенство треугольника)

Определение5(ортогональной системы функций). Бесконечная система действительных функций {ϕn } из рассматриваемого пространства

L[a;b]называется ортогональной, если все функции этой системы попарно

ортогональны,

b

 

т.е. ϕi (xk (x)dx = 0 , при i ¹ k , т.е (ϕi ,ϕk )=0 i ¹ k

(3)

a

Определение 6(ортонормированной системы функций).

 

 

Система { ( )} называется ортонормированной, если

( , ) = 0, ≠

( ) ( ) = 0, ≠

 

т. е.

( , ) = 1, =

 

( ) = 1, =

 

 

(4)

 

Примеры ортогональных систем функций:

 

 

Тригонометрические системы функций:

 

1)

1,cos(nx), sin(nx) - ортогональна на (-π; π)

ненорми-

 

(полная, в лин.пространстве функций с интегрируемым

рованные

 

квадратом на(-π; π)

 

системы ф-ий

2)

1,cos

nπx

, sin

nπx

 

- ортогональна на (-ℓ;ℓ)

 

 

l

l

 

 

 

 

 

 

 

 

 

(полная на (-ℓ;ℓ))

 

 

3)

sin

nπx

 

- ортогональна на [0;ℓ]

 

L

 

 

 

 

 

 

 

 

 

Нормы функций, входящих в систему 1) имеют вид:

 

 

||1||2=2L , ||cos!"#%

$&||2 =L, ||sin!"#% $&||2 =L

 

 

||1||2=2π , ||cos(nx)||2 =π, ||sin(nx)||2 =π, а в системе 2)

 

Примеры других систем:

полином Лежандра – орт. на (-1;1) функции Чебышева – орт. на (-1;1) функции Эрмита – орт. на (-∞;+∞) функции Лагерра – орт. на [0;+∞) и другие

Изучить по указанной ранее книге Приложение 1 стр. 82-96

Определение 7(полной системы):

 

 

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

 

называется полной, в

 

пространстве"

функций с интегрируемым

соответствующем линейном

 

{ }

 

 

квадратом, если не существует в этом пространстве ненулевого{ } элемента такого, который был бы ортогонален каждой функции " из системы. Пример1. Базис в конечномерном евклидовом пространстве является

полной системой. Здесь только 0-вектор ортогонален каждому вектору базиса Пример2. Тригонометрические системы 1),2),3) каждая в своём

пространстве, на своём отрезке.

Определение 8(сходимость по норме):

{fn}Є L сходится по норме к элементу f, если выполняется следующее условие:

 

 

 

 

 

 

 

 

 

 

"

 

 

0

 

"

 

"→,1223 0

||) * )|| "→,-./

 

 

()( ) * ) ( ))

 

 

 

 

 

 

 

 

 

 

(если это выполняется, то говорят, что {)"}сходится по норме)

Вопрос:

{)"}

 

 

 

 

 

 

Замечание 2:

сходимость по норме называют еще сходимостью в среднем

Замечание 3:

 

 

 

Є L [a,b]

 

 

 

 

Всегда ли

 

из[a,b] будет сходиться по норме к f Є L[a,b]?

всегда,"

так как существуют последовательности, которые могут

Ответ: не {) }

 

 

 

 

 

 

 

 

сходиться к элементу f Є другому пространству, а не пространству L[a,b].

Пример 1:

Частичные суммы Sn(x) Є C[a,b] не всегда сходятся к непрерывной функции f Є C[a,b], иногда сходится к f Є [a,b], т.е. к кусочно-непрерывной функции. В результате мы можем получить предельные элементы, которые не принадлежат4 (исходному) 6 78; :;пространству.

Пример 2: 5 сходимость не к кусочно-непрерывным функциям, а к обобщенным функциям (дельта функция Дирака).

45

( ) = < 1

, | | ? >

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

 

 

2>

| | @ >

 

0,

При ε =

1

 

1

 

 

 

45

 

 

 

28

( ) = A 1, | | ? 21

 

 

2

 

 

 

 

0,

| | @ 2

 

Пределом этой последовательности кусочнонепрерывных функции является обобщенная4 * функция

функция Дирака –

 

 

 

 

 

 

 

,

 

 

Определение 9: 4( ) обобщенная функция Дирака(английский физик)

4( )

4( ) = 0, при ≠ 0,

 

4( ) = 1

 

 

 

 

 

 

 

 

- плотность единичной массы,

 

С физической точки зрения

L,

 

 

 

 

 

 

 

 

“0”.

 

 

сосредоточенной в точке

 

4( )

 

 

Свойства:

 

 

 

 

 

 

M

N

,

O = 4(

 

 

 

 

 

 

1) интеграл Фурье имеет вид :

 

 

#

L,

P$

)

 

 

 

 

 

4( ) * производная

4( ) = VW( ), где

h(x)

f(x) =

1, x>0

 

 

функции Хевисайда

2) функция Хевисайда:

0, x<0

Замечание4. Теория обобщенных функций была построена советским математиком С.Л. Соболевым и французским математиком Шварцем.

Определение10(полное линейное пространство).

Если любая последовательность {fn} сходится по норме к f из этого же пространства, то такое пространство является полным.