Preview only show first 10 pages with watermark. For full document please download

P07 [246,46 Kib]

   EMBED


Share

Transcript

verzija 0.1 – 23.11.2015. 7. Fourierov red po ortogonalnim sustavima U ovom c´emo poglavlju promatrati op´cenite vektorske prostore u kojima je definiran skalarni umnoˇzak. Postojanje skalarnog umnoˇska omogu´cava promatranje geometrijskih pojmova: okomitosti, duljine vektora, udaljenosti dvaju vektora na prirodan naˇcin, poop´cenjem tih pojmova iz trodimenzionalnog prostora. Cilj nam je izvesti op´cenite postupke za prikaz vektora preko elemenata ortogonalne baze i primjeniti do bivene formule na prostor funkcija. 1. Fourierovi redovi u unitarnom prostoru 1.1. Skalarni produkt. Unitarni prostori. Neka je X vektorski prostor, x , y , z ,. . . njegovi elementi. Skalarni umnoˇzak Skalarni umnoˇzak na vektorskom prostoru X je svaka funkcija (· | ·) : X × X → R koja ima sljede´ca cˇetiri svojstva: 1) Pozitivnost. Za svaki x ∈ X vrijedi (x | x) > 0 . Nadalje, (x | x) = 0 ako i samo ako je x = 0 . 2) Homogenost. Za svaki skalar α ∈ R vrijedi (α x | y) = α (x | y). 3) Komutativnost. Za svake x, y ∈ X vrijedi (x | y) = (y | x). 4) Aditivnost. Za svaka tri vektora x, y, z ∈ X vrijedi (x + y | z) = (x | z) + (y | z). Vektorski prostor X na kojem postoji skalarni umnoˇzak zovemo unitarni prostor. 2 7. FOURIEROV RED PO ORTOGONALNIM SUSTAVIMA Iz skalarnog umnoˇska izvodimo i druge vaˇzne geometrijske pojmove. Tu je najprije duljina vektora (norma): p kxk := (x | x), te udaljenost dvaju vektora: d(x, y) := kx − yk. Pojam ortogonalnosti bit c´e u nastavku kljuˇcan. Za vektore x i y kaˇzemo da su okomiti (ortogonalni) ako vrijedi (x | y) = 0 . 1.2. Primjer. Skalarni umnoˇzak u Rn . Elementi prostora Rn uredene su n torke: x = (x1 , . . . , xn ), y = (y1 , . . . , yn ). Skalarni umnoˇzak vektora x i y definiramo formulom (x | y) := x1 y1 + x2 y2 + . . . + xn yn . (7.1) Norma vektora raˇcuna se formulom q p kxk = (x | x) = x21 + . . . + x2n . (7.2) Vektori x i y su okomiti ako vrijedi x1 y1 + x2 y2 + . . . + xn yn = 0. (7.3) 1.3. Prostor L2 (a, b) . Neka je I = (a, b) interval u R , pri cˇemu moˇze biti a = −∞ i b = +∞ . S L2 (a, b) oznaˇcavamo prostor svih kvadratno integrabilnih funkcija. Funkcija f definirana na intervalu [a, b] pripada tom prostoru ako vrijedi Z b |f (x)|2 dx < ∞. a Ukoliko za neku funkciju f vrijedi Z b |f (x)|2 dx = 0 a tada je ona jednaka nuli svuda, osim moˇzda u nekoliko toˇcaka. Takvu funkciju c´emo u okviru ove teorije poistovjetiti s nul-fukcijom, zato sˇ to nam ponaˇsanje funkcije u jednoj ili nekolicini toˇcaka ne´ce biti vaˇzno. U tom smislu, ako za dvije funkcije f i g vrijedi Z b |f (x) − g(x)|2 dx = 0, a smatrat c´emo da su te dvije funkcije jednake. 7. FOURIEROV 3 RED PO ORTOGONALNIM SUSTAVIMA U prostoru L2 (a, b) definiran je skalarni umnoˇzak formulom (f | g) = Z b f (x)g(x) dx. (7.4) a Lako je provjeriti da ovako definiran skalarni produkt zaista ima cˇetiri traˇzena svojstva. Zato je L2 (a, b; ρ) unitarni prostor. Norma elementa prostora L2 (a, b) glasi kf k = (f | f )1/2 = Z a b 1/2 |f (x)|2 dx . Udaljenost dviju funkcija u ovom prostoru raˇcunamo na naˇcin d(f , g) = kf − gk2 = Z a b 1/2 |f (x) − g(x)| dx . 2 ˇ Udaljenost d jest prikladna mjera za udaljenost funkcija u mnogim primjenama. Cinjenica da je ta udaljenost inducirana skalarnim produktom omogu´cava da se koriste sva svojstva unitarnih prostora u prouˇcavanju prostora funkcija opskrbljenog ovom udaljenoˇsc´u. Postojanje skalarnog produkta omogu´cava razmatranje pojmova kao sˇ to su ortogonalnost i ortogonalna projekcija. Mjerenje udaljenosti na prostoru funkcija pomo´cu skalarnog produkta ima odredene prednosti u odnosu na udaljenost koja promatra razlike funkcijskih vrijednosti. 1.4. Primjer. Promotrimo funkcije f n (x) = xn definirane na intervalu (0, 1) . Sa slike nam je vidljivo da porastom broja n ove funkcije se “pribliˇzavaju” funkciji f (x) = 0, 0 < x < 1. Medutim, ako udaljenost dviju funkcija definiramo na uobiˇcajeni naˇcin, preko supremuma razlika njihovih funkcijskih vrijednosti, za njihovu udaljenost c´e vrijediti za svaki n d∞ (f n , f ) = sup |f n (x) − f (x)| = sup |xn | = 1. 0 kxk2 − n X k=1 n X k=1 n X k=1 n X k=1  x− n X k=1  n X ak ek ak ek x − k=1  X n n n X  X aj ej ak ek ) + ak ek ak ek x − (x | k=1 ak (ek | x) + ak ck kek k2 + n X k=1 k=1 (ej | ek ) j,k=1 n X k=1 a2k kek k2 n X (ak − ck )2 kek k2 c2k kek k2 + k=1 c2k kek k2 . Moˇzemo zakljuˇciti sljede´ce: 1) Ova je udaljenost P najmanja kad je ak = ck . 2) Uvijek vrijedi n1 c2k kek k2 6 kxk2 . Time smo pokazali sljede´ci teorem. Svojstvo najbolje aproksimacije Fourierovog reda Teorem 1.4. Fourierov red (??) ima svojstvo da najbolje aproksimira element x , za svaki n je udaljenost n X ak ek || kx − k=1 (x | ek ) . kek k2 Nadalje, za koeficijente ck Fourierovog reda vrijedi Besselova nejednakost najmanja, ukoliko je ak = ck = ∞ X k=1 c2k kek k2 6 kxk2 . 1.15. Uvjeti za ortogonalnu bazu. P Niz sin x , sin 2x ,. . . , sin nx je ortogonalan. Medutim, kako je bk sin kx uvijek neparna funkcija, niti jedna drukˇ c ija funkcija ne moˇ z e se napisati u obliku P f (x) = bk sin kx . Ovaj je niz presiromaˇsan da bi cˇinio bazu u prostoru L2 (−π , π ) . 7. FOURIEROV RED PO ORTOGONALNIM SUSTAVIMA 11 Postavlja se pitanje: kad neki ortogonalni skup cˇini bazu? S tim u vezi bit c´e sljede´ca definicija: Svojstva ortogonalnog niza Neka je (en ) ortogonalan niz u unitarnom prostoru X . Za njega kaˇzemo da je • zatvoren, ako za koeficijente Fourierovog reda vrijedi Parsevalova jednakost n X c2k kek k2 = kf k2 ; k=1 • potpun, ako je svaki element koji je okomit na sve svaki ek identiˇcki jednak nuli. Teorem 1.5. Neka je (ek ) ortogonalan sistem u L2 (a, b) . Slijede´ce su tvrdnje ekvivalentne 1) (ek ) je baza. 2) (ek ) je zatvoren. 3) (ek ) je potpun. 1.16. Ortogonalni sustavi. Istaknutu ulogu u prostoru L2 (a, b, ρ) cˇine polinomi. Sistem polinoma 1 , x , x2 , . . . linearno su nezavisni, jer identitet λ0 + λ1 x + λ2 x2 + . . . λn xn = 0 moˇze vrijediti samo onda ako su svi skalari jednaki nuli. Konaˇcan niz polinoma 1 , x ,. . . , xn razapinje prostor Pn svih polinoma stupnja n . Naime svaki se polinom P stupnja n moˇze prikazati u obliku linearne kombinacije elemenata ove baze. Ta baza nije ortogonalna. Npr., vrijedi (x | x3 ) = Z a b x · x3 ρ(x)dx > 0. Pokazat c´emo kako se Gramm-Schmidtovim postupkom ortogonalizacije ta baza moˇze zamijeniti sistemom polinoma Q0 , Q1 ,. . . , Qn . . . koji su ortogonalni na intervalu [a, b] , s teˇzinskom funkcijom ρ . Pri tom je Qk polinom stupnja k . Neka je Pk bilo koji polinom stupnja k . On se moˇze prikazati u obliku linearne kombinacije polinoma Q0 , Q1 , . . . , Qk . Kako je Qn okomit na sve polinome Qj , j < n , zakljuˇcujemo da c´e Qn biti okomit na po volji odabran polinom Pk , cˇim je k < n . Dakle, Qn je okomit na prostor Pn−1 . 12 7. FOURIEROV RED PO ORTOGONALNIM SUSTAVIMA 1.17. Nul-toˇcke ortogonalnih polinoma. U mnogim primjenama od osnovne s vaˇznosti nul toˇcke ortogonalnih polinoma. U sljede´cem teoremu c´emo dokazati neke osnovne cˇinjenice o karakteru tih nul-toˇcaka. Teorem 1.6. Neka je n > 1 . Sve nul-toˇcke polinoma Qn su realne razliˇcite i leˇze unutar intervala [a, b] . DOKAZ. Dokaˇzimo najprije da realne nul toˇcke ne mogu biti viˇsestruke. Pretpostavimo obratno: x0 je viˇsestruka realna nul-toˇcka. Tada je Rn−2 (x) := Qn (x) (x − x0 )2 polinom stupnja n − 2 . On je stoga okomit na Qn . Medutim, vrijedi Z b Qn (x)2 ρ(x)dx > 0 (Qn | Rn−2 ) = 2 a (x − x0 ) jer je podintegralna funkcija pozitivna. Time smo dobili kontradikciju. Stoga su sve realne nul-toˇcke jednostruke. Na sliˇcan naˇcin c´emo pokazati da realnih nul-toˇcaka ima n te da sve leˇze unutar intervala [a, b] . kad bi ih bilo manje, recimo x1 , . . . , xk , ( k < n ), tada bismo mogli napisati Qn (x) = (x − x1 ) · · · (x − xk )Rn−k (x). Tu je Rn−k polinom stupnja n − k koji nema realnih nul-toˇcaka unutar intervala [a, b] i stoga je na tom intervalu stalnog znaka. Nadalje, polinom Qn je okomit na polinom Qn (x) = (x − x1 ) · · · (x − xk ) , jer je ovaj stupnja k < n . Stoga, Rn−k (x) Z b 2 Qn Qn (x) ρ(x)dx = 0 (Qn | )= R n −k a Rn−k (x) sˇ to je ponovo kontradikcija, jer je podintegralna funkcija stalnog znaka na [a, b] . 2. Gramm-Schmidtov postupak ortogonalizacije 2.1. Postupak ortogonalizacije. Postupak nalaˇzenja niza ortogonalnih vektora u nekom vektorskom prostoru vaˇzan je postupak. Opisat c´emo op´ceniti algoritam za zamjenjivanje linearno nezavisnih vektora medusobno ortogonalnim vektorima. Taj se postupak naziva Gramm-Schmidtov postupak ortogonalizacije. Pretpostavimo da su x1 , . . . , xn bilo koji linearno nezavisni vektori. Definirajmo e1 := x1 . Vektor e2 c´emo izabrati tako da on zadovoljava sljede´ca svojstva: • okomit je na e1 , • potprostor L(e1 , e2 ) podudara se s potprostorom L(x1 , x2 ) , 7. FOURIEROV RED PO ORTOGONALNIM SUSTAVIMA 13 • sustav (O; x1 , x2 ) ima istu orijentaciju kao i sustav (O; e1 , e2 ) . Sva ta svojstva zadovoljava sljede´ci vektor: e2 := x2 − (x2 | e1 ) e1 . ke1 k2 Tre´ci vektor dobivamo formulom e3 := x3 − (x3 | e1 ) (x3 | e2 ) e1 − e2 , 2 ke1 k ke2 k2 i op´cenito, ek := xk − (xk | ek−1 ) (xk | e1 ) e1 − . . . − ek−1 . ke1 k2 kek−1 k2 Ukoliko nakon svakog koraka vektor ek zamijenimo s ek /kek k , dobit c´emo ortonormirani niz istih svojstava. x3 e3 e2 e1 Da je ovaj sistem zaista ortogonalan, lako provjeravamo indukcijom. Vrijedi (e2 | e1 ) = (x2 | e1 ) − (x2 | e1 ) (e1 | e1 ) = 0. ke1 k2 ˇ Neka je k < n . Pretpostavimo da je (ej | ek ) = 0 za sve j < n , j = 6 k . Zelimo pokazati da vrijedi i (en | ek ) = 0 . Mnoˇze´ci relaciju (??) skalarno s ek dobivamo (en | ek ) = (xn | ek ) − n −1 X (xn | ej ) (ej | ek ) kej k2 j=0 (xn | ek ) (ek | ek ) = 0. = (xn | ek ) − kek k2 2.2. Primjer. Ortogonalizirajmo skup {x1 , x2 , x3 } ⊂ R3 ako je x1 = (1, 1, 1) , x2 = (1, 1, 0) , x3 = (1, 0, 0) . 14 7. FOURIEROV RED PO ORTOGONALNIM SUSTAVIMA Prema navedenom algoritmu imamo e1 = x1 = (1, 1, 1), (x2 | e1 ) e1 e2 = x2 − ke1 k2 = (1, 1, 1) − 32 (1, 1, 1) = ( 13 , 31 , − 23 ) (x3 | e1 ) (x3 | e2 ) e3 = x3 − e1 − e2 , 2 ke1 k ke2 k2 = (1, 0, 0) − 31 (1, 1, 1) − 1 3 2 3 ( 31 , 31 , − 23 ) = ( 21 , − 12 , 0). Dobiveni vektori su okomiti. Pritom vrijedi L(e1 ) = L(x1 ) , L(e1 , e2 ) = L(x1 , x2 ) , - svi sustavi imaju istu orijentaciju. L(e1 , e2 , e3 ) = L(x1 , x2 , x3 ) . Takoder, Rezultat ortogonalizacije bitno ovisi o poretku vektora koje ortogonaliziramo. Razlog tome je upravo zahtjev da u svakom koraku k vrijedi L(e1 , . . . , ek ) = L(x1 , . . . , xk ) . Tako na primjer, ako u gornjem primjeru izaberemo poredak x2 , x3 , x1 dobit c´emo za rezultat niz (1, 1, 0) , ( 12 , − 12 , 0) , (0, 0, 1) . Uzmemo li pak poredak x3 , x2 , x1 , ortogonalizacijom dobivamo kanonsku bazu u R3 . Provjerite! ∗*∗ Ovaj postupak moˇzemo zamisliti i u poneˇsto izmjenjenim okolnostima. Zamislimo da imamo nekoliko ortogonalnih vektora {e1 , . . . , ek } koji ne cˇine bazu. Lako je zamisliti postupak koji c´e te vektore nadopuniti do bilo koje baze prostora: dovoljno je dodavati jedan po jedan vektor {xk+1 , . . . , xn } koji se ne nalazi u potprostoru razapetom s ve´c odabranim vektorima. Sad Gramm-Schmidtovim postupkom ortogonaliziramo zadani skup. Pri tom prvih k vektora ostaje nepromijenjeno (jer su ve´c ortogonalni). Istaknimo ovaj zakljuˇcak. Pri tom rijeˇc ortogonalni moˇzemo po potrebi zamijeniti s ortonormirani. Teorem 2.1. Svaki se ortogonalni (ortonormirani) skup e1 , . . . , ek moˇze nadopuniti do ortogonalne (ortonormirane) baze. 2.3. Legendreovi polinomi. Neka je P vektorski prostor svih polinoma, promatranih na intervalu [−1, 1] . Istaknimo u njemu sljede´ce polinome f 0 (x) = 1 , f 1 (x) = x , f 2 (x) = x2 .. . Tada se svaki polinom f (x) = an xn + . . . + a1 x + a0 moˇze napisati u obliku f = an f n + . . . + a1 f 1 + a0 f 0 . Kaˇzemo da smo polinom f rastavili po bazi (f n ) 7. FOURIEROV 15 RED PO ORTOGONALNIM SUSTAVIMA u prostoru P . Medutim, taj sistem nije ortogonalan na intervalu [−1, 1] . Vrijedi npr. Z 1 2 (f 1 | f 3 ) = x · x3 dx = . 5 −1 Sistem ortogonalnih polinoma moˇzemo dobiti od sistema 1 , x2 ,. . . s pomo´cu Gramm–Schmidtovog postupka ortogonalizacije. Stavimo P0 (x) = 1 P1 (x) = x . R1 Vrijedi (P0 | P1 ) = −1 x dx = 0 . Polinom P2 potraˇzimo u obliku: P2 (x) = ax2 + bx + c . Tada mora biti (P2 | P0 ) = Z 1 −1 (ax2 + bx + c)dx = 32 a + 2c = 0 =⇒ a = −3c (P2 | P1 ) = Z 1 −1 (ax3 + bx2 + cx)dx = 32 b = 0 . Dobili smo P2 (x) = c(−3x2 + 1) . Konstantu c c´emo izabrati tako da bude P2 (1) = 1 . Dobivamo 3 1 P2 (x) = x2 − . 2 2 Nastavimo ovakav postupak s polinomom P3 (x) = ax3 + bx2 + cx + d . Dobiti c´emo: P3 (x) = 5 3 3 x − x, 2 2 itd. Na ovaj naˇcin dobivamo sistem (Pn ) polinoma koji su ortogonalni na intervalu [−1, 1] . Polinom Pn naziva se n–ti Legendreov polinom. Zadatak 2.1. Gramm-Schmidtovim postupkom ortogonaliziraj funkcije f 0 (x) = 1 , f 2 (x) = x2 , f 1 (x) = x , f 3 (x) = x3 na intervalu [0, 1] , uz teˇzinsku funkciju ρ(x) = x . ⊲ Raˇcunat c´emo ortogonalan niz (wn) po formuli . Izraˇcunajmo potrebne skalarne produkte i norme. w0 (x) = f 0 (x) = 1, Z 1 x dx = 12 , kw0 k2 = 0 (f 1 | w0 ) = Z 0 1 x · 1 · x dx = 13 . 16 7. FOURIEROV RED PO ORTOGONALNIM SUSTAVIMA Zato w1 (x) = f 1 (x) − (f 1 | w0 ) w0 (x) = x − kw0 k2 1 3 1 2 · 1 = x − 23 , Da bi odredili w2 , raˇcunamo 2 kw1 k = (f 2 | w0 ) = (f 2 | w1 ) = Z 1 (x − 32 )2 x dx = 0 Z 1 x2 · 1 x dx = 0 Z 1 36 , 1 36 , 1 0 x2 (x − 23 ) x dx = 1 30 . i odavde w2 = f 2 (x) − (f 2 | w1 ) (f 2 | w0 ) w0 (x) − w1 (x) kw0 k2 kw1 k2 1 4 1 2 = x2 − ·1− 1 30 1 36 (x − 23 ) = x2 − 65 x + 3 10 , Joˇs jedan krug raˇcunanja 2 kw2 k = (f 3 | w0 ) = (f 3 | w1 ) = (f 3 | w2 ) = Z 1 (x2 − 56 x + 0 Z 0 1 600 , x3 · 1 · x dx = 15 , 1 x3 (x − 23 ) x dx = 0 Z = 1 0 Z 3 2 10 ) x dx 1 30 , 1 x3 (x2 − 65 x − 3 10 ) x dx = 1 350 , i ponovo po dobivamo w3 (x) = x3 − = x3 − 1 5 1 2 12 7 ·1− 1 30 1 36 x2 + 65 x − (x − 23 ) − 4 35 . ⊳ 1 350 1 600 x2 − 65 x + 3 10  2.4. Rademacherov sistem. Uz klasiˇcne ortogonalne sisteme, sastavljene od klasa lijepih’, beskonaˇcno-diferencijabil funkcija poput sinusoida ili polinoma, vaˇzni su i ortogonalni sistemi sasvim drukˇcijih svojstava, koje saˇcinjavaju diskontinuirane funkcije. Jedan od tih sistema je i Rademacherov sistem. 7. FOURIEROV 17 RED PO ORTOGONALNIM SUSTAVIMA Definirajmo niz funkcija {rn } na intervalu [0, 1] na naˇcin r0 (x) = 1, ( r1 (x) = r2 (x) = ( 0 n . Neka je In bilo koji od 2n intervala na kojima je rn konstantna (i iznosi +1 ili −1 ). Tada je Z Z rm (x)dx. rn (x)rm (x)dx = ±1 · In In Kako je m > n , to funkcija rm poprima na intervalu In vrijednosti +1 ili −1 , na svakom od 2m−n dijelova na koje se raspada taj podinterval. Funkcija rm c´e imati vrijednost +1 na polovini tih podintervala, te −1 na preostaloj polovini. Zato je Z rm (x)dx = 0 In te je Z 0 n 1 rn (x)rm (x)dx = 2 Z X k=1 In rn (x)rm (x) = 0 18 7. FOURIEROV RED PO ORTOGONALNIM SUSTAVIMA S druge strane je 1 Z Z rn (x)2 dx = 1 1 · dx = 1, 0 0 te je {rn } ortonormirani sistem. ⊳ Zadatak 2.3. Funkciju f (x) = x razvij u red po Rademacherovom sistemu. ⊲ Kako je {rn } ortonormirani sistem, to c´e Fourierov red glasiti f = ∞ X (f | rn )rn . n=0 Nulti koeficijent dobivamo iz (f | r0 ) = Z 1 0 x · 1 · dx = 21 . Pri raˇcunanju ostalih koeficijenata koristimo prikaz rn (x) = (−1)j−1 , x∈ Zato je (f | rn ) = = Z n 1 f (x)rn (x)dx = 0 2n X (−1) j−1 1 22n+1  1  2 j=1 = j − 1 j  , j = 1, 2, . . . , 2n . , 2n 2n n 2 Z X n j2 j/2n x · (−1)j−1 dx j=1 (j−1)/2n 2  j − 1 2  − 2n 2 X (−1)j−1 (2j − 1) j=1 1 [1 − 3 + 5 − 7 + . . . − (2n − 1)] 22n+1 1 1 = 2n+1 [(−2) · 2n−1 ] = − n+1 2 2 = Dobili smo: ∞ f (x) = X 1 1 − rn (x). ⊳ r0 2n+1 n=1 7. FOURIEROV 19 RED PO ORTOGONALNIM SUSTAVIMA ˇ sevljevi polinomi 3. Cebiˇ ˇ sevljevih polinoma definiran je formulom Niz (Tn ) Cebiˇ Tn (x) = cos n(arc cos x). Ovom su formulom uistinu definirani polinomi. Stavimo t := arc cos x . Tada je x = cos t i Tn (x) = cos nt . Znamo da se cos nt uvijek moˇze prikazati u obliku polinoma po cos t . Npr. T0 (x) = cos 0 · t = 1 T1 (x) = cos t = x T2 (x) = cos 2t = cos2 t − sin2 t = 2 cos2 t − 1 = 2x2 − 1 T3 (x) = cos 3t = cos3 t − 3 cos t sin2 t = 4 cos3 t − 3 cos t = 4x3 − 3x ˇ sevljevi polinomi zadovoljavaju rekurzivnu formulu, koju c´emo itd. Op´cenito, Cebiˇ izvesti iz identiteta cos(n + 1)t + cos(n − 1)t = 2 cos t cos nt . Supstitucija t = arc cos x daje Tn+1 (x) = 2x Tn (x) − Tn−1 (x) , T0 (x) = 1, T1 (x) = x. n > 1, Odavde zakljuˇcujemo da je Tn polinom stupnja n , s vode´cim koeficijentom 2n−1 . ˇ sevljevih polinoma Iz rekurzivne formule lako dobivamo nekoliko prvih Cebiˇ T0 (x) = 1 T1 (x) = x T2 (x) = 2x2 − 1 T3 (x) = 4x3 − 3x T4 (x) = 8x4 − 8x2 + 1 T5 (x) = 16x5 − 20x3 + 5x T6 (x) = 32x6 − 48x4 + 18x2 − 1 .. . 20 7. FOURIEROV T0 RED PO ORTOGONALNIM SUSTAVIMA 1 T1 1 T2 T4 T3 ˇ sevljevih polinoma iskazano je u sljede´cem teoremu Temeljno svojstvo Cebiˇ Teorem 3.1. Polinom 1 2n−1 - svim polinomima stupnja Tn (x) ima svojstvo da medu n s vode´cim koeficijentom 1 najmanje odstupa od nule na intervalu [−1, 1] . 1 Maksimalno odstupanje tog polinoma je ± n−1 i ono se dostiˇze u n + 1 toˇcki. 2 Dokaz. Neka je Pn (x) = xn + . . . bilo koji polinom s vode´cim koeficijentom jednakim 1 i neka vrijedi |Pn (x) 6 1 2n−1 , za |x| 6 1. 1 Trebamo dokazati da je Pn (x) = 2n−1 Tn (x) . ˇ sevljevih polinoma: ako je Pri tom c´emo koristiti samo jedno svojstvo Cebiˇ kπ xk = cos( n ) , onda vrijedi Tn (xk ) = cos kπ = (−1)k , k = 0, 1, . . . n. Promotrimo polinom Q(x) = 1 Tn (x) − Pn (x). 2n−1 7. FOURIEROV 21 RED PO ORTOGONALNIM SUSTAVIMA 1 Budu´ci je vode´ci koeficijent polinoma 2n−1 Tn (x) i P(x) jednak 1, stupanj polinoma Q(x) ne prelazi (n − 1) . Nadalje, vrijedi 1 2n−1 T( xk ) = (−1)k , |Pn (xk )| 6 1. Zato se predznak od Q(xk ) podudara s predznakom od Tn (xk ) , ili je Q(xk ) = 0 . Dakle, u krajnjim toˇckama svakog intervala [xk+1 , xk ] , polinom Q(x) poprima vrijednosti suprotnih predznaka, pa zato Q(x) ima nultoˇcku u svakom od tih intervala. Ako je Q(xk ) = 0 , trebamo uˇciniti malo detaljniju analizu. U tom sluˇcaju, ili je nultoˇcka dvostruka, pa taj polinom opet ima dvije nultoˇcke u intervalu [xk+1 , xk−1 ] , ili u jednom od susjednih intervala polinom ima joˇs jednu nultoˇcku (vidi sliku). Razlog za to je sˇ to u toˇckama xk+1 i xk−1 polinom Q(x) ima vrijednosti istog predznaka 1 . Budu´ci promatranih intervala ima n , zakljuˇcujemo da polinom Q(x) ima barem n nultoˇcaka na intervalu [−1, 1] . Njegoc je stupanj (n − 1) , zato je on identiˇcki jednak 1 Tn (x) . nuli. Dakle, vrijedi P(x) = 2n−1 ˇ sevljevih polinoma koristi u postupku telesPokaˇzimo kako se ovo svojstvo Cebiˇ kopiranja reda potencija, bitnom postupku za dobivanje dobrih polinomijalnih aproksimacija. Zadatak 3.4. Odredi polinom stupnja 4 koji najmanje odstupa od polinoma P(x) = x6 , na intervalu [−1, 1] . - svim polinomima stupnja 6 , s vode´cim koeficijentom 1 , najmanje od ⊲ Medu 1 nule odstupa polinom 32 T6 (x) . Vrijedi 1 32 T6 (x) = = 6 4 2 1 32 (32x − 48x + 18x 9 2 1 x6 − 32 x4 + 16 x − 32 − 1) 9 2 1 x + 32 najmanje odstupa od polinoma P(x) = x6 . Zbog toga polinom Q(x) = 32 x4 − 16 Polinom Q moˇzemo dobiti na joˇs jedan naˇcin, koji bolje opisuje postupak aproksimiranja funkcije s njenim Fourierovim redom. Prikaˇzimo u tu svrhu x6 pomo´cu ˇ sevljevih polinoma (vidi Tablicu 2.1). Cebiˇ x6 = 1 32 (10 + 15 T2 + 6 T4 + T6 ) 1 Ako je vrijednost u nekoj of tih toˇcaka opet jednaka nuli, razmatranje bi se trebalo nastaviti na analogan naˇcin, s istim konaˇcnim zakljuˇckom 22 7. FOURIEROV RED PO ORTOGONALNIM SUSTAVIMA a zatim zanemarimo cˇlan sa T6 : x6 ≈ 1 32 (10 + 15 T2 + 6 T4 2 1 32 10 + 15(2x − 1) 3 4 9 2 1 2 x − 16 x + 32 ⊳ = = ) + 6(8x4 − 8x2 + 1)  Zadatak 3.5. Odredi polinom cˇ etvrtog stupnja koji dobro aproksimira funkciju f (x) = cos x na intervalu [−1, 1] . ⊲ Funkciju f moˇzemo aproksimirati Taylorovim polinomom stupnja 4 : cos ≈ 1 − x4 x2 + 2 24 Pri tom je uˇcinjena pogreˇska koja iznosi R4 = cos(5) (ϑ ) x5 , 5! ϑ ∈ (−1, 1) ˇ i manja je od 1/5! = 0.0083 . Zelimo li posti´ci ve´cu toˇcnost, aproksimirati c´emo funkciju f Taylorovim polinomom ve´ceg stupnja: x2 x4 x6 + + 2 24 720 6 ˇ a zatim polinom x aproksimiramo pomo´cu Cebiˇsevljevih polinoma polinomom cˇevtrtog stupnja, kao u proˇslom zadatku: cos x ≈ 1 − x6 ≈ 1 32 (1 − 18x2 + 48x4 ) . Tako dobivamo x4 1 1 x2 + − (1 − 18x2 + 48x4 ) 2 24 720 32 = 0.99996 − 0.49922 x2 + 0.03958 x4 cos x ≈ − Uˇcinjena pogreˇska je manja od R6 + 1 1 1 1 1 |T6 | < + = 0.00024 . 720 32 7! 720 32 Joˇs ve´cu toˇcnost postiˇzemo ako dodamo i cˇlan 1 x8 = 8! 40 320 1 ≈ 40 320 1 ≈ 40 320 1 (35T)0 + 56T2 + 28T4 + T8 ) 128 1 (35T0 + 56T2 + 28T4 ) 128 1 (7 − 112x2 + 224x4 ) 128 7. FOURIEROV 23 RED PO ORTOGONALNIM SUSTAVIMA Tada c´emo imati cos x ≈ 0.9999580 − 0.4992405 x2 + 0.0396267 x4 Uˇcinjena pogreˇska je sada sigurno manja od 1 1 |T6 | + (8|T6 | + |T8 |) + R8 < 0.000048 ⊳ 720 · 32 40320 · 128 ˇ sevljevih polinoma pomo´cu polinoma x 7→ xn i obratne veze. Tablica 2.1. Prikaz Cebiˇ T0 T1 T2 T3 T4 T5 T6 T7 T8 T9 T10 1 1 x x2 x3 x4 x5 x6 x7 x8 x9 x10 1 −1 1 2 −3 5 −1 1 4 −8 8 −20 18 −7 9 −1 56 −32 50 1 x x2 x3 x4 x5 x6 x7 x8 x9 x10 160 −120 T0 1 1 2 1 4 1 8 1 16 1 32 1 64 1 128 1 256 1 512 16 −48 T1 32 −112 432 −400 64 −256 1120 128 −576 256 −1280 512 T2 T3 T4 T5 T6 T7 T8 T9 T10 1 1 1 3 3 1 4 10 10 1 5 15 35 35 6 21 56 126 126 1 28 84 210 1 7 1 8 36 120 1 9 45 1 10 1 ˇ sevljevih polinoma. 3.1. Svojstva Cebiˇ ˇ sevljevi polnomi zadovoljavaju identitete Teorem 3.2. Cebiˇ (1 − x2 )Tn′ (x) = n(Tn−1 (x) − xTn (x)), (1 − x )(Tn′ (x))2 = n2 (1 − Tn (x))2 , (1 − x2 )Tn′′ (x) − xTn′ (x) + n2 Tn (x) = 0. 2 (7.10) (7.11) (7.12) 24 7. FOURIEROV RED PO ORTOGONALNIM SUSTAVIMA ⊲ Ove se relacije lako dokazuju koriste´ci se vezama x = cos t , Tn (x) = cos nt . Odavde imamo d cos nt n sin nt dt = , = dcost sin t dt   −1 n cos t sin nt − n2 cos nt sin t d n sin nt · = Tn′′ (x) = dt sin t sin t sin3 t Tn′ (x) S ovim supstitucijama, formula (7.10) postaje (1 − cos2 t) · n sin nt = n(cos(n − 1)t − cos t cos nt) sin t odnosno sin t sin nt + cos t cos nt = cos(n − 1)t sˇ to je istina, po adicijskom teoremu. Na sliˇcan naˇcin se dokazuju preostala dva identiteta. ⊳ Fourierov red funkcije f po ortogonalnom sistemu (Tn ) glasi f (x) = c0 T0 (x) + . . . + cn Tn (x) + . . . , Z 1 1 f (x) √ c0 = dx π −1 1 − x2 Z 2 1 f (x) Tn (x) √ cn = dx , n > 1 . π −1 1 − x2 Ovi se integrali raˇcunaju supstitucijom x = cos t , jer je tada Tn (x) = cos nt . Tako gornje formule prelaze u Z 1 π f (cos t)dt, c0 = π 0 Z 2 π cn = f (cos t) cos ntdt π 0 Zadatci za vjeˇzbu 7.1. Vrijedi 1 z + 2 = z 2  2 1 −2 z+ z 1 i z4 + z14 u obliku polinoma po z + . z x . 7.2. Neka je Pn (x) = 2Tn 2 Prikaˇzi z3 + 1 z3 7. FOURIEROV RED PO ORTOGONALNIM SUSTAVIMA 25 1) Odredi rekurzivnu relaciju za ove polinome i pokaˇzi da su im koeficijenti cjelobrojni. 2) Dokaˇzi da vrijedi 1 1 Pn (z + ) = zn + n . z z 7.3. Neka je Un funkcija definirana formulom sin nt = Un (cos t) sin t. ˇ sevljev Dokaˇzi da je Un polinom s cjelobrojnim koeficijentima. Un se naziva Cebiˇ polinom druge vrste. 7.4. Dokaˇzi da polinomi (Un ) zadovoljavaju identiˇcnu rekurzivnu relaciju Un (x) = 2xUn−1 (x) − Un−2 (x), uz poˇcetne uvjete U0 (x) = 1 , U1 (x) = 2x . Dokaˇzi sljede´ce relacije A. Un (x) = 2Tn (x) + Un−2 (x) ; B. Tn (x) = Un (x) − xUn−1 (x) ; Pn C. Un (x) = k=0 xk Tn−k (x) ; 7.6. Izvedi eksplicitne prikaze: ⌊n/2⌋   X n−k (2x)n−2k ; A. Un (x) = (−1)k k 7.5. k=0 ⌊n/2⌋ B. Un (x) = X k=0 (−1)k   n + 1 n−2k 2 x (x − 1)k . 2k + 1