Continuous-time Dynamical System - Periodic Trajectory

Calculus Level 5

Let v ( t ) = [ x ( t ) y ( t ) z ( t ) ] \mathbf{v}(t) = \begin{bmatrix} x(t) \\ y(t) \\ z(t) \end{bmatrix}

If

d v ( t ) d t = A v ( t ) \dfrac{d\mathbf{v}(t)}{dt} =A \mathbf{v}(t)

where

A = [ 0 0.809016994 0.50903696 0.809016994 0 0.293892626 0.50903696 0.293892626 0 ] A = \begin{bmatrix} 0 && -0.809016994 && 0.50903696 \\ 0.809016994 && 0 &&-0.293892626 \\ -0.50903696 && 0.293892626 && 0 \end{bmatrix}

And v ( 0 ) = [ 5 0 0 ] \mathbf{v}(0) = \begin{bmatrix} 5 \\ 0 \\0 \end{bmatrix}

Then v ( t ) \mathbf{v}(t) will trace a circle in 3D space. Find the radius of this circle.


The answer is 4.779.

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.

1 solution

Karan Chatrath
Sep 22, 2020

Another good problem. Consider the matrix A, which is of the form:

A = [ 0 a 1 a 2 a 1 0 a 3 a 2 a 3 0 ] A = \left[\begin{matrix} 0&-a_1&a_2\\a_1&0&-a_3\\-a_2&a_3&0\end{matrix}\right]

The characteristic equation of this matrix is:

λ 3 + ( a 1 2 + a 2 2 + a 3 2 ) λ = 0 \lambda^3 + (a_1^2 + a_2^2+a_3^2)\lambda =0

It can be seen that:

a 1 2 + a 2 2 + a 3 2 = 1 a_1^2 + a_2^2+a_3^2=1

So the eigen values of A A are ( 0 , ± i ) (0,\pm i) .

The eigenvalue decomposition of A A reads:

A = V D V 1 A = V D V^{-1}

Where D D is:

D = [ i 0 0 0 i 0 0 0 0 ] D = \left[\begin{matrix} i&0&0\\0& -i&0\\0&0&0\end{matrix}\right]

The general solution of the given system is:

v ( t ) = V e D t V 1 v ( 0 ) \mathrm{v}(t) = V \mathrm{e}^{Dt} V^{-1} \mathrm{v}(0) v ( t ) = V [ e i t 0 0 0 e i t 0 0 0 1 ] V 1 v ( 0 ) \mathrm{v}(t) = V \left[\begin{matrix} \mathrm{e}^{it}&0&0\\0& \mathrm{e}^{-it}&0\\0&0&1\end{matrix}\right] V^{-1} \mathrm{v}(0)

The solution of this system is periodic when:

v ( t + T ) = v ( t ) = V [ e i ( t + T ) 0 0 0 e i ( t + T ) 0 0 0 1 ] V 1 v ( 0 ) \mathrm{v}(t +T) =\mathrm{v}(t) = V \left[\begin{matrix} \mathrm{e}^{i(t+T)}&0&0\\0& \mathrm{e}^{-i(t+T)}&0\\0&0&1\end{matrix}\right] V^{-1} \mathrm{v}(0)

The above is only possible when T = 2 π T = 2\pi . This means that a point on the trajectory which would be diametrically opposite to the initial condition v ( 0 ) \mathrm{v}(0) would be:

v ( π ) = V [ e i π 0 0 0 e i π 0 0 0 1 ] V 1 v ( 0 ) \mathrm{v}(\pi) = V \left[\begin{matrix} \mathrm{e}^{i \pi}&0&0\\0& \mathrm{e}^{-i \pi}&0\\0&0&1\end{matrix}\right] V^{-1} \mathrm{v}(0)

Therefore, the radius is:

R = v ( π ) v ( 0 ) 2 = 4.779 R = \frac{\lvert \mathrm{v}(\pi) - \mathrm{v}(0) \rvert}{2} = 4.779

What you see here is a bug in the Brilliant website. I am not deleting this for now as I will not be able to post another solution later. I will, however, update this solution at a later time.

Karan Chatrath - 8 months, 3 weeks ago

Log in to reply

Solution edited. This is the second time I have encountered a bug where all latex syntax becomes HTML-like code which is unreadable and impossible to work with.

Karan Chatrath - 8 months, 3 weeks ago

What is the matrix V V , that is, how did you calculate R R ? This method is new to me and is interesting.

A Former Brilliant Member - 8 months, 3 weeks ago

@Foolish Learner The eigen value of a square matrix λ \lambda and its corresponding eigen vector v v are defined as:

A v = λ v Av = \lambda v

Since A is 3X3, it has three eigenvalues and three corresponding eigenvectors:

A v 1 = λ 1 v 1 Av_1 = \lambda_1 v_1 A v 2 = λ 2 v 2 Av_2 = \lambda_2 v_2 A v 3 = λ 3 v 3 Av_3 = \lambda_3 v_3

Now, if we were to stack each of the eigenvectors (which are column matrices) in a matrix as such:

V = [ v 1 v 2 v 3 ] V = \left[\begin{matrix} v_1&v_2&v_3\end{matrix}\right] A V = [ A v 1 A v 2 A v 3 ] \implies AV = \left[\begin{matrix} Av_1&Av_2&Av_3\end{matrix}\right] Using the definiton of eigen values:

A V = [ λ 1 v 1 λ 2 v 2 λ 3 v 3 ] \implies AV = \left[\begin{matrix} \lambda_1v_1&\lambda_2 v_2&\lambda_3 v_3\end{matrix}\right]

A V = [ v 1 v 2 v 3 ] [ λ 1 0 0 0 λ 2 0 0 0 λ 3 ] \implies AV = \left[\begin{matrix} v_1&v_2&v_3\end{matrix}\right] \left[\begin{matrix} \lambda_1&0&0\\ 0&\lambda_2&0\\0&0&\lambda_3\end{matrix}\right]

A V = V D \implies AV = VD A = V D V 1 \implies A = VDV^{-1}

So the deifnition of V V and D D becomes clear from the above steps. As for the radius calculation, I just computed a point which is a solution to the given ODE system such that it lies diametrically opposite to the initial point. Having found that point, I just computed the half of the Euclidean distance between them to find the radius.

Karan Chatrath - 8 months, 3 weeks ago

Log in to reply

@Karan Chatrath Sir please take a look on my new note. If you solved any problem from that, please tell me I will show my attempt. Thanks in advance.

Talulah Riley - 8 months, 3 weeks ago

0 pending reports

×

Problem Loading...

Note Loading...

Set Loading...