In a recent post, we looked at which Toeplitz real matrices gave useful quadrature oscillators by the recurrence . There, we computed their eigenvalues and solved the recurrence in terms of them. Of course, there are many other ways to approach this problem. Here we look at another approach that gives a good geometric picture of what happens, can be applied to general matrices, and may be used as a starting point for the case.
We are interested in characterizing the matrices such that if we consider the sequence defined by , fixed, then its two components are sinusoids of constant phase difference 90º (or perhaps just approximately 90º).
We will denote by the matrix of rotation by in :
Reasoning as in the previous post, we see that the eigenvalues of need to have modulus and be complex conjugates. It follows that is similar to a rotation matrix:
In fact it makes sense and is worthy to consider the "continuous version" of the recurrence, which has the solution
We denote by the singular value decomposition of . is a diagonal matrix with strictly positive entries on the diagonal, and and are real orthogonal matrices. We note that depending on whether preserves or reverses the orientation. Hence,
We put and make the change of variables . We have
Now we can clearly see geometrically what happens. As varies along , traces the circumference of centre and radius . Then stretches this circumference to produce an ellipse (except in the trivial case when is a multiple of the identity, which implies that the we started with was just a rotation matrix). The axes of this ellipse are the coordinate axes. Then , which is a rotation and perhaps an orthogonal symmetry, just rotates this ellipse, so that its axes may no longer be the coordinate axes.
Therefore, for any , the orbit traces an ellipse. Moreover, if we take two nonzero initial data and , their corresponding orbits and , seen as sets, are proportional: . So the initial data that we choose is not important.
Of course, the components of are sinusoids of angular frequency , because they are linear combinations of the the components of , which are sinusoids of angular frequency themselves. Hence,
for some , , and . The curve traced by is an ellipse whose axes are rotated with respect to the coordinate axes.
It follows that if we want the components of to be in quadrature, which amounts to , then the axes of the ellipse traced by should be the coordinate axes. This implies that performs no rotation at all. It can only permute and change the sign of the vectors in the standard basis. Therefore, we may as well assume that is the identity, because these simple transformations of the vectors in the standard basis don't make much difference.
We recall that . Putting , denoting the diagonal elements of by , and computing the matrix product, we see that
where . We see that the only matrices that produce a quadrature oscillator are Toeplitz: precisely those which we studied in the previous post.
Perhaps one is interested in matrices that give oscillators which are almost in quadrature, meaning that is very small. Of course, this just means that is essentially a rotation of very small angle, so the matrices that satisfy this condition may be computed in a similar fashion.