Metodi 1

From Qualcuno?
Jump to navigation Jump to search

Come affrontare un'equazione differenziale a variabili separabili?

Partiamo da un'equazione come $$ \frac{\partial u}{\partial t} = \frac{\partial^2 u}{\partial x^2}\text{.}$$

Ottenere le soluzioni generiche

Prima di tutto, separare le variabili e assumere che le soluzioni separate siano somme di esponenziali complesse. Assumiamo che <math>u(x,t) = X(x)T(t)</math>.

Da ciò, deriviamo. Ottenendo: $XT^\prime = X^{\prime\prime}T$, e quindi $$\frac{T^\prime}{T} = \frac{X^{\prime\prime}}{X} = -\lambda^2 \text{;}\qquad \begin{cases}T^\prime + \lambda^2 T = 0 \\ X^{\prime\prime} +\lambda^2 X =0 \end{cases}$$

Assumiamo anche che <math>X(x) = b e^{\beta x}</math>, e che <math>T(t) = a e^{\alpha t}</math>.

$$\begin{cases}a\alpha e^{\alpha t} + \lambda^2 a e^{\alpha t} = (\alpha + \lambda^2) a e^{\alpha t} = 0 \\ b \beta^2 e^{\beta x} + \lambda^2 b e^{\beta x} = (\beta^2+\lambda^2) b e^{\beta x} = 0 \end{cases} $$

Una volta impostato il sistema con entrambe le equazioni per $X$ e $T$, cerchiamo di ottenere $\alpha$ e $\beta$ in funzione della costante $\lambda$ che collega le due equazioni.

$$\begin{cases} \alpha + \lambda^2 = 0 \\ \beta^2+\lambda^2 = 0 \end{cases}$$ $$\begin{cases} \alpha = - \lambda^2 \\ \beta = \sqrt{\lambda^2} \implies \beta_1 = \lambda \text{;} \beta_2 = i^2 \lambda \text{.} \end{cases}$$

Abbiamo quindi ottenuto due equazioni, $$\begin{cases}T(t) = a e^{-\lambda^2 t}\\ X(x) = b_1 e^{\lambda x} + b_2 e^{-\lambda x}\end{cases}\tag{generica}$$

Condizioni al contorno iniziali

Adesso possiamo applicare le condizioni al contorno. Nell'esame che stiamo seguendo, queste erano: $$\begin{cases}u(-\frac{\pi}{2},t) = 0 \\ u_x(\frac{\pi}{2}, t) = 0\end{cases}$$

Sostituiamo quindi quanto ottenuto in precedenza nel sistema delle condizioni al contorno: $$\begin{cases} b_1 e^{-\frac{\pi}{2} \lambda} + b_2 e^{\frac{\pi}{2}\lambda } = 0 \\ b_1 \lambda e^{\frac{\pi}{2}\lambda} - b_2 \lambda e^{-\frac{\pi}{2}\lambda } = 0 \end{cases}$$ Notiamo che $\lambda = 0$ ci conduce alla soluzione banale... $$\begin{cases} \frac{b_1}{b_2} e^{-\frac{\pi}{2} \lambda} = - e^{\frac{\pi}{2}\lambda } \\ \frac{b_1}{b_2} e^{\frac{\pi}{2}\lambda} = e^{-\frac{\pi}{2}\lambda } \end{cases}$$ Sostituiamo $-1 = e^{(1+2n)i\pi}$, $n \in \mathbb{Z}$: $$\begin{cases} \frac{b_1}{b_2} e^{-\frac{\pi}{2} \lambda} = e^{(1+2n)i\pi} e^{\frac{\pi}{2}\lambda } \\ \frac{b_1}{b_2} e^{\frac{\pi}{2}\lambda} = e^{-\frac{\pi}{2}\lambda } \end{cases}$$ Divido la prima equazione per $e^{-\frac{\pi}{2} \lambda}$ e la seconda per $e^{\frac{\pi}{2}\lambda}$: $$\begin{cases} \frac{b_1}{b_2}= e^{(1+2n)i\pi} e^{\pi\lambda } \\ \frac{b_1}{b_2}= e^{-\pi\lambda } \end{cases}$$ Sostituisco la seconda equazione nella prima, e divido l'equazione così ottenuta per $e^{-\pi\lambda}$, inoltre, scrivo $1$ come $e^0$: $$\begin{cases} e^0 = e^{(1+2n)i\pi} e^{2\pi\lambda } \\ \frac{b_1}{b_2}= e^{-\pi\lambda } \end{cases}$$ Prendo il logaritmo della prima equazione: $$\begin{cases} 0 = (i + 2ni + 2\lambda)\pi \\ \frac{b_1}{b_2}= e^{-\pi\lambda } \end{cases}$$ Da cui, $$\begin{cases} \lambda = i\frac{1+2n}{2} \\ b_1= b_2 e^{-i\pi\frac{1+2n}{2}} \end{cases}$$

Le soluzioni, date le condizioni iniziali, sono quindi, ricordando $ n\in \mathbb Z$

$$\left\{\begin{array}{ll} T_n(t) &= a e^{(\frac{1+2n}{2})^2 t}\\ X_n(x) &= b (e^{-i\pi\frac{1+2n}{2}} e^{i\frac{1+2n}{2} x} + e^{-i\frac{1+2n}{2} x})\\ &= b (e^{i\frac{1+2n}{2}(x-\pi)} + e^{-i\frac{1+2n}{2} x}) \end{array}\right.$$

Parità e condizioni al contorno generiche

La consegna chiede poi di verificare che perché le $X_n$ abbiano una parità definita, la condizione $AD=-BC$ sia verificata nelle seguenti condizioni al contorno:

$$\begin{cases} A u(-\frac{\pi}{2},t) + B u_x(-\frac{\pi}{2}, t) = 0\\ C u(\frac{\pi}{2},t) + D u_x(\frac{\pi}{2}, t) = 0 \end{cases}$$

Per studiare la parità è conveniente scrivere $X$ e $X^\prime$ in termini di seni e coseni, usando la relazione $e^{i\theta} =\cos{\theta} + i\sin{\theta}$, ovvero $e^z = \cos{(iz)} - i\sin{(iz)}$:

$$\begin{array}{rcl} X(x) &=& b_1 e^{\lambda x} + b_2 e^{-\lambda x}\\ &=& b_1 (\cos (i \lambda x) - i \sin (i \lambda x) ) + b_2 (\cos (i\lambda x) + i \sin (i\lambda x) ) \\ &=& (b_1 + b_2) \cos (i\lambda x) - i(b_1 - b_2) \sin (i\lambda x) \\ \\ X^\prime (x) &=& \lambda (b_1 e^{\lambda x} - b_2 e^{-\lambda x}) \\ &=& \lambda \left[ b_1 \left(\cos (i \lambda x) - i \sin(i \lambda x)\right) - b_2 \left(\cos (i \lambda x) + i \sin (i \lambda x) \right)\right]\\ &=& \lambda \left[ (b_1-b_2)\cos (i \lambda x) - i(b_1+b_2) \sin (i \lambda x) \right] \end{array}$$

Da ciò sembra evidente che le $X$ sono una somma di oggetti pari (i coseni) e oggetti dispari (i seni), per imporre una parità definita dovremmo annullare una delle due parti. Per $X$ pari, dovremmo avere $b_1=b_2$, e per $X$ dispari si avrà $b_1 = - b_2$: $$\begin{array}{rcl} X_p(x) &=& 2 b_1 \cos (i \lambda x) \\ X_p^\prime(x) &=& - 2 i \lambda b_1 \sin (i \lambda x) \\ X_d(x) &=& - 2i b_1 \sin (i \lambda x)\\ X_d^\prime(x) &=& 2 \lambda b_1 \cos (i \lambda x) \end{array}$$

A ognuna di queste vanno applicate le nuove condizioni al contorno: $$\begin{cases} AX_p(-\frac{\pi}{2}) + BX_p^\prime(-\frac{\pi}{2}) = 0 \\ CX_p(\frac{\pi}{2}) + DX_p^\prime(\frac{\pi}{2}) = 0 \\ \end{cases}$$ $$\begin{cases} AX_d(-\frac{\pi}{2}) + BX_d^\prime(-\frac{\pi}{2}) = 0 \\ CX_d(\frac{\pi}{2}) + DX_d^\prime(\frac{\pi}{2}) = 0 \\ \end{cases}$$

Common pitfalls

$$\ln (a e^b) = \ln (e^{\ln{a}} e^b) = \ln (e^{\ln{(a)} + b}) = \ln (a) + b \neq b \ln (a) = \ln (a^b) $$