Search
 
 
  Engleski
 
 
 
Open in this window (click to change)
Forum@DeGiorgi: Početna
Forum za podršku nastavi na PMF-MO
Login Registracija FAQ Smajlići Članstvo Pretražnik Forum@DeGiorgi: Početna

ODJ - metoda varijacije konstanti
WWW:

Moja sarma
 
Započnite novu temu   Odgovorite na temu   printer-friendly view    Forum@DeGiorgi: Početna -> Kolegiji 3. godine -> Obične diferencijalne jednadžbe
Prethodna tema :: Sljedeća tema  
Autor/ica Poruka
ZELENIZUBNAPLANETIDO
SADE

Forumaš(ica)
Forumaš(ica)


Pridružen/a: 04. 03. 2004. (19:56:15)
Postovi: (54F)16
Sarma = la pohva - posuda
= 12 - 5
Lokacija: hm?

PostPostano: 16:52 pet, 24. 9. 2004    Naslov: ODJ - metoda varijacije konstanti Citirajte i odgovorite

Da li bi netko mogao ukratko opisati taj postupak u nadi da cu moci sam izvesti medjukorake?
Moze opcenito, a moze i direktna primjena na izvod rjesenja linearne diferencijalne jednadjbe?
Da li bi netko mogao ukratko opisati taj postupak u nadi da cu moci sam izvesti medjukorake?
Moze opcenito, a moze i direktna primjena na izvod rjesenja linearne diferencijalne jednadjbe?



_________________

Pupoljak nije negiran. Rekao sam to i ponovit cu to jos jedanput. Pupoljak NIJE negirAn.
MADD
(Mothers Against Dirty Dialectics)
Based on a true story. NOT.
Ko ih sljivi, mi sviramo punk Wink
[Vrh]
Korisnički profil Pošaljite privatnu poruku Posjetite Web stranice
nenad
Moderator
Moderator


Pridružen/a: 08. 10. 2002. (14:08:30)
Postovi: (355)16
Sarma = la pohva - posuda
92 = 106 - 14

PostPostano: 18:28 pet, 24. 9. 2004    Naslov: Citirajte i odgovorite

Najprije treba odrediti OPĆE rješenje HOMEGENE linearne jednadžbe.
Na primjer, za y'+py=g, homogena je jednadžba y'+py=0, i rješava se
(obično) separacijom varijabli. Rješenje je y(x) = c exp(-int p) , pri čemu
je c proizvoljna konstanta.
Lagrangeova metoda sastoji se u tome da se kao [i]ansatz[/i] (pretpostavljeni
oblik rješenja) za NEHOMOGENU jednadžbu uzme funkcija (konstantu c
smo zamijenili funkcijom c - VARIRALI smo konstantu)
y_p(x) = c(x) exp(-int p)
koja se uvrštava u nehomogenu jednadžbu i dobiva se jednadžba za c
koju možemo lako riješiti (izbor konstante kod određivanja rješenja c je
proizvoljan - trebamo samo jedno PARTIKULARNO rješenje).
Kako je jednadžba linearna, to kao u Linearnoj algebri za linearne
sustave algebarskih jednadžbi dobivamo da je opće rješenje nehomogene
jednadžbe zbroj partikularnog rješenja nehomogene jednadžbe i općeg
rješenja homogene.

Ukoliko se radi o linearnoj jednadžbi višeg (n-tog) reda, imamo fundamentalni
skup rješenja homogene jednadžbe (n linearno neovisnih funkcija koje razapinju
prostor rješenja, koji je n-dim vektorski prostor), i opće rješenje je njihova
linearna kombinacija:
y(x) = c_1 y_1(x) + ... + c_n y_n(x).
Partikularno rješenje tražimo u obliku:
y_p(x) = c_1(x) y_1(x) + ... + c_n(x) y_n(x),
gdje su c_i sada nepoznate funkcije. Uvrštavamo, deriviramo, sređujemo ...
Za n=2 trik je da se iskoristi sloboda, te se proizvoljno stavi u izrazu za
y_p' da je c_1'y_1+c_2'y_2 = 0, što bitno pojednostavljuje račun y_p''.
Slično se radi i za n>2.

Nadam se da će ovo, uz malo papira i olovku za vježbu, pomoći.

Naravno, posve je druga priča kako odrediti fundamentalni skup
rješenja (lako za n=1 - v. gore, ili za konstantne koeficijente).

- Nenad.
Najprije treba odrediti OPĆE rješenje HOMEGENE linearne jednadžbe.
Na primjer, za y'+py=g, homogena je jednadžba y'+py=0, i rješava se
(obično) separacijom varijabli. Rješenje je y(x) = c exp(-int p) , pri čemu
je c proizvoljna konstanta.
Lagrangeova metoda sastoji se u tome da se kao ansatz (pretpostavljeni
oblik rješenja) za NEHOMOGENU jednadžbu uzme funkcija (konstantu c
smo zamijenili funkcijom c - VARIRALI smo konstantu)
y_p(x) = c(x) exp(-int p)
koja se uvrštava u nehomogenu jednadžbu i dobiva se jednadžba za c
koju možemo lako riješiti (izbor konstante kod određivanja rješenja c je
proizvoljan - trebamo samo jedno PARTIKULARNO rješenje).
Kako je jednadžba linearna, to kao u Linearnoj algebri za linearne
sustave algebarskih jednadžbi dobivamo da je opće rješenje nehomogene
jednadžbe zbroj partikularnog rješenja nehomogene jednadžbe i općeg
rješenja homogene.

Ukoliko se radi o linearnoj jednadžbi višeg (n-tog) reda, imamo fundamentalni
skup rješenja homogene jednadžbe (n linearno neovisnih funkcija koje razapinju
prostor rješenja, koji je n-dim vektorski prostor), i opće rješenje je njihova
linearna kombinacija:
y(x) = c_1 y_1(x) + ... + c_n y_n(x).
Partikularno rješenje tražimo u obliku:
y_p(x) = c_1(x) y_1(x) + ... + c_n(x) y_n(x),
gdje su c_i sada nepoznate funkcije. Uvrštavamo, deriviramo, sređujemo ...
Za n=2 trik je da se iskoristi sloboda, te se proizvoljno stavi u izrazu za
y_p' da je c_1'y_1+c_2'y_2 = 0, što bitno pojednostavljuje račun y_p''.
Slično se radi i za n>2.

Nadam se da će ovo, uz malo papira i olovku za vježbu, pomoći.

Naravno, posve je druga priča kako odrediti fundamentalni skup
rješenja (lako za n=1 - v. gore, ili za konstantne koeficijente).

- Nenad.


[Vrh]
Korisnički profil Pošaljite privatnu poruku
Crni
Forumaš(ica)
Forumaš(ica)


Pridružen/a: 15. 12. 2003. (01:20:43)
Postovi: (23C)16
Spol: muško
Sarma = la pohva - posuda
= 29 - 25
Lokacija: Zagreb

PostPostano: 18:58 pet, 24. 9. 2004    Naslov: Re: ODJ - metoda varijacije konstanti Citirajte i odgovorite

[quote="ZELENIZUBNAPLANETIDOSADE"]Da li bi netko mogao ukratko opisati taj postupak u nadi da cu moci sam izvesti medjukorake?
Moze opcenito, a moze i direktna primjena na izvod rjesenja linearne diferencijalne jednadjbe?[/quote]

Pogledaj si u [i]Demidoviču[/i]. Tam' ti je to dobro objašnjeno.
ZELENIZUBNAPLANETIDOSADE (napisa):
Da li bi netko mogao ukratko opisati taj postupak u nadi da cu moci sam izvesti medjukorake?
Moze opcenito, a moze i direktna primjena na izvod rjesenja linearne diferencijalne jednadjbe?


Pogledaj si u Demidoviču. Tam' ti je to dobro objašnjeno.


[Vrh]
Korisnički profil Pošaljite privatnu poruku
ZELENIZUBNAPLANETIDO
SADE

Forumaš(ica)
Forumaš(ica)


Pridružen/a: 04. 03. 2004. (19:56:15)
Postovi: (54F)16
Sarma = la pohva - posuda
= 12 - 5
Lokacija: hm?

PostPostano: 19:19 pet, 24. 9. 2004    Naslov: Citirajte i odgovorite

Hvala Nenad :D i tebi crni :)

na posao :) pa ako ispadne neki tex iz toga zavrsi ovdje :)
Hvala Nenad Very Happy i tebi crni Smile

na posao Smile pa ako ispadne neki tex iz toga zavrsi ovdje Smile



_________________

Pupoljak nije negiran. Rekao sam to i ponovit cu to jos jedanput. Pupoljak NIJE negirAn.
MADD
(Mothers Against Dirty Dialectics)
Based on a true story. NOT.
Ko ih sljivi, mi sviramo punk Wink
[Vrh]
Korisnički profil Pošaljite privatnu poruku Posjetite Web stranice
ZELENIZUBNAPLANETIDO
SADE

Forumaš(ica)
Forumaš(ica)


Pridružen/a: 04. 03. 2004. (19:56:15)
Postovi: (54F)16
Sarma = la pohva - posuda
= 12 - 5
Lokacija: hm?

PostPostano: 11:20 sub, 25. 9. 2004    Naslov: Citirajte i odgovorite

Polazimo od C-problema sa linearnom diferencijalnom jednadjbom:
[latex]\displaystyle \frac{dx}{dt}=a(t)x+b(t) \\ x(t_0)=x_0[/latex]
[quote="nenad"]Najprije treba odrediti OPĆE rješenje HOMEGENE linearne jednadžbe.[/quote]
Dakle:
[latex]\displaystyle $(1) $\frac{dx}{dt}=a(t)x
\\$(2) $ x(t_0)=c~ , ~c \in R[/latex]
Diferencijalna jednadjba (1) je diferencijalna jednadjba sa separiranim varijablama:
[latex]\frac{dx}{dt}=g(t)h(x)$, za $\\
g(t):=a(t) ~ , ~t \in I$ i $ h(x)=x ~ , ~x \in R[/latex]
Pri cemu su obije fje klase C1, dakle mozemo primijeniti zakljucke iz poglavlja o SEP-ODJ. Tocnije, radi se o teoremu 2.4., koji kaze da tada nas (C)-problem (1)+(2) ima jedinstveno, neprosirivo rijesenje i za c=0 je to rjesenje globalno i dano formulom u(t)=0, za t iz I.
Buduci da je h(0)=0 imamo dvije komponente povezanosti rjesenja (|R+ i |R-).
Zbog antisimetricnosti toga preslikavanja sa obzirom na h(x)=x, imamo skoro-pa-isti racun:
U skladu sa oznakama iz poglavlja o SEP-ODJ, def fju [latex]G : I \rightarrow R[/latex] formulom:
[latex]\displaystyle \int_{t_0}^t a(\tau)d\tau~ , ~t \in I[/latex]
te za c>0 fju H:
[latex]\displaystyle \int_c^x \frac{ d \xi }{ h( \xi ) } = \int_c^x \frac{ d \xi }{ \xi }=\ln x-\ln c = \ln \frac xc~ , ~x \in R^+ \\
\Rightarrow H^{-1}(y)=ce^y~ , ~y \in R[/latex]
Dakle (po tm2.4):
[latex]\displaystyle u(t)=H^{-1}(G(t))~ , ~t \in I \Rightarrow u(t)=ce^{\int_{t_0}^t a(\tau)d\tau}~ , ~t \in I[/latex]

Za slucaj c<0 imamo potpuno isti slucaj, dakle, uz oznaku:
[latex]U(t,s)=\int^t_s a(\tau)d\tau~ , ~t,s \in I[/latex]
slijedi da je za bilo koji c iz R rijesenje C-problema dano formulom:
[latex]u(t)=c \cdot U(t, t_0)~ , ~t \in I[/latex]

Izvod formule za rijesenje gornjeg C-problema bazira se na dosjetki zvanoj [b]METODA VARIJACIJE KONSTANTI[/b] koja se bazira na pretpostavci da je c funkcija od t, tj:
[latex](15)$ $u(t)=c(t) \cdot U(t, t_0)~ , ~t \in I[/latex]
Da bi odredili tu fju c(t) treba osigurati da u(t) zadovoljava uvjete C-problema, tj:
[latex]\displaystyle \frac d{dt}c(t)U(t, t_0)=a(t)c(t)U(t,t_0)+b(t) $ i $\\
c(t_0)U(t,t_0)=x_0[/latex]
Raspisemo diferencijal:
[latex]\displaystyle \frac d{dt}c(t)U(t,t_0)=\frac{dc}{dt}(t)U(t,t_0)+c(t)\partial_1U(t,t_0)[/latex]
..koristeci se rezultatima leme 3.3 o svojstvima fje U(t,s) dobivamo jednadjbe:
[latex]\displaystyle \frac{dc}{dt}(t)U(t,t_0)=b(t) $ i $ c(t_0)=x_0[/latex]
mnozimo fju sa [latex]\frac 1{U(t,t_0)}=U(t_0,t)[/latex] i dobivamo:
[latex]\displaystyle \frac{dc}{dt}(t)=U(t_0,t)b(t)~ , ~t \in I[/latex]
Integracijom:
[latex]\displaystyle c(t)-c(t_0)=\int_{t_0}^tU(t_0, \tau)b(\tau)d\tau \\
\Rightarrow c(t)=\int_{t_0}^tU(t_0, \tau)b(\tau)d\tau+x_0[/latex]

Uvrstavajuci u pocetnu pretpostavku (15) dobivamo:
[latex]\displaystyle u(t)=c(t)U(t,t_0)=U(t,t_0)\left[ x_0 + \int_{t_0}^t U(t_0, \tau)b(\tau)d\tau \right][/latex]
..pa mnozeci elemente u zagradama, opet na osnovu leme 3.3 dobivamo:
[latex]\displaystyle u(t)=U(t,t_0)x_0+\int^t_{t_0}U(t,\tau)b(\tau)d\tau~ , ~t \in I[/latex]
...sto je rjesenje iz zakljucka teorema 3.2 o rijesenju c-problema sa L-ODJ

phew :) eto :)

[quote="nenad"][i]ansatz[/i][/quote]
samo ovo ne razumijem :shock: :?
Polazimo od C-problema sa linearnom diferencijalnom jednadjbom:

nenad (napisa):
Najprije treba odrediti OPĆE rješenje HOMEGENE linearne jednadžbe.

Dakle:

Diferencijalna jednadjba (1) je diferencijalna jednadjba sa separiranim varijablama:

Pri cemu su obije fje klase C1, dakle mozemo primijeniti zakljucke iz poglavlja o SEP-ODJ. Tocnije, radi se o teoremu 2.4., koji kaze da tada nas (C)-problem (1)+(2) ima jedinstveno, neprosirivo rijesenje i za c=0 je to rjesenje globalno i dano formulom u(t)=0, za t iz I.
Buduci da je h(0)=0 imamo dvije komponente povezanosti rjesenja (|R+ i |R-).
Zbog antisimetricnosti toga preslikavanja sa obzirom na h(x)=x, imamo skoro-pa-isti racun:
U skladu sa oznakama iz poglavlja o SEP-ODJ, def fju formulom:

te za c>0 fju H:

Dakle (po tm2.4):


Za slucaj c<0 imamo potpuno isti slucaj, dakle, uz oznaku:

slijedi da je za bilo koji c iz R rijesenje C-problema dano formulom:


Izvod formule za rijesenje gornjeg C-problema bazira se na dosjetki zvanoj METODA VARIJACIJE KONSTANTI koja se bazira na pretpostavci da je c funkcija od t, tj:

Da bi odredili tu fju c(t) treba osigurati da u(t) zadovoljava uvjete C-problema, tj:

Raspisemo diferencijal:

..koristeci se rezultatima leme 3.3 o svojstvima fje U(t,s) dobivamo jednadjbe:

mnozimo fju sa i dobivamo:

Integracijom:


Uvrstavajuci u pocetnu pretpostavku (15) dobivamo:

..pa mnozeci elemente u zagradama, opet na osnovu leme 3.3 dobivamo:

...sto je rjesenje iz zakljucka teorema 3.2 o rijesenju c-problema sa L-ODJ

phew Smile eto Smile

nenad (napisa):
ansatz

samo ovo ne razumijem Shocked Confused



_________________

Pupoljak nije negiran. Rekao sam to i ponovit cu to jos jedanput. Pupoljak NIJE negirAn.
MADD
(Mothers Against Dirty Dialectics)
Based on a true story. NOT.
Ko ih sljivi, mi sviramo punk Wink
[Vrh]
Korisnički profil Pošaljite privatnu poruku Posjetite Web stranice
veky
Forumaš(ica)
Forumaš(ica)


Pridružen/a: 09. 12. 2002. (19:59:43)
Postovi: (5B0)16
Sarma = la pohva - posuda
22 = 24 - 2
Lokacija: negdje daleko...

PostPostano: 16:33 pon, 27. 9. 2004    Naslov: Citirajte i odgovorite

[quote="ZELENIZUBNAPLANETIDOSADE"]
[quote="nenad"][i]ansatz[/i][/quote]
samo ovo ne razumijem :shock: :?[/quote]

Nešto kao "za početak", odnosno kao početna pretpostavka. IIRC.
ZELENIZUBNAPLANETIDOSADE (napisa):

nenad (napisa):
ansatz

samo ovo ne razumijem Shocked Confused


Nešto kao "za početak", odnosno kao početna pretpostavka. IIRC.


[Vrh]
Korisnički profil Pošaljite privatnu poruku Pošaljite e-mail Posjetite Web stranice
Prethodni postovi:   
Započnite novu temu   Odgovorite na temu   printer-friendly view    Forum@DeGiorgi: Početna -> Kolegiji 3. godine -> Obične diferencijalne jednadžbe Vremenska zona: GMT + 01:00.
Stranica 1 / 1.

 
Forum(o)Bir:  
Možete otvarati nove teme.
Možete odgovarati na postove.
Ne možete uređivati Vaše postove.
Ne možete izbrisati Vaše postove.
Ne možete glasovati u anketama.
You cannot attach files in this forum
You can download files in this forum


Powered by phpBB © 2001, 2002 phpBB Group
Theme created by Vjacheslav Trushkin
HR (Cro) by Ančica Sečan