Similarly as time proceeds to negative infinity the contribution of the first vector goes to 0, and so vt begins to look like a multiple of the second eigenvector. Hence the vector vt begins to look more and more like a multiple of the first eigenvector.
This means as time proceeds to infinity, the contribution of the second vector goes to 0 (since (1/2)^t -> 0 as t->infinity). In this case the first eigenvalue is greater than 1 and the second is less than 1. Then the state of the system at time t is given by This happens for a simple algebraic reason: if we express our vector as a linear combination of the eigenvectors Notice that as time progresses into the future the vectors in our system are starting to hug the line spanned by the eigenvector with eigenvalue 11/10, whereas as we move into the past the vectors in our system hug the line spanned by the eigenvector with eigenvalue 1/2. Now I'd like to see all these images together at once: we'll watch the path our initial condition defines from time t=-7 to time t=55. The picture below (which should be animated) starts at time t=-7 and progresses to time t=-1. I'll do this by just taking negative powers of my matrix. Now I'd like to see where my initial condition came from, which is to say I'd like to see where the point was in the past. I'll get Mathematica to do this by producing a table of images, and it will even animate the images for me! The animation is below. What I'd like to do now is plot the trajectory of the initial condition as time progresses forward. You can see that the point looks like its right on the origin, but this is only because my picture is blown up pretty large. I'll define v0=(1,1) and will plot it on the axes above. Now I'm ready to give Mathematica my initial vector. Notice that the plot Mathematica generates has a weird scale so that the graph appears wider than it is long (even though the y values go between -300 and 300 whereas the x values only go between -100 and 100). I'll go ahead and plot the lines spanned by these vectors. The output tells me that the two eigenvalues are 11/10 and 1/2, with a basis for the corresponding eigenspaces given by the vectors (-1,1) and (1/3,1). Mathematica will give me this information when I run the command Eigensystem. To plot the trajectory of my initial vector I'll need to know eigenvalues and eigenvectors. We'll start with a matrix A and see what it does to an initival vector in. Plotting Trajectories in Linear Dynamical Systems
0 kommentar(er)
