analiza

O spektru nelinearnih operatora

Sanela Halilović
Odsjek Matematika, Prirodno-matematički fakultet
Univerzitet u Tuzli
[email protected]
Samra Pirić
Odsjek Matematika, Prirodno-matematički fakultet
Univerzitet u Tuzli
[email protected]






Sažetak
U ovom radu pokazujemo da preslikavajući spektar nelinearnih operatora ne zadržava neke bitne osobine koje ima spektar linearnih operatora. To ilustriramo nizom primjera.

Ključne riječi: spektar, nelinearni operator, bijekcija

1Uvod

Dobro je poznata važnost spektralne teorije za linearne operatore. Podsjetimo se nekih najznačajnijih osobina spektra linearnih operatora. Neka je X Banachov prostor nad poljem \mathbb{K} realnih ili kompleksnih brojeva.

Definicija 1. Spektar ograničenog linearnog operatora L:X\rightarrow X je skup
(1)
\sigma(L)=\lbrace \lambda\in \mathbb{K}:\lambda I- L\text{ nije bijekcija}\rbrace .

Za svako \lambda \in \mathbb{K}\setminus \sigma (L) postoji rezolventni operator (\lambda I- L)^{-1} koji je ograničen.

Definicija 2. Neka je X kompleksan Banachov prostor, \sigma(L) spektar ograničenog linearnog operatora L i \lambda\in \sigma(L).
1) Kažemo da \lambda pripada točkovnom spektru ako \lambda I-L nije injekcija. Skup svih takvih \lambda označavamo sa \sigma_{p}(L) i nazivamo točkovnim spektrom operatora L.
2) \lambda je element neprekidnog spektra ako je \lambda I-L injekcija i (\lambda I-L) (X) gust potprostor od X. Skup svih takvih \lambda označavamo sa \sigma_{c}(L) i nazivamo neprekidnim spektrom operatora L.
3) \lambda je element rezidualnog spektra ako je \lambda I-L injekcija, ali (\lambda I-L) (X) nije gust potprostor od X. Skup svih takvih \lambda označavamo sa \sigma_{r}(L) i nazivamo rezidualnim spektrom operatora L.

U slučaju da dim X\lt \infty onda je \sigma(L)=\sigma_{p}(L). Primijetimo da \lambda_{0}\in \sigma_{c}(L) povlači da operator (\lambda_{0} I-L) ima inverzni operator (\lambda_{0} I- L)^{-1} ali da taj nije ograničen. Situacija \lambda_{0}\in \sigma_{r}(L) znači da rezolventni operator postoji, ali njegovo područje definicije nije gusto u X; u tom slučaju rezolventni operator može biti ograničen ili neograničen.

Slobodno govoreći, elementi \lambda u subspektru \sigma_{p}(L) karakteriziraju neki gubitak injektivnosti, oni iz \sigma_{r}(L) neki gubitak surjektivnosti, a oni iz \sigma_{c}(L) neki gubitak stabilnosti operatora \lambda I-L.

Ovi dijelovi spektra formiraju disjunktnu podjelu spektra
\sigma (L)=\sigma_{p}(L)\cup \sigma_{c}(L)\cup \sigma_{r}(L).
Vrijedi i sljedeći teorem koji daje tzv. formulu spektralnog preslikavanja polinoma.

Teorem 3. Neka je L linearan operator u Banachovu prostoru X nad poljem \mathbb{K }. Za svaki polinom p:\mathbb{K}\rightarrow \mathbb{K}, p(\lambda)=a_{n}\lambda^{n}+a_{n-1} \lambda^{n-1}+\cdots +a_{1}\lambda+a_{0} vrijedi
(2)
\sigma(p(L))=p (\sigma(L)),
gdje je p(L)=a_{n}L^{n}+\cdots +a_{1}L+a_{0} I i p(\sigma(L))=\lbrace p(\lambda):\lambda\in\sigma(L)\rbrace.

Spektar linearnog operatora \sigma (L) ima sljedeće važne osobine:
\bullet zatvoren je i ograničen skup (dakle kompaktan)
\bullet neprazan je skup kad je \mathbb{K} polje kompleksnih brojeva
\bullet vrijedi formula spektralnog preslikavanja (2).


2Preslikavajući spektar

Kod definiranja spektra nelinearnih operatora, cilj je, po mogućnosti:
\bullet u slučaju linearnog operatora da se svodi na poznati spektar (1),
\bullet da ima bar neke zajedničke osobine s linearnim spektrom (npr. zatvorenost, kompaktnost),
\bullet da sadrži svojstvene vrijednosti operatora.
S obzirom na definiciju (1) predstavljalo bi izazov definirati spektar neprekidnog linearnog operatora F jednostavno s
(3)
\Sigma(F):=\lbrace \lambda\in \mathbb{K}:\lambda I- F \text{ nije bijekcija}\rbrace .
Ovo ćemo nazivati preslikavajući spektar operatora F. Preciznije mogli bismo proučavati spektar injektivnosti
(4)
\Sigma_{i}(F):=\lbrace \lambda\in \mathbb{K}:\lambda I- F \text{ nije injekcija}\rbrace
i spektar surjektivnosti
(5)
\Sigma_{s}(F)=\lbrace \lambda\in \mathbb{K}:\lambda I- F \text{ nije surjekcija}\rbrace ,
pri čemu je \Sigma(F)=\Sigma_{i}(F)\cup \Sigma_{s}(F). Međutim, u nelinearnom slučaju se pokazuje da ovaj pristup nije od neke koristi. Zapravo, ove jednostavne definicje imaju smisla samo u linearnom slučaju, jer tada imamo veoma rigidnu strukturu linearnosti, a također i tako moćno oruđe kao teorem o zatvorenm grafu koji garantira ograničenost inverza ograničenog operatora.

Pokazat ćemo primjerima da preslikavajući spektar (3) ne mora imati ni jednu od poznatih osobina linearnog spektra.

Primjer 2.1 Neka je operator F:\mathbb{R}\rightarrow \mathbb{R} definiran s
F(x)=\sqrt{|x|}.
Odredimo spektre injektivnosti i surjektivnosti, odnosno preslikavajući spektar operatora F. Označimo G(x)=(\lambda I-F)(x)=\lambda x-\sqrt{|x|}.
a) Ispitajmo kad je ovo preslikavanje injektivno.
Za \lambda =0 je G(x)=-\sqrt{|x|}, a to nije injekcija (jer za x\neq 0 vrijedi G(x)=G(-x)). Prema tome, 0 \in \Sigma_{i}(F).
Slika 1: G(x)=-\sqrt{|x|},(\lambda=0)


Neka je sada \lambda \neq 0 i promatrajmo jednakosti
G(x_{1})= \lambda x_{1}-\sqrt{|x_{1}|}=\lambda x_{2}-\sqrt{|x_{2}|}=G(x_{2})
\lambda (x_{1}-x_{2})=\sqrt{|x_{1}|}-\sqrt{|x_{2}|}.
U slučaju da su x_{1} i x_{2} pozitivni i x_{1}\neq x_{2}, vrijedi
\lambda (\sqrt{x_{1}}-\sqrt{x_{2}})(\sqrt{x_{1}}+\sqrt{x_{2}})=\sqrt{x_{1}}-\sqrt{x_{2}}
\lambda (\sqrt{x_{1}}+\sqrt{x_{2}})=1 \Rightarrow \sqrt{x_{2}}=\frac{1}{\lambda}-\sqrt{x_{1}}.
Iz uvjeta \frac{1}{\lambda}-\sqrt{x_{1}}\gt 0 dobivamo \lambda \gt 0 i x_{1} \in (0,\frac{1}{\lambda^{2}}).
U slučaju da su x_{1} i x_{2} negativni i x_{1}\neq x_{2}, vrijedi
\lambda (-\sqrt{|x_{1}|}+\sqrt{|x_{2}|})(\sqrt{|x_{1}|}+\sqrt{|x_{2}|})=\sqrt{|x_{1}|}-\sqrt{|x_{2}|}
\lambda (\sqrt{|x_{1}|}+\sqrt{|x_{2}|})=-1 \Rightarrow \sqrt{|x_{2}|}=-\frac{1}{\lambda}-\sqrt{|x_{1}|}.
Rješavanjem nejednakosti -\frac{1}{\lambda}-\sqrt{|x_{1}|}\gt 0 dobivamo \lambda \lt 0 i x_{1} \in (-\frac{1}{\lambda^{2}},0). Ovim smo pokazali da za proizvoljno \lambda \gt 0, ako uzmemo x_{1} \in (0,\frac{1}{\lambda^{2}}) i
x_{2}=(\frac{1}{\lambda}-\sqrt{x_{1}})^{2}, onda dobivamo G(x_{1})=G(x_{2}). To znači da za \lambda \in (0,\infty) preslikavanje G nije injektivno, odnosno (0,\infty)\subseteq \Sigma_{i}(F). S druge strane, za proizvoljno \lambda \lt 0, ako uzmemo x_{1} \in (-\frac{1}{\lambda^{2}},0) i x_{2}=-(\frac{1}{\lambda}+\sqrt{|x_{1}|})^{2}, onda opet dobivamo G(x_{1})=G(x_{2}). Dakle, (-\infty,0)\subseteq \Sigma_{i}(F). Sveukupno, našli smo spektar injektivnosti
\Sigma_{i}(F)=(-\infty,0)\cup \lbrace 0\rbrace \cup (0,\infty)=\mathbb{R}
Slika 2: G(x) za \lambda=1 i \lambda=0.8


b) Ispitajmo kad je G surjektivno preslikavanje.
Za \lambda =0 je G(x)=-\sqrt{|x|}, a ovo nije surjekcija jer je G(\mathbb{R})=(-\infty,0]. Prema tome, 0\in \Sigma_{s}(F). Neka je sad \lambda \neq 0 i y \in \mathbb{R} proizvoljno. Ispitajmo rješenja jednadžbe
\lambda x-\sqrt{|x|}=y.
Nalazimo: za \lambda \gt 0 je
x=\begin{cases} \frac{1+2\lambda y + \sqrt{1+4\lambda y}}{2\lambda ^{2}} & y \in [-\frac{1}{4\lambda},\infty) \\ \frac{1+2\lambda y - \sqrt{1+4\lambda y}}{2\lambda ^{2}} & y \in [-\frac{1}{4\lambda},0] \\ \frac{-1+2\lambda y+ \sqrt{1-4\lambda y}}{2\lambda ^{2}} & y \in (-\infty,0] \\ \end{cases},
a za \lambda \lt 0 je
x= \begin{cases} \frac{1+2\lambda y - \sqrt{1+4\lambda y}}{2\lambda ^{2}} & y \in (-\infty,0] \\ \frac{-1+2\lambda y + \sqrt{1-4\lambda y}}{2\lambda ^{2}} & y \in [\frac{1}{4\lambda},0] \\ \frac{-1+2\lambda y - \sqrt{1-4\lambda y}}{2\lambda ^{2}} & y \in [\frac{1}{4\lambda}, \infty) \\ \end{cases}.
Tako da za \lambda \neq 0 i proizvoljno y, postoji x \in \mathbb{R} takav da je G(x)=y; odnosno G je surjekcija. Ostaje samo \Sigma_{s}(F)=\lbrace 0\rbrace . Preslikavajući spektar je
\Sigma (F)=\Sigma_{i}(F)\cup \Sigma_{s}(F)=\mathbb{R},
pa vidimo da nije ograničen skup.
Primjer 2.2 Neka je F:\mathbb{R}\rightarrow \mathbb{R} definiran s
(6)
F(x)= \left\lbrace \begin{array}{ll} x & \text{ako je }|x|\gt 1\text{}, \\ x^{2} & \text{ako je }0\leq x\leq 1\text{}, \\ -x^{2} & \text{ako je }-1\leq x\leq 0.\text{} \end{array} \right.
U slučaju da je \lambda =0, promatrajmo preslikavanje
G(x)=-F(x)=\begin{cases} -x & \text{ ako je } |x|\gt 1,\\ -x^{2} & \text{ ako je } 0\leq x\leq 1,\\ x^{2} & \text{ ako je } -1\leq x\leq 0 \end{cases}
i pokažimo da je u pitanju bijekcija (vidi sliku 3).
Slika 3: G(x) za \lambda=0


Za y\in (-\infty, -1)\cup (1,\infty),\quad \exists x=-y, tako da G(x)=y. Za y\in [-1,0],\quad \exists x=\sqrt{-y}, tako da G(x)=y. I za y\in [0,1], \exists x=-\sqrt{y},\quad G(x)=y. Dakle, G(\mathbb{R})=\mathbb{R}, pa je G surjekcija i 0\notin \Sigma_{s}(F). Jasno je i da je injekcija jer iz svake jednadžbe G(x_{1})=G(x_{2}) slijedi da je x_{1}=x_{2}. Injektivnost se može dokazati i činjenicom da je G neprekidna i stalno opadajuća funkcija od +\infty do -\infty na čitavoj realnoj osi. Dakle, 0\notin \Sigma_{i}(F). Kad je \lambda =1 imamo:
G(x)=\begin{cases} 0 & \text{ ako je } |x|\gt 1,\\ x-x^{2} & \text{ ako je } 0\leq x\leq 1,\\ x+x^{2} & \text{ ako je } -1\leq x\leq 0. \end{cases}
G_{max}=G(\frac{1}{2})=\frac{1}{4}, a G_{min}=G(-\frac{1}{2})=-\frac{1}{4} (vidi sliku 4).
Slika 4: G(x) za \lambda=1


Budući da je G neprekidna funkcija, vrijedi G(\mathbb{R})=[-\frac{1}{4}, \frac{1}{4}], pa nije surjekcija i 1\in \Sigma_{s}(F). Jasno je da G nije ni injekcija jer npr. G(2)=G(3)=0. Prema tome, 1\in \Sigma_{i}(F). Neka je \lambda \in (0,1) (vidi sliku 5). Tada
(7)
G(x)=\begin{cases} x(\lambda-1) & \text{ ako je } |x|\gt 1,\\ \lambda x-x^{2} & \text{ ako je } 0\leq x\leq 1,\\ \lambda x+x^{2} & \text{ ako je } -1\leq x\leq 0. \end{cases} G^{\prime} (x)= \begin{cases} \lambda-1 & \text{ ako je } |x|\gt 1,\\ \lambda -2x & \text{ ako je } 0\leq x\leq 1,\\ \lambda +2x & \text{ ako je } -1\leq x\leq 0. \end{cases}.



x (-\infty,-1) (-1,-\lambda/2) (-\lambda/2,0) (0,\lambda/2) (\lambda/2,1) (1,\infty)
G'(x) - - + + - -
G(x) \searrow \searrow \nearrow \nearrow \searrow \searrow

Budući da je G neprekidna funkcija na \mathbb{R} i dostiže lokalni minimum za x=-\frac{\lambda}{2}, a lokalni maksimum za x=\frac{\lambda}{2}, slijedi da nije injekcija. Tako npr. za x_{1}\in (0,\frac{\lambda}{2}) i x_{2}\in (\frac{\lambda}{2}, 1), iz jednadžbe G(x_{1})=G(x_{2}) dobivamo x_{2}=\lambda -x_{1}. Ovim smo pokazali (0,1)\subseteq \Sigma_{i}(F).
Slika 5: G(x) za \lambda=\frac{1}{2}


Jasno je da je G surjekcija jer je neprekidna funkcija i
\lim_{x \to -\infty} G(x)= + \infty,\quad \lim_{x \to +\infty} G(x)= - \infty.
U slučaju \lambda \in (1,2) je



x (-\infty,-1) (-1,-\lambda/2) (-\lambda/2,0) (0,\lambda/2) (\lambda/2,1) (1,\infty)
G`(x) + - + + - +
G(x) \nearrow \searrow \nearrow \nearrow \searrow \nearrow

tako da opet nemamo injekciju, pa (1,2)\subseteq \Sigma_{i}(F). Sada je
\lim_{x \to -\infty} G(x)= - \infty,\quad \lim_{x \to +\infty} G(x)= +\infty
pa G jest surjekcija. Iz (7) slijedi: za \lambda \geq 2 je G^{\prime}(x)\gt 0, te G stalno raste od -\infty do +\infty; a za \lambda \lt 0 je G^{\prime}(x)\lt 0, pa G stalno opada od +\infty do -\infty na čitavoj realnoj osi. Ovo znači, zbog neprekidnosti, da je za ovakve \lambda preslikavanje G bijekcija. Našli smo, konačno, spektre:
\Sigma_{i}(F)=\Sigma (F)=(0,2), \qquad \Sigma_{s}(F)=\lbrace 1\rbrace .
Dakle, \Sigma (F) nije zatvoren skup.


Jedna od najvažnijih osobina linearnog spektra je ta da je on uvijek neprazan u slučaju kad je \mathbb{K} polje kompleksnih brojeva. Međutim, pokazuje se da ovo više ne vrijedi kad je riječ o nelinearnom operatoru.

Primjer 2.3 Neka je operator F:\mathbb{C}^{2}\rightarrow \mathbb{C}^{2} definiran s
(8)
F(z,w)=(\overline{w},i\overline{z}).
Tada je preslikavanje (\lambda I-F)(z,w)=(\lambda w-i\overline{z}), za svako \lambda \in \mathbb{C}, bijekcija na \mathbb{C}^{2} s inverzom
(\lambda I-F)^{-1} (\zeta,\omega)= \left( \frac{\overline{\lambda}\zeta + \overline{\omega}}{i + |\lambda |^{2}}, -\frac{\overline{\lambda}\omega + i \overline{\zeta}}{i - |\lambda |^{2}} \right).
Slijedi da je:
\Sigma_{i}(F)= \Sigma_{s}(F)=\Sigma(F)=\emptyset.
Prema tome preslikavajući spektar ovog nelinearnog operatora (9) je prazan skup.
Primjer 2.4 Neka je operator F:\mathbb{R}\rightarrow \mathbb{R} definiran s
(9)
F(x)= \begin{cases} 0 & \text{ ako je } x\leq 1, \\ x-1 & \text{ ako je } 1 \lt x \lt 2, \\ 1 & \text{ ako je } x \geq 2. \end{cases}
Nije teško pokazati da je spektar \Sigma (F)=[0,1]. Za polinom p(z)=z^{2} vrijedi p(\Sigma (F))=[0,1]. S druge strane, iz činjenice da je F^{2} (x)\equiv 0 slijedi da je \Sigma (p(F))=\Sigma (F^{2})=\lbrace 0\rbrace. Ovaj primjer pokazuje da ne vrijedi formula spektralnog preslikavanja (2).

Spektar injektivnosti (4) tijesno je povezan s točkovnim spektrom
\sigma_{p} (F):=\lbrace \lambda \in \mathbb{K}: F(x)= \lambda x \text{ za neko } x\neq 0\rbrace .
Kao i u linearnom slučaju, elemente \lambda \in \sigma_{p} (F) nazivat ćemo svojstvenim vrijednostima operatora F. U slučaju da je F(0)=0, vrijedi inkluzija
\sigma_{p} (F)\subseteq \Sigma_{i} (F)
koja može biti i stroga. U Primjeru 2.1 imamo da je \sigma_{p} (F)=\mathbb{R}\setminus \lbrace 0\rbrace \subset \mathbb{R}= \Sigma_{i}(F). Naravno, za linearne operatore L uvijek vrijedi da je \sigma_{p} (L)=\sigma_{i} (L), po definiciji.

Pogledajmo sad spektar surjektivnosti (5). Za z\in X definiramo translaciju F_{z} operatora F s
(10)
F_{z}(x)=F(x)+z
Sljedeći rezultat daje nam vezu između spektra surjektivnosti \Sigma_{s}(F) i točkovnog spektra \sigma_{p}(F_{z}) svih translacija (10).

Propozicija 4. Za neprekidan operator F:X \to X vrijedi jednakost
(11)
\mathbb{K}\setminus \Sigma_{s}(F)=\bigcap_{z\in X \setminus \lbrace -F(0)\rbrace } \sigma_{p}(F_{z}).

Dokaz.
(i) Pokažimo najprije da
\bigcap_{z\in X \setminus \lbrace -F(0)\rbrace } \sigma_{p}(F_{z})\subseteq \mathbb{K}\setminus \Sigma_{s}(F).
Neka za svako z\neq -F(0) vrijedi \lambda \in \sigma_{p}(F_{z}). Tada:
(\exists x_{z}\neq 0) \lambda x_{z}=F_{z} (x_{z})=F(x_{z})+z \Rightarrow (\lambda I-F)(x_{z})=z,
što znači da z \in R(\lambda I-F). Dakle, \forall z \in X \setminus {-F(0)}, z \in R(\lambda I-F), pa R(\lambda I-F)\supseteq X \setminus {-F(0)}. Budući da još, očito, vrijedi i -F(0) \in R(\lambda I-F), imamo: R(\lambda I-F)\supseteq X. Svakako je R(\lambda I-F)\subseteq X, pa R(\lambda I-F)=X. Dakle \lambda I-F je surjektivno preslikavanje, tj. \lambda \in \mathbb{K}\setminus \Sigma_{s}(F).
(ii) Neka je sada \lambda \in \mathbb{K}\setminus \Sigma_{s}(F), tj. \lambda I-F je surjektivno, te vrijedi:
(\forall z \in X) (\exists x_{z} \in X) \lambda x_{z}-F(x_{z})=z. Ako z\neq -F(0), onda x_{z}\neq 0, a to znači da je x=x_{z} netrivijalno rješenje jednadžbe F_{z}(x)=\lambda x. Prema tome, \lambda \in \sigma_{p}(F_{z}). Budući da ovo vrijedi \forall z \in X \setminus {-F(0)}, onda je
\lambda \in \bigcap_{z\in X \setminus \lbrace -F(0)\rbrace } \sigma_{p}(F_{z}).
Dakle,
\mathbb{K}\setminus \Sigma_{s}(F)\subseteq \bigcap_{z\in X \setminus \lbrace -F(0)\rbrace } \sigma_{p}(F_{z}).
Na osnovi pokazanog u (i) i (ii) slijedi tražena jednakost.
\ \blacksquare


Za F(x)=\sqrt{|x|} iz Primjera 2.1, translatirana funkcija je F_{z}(x)=\sqrt{|x|}+z. Budući da je za svako z\in \mathbb{R} točkovni spektar \sigma_{p}(F_{z})=\mathbb{R}\setminus \lbrace 0\rbrace, na osnovi (11) slijedi da je \Sigma_{s}(F)=\lbrace 0\rbrace.

Primjer 2.5 U prostoru neprekidnih funkcija C[0,1] dan je Hammersteinov integralni operator
H(x)(s)=s^{\beta +1} \int_{0}^{1}t^{\beta} \sin x(t)dt \qquad (0\leq s\leq 1, \beta \geq 0).
Operator H je kompozicija H=KF, nelinearnog operatora F definiranog s
F(x)(t)=\sin x(t)
i linearnog Fredholmova integralnog operatora
Ky(s)=\int_{0}^{1} s^{\beta +1}t^{\beta}y(t)dt.
Odredimo točkovni spektar \sigma_{p}(H). Operator H je kompaktan. Za neprekidnu funkciju x_{n}(t)\equiv n\pi\equiv 0, pa H(x_{n})=0=0x_{n}. To znači da 0 \in \sigma_{p}(H). Razmotrimo sad jednadžbu H(x)=\lambda x, za \lambda \neq 0. Imamo
H(x)(s)=s^{\beta +1} \int_{0}^{1}t^{\beta} \sin x(t)dt=\lambda x(s)
x(s)=cs^{\beta +1} \text{ za neko } c\neq 0.
(12)
\lambda =\frac{1}{c}\int_{0}^{1}t^{\beta} \sin (ct^{\beta +1})dt=:\psi (c).
Vrijedi i obrnuto, svaka funkcija x(t)=ct^{\beta +1}, c \neq 0 je svojstvena funkcija operatora H koja odgovara svojstvenoj vrijednosti \lambda =\psi(c). Prema tome,
\sigma_{p}(H)=\lbrace \psi(c):c \in \mathbb{R}\setminus \lbrace 0\rbrace \rbrace .
Uvedimo u (12) zamjenu varijabli ct^{\beta +1}=\taudt=\frac{1}{\beta +1} c^{-\frac{1}{\beta +1}} \tau ^{-\frac{\beta}{\beta +1}}d\tau, pa dobivamo
\lambda =\psi (c)=\frac{1}{(\beta +1) c^{2}}\int_{0}^{c} \sin \tau d\tau =\frac{1-\cos c}{(\beta +1) c^{2}} \quad (c\neq 0).
Sada je
\lim_{c \rightarrow 0} \psi (c)=\lim_{c \rightarrow 0}\frac{1-\cos c}{(\beta +1) c^{2}}= \lim_{c \rightarrow 0}\quad \frac{\sin c}{2(\beta +1) c}=\frac{1}{2(\beta +1)}.
Možemo dodefinirati funkciju \psi u točki 0 tako da bude neprekidna na \mathbb{R}:
\widetilde{\psi}(c)= \begin{cases} \frac{1-\cos c}{(\beta +1) c^{2}} & c\neq 0, \\ \frac{1}{2(\beta +1)} & c =0. \end{cases}
Budući da je funkcija \widetilde{\psi} neprekidna, vrijedi
0\leq \widetilde{\psi}(c) \leq \frac{1}{2(\beta +1)}, \quad (-\infty \lt c\lt \infty),
pri čemu se dostižu sve vrijednosti između lijeve i desne granice. Lijeva strana nejednakosti dostiže se u točkama c=2k \pi, k \in \mathbb{Z}, tj. \psi(2k\pi)= 0. Desna strana nejednakosti ne dobiva se ni za jedno c\neq 0, jer iz \psi (c)=\frac{1}{2(\beta +1)} slijedi
\frac{1}{2(\beta +1)}=\frac{1-\cos c}{(\beta +1) c^{2}}=\frac{2\sin^{2}(\frac{c}{2})}{(\beta +1) c^{2}}= \frac{1}{2(\beta +1)}\Big (\frac{\sin \frac{c}{2}}{\frac{c}{2}}\Big )^{2},
\sin \frac{c}{2}=\frac{c}{2}\Rightarrow c = 0.
Prema tome \sigma_{p} (H)=\Big [0,\frac{1}{2(\beta +1)}\Big ).

3Zaključak

Ovim primjerima pokazali smo da nam treba drugačiji pristup pri definiranju spektra nelinearnih operatora. Za nelinearni neprekidni operator F i neku klasu neprekidnih operatora \mathcal{M}(X) koja sadržava F možemo definirati rezolventni skup
\rho (F)=\Big \lbrace \lambda \in \mathbb{K}: \lambda I-F \text{ je bijekcija i } (\lambda I-F)^{-1} \in \mathcal{M}(X)\Big \rbrace
i spektar
\sigma (F)=\mathbb{K} \setminus \rho (F).
Ovisno o tome što uzmemo za klasu \mathcal{M}(X) (npr. neprekidno diferencijabilni, Lipschitz neprekidni, stabilno rješivi ili epi operatori) dobivamo razne spektre. Nazive su dobivali po matematičarima koji su ih prvi uveli: Rhodius, Neuberger, Kachurovski, Feng, itd. Na ovaj način dobivaju se spektri koji imaju samo neke dobre osobine koje imaju spektri linearnog operatora. Nelinearne spektralne teorije i dalje su u razvoju, a opseg njihove primjene vrlo je širok.
Bibliografija
[1] S. Kurepa, Funkcionalna analiza Elementi teorije operatora, Školska knjiga, Zagreb, 1990.
[2] J. Appell, E. Pascale and A. Vignoli, Nonlinear Spectral Theory, Walter de Gruyter Berlin, New York, 2004
[3] S. Kurepa, Matematička analiza 1, Školska knjiga, Zagreb, 1997.

Fourierov red i Fourierova transformacija

 

Domagoj Matijević
Odjel za matematiku
Sveučilišta u Osijeku
Stjepan Poljak
Odjel za matematiku
Sveučilišta u Osijeku

 
Sažetak
U ovom radu objasniti ćemo kako svaku periodičku funkciju možemo napisati kao sumu (ne nužno konačnu) sinusa različitih amplituda, faza i frekvencija – Fourierov red. U drugom dijelu rada motivirati ćemo formulu za Fourierovu transformaciju preko Fourierova reda te objasniti formulu za diskretnu Fourierovu transformaciju.

1Uvod

U današnjoj se literaturi Fourierova transformacija najčešće prezentira čitatelju kao cjelovit samostalan izraz te se rijetko motivira čitatelja da shvati kako je sam izraz nastao (vidi [1], [3], [2]). U ovom ćemo članku pokušati objasniti pojam Fourierova reda te preko njega motivirati nastanak formule za Fourierovu transformaciju, slijedeći na taj način i sam povijesni tijek njihovih nastanaka.

Fourierova analiza proizlazi iz glavne ideje da svaku periodičku funkciju možemo zapisati kao sumu (ne nužno konačnu) sinusa različitih amplituda, faza i frekvencija. Takva suma naziva se Fourierov red (vidi npr. [6]).


Slika 1: Joseph Fourier (1768.-1830.)


Sama motivacija koju je Jean Baptiste Joseph Fourier (Slika 1) postavio je problem rješavanja parcijalne diferencijalne jednadžbe za širenja topline. Fourier je teorije o širenju topline i Fourierovu redu razvio u razdoblju od 1804. do 1807. godine te poslije objavio u svom bitnom radu O širenju topline u čvrstim tijelima objavljenom 1822. godine. Taj rad naišao je na određen otpor, ponajviše matematičara (među njima su bili i Lagrange i Laplace); razlog tomu je Fourierova tvrdnja da se sve (pa i nederivabilne) funkcije mogu razviti u trigonometrijski red. Iako derivabilnost neće biti nužan uvjet za zapisivanje funkcije kao sume sinusa, poslije će se ipak pokazati da je Fourierova tvrdnja ipak bila preambiciozna. Danas su Fourierove teorije dalje razvijene; među njima je najvažnija diskretna Fourierova transformacija (DFT) i algoritam za brzo i efikasno izračunavanje DFT-a u O(n\log n) vremenu - brza Fourierova transformacija (FFT), koja je temelj velikog dijela modernih multimedijskih primjena (MP3, JPEG).

Iako je Fourierova transformacija proizašla kao rezultat proučavanja Fourierova reda, u današnjoj literaturi često se servira kao „gotov alat” i vrlo je teško shvatiti iz same formule Fourierove transformacije motivaciju za njeno nastajanje.

Struktura ovog rada sastoji se od dvaju osnovnih poglavlja. U 2. poglavlju formalno ćemo definirati Fourierov red, dok ćemo u 3. poglavlju motivirati Fourierovu transformaciju preko Fourierova reda te objasniti formulu za diskretnu Fourierovu transformaciju.

2Fourierov red

Pri razvoju funkcije u Fourierov red, uvjet je da je funkcija periodična. Prvo pitanje koje se postavlja - možemo li to napraviti i za neperiodične funkcije na nekom intervalu? Odgovor je: da, ako ih učinimo periodičnima! Tada taj interval postaje period te funkcije i ponavlja se beskonačno mnogo puta. Iako sva razmatranja koja ćemo ovdje iznijeti vrijede za periodične funkcije nekog općenitog perioda [-L/2, L/2 ], u ostatku teksta ćemo zbog jednostavnosti pretpostaviti da je promatrana funkcija periodična na nekom intervalu [-\pi, \pi ]. Dakle, za danu funkciju (ili signal, ako je funkcija ovisna o vremenu) f:\mathbb{R}\rightarrow\mathbb{R}, periodičnu na [-\pi,\pi], želimo pronaći brojeve a_{0}, a_{k}, \phi_{k}\in\mathbb{R}, za k=1,\ldots,n tako da vrijedi
f(x)=\frac{a_{0}}{2}+\sum_{k=1}^{n}{a_{k}\sin{(k x+\phi_{k})}}, \quad n\ge 0.


U ovom izrazu a_{k} ima ulogu amplitude, a \phi_{k} ulogu faze za sinus funkciju frekvencije k. Član \frac{a_{0}}{2} je poseban i on služi za translaciju funkcije duž y-osi (često se naziva „DC komponenta”). Ako u gornjem izrazu primijenimo adicijsku formulu za sinus funkciju,
\sin{(\alpha+\beta)}=\sin{\alpha}\cos{\beta}+\cos{\alpha}\sin{\beta},


vidimo da možemo zapisati Fourierov red i drugačije:
f(x)=\frac{a_{0}}{2}+\sum_{k=1}^{n}{(a_{k}\sin{(k x)}\cos{\phi_{k}}+a_{k}\cos{(k x)}\sin{\phi_{k}})}.


Primijetimo da u ovoj formuli \sin{\phi_{k}} i \cos{\phi_{k}} ne ovise ni o varijabli x ni o k - to su brojevi koji su dio koeficijenata uz \sin{(k x)} i uz \cos{(k x)}. Stoga, možemo definirati nove koeficijente a_{k}, b_{k}\in\mathbb{R} u kojima će implicitno biti uključeni i \phi_{k}. Tako da nova formula izgleda ovako:
(1)
f(x)=\frac{a_{0}}{2}+\sum_{k=1}^{n}{(a_{k}\sin{(k x)}+b_{k}\cos{(k x)})}, \quad n\ge 0.

Traženje Fourierova reda sada se svodi na traženje koeficijenata a_{k} i b_{k}. Primijetite da više eksplicitno ne moramo promatrati fazu pojedine sinus funkcije, već je dovoljno tražiti koeficijente (amplitude) uz sinus i kosinus.

2.1Uvjeti na funkciju f(x)

Vrijedi li jednakost (1) za svaku periodičku funkciju i je li uistinu istina da baš svaku periodičku funkciju f(x) možemo zapisati kao konačnu sumu sinusa i kosinusa različitih amplituda i frekvencija? Da bismo bolje razumjeli odgovor na ovako postavljeno pitanje, prisjetimo se elementarnih tvrdnji iz diferencijalnog računa kao što je tvrdnja da {\it zbroj svakih dviju neprekidnih funkcija mora opet biti neprekidna funkcija}, ili da je {\it zbroj dvije derivabilne funkcije opet derivabilna funkcija}, ili još jače tvrdnje da {\it zbroj svakih dviju m-puta derivabilnih funkcija, m\in\mathbb{N}, mora opet biti m-puta derivabilna funkcija.} I sinus i kosinus su i neprekidne i m-puta derivabilne funkcije pa vrijedi i da je svaka njihova konačna suma takva. Iako je dovoljno uzeti uvjet derivabilnosti, htjeli smo posebno naglasiti i slabije svojstvo neprekidnosti. Uzevši u obzir ove činjenice, čini se da je nemoguće razviti nederivabilnu funkciju u ovako definiran Fourierov red!

No, stvar prima drugačiju perspektivu kada iskoristimo činjenicu da veće frekvencije sinusa i kosinusa u sumi (1) uzrokuju oštrije rubove u rezultatnoj krivulji. Primjer toga prikazan je na slici 2 za tzv. rectangle funkciju koja je ovdje definirana na sljedeći način:
f(x)= \begin{cases} 1, & x\in(-1,1), \\ 0, & \text{inače.} \\ \end{cases}.




Slika 2: Na slici je prikazana aproksimacija tzv. rectangle funkcije na intervalu [-2,2] s pomoću Fourierova reda za: a) n=2; b) n=4; c) n=8; d) n=16.


Stoga, kada bismo dodavali sve više i više funkcija, brojač k bio bi sve veći i veći i naša rezultatna krivulja sve bi bolje aproksimirala vrhove u kojima funkcija nije derivabilna. Ovaj argument iskoristiti ćrmo za motivaciju da sumu (1) preoblikujemo tako da ide u beskonačnost, na sljedeći način:
(2)
f(x)=\frac{a_{0}}{2}+\sum_{k=1}^{\infty}{(a_{k}\sin{(k x)}+b_{k}\cos{(k x)})}.

I ova se suma može bolje zapisati, prisjetimo li se Eulerove formule
r e^{i\phi}=r(\cos{\phi}+i\sin{\phi}).


Formulu (2) možemo zapisati na nov način tako da uvedemo nove koeficijente c_{k}\in\mathbb{C}, k\in\mathbb{Z}. No, želimo biti sigurni da ćemo na kraju dobiti realne koeficijente. Stoga ćemo iskoristiti svojstvo da je z+\overline{z}=2Re(z), gdje je z\in\mathbb{C}, a \overline{z} oznaka za konjugirano kompleksni broj. Primijetimo da je s pomoću Eulerove formule, da bismo dobili konjugirano kompleksni broj, dovoljno za \phi uzeti -\phi (svojstvo neparnosti funkcije sinus). Imajući to na umu, želimo pronaći koeficijente c_{k} takve da vrijedi
c_{k} e^{ikx} + \overline{ c_{k}e^{ikx}} = a_{k} \sin(kx) + b_{k} \cos(kx),
što je ekvivalentno
(\cos{(k x)}-i\sin{(k x)})\overline{c_{k}}+(\cos{(k x)}+i\sin{(k x)})c_{k}=a_{k}\sin{(k x)}+b_{k}\cos{(k x)}.


Nadalje, vrijedi
\cos{(k x)}(\overline{c_{k}}+c_{k})+\sin{(k x)}(c_{k} i-\overline{c_{k}} i)=a_{k}\sin{(k x)}+b_{k}\cos{(k x)}.


Primijetimo da je c_{k} i-\overline{c_{k}}i=-2Im(c_{k}). Sada iz gornjeg izraza dobivamo sljedeće dvije jednadžbe:
\begin{eqnarray*} a_{k}&=&-2Im(c_{k}),\\ b_{k}&=&2Re(c_{k}). \end{eqnarray*}


Ako prvu jednadžbu pomnožimo s imaginarnom jedinicom i, tj. a_{k} i=-2Im(c_{k})i te joj pribrojimo drugu jednadžbu, dobivamo
\begin{eqnarray*} b_{k}+a_{k} i&=&2Re(c_{k})-2Im(c_{k})i,\\ \overline{c_{k}}&=&\frac{1}{2}(b_{k}+i a_{k}). \end{eqnarray*}
Primijetite da smo se gore koristili jednakošću Re(c_{k})-Im(c_{k})i=\overline{c_{k}}. Nadalje, sada je lako dobiti samu vrijednost c_{k}.

Ako uvedemo oznaku c_{-k}=\overline{c_{k}}, općenitiji razvoj u Fourierov red možemo zapisati sljedećom formulom
f(x)=\sum_{k=-\infty}^{\infty}{c_{k} e^{i k x}},


gdje je c_{-k}=\frac{1}{2}(b_{k}+i a_{k}) i c_{k}=\frac{1}{2}(b_{k}-i a_{k}). Nulti član, „DC komponenta”, u ovom je slučaju c_{0}=\overline{c_{0}}=\frac{a_{0}}{2}.

2.2Određivanje Fourierovih koeficijenata c_{k}

U ovom poglavlju pokazat ćemo kako za danu periodičku funkciju f(x) izračunati k-ti koeficijent c_{k}. Budući da je
f(x)=\ldots+c_{k-1}e^{(k-1)x i}+c_{k} e^{k x i}+c_{k+1}e^{(k+1)x i}+\ldots,


dijeleći ovu jednadžbu s e^{k x i} (koji je uvijek različit od nule) dobivamo
\begin{eqnarray*} c_{k}&=&\frac{f(x)}{e^{k x i}}-(\ldots+c_{k-1}e^{-x i}+c_{k+1}e^{x i}+\ldots)\\ &=&f(x)e^{-k x i}-\sum_{j=-\infty,j\neq k}^{\infty}{c_{j}e^{(j-k) x i}}, \end{eqnarray*}


za k\in\mathbb{Z}. Sada iskoristimo činjenicu da je c_{k} konstanta, tako da integrirajući obje strane po varijabli x na [-\pi,\pi] dobivamo
\begin{eqnarray*} \int_{-\pi}^{\pi}{c_{k} dx}&=&\int_{-\pi}^{\pi}{f(x)e^{-k x i}dx}-\int_{-\pi}^{\pi}{\sum_{j=-\infty,j\neq k}^{\infty}{c_{j}e^{(j-k) x i}}dx},\\ c_{k}\pi-(-c_{k}\pi)&=&\int_{-\pi}^{\pi}{f(x)e^{-k x i}dx}-\sum_{j=-\infty,j\neq k}^{\infty}{\int_{-\pi}^{\pi}{c_{j}e^{(j-k) x i}dx}},\\ 2\pi c_{k}&=&\int_{-\pi}^{\pi}{f(x)e^{-k x i}dx}-\sum_{j=-\infty,j\neq k}^{\infty}{\frac{1}{(j-k)i}c_{j}(e^{(j-k) \pi i}-e^{-(j-k) \pi i})}, \end{eqnarray*}


za k\in\mathbb{Z}. Čini se da smo dobili još kompliciraniji izraz za c_{k}, ali ako primijetimo da je, po Eulerovoj formuli, te jer je (j-k)\in\mathbb{Z},
e^{(j-k)\pi i}-e^{-(j-k)\pi i}=i\sin{(j-k)\pi}=0,


vidimo da se cijela suma poništava. Dijeleći jednadžbu s 2\pi dobivamo formulu za k-ti koeficijent Fourierova reda:
c_{k}=\frac{1}{2\pi}\int_{-\pi}^{\pi}{f(x)e^{-k x i}dx}.


Slično se može pokazati da vrijedi i za općenitiji interval, odnosno period, [-L/2,L/2], te da tada Fourierov red glasi
(3)
f(x)=\sum_{k=-\infty}^{\infty}{c_{k} e^{(2\pi i k x)/L}},

s pridruženim koeficijentima
(4)
c_{k}=\frac{1}{L}\int_{-L/2}^{L/2}{f(x)e^{-(2\pi k x i)/L}dx},

gdje je k\in\mathbb{Z}. Primijetimo da koeficijent \frac{2\pi}{L} služi da bismo sveli sinus i kosinus funkcije na interval [0, L] s periodom L. Time je možemo promatrati i na [-L/2, L/2] s periodom L.

2.3Konvergencija Fourierovog reda

Ostaje nam samo pitanje konvergencije, koje je stoljećima predstavljalo poveći problem. Razlog tomu je alternirajuća priroda sinus i kosinus funkcije; tim svojstvom znatno je otežano promatranje konvergencije ovako definiranog reda. No, stvar otežava i promatranje Fourierove tvrdnje - da se ovako može razviti svaka, pa i nederivabilna funkcija! Rezultati promatranja ovog reda govore da je najvažniji uvjet zapravo integrabilnost funkcije koju razvijamo u Fourierov red. Što se tiče konvergencije, navest ćemo jedan važniji rezultat, a to je teorem koji govori o konvergenciji po L^{2} normi:

Teorem 1. [Riesz-Fischer] Funkcija f\in L^{2}([-P/2,P/2]) s periodom P je L^{2}-integrabilna, tj. vrijedi
\int_{-P/2}^{P/2}{|f(x)|^{2}dx}\lt \infty,


ako i samo ako je Fourierov red te funkcije L^{2}-konvergentan, tj. ako
\int_{-P/2}^{P/2}{|f(x)-\sum_{k=-n}^{n}{c_{k} e^{(2\pi i k x)/P}}|^{2}}dx\rightarrow 0,

kada n\rightarrow\infty.



Dakle, kvadrat razlike funkcije f(x) i njene konačne aproksimacije na nekom intervalu teži u nulu. Ovakve rezultate o konvergenciji Fourierova reda dobivamo ako se koristimo općenitijim, odnosno Lebesgueovim integralom. U analizi signala L^{2} integral funkcije, odnosno signala, predstavlja energiju signala. Ako je taj integral konačan, kaže se da signal ima konačnu energiju. Ovo rješenje postignuto je tek 1907. godine, što je cijelo stoljeće od prvotnog utemeljenja Fourierove analize. Ovaj teorem usko je vezan s Parsevalovim teoremom, pa se nerijetko tako i naziva.

3Fourierova transformacija


Promatrali smo kako razviti periodičnu funkciju perioda L definiranu na intervalu [-L/2,L/2] u Fourierov red. No, ako pokušamo pustiti da L teži u beskonačnost, tada je funkcija periodična na intervalu (-\infty,\infty) s beskonačno velikim periodom (pretpostavljamo da je funkcija definirana svugdje). Ovakvim razmišljanjem čini se da možemo dobiti koeficijente i za funkciju koja je potpuno neperiodična te na taj način proširiti teoriju koja vrijedi za Fourierov red i na neperiodične funkcije. Međutim, ima li to smisla? Primijetimo da ako postavimo L\rightarrow\infty, ako integral ima konačnu vrijednost na (-\infty,\infty), izraz (4) će otići u nulu (zbog 1/L ispred integrala). Stoga se čini da naivan pristup u kojemu želimo simulirati neperiodičnost s beskonačno dugim periodom ipak nema smisla, budući da će vrijednost koeficijenata težiti u nulu.

Međutim, ovakvo razmišljanje vodit će nas prema transformaciji izraza (4) za k-ti Fourierov koeficijent ka poznatoj formuli za Fourierovu transformaciju.

Kao prvi korak u toj transformaciji zapisat ćemo formulu (4) u obliku funkcije F ovisne o varijabli k/L na sljedeći način:
(5)
F\left(\frac{k}{L}\right)=\int_{-L/2}^{L/2}{f(x)e^{-(2\pi k x i)/L}dx}.

Tada nova formula za Fourierov red glasi:
(6)
f(x)=\sum_{k=-\infty}^{\infty}{\frac{1}{L}F\left(\frac{k}{L}\right) e^{(2\pi i k x)/L}}.

Sada pustimo da L\rightarrow\infty. Primijetimo da sada varijabla k/L, koja je bila diskretna, prelazi u kontinuiranu. Iako bismo prvo mislili da za fiksni k ovaj izraz prelazi u nulu, ideja je da je -\infty\lt k\lt \infty. Stoga se povećavanjem broja L zapravo samo smanjuje „razmak” između vrijednosti koje varijabla prima. Drugim riječima, varijabla k/L prelazi iz diskretne u kontinuiranu varijablu s. Sada zapišemo (5) na novi način:
F(s)=\int_{-\infty}^{\infty}{f(x)e^{-2\pi s x i}dx}.


Dobivena funkcija F:\mathbb{C}\rightarrow\mathbb{C} naziva se Fourierovom transformacijom funkcije f:\mathbb{C}\rightarrow\mathbb{C}; budući da izrazi za F i f sadržavaju kompleksne brojeve, i jedan i drugi daju kompleksne vrijednosti, a samim time i primaju. Promotrimo što se dogodi s formulom za Fourierov red. Budući da smo stavili L\rightarrow\infty, možemo primijetiti da suma u (6) prelazi u integral; 1/L zapravo predstavlja beskonačno malu vrijednost ds:
f(x)=\int_{-\infty}^{\infty}{F(s) e^{2\pi i s x}ds}.


U ovom izrazu f(x) predstavlja inverznu Fourierovu transformaciju - transformaciju koja temeljeno na koeficijentima F(s) za svaki s\in\mathbb{R} vraća originalnu funkciju.

Bitno je primijetiti da smo Fourierovu transformaciju dobili iz jednadžbe za koeficijente Fourierova reda, dok smo inverznu Fourierovu transformaciju dobili iz same jednadžbe Fourierova reda. U sljedećem primjeru pokazat ćemo Fourierovu transformaciju za rectangle funkciju.

Primjer 2. Neka je zadana funkcija izrazom
f(x)=\left\lbrace \begin{array}{ll} 1, & \text{}x\in(-\frac{1}{2},\frac{1}{2})\text{;} \\ 0, & \text{inače.} \\ \end{array} \right.


Ovakva se funkcija, kao što je prije napomenuto, često naziva rectangle funkcija (vidi sliku 3) i ima bitnu primjenu pri obradi signala. Nekada se definira i vrijednost od \frac{1}{2} u točkama \frac{1}{2} i -\frac{1}{2}, ali račun ostaje isti. Na\dj imo Fourierovu transformaciju ove funkcije.


Slika 3: Graf rectangle funkcije.


Rješenje 3. Fourierova transformacija ove funkcije dobiva se preko već opisane formule:
F(s)=\int_{-\infty}^{\infty}{f(x)e^{-2\pi s x i}dx}.


Da bismo mogli izračunati ovaj integral za zadani f(x), moramo ga rastaviti na tri dijela i uvrstiti odgovarajuće slučajeve za f(x) na sljedeći način:
F(s)=\int_{-\infty}^{-\frac{1}{2}}{0\cdot e^{-2\pi s x i}dx}+\int_{-\frac{1}{2}}^{\frac{1}{2}}{1\cdot e^{-2\pi s x i}dx}+\int_{\frac{1}{2}}^{\infty}{0\cdot e^{-2\pi s x i}dx}.


Iz gornjeg izraza vidimo da iščezavaju svi integrali osim srednjeg te je dovoljno pronaći:
F(s)=\int_{-\frac{1}{2}}^{\frac{1}{2}}{e^{-2\pi s x i}dx}.


Tada po Eulerovoj formuli slijedi:
\begin{eqnarray*} F(s)&=&\int_{-\frac{1}{2}}^{\frac{1}{2}}{(\cos{(-2\pi s x)}+i \sin{(-2\pi s x)})dx}\\ &=&\int_{-\frac{1}{2}}^{\frac{1}{2}}{\cos{(-2\pi s x)}dx}+i \int_{-\frac{1}{2}}^{\frac{1}{2}}{\sin{(-2\pi s x)}dx}\\ &=&-\frac{1}{2\pi s}\left(\sin{\left(-2\pi s\cdot \frac{1}{2}\right)}-\sin{\left(-2\pi s\cdot\left(-\frac{1}{2}\right)\right)}\right)\\ &&-\frac{i}{2\pi s}\left(\cos{\left(-2\pi s \cdot \frac{1}{2}\right)}-\cos{\left(-2\pi s \cdot\left(-\frac{1}{2}\right)\right)}\right). \end{eqnarray*}


Koristeći se svojstvom parnosti funkcije kosinus i neparnosti funkcije sinus dobivamo
\begin{eqnarray*} F(s)&=&-\frac{1}{2\pi s}(-\sin{(\pi s)}-\sin{(\pi s)})-\frac{i}{2\pi s}(\cos{(\pi s)}-\cos{(\pi s)})\\ &=&-\frac{1}{2\pi s}(-2\sin{(\pi s)})\\ &=&\frac{\sin{(\pi s)}}{\pi s}. \end{eqnarray*}


Gornjim izrazom dana je Fourierova transformacija tražene funkcije.


{\bf Napomena.} Funkcija dobivena Fourierovom transformacijom rectangle funkcije često se u analizi signala naziva sinc funkcija (vidi sliku 4) i definira na sljedeći način:
\text{sinc}(x)=\frac{\sin{(\pi x)}}{\pi x}.


Slika 4: Graf sinc funkcije.


3.1Diskretna Fourierova transformacija (DFT)

Ako imamo diskretan skup kompleksnih brojeva (koji smo mogli dobiti i uzorkovanjem funkcije na nekom intervalu), Fourierova transformacija koja je definirana za kontinuiran skup nam ne odgovara. Tada definiramo DFT na sljedeći način:
(7)
F_{j} = \sum_{k=0}^{n-1} f_{k} e^{\frac{2\pi ikj}{n}},
koja niz f_{0}, \ldots , f_{n-1} kompleksnih brojeva pretvara u niz F_{0}, \ldots , F_{n-1} kompleksnih brojeva s pripadajućim inverzom
(8)
f_{k} =\frac{1}{n}\sum_{j=0}^{n-1}F_{j} e^{\frac{-2\pi ikj}{n}},
koja kompleksne brojeve F_{0}, \ldots , F_{n-1} pretvara u kompleksne brojeve f_{0}, \ldots , f_{n-1}.

Primijetite da su u izrazima (7) i (8) diskretne vrijednosti u kojima evaluiramo funkciju f(x), npr. x_{0},\ldots,x_{n-1}, te funkciju F(s), npr. s_{0},\ldots,s_{n-1}, potpuno zamijenjene indeksima k = 0,\ldots,n-1 i j=0,\ldots,n-1 na uzorku veličine n na način da definiramo F(s_{j})=F_{j} te f(x_{k})=f_{k}. Nadalje, u eksponentu Eulerove konstante e umnožak dvaju odgovarajućih uzoraka x_{k}s_{j} bit će jednak kj/n, budući da bez smanjenja općenitosti možemo pretpostaviti da uzorkujemo funkcije f(x) s nekog intervala uzorkovanja oblika [0,L] duljine L s uzorcima x_{k}=k/B, te F[s] s nekog intervala uzorkovanja oblika [0,B] duljine B s uzorcima s_{j}=j/L. Uzorci x_{k} i s_{j} uzimaju se tako da vrijedi
n\cdot \frac{1}{B} = L \text{ tj. } n\cdot \frac{1}{L} = B.

Budući da su ovako definirane dvije Fourierove transformacije inverz jedna drugoj, nije nam nužno strogo definirati jednu od tih transformacija izričito kao inverz. Stoga se u literaturi (kao i za općenitu Fourierovu transformaciju), uglavnom radi jednostavnosti, inverzi označavaju različito. Tako se (8) u analizi signala češće koristi kao DFT, a (7) kao inverz (uz promjenu predznaka u eksponentu od e). To ne predstavlja problem, jer se ionako svi računi obavljaju slično.

4Zaključak i diskusija

Fourierov red, a posebice Fourierova transformacija u današnje su vrijeme jedan od najčešće korištenih matematičkih alata, a time i predmet brojnih šala na račun matematičara i inženjera (vidi sliku 5).


Slika 5: Slika preuzeta s http://imgs.xkcd.com/comics/fourier.jpg.


Značenje Fourierove transformacije je višestruko i ovisno o primjeni. Česta je primjena kada se funkcija - signal, koja je definirana na cijelom realnom pravcu, pretvara u funkciju koja za određenu frekvenciju k vraća F(s), vrijednost iz koje možemo saznati informacije o amplitudi i fazi za frekvenciju s. Mnogo je operacija lakše obavljati preko Fourierove transformacije funkcije, među njima je najbitnija konvolucija. Jednako tako, Fourierova transformacija ima puno svojstava koja olakšavaju analiziranje funkcije - signala (vidi [1], [5]).

Ovaj rad je rezultat uvodne diskusije diplomskog rada „Brza Fourierova transformacija i primjene” studenta S. Poljaka napisanog pod mentorstvom D. Matijevića.



Bibliografija
[1] Brigham, E. Oran. The fast Fourier transform and its applications, Englewood Cliffs, N.J.: Prentice Hall. ISBN 0-13-307505-2, 1988
[2] Cormen, Thomas H.; Leiserson, Charles E. ; Rivest, Ronald L.; and Stein, Clifford. "Chapter 30: Polynomials and the FFT". Introduction to Algorithms (Second ed.), MIT Press and McGraw-Hill. pp. 822–848 ISBN 0-262-03293-7, 2001
[3] Oppenheim, Alan V.; Schafer, R. W.; and Buck, J. R., Discrete-time signal processing, Upper Saddle River, N.J.: Prentice Hall. ISBN 0-13-754920-2, 1999
[4] Osgood, Brad G. The Fourier Transform and its Applications, on-line lectures by Stanford School of Engineering, at http://academicearth.org/courses/the-fourier-transform-and-its-applications.
[5] Poljak, Stjepan. Brza Fourierova transformacija i primjene, diplomski rad, Odjel za matematiku Sveučilište J.J. Strossmayer u Osijeku, Java implementacija: http://www.mathos.hr/FFT/fourier.html, 2010.
[6] Tolstov, Georgi P. Fourier Series. Courier-Dover. ISBN 0486633179, 1976

 

Diferenciranje i integriranje pod znakom integrala

Irfan Glogić , Harun Šiljak

When guys at MIT or Princeton had trouble doing a certain integral, it was because they couldn't do it with the standard methods they had learned in school. If it was contour integration, they would have found it; if it was a simple series expansion, they would have found it. Then I come along and try differentiating under the integral sign, and often it worked. ([1])


1Teorijski uvod

Diferenciranje i integriranje pod znakom integrala je tehnika koja je često korisna u izračunavanju integrala funkcija jedne realne varijable. Prije nego što krenemo s primjerima, navedimo osnovne teoreme kojima ćemo se koristiti.

Teorem 1. Neka je funkcija f(x,y) definirana na pravokutniku [a,b]\times[c,d] i neka je neprekidna po x na [a,b] za proizvoljan y. Pretpostavimo također da postoji parcijalna derivacija \frac{\partial}{\partial y}f(x,y) i da je neprekidna kao funkcija dviju varijabli. Tada za svaki y\in[c,d] vrijedi

\frac{d}{dy}(\int_{a}^{b}f(x,y)\, dx)=\int_{a}^{b}\frac{\partial}{\partial y}f(x,y)\, dx.


Što se tiče diferenciranja pod znakom integrala koji je neodređen, osobitu ulogu ima pojam uniformne konvergencije integrala.

Naime, ako postoji integral I(y)=\int_{a}^{\infty}f(x,y)\, dx (definiran kao \lim_{b\to\infty}\int_{a}^{b}f(x,y)dx) za y\in Y i za svaki \varepsilon\gt 0 postoji b_{0}\ge a koji ne ovisi o y, takav da za b\gt b_{0} vrijedi

\left|\int_{a}^{\infty}f(x,y)\, dx-\int_{a}^{b}f(x,y)\, dx\right|=\left|\int_{b}^{\infty}f(x,y)\, dx\right|\lt \varepsilon\text{ za sve } y\in Y,

tada kažemo da integral I(y) konvergira uniformno po y\in Y.

Za dokazivanje uniformne konvergencije integrala koriste se razni kriteriji. Mi ćemo navesti dva.



Kriterij 2. Pretpostavimo da je funkcija f(x,y) integrabilna po x na svakom konačnom segmentu [a,\eta] (\eta\ge a). Ako postoji funkcija \varphi(x) koja ovisi samo o x, integrabilna na [a,\infty\rangle takva da za svaki y\in Y vrijedi |f(x,y)|\le\varphi(x) (za x\ge a), onda integral I(y) konvergira uniformno po y.

Kriterij 3. Ako je integral \int_{a}^{\infty}f(x)\, dx konvergentan, a funkcija g(x,y) monotona po x i uniformno ograničena, onda integral \int_{a}^{\infty}f(x)g(x,y)\, dx konvergira uniformno po y.


Sljedeći teorem daje dovoljne uvjete za prijelaz limesa pod znak integrala.

Teorem 4. Neka je funkcija f(x,y) za y\in Y integrabilna po x na segmentu [a,A] za sve A\gt a i neka na svakom segmentu konvergira uniformno po x graničnoj funkciji \varphi(x) kada y\to y_{0}. Ako pored toga integral I(y)=\int_{a}^{\infty}f(x,y)\, dx konvergira uniformno po y\in Y, onda vrijedi

\lim_{y\to y_{0}}\int_{a}^{\infty}f(x,y)\, dx=\int_{a}^{\infty}\varphi(x)\, dx,

gdje y_{0} može biti i \infty.

Što se tiče diferenciranja pod znakom neodređenog integrala, pokazuje se da i u ovom slučaju vrijedi Leibnizovo pravilo.

Teorem 5. Neka je funkcija f(x,y) definirana i neprekidna po x za x\ge a i y iz segmenta [c,d] te neka za x\ge0 i y\in[c,d] ima derivaciju \frac{\partial}{\partial y}f(x,y) koja je neprekidna funkcija po obje varijable. Pretpostavimo također da integral I(y)=\int_{a}^{\infty}f(x,y)\, dx konvergira za sve y\in[c,d], a integral \int_{a}^{\infty}\frac{\partial}{\partial y}f(x,y)\, dx konvergira uniformno po y na tom istom segmentu. Tada za proizvoljni y\in[c,d] vrijedi

I'(y)=\int_{a}^{\infty}\frac{\partial}{\partial y}f(x,y)\, dx.


Za potrebe integracije pod znakom integrala navodimo sljedeće teoreme.

Teorem 6. Ako je funkcija f(x,y) neprekidna na pravokutniku [a,b]\times[c,d], tada vrijedi

\int_{c}^{d}\int_{a}^{b}f(x,y)\, dx\, dy=\int_{a}^{b}\int_{c}^{d}f(x,y)\, dy\, dx.


Teorem 7. Neka je funkcija f(x,y) definirana i neprekidna za x\ge a i y\in[c,d]. Ako integral I(y)=\int_{a}^{\infty}f(x,y)\, dx konvergira uniformno po y na segmentu [c,d], tada vrijedi

\int_{c}^{d}I(y)\, dy=\int_{c}^{d}\int_{a}^{\infty}f(x,y)\, dx\, dy=\int_{a}^{\infty}\int_{c}^{d}f(x,y)\, dy\, dx.


Teorem 8. Neka je funkcija f(x,y) definirana i neprekidna za x\ge a i y\ge c. Pretpostavimo također da oba integrala \int_{a}^{\infty}f(x,y)\, dx i \int_{c}^{\infty}f(x,y)\, dy konvergiraju uniformno, prvi po y, a drugi po x. Tada, ako postoji bar jedan od integrala \int_{c}^{\infty}\int_{a}^{\infty}|f(x,y)|\, dx\, dy i \int_{a}^{\infty}\int_{c}^{\infty}|f(x,y)|\, dy\, dx, onda postoje i jednaki su integrali

\int_{c}^{\infty}\int_{a}^{\infty}f(x,y)\, dx\, dy\text{ i }\int_{a}^{\infty}\int_{c}^{\infty}f(x,y)\, dy\, dx.


Dokazi navedenih tvrdnji lako se mogu izvesti iz rezultata iz naprimjer [2, 3].

2Riješeni primjeri

U sljedećim primjerima prikazana je primjena diferenciranja i integriranja pod znakom integrala. Pritom su neki zadaci detaljno riješeni dok su ponegdje dijelovi rješenja ostali neprikazani - čitatelju se savjetuje da sam pokuša napraviti potrebne dopune.

Ako nije naveden izvor zadatka, on je zajedno s rješenjem preuzet iz [3].

Primjer 9.([4]) Nađite vrijednost integrala

\int_{0}^{1}\frac{\ln(1+x)}{1+x^{2}}\, dx.

Dokaz. Definirajmo funkciju f(x,t)=\frac{\ln(xt+1)}{x^{2}+1}. Ova funkcija neprekidna je na pravokutniku P=[0,1]\times[0,1] i parcijalna derivacija \frac{\partial}{\partial t}f(x,t)=\frac{x}{(xt+1)(x^{2}+1)} postoji i neprekidna je na P. Tada prema Teoremu 1 za svaki t\in[0,1] i I(t)=\int_{0}^{1}f(x,t)\, dx vrijedi I'(t)=\int_{0}^{1}\frac{x\, dx}{(xt+1)(x^{2}+1)}. Rastavljanjem na parcijalne razlomke i integracijom dobijamo
I'(t)=\frac{2t\mathop{\text{arctg}} x-2\ln(tx+1)+\ln(x^{2}+1)}{2(t^{2}+1)}\Big|_{x=0}^{x=1}=\frac{\pi t+2\ln2-4\ln(t+1)}{4(t^{2}+1)}.
Odavde slijedi da je
I(t)=\frac{\ln2\mathop{\text{arctg}} t}{2}+\frac{\pi\ln(t^{2}+1)}{8}-\int_{0}^{t}\frac{\ln(t+1)}{t^{2}+1}dt,
pa je
I(1)=\frac{\pi\ln2}{4}-\int_{0}^{1}\frac{\ln(t+1)}{t^{2}+1}dt.
Budući da je integral na desnoj strani jednakosti upravo traženi integral I(1), slijedi da je tražena vrijednost f(1)=\frac{\pi\ln2}{8}.
\ \blacksquare

Primjer 10.([5]) Nađite vrijednost integrala

\int_{0}^{\infty}\frac{\mathop{\text{arctg}}\pi x-\mathop{\text{arctg}} x}{x}\, dx.

Dokaz. Neka je f(x,t)=\frac{\mathop{\text{arctg}} tx-\mathop{\text{arctg}} x}{x}. Ova funkcija je neprekidna za x\ge0 i t\in[1,\pi], pri čemu je na istom skupu neprekidna i derivacija \frac{\partial}{\partial t}f(x,t)=\frac{1}{1+t^{2}x^{2}}. Integral I(t)=\int_{0}^{\infty}f(x,t)\, dx konvergira za sve t\in[1,\pi], što je moguće pokazati npr. Lagrangeovim teoremom o srednjoj vrijednosti, jer iz njega slijedi da za t\gt 1 vrijedi

\frac{\mathop{\text{arctg}} tx-\mathop{\text{arctg}} x}{x(t-1)}=(\mathop{\text{arctg}} x)'_{x=c},

gdje je c neki broj iz segmenta [x,tx], pa imamo

\frac{\mathop{\text{arctg}} tx-\mathop{\text{arctg}} x}{x(t-1)}=\frac{1}{1+c^{2}}\le\frac{1}{1+x{}^{2}}.

Odavde je konvergencija integrala očita. Sada možemo, koristeći se Kriterijem 1, pokazati da \int_{0}^{\infty}\frac{\partial}{\partial t}f(x,t)\, dx konvergira uniformno po t, jer je \frac{1}{1+t^{2}x^{2}}\le\frac{1}{1+x^{2}} za t\ge1. Tada prema Teoremu 5 vrijedi

I'(t)=\int_{0}^{\infty}\frac{\partial}{\partial t}f(x,t)\, dx=\int_{0}^{\infty}\frac{1}{x}\left(\frac{x}{1+t^{2}x^{2}}\right)dx=\frac{1}{t}\int_{0}^{\infty}\frac{dy}{1+y^{2}}=\frac{\pi}{2t}.

Stoga je I(t)=\frac{\pi}{2}\ln t+C. Budući da je I(1)=0, slijedi C=0. Dakle, I(\pi)=\frac{\pi}{2}\ln\pi.
\ \blacksquare

Primjer 11. Nađite vrijednost integrala

\int_{0}^{\infty}\frac{\sin x}{x}\, dx.

Dokaz. Neka je f(x,t)=\frac{\sin x}{x}e^{-tx}. Ova funkcija je neprekidna na P=[0,\infty\rangle\times[0,a] za svako a\gt 0, pri čemu je na istom skupu neprekidna i derivacija \frac{\partial}{\partial t}f(x,t)=-\sin x\, e^{-tx}. Integral I(t)=\int_{0}^{\infty}f(x,t)\, dx konvergira za svako t jer za t=0 imamo \int_{0}^{\infty}\frac{\sin x}{x}e^{-tx}\, dx=\int_{0}^{\infty}\frac{\sin x}{x}\, dx\lt \infty (posljednja nejednakost slijedi npr. iz Dirichletova kriterija konvergencije, jer \left|\int_{0}^{a}\sin x\, dx\right|\lt 2 za svaki a, dok \frac{1}{x} monotono opada k nuli za x\gt 0), a za t\ge t_{0}\gt 0 imamo \left|\int_{0}^{\infty}\frac{\sin x}{x}e^{-tx}\, dx\right|\le\int_{0}^{\infty}e^{-tx}\, dx=\frac{1}{t}\lt \infty.

Za t\gt 0 vrijedi \left|\frac{\partial}{\partial t}f(x,t)\right|=\left|-\sin x\, e^{-tx}\right|\le e^{-tx}\le e^{-t_{0}x}=\varphi(x), a budući da je \varphi(x) integrabilna na [0,\infty\rangle, po Kriteriju 1 integral \int_{0}^{\infty}\frac{\partial}{\partial t}f(x,t)\, dx konvergira uniformno na svakom skupu oblika \lbrace t\in\mathbb{R}|t\ge t_{0}\gt 0\rbrace. Prema Teoremu 5 vrijedi

I'(t)=-\int_{0}^{\infty}\sin x\, e^{-tx}\, dx.

Primjenom parcijalne integracije dva puta dobijamo

I'(t)=-\frac{1}{1+t^{2}}\text{ za }t\gt 0.

Odatle slijedi da je I(t)=-\mathop{\text{arctg}} t+C za t\gt 0. Pokažimo da je \lim_{t\to\infty}I(t)=0.

Neka je \varepsilon\gt 0 proizvoljan i odaberimo x_{0}\gt 0 takvav da je \sin x\ge0 za x\in[0,x_{0}] i
0\lt \int_{0}^{x_{0}}\frac{\sin x}{x}e^{-tx}\, dx\lt \int_{0}^{x_{0}}\frac{\sin x}{x}\lt \frac{\varepsilon}{2}

za sve t\gt 0. Funkcija f(x,t) je integrabilna na segmentu [x_{0},A] za svaki A\gt x_{0} i na svakom takvom segmentu konvergira uniformno po x k 0 kada t\to\infty. Integral \int_{x_{0}}^{\infty}\frac{\sin x}{x}e^{-tx}\, dx konvergira uniformno po t\gt 0 (prema Kriteriju 2) pa koristeći se Teoremom 4 dobijamo \lim_{t\to\infty}I(t)=0. Dakle, C=\frac{\pi}{2}, a vrijednost traženog integrala je I(0)=\frac{\pi}{2}.

Dokaz da je \lim_{t\to\infty}I(t)=0 mogao je biti izveden i jednostavnije, kao što ćemo vidjeti u sljedećem primjeru.
\ \blacksquare

Primjer 12. Dokažite jednakost

\int_{0}^{\infty}\frac{2-2\cos x}{xe^{x}}\, dx=\ln2.

Dokaz. Neka je f(x,t)=\frac{2-2\cos x}{xe^{x}}e^{-tx}. Koristeći se idejama iz prethodnih zadataka, može se pokazati da integral I(t)=\int_{0}^{\infty}f(x,t)\, dx konvergira za sve t\in\lbrace t\in\mathbb{R}|t_{0}\le t\lt \infty\rbrace gdje je 0\lt t_{0}\lt 1 te da na istom skupu integral \int_{0}^{\infty}f(x,t)\, dx konvergira uniformno. Stoga prema Teoremu 5 imamo
I'(t)=\int_{0}^{\infty}f'_{t}(x,t)\, dx=\int_{0}^{\infty}(2-2\cos x)e^{-tx}dx=-\frac{2}{t}+\frac{2t}{t^{2}+1},
I(t)=\ln\left(1+\frac{1}{t^{2}}\right)+C.
Primijetimo da je \frac{2-2\cos x}{x} neprekidna ograničena funkcija na [0,\infty\rangle i da je
\lim_{t\to\infty}\int_{0}^{\infty}e^{-tx}dx=0.
Stoga je \lim_{t\to\infty}I(t)=0, pa je C=0. Dakle, traženi integral je I(1)=\ln2.
\ \blacksquare

U nastavku slijede primjeri primjene integracije pod znakom integrala u rješavanju zadataka. Najprije ćemo ovu tehniku primijeniti na već dane primjere 10 i 12, a zatim i na jedan poznati rezultat iz analize.

Primjer 13.([5]) Nađite vrijednost integrala

\int_{0}^{\infty}\frac{\mathop{\text{arctg}}\pi x-\mathop{\text{arctg}} x}{x}\, dx.

Dokaz.
I=\int_{0}^{\infty}\frac{\mathop{\text{arctg}}\pi x-\mathop{\text{arctg}} x}{x}\, dx=\int_{0}^{\infty}\frac{1}{x}\mathop{\text{arctg}} tx\big|_{t=1}^{t=\pi}\, dx=\int_{0}^{\infty}\int_{1}^{\pi}\frac{1}{1+(xt)^{2}}dt\, dx.

Budući da je \frac{1}{1+(tx)^{2}} neprekidna funkcija za x\ge0 i t\in[1,\pi] i integral \int_{0}^{\infty}\frac{1}{1+(tx)^{2}}\: dx prema Kriteriju 1 konvergira uniformno za svako t\in[1,\pi], iz Teorema 7 slijedi opravdanost sljedećeg postupka:
I=\int_{0}^{\infty}\int_{1}^{\pi}\frac{1}{1+(xt)^{2}}dt\, dx=\int_{1}^{\pi}\int_{0}^{\infty}\frac{1}{1+(xt)^{2}}dx\, dt,
I=\int_{1}^{\pi}\frac{1}{t}\cdot\frac{\pi}{2}dt=\frac{\pi}{2}\ln\pi.
\ \blacksquare

Primjer 14. Dokažite jednakost

\int_{0}^{\infty}\frac{2-2\cos x}{xe^{x}}\, dx=\ln2.

Dokaz. Budući da je \int_{0}^{1}\sin xt\ dt=\frac{1-\cos x}{x}, imamo
I=\int_{0}^{\infty}\frac{2-2\cos x}{x}\cdot e^{-x}dx=2\int_{0}^{\infty}\int_{0}^{1}\sin xt\ e^{-x}dt\ dx.

Budući da je \sin xt\ e^{-x} neprekidna funkcija za x\ge0 i t\in[0,1] i integral \int_{0}^{\infty}\sin xt\ e^{-x}\, dx konvergira uniformno za svako t\in[0,1] po prvom kriteriju, iz Teorema 7 slijedi opravdanost sljedećeg postupka:

I=2\int_{0}^{1}\int_{0}^{\infty}\sin xt\ e^{-x}dx\ dt=2\int_{0}^{1}\frac{t}{1+t^{2}}dt=\ln2.
\ \blacksquare

Primjer 15. Nađite vrijednost integrala

\int_{0}^{\infty}e^{-x^{2}}\, dx.

Dokaz. Iako ovo rješenje ne predstavlja klasičnu integraciju pod znakom integrala, ipak ga vrijedi prikazati.

Neka je I=\int_{0}^{\infty}e^{-x^{2}}dx. Vrijedi
\int_{0}^{\infty}e^{-tx^{2}}dx=\frac{I}{\sqrt{t}}.
Množenjem zadnje jednakosti s e^{-t} i integriranjem obiju strana na intervalu [0,\infty\rangle dobijamo
2I^{2}=\int_{0}^{\infty}\int_{0}^{\infty}e^{-tx^{2}}dx\ e^{-t}dt=\int_{0}^{\infty}\int_{0}^{\infty}e^{-tx^{2}}e^{-t}dt\ dx.

Promjena reda integriranja opravdana je prema Teoremu 8, budući da funkcija e^{-t(x^{2}+1)} zadovoljava uvjete ovog teorema: uniformnu konvergenciju integrala \int_{0}^{\infty}e^{-t(x^{2}+1)}dx i \int_{0}^{\infty}e^{-t(x^{2}+1)}dt lako je pokazati, a egzistencija jednog od dva tražena uzastopna integrala slijedi iz sljedećeg računa:

2I^{2}=\int_{0}^{\infty}\frac{1}{1+x^{2}}dx=\frac{\pi}{2}.
Prema tome je I=\frac{\sqrt{\pi}}{2}.
\ \blacksquare

3Zadaci za samostalan rad

Čitatelj može predstavljene metode primijeniti na sljedeće zadatke koji su ostavljeni za vježbu.

Zadatak 16. Dokažite jednakost

\int_{0}^{\infty}\frac{\ln(1+ax)}{\cos x}dx=\frac{1}{2}\mathop{\text{arctg}} a\log(1+a^{2}).


Zadatak 17. Dokažite jednakost

\int_{0}^{\infty}\frac{\mathop{\text{arctg}}\sqrt{x^{2}+2}}{(x^{2}+1)\sqrt{x^{2}+2}}dx=\frac{5\pi^{2}}{96}.


Zadatak 18. Nađite vrijednost integrala

\int_{0}^{\infty}\frac{x^{p-1}}{1+x}\ln x\, dx


za 0\lt p\lt 1.

Zadatak 19. Dokažite jednakost

\int_{0}^{\infty}2\sec x\ln\left(\frac{1+\beta\cos x}{1+\alpha\cos x}\right)dx=\arccos^{2}\alpha-\arccos^{2}\beta.


Zadatak 20.([6]) Dokažite jednakost

\int_{0}^{\infty}\left(x-\frac{x^{3}}{2}+\frac{x^{5}}{2\cdot4}-\frac{x^{7}}{2\cdot4\cdot6}+\cdots\right)\left(1+\frac{x^{2}}{2^{2}}+\frac{x^{4}}{2^{2}\cdot4^{2}}+\frac{x^{6}}{2^{2}\cdot4^{2}\cdot6^{2}}+\cdots\right)dx=\sqrt{e}.


Zadatak 21. Nađite vrijednost integrala

\int_{-\infty}^{\infty}\frac{\cos nx}{(x^{2}+1)^{2}}dx.

Bibliografija
[1] Feynman, R. P., Surely You're Joking, Mr. Feynman by Richard P. Feynman, Bantam Books, 1989
[2] Pandžić, P., Tambača J., Integrali funkcija više varijabli, Sveučilišna skripta, Sveučilište u Zagrebu, 2011.
[3] Ungar Š., Matematička analiza 3, Sveučilište u Zagrebu, 2002.
[4] Mathlinks forum: www.mathlinks.ro
[5] Kedlaya, K. S., Ng, L., Solutions to the 63rd William Lowell Putnam Mathematical Competition, 2002
[6] Alexanderson G. L., Klosinski L. F., Larson L. C., The William Lowell Putnam Mathematical Competition Problems and Solutions: 1965-1984, MAA, 1985
[7] Kedlaya, K. S., Poonen B., Vakil R., The William Lowell Putnam Mathematical Competition 1985–2000: Problems, Solutions, and Commentary, MAA, 2002


Alternativna definicija limesa funkcije

 


Ozren Perše, Ana Zeman

 



Sažetak
U ovom preglednom radu prezentiramo alternativnu definiciju limesa funkcije, danu u radu B.M. Baishanski, arXiv:0805.3671. Pokazujemo da osnovna svojstva limesa funkcije lagano slijede iz te definicije, te da je ta definicija ekvivalentna tradicionalnoj „\epsilon-\delta” definiciji.

1Uvod

Tradicionalna „\epsilon-\delta” definicija limesa realne funkcije jedne realne varijable (vidi npr. [2], [3]) dosta je apstraktna za mnoge studente koji se prvi put s njom susreću. Iz tog se razloga u radu [1] predlaže alternativna definicija limesa funkcije, koja bi trebala biti intuitivnija za početnike. Naime, autor tog članka mišljenja je da je jednostavnije geometrijski interpretirati Definicije 1 i 8 (navedene u ovom radu) nego tradicionalnu „\epsilon-\delta” definiciju. Alternativna definicija limesa funkcije dobiva se kombinacijom dviju dobro poznatih činjenica, koje su u sadašnjoj aksiomatici jednostavne posljedice tradicionalne „\epsilon-\delta” definicije (vidi Definiciju 1). Jedna od tih činjenica (svojstvo (2) iz Definicije 1) u novoj aksiomatici postaje fundamentalno svojstvo limesa funkcije.

U ovom radu promatramo realne funkcije realne varijable i, slijedeći [1], definiramo \lim f(x) kada x teži \infty. Potpuno analogno može se definirati \lim f(x) kada x teži c+ i c-, za c \in \mathbb{R}, te kada x teži -\infty. Također, uz manje modifikacije, moguće je dati i definiciju beskonačnog limesa (tj. \lim f(x)=\pm \infty).

Limes funkcije možemo shvatiti kao preslikavanje s klase realnih funkcija u realne brojeve. Slijedeći [1], definiramo to preslikavanje i maksimalnu klasu funkcija na kojoj je to preslikavanje definirano, tako da je osnovno svojstvo limesa (svojstvo (2) iz Definicije 1) zadovoljeno. To radimo u tri koraka:
1) definiramo limes za monotone ograničene funkcije,
2) definiramo klasu konvergentnih funkcija (tj. dajemo novo značenje pojmu konvergencije - kasnije pokazujemo da se podudara s tradicionalnim pojmom konvergencije (vidi teoreme 16 i 17)),
3) proširujemo definiciju limesa na sve konvergentne funkcije.


Nadalje, pokazujemo da iz te definicije limesa jednostavno slijede uobičajena svojstva limesa (teorem 15), te da je ova definicija limesa ekvivalentna tradicionalnoj „\epsilon-\delta” definiciji (teoremi 16 i 17).

Za a \in \mathbb{R}, u ovom radu s (a,\infty) označujemo otvoreni interval \lbrace x \in \mathbb{R} | \ x\gt a \rbrace.



2Limes monotonih i ograničenih funkcija

U ovom poglavlju promatramo klasu funkcija koje su monotone i ograničene na nekoj okolini beskonačnosti:
BM(\infty) = \lbrace f| \text{ postoji } a \in \mathbb{R} \text{ takav da je } f \text{ monotona i ograničena na intervalu } (a,\infty)\rbrace
i definiramo pojam limesa na toj klasi.

Definicija 1. Kažemo da je preslikavanje
L: BM(\infty) \to \mathbb{R}
limes na BM(\infty) ako su zadovoljena sljedeća svojstva:
(1) ako je f(x) = c za sve x \in \mathbb{R}, tada L(f) = c,
(2) ako je L(f) \lt L(g), tada postoji a \in \mathbb{R} takav da je
f(x) \lt g(x) \text{ za } x \gt a.


Umjesto L(f) = \lambda, često ćemo se koristiti uobičajenom notacijom „\lim_{x \to \infty} f(x) =\lambda” ili „f(x)\to\lambda \text{ kada } x\to \infty ”. Ako je f rastuća (odnosno padajuća), pišemo f(x)\nearrow \lambda kada x \to \infty (odnosno f(x)\searrow \lambda kada x\to \infty).

Sada pokazujemo da postoji jedinstveno preslikavanje sa svojstvima iz Definicije 1.

Teorem 2. Postoji limes L na BM(\infty).

Dokaz. Ako je f rastuća i ograničena za x \gt a, definiramo L(f) = \sup \lbrace f(x)| \ x \gt a \rbrace, a ako je f padajuća i ograničena za x \gt a, definiramo L(f) = \inf \lbrace f(x)| \ x \gt a \rbrace. Očito, trebamo samo provjeriti je li zadovoljen uvjet (2) iz Definicije 1. Postoje četiri slučaja s obzirom na to jesu li funkcije f i g rastuće ili padajuće. Dokaz dajemo samo za slučaj kada je f padajuća za x \gt a' i g rastuća za x \gt a'' (preostala tri slučaja su slična). Neka je \gamma \in \mathbb{R} takav da je L(f)\lt \gamma \lt L(g). Budući da je L(f) = \inf \lbrace f(x)| \ x \gt a' \rbrace, L(g) = \sup \lbrace g(x)| \ x \gt a'' \rbrace, dobivamo da postoje b' i b'' iz \mathbb{R} takvi da f(b') \lt \gamma \lt g(b''). Budući da je f padajuća, a g rastuća, slijedi da je f(x) \lt g(x) za x \gt \max \lbrace b',b'' \rbrace.
\ \blacksquare

Teorem 3. Limes L na BM(\infty) je jedinstven.

Dokaz. Pretpostavimo da postoje dva limesa, L' i L'', te funkcija f iz BM(\infty) takva da je L'(f) različit od L''(f). Bez smanjenja općenitosti možemo pretpostaviti da je L'(f) \lt L''(f). Neka je \gamma \in \mathbb{R} takav da je
(1)
\begin{eqnarray} L'(f) \lt \gamma \lt L''(f). \end{eqnarray}
Ako označimo s \gamma i konstantnu funkciju s vrijednosti \gamma, iz svojstva (1) iz Definicije 1 slijedi da je L'(\gamma) = L''(\gamma) = \gamma. Sada relacija (1) povlači da je L'(f) \lt L'(\gamma) i L''(\gamma) \lt L''(f), pa iz svojstva (2) slijedi da postoje a' i a'' iz \mathbb{R} takvi da je f(x) \lt \gamma za x\gt a' i f(x) \gt \gamma za x\gt a''. Dakle, za x \gt \max\lbrace a',a''\rbrace, dobivamo f(x) \lt f(x), što je kontradikcija.
\ \blacksquare

Teorem 4. Ako su funkcije f i g iz klase BM(\infty) i ako postoji a \in \mathbb{R} takav da
f(x)\le g(x) \text{ za } x\gt a, \text{ tada } L(f) \le L(g).

Dokaz. Pretpostavimo suprotno, odnosno L(f) \gt L(g). Iz svojstva (2) iz Definicije 1 slijedi da postoji a' \in \mathbb{R} takav da je f(x)\gt g(x) za x\gt a'. Sada za x \gt \max \lbrace a,a' \rbrace dobivamo kontradikciju.
\ \blacksquare

Teorem 5. Neka je N pozitivna padajuća funkcija na (a,\infty). Tada su sljedeće tvrdnje ekvivalentne:
(i) L(N) = 0
(ii) za svaki n \in \mathbb{N} postoji x_{n} \gt a takav da je N(x_{n}) \lt 1/n.

Dokaz. (i)\Rightarrow(ii). Budući da je L(N) \lt L(1/n) (pri čemu s 1/n označavamo pripadnu konstantnu funkciju), iz svojstva (2) iz Definicije 1 slijedi da postoji c \in \mathbb{R} takav da je N(x) \lt 1/n za x\gt c.

(ii)\Rightarrow(i). Budući da je N padajuća na (a,\infty), vrijedi 0\lt N(x)\lt 1/n za x\gt x_{n}. Iz Teorema 4 i svojstva (1) dobivamo 0 \le L(N) \le 1/n za svaki n \in {\mathbb{N}}. Dakle, L(N) = 0.
\ \blacksquare


Tvrdnje sljedećeg korolara jednostavno slijede iz definicije limesa na klasi BM(\infty), dane u dokazu Teorema 2:

Korolar 6. Neka je f(x) = \lambda + N(x), g(x) = \lambda - N(x). Sljedeće su tvrdnje ekvivalentne:
(i) N(x) \searrow 0,
(ii) f(x)\searrow \lambda,
(iii) g(x) \nearrow \lambda.


Tvrdnje sljedećeg korolara jednostavno slijede iz Teorema 5:

Korolar 7. Pretpostavimo da N'(x)\searrow 0, N''(x)\searrow 0 kada x\to\infty. Tada vrijedi:
(i) Ako je N(x) = N'(x)+ N''(x), tada N(x)\searrow 0 kada x\to\infty.
(ii) Ako je N(x) = C N'(x), pri čemu je C pozitivna konstanta, tada N(x)\searrow 0 kada x\to\infty.

3Konvergentne funkcije i limes

U ovom poglavlju uvodimo pojam limesa na širu klasu funkcija.

Definicija 8. Kažemo da funkcija f konvergira kada x teži \infty ako postoji a \in \mathbb{R} i funkcije m i M, koje su ograničene i monotone na intervalu (a, \infty) i da vrijedi
(1) f je definirana na (a, \infty)
(2) m(x) \le f(x) \le M(x), za x\gt a
(3) L(m) = L(M).
Klasu funkcija f koje konvergiraju kada x teži \infty označavamo s C(\infty).

Napomena 9. Iz Definicije 8 slijedi da je svaka konvergentna funkcija ujedno i ograničena na nekom intervalu (a, \infty).

Teorem 10. Limes L na BM(\infty) može se proširiti na klasu C(\infty), tako da svojstvo (2) iz Definicije 1 ostaje zadovoljeno.

Dokaz. Ako je f iz klase C(\infty), stavimo L(f) = L(m) = L(M). Prvo moramo provjeriti je li ova definicija dobra, odnosno da, ako imamo funkcije m',M', m'', M'' iz klase BM(\infty) takve da m'(x) \le f(x) \le M'(x) za x\gt a' i m''(x) \le f(x) \le M''(x) za x\gt a'', te L(m') = L(M') = L' i L(m'') = L(M'') = L'', onda vrijedi L' = L''. Iz gornjih pretpostavki slijedi da je m'(x) \le M''(x) za x\gt \max \lbrace a',a'' \rbrace. Sada Teorem 4 povlači L' = L(m') \le L(M'') = L''. Slično se pokazuje L'' \le L'. Dakle, L'= L''.

Preostaje provjeriti je li svojstvo (2) iz Definicije 1 zadovoljeno na C(\infty). Neka su sada m',M', m'', M'' iz klase BM(\infty) takve da m'(x) \le f(x) \le M'(x) za x\gt a', m''(x) \le g(x) \le M''(x) za x\gt a'', te L(m') = L(M') = L(f) i L(m'') = L(M'') = L(g). Budući da je L(M') = L(f) \lt L(g) = L(m''), iz svojstva (2) za funkcije iz klase BM(\infty) slijedi da postoji a \in \mathbb{R} takav da je M'(x)\lt m''(x) za x\gt a. Slijedi f(x) \le M'(x) \lt m''(x) \le g(x) za x\gt \max \lbrace a,a',a'' \rbrace.
\ \blacksquare


Vrijede sljedeća poopćenja teorema 3 i 4 na klasu C(\infty):

Teorem 11. Limes L na C(\infty) je jedinstven.

Teorem 12. Ako postoji a \in \mathbb{R} takav da je
f(x)\le g(x) \text{ za } x\gt a,
i ako f i g konvergiraju kada x teži \infty, tada
L(f) \le L(g).


Dokazi teorema 11 i 12 identični su dokazima teorema 3 i 4.

Tvrdnje sljedećih lema jednostavno slijede iz Definicije 8 i odgovarajućih tvrdnji na klasi BM(\infty) (korolari 6 i 7):

Lema 13. a) Pretpostavimo da je |f(x)-\lambda | \lt N(x) za x \gt a. Tada, ako N(x) \searrow 0 kada x \to \infty, onda f(x) \to \lambda kada x\to\infty.

b) Neka je f(x) = \lambda + z(x). Sljedeće su tvrdnje ekvivalentne:
(i) f(x)\to \lambda kada x\to \infty
(ii) z(x)\to 0 kada x\to\infty.

Lema 14. (i) Ako z'(x)\to 0, z''(x)\to 0 kada x\to\infty i z(x)=z'(x)+z''(x), tada z(x)\to 0 kada x\to \infty.

(ii) Ako z(x)\to 0 kada x\to\infty i w(x)=b(x)z(x), pri čemu je funkcija b ograničena na nekom intervalu (a,\infty), tada w(x)\to 0 kada x\to\infty.


U sljedećem teoremu navodimo uobičajena svojstva limesa s obzirom na zbrajanje, množenje i dijeljenje funkcija:

Teorem 15. Pretpostavimo da f(x)\to \alpha i g(x)\to\beta kada x\to\infty. Neka je s = f + g, p = f g i q = 1/g. Tada s(x)\to\alpha + \beta, p(x)\to \alpha \beta kada x\to\infty. Ako je i \beta \neq 0, tada q(x) \to 1/\beta kada x\to\infty.

Dokaz. Slijedi iz lema 13 i 14.
\ \blacksquare

Teorem 16. Pretpostavimo da f konvergira kada x teži \infty i da je L(f) = \lambda. Tada vrijedi:
(i) Ako je \alpha \lt \lambda \lt \beta, tada postoji a \in \mathbb{R} takav da je \alpha \lt f(x) \lt \beta za x \gt a,
(ii) Za svaki \epsilon \gt 0 postoji X=X(\epsilon) \in \mathbb{R} takav da je
|f(x)-\lambda| \lt \epsilon \text{ za } x \gt X.

Dokaz. (i) Kao i do sada, koristimo se istom oznakom za realan broj i za konstantnu funkciju čija je jedina vrijednost taj realan broj, dakle L(\alpha)= \alpha, L(\beta)= \beta. Pretpostavka je da
L(\alpha) \lt L(f) \lt L(\beta).
Koristeći se svojstvom (2) dobivamo
\alpha \lt f(x) \text{ za } x \gt a', f(x) \lt \beta \text{ za } x \gt a'',
pa tvrdnja (i) slijedi za a=\max\lbrace a',a''\rbrace.

(ii) Slijedi iz (i) za \alpha = \lambda - \epsilon, \beta = \lambda + \epsilon.
\ \blacksquare

Teorem 17. Pretpostavimo da za svaki \epsilon \gt 0 postoji X = X(\epsilon) \in \mathbb{R} takav da
(2)
\begin{eqnarray} |f (x)-\lambda | \lt \epsilon \text{ za } x \gt X. \end{eqnarray}
Tada f konvergira kada x teži \infty (u smislu Definicije 8) i L(f) = \lambda.

Dokaz. Moramo provjeriti jesu li zadovoljeni uvjeti Definicije 8. Neka je M(x) = \sup \lbrace f(t)| \ t \ge x \rbrace, m(x) = \inf \lbrace f(t)| \ t \ge x \rbrace. Lagano se pokazuje da su funkcije M i m iz klase BM(\infty) i da je m(x) \le f(x) \le M(x), recimo za x \gt X(1). Preostaje provjeriti je li L(m) = L(M) =\lambda. Pokazat ćemo da je L(M) =\lambda, analogno se pokazuje L(m) = \lambda.

Za proizvoljan \epsilon \gt 0, iz relacije (2) slijedi da je \lambda - \epsilon \lt f(t) \lt \lambda + \epsilon za t \gt X(\epsilon). Dakle, \lambda - \epsilon \le M(x) \le \lambda + \epsilon za x \gt X(\epsilon). Iz Teorema 4 sada slijedi da za svaki \epsilon \gt 0 vrijedi \lambda - \epsilon \le L(M) \le \lambda + \epsilon. Odavde dobivamo L(M) = \lambda.
\ \blacksquare


Teoremi 16 i 17 pokazuju da je definicija limesa funkcije iz ovog rada ekvivalentna uobičajenoj definiciji limesa funkcije.
Bibliografija
[1] B. M. Baishanski, A more intuitive definition of limit, arXiv:0805.3671
[2] S. Kurepa, Matematička analiza 1, Tehnička knjiga, Zagreb, 1989.
[3] S. Kurepa, Matematička analiza 2, Tehnička knjiga, Zagreb, 1990.