Equazioni differenziali
Equazioni la cui incognita è una funzione. La funzione compare con le sue derivate.
Esempio 1
y'(x) = xy(x)? Incognita y(x) = ∫ g(x) c(x) x dx = x2/2 + cy'(0) Pr problema di Cauchy Condizione iniziale y(0) = 0 y(x) = x2/2 + c1 = 02 / + c → c = 1 y'(x)=f(x) → y(x) = ∫ f(x) dx
Esempio 2
y'(x) + tg(x) = xy(x) = x - y(x). Non si può procedere come prima perché non possiamo trovare una primitiva di g(x). Le frecce rappresentano la direzione delle derivate in ogni punto. Una soluzione y = x - 1. Verifica (x - 1) = x.
Equazioni differenziali
Equazione la cui incognita è una funzione. La funzione compare con le sue derivate.
Es. 1
y'(x)=x. Incognita: y(x)=∫g(x)dx=∫xdx= x2/2 + Cy'(x0) = 0. Problema di Cauchy. Condizione iniziale y(x)= x2/2 + Cy(x0) = y01= 02/2 + C→ C = 1.
Es. 2
y'(x)+tg(y(x))=xy(x)= x-y(x). Non si può procedere come prima perché non possiamo trovare una primitiva di g(x). Le frecce rappresentano la direzione delle derivate in ogni punto. Una soluzione y= x-1. Verifica (x-1)= xy(x) + y(x) = x.
Equazione differenziale lineare
Eq. differenziale lineare (equazione lineare ordinaria di 1°) del primo ordine (grado di derivazione) y(x) + y(x) = x.
Si moltiplica per un fattore integrante, in questo caso ex. exy(x) + exy(x) = x ex.
D(ex y(x)) = ex x: Ora abbiamo, con condizione ideale: derivata della f.ne uguale alla k, possiamo quindi integrare ∫D(exy(x)) dx = ∫ex x dx => exy(x) = ∫ex x dx + c: g(x) = ex(∫ ex x dx + c)e-x.
∫ex x dx: exx|ex-e-x. g(x) = x - 1 + c e-x. a) g(0) = 0 y0(0) = -1|c = 0 => c = 1.
In generale
y + a(x) g(x) = b(x). Eq. differenziale lineare del 1° ordine (a coefficienti non costanti).
Fattore integrante: A(x) e∫a(x). A(x) = una primitiva di a(x). D(eA(x) g(x)) = e∫A(x) f(x).
eA(x)y(x) + ∫a(x) g(x) = e∫A(x) b(x). eA(x)g(x) = ∫A(x) g(x) si integra.
y(x) = e-∫a(x)dx [∫ e∫a(x)dxc dx + C. Formula risolutiva y'(x) + g(x) x = ⇒ a(x) = 1 f(x) = x A(x) = x y(x) = e-∫x |e∫xdx.
L'importante é ricordare: la costante arbitraria basta il risultato é sbagliato.
Esempio
y(x) y'(x) = f(x)=4x2 a(x) = -1 ⇒ A(x) = log x ⇒ eA(x) = |x| g(t) = -∫xg(t) = 0 g(x) = 1/|x| [| | 4x2 dx = x = x0 y(x)= x ∫ 1/x [| 4x x4 dx= = xx 1/x3 (y'= x(sub>-4+c)= = |^(3)(C/x 0 = g(1)= =(-1)3 = x0 diff non la condizione basi che per x0 l'eq integrale che x=1 e dei dei per 'x=1' ) ε ε.
Verifica g(x) = 3 - 1/x1/2 | x 2/x| ety.
ES.{( ex + 1 ) ẏ(x) - ḡ(x) = ( ex + 1 ) exy(b) = 1} ẋ = ( ex + 1 )e A(x) = ∫ dx = ∫ 1A(x) = ∫ 1 dx = ∫∫ dx =e = ( ex + 1 )= ln ( ex + 1 )e A(x) = e-xg(x) =.
In fine g(x) =X + ex -3ex + 1C = 3.
ES.{ẏ(x) + 2x ẏ(x) = 3x + 4 g(0) = 5g(x) espisto a ∞g(x) = e-x2 B(x3 + 4) g(x) =lim e-x2x g -d B(x3 + 4)lim ex2 + 4.
Ta qui la primitiva non si prova a inserire 2 estremi dell'integrale {g(0) = ∫[ 2 (B(x4 + 3 lim g(0) = lim- x 2 e ... (3x + 4)C = 5.
y''(x) + 2y(x) = 3e-2xg(x) = 3xA: x = 2, eA(x) = e-2xeA(x) = x e-2 = A = 2xg(b) = 1g(x) = 1/e2 ∫e2x e-2x3 = x = A.
In fine g(x) = BX + 1/e3x e-2xy = | yf |.
Equazioni differenziali lineari a coefficienti costanti
an y(n)(x) + am y(m-1)(x) + ... + a0 y(x) = f(x) ordine n. y = e3.
ESE.
y''−3y' + 2y = x2.
- Si prende prima la parte omogenea e si risolve.
y(x) = yo(x) + yg(x). La soluzione sarà data dalla somma di: quella omogenea e quella particolare.
y''−2y' + y = 0. y = e2xy2,3 y = 2y.
Polinomio caratteristico: λ2−3λ+2=0, (x−2)(x−1)=0. y' = 0 y = xyo = C1 + C2ex2xλ2=1.
Si prendono le soluzioni e si costruiscono gli esponenziali in questo modo. yo = C1 + C2xy' y − e2x = 0. yo(x) = C1ex + C2x ex.
Se le soluzioni sono distinte si hanno un numero di esponenziali uguale al numero delle soluzioni con ciascuno il rispettivo esponente.
In generale
anzn + an-1zn-1 + ... + a0 = 0.
Fattorizzazione
- (zn - 1) ⇾ mdx ⇾ 2dx ⇾ ... ⇾ z...z ⇾ n fattori
- (z2 + 1)n ⇾ grado complessivo: 2m - auto
- Δko ↘ dx eiβx ⇾ el ecosβx
- ↘ x ex senβx x ex dx
ES. 1 - Parte omogenea
y(x) - 2y(x) - 3y(x) = 0 z2 - 3z - 3 = 0 z1,2 = 2 ± √16/2 (z+1)(z-3) = 0.
g(1)0(x) = c1e3x + c2e-x.
ES. 2
y(x) + 4y(x) + 6y(x) + 10 = 0 z2 + 4z + 6 = 0 (z+2)2 = 0 ⇒ z = -2 ordine 2.
g0(x) = c1e-2x + c2xe-2x.
ES. 3
y + 2y1 + 5y = 0 z2 + 2z + 5 = 0 -2 ± √4 - 20 / 2 = -1 ± 2iyg(x) = yg(x) = e-x(c1 cos(2x) + c2 sen(2x)).
ES. 6.1
g(x) - g'(x) = 0. g(0) = 1. g(0)' = 6. g(x) = c1e3x + c2ex.
y(x) = c1e3x + c2ex. g(x) = 3c1e3x + c2ex.
c1 + c2 = 13. c1 + c2 = 6. c1 = 5/2. c2 = -3/2.
ES. 5
f(x) = 16 g(x) lim (x->4)= g(x) = 3G(x) = 0. Se il limite v. +/- oo di convergenza dell'integro = 0.
ES. 6
f(x) = gg(0) = h0 g(6) = v0. Ora ci occupiamo di trovare la soluzione particolare. Questa è la forma della f.n.e. a queste fa si associa un numero z1+jβ si determina la molteplicità m di z1β nel pol. caratteristico.
yp(x) = eβx [ Pm(x) cos βx + Qm(x) sen βx ] ← Polinomi generici di grado ≤ a quello di P(yk)(g) = gyyh(0) = h0(y)g(0) = v00+i0=0.
y(x) = C1 ⋅ iχ + Ch + ( t eβx xi )2ᴧῸᴿᴇᴿᴵᴬm = 2z2 = 0.
(yR)(g)(x) = x ⋅ Agx(x) = 2x ⋅ A gx(x) = 2A2A = gA = g/₂gk(x) = ½ g x2(y(x)) + 3(y(x)) + 8xz3 + Β2z2 = 0z[(z+3)-0z = 0ord. 1z-3-ᓘᓗ-3ᓛ.
yh(x) = C4 ex + Cx gx e-3xa(x) = C4 x + eβx = A)bx(x) = (A - (Cxf) = A3+b χa(x) => b(x) => h = 0.
Ha multiplicítà δ nel polinomio caratteristico m = zgχ(c) = z (AХ+B) = A0+Z3 .....y(x) = (v aA+B) x6Α+3 (6Aχ+Β) = 8χ18Αx=8χ A=1/2C(A+ḐAχ+⌑∃Σ∃χ3χCA +5(6Ax +Ƀ) 9χ8Α+5Α+Ƀ(J(x)) = 6A.
g(x) = C1 ⋅ iχ + Cx + e-3x− 3/₂ χ3 −⃡1−1/3χ2.
y(x) = e2x ⋅ (P1(x) cos βx + Q1(x) sen βx)θ = β. Bisogna controllare come si comporta questo numero rispetto al polinomio caratteristico (molteplicità).
ym(x)= xm ⋅ eαx (Q1(x) cos βx + Q2(x) sen βx) 3:y'{ y″(x) -2y′(x)-y(x)= x - e3x(x) = C1 + C2e2x.
ES.2− 2−0 ℤ(−2)−0 −0. Sapere per una possa per l'(. La molteplicità a nel pol. caratteristico g{(x)= X ((x+ ){'A1Ax + BxG2(x)= 2A + B2A-2(ax+5)xg(x)=gek(x)+te-3x3.
Molteplicità a nel pol. caratteristico+ em =-2 e{q}= paralel01 g(x) = S A e3x∫excos x dxz=0 =>1y0(x)=C1g(x)=ex∙cos(x)gx(x) = xy∞ cos y∞ mod. g(x) = 0.
Continua...
Eq. diff. a variabili seperabili (non lineare)
y-1(x) = g1(x)z:=0 g(x)=±x7y2(x) g′(x) − g2(x).
- y0(x)−g1(x)0 ≠ +∫ y′(x) (y(x)− g−1(x))−2X∞ ∫y′(x)−g(x)g(y) ∫ (dy/y−g)y−g = x−x0
1/(1−y)A/(g−y) + B/(g) + 1/(−g)⇒ ∫ dy/y + ∫dy/−g = ln|g|−ln|y|=ln|y/g|x = loglog (g(x))y: = log (g(x)/(1−g(x)))x=log g(x)xg(x)/(1−g(x))ex/(1+ex)x=0g(x)=ex = exg(x)g(x)+exy2(x)=exy(x)= ex/(+e).
In generale
- y'(x)=a x : b y(x) (y0 - ʸ₀) / (b(y₀)) x = ∫xx₀ - a(x) dx →
- ∫yy₀ dy / b(y) = ∫xx₀ a(x) dx
- H(g(x)) - H(g₀) = A(x) - A(x₀)
- H(g(x)) = A(x) - A(x₀) + H(g₀)
g(x)=h-1(A(x) - A(x₀) + H(g₀)).
ES.
{ y'(x) = 2 x (g(x))2y(0) = 1y / (g(x)) = 2x x₀ ed. potenziale ∫x y / (g(x))2 dx = ∫x 2x dx →∫x dy y2 x = ∫x 2x dx→1 / g(x) - 1 / x2→1 / y = x2.
(x+1) y'(x) = x (y(x)-1)y(0) = ∞g(x) =g(y) =1 dxg(y(x)-1)g(y)x / x2 +1dy / y-1 = yx / x+1ln (y(x)-1) 0 y(x)½ ln (x2 +1)ln (g(x)-1) =½ ln (1 + x2)y(x) = 1 - √(1+x2)g(x) =g(x) = 1 - √(1+x2)∞ ẟ S.I.
g(x) = 2x (cos (x))xg(y(x)-1) =cod x=0g(cos x2+1)lim g(x)(x→+∞)dy = ∫ x dxtg (g(y)) = [x2]cos (g)1-0tg g(x) = x2g(x) = arctg (x)2 |x| > 1y(x) = arctg (x) + π/2lim g(x-y) gt x→-∞ y→1-x-2-1