Метод Фур'є (відокремлення змінних)

Матеріал з testwiki
Перейти до навігації Перейти до пошуку

Метод Фур'є - один з основних методів розв'язування диференціальних рівнянь з частинними похідними. І незважаючи на те, що він доволі громіздкий його знання в будь-якому разі необхідне щоб здати СРП.

На прикладі задачі про вільні коливання струни, кінці якої жорстко закріплені.

Ця умова означає, що на струну не діють зовнішні сили, і вона рухається під дією власної сили натягу, яка з'явилась через початкове відхилення, та початкової швидкості. Кінці не рухаються. В рівняннях це означає: Шаблон:Eq

вільні коливання Шаблон:Eq початкове положення та швидкість Шаблон:Eq закріплені кінці.

Відокремлення змінних

Спочатку шукаємо частинні розв'язки рівняння Шаблон:Er, які нерівні тотожньо нулеві, і задовольняють умови Шаблон:Er, у вигляді: Шаблон:Eq

Через формулу Шаблон:Er метод назвали методом відокремлення змінних. Ця формула - один з ключових моментів.

Підставляємо Шаблон:Er в Шаблон:Er, і отримуємо Шаблон:Eq Ділимо обидві частини на a2X(x)T(t), дістанемо Шаблон:Eq

Для того, щоб функція u(x,t)=X(x)T(t) була розв'язком рівняння Шаблон:Er рівність Шаблон:Er повинна задовольнятися тотожно, тобто для всіх значень незалежних змінних 0<x<l. Права частина рівності Шаблон:Er є функцією тільки змінної t, а ліва - тільки x. Фіксуючи деяке значення x, і змінюючи t (або навпаки), дістанемо, що права і ліва частини рівності Шаблон:Er при зміні своїх аргументів зберігають стале значення. Шаблон:Eq

де λ - стала, яка для зручності вибирається зі знаком "мінус", при цьому нічого не припускається про її знак.

Із співвідношення Шаблон:Er дістаємо рівняння для визначення функцій X(x), та T(t):

Шаблон:Eq Шаблон:Eq


Крайові умови Шаблон:Er дають: Шаблон:Eq Шаблон:Eq

Звідси випливає, що X(0)=X(l)=0, бо інакше було б Шаблон:Eq але ж задача полягає у знаходженні нетривіальних розв'язків.

Таким чином для функції X(x) одержали так звану задачу Штурма-Ліувілля (задачу про власні значення): знайти ті значення параметра λ, при яких існують нетривіальні розв'язки задачі: Шаблон:Eq

а також знайти ці розв'язки. Такі значення параметра λ називаються власними значеннями, а відповідні їм нетривіальні розв'язки - власними функціями задачі Шаблон:Er.

Розглянемо окремо випадки, коли параметр λ від'ємний, дорівнює нулю, або додатній.

Випадок λ<0

Загальний розв'язок рівняння Шаблон:Er має вигляд

Шаблон:Eq

Крайові умови дають Шаблон:Eq Шаблон:Eq тобто C1=C2 і C1(eα+eα)=0. Але в розглянутому вище випадку α - дійсне і додатне, а тому eαeα0 Тому C1=0 і C2=0. Отже, X(x)0 і, значить в цьому випадку задача не має нетривіальних розв'язків.

Випадок λ=0

В цьому випадку загальний розв'язок рівняння Шаблон:Er має вигляд: Шаблон:Eq

Крайові умови дають Шаблон:Eq Шаблон:Eq тобто C1=0. Отже, C1=C2=0, і значить X(x)0, тобто також не існує нетривіальних розв'язків.

Випадок λ>0

І нарешті останнє. Загальний розв'язок рівняння Шаблон:Er в такому разі виглядає: Шаблон:Eq

Крайові умови дають Шаблон:Eq Шаблон:Eq Якщо C1=0, то і в цьому випадку буде X(x)0. Ми шукаємо нетривіальні розв'язки, тому C10 і Шаблон:Eq або λ=πnl, де n=±1,±2,,n0, тому що за умовою λ>0.

Отже, нетривіальні розв'язки задачі Шаблон:Er можливі лише при значеннях

Шаблон:Eq

Цим власним значенням відповідають власні функції

Шаблон:Eq де Dn - довільна стала.

Надалі будемо надавати n тільки додатні значення, оскільки при від'ємних n дістанемо розв'язки, які мають такий самий вигляд (адже Dn - довільні сталі, які можуть мати будь-який знак).


Знаходимо T(t)

Підставивши знайдені значення λn в Шаблон:Er, дістанемо рівняння Шаблон:Eq загальний розв'язок якого має вигляд Шаблон:Eq де Cn та Kn - довільні сталі.

Частинні розв'язки

Підставивши Шаблон:Er і Шаблон:Er в Шаблон:Er знайдемо частинні розв'язки рівняння Шаблон:Er, які задовольняють крайові умови Шаблон:Er. При цьому кожному значенню n=1,2, буде відповідати розв'язок

Шаблон:Eq

Введемо позначення An=CnDn,Bn=KnDn. Тоді un(x,t) можна записати у вигляді

Шаблон:Eq

Функції un(x,t) задовольняють початкові умови Шаблон:Er вихідної задачі тільки для окремих випадків початкових функцій φ(x) і ψ(x).

За допомогою розв'язків Шаблон:Er побудуємо розв'язок, який задовольняє початкові умови Шаблон:Er. Розв'язок рівняння Шаблон:Er, який задовольняє умови Шаблон:Er і Шаблон:Er шукатимемо у вигляді ряду:

Шаблон:Eq

Якщо ряд Шаблон:Er збігається рівномірно, і його можна двічі продиференціювати по x і по t почленно, то сума ряду буде задовольняти рівняння Шаблон:Er і крайові умови Шаблон:Er. Поставимо вимогу, щоб функція u(x,t), яка визначена рядом Шаблон:Er задовольняла початкові умови Шаблон:Er:

Шаблон:Eq

Ряд Фур'є

Припустимо, що функції φ(x) і ψ(x) задовольняють умови розкладу в w:ряд Фур'є за синусами на проміжку (0,l). Тоді

Шаблон:Eq Шаблон:Eq

Порівняння рядів Шаблон:Er і Шаблон:Er з формулами Шаблон:Er показує, що для виконання початкових умов треба покласти

Шаблон:Eq

Таким чином розв'язок задачі Шаблон:Er-Шаблон:Er, одержаний методом відокремлення змінних, має вигляд

Шаблон:Eq

Це був розв'язок однорідного рівняння гіперболічного типу. І хоча це не всі класи рівнянь, але спосіб яким розв'язано це - найзагальніший, і входить у розв'язки багатьох інших. Проте людина, яка оцифровувала цю методичку вже здала сесію, і іде готуватись до наступної, тому продовження читайте з тридцятої сторінки методички зі списку літератури. Формула (19).

Література

  1. Методичні розробки до вивчення нормативного курсу "Рівняння математичної фізики" (класифікація рівнянь другого порядку, формули Даламбера і Пуассона, метод відокремлення змінних) для студентів факультету кібернетики. Упоряд. С.О.Войцеховський, В.І.Гаркуша, М.П.Копистра та ін.- Київ, 2001 р. 65 с.