„Anal2-magic” változatai közötti eltérés
Ékezetek, könyörgöm |
a ʎ -> λ |
||
121. sor: | 121. sor: | ||
=== Magasabbrendű DE-k === | === Magasabbrendű DE-k === | ||
=== Homogén lineáris, állandó együtthatós DE === | === Homogén lineáris, állandó együtthatós DE === | ||
Megoldás: C * e<sup> | Megoldás: C * e<sup>λ*x</sup> alakban kell keresni. De néha bejön cos/sin is melle, hogy ne legyen egyszerű.<br /> | ||
'''Példa:'''<br /> | '''Példa:'''<br /> | ||
y<sup>(3)</sup> + 2 * y<sup>(2)</sup> + y' = 0<br /> | y<sup>(3)</sup> + 2 * y<sup>(2)</sup> + y' = 0<br /> | ||
λ<sup>3</sup> + 2 * λ<sup>2</sup> + λ = 0<br /> | |||
λ * ( λ<sup>2</sup> + 2 * λ + 1 ) = 0 // annyit emelsz ki amennyi a legkisebb lambda hatványa (itt 1)<br /> | |||
λ * ( λ + 1 )<sup>2</sup> = 0<br /> | |||
első feléből | első feléből λ<sub>1</sub> = 0<br /> | ||
második feléből | második feléből λ<sub>2</sub> = -1 <br /> | ||
DE 3 megoldás kell!!!<br /> | DE 3 megoldás kell!!!<br /> | ||
ilyenkor a homogén megoldáshoz hozzárakunk meg x-el beszorzott tagokat<br /> | ilyenkor a homogén megoldáshoz hozzárakunk meg x-el beszorzott tagokat<br /> | ||
136. sor: | 136. sor: | ||
'''Példa 2:'''<br /> | '''Példa 2:'''<br /> | ||
y<sup>(3)</sup> + 4 * y<sup>(2)</sup> + 13 * y' = 0<br /> | y<sup>(3)</sup> + 4 * y<sup>(2)</sup> + 13 * y' = 0<br /> | ||
λ<sup>3</sup> + 4 * λ<sup>2</sup> + 13 * λ = 0 // kiemelsz λ-t<br /> | |||
λ( λ<sup>2</sup> + 4 * λ + 13 ) = 0<br /> | |||
λ( (λ + 2)<sup>2</sup> + 9 ) = 0<br /> | |||
első feléből | első feléből λ<sub>1</sub> = 0<br /> | ||
második feléből:<br /> | második feléből:<br /> | ||
-9 = ( | -9 = (λ + 2)<sup>2</sup><br /> | ||
-9<sup>1/2</sup> = | -9<sup>1/2</sup> = λ + 2<br /> | ||
-9<sup>1/2</sup> - 2 = | -9<sup>1/2</sup> - 2 = λ<br /> | ||
3*i - 2 = | 3*i - 2 = λ // két darab komplex megoldás lesz, ezek miatt be kell rakni sin/cos-t is a megoldásba<br /> | ||
-3*i - 2 = | -3*i - 2 = λ<br /> | ||
<br /> | <br /> | ||
y<sub>h</sub> = C1 * e<sup>0 * x</sup> + C2 * e<sup>-2 * x</sup> * cos(3 * x) + C3 * e<sup>-2 * x</sup> * sin(3 * x)<br /> | y<sub>h</sub> = C1 * e<sup>0 * x</sup> + C2 * e<sup>-2 * x</sup> * cos(3 * x) + C3 * e<sup>-2 * x</sup> * sin(3 * x)<br /> | ||
tehát a valós rész lesz a | tehát a valós rész lesz a λ, a képzetes rész pedig a cos/sin (pozitív/negatív) belseje<br /> | ||
<br /> | <br /> | ||
'''Példa 3:'''<br /> | '''Példa 3:'''<br /> | ||
adott egy megoldás: 2 * e<sup>5 * x</sup> - e<sup>-3 * x</sup><br /> | adott egy megoldás: 2 * e<sup>5 * x</sup> - e<sup>-3 * x</sup><br /> | ||
ebből kell a DE-et felírni.<br /> | ebből kell a DE-et felírni.<br /> | ||
ebből rögtön latjuk is, hogy | ebből rögtön latjuk is, hogy λ<sub>1</sub> = 5<br /> | ||
λ<sub>2</sub> = -3<br /> | |||
tehát ebből következtethetünk a karakterisztikus egyenletre:<br /> | tehát ebből következtethetünk a karakterisztikus egyenletre:<br /> | ||
( | (λ - 5) * (λ + 3) = 0<br /> | ||
innentől 'csak' át kell rendezni, és megoldani<br /> | innentől 'csak' át kell rendezni, és megoldani<br /> | ||
λ<sup>2</sup> + 3 * λ - 5 * λ - 15 = 0<br /> | |||
λ<sup>2</sup> - 2 * λ - 15 = 0<br /> | |||
y<sup>(2)</sup> - 2 * y - 15 = 0<br /> | y<sup>(2)</sup> - 2 * y - 15 = 0<br /> | ||
<br /> | <br /> | ||
168. sor: | 168. sor: | ||
Ebből kell a homogén DE megoldása.<br /> | Ebből kell a homogén DE megoldása.<br /> | ||
a(x) * y<sup>(2)</sup> + b(x) * y' + c(x) * y = 0<br /> | a(x) * y<sup>(2)</sup> + b(x) * y' + c(x) * y = 0<br /> | ||
y<sub>h</sub> = C1 * y1(x) + C2 * y2(x) // ez ugye az e<sup> | y<sub>h</sub> = C1 * y1(x) + C2 * y2(x) // ez ugye az e<sup>λ*x</sup> -os alak<br /> | ||
Az inhomogén partikuláris megoldása ebből az alábbi alakban lesz:<br /> | Az inhomogén partikuláris megoldása ebből az alábbi alakban lesz:<br /> | ||
itt C1 és C2 ismeretlenek, y1, és (opcionálisan, ha kell) y2, y3 stb. pedig f(x)-ből 'generálódnak' (lásd alább)<br /> | itt C1 és C2 ismeretlenek, y1, és (opcionálisan, ha kell) y2, y3 stb. pedig f(x)-ből 'generálódnak' (lásd alább)<br /> | ||
187. sor: | 187. sor: | ||
'''Konkrét példa:'''<br /> | '''Konkrét példa:'''<br /> | ||
y<sup>(2)</sup> - 5 * y' + 6 * y = 2 * sin(2 * x)<br /> | y<sup>(2)</sup> - 5 * y' + 6 * y = 2 * sin(2 * x)<br /> | ||
λ<sup>2</sup> - 5 * λ + 6 = 0<br /> | |||
λ<sub>1</sub> = 2<br /> | |||
λ<sub>2</sub> = 3<br /> | |||
y<sub>h</sub> = C1 * e<sup>2 * x</sup> + C2 * e<sup>3 * x</sup><br /> | y<sub>h</sub> = C1 * e<sup>2 * x</sup> + C2 * e<sup>3 * x</sup><br /> | ||
y<sub>ip</sub> = A * f(x) + B * f'(x)<br /> | y<sub>ip</sub> = A * f(x) + B * f'(x)<br /> |
A lap 2017. március 8., 23:02-kori változata
Fontos
Ezek a 2 félévnyi Analízis 2 (sima/kereszt) alatt gyűltek össze, többnyire típuspéldákra mennek rá, 2-est (elvileg) simán össze lehet vele szedni.
BTW, a kereszt nem azért jött össze, mert a sima nem ment, hanem mert már nem volt időm tanulni a vizsgára. A gyakorlást NEM helyettesíti. Tehát ezt bemagolod, és utána megoldasz sok zh-t / vizsgát, úgy már jó (elvileg :D ).
Keresztet nem ajánlom :D ua. az anyag, de máshogy kérdezik.
Ha nem mész át ezzel, az a TE hibád :P
A pontosításoknak természetesen mindenki örül
Deriválttábla, számológép nem art :P
Alapok
Azonosságok, amiket jó, ha tudsz
sin2(x) + cos2(x) = 1
cosh2(x) - sinh2(x) = 1
sin(2x) = 2 * sin(x) * cos(x)
sin(x + y) = sin(x) * cos(y) + cos(x) * sin(y)
sin(x ‒ y) = sin(x) * cos(y) ‒ cos(x) * sin(y)
cos(x ‒ y) = cos(x) * cos(y) + sin(x) * sin(y)
cos(x + y) = cos(x) * cos(y) ‒ sin(x) * sin(y)
sinh(a+b) = sinh(a) * cosh(b) + cosh(a) * sinh(b) // sincos-cossin (h)
cosh(a+b) = cosh(a) * cosh(b) + sinh(a) * sinh(b) // coscos-sinsin (h)
limx->0 sin(x) / x = 1 // ezek talán meg anal1-ről :P
limx->0 x / sin(x) = 1
f'(x0) = limx->0 ( f(x) - f(x0) ) / (x - x0)
f'(x0) = limdeltax->0 ( f(x0 + deltax) - f(x0) ) / deltax
cosh(x) = ( ex + e-x ) / 2
sinh(x) = ( ex - e-x ) / 2
Deriválás
f'(c * x) = c * f'(x) // konstanssal szorzás
(f + g)'(x) = f'(x) + g'(x) // összeadás
(f * g)'(x) = f'(x) * g(x) + g'(x) * f(x) // szorzás
(f / g)'(x) = ( f'(x) * g(x) - g'(x) * f(x) ) / g2(x) // osztás
f'( g(x) ) = f'( g(x) ) * g'(x) // összetett fv
(f-1)'(x) = 1 / ( f'( f-1(x) ) ) // inverz fv
Érintő egyenes egyenlete: f(x) = f(x0) + f'(x0) * (x - x0)
Integrálás
ʃ f(x) dx = F(x) + C
ʃ f( fi(x) ) * fi'(x) dx = F( fi(x) ) + C
ʃ fa(x) * f'(x) dx = ( f(x)a + 1 ) / (a + 1) + C // a != -1
ʃ ef(x) * f'(x) dx = ef(x) + C
ʃ f'(x) / f(x) dx = ln| f(x) | + C
ʃ f' * g dx = f * g - ʃ f*g' // parciális integrálás
ʃ (a * x + b) dx = F(a * x + b) / a + C // a != 0
Helyettesítéses integrál:
ʃ f(x) dx // ez vmi bonyolult integrált akar lenni :P
u = f(x) // ez lesz a helyettesítés
du = f'(x) //lederiválod f(x)-et, mert le kell vele osztani
ʃ u / f'(x) du = kijön vmi --> visszahelyettesítesz
Parciális törtekre bontás integrálás
EZT VKI LEÍRHATNÁ IDE
Diffegyenletek (DE)
Elsőrendű DE-k
Szeparábilis DE
y'(x) = g(y) * f(x) // ilyen alakban kell keresni. Ha nincs f(x), akkor beszorzol 1-el, az lesz az f(x) !!!
Meg kell nézni, hogy g(y) mikor lesz 0
g(y) = 0
Megoldod, ha van megoldás, akkor az egy megoldás lesz!
ʃ 1 / g(y) dy = ʃ f(x) dx
ebből kijön: y = K * h(x) // itt a K = eC ; C az integrálás során keletkezik
néha nem kell pont ilyen alakra hozni, azt jelzik.
Lineáris DE
y'(x) + g(x) * y = f(x) // ilyen alakban kell keresni
y'(x) + g(x) * y = 0 --> homogén lineáris DE --> innen szeparábilis, megoldható
y = K * h(x) --> az inhomogén általánoshoz kell K(x) is
K(x) = ʃ f(x) / h(x) dx
//inhomogén általános megoldása
yia = K * h(x) + K(x) * h(x) // homogén + inhomogén partikuláris megoldás
Kezdeti érték probléma: behelyettesítesz, kijön: K = valami
K-t visszahelyettesíted yia-ba --> megkapod: ykonkrét
DE helyettesítéssel
Példán keresztül bemutatva:
y' = 1 / (x + y)
ezt nehéz lenne bármelyik kategóriába besorolni (lineáris, szeparábilis), így valami helyettesítést kell alkalmazni.
Simán megadták, hogy mik lehetnek a helyettesítések, azokból kellett az egyiket alkalmazni.
Lehetséges helyettesítések:
u = x + y
u = y / x
Ehhez a feladathoz az elsőt választjuk. A célunk az, hogy az egyenletben ne legyen csak u alapú változó. Tehát:
u = x + y
kifejezzük y-t:
y = u - x
lederiváljuk:
y' = u' - 1
Tehát most már minden változó y', x+y megvan, behelyettesítünk:
u' - 1 = 1 / u
kicsit rendezzük:
u' = 1 + 1 / u = (u + 1) / u
Ez tehát szeparábilis, g(y) helyett g(u) van, f(x)-et pedig 1 fogja jelképezni.
Megnézzük a 0-re vonatkozó megoldást:
g(u) = (u + 1) / u = 0
u = -1
Tehát visszahelyettesítve: y = -1 - x egy megoldása lesz a DE-nek.
Tovább haladunk a megoldással:
ʃ u / (u + 1) du = ʃ 1 dx
A második fele: x + C
Az első fele:
ʃ u / (u + 1) du = ʃ (u + 1 - 1) / (u + 1) du = ʃ 1 - ( 1 / (u + 1) ) du = u - ln| u + 1 | + C
Ezekből:
u - ln| u + 1 | = x + C --> visszahelyettesítünk
x + y - ln| x + y + 1 | = x + c
Magasabbrendű DE-k
Homogén lineáris, állandó együtthatós DE
Megoldás: C * eλ*x alakban kell keresni. De néha bejön cos/sin is melle, hogy ne legyen egyszerű.
Példa:
y(3) + 2 * y(2) + y' = 0
λ3 + 2 * λ2 + λ = 0
λ * ( λ2 + 2 * λ + 1 ) = 0 // annyit emelsz ki amennyi a legkisebb lambda hatványa (itt 1)
λ * ( λ + 1 )2 = 0
első feléből λ1 = 0
második feléből λ2 = -1
DE 3 megoldás kell!!!
ilyenkor a homogén megoldáshoz hozzárakunk meg x-el beszorzott tagokat
yh = C1 * e0 * x + C2 * e-1 * x + C3 * x * e-1 * x
Példa 2:
y(3) + 4 * y(2) + 13 * y' = 0
λ3 + 4 * λ2 + 13 * λ = 0 // kiemelsz λ-t
λ( λ2 + 4 * λ + 13 ) = 0
λ( (λ + 2)2 + 9 ) = 0
első feléből λ1 = 0
második feléből:
-9 = (λ + 2)2
-91/2 = λ + 2
-91/2 - 2 = λ
3*i - 2 = λ // két darab komplex megoldás lesz, ezek miatt be kell rakni sin/cos-t is a megoldásba
-3*i - 2 = λ
yh = C1 * e0 * x + C2 * e-2 * x * cos(3 * x) + C3 * e-2 * x * sin(3 * x)
tehát a valós rész lesz a λ, a képzetes rész pedig a cos/sin (pozitív/negatív) belseje
Példa 3:
adott egy megoldás: 2 * e5 * x - e-3 * x
ebből kell a DE-et felírni.
ebből rögtön latjuk is, hogy λ1 = 5
λ2 = -3
tehát ebből következtethetünk a karakterisztikus egyenletre:
(λ - 5) * (λ + 3) = 0
innentől 'csak' át kell rendezni, és megoldani
λ2 + 3 * λ - 5 * λ - 15 = 0
λ2 - 2 * λ - 15 = 0
y(2) - 2 * y - 15 = 0
Inhomogén lineáris, állandó együtthatós DE
absztrakt példa:
a(x) * y(2) + b(x) * y' + c(x) * y = f(x)
Ebből kell a homogén DE megoldása.
a(x) * y(2) + b(x) * y' + c(x) * y = 0
yh = C1 * y1(x) + C2 * y2(x) // ez ugye az eλ*x -os alak
Az inhomogén partikuláris megoldása ebből az alábbi alakban lesz:
itt C1 és C2 ismeretlenek, y1, és (opcionálisan, ha kell) y2, y3 stb. pedig f(x)-ből 'generálódnak' (lásd alább)
c * | yip = C1 * y1(x) + C2 * y2(x)
b * | y'ip = C1' * y1(x) + C2 * y2'(x) + C1' * y1(x) + C2 * y2'(x)
// note: C1' * y1(x) + C2 * y2'(x) = 0
a * | y(2)ip = C1' * y1'(x) + C1 * y1(2)(x) + C2' * y2'(x) + C2 * y2(2)(x)
ezt C1, C2-re kell megoldani.
ezután az inhomogén általános megoldás = homogén megoldás + inhomogén partikuláris megoldás
Speciális f(x) esetek:
itt A, Bi ismeretlenek
f(x) = K * ea * x --> yip = A * ea * x
f(x) = amxm+ ... + a0 --> yip = Bmxm+ ... + B0
f(x) = K1 * sin(a * x) --> yip = A * sin(a * x) + B * cos(a * x) // tehát bejön egy cos(a * x) is!
f(x) = K2 * cos(b * x) --> yip = A * cos(b * x) + B * sin(b * x) // tehát bejön egy sin(b * x) is!
Konkrét példa:
y(2) - 5 * y' + 6 * y = 2 * sin(2 * x)
λ2 - 5 * λ + 6 = 0
λ1 = 2
λ2 = 3
yh = C1 * e2 * x + C2 * e3 * x
yip = A * f(x) + B * f'(x)
// annyiszor kell deriválni yip-t, amennyi fokú az eredeti DE is (itt 2)
// ha a homogének között szerepel az yip, akkor külső rezonancia van!
// tehát yip *= x, és utána már lehet deriválni --> ezt kell gyakorolni
// magic: be kell szorozni a deriváltakat az együtthatókkal
6 * | yip = A * sin(2 * x) + B * cos(2 * x)
-5 * | y'ip = 2 * A * cos(2 * x) - 2 * B * sin(2 * x)
1 * | y(2)ip = -4 * A * sin(2 * x) - 4 * B * cos(2 * x)
// magic: Ha megnézed, akkor beszoroztam az elején levő számokkal ott ahol kellett.
sin(2 * x)-ből 2 volt az eredeti DE-ben, tehát:
2 = 6 * A + 5 * 2 * B - 4 * A
cos(2 * x)-ből 0 volt az eredeti DE-ben, tehat:
0 = 6 * B - 5 * 2 * A - 4 * B
ezekből:
A = 1 / 26
B = 5 / 26
yia = C1 * e2 * x + C2 * e3 * x + (1 / 26) * sin(2 * x) + (5 / 26) * cos(2 * x)
Izoklinák
példa:
y' = ey + 2 - x
ebből magic: K = ey + 2 - x
kifejezzük y-t:
y = ln( x + K ) - 2
Ha kérdeznek lokális szélsőértéket, akkor y'-at kell megvizsgálni helyettesítéssel
Az inflexiós ponthoz y(2)-at kell megnézni:
y(2) > 0 --> lokális minimum
y(2) < 0 --> lokális maximum
Ha párhuzamosságot kérdeznek, akkor a meredekség = K-val.
Ezekhez ajánlott megnézni par feladatot, és azon értelmezni :D
Lineáris rekurzió
(ez nagyon magic)
megoldás alakja: f(n) = qn // q != 0
pelda:
f(n) = 4 * f * (n - 1) - 3 * f * (n - 2)
ebből:
qn = 4 * qn - 1 - 3 * qn - 2
a legalacsonyabb hatványú q-val osztunk.
q2 = 4 * q - 3 --> másodfokú
q1 = 1
q2 = 3
ebből:
f(n) = C1 * 1n + C2 * 3n
Ha O(1) típusú megoldások kellenek:
f(n) = O(1): létezik olyan K, hogy |f(n)| <= K * 1, n > N (veges sok kivétel)
tehát: f(n)-nek korlatosnak kell lenni:
C2 = 0
Taylor sorok
// easter egg :D, t.g.o.d.: JrVt10PTD8I
A Taylor sorok arra jók, hogy egy függvényt közelítsünk a deriváltjai segítségével.
Fun fact: ezt régebben arra is használták, hogy a 'drága' sin/cos és hasonló fv-eket helyettesítsek egy 'olcsó' változattal.
f(x) függvény x0 bázispontú n-ed fokú Taylor polinomja:
tehát ahhoz, hogy felírjuk a T-sorát egy függvénynek n db deriváltra lesz szükség.
Analitikus függvény: egy intervallumon analitikus egy függvény, ha ott előállítja a T-sora
Nevezetes függvények T-sorai
Értelmezés sikertelen (formai hiba): {\displaystyle \frac {x^m} {1-x} = \sum_{k=m}^{\infty} x^k \rightarrow Konvergencia tartomány: |x| < 1 }
Lagrange-hiba becsles
Tehát a hibát meg lehet becsülni az n+1-ik T-sor taggal.
xi eleme lesz az [x ; x0] tartománynak, érdemes úgy választani, hogy egyszerű legyen számolni (pl x0 általában jó)
Lagrange-hiba: ( fn + 1(xi) / (n + 1)! ) * (x - xi)n + 1
Példa (keresztről):
y' = sin( y ) + 2 + x
y( x = pi ) = 1
y( x = 3 ) = kb mennyi ? (becsles kell)
felső becsles a hibára?
y'( x = pi ) = sin( 1 ) + 2 + pi // itt az 1 elvileg radiánban van --> számológép!
y(2)( x = pi ) = cos( y ) * y' + 1 = cos( 1 ) * ( sin( 1 ) + 2 + pi ) + 1
T( x0 = pi ) = y( pi ) + y'( pi ) * (x - pi) = 1 + ( sin( 1 ) + 2 + pi ) * (x - pi)
y(3) ~= T( x0 = pi, x = 3 ) = 1 + ( sin( 1 ) + 2 + pi ) * (3 - pi) ~= -0.2 // ezt a tanár nagyon becsülte!
létezik olyan xi, hogy [3 ; pi] tartományban van, mivel felső becslest csinálunk, ezert pi-t választjuk xi-nek.
hiba = | y(3) - T( x0 = pi, x = 3 ) | = Lagrange = ( f(2)(xi) / 2! ) * (3 - pi)2 ~= 0.1 // meg ezt is!
tehát a megoldás: -0.2 +- 0.1
Konvergencia tartomány (KT)
Általában meg van adva vmi T-sor, szummás alakban. Erre alkalmazzuk a hányados / gyökkritériumot.
|an|1/n vagy | (an + 1) / an |
ezután kijön vmi, ami egyenlő 1 / R-el, kifejezzük R-t.
az (x - x0) = 0 egyenletből megkapjuk x-et, ez lesz a KT középpontja.
tehát KT = (x - R, x + R)
végpontokban külön meg kell nézni:
ha divergens --> (
ha konvergens --> [
kell.
Ha x2 van (már nem tudom hol, nézz rá feladatot :D ), akkor u = x2 (helyettesítünk)
a végén meg kell nézni, hogy a KT jó-e.
a <= u=x2 <= b
ez minden x-re teljesül --> |x| < sqrt(a) --> KT: ( -sqrt(a), sqrt(a) )
Fourier-sorok
Megoldás lepései:
- fel kell rajzolni a függvényt
- ha a függvény páros --> bk = 0
- ha a függvény páratlan --> ak = 0, a0 = 0
- fi(x) = a0 / 2 + sum( akcos(k * x) + sin(k * x) )
- ak = 1 / pi * ʃpi-pi f(x) * cos(k * x) dx // itt k = 0, 1, 2, ...
- bk = 1 / pi * ʃpi-pi f(x) * sin(k * x) dx // itt k = 0, 1, 2, ...
- páratlan * páratlan = páros // ezt már nem tudom miért volt a lapon :D --> nézzél feladatot
- ha [-pi ; 0] és [0 ; pi] között ugyanaz a függvény, akkor ez páros függvény lesz
- ekkor elég [0 ; pi] -ig integrálni. A 0 értékű tartományokat ezután ki lehet venni.
- ki kell integrálni a függvényt
- vissza kell helyettesíteni fi(x)-be
- fi(x) = f(x) --> be kell helyettesíteni, a szakadási helyeknél: ( f(x+) + f(x-) ) / 2
Gradiens (aka többváltozós fv-ek deriválása)
Általában adott P0 = (a, b) vektor.
gradf(P0) = f 'x(P0) * i + f 'y(P0) * j = (f 'x, f 'y) // itt i, j egységvektorok
f 'x illetve f 'y úgy jön ki, hogy x illetve y szerint deriválsz.
Pl ha x szerint deriválsz, akkor y csak egy konstans lesz, nem kell vele foglalkozni.
df/de|P0 = gradf(P0) * e // itt e az egységvektor, amit általában megadnak, néha normalizálni kell, a szorzás a két vektor komponensek szerinti szorzása (tehát nem skalár vagy vektor szorzás)
max df/de|P0 = |gradf(P0)| = sqrt((f 'x)2 + (f 'y)2) // azaz a vektor hossza
a maximum iránya: v = gradf(P0) / |gradf(P0)| // normalizálod
Miért létezik gradf? mert a parciális deriváltak f 'x és f 'y léteznek és f(x,y) folytonos P0-ban
Akkor totálisan differenciálható, ha a parciális deriváltak folytonosak P0 pontban, tehát létezik a határértékük // vagy mi :D
f(x,y) P0(x0,y0) érintősík egyenlete: f 'x(x0,y0) * (x - x0) + f 'y(x0,y0) * (y - y0) + f(x0,y0) = z
Megoldási menet kétváltozós szélsőértékszámítás:
1)
Adott f( x,y) kétváltozós függvény
f 'x = ....... = 0
f 'y = ....... = 0
lehetséges szélsőérték
p1(..,..)
p2(..,..)
p3(..,..)
2)
f ′′xx = ...
f ′′xy = ...
f ′′yy = ...
3)
Ha D(p1) =......> 0 akkor szélsőérték
Ha D(p1) =......< 0 akkor nem szélsőérték!
Vagy
D = f ′′xx f ′′yy − (f ′′xy)2 = ...
4)
f xx(p1) =.... ha > 0 akkor min vagy ha < 0 akkor max
Körintegrál
Ebből en két fajtával találkoztam:
- amikor az alakzat egy kor
- amikor az alakzat egy ellipszis
Az integrál alakja általában:
ʃ f(z) / (z - z0)n + 1 dz
A tartomány alakja:
- |z - a| = x // x sugarú, a középpontú kor
- |z - a| + |z - b| = x // ez meg egy ellipszis
Fel kell rajzolni az alakzatot (az ellipszisnek nézz utána!)
Itt négy eset jöhet szoba:
- ha z0 a koron kívülre esik, akkor a függvény reguláris, ʃ f(z) ds = 0
- ha z0 pont a körön van --> nem értelmezett az integrál
- ha z0 a körben van --> ʃ f(z) / (z - z0)n + 1 dz = (2 * pi * i) / n! * f(n)(z0)
- ha több z0 is van (asszem az ellipszis ilyen) akkor kettő integrál összegére kell bontani, majd mind a kettőre |z|-t úgy kell választani, hogy egyszerre csak egy z0-t tartalmazzon a kor (nézz feladatot!), természetesen ilyenkor előbb meg kell nézni az előző három pontot, hogy esetleg nem értelmezzük, kívül esik-e stb. mert akkor nem kell integrálni...
Példa:
ʃ ( f(z) + 1 / (z + 8) ) dz
tartomány: |z - 2 * i| = 2
tehát a középpont = 2 * i
z0 = -8
r = 2
ebből felrajzoljuk a kort, és akkor latjuk, hogy z0 kivul esik a koron, tehát ʃ f(z) dz = 0
Példa 2:
ʃ cos( z ) / ( z4 + 8 * z2 + 16 ) dz
tartomány: |z + 2| + |z - 2| = 5
tehát ez egy ellipszis lesz, több z0 is van.
r = 5
A z0-ok kiszámolása:
z4 + 8 * z2 + 16 = (z2 + 4) * (z2 + 4)
sqrt(z2) = -4
z1 = 2 * i
z2 = -2 * i
felrajzoljuk:
http://i.imgur.com/oon9cwS.png
Ki kell kiszámolni b hosszat, hogy megtudjuk, hogy a z0-ok merre vannak.
Azt tudjuk, hogy A = -2, B = 2.
Tehát (Pitagorasz-tétel, huh?):
b = sqrt( (R / 2)2 - 22 ) = 1.5 // a 2 az A-ból jött, R = 5 ugye
tehát a két z0 kívül esik, emiatt ʃ f(z) dz = 0
Érdemes a többi típusra is nézni feladatot!
Alternatív koordinátarendszerek
Polárkoordináták
Eredetileg ugye egy 2D-s vektor: v = (x, y)
polárban: v = (r, fi)
Átváltás:
x = r * cos( fi )
y = r * sin( fi )
itt r a vektor hossza, fi az x tengellyel bezárt szög
r = sqrt( x2 + y2 )
fi eleme [0 ; 2 * pi]
Jakobi determináns |J|:
|matrix| = r // azaz az alábbi matrix determinánsa
A matrix első oszlopa az r szerinti deriváltakból all, a második pedig a fi szerintiekből. // HF: számold ki ;)
Ha pl egy integrálnál at kell váltani a koordinátarendszert, akkor a függvényt az átváltás után be kell szorozni |J|-vel.
ez a típus hasznos x2 + y2 esetben (amikor ilyesmi van az integrálban)
Hengerkoordináták
ugyanaz mint a polar csak térben, hozzájön z = z is (nem változik)
ez a típus hasznos x2 + y2 + z2 esetben
|J| ugyanaz mint a polárnal.
Gömbkoordináták
ugyanaz mint a henger, csak itt egy gömb felületen van az egész.
átváltás:
x = r * sin( b ) * cos( fi )
y = r * sin( b ) * sin( fi )
z = r * cos( b )
r = sqrt( x2 + y2 + z2 )
fi eleme [0 ; 2 * pi]
b eleme [0 ; pi] // am a b az beta, de mindegy hogy hívod :D
|J| = r2 * sin( b )
A matrix első oszlopa az r szerinti deriváltak (már 3 elem), a második a fi szerinti deriváltak, a harmadik a b szerintiek. // HF: számold ki ;)
Példa:
ʃʃ (2 * x2 + 2 * y2 + 4)7 dT = ?
T: x2 + y2 <= 9, x <= 0, y >= 0
Itt kérdés a tartomány amin integrálni kéne.
Jah és van amikor két alakzat által bezárt területet/térfogatot kérdeznek, ekkor ahhoz, hogy megkapd r-t meg kell nézni, hogy hol metszenek ezek (egyenletmegoldás)
T első részéből megtudjuk, hogy ez egy kor lesz, illetve, hogy r = 3
A második, harmadik részből megtudjuk, hogy ennek a körnek a második negyedének a területe kell.
Tehát fi eleme [pi / 2 ; pi] tartománynak (itt kell majd integrálni)
x = r * cos( fi ) = 3 * cos( fi )
y = r * sin( fi ) = 3 * sin( fi )
Ne felejts el beszorozni |J|-vel!
átváltás után:
ʃʃ r * ( 2 * r2 + 4 )7 dfidr // tartomány: r: [0 ; 3], fi: [pi / 2 ; pi]
ezt már ki tudjuk integrálni, a megoldás: pi / 64 * ( 228 - 48 )
Példa 2:
Térfogatszámítasos integrál.
T: sqrt( x2 + y2 ) <= z <= 6 - ( x2 + y2 )
Ilyenkor az integrált ʃʃ 1 dT-nek kell tekinteni, és itt ki kell találni, hogy hol integráljunk, illetve, hogy mit ( |J| )
T bal és jobb oldalából, ha felrajzoljuk, akkor latszik, hogy ez egy két görbe közötti terület lesz.
Mivel x2 és y2 illetve z van, ezert hengerkoordinátákat fogunk használni. (azért nem gömbit, mert az bonyolultabb)
T polárral: R <= z <= 6 - R2
amint az előző példánál említettem, itt meg kell nézni, hogy hol találkozik a két görbe.
R = 6 - R2 --> másodfokú, R1 = -3, R2 = 2, ebből a -3 nem valószínű, hogy jó.
Tehát az integrál a következő lesz:
ʃʃʃ r dzdrdfi, a tartomány:
z: [r ; 6 - r2 // ezzel nem tudunk mit csinálni, viszont meg így is számolható lesz.
r: [0 ; 2]
fi: [0 ; 2 * pi] // mivel nem mondtak meg, hogy ezzel mi legyen, ezert az alapértelmezett tartományt használjuk.
Innen ez már simán kiintegrálható, nekem 33.51 jött ki.
Példa 3:
Tartománycserés integrál.
ʃʃ (1 + x3)1 / 5 dxdy
T: x: [sqrt(y) / 2 ; 2], y: [0 ; 16]
Ha felrajzoljuk a tartományt, akkor láthatjuk, hogy ez valami vitorla alakú izé lesz.
Mikmakról tanult GTK-s (elfordítod a koordinátarendszert, mert az milyen jó...) módszerrel a tartomány első felénél:
kifejezzük y-t: y = (2 * x)2
Tehát ami történt az az, hogy x(y)-ból áttranszformáltuk y(x)-re (tehát GTK-s ból a normálisra)
Ha ránézünk a rajzra, akkor láthatjuk, hogy x: [0 ; 2] lesz.
Tehát az integrál a következő lesz:
ʃʃ (1 + x3)1 / 5 dydx
T: x: [0 ; 2], y: [0 ; (2 * x)2]
Innen ez is simán kiintegrálható.
Komplex függvénytan
Komplex számok
z = x + i * y // itt x a valós rész, y a képzetes, i = sqrt(-1)
f komplex differenciálható mindenhol --> u,v harmonikus fv-ek --> delta u = u ' 'xx + u ' 'yy = 0
Azonosságok:
|z| = sqrt( x2 + y2 )
/z = x - i * y // konjugált
|z1 * z2| = |z1| * |z2|
|z1 / z2| = |z1| / |z2|
|z|2 = z * /z
|z| = |/z|
arg(z) = -arg(/z) // arg(z): az x tengellyel bezárt szög, trigonometrikusnál fi
/(z1 + z2) = /z1 + /z2
Trigonometrikus alak:
z = r * ( cos(fi) + i * sin(fi) ) // itt r a vektor hossza, fi az x tengellyel bezárt szög
r = |z|
fi = arg(z) // fi: [-pi ; pi]
Exponenciális alak:
z = r * efi * i // úgy lehet megjegyezni, hogy Réffy J. (már akit tanított)
Komplex szorzás:
z1 * z2 = r1 * r2 * ( cos(fi + b) + i * sin(fi + b) ) = r1 * r2 * e(fi + b) * i
Osztás:
z1 / z2 = r1 / r2 * ( cos(fi - b) + i * sin(fi - b) ) = r1 / r2 * e(fi - b) * i
Hatványozás:
zn = rn * ( cos(fi * n) + i * sin(fi * n) )
Gyökvonás:
z1 / n = r1 / n * e( (fi + 2 * k * pi) / n ) * i = r1 / n * ( cos( (fi + 2 * k * pi) / n ) + i * sin( (fi + 2 * k * pi) / n ) )
Euler-formula:
ei * fi = cos(fi) + i * sin(fi) // erre nézz feladatot!
Harmonikus függvények
f(z) = u(x, y) + i * v(x, y) //azaz u a valós rész, v a képzetes rész (függvény)
Azonosságok:
u'x = v'y
u'y = -v'x
u(2)xx = v(2)yx
u(2)yy = -v(2)xy
u(2)xy = v(2)yy
u(2)yx = -v(2)xx
Young tétel: ha egy pont környékén a >= 2 fokú parciális deriváltak léteznek és folytonosak, akkor függetlenek a deriválás sorrendjétől. Tehát xy és yx ugyanaz lesz.
deltau = u(2)xx + u(2)yy
Lokális szélsőértékek:
van, ha f'x = f'y = 0, és // szükséges feltétel
|f(2)xx f(2)xy|
|f(2)yx f(2)yy|
|det| > 0
Ha f(2)xx > 0 --> lokális minimum
Ha f(2)xx < 0 --> lokális maximum
// note: néha a valós részből kell a képzetest kiszámolni. Ilyenkor kiszámolod az elsőfokú deriváltakat, abból ugye megkapod a képzetes elsőfokú deriváltjait, ezt viszont vissza lehet integrálni. --> erre nézz feladatot