|
|
308. sor: |
308. sor: |
|
| |
|
| <math>u(0, t) = u(3, t) = 0,~u(x,0)=sin\frac{4\pi}{3}x,~\frac{\partial u}{\partial t}(x, 0) = 2\sin\frac{\pi}{3}x</math> | | <math>u(0, t) = u(3, t) = 0,~u(x,0)=sin\frac{4\pi}{3}x,~\frac{\partial u}{\partial t}(x, 0) = 2\sin\frac{\pi}{3}x</math> |
| | |
| | {{Rejtett |
| | |mutatott='''Megoldás:''' |
| | |szöveg= |
| | * Az <math>U(x, t)</math>-t keressük szorzat alakban: <math>U(x, t) = X(x)T(T)</math> |
| | * A diffegyenlet így átírva: <math>X(t)\ddot{T}(t) = 4*X''(x)T(T)</math> |
| | * Ez így már szeparálható (figyeljünk arra, hogy a deriváltak a számlálóban legyenek): |
| | <math>4 \cdot \frac{X''(x)}{X(x)} = \frac{\ddot{T}(t)}{T(T)} = -b^2</math> |
| | * Nézzük meg, hogy melyik változóra van feltételünk, aminek a jobb oldalán konstans szerepel. |
| | ** Az első két féltétel átírva: X(0)T(t) = X(3)T(t) = 0, minden t-re, vagyis X(0) = X(3) = 0 |
| | ** Tehát az X-re van a T-től nem függő feltételünk, ezért először az X-re oldjuk meg a diffegyenletet! |
| | * Oldjuk meg a diff-egyenletet: |
| | <math>4 \cdot \frac{X''(x)}{X(x)} = -b^2</math> |
| | |
| | <math>4 \cdot X''(x) + b^2 \cdot X(x) = 0</math> |
| | |
| | * Írjuk fel a karakterisztikus függvényt! |
| | <math>4 \cdot \lambda^2 + b^2 = 0</math> |
| | |
| | <math>\lambda^2 = -\frac{b^2}{4}</math> |
| | |
| | <math>\lambda = \pm i \frac{b}{2}</math> |
| | |
| | * Vagyis a diff-egyenlet megoldása: |
| | <math>X(x) = c_1 \cos{\frac{b}{2}x} + c_2 \sin{\frac{b}{2}x}</math> |
| | |
| | * Vizsgáljuk meg a kezdeti feltételeket: |
| | <math>X(0) = c_1 \cos{0} + c_2 \sin{0} = c_1 = 0</math> |
| | |
| | <math>X(3) = c_2 \sin{\frac{b}{2}x} = 0</math> |
| | |
| | Ami csak olyan egész k értékekre teljesülhet, amikre: <math>\frac{b}{2} = k \pi,~b = 2 k \pi</math> |
| | |
| | * Most oldjuk meg a diff-egyenletet T(t)-re, de a b helyére az újonnan kapott képletet írjuk be. |
| | |
| | <math>\frac{\ddot{T}(t)}{T(t)} = -(2 k \pi)^2</math> |
| | |
| | <math>\lambda^2 = -(2 k \pi)^2</math> |
| | |
| | <math>\lambda = \pm 2 i k \pi</math> |
| | |
| | * A T-re vonatkozó (k-tól függő) diff-egynelet: |
| | <math>T_k(t) = a_k \cos{2 k \pi t} + b_k \sin{2 k \pi t}</math> |
| | |
| | * Az <math>U(x, t)</math>-re vonatkozó k-tól függő egyenlet tehát: |
| | <math>U_k(x, t) = c_2 \sin{k \pi x} (a_k \cos{2 k \pi t} + b_k \sin{2 k \pi t})</math> |
| | |
| | * Vezessük be az <math>A_k = c_2 \cdot a_k</math> és <math>B_k = c_2 \cdot b_k</math> konstansokat! |
| | <math>U_k(x, t) = A_k \sin{k \pi x} \cos{2 k \pi t} + B_k \sin{k \pi x} \sin{2 k \pi t}</math> |
| | |
| | * Az <math>U(x, t)</math> pedig felírható az <math>U_k(x, t)</math>-k összegeként az összes k-ra. |
| | <math>U(x, t) = \sum_0^\infty U_k(x, t)</math> |
| | |
| | * A maradék két feltétel segítségével számoljuk ki az <math>A_k</math> és <math>B_k</math> konstansok értékeit. |
| | }} |
|
| |
|
| <hr> | | <hr> |
Integrál trafók témakör
Laplace trafó diff-egyenlet
1) [2015ZH1] Laplace transzformáció segítségével számítsuk ki x(t)-t, ha
Megoldás:
- Vegyük mindkét egyenlet Laplace trafóját ():
- Az egyenleteket átrendezve, és x(0), y(0)-t behelyettesítve:
- Megoldás X-re (a számlálóban a mátrix első oszlopa le lett cserélve az egyenlet jobb oldalára. Ha y-t számolnánk, akkor a második oszlopot kéne lecserélni):
- Az inverz laplacehoz bontsuk parciális törtekre:
- Együtthatókat összehasonlítva:
- Vagyis
- Tehát a táblázat alapján
2) [2016ZH1] Laplace transzformáció segítségével számítsuk ki x(t)-t, ha
Megoldás:
- Vegyük mindkét egyenlet Laplace trafóját:
- Átrendezve és mátrixos alakra hozva:
- Inverz Laplace után:
3) [2016ZH1] Transzformáljuk elsőrendűvé a differenciálegyenlet Laplace transzformációval (Nem kell megoldani!)!
Megoldás:
- Számítsuk ki a tagok Laplace trafóját (x szerint):
- Tehát az egyenlet Laplace transzformáltja (elsőrendű Y-ban):
Laplace trafó szabályok alkalmazása
1) [2016PZH] Számítsuk ki az alábbi jobboldali határétrékeket:
Fourier diff-egyenlet
1) [2015ZH1] Oldjuk meg Fourier transzformáció segítségével!
Megoldás:
- Vegyük az egyenlet Fourier trafóját (a táblázatban a Fourier trafó y függvénye, de az y itt mást jelent, a táblázatbeli y-ok helyére írjuk s-t, illetve vezessük be az alábbi jelölést: )!:
- Aminek a disztribúció értelemben vett megoldás Y-ra:
- Ha , akkor leoszthatunk vele.
- Ha , akkor , vagyis bármilyen konstans lehet, ezt jelöljük pl c-vel.
- Az összeg jobboldali tagja egyszerűsíthető, ha kihasználjuk, hogy az egy disztribúció (a a nevezőben lévő s-be is nullát helyettesít):
- Aminek vegyük az inverz Fourier transzformáltját:
- Megjegyzés: A táblázatban szerepel , de nekünk inverz trafó kell
2) [2016ZH1] Transzformáljuk elsőrendűvé a differenciálegyenlet Fourier transzformációval (Nem kell megoldani!)!
Megoldás:
- Számítsuk ki az egyenlet tagjainak Fourier trafóját (x szerint):
- Vagyis az egyenlet Fourier trafója (elsőrendű diff-egyenlet -ra):
Fourier trafó szabályok alkalmazása
1) [2015ZH1] Számítsuk ki az Fourier transzformáltját, ha tudjuk, hogy
Megoldás:
Vezessük be a jelölést!
Disztribúciók
1) [2015ZH1] Adjuk meg és lineáris kombinációjaként az disztribúciót!
2) [2016ZH1] Számítsuk ki a reguláris disztribúcuó és a disztribúció konvolúciójának hatását a függvényre:
Megoldás:
- Elődáson volt, hogy
- Ezt felasználva alkalmazzuk a disztribúciót a függvényre:
3) [2016ZH1] Mi az disztribúció értelemben vett egyenelet összes megoldása? (+1 miért?)
Megoldás:
- Ha , akkor leoszthatunk vele, és azt kapjuk, hogy .
- Ha , akkor , vagyis bármilyen konstans értéket felvehet, ezt jelöljük pl c-vel.
- Tehát ha , akkor , ha , akkor tetszőleges értékű, ez röviden:
4) [2016ZH1] Adjuk meg az disztribúciót a eltolt deriváltjainak lineáris kombinációjaként!
Megoldás:
5) [2016PZH] Legyen u az által generált reguláris disztribúció, . Számítsuk ki -t!
Megoldás:
- Először szabaduljunk meg a konvulúciótól:
- Majd értékeljük ki a disztribúciót (ez egy közismert integrál, de viszonylag nehéz kiszámolni):
Wavelet trafók
Megjegyzés: a kitevőbe írt törtek (pl: ) sok böngészőben hibásan jelennek meg, ezért ezekben az esetekben törtek helyett osztás jelet fogok használni.
1) [2015ZH1] Legyen , a mexikói kalap wavelet.
a) Legyen .
b) Legyen . Tudjuk, hogy
Megoldás:
a)
A wavelet Fourier trafóját közvetlenül megkaphatjuk a wavelet kiértékelése nélkül:
A táblázatban nincs benne, de közismert, hogy
A táblázatból kiolvasott képletbe behelyettesítve:
b)
Helyettesítésel integrállal tegyük egyszerűbbé a fenti képletet:
Használjuk ki, hogy korábban már kiszámoltuk, hogy
Amit kétszer parciálisan integrálva meg is kapjuk az eredményt:
2) [2016ZH1] A Poisson wavelet a következő:
a) Mutassuk meg, hogy , ha
b) Mutassuk meg, hogy
c)
3) [2016PZH] Legyen . Adjuk meg f által generált wavelet transzformáltjának Fourier transzformáltját!
Numerikus módszerek témakör
Parcdiff egyenletek (Fourier)
1) [2015ZH2] Oldjuk meg Fourier módszerrel az alábbi parciális differenciálegyenletet!
Megoldás:
- Az -t keressük szorzat alakban:
- A diffegyenlet így átírva:
- Ez így már szeparálható (figyeljünk arra, hogy a deriváltak a számlálóban legyenek):
- Nézzük meg, hogy melyik változóra van feltételünk, aminek a jobb oldalán konstans szerepel.
- Az első két féltétel átírva: X(0)T(t) = X(3)T(t) = 0, minden t-re, vagyis X(0) = X(3) = 0
- Tehát az X-re van a T-től nem függő feltételünk, ezért először az X-re oldjuk meg a diffegyenletet!
- Oldjuk meg a diff-egyenletet:
- Írjuk fel a karakterisztikus függvényt!
- Vagyis a diff-egyenlet megoldása:
- Vizsgáljuk meg a kezdeti feltételeket:
Ami csak olyan egész k értékekre teljesülhet, amikre:
- Most oldjuk meg a diff-egyenletet T(t)-re, de a b helyére az újonnan kapott képletet írjuk be.
- A T-re vonatkozó (k-tól függő) diff-egynelet:
- Az -re vonatkozó k-tól függő egyenlet tehát:
- Vezessük be az és konstansokat!
- Az pedig felírható az -k összegeként az összes k-ra.
- A maradék két feltétel segítségével számoljuk ki az és konstansok értékeit.
2) [2016ZH2] Oldjuk meg Fourier módszerrel az alábbi parciális differenciálegyenletet!
Parcdiff egyenletek (véges differenciák)
1) [2015ZH2] Véges differenciák segítségével, felosztás mellett adjuk meg az értékét, ha
2) [2016ZH2] Vázoljuk fel az alábbi feladat megoldását véges differenciák módszerével, ha , az x irányú távolság, h = 1. Mennyi lesz ?
Jordan normál-forma
1) [2016ZH2] Adjuk meg az egyenlet megoldását, ha
Nem lineáris egyenletek numerikus megoldása
1) [2015ZH2] Keressük a egyenlet megoldását. Tudjuk, hogy a gyök a [4, 5] intervallumban van.
a) A gyökhöz milyen közel kell indítani a húrmódszert, hogy az eljárás konvergáljon?
b) Használható-e a [4, 5] intervallumon az iteráció?
2) [2016ZH2] Tekintsük az egyenletet az [1, 2] intervallumon! Megoldható-e iterációval az [1, 2] valamely részintervallumán? Ha igen, milyen rövid legyen? Megoldható-e húrmódszerrel az [1, 2] valamely részintervallumán? Ha igen, milyen rövid legyen?
3) [2016PZH] Az egyenlet esetében az intervallum felezés, vagy az iteráció a célravezetőbb az [1, 2] intervallumon? És a [2, 3]-n?
Lagrange multiplikátor módszer
1) [2015ZH2] Keressük meg az szélsőértékét az feltétel mellett! Vizsgáljuk meg a feltételes definitséget a kapott pontban!
2) [2016ZH2] Hol lehet feltételes szélsőértéke a függvénynek az feltétel mellett? (+3 pontért: Az egyik lehetséges pontban nézzük meg, hogy van-e!)
3) [2016PZH] Hol lehet feltételes szélsőértéke a függvénynek az feltétel mellett? Állapoítsuk meg a szélsőértékek jellegét!
Variáció számítás
1) [2015ZH2] Keressük meg az funkcionálhoz tartozó extremális y függvényt!
2) [2015ZH2] Keressük meg az funkcionálhoz tartozó extremális y függvényt!