2 svar
110 visningar
Minounderstand är nöjd med hjälpen
Minounderstand 154
Postad: 8 aug 2017 19:06 Redigerad: 8 aug 2017 19:13

Lösa ODE av grad 3 med Euler framåt

Försöker lösa denna differentialekvation med Euler framåt men vet inte riktigt hur jag ska bära mig åt algebraiskt när den är av grad >1

y''-xy=2+y' y''=xy+y'+2 med initialvärden:

y(2)=3y'(2)=2och steglängd h=0.5

De frågar alltså efter y(2.5) och y'(3)

 

Tänker att man skriver om den som ett system av första graden:

y'=vv'=xy+v+2ddxyv=01x1yv+02

y(2.5) kan jag då räkna ut med Eulers metod:
y(2.5)=y(2)+h·y'(2)=3+0.5·2=4
och y'(2.5)=v(2.5)=2+0.5(2·3+2+2)=7

men hur bär jag mig åt för att beräkna y'(3)?

 

Edit: råkade skriva grad 3, men menar förstås grad 2.

Stokastisk 3597 – Fd. Medlem
Postad: 8 aug 2017 19:17 Redigerad: 8 aug 2017 19:18

Om du har att

z˙=Az+b \mathbf{\dot{z}} = \mathbf{A}\mathbf{z} + \mathbf{b}

där z=yv \mathbf{z} = \begin{pmatrix}y \\ v\end{pmatrix} , så innebär Euler framåt att du har steget

zn+1=zn+hAzn+b \mathbf{z}_{n + 1} = \mathbf{z}_n + h \left(\mathbf{A}\mathbf{z}_n + \mathbf{b}\right) .

Så du ska alltså stega med hela vektorn yv \begin{pmatrix}y \\ v\end{pmatrix} och inte bara med y y enskilt.

Minounderstand 154
Postad: 8 aug 2017 19:41

Ah, vackert.

Tack!

Svara Avbryt
Close