киб_3_семестр_матан_лекции
.pdfЕсли к правой части равенства (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 |
|
kπ |
|
|
|
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 (x)ϕk (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 из этого же пространства, то такое пространство является полным.