Metodo del nucleo risolvente per equazione differenziale

Prima di postare leggi le regole del Forum. Puoi anche leggere le ultime discussioni.
Questa sezione è un contenitore temporaneo per i topic speciali a pagamento (One-Shot)

Se hai acquistato uno o più Topic speciali, puoi pubblicarli qui cliccando su "Apri un Topic".

Il registro completo dei Topic risolti e in corso è disponibile sul proprio profilo.

Metodo del nucleo risolvente per equazione differenziale #78700

avt
ferdcip
Punto
Avrei bisogno di aiuto per risolvere un problema di Cauchy con il metodo del nucleo risolvente:

y''+4y = sin(2x) ; y(0) = 0 ; y'(0) = 1
 
 

Metodo del nucleo risolvente per equazione differenziale #78715

avt
Galois
Amministratore
Ciao ferdcip emt

Abbiamo il seguente PdC

y''+4y = sin(2x) ; y(0) = 0 ; y'(0) = 1

che dobbiamo risolvere applicando il metodo del nucleo risolvente. Iniziamo dal risolvere l'equazione differenziale del secondo ordine omogenea a coefficienti costanti

y''+4y = 0

Il polinomio caratteristico ad essa associato è

P(λ) = λ^2+4

Troviamone gli zeri che si ottengono risolvendo l'equazione di secondo grado

λ^2+4 = 0

che ha come soluzioni

λ_1 = 2i, λ_2 = -2i

Ovvero abbiamo due soluzioni complesse della forma

λ_(1,2) = α±iβ

con α = 0 e β = 2

Di conseguenza l'integrale generale dell'omogenea sarà

y_O(x) = c_1e^(α x)cos(β x)+c_2e^(α x)sin(β x), con α = 0 e β = 2

ossia

y_O(x) = c_1cos(2x)+c_2sin(2x)

Passiamo ora a trovare una soluzione particolare, che indicherò con bary(x) della non omogenea.

Personalmente, visto che il termine noto f(x) = sin(2x) della nostra equazione differenziale di partenza è di tipo particolare, ritengo che sarebbe stato molto più semplice procedere con il metodo di somiglianza. Visto però che il testo dell'esercizio richiede espressamente di utilizzare il metodo del nucleo risolvente procediamo con esso.

Vediamo innanzitutto come si applica tale metodo

Metodo del nucleo risolvente per edo del secondo ordine non omogenee a coefficienti costanti

Data una generica equazione differenziale non omogenea del secondo ordine a coefficienti costanti avente come termine noto f(x) sappiamo che l'integrale generale dell'omogenea ad essa associata sarà della forma

y_O(x) = c_1 y_1(x)+c_2 y_2 (x)

Il metodo del nucleo risolvente ci dice che una soluzione particolare della non omogenea è data da

bary(t)∫_(x_0)^(x)(K(x,p)·f(p))dp

Dove:

x_0 è un generico punto (che scegliamo noi) appartenente al dominio della nostra equazione differenziale. Ovviamente sceglieremo il più semplice possibile.

f(p) è il termine noto dell'equazione differenziale dove, semplicemente, si sostituisce la variabile x con la variabile p

K(x,p) è quello che si dice nucleo risolvente ed è dato da:

K(x,p) = (det[ y_1(p) y_2(p) ; y_1(x) y_2(x) ])/(det[ y_1(p) y_2(p) ; y_1'(p) y_2'(p) ])

ossia è il rapporto tra due determinanti.

------------

Questa la teoria. Applichiamo quanto detto alla nostra equazione differenziale

y''+4y = sin(2x)

il cui termine noto è

f(x) = sin(2x) da cui f(p) = sin(2p)

Abbiamo già trovato l'integrale generale dell'omogenea associata

y_O(x) = c_1cos(2x)+c_2sin(2x)

ossia si ha

y_1(x) = cos(2x) → y_1(p) = cos(2p) → y_1'(p) = -2sin(2p)

y_2(x) = sin(2x) → y_2(p) = sin(2p) → y_2'(p) = 2cos(2p)

Pertanto il nucleo risolvente è dato da

K(x,p) = (det[ y_1(p) y_2(p) ; y_1(x) y_2(x) ])/(det[ y_1(p) y_2(p) ; y_1'(p) y_2'(p) ]) = (det[ cos(2p) sin(2p) ; cos(2x) sin(2x) ])/(det[ cos(2p) sin(2p) ;-2sin(2p) 2cos(2p) ])

Ricordando come si calcola il determinante di una matrice abbiamo

det[ cos(2p) sin(2p) ; cos(2x) sin(2x) ] = cos(2p)sin(2x)-cos(2x)sin(2p)

det[ cos(2p) sin(2p) ;-2sin(2p) 2cos(2p) ] = 2cos^2(2p)+2sin^2(2x) = 2[cos^2(2p)+sin^2(2p)] = 2

Ricordiamo infatti che, per la relazione fondamentale della trigonometria

cos^2(2p)+sin^2(2p) = 1

Il nostro nucleo risolvente è quindi

K(x,p) = (cos(2p)sin(2x)-cos(2x)sin(2p))/(2) = (1)/(2)[cos(2p)sin(2x)-cos(2x)sin(2p)]

Di conseguenza la soluzione particolare è data da

bary(t) = ∫_(x_0)^(x)(K(x,p)·f(p))dp =

= ∫_(x_0)^(x)[(1)/(2)[cos(2p)sin(2x)-cos(2x)sin(2p)]]·sin(2p)dp

Scegliamo x_0 = 0 che appartiene al dominio dell'equazione differenziale data e svolgiamo il prodotto. Sfruttando poi l'additività e la linearità dell'integrale abbiamo

bary(x) = (1)/(2)sin(2x)∫_(0)^(x)[cos(2p)sin(2p)]dp (•)-(1)/(2)cos(2x)∫_(0)^(x)[sin^2(2p)]dp (• •)

Armiamoci ora di tanta pazienza e buona volontà e procediamo con il calcolo dei due integrali.

∫_(0)^(x)[cos(2p)sin(2p)]dp

Procediamo con il metodo di sostituzione per gli integrali ponendo

sin(2p) = t

da cui

2cos(2p) dp = dt

Inoltre per

p = 0, t = 0 e per p = x, t = sin(2x)

Allora

∫_(0)^(x)[cos(2p)sin(2p)]dp = (1)/(2)∫_(0)^(x)[2cos(2p)sin(2p)]dp = (1)/(2)∫_(0)^(sin(2x))[t]dt

Ci siamo così ricondotti al calcolo di un integrale fondamentale

= (1)/(2)[(t^2)/(2)]_(0)^(sin(2x)) = (1)/(4)sin^2(2x)

Ossia

∫_(0)^(x)[cos(2p)sin(2p)]dp = (1)/(4)sin^2(2x)

Allora il primo pezzo della soluzione particolare, quello che abbiamo indicato con (•) è:

(1)/(2)sin(2x)∫_(0)^(x)[cos(2p)sin(2p)]dp = (1)/(2)sin(2x)·(1)/(4)sin^2(2x) = (1)/(8)sin^3(2x)

Passiamo ora al calcolo del seguente integrale

∫_(0)^(x)[sin^2(2p)]dp

Poniamo 2p = t, da cui

2dp = dt ossia dp = (1)/(2)t

Sostituendo nell'integrale di partenza avremo

∫_(0)^(x)[sin^2(2p)]dp = (1)/(2)∫_(0)^(2x)[sin^2(t)]dt

Ricordando che l'integrale del seno al quadrato - (click per i due metodi di calcolo) è dato da

∫[sin^2(t)]dt = (t)/(2)-(1)/(2)sin(t)cos(t)+c

Abbiamo che

(1)/(2)∫_(0)^(2x)[sin^2(t)]dt = (1)/(2)[(2x)/(2)-(1)/(2)sin(2x)cos(2x)] = (1)/(2)x-(1)/(4)sin(2x)cos(2x)

Ossia

∫_(0)^(x)[sin^2(2p)]dp = (1)/(2)x-(1)/(4)sin(2x)cos(2x)

Ne segue che il secondo pezzo della soluzione particolare, quello che abbiamo indicato con (• •) sarà

(1)/(2)cos(2x)∫_(0)^(x)[sin^2(2p)]dp = (1)/(2)cos(2x)[(1)/(2)x-(1)/(4)sin(2x)cos(2x)] = (1)/(4)xcos(2x)-(1)/(8)sin(2x)cos^2(2x)

Possiamo quindi concludere che

bary(x) = (1)/(2)sin(2x)∫_(0)^(x)[cos(2p)sin(2p)]dp-(1)/(2)cos(2x)∫_(0)^(x)[sin^2(2p)]dp =

= (1)/(8)sin^3(2x)-(1)/(4)xcos(2x)+(1)/(8)sin(2x)cos^2(2x)

Per facilitare i conti che verranno osserviamo che

(1)/(8)sin^3(2x)+(1)/(8)sin(2x)cos^2(2x) =

(effettuando un raccoglimento totale e ricordando nuovamente la relazione fondamentale della trigonometria)

= (1)/(8)sin(2x)[sin^2(2x)+cos^2(2x)] = (1)/(8)sin(2x)

Morale della favola

bary(x) = (1)/(8)sin(2x)-(1)/(4)xcos(2x)

L'integrale generale dell'equazione differenziale di partenza è quindi

y(x) = y_O(x)+ bary(x) = c_1cos(2x)+c_2sin(2x)+(1)/(8)sin(2x)-(1)/(4)xcos(2x)

che possiamo riscrivere come

y(x) = (c_1-(1)/(4)x)cos(2x)+(c_2+(1)/(8))sin(2x)

Calcoliamo ora la derivata prima. Applicando le regole di derivazione e ricordando le derivate fondamentali otteniamo

y'(x) = -(1)/(4)cos(2x)-2(c_1-(1)/(4)x)sin(2x)+2(c_2+(1)/(8))cos(2x)

Imponiamo le condizioni iniziali. Da y(0) = 0 si ottiene

0 = c_1cos(0)+(c_2+(1)/(8))sin(0) = c_1

ossia c_1 = 0. Da y'(0) = 1 abbiamo invece

1 = -(1)/(4)cos(0)-2(c_1-(1)/(4))sin(0)+2(c_2+(1)/(8))cos(0)

1 = -(1)/(4)+2c_2+(1)/(4) ⇔ c_2 = (1)/(2)

Abbiamo finito! La soluzione del Problema di Cauchy di partenza è

y(x) = -(1)/(4)xcos(2x)+(5)/(8)sin(2x)

emt
Ringraziano: Omega, Ifrit, CarFaby, ferdcip, castro87
  • Pagina:
  • 1
Os