Let v ( k ) = ⎣ ⎡ x ( k ) y ( k ) z ( k ) ⎦ ⎤
If
v ( k + 1 ) = A v ( k )
where
A = ⎣ ⎡ 0 . 9 4 4 9 0 1 5 4 3 0 . 2 8 5 7 2 2 2 2 6 − 0 . 1 5 9 7 6 1 9 6 3 − 0 . 2 6 7 6 7 7 9 9 1 0 . 9 5 5 3 1 9 3 8 6 0 . 1 2 5 3 5 2 9 5 6 0 . 1 8 8 4 3 9 8 2 6 − 0 . 0 7 5 6 8 1 4 4 0 . 9 7 9 1 6 4 3 1 3 ⎦ ⎤
And v ( 0 ) = ⎣ ⎡ 3 4 5 ⎦ ⎤
Then sequence { v ( k ) } is a periodic sequence that lies on a circle in 3D space. Find the radius of this circle.
This section requires Javascript.
You are seeing this because something didn't load right. We suggest you, (a) try
refreshing the page, (b) enabling javascript if it is disabled on your browser and,
finally, (c)
loading the
non-javascript version of this page
. We're sorry about the hassle.
Computations such as the eigenvalue decomposition and and matrix power calculations are done using a computer. It can also be done using online open source tools.
Problem Loading...
Note Loading...
Set Loading...
The general solution of this dynamic system is:
v ( k ) = A k v ( 0 )
If the solution is a circle and is periodic, then there must be a positive integer k such that:
v ( k ) = v ( 0 ) = A k v ( 0 ) ⟹ A k = I
The eigenvalue decomposition of A is can be done as such:
A = V D V − 1
A k = V D k V − 1
For the above matrix to be the identity matrix, basically:
D k = I
Here, D is already a diagonal matrix with entries:
D = ⎣ ⎡ 0 . 9 3 9 7 + 0 . 3 4 2 i 0 0 0 0 . 9 3 9 7 − 0 . 3 4 2 i 0 0 0 1 ⎦ ⎤
Let:
e i θ = 0 . 9 3 9 7 + 0 . 3 4 2 i
⟹ D = = ⎣ ⎡ e i θ 0 0 0 e − i θ 0 0 0 1 ⎦ ⎤ ⟹ D k = ⎣ ⎡ e i k θ 0 0 0 e − i k θ 0 0 0 1 ⎦ ⎤
For the above matrix to be the identity matrix:
k θ = 2 π k = 1 8
So the 19th term which is v ( 1 8 ) is when the trajectory completes a complete circle. The term in between the first and 19th term is the 10th term which is v ( 9 ) . Therefore, the radius is:
R = 2 ∣ v ( 9 ) − v ( 0 ) ∣ = 1 . 2 3 2
v ( 9 ) corresponds to
k θ = π
It is essentially the point diametrically opposite to the initial vector, on the trajectory.