06.03.2013 Views

Capitolul 1 Ecuatii diferentiale de ordinul ˆıntâi rezolvabile prin ...

Capitolul 1 Ecuatii diferentiale de ordinul ˆıntâi rezolvabile prin ...

Capitolul 1 Ecuatii diferentiale de ordinul ˆıntâi rezolvabile prin ...

SHOW MORE
SHOW LESS

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

Meto<strong>de</strong> numerice 135<br />

ceea ce <strong>de</strong>monstreazǎ cǎ funct¸iile X ε (t) sunt echicontinue pe Iα.<br />

Cu teorema lui Arzela-Ascoli rezultǎ cǎ existǎ un ¸sir εn → 0, astfel ca<br />

¸sirul {X εn }εn sǎ fie uniform convergent pe intervalul Iα la o funct¸ie continuǎ<br />

X pe intervalul Iα, ¸si aceasta satisface X(t) −X 0 ≤ b, pentru orice t ∈ Iα.<br />

Continuitatea uniformǎ a funct¸iei F pe cilindrul ∆ ¸si convergent¸a uniformǎ<br />

a ¸sirului <strong>de</strong> funct¸ii X εn la funct¸ia X asigurǎ convergent¸a uniformǎ a<br />

¸sirului <strong>de</strong> funct¸ii F(τ, X εn (τ)) la funct¸ia F(τ, X(τ)) pe intervalul Iα.<br />

Trecem la limitǎ în egalitatea:<br />

X εn (t) = X 0 +<br />

t<br />

¸si obt¸inem cǎ funct¸ia X(t) verificǎ<br />

X(t) = X 0 +<br />

t0<br />

t<br />

t0<br />

F(τ, X εn (τ))dτ +<br />

t<br />

t0<br />

θ εn (τ)dτ<br />

F(τ, X(τ))dτ, (∀) t ∈ Iα.<br />

Aceasta <strong>de</strong>monstreazǎ cǎ limita X = X(t) este solut¸ia problemei cu date<br />

init¸iale<br />

˙X = F(t, X), X(t0) = X 0 .<br />

Din teorema <strong>de</strong> unicitate rezultǎ cǎ funct¸ia X(t) coinci<strong>de</strong> cu solut¸ia saturatǎ<br />

X(t; t0, X0) pe intervalul Iα:<br />

X(t) = X(t; t0, X 0 ), (∀) t ∈ Iα.<br />

Se obt¸ine în acest fel cǎ funct¸ia X ε (t) aproximeazǎ solut¸ia<br />

neprelungibilǎ X(t; t0, X 0 ) pe intervalul Iα.<br />

Valorile funct¸iei X ε (t) în punctele ti se obt¸in cu formula <strong>de</strong> recurent¸ǎ:<br />

X ε (ti) = X ε (ti−1) + h · F(ti−1, X ε (ti−1)), i = 1, q<br />

iar în punctele t−i cu formula <strong>de</strong> recurent¸ǎ:<br />

X ε (t−i) = X ε (t−i+1) − h · F(t−i+1, X ε (t−i+1)), i = 1, q<br />

Aceste proceduri <strong>de</strong> trecere <strong>de</strong> la (ti−1, Xε i−1 ) la (ti, Xε i<br />

(t−i+1, Xε −i+1 ) la (t−i, Xε −i ) sunt u¸sor <strong>de</strong> programat.<br />

) ori <strong>de</strong> la

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!