Disambiguazione – Se stai cercando il Teorema di Picard nell'ambito dell'analisi complessa, vedi Teorema di Picard. In matematica, il teorema di esistenza e unicità per un problema di Cauchy, detto anche teorema di Picard-Lindelöf, teorema di esistenza di Picard o teorema di Cauchy-Lipschitz, stabilisce le condizioni di esistenza e unicità della soluzione di un'equazione differenziale ordinaria.
Il teorema dice che dato il problema ai valori iniziali:
![{\displaystyle y'(t)=f(t,y(t)),\qquad y(t_{0})=y_{0},\qquad t\in [t_{0}-\varepsilon ,t_{0}+\varepsilon ],}](https://wikimedia.org/api/rest_v1/media/math/render/svg/d85ecb44a9bd33aeb471c0818c42b73fc769349e)
se
è una funzione lipschitziana in
e continua in
allora per qualche
esiste un'unica soluzione
al problema ai valori iniziali sull'intervallo
Il teorema
Sia
una funzione definita in un intorno del punto
della forma:
![{\displaystyle I\times J=\{(x,y)\in \mathbb {R} \times \mathbb {R} ^{n}:|x-x_{0}|\leq a,\|y-y_{0}\|\leq b\},}](https://wikimedia.org/api/rest_v1/media/math/render/svg/a4ca4d259262558a74de2cc1b91fa0c3bc23906e)
con
,
reali positivi, e si ponga che
è almeno di classe
in tale intorno. Si supponga inoltre che
sia lipschitziana rispetto alla variabile
(uniformemente nella variabile
, cioè la costante di Lipschitz non dipende da
):
![{\displaystyle \|f(x,y_{1})-f(x,y_{2})\|\leq L\cdot \|y_{1}-y_{2}\|\quad \forall x\in I\quad \forall y_{1},y_{2}\in J,}](https://wikimedia.org/api/rest_v1/media/math/render/svg/d72ee0892a107eddbdf09f97145e8d68b3ece826)
con
costante di Lipschitz. Allora il problema di Cauchy:
![{\displaystyle \Theta =\left\{{\begin{array}{ll}y'&=f(x,y)\\y(x_{0})&=y_{0}\end{array}}\right.}](https://wikimedia.org/api/rest_v1/media/math/render/svg/5164e0f17d8cf066b271e13f327d91b15f229eb1)
possiede una soluzione unica.[1]
Sotto l'ipotesi di continuità della funzione è possibile dimostrare l'equivalenza tra il problema di Cauchy e la seguente equazione integrale, detta equazione di Volterra:
![{\displaystyle y(x)=y_{0}+\int _{x_{0}}^{x}f(t,y(t))\mathrm {d} t,\qquad \forall x\in I_{\delta },}](https://wikimedia.org/api/rest_v1/media/math/render/svg/bc08c7986f95b1a774dc0f9f4514b63b9f975e36)
dove
è un intorno di
, con
valore opportuno. L'esistenza di una funzione
che soddisfa al sistema
si verifica se e solo se tale equazione ammette soluzione.
Dimostrazioni
Nel seguito sono elencate due diverse dimostrazioni del teorema. La prima sfrutta concetti basilari di analisi funzionale, mentre la seconda utilizza argomenti di analisi reale e ha il pregio di mostrare come costruire operativamente una soluzione attraverso approssimazioni successive, e di dare una stima generalmente più accurata dell'ampiezza
dell'intervallo di definizione della soluzione.
Prima dimostrazione
Sia
con
. Si noti che
per il teorema di Weierstrass (poiché
è compatto). Nel caso in cui
, ovvero qualora
sia identicamente nulla, il sistema ammette come unica soluzione la funzione costante
, quindi si può supporre
.
Sia
Si può considerare lo spazio metrico
delle funzioni
continue con la norma dell'estremo superiore, e una palla al suo interno, definita da:
![{\displaystyle B=\{g\in X:\|g-y_{0}\|_{C^{0}}\leq b\}.}](https://wikimedia.org/api/rest_v1/media/math/render/svg/282ed8451609a9cb7ff5047038eb67af8d785ed5)
Essendo lo spazio
completo, e
chiuso, allora anche quest'ultimo risulta essere uno spazio completo rispetto alla metrica indotta.
Si procede quindi definendo l'operatore
, detto "operatore di Volterra", tale che
, dove:
![{\displaystyle {\widehat {y}}=y_{0}+\int _{x_{0}}^{x}f(t,y(t))\mathrm {d} t.}](https://wikimedia.org/api/rest_v1/media/math/render/svg/a6904b44d561032d560b75fd01936a6fa3dc2350)
Si nota innanzitutto che
è ben definito, ossia che
si ha
. Infatti:
![{\displaystyle |{\widehat {y}}-y_{0}|=\left|\int _{x_{0}}^{x}f(t,y(t))\mathrm {d} t\right|\leq \int _{x_{0}}^{x}|f(t,y(t))|\mathrm {d} t,}](https://wikimedia.org/api/rest_v1/media/math/render/svg/e9fbf3e09de18eed1dabc37ca274b2062396e53b)
per ogni
. Ma per ipotesi
, da cui si deduce che:
![{\displaystyle |{\widehat {y}}-y_{0}|\leq \int _{x_{0}}^{x}|f(t,y(t))|\mathrm {d} t\leq M|x-x_{0}|\leq M\delta \leq b.}](https://wikimedia.org/api/rest_v1/media/math/render/svg/51ecd26e365d8ab7454ae3428efa19351d462c8d)
Una volta assicurata la buona definizione di
è sufficiente dimostrare che questa è una contrazione su
per completare il teorema. Il teorema delle contrazioni infatti ci assicura l'esistenza di un unico punto fisso di
in
, quindi nel nostro caso di una funzione
tale che
, cioè
![{\displaystyle y(x)=y_{0}+\int _{x_{0}}^{x}f(t,y(t))\mathrm {d} t,}](https://wikimedia.org/api/rest_v1/media/math/render/svg/0f6f261c56d33dfcd18f951685bef9ef998d455b)
definita sull'intervallo
, e risolvente dunque il sistema
. Tenendo conto delle ipotesi su
(in particolare la lipschitzianità) si può scrivere:
![{\displaystyle {\begin{aligned}|F(y_{1})-F(y_{2})|&=\left|\int _{x_{0}}^{x}{\big (}f(t,y_{1}(t))-f(t,y_{2}(t)){\big )}\mathrm {d} t\right|\leq \int _{x_{0}}^{x}|f(t,y_{1}(t))-f(t,y_{2}(t))|\mathrm {d} t\\&\leq \int _{x_{0}}^{x}L|y_{1}(t)-y_{2}(t)|\mathrm {d} t\leq L\delta \|y_{1}-y_{2}\|_{C^{0}}.\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/4e8cfcb4b2663491190890c066d77591fb399f81)
e prendendo l'estremo superiore al variare di
si ottiene:
![{\displaystyle \|F(y_{1})-F(y_{2})\|_{C^{0}}\leq L\delta \|y_{1}-y_{2}\|_{C^{0}},}](https://wikimedia.org/api/rest_v1/media/math/render/svg/d38a0cf82229c624740c936d3f6d52d659386566)
e poiché
,
è una contrazione.
Seconda dimostrazione (Picard-Lindelöf)
Nel corso della seguente dimostrazione si giunge ad una stima generalmente più accurata del numero reale
. Inizialmente, si ponga
. Il passo successivo consiste nel definire per ricorrenza una successione di funzioni
come:
![{\displaystyle \left\{{\begin{array}{ll}y_{0}(x)=y_{0}\\y_{k+1}(x)=y_{0}+\int _{x_{0}}^{x}f(t,y_{k}(t))\mathrm {d} t.\end{array}}\right.}](https://wikimedia.org/api/rest_v1/media/math/render/svg/6ee6e29d54ed54ee189f74d32d12e28ceef0796a)
È necessaria quindi una verifica della buona definizione della successione, più precisamente bisogna mostrare (ad esempio tramite induzione) che
; il passo base è immediato per come è stato definito
, mentre per il passo induttivo si supponga
, da cui banalmente
. Per le ipotesi preliminarmente fatte su
si può quindi maggiorare il valore assoluto di
con
. È dunque di immediata verifica che:
![{\displaystyle |y_{k+1}(x)-y_{0}|=\left|\int _{x_{0}}^{x}f(t,y_{k}(t))\mathrm {d} t\right|\leq \left|\int _{x_{0}}^{x}|f(t,y_{k}(t))|\mathrm {d} t\right|\leq M|x-x_{0}|\leq M\delta \leq b.}](https://wikimedia.org/api/rest_v1/media/math/render/svg/a1781c100e924ac5a4ef958b555fc01e82034c9d)
Si procede nella dimostrazione stimando ricorsivamente la distanza tra due termini consecutivi della successione puntualmente in
con un metodo analogo a quello induttivo appena usato.
Inizialmente si ha:
![{\displaystyle |y_{1}(x)-y_{0}|=\left|\int _{x_{0}}^{x}f(t,y_{0}(t))\mathrm {d} t\right|\leq M\left|x-x_{0}\right|,}](https://wikimedia.org/api/rest_v1/media/math/render/svg/cb67fbdf81ea0e55363f1a3d1b303e8c1c6e7c1a)
mentre per i passi seguenti bisogna usare anche l'ipotesi di lipschitzianità di cui gode
:
![{\displaystyle {\begin{aligned}|y_{2}(x)-y_{1}(x)|&=\left|\int _{x_{0}}^{x}[f(t,y_{1}(t))-f(t,y_{0})]\mathrm {d} t\right|\leq \left|\int _{x_{0}}^{x}|f(t,y_{1}(t))-f(t,y_{0})|\mathrm {d} t\right|\\&\leq \left|\int _{x_{0}}^{x}L|y_{1}(t)-y_{0}|\mathrm {d} t\right|\leq ML\left|\int _{x_{0}}^{x}|t-x_{0}|\mathrm {d} t\right|={\frac {ML}{2}}|x-x_{0}|^{2}.\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/f8dd5d95fdc199d7bf021b7b772c6c40442cbfed)
Per avere una migliore comprensione della formula generale per la stima che verrà data tra poco è consigliabile sviluppare almeno un altro passo dell'induzione:
![{\displaystyle {\begin{aligned}|y_{3}(x)-y_{2}(x)|&=\left|\int _{x_{0}}^{x}[f(t,y_{2}(t))-f(t,y_{1}(t))]\mathrm {d} t\right|\leq \left|\int _{x_{0}}^{x}|f(t,y_{2}(t))-f(t,y_{1}(t))|\mathrm {d} t\right|\\&\leq \left|\int _{x_{0}}^{x}L|y_{2}(t)-y_{1}(t)|\mathrm {d} t\right|\leq {\frac {ML^{2}}{2}}\left|\int _{x_{0}}^{x}|t-x_{0}|^{2}\mathrm {d} t\right|={\frac {ML^{2}}{3!}}|x-x_{0}|^{3}.\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/585c5dcd30ccbcee33dbc042274a6e9d2a5acfed)
Risulta a questo punto chiara la seguente stima generale, alla quale si può giungere tramite processo induttivo:
![{\displaystyle |y_{k+1}(x)-y_{k}(x)|\leq {\frac {ML^{k}}{(k+1)!}}|x-x_{0}|^{k+1}\quad \forall x\in I_{\delta },}](https://wikimedia.org/api/rest_v1/media/math/render/svg/2d9b43e7278c1d307db1bb047219784a12ec1788)
da cui si può dedurre la convergenza uniforme di questa successione di funzioni nell'intervallo
, dato che maggiorando ulteriormente con:
![{\displaystyle {\frac {M}{L}}{\frac {L^{k+1}\delta ^{k+1}}{(k+1)!}}.}](https://wikimedia.org/api/rest_v1/media/math/render/svg/bcb25dc7ad240bf175ceca8b3c7ef7a2eed652d7)
si ottiene chiaramente la ridotta della serie esponenziale numerica:
![{\displaystyle \sum _{k=1}^{\infty }{\frac {M}{L}}{\frac {L^{k+1}\delta ^{k+1}}{(k+1)!}}={\frac {M}{L}}(e^{L\delta }-1).}](https://wikimedia.org/api/rest_v1/media/math/render/svg/acdf01d8f953a04744724e1b8c0ae20e3bf6c10b)
Passando al limite per
e sfruttando nuovamente la lipschitzianità di
rispetto a
, si ottiene la convergenza totale e quindi uniforme della serie telescopica
(maggiorata dalla serie
convergente) alla funzione
, mentre per quanto riguarda il secondo membro della successione definita all'inizio
, la sua funzione integranda converge a
.
Si può a questo punto utilizzare il teorema di passaggio al limite sotto il segno di integrale per ottenere:
![{\displaystyle y(x)=y_{0}+\int _{x_{0}}^{x}f(t,y(t))\mathrm {d} t.}](https://wikimedia.org/api/rest_v1/media/math/render/svg/3d393420e2469a3803c1c340faff0c6b0398aada)
Ma questa è la formulazione integrale (ed equivalente) del problema di Cauchy, quindi per concludere la dimostrazione non resta altro che mostrare l'unicità di tale soluzione. Il modo migliore è procedere per assurdo: si supponga che esista un'altra funzione
(soluzione del PdC) definita in un nuovo intorno
(la notazione rimane coerente con quanto esposto in precedenza) della condizione iniziale
(quindi con lo stesso centro) e tale che esiste
per cui
. Definito
si consideri la relazione (valida per ipotesi di assurdo):
![{\displaystyle |g(x)-y_{0}|=\left|\int _{x_{0}}^{x}f(t,g(t))\mathrm {d} t\right|\leq M|x-x_{0}|\quad \forall x\in I_{\tilde {\delta }}.}](https://wikimedia.org/api/rest_v1/media/math/render/svg/f4232808ed0e4b42003c62196a4b23bad3479077)
Con un procedimento completamente analogo al precedente si giunge però alla stima:
![{\displaystyle |g(x)-y_{k}(x)|\leq {\frac {ML^{k}}{(k+1)!}}|x-x_{0}|^{k+1}\quad \forall x\in I_{\tilde {\delta }}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/275dbb2a345fa439639be363ea55ddce942d02de)
Dato che il secondo membro della disuguaglianza tende a 0 al tendere di
all'infinito, si può dedurre che:
![{\displaystyle g(x)=y(x)=\lim _{k\to +\infty }y_{k}\quad \forall x\in I_{\tilde {\delta }},}](https://wikimedia.org/api/rest_v1/media/math/render/svg/8d19b54eaa4ebc7980d43a6f85c036d628989e8d)
e ciò contraddice l'ipotesi se
, mentre se
non contraddice l'ipotesi ma dimostra che la funzione
quando
che è il nostro intervallo di partenza.
Generalizzazioni
Il teorema è un valido strumento nello studio delle equazioni differenziali, ma a priori garantisce unicamente l'esistenza della soluzione localmente, ossia in un intorno delle condizioni iniziali. Non è assicurata invece l'esistenza di un'unica funzione risolvente
in un intervallo arbitrario (eventualmente tutto
), sotto ipotesi più strette (ad esempio la sublinearità rispetto a
di
) rispetto a quelle richieste per la versione locale. Se
soddisfa queste ulteriori richieste si può dimostrare inoltre che la soluzione ammette un prolungamento massimale sul suo intervallo di definizione.
Un altro enunciato, il teorema di esistenza di Peano, mostra invece soltanto l'esistenza della soluzione (non l'unicità), ma considera una funzione che è solamente una funzione continua, e non lipschitziana. Ad esempio, il secondo membro dell'equazione
con la condizione iniziale
è continuo, ma non secondo Lipschitz. Difatti, l'equazione ha tre soluzioni, di cui la prima è
e le altre due sono:
![{\displaystyle y(t)=\pm {\big (}{\tfrac {2}{3}}t{\big )}^{3/2}.}](https://wikimedia.org/api/rest_v1/media/math/render/svg/bab86cc92777c27bc5752345d0024e2af4e8cc53)
Più in generale, il teorema di esistenza di Carathéodory dimostra l'esistenza per condizioni più deboli per
. Si nota che nonostante tali condizioni siano soltanto sufficienti ci sono risultati, come quello di Okamura, che forniscono condizioni necessarie e sufficienti affinché il problema ai valori iniziali abbia soluzione unica.[2]
Esempi
- Sia dato il problema di Cauchy:
![{\displaystyle \left\{{\begin{array}{ll}y'=\lambda y\\y(0)=1\end{array}}\right.}](https://wikimedia.org/api/rest_v1/media/math/render/svg/64d308247e09890ea032bdd7aa5dd12cb6bf3a80)
- La funzione
soddisfa tutte le ipotesi, quindi localmente la soluzione è unica (in realtà si potrebbe osservare che poiché
per una certa costante reale
la soluzione è globalmente unica al variare di
). La soluzione è quindi (tenendo conto della condizione iniziale
) la funzione ![{\displaystyle y=e^{\lambda x}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/2220f8fdc9ea3f9f8223d28182ed379396a4c556)
- Un tipico esempio di un problema che non rispetta le ipotesi è:
![{\displaystyle \left\{{\begin{array}{ll}y'=3y^{\frac {2}{3}}\\y(0)=0\end{array}}\right.}](https://wikimedia.org/api/rest_v1/media/math/render/svg/35bdfd11d83083433e8bb62e5b59a2b5753cfa11)
- La funzione
non è localmente lipschitziana rispetto a
in nessun intorno dell'origine e infatti non si ha un'unica soluzione con questa condizione iniziale (anzi, se ne possono trovare infinite: è il fenomeno del pennello di Peano), quali ad esempio
oppure
.
![{\displaystyle y''+\omega ^{2}y=0,\quad \omega \in \mathbb {R} }](https://wikimedia.org/api/rest_v1/media/math/render/svg/426142081024bfee201ced522e958f04c98c0315)
- riconducibile mediante sostituzione al sistema:
![{\displaystyle \left\{{\begin{array}{ll}y'=z\\z'=-\omega ^{2}y\end{array}}\right.}](https://wikimedia.org/api/rest_v1/media/math/render/svg/a47b9e7f247be297b67630b8289fdef521e7d18f)
- Aggiungendo le condizioni iniziali (la scelta
è arbitraria)
e
si ottiene come unica soluzione:
![{\displaystyle y(x)=y_{0}\cos(\omega x)+{\frac {z_{0}}{\omega }}\sin(\omega x)}](https://wikimedia.org/api/rest_v1/media/math/render/svg/06ac076fad3f9bfce3c29f27bd7b60de65ba4637)
Note
- ^ Mathworld - Picard's Existence Theorem, su mathworld.wolfram.com. URL consultato il 07-01-2013.
- ^ Ravi P. Agarwal e V. Lakshmikantham, Uniqueness and Nonuniqueness Criteria for Ordinary Differential Equations, World Scientific, 1993, ISBN 978-981-02-1357-2., page 159
Bibliografia
- (EN) Vladimir Igorevich Arnold (1988): Geometrical Methods in the Theory of Ordinary Differential Equations, 2nd ed., Springer, ISBN 0-387-96649-8
- (EN) Vladimir Igorevich Arnold (1992): Ordinary Differential Equations, Springer, ISBN 3-540-54813-0
- (FR) G. Peano, Démonstration de l'intégrabilité des équations différentielles ordinaires Math. Ann. , 37 (1890) pp. 182–228
- (EN) I.G. Petrovskii, Ordinary differential equations , Prentice-Hall (1966) (Translated from Russian)
- (EN) P. Hartman, "Ordinary differential equations" , Birkhäuser (1982)
Voci correlate
Collegamenti esterni
- (EN) Encyclopedia of Mathematics, Cauchy-Lipschitz theorem, in Encyclopaedia of Mathematics, Springer e European Mathematical Society, 2002.
- (EN) Fixed Points and the Picard Algorithm, su krellinst.org. URL consultato il 6 gennaio 2013 (archiviato dall'url originale il 16 giugno 2010).
- (EN) Picard Iteration, su math.fullerton.edu. URL consultato il 6 gennaio 2013 (archiviato dall'url originale il 26 settembre 2013).
- (EN) Proof of the Picard–Lindelöf theorem (PDF), su math.byu.edu. URL consultato il 6 gennaio 2013 (archiviato dall'url originale il 26 settembre 2020).
Portale Matematica: accedi alle voci di Wikipedia che trattano di matematica