Laplace getransformeerde van de matrix-vergelijking geeft:
(sX[s] - x[0] ; sY[s] - y[0]) = (5 3;-6 -4)(X[s] ; Y[s]) + (1/s;1/s)
Beginvoorwaarden invullen:
(sX[s] ; sY[s]) = (5 3;-6 -4)(X[s] ; Y[s]) + (1/s ; 1/s)
Wat links staat is te schrijven als een diagonaalmatrix*vector:
(sX[s] ; sY[s]) = (s 0; 0 s)(X[s] ; Y[s])
Hierdoor kan je de matrices in elkaar schuiven:
(s-5 -3; 6 s+4)(X[s] ; Y[s]) = (1+1/s ; 1+1/s)
Je ziet dat de twee rijen aan elkaar gelijk zijn, dus:
(s-5)X[s] - 3Y[s] = 6X[s] + (s+4)Y[s]
Oplossen naar Y[s]:
(s-11)X[s] = (s+7)Y[s]
Y[s] = (s-11)/(s+7) X[s]
Invullen in de eerste rij, en gelijk stellen aan 1+1/s (=(s+1)/s):
(s-5)X[s] - 3(s-11)/(s+7) X[s] = 1 + 1/s
X[s] *( s-5 - (3s-33)/(s+7) ) = 1 + 1/s
X[s] *( (s2 - 35 +2s -3s + 33)/(s+7) ) = 1 + 1/s
X[s] *( (s2 -s + 2)/(s+7) ) = 1 + 1/s
X[s] *( (s-2)(s+1)/(s+7) ) = (s+1)/s
X[s] = (s+7)/(s(s-2))
X[s] = 1/(s-2) + 7/(s(s-2))
Het kan nog iets mooier, want 1/(s(s-2)) = 1/2( 1/(s-2) -1/s ), waardoor:
X[s] = 1/(s-2) + +7/(2(s-2)) -7/(2s)
X[s] = 9/(2(s-2)) - 7/(2s)
Terugtransformeren geeft inderdaad het goede antwoord voor x[t]. y[t] kan je nu ook wel vinden denk ik. Misschien dat het fout is gegaan bij het oplossen naar X[s]?
Laatst gewijzigd op 08-11-2005 om 16:50.
|