Loading [MathJax]/jax/element/mml/optable/GeneralPunctuation.js

Vandermonde

Determinante naše svagdašnje

 

Ivana Grgić,
Sveučilište u Splitu, Fakultet elektrotehnike, strojarstva i brodogradnje, Ruđera Boškovića 32, 21 000 Split, Hrvatska; e-mail address: ipletiko@fesb.hr
Mirjana Strukan
Osnovna škola "Blatine škrape" Split; e-mail address: mirjana.strukan@gmail.com

Sažetak
U ovom radu ćemo opisati i objasniti razliku izmeđ u nekoliko važnih i poznatih matrica i determinanti nazvanih po slavnim matematičarima. Definirani su Jakobijan, Wronskijan i Hesijan, Vandermondeova i Cayley-Mengerova determinanta, njihovo značenje i područje primjene. Cilj rada je na jednom mjestu objediniti osnove tako važnih determinanti koje profesionalni matematičari, a i mnogi studenti, barem u nekom dijelu života, koriste gotovo svakodnevno.




1Uvod

U linearnoj algebri determinante igraju ključnu ulogu. Iz tog razloga su jedan od temeljnih matematičkih pojmova s kojim se susreću studenti na prvim godinama ne samo matematičkog, tehničkog, nego i društvenog usmjerenja. Primjena determinanti je zaista široka. Determinante opisuju prirodu rješenja sustava linearnih jednadžbi, ukazuju na linearnu zavisnost/nezavisnost vektora, utjelovljuju određ ena geometrijska svojstva linearnih transformacija te su najvažniji alat kod pronalaženja svojstvenih vrijednosti matrica. U ovom radu pretpostavljamo da je čitatelju poznat pojam determinante, kao i derivacije funkcija jedne i više varijabli. Cilj rada je na jednom mjestu objediniti osnove važnih determinanti koje profesionalni matematičari, a i mnogi studenti, barem u nekom dijelu života, koriste gotovo svakodnevno.

2Gradijent

Prije nego opišemo Jakobijan i Hessijan, zgodno je uvesti pojam gradijenta. Gradijent je generalizacija pojma derivacije na skalarne funkcije više varijabla. često umjesto riječi funkcija koristimo riječ polje. Pojam polja je uveo irski matematičar Hamilton, u svezi s promatranjem fizikalnih veličina u električnim, magnetnim i drugim poljima. Skalarna i vektorska polja su u stvari drugi naziv za skalarne i vektorske funkcije.


Definition 1. Neka je f:DR,DR3 skalarno polje. Gradijent skalarnog polja f, u oznaci gradf, je vektorsko polje koje definiramo na sljedeći način gradf=fxi+fyj+fzk.

Definition 2. Hamiltonov diferencijalni operator (nabla) glasi =ix+jy+kz.

Diferencijalni operator istovremeno ima svojstva i vektora i derivacije. Vrijedi gradf=f. Gradijent možemo smatrati retčanom matricom, gradf=[fxfyfz].

Gradijent označava smjer najbržeg rasta funkcije f u točki T0(x0,y0,z0). Iznos tog rasta je |gradf(T0)|. Primjerice, ako skalarno polje f opisuje temperaturu u proizvoljnoj točki prostora, onda će gradijent od f u točki T(x,y,z) pokazivati u smjeru u kojem temperatura najbrže raste.

3Jakobijan

U literaturi, termin Jakobijan1 se koristi naizmjenično za Jacobijevu matricu i determinantu. I matrica i determinanta imaju korisnu i važnu primjenu: matrica sadrži parcijalne derivacije prvog reda, a determinanta je korisna u procesu promjene varijabla u integralnom računu.


Neka je F:RnRm vektorsko polje. Jakobijan možemo smatrati generalizacijom derivacije na vektorska polja. Potrebna nam je brzina promjene svake komponente funkcije F za ulaznu varijablu x=(x1,x2,...,xn), a upravo to je uhvaćeno u matrici koju nazivamo Jacobijeva matrica J.
Dakle,

J=[f1x1f1x2f1xnf2x1f2x2f2xnfmx1fmx2fmxn].
Jakobijan je determinanta Jacobijeve (kvadratne) matrice. U slučaju kada je m=1 i n=3 Jacobijeva matrica je isto što i gradijent (retčana matrica). Jedna od važnih primjena Jacobijeve matrice, odnosno njene determinante je u računanju dvostrukih integrala, kad nam je potrebna zamjena varijabli kojom transformiramo originalnu funkciju u neku mnogo jednostavniju koju ćemo lakše integrirati. Primjerice, često pri računanju koristimo vezu izmeđ u kartezijevih i polarnih koordinata: x(r,φ)=rcosφ, y(r,φ)=rsinφ. Općenito, pojam Jakobijana za funkcije x=x(u,v), y=y(u,v) definiramo jednakošću detJ=|xuxvyuyv|. Dakle, Jakobijan za promjenu koordinata u polarne je jednak detJ=|xrxφyryφ|=|cosφrsinφsinφrcosφ|=r.
Formula za prijelaz iz kartezijevih u polarne koordinate u dvostrukom integralu glasi Df(x,y)dxdy=Df(x(r,φ),y(r,φ))|detJ|drdφ=Df(rcosφ,rsinφ)rdrdφ.


Example 3. Izračunati površinu lika D u xy ravnini omeđ enog kružnicama x2+y2=1, x2+y2=4 i pravcima y=33x, y=x za x0.


Kako je područje integracije omeđ eno kružnicama prikladno je prijeći u polarne koordinate. Jakobijan je tada jednak detJ=r.


Rubne krivulje i pravci područja D u polarnim koordinatama imaju jednadžbe r1=1, r2=2, tanφ1=1, tanφ2=33. Tim jednadžbama je određ eno područje D u rφ ravnini. Područje D je slika područja D u odnosu na zadanu transformaciju. Računamo: P=Ddxdy=Drdrdφ=π4π6dφ21rdr=π4π6r22|21dφ=π4π6(212)dφ=32φ|π4π6=32π12=π8.


Apsolutna vrijednost Jakobijana |detJ| u slučaju kada je detJ0 je faktor kojim množimo površinu pravokutnika da bismo dobili površinu krivocrtnog lika. Jakobijan djeluje kao faktor skaliranja izmeđ u dva koordinatna sustava.

4Hessijan

Vidjeli smo da je gradijent retčana matrica koja sadrži prve derivacije funkcije više varijabli. Hesseova2 matrica je matrica parcijalnih derivacija drugog reda funkcije više varijabli.

Definition 4. Hesseova matrica je kvadratna matrica dimenzija n×n parcijalnih derivacija drugog reda funkcije f:RnR tj. funkcije koja n-dimenzionalni vektor preslikava u skalar. Element Hesseove matrice definiran je s Hi,j=2fxixj, a sama matrica ovako:

H=f(x)=[2fx212fx1x22fx1xn2fx2x12fx222fx2xn2fxnx12fxnx22fx2n].


Determinantu Hesseove matrice nazivamo Hessijan.

Zbog Schwarzovog teorema, za funkcije s neprekidnim parcijalnim derivacijama drugog reda, Hesseova matrica je simetrična. Takođ er, uočimo vezu Hesseove matrice i Jacobijeve matrice. Vrijedi H(f(x))=J(f(x)), tj. Hesseova matrica funkcije f je Jacobijeva matrica gradijenta funkcije f.


Najčešća primjena Hessijana je pri određ ivanju ekstrema funkcije više varijabli.

Definition 5. Za kvadratnu matricu AMn njene glavne minore su determinante kvadratnih matrica A1,A2,,An1,An=A, gdje je Ak matrica koja se iz A dobije tako da uzmemo njenih prvih k redaka i k stupaca, tj. Ai=[aij]Mk.
 


Funkcija f u stacionarnoj točki x0 ima lokalni minimum ako su sve minore Hesseove matrice H(f)(x0) pozitivne, a lokalni maksimum ako predznaci minora Hesseove matrice H(f)(x0) alterniraju počevši s negativnim.

Ukoliko su neke od minora nula, ali nema negativnih ili pak alterniraju tako da neparne po redu nisu pozitivne, a parne nisu negativne, potrebne su druge metode provjere. U preostalim slučajevima riječ je o sedlastoj točki, tj. stacionarnoj točki koja nije ekstrem.



Example 6. Odrediti lokalne ekstreme funkcije zadane formulom f(x,y)=4xyx4y4.

Gradijent funkcije f jednak je f(x,y)=(4y4x3)i+(4x4y3)j,
a Hesseova matrica H(f)(x,y)=[12x24412y2].
Stacionarne točke funkcije f su (0,0), (1,1) i (1,1).

Računamo: 2fx2(0,0)=0, 2fxy(0,0)=4, 2fy2(0,0)=0.

Sada je detH(f)(0,0)=16 pa je (0,0) sedlasta točka.


Zatim: 2fx2(1,1)=12, 2fxy(1,1)=4, 2fy2(1,1)=12.

Vrijedi detH(f)(1,1)=128 pa je (1,1) lokalni maksimum i vrijednost lokalnog maksimuma u (1,1) iznosi f(1,1)=2.


I konačno: 2fx2(1,1)=12, 2fxy(1,1)=4, 2fy2(1,1)=12

te je detH(f)(1,1)=128 pa je i (1,1) lokalni maksimum i vrijednost lokalnog maksimuma u (1,1) opet iznosi f(1,1)=2.

 

5Wronskijan

Do sada smo upoznali Jacobijevu matricu i vidjeli njenu primjenu u rješavanju zadataka s dvostrukim integralima te Hesseovu matricu primijenjenu u računanju ekstrema funkcija više varijabli. U ovom poglavlju kratko ćemo opisati determinantu matrice Wronskog3 koju nazivamo Wronskijan te pokazati njenu primjenu u rješavanju diferencijalnih jednadžbi drugog reda.
Iz linearne algebre znamo definiciju linearne nezavisnosti vektora.

Definition 7. Vektori a1,a2,,ak su linearno nezavisni ako za sve skalare λ1,,λkR λ1a1+λ2a2++λkak=0λ1==λk=0.
U protivnom su vektori a1,,ak linearno zavisni.

Slično želimo definirati i linearnu nezavisnost derivabilnih funkcija jedne varijable.

Definition 8. Dvije funkcije y1 i y2 su linearno nezavisne na intervalu IR ako za skalare k1,k2R k1y1+k2y2=0k1=0,k2=0. U protivnom su funkcije y1,y2 linearno zavisne.

Iz definicije vidimo da kada je k10 ili k20 možemo dobiti y1=k2k1y2 ili y2=k1k2y1, odnosno y1=my2 ili y2=ny1 (y1 i y2 su proporcionalne funkcije).



Example 9. Funkcije f(x)=2sin2x i g(x)=1cos2x su linearno zavisne jer (1)(2sin2x)+(2)(1cos2x)=0.



Osim ovoga, postoji i sistematičniji način ispitivanja linearne nezavisnosti funkcija. O tome govori sljedeća definicija i teorem.

Definition 10. Neka su y1,y2:IR derivabilne funkcije. Funkcija W(x)=|y1(x)y2(x)y1(x)y2(x)|=y1(x)y2(x)y1(x)y2(x)
je determinanta Wronskog ili Wronskijan funkcija y1 i y2 .
 

Theorem 11. Ako su funkcije y1 i y2 linearno zavisne na intervalu I , onda je njihov Wronskijan identično jednak nula.4



Example 12. Pokazati da su funkcije f(x)=ex i g(x)=xex linearno nezavisne.


Izračunajmo Wronskijan W(x)=|exxexex(x+1)ex|=(x+1)e2xxe2x=e2x0. Kako Wronskijan nije nula, zaključujemo da su ex i xex linearno nezavisne funkcije na svakom intervalu.


Vidjeli smo definiciju Wronskijana u slučaju n=2. Za općeniti n, Wronskijan je determinanta W(y1,y2,...,yn)=|y1yny1ynyn11yn1n|.

Studenti matematičkih i tehničkih fakulteta Wronskijan najčešće susreću pri rješavanju homogenih linearnih diferencijalnih jednadžbi (DJ) drugog (ili višeg) reda. Homogene linearne DJ drugog reda su DJ oblika y. Znamo da ako su y_{1} i y_{2} dva rješenja takve jednadžbe, onda je i svaka njihova linearna kombinacija c_{1}y_{1}+c_{2}y_{2} takođ er rješenje te jednadžbe. No ako su funkcije y_{1} i y_{2} linearno nezavisne funkcije na nekom intervalu \mathcal{I} onda je opće rješenje jednadžbe dano s c_{1}y_{1}+c_{2}y_{2}.


Example 13. Riješiti DJ y''-y=0, y(0)=4, y'(0)=-2.


Funkcije y_{1}=e^{x} i y_{2}=e^{-x} su rješenja ove homogene linearne DJ drugog reda za sve x\in\mathbb{R}. Naime, direktnom provjerom za y_{1}=e^{x} dobijemo (e^{x})''-e^{x}=e^{x}-e^{x}=0, a slično i za y_{2}=e^{-x}. Njihova linearna kombinacija y=c_{1}e^{x}+c_{2}e^{-x} je takođ er rješenje DJ. Iz početnih uvjeta dobijemo: y(0)=c_{1}+c_{2}=4 \quad y'(0)=c_{1}-c_{2}=-2 . Rješenje ovog sustava je c_{1}=1, c_{2}=3. Time smo dobili rješenje DJ koje glasi y=e^{x}+3e^{-x} koje je ujedno i opće rješenje jer je W(x)=\begin{vmatrix} e^{x} & 3e^{-x} \\ e^{x} & -3e^{-x} \end{vmatrix}=-3-3=-6\neq 0. Primjedba. Da smo umjesto e^{x} i e^{-x} uzeli funkcije y_{1}=e^{x} i y_{2}=ke^{x}, \,k\in\mathbb{R}, te stavili y=c_{1}e^{x}+kc_{2}e^{x} naše rješenje ne bi bilo opće jer su ovakvi y_{1} i y_{2} linearno zavisne funkcije (proporcionalne su).

6Vandermondeova determinanta

Vandermondeova5 determinanta je jedna od najpoznatijih eksplicitnih formula za neku determinantu u matematici. Retci Vandermondeove matrice su geometrijski nizovi. Neka su x_{0},x_{1},...,x_{n} realni brojevi, x=( x_{0},x_{1},...,x_{n}). Matricu u oznaci V_{n}(x) definiranu na način V_{n}(x)=\begin{bmatrix} 1 & x_{0} & x_{0}^{2} & \cdots & x_{0}^{n} \\ 1 & x_{1} & x_{1}^{2} & \cdots & x_{1}^{n} \\ \vdots & \vdots & \vdots & \vdots & \vdots \\ 1 & x_{n} & x_{n}^{2} & \cdots & x_{n}^{n} \end{bmatrix} nazivamo Vandermondeova matrica. Determinanta kvadratne matrice V_{n}(x) je Vandermondeova determinanta. Njena vrijednost je polinom \det(V_{n}(x))=\prod_{1\leq i\lt j\leq n} (x_{j}-x_{i}). Jedna od primjena Vandermondeove determinante je u određ ivanju interpolacijskog polinoma p_{n}(x), n-tog stupnja koji će dovoljno dobro aproksimirati funkciju f za koju je poznato njeno djelovanje na konačnom skupu točaka, tj. znamo da je f(x_{i})=y_{i},\, i=0,1,...,n.

Example 14. Odrediti interpolacijski polinom čiji graf prolazi točkama (1,2), (5,3), (-2,0), (7,-1).

Ako p_{3}(x) prikažemo u kanonskom obliku p_{3}(x)=a_{0}+a_{1}x+a_{2}x^{2}+a_{3}x^{3}, tada treba odrediti koeficijente a_{0}, a_{1}, a_{2}, a_{3} tako da vrijede uvjeti interpolacije p_{3}(x_{i})=y_{i}, \, i=0,1,2,3, tj. \begin{align*} a_{0}+a_{1}\cdot 1+a_{2}\cdot 1^{2}+a_{3}\cdot 1^{3} &= 2 \\ a_{0}+a_{1}\cdot 5+a_{2}\cdot 5^{2}+a_{3}\cdot 5^{3} &= 3 \\ a_{0}+a_{1}\cdot (-2)+a_{2}\cdot (-2)^{2}+a_{3}\cdot (-2)^{3} &= 0 \\ a_{0}+a_{1}\cdot 7+a_{2}\cdot 7^{2}+a_{3}\cdot 7^{3} &= -1 . \end{align*}
U matričnom obliku sustav možemo zapisati:

\underbrace{\begin{bmatrix} 1 & 1 & 1 & 1 \\ 1 & 5 & 5^{2} & 5^{3} \\ 1 & -2 & (-2)^{2} & (-2)^{3} \\ 1 & 7 & 7^{2} & 7^{3} \end{bmatrix}}_{V} \cdot \underbrace{\begin{bmatrix} a_{0} \\ a_{1} \\a_{2} \\ a_{3} \end{bmatrix}}_{\vec{a}}=\underbrace{\begin{bmatrix} 2 \\ 3 \\ 0 \\ -1 \end{bmatrix}.}_{\vec{y}}
Cramerovim pravilom riješimo sustav,

a_{i}=\frac{D_{i}}{\det V}, \, i=0,1,2,3,
gdje je \det V Vandermondeova determinanta, a D_{i} determinanta matrice koju dobijemo tako da se u matrici V(i+1)-vi stupac zamijeni vektorom \vec{y}. Odredimo Vandermondeovu determinantu

\det V= \prod_{1\leq i\lt j\leq 3} (x_{j}-x_{i}) = (7-1)(7-5)(7-(-2))(-2-1)(-2-5)(5-1)=9072.
Odredimo vrijednosti determinanti D_{i}, te zapišemo rješenje sustava

a_{0}=\frac{D_{0}}{\det V}=\frac{9996}{9072}, \,a_{1}=\frac{D_{1}}{\det V}=\frac{7734}{9072}, \, a_{2}=\frac{D_{2}}{\det V}=\frac{732}{9072}, \, a_{3}=\frac{D_{3}}{\det V}=-\frac{318}{9072}.
Sada možemo zapisati interpolacijski polinom

p_{3}(x)=\frac{1}{9072}\left(9996+7734x+732x^{2}-318x^{3}\right).



7Cayley-Mengerova determinanta

Studentima vjerojatno manje poznata determinanta od prethodno opisanih je Cayley-Mengerova6 determinanta koja se koristi u linearnoj algebri, geometriji i trigonometriji. Predstavlja formulu za sadržaj (duljinu/površinu/volumen) n-dimenzionalnog simpleksa izraženu preko kvadrata udaljenosti svih parova vrhova simpleksa. U geometriji, simpleks je najjednostavniji politop koji se može formirati, odnosno generalizacija pojma trokuta ili tetraedra na proizvoljne dimenzije. Na primjer 0-dimenzionalni simpleks je točka, 1-dimenzionalni simpleks je segment (dio pravca), 2-dimenzionalni simpleks je trokut, a 3-dimenzionalni tetraedar.

Definition 15. Neka su A_{1},A_{2},...,A_{n+1} vrhovi n-dimenzionalnog simpleksa u n-dimenzionalnom Euklidskom prostoru. Neka je d_{ij} euklidska udaljenost izmeđ u vrhova A_{i} i A_{j}, d_{ij}=\Vert A_{i}-A_{j}\Vert _{2}. Tada se n-dimenzionalni sadržaj simpleksa v_{n} može izračunati iz sljedeće formule: v^{2}_{n}=\frac{(-1)^{n+1}}{(n!)^{2} 2^{n}} \begin{vmatrix} 0 & d_{12}^{2} & \cdots & d_{1(n+1)}^{2} & 1 \\ d_{12}^{2} & 0 & \cdots & d_{2(n+1)}^{2} & 1 \\ d_{13}^{2} & d_{23}^{2} & \cdots & d_{3(n+1)}^{2} & 1 \\ \vdots & \vdots & & \vdots & \vdots \\ d_{1(n+1)}^{2} & d_{2(n+1)}^{2} & \cdots & 0 & 1 \\ 1 & 1 & \cdots & 1 & 0 \end{vmatrix}.


Determinanta u gornjoj definiciji je Cayley-Mengerova determinanta \Delta. Vezu izmeđ u Cayley-Mengerove determinante i sadržaja n-dimenzionalnog simpleksa lakše uočavamo ako raspišemo determinantu za n=1, n=2 i n=3.
Kako smo već napisali 1-dimenzionalni simpleks je dio pravca izmeđ u dvije točke A_{1} i A_{2}. Tada je \Delta =\begin{vmatrix} 0 & d_{12}^{2} & 1 \\ d_{12}^{2} & 0 & 1 \\ 1 & 1 & 0 \end{vmatrix}=2d_{12}^{2}=2v^{2}_{1},
odnosno Cayley-Mengerova determinanta \Delta je jednaka dvostrukom kvadratu udaljenosti izmeđ u dviju točaka.

Ako za n=2 označimo d_{12}=a, d_{13}=b i d_{23}=c onda je \Delta =\begin{vmatrix} 0 & a^{2} & b^{2} & 1 \\ a^{2} & 0 & c^{2} & 1 \\ b^{2} & c^{2} & 0 & 1 \\ 1 & 1 & 1 & 0 \end{vmatrix}=a^{4}-2a^{2}b^{2}+b^{2}-2a^{2}c^{2}-2b^{2}c^{2}+c^{4}.
Stavimo li s=\frac{1}{2}(a+b+c) vrijedi \Delta=-16s(s-a)(s-b)(s-c).
Iz Heronove formule znamo da je površina trokuta sa stranicama a, b i c jednaka \displaystyle\sqrt{s(s-a)(s-b)(s-c)}. Dakle, za n=2 Cayley-Mengerova determinanta je proporcionalna kvadratu poršine trokuta; \displaystyle\Delta=-16v^{2}_{2}. Za n=3 sadržaj 3-simpleksa (tj. volumen tetraedra) je dan s

v^{2}_{3}=\frac{1}{288}\begin{vmatrix} 0 & d_{12}^{2} & d_{13}^{2} & d_{14}^{2}& 1 \\ d_{12}^{2} & 0 & d_{23}^{2} & d_{24}^{2} & 1 \\ d_{13}^{2} & d_{23}^{2} & 0 & d_{34}^{2} & 1 \\ d_{14}^{2} & d_{24}^{2} & d_{34}^{2} & 0 & 1 \\ 1 & 1 & 1 & 1 & 0 \end{vmatrix}.
U ovom radu smo opisali matrice i determinante s kojima se najčešće susreću studenti prve i druge godine i predavači matematike na fakultetima. Bez njih je teško i zamisliti bilo koje fakultetsko gradivo. Osim navedenih, brojne su druge zanimljive determinante kao što su Toeplitzova, Sylvesterova, Slaterova, Diudonneova, Hankelova, Pascalova, Cauchyjeva, ali one izlaze iz okvira ovog rada.
Bibliografija
[1] E. Kreyszig: Advanced engineering mathematics, 8th edition, 1999.
[2] T. Burić, L. Korkut, J.P. Milišić, M. Pašić, I. Velčić: Vektorska analiza, Element, Zagreb, 2014.
[3] I. Brnetić, V. županović: Višestruki integrali, Element, Zagreb, 2014.
[4] https://math.libretexts.org/Bookshelves/Analysis/Supplemental
[5] http://www.mathematics.digital/matematika1/index.html
[6] http://lavica.fesb.unist.hr/matematika2/
[7] https://najeebkhan.github.io/blog/VecCal.html
[8] S. Kurepa: Matematička analiza III, Tehnička knjiga, Zagreb, 1975.
[9] https://www.pmf.unizg.hr/_download/repository/mat2-pred9-novo.pdf
[10] https://en.wikipedia.org/wiki/Cayley-Menger_determinant