Scarica il documento per vederlo tutto.
Scarica il documento per vederlo tutto.
Scarica il documento per vederlo tutto.
Scarica il documento per vederlo tutto.
vuoi
o PayPal
tutte le volte che vuoi
Equazioni Differenziali
y' = f(t)g(y)
f(t) sia definite su un intervallo I ⊆ ℝ
g(y) ···································· J ⊆ ℝ
Soluzioni costanti:
Prop. ϕ(t)=K ∀ t ∈ I è sol. dell'eq. diff. g(K)=0y=K punto di equilibrio o stazionario per l'eq. diff.
Teorema di Peano (esistenza locale delle soluzioni)
(y' = f(t)g(y))f continua in un int I che contiene t₀,g continua in un int J che contiene y₀.
∃ un intervallo aperto I₀ ⊂ I, ∃ ϕ: I₀ -> ℝ derivabile etc(ϕ'(t)=f(t)g(ϕ(t)) soluzione dell'eqϕ(t₀)=y₀ soddisfa la condizione
y' = ³√yy(0)=0f(t)=1 g(y)=³√y è continue in y=0Questo problema di Cauchy ha infinite soluzioni
Oss. ϕ'(t), f(t), g(ϕ(t))Cont. per ip.Cont. su I=> ϕ(t) è continuse su I₀=> ϕ ∈ C¹ (I₀)
Definizione
f : I ⊂ ℝ ⟶ ℝ è lipschitziana con costante di Lipschitz h > 0 se
- ∃ h > 0 t.c. ∀ x1, x2 ∈ I |f(x1) - f(x2)| ≤ h |x1 - x2|
- se x1 ≠ x2 |(f(x1) - f(x2))/(x1 - x2)| ≤ h
- f(x)=|x| non è lipschitziana in un intorno di 0 perché
- limx⟶x0 (f(x) - f(x0))/(x - x0) = +∞
- ∀ h : ∃ B(0) x ∈ B(0) => |(f(x) - f(0))/(x - 0)| > h
- 1) f è lipschitziana in I ⇒ f è continua su I
- 2) f è di classe C1 in I ⇒ f è lipschitziana su I
Teorema di Cauchy/Picard-Lindelof
esistenza e unicità locale della soluzione
- { y' = f(t) g(y) / y(t0) = y0
- f continua in I x Iy c ℝ cont t0
- g lipschitziana in Iy c ℝ cont t0
=> ∃ I0 ⊂ I, ∃ φ ∈ C1 su I0 φ : I0 ⇒ ℝ
- tale che φ'(t) = f(t) g(φ(t)) ∀ t ∈ I0
- φ(t0) = y0
- e φ è l'unica soluzione del prob. di Cauchy su I0
Sol. globali = sol. definite su I1
Pannello di Peano:
1 4 1() = 3√ 1(0) = 0
() = 1
Cont. su
3√
Non è lipschitziana in
alcun intorno di 0
Non vale il teo di esistenza e unicità delle soluzioni
1)
41 = 0 4 = 0 1(t) = 0
l'unica sol. costante def. su
2)
'1 = ( 3√ ) · 1
t ≤ 0
1/3
'1 = 1
=> 1/3∫ d = ∫ dt = t + c
∫ d = ∫ 3/2 4/3
3/2 = 2/3 [t + c]3
= ± [2/3 (t + c)]3
'1(t):
[2/3 (t + c)]3
2√ (2/3 (t + c))3
(0) = 0
dom {t ∈ : t > c}
dom =
dom :
γ(t):
[1(t)
0
t < c
ψ(t):
[2(t)
0
t < 0
Sono soluzioni def su D se sono derivabili e se
' 1(t) = 3√
√(2/33√) t > 0
' 1(t)
[([2/3 t]3/2) ]1/2 : 1/2 1/√ 3/2 t√ 1/6 √/ 1/25 √/ 1/6 1/5
Equazioni differenziali lineari del secondo ordine
(*) y'' + a(t)y' + b(t)y = f(t)
a(t), b(t), f(t) def. su uno stesso intervallo I ⊂ ℝ
eq. complete (o affine) se f(t) ≠ 0 per qualche t ∈ I
eq. omogenee se f(t) = 0 cioè f(t) = 0 ∀ t ∈ I
f(t) termine forzante y'' ::= a(t)y' + b(t)y + f(t)
Esempio: y'' - w1y
y'' + w1y = 0 oscillatore armonico
eq. diff. lineare a(t):=0 b(t):= w12 omogenea
Oscillatore armonico smorzato (con attrito)
y'' - w1y - b1y' (b > 0) => y'' + b1y' + w12y = 0
- q(t) := b
- b(t) := w12
- omogenea
- f(t) := g
- [termino forzante]
Caduta di un grave in un mezzo viscoso
Teorema di esistenza e unicità della soluzioni
(y'' + a(t)y' + b(t)y = f(t))
PC { y(t0) = y0, y'(t0) = y'0}
se a(t), b(t), f(t) continue su un intervallo I contenente t0
=> ∃ unica soluzione di PC, definita su I
Cioè y: I → ℝ derivabile 2 volte su I e tale che
{ y''(t) + a(t)y'(t) + b(t)y(t) = f(t) ∀ t ∈ I
y(t0) = y0, y'(t0) = y'0}
Cerco le φ(t) = ezt che sono soluzioni dell'eq(1)
φ'(t) = zezt φ''(t) = z2 ezt
=> z2 ezt + aezt + bezt = 0 ∀t ∈ ℝ
=> ezt (z2 + az + b) = 0
ezt ≠ 0 ∀t
=> posso dividere per ezt
(z2 + az + b) = 0 eq. algebrica
Quindi dati φ'' + aφ' + bφ = 0 a,b,c ∈ ℝ
φ(t) = ezt è soluz. dell'eq. (=> z è una soluzione dell'eq.
(z2 + az + b = 0)
- Se Δ > 0
z2 + az + b = 0 ha 2 sol. reali z₁, z₂ ∈ ℝ
=> ho 2 soluzioni reali ez₁t e ez₂t
e linearmente indipendenti
=> { C1 ez₁t + C2 ez₂t , C1, C2 ∈ ℝ }
integrare generale
Esempio. φ' + 3φ' + 2φ = 0
z2 + 3z + 2 = 0 = (z+2)(z+1)=0
=> {C1 e-2t + C2 e-t , C1, C2 ∈ ℝ }
- Se Δ = 0
z2 + az + b = 0 => una sol. reale doppia (z - z₁)2 = 0
Si è scoperto che anche φ(t)= t ez₁t è soluzione
infatti φ'(t)=z1 ez1t + t z2 ez1t e= ez1t (1+t z₁)
φ''(t)=z12ez1t + z1ez1t + t z12ez1t e= t ez1t (t z2 + t z1)
=> e-z1t (z2 t + z1 > z2 ezzt + b = e = 0 ∀t
ezt (z1 + t2 + 1 + t ( z₁ + b )) =e t (zzz,t + 1) (z -z2 t z1+ b))
=> z₁, z₂ = ±√Δ
=> a,z = z%^2=0%^2=0
αlibiamo soluzioni φ1(t) = ez1t φ2(t) = tez1t
(c a/c - 3d 1/c)d t + (2 c 1 2d 3d) = c 1 t 1
- 5c t (2 c - d) = 4 t 1
{c a e^2t + c2e-4t - 1/5
fh(t) smt Pm(t)e d t cosp t
pm(t) = 0 Qm(t)^n grado 0 β= n , 2 = v
Ψ(t) = c e2t cosp t + d e2 t sinp t
Ψ'(t) = - c sin t d cos t
Ψ''(t) c cos t d sin t
Ψ(t) = c2 e-4t + 1/5
∴Ψ(t) 3/10 cos t - sin t
y'' - 3y 2y = fh(t) smt
cos t C-3d 2c) e-4t/3
∫ pdt sin t sint
(c-3d=0 ⇒ c 3d:d
Ψ(t) = 3/10 cost 440 sint
c2e2t + c1e3t
3/10 cost 41 sint
y'' - 34 + 24 - 32t + sint
se lo ho la same lasciar
delle forzono
y4(t) 3 t2 t2 10 t sint