Uitschrijven euler (direct)

Moderators: ArcherBarry, Fuzzwood

Reageer
Gebruikersavatar
Berichten: 94

Uitschrijven euler (direct)

Goedendag,

Ik probeer te uit te zoeken hoe ik de Euler (direct) kan uitschrijven om de grafiek te kunnen tekenen hiervoor.

Hieronder heb ik de gegevens

Situatieschets:

Een massa van 10kg ligt op een wrijvingsloze vlak, een veer
\( k_{veer} = 1 [N/m] \)
en een demper
\( k_{demper} = 1/2 [N/m^{2}] \)
die naast elkaar zijn aan de massa bevestigd. Men trekt de massa
\( y = 1 [m] \)
.
\(u(t)=0\)
\( \forall \)
\( t<1 \)
\(u(t)=1\)
\( \forall \)
\( t \geq 1 \)
Ik kom dmv een simulatiediagram uit op:
\( x[t+h] = \phi (h) * \Delta (t) + \rho (h) * u (t) \)
\( y [h] = C * X(t) + D * u (t) \)
\( h=0,1 \)
\( \phi = \begin{array}{cc} 1 & 0.1 \\ -0.01 & 0.995 \end{array} \)
\( \rho = \begin{array}{cc} 0.01 \\ 0.01 \end{array} \)
\( C = \begin{array}{cc} 1 & 0 \end{array} \)
\( D = 0 \)
\( X_1 = 1*x_{1}(t) + 0.1 x_{2} (t) \)
\( X_2 = -0.01*x_{1}(t) + 0.1 x_{2} (t) \)
Zou iemand mij op het goede spoor kunnen zetten? Ik wil per tijdsperiode (h=0,1) uitrekenen wat waarde u(t) is, zodat ik deze in een grafiek kan zetten, alleen heb ik zoveel uitgerekend dat ik door de bomen het bos niet meer zie.

Alvast bedankt.

Gebruikersavatar
Berichten: 94

Re: Uitschrijven euler (direct)

Ik heb nog even nagerekend met Euler(exact).
\( e^{ \left[ \begin{array}{rc1} 0 & 1 & 0 \\ -1/10 & 1/20 & 1/10 \\ 0 & 0 & 0 \end{array} \right] * 0,1 } = \left[ \begin{array}{rc1} 0,9995 & 0,0997 & 0,0004 \\ -0,0099 & 0,9945 & 0,009 \\ 0 & 0 & 1 \end{array} \right] \)
Het komt overeen met de gevonden antwoorden, dus de methode klopt in ieder geval.

Nu alleen nog de vraag hoe ik de "directe" methode kan omzetten naar een grafiek.

Reageer