Abstract
This paper investigates the characteristics of the Kalman filter for a broad class of complex Fibonacci systems and represents an extension to the complex domain of the state estimation problem for the real-valued Fibonacci system. Complex Fibonacci systems are obtained by modifying the real-valued Fibonacci recurrence relation to include complex coefficients, control and noise inputs, and a noisy output-measurement equation. Analytic expressions for the Kalman filter’s steady-state gain and error covariance matrices are obtained, and it is found that for a broad subclass of these complex systems the elements of the matrices are functions of the golden ratio.
1. Introduction
The th number in the real-valued Fibonacci sequence can be generated either recursively from the recurrence relation (starting with and ) or directly using Binet’s formula, where , the golden ratio, is given by . The Fibonacci numbers and the golden ratio are well known for their number-theoretic properties and for their appearances in natural phenomena and art, and they have also appeared more recently in a state estimation context [1–3]. In [1, 2], connections were made between the Fibonacci numbers, the golden ratio, and properties of the Kalman filter for certain classes of one- and multidimensional systems. In [3], the properties of the filter for the Fibonacci system itself—defined by adding a control input and a noisy measurement equation to the Fibonacci recurrence relation—were also found to involve the Fibonacci numbers and the golden ratio.
This paper continues in the vein of [3] and extends the results presented there to a broad class of suitably defined complex Fibonacci systems (which include, as a special case, the real-valued Fibonacci system studied in [3]). Complex systems are obtained by introducing complex coefficients and a noise and control input into the classical Fibonacci recurrence relation (1), along with a noisy measurement equation. The resulting systems’ impulse responses—for different complex coefficients—are different complex-valued Fibonacci sequences. Analytic expressions for the Kalman steady-state gain and a priori error covariance matrices are developed, and for a broad subclass of these systems it is shown that the elements of these matrices are functions of the golden ratio.
While the real-valued Fibonacci sequence is very well known, complex Fibonacci sequences are lesser known, although their number-theoretic properties have been investigated by several authors. Studies of complex sequences generated from the classical recurrence relation (1) either by starting with complex numbers rather than with 0 and 1 or by introducing complex coefficients into the recurrence can be found in [4–8], and studies of the properties of two-dimensional complex sequences can be found in [4–6, 9]. In this paper, the number-theoretic properties of complex sequences are not studied. Rather, by modifying the classical recurrence relation to create complex Fibonacci systems whose impulse responses are complex Fibonacci sequences, the recurrence relations themselves are cast into systems context and the properties of their Kalman filters are investigated.
2. Representations of Complex Fibonacci Systems
The complex Fibonacci systems considered here are obtained by modifying the Fibonacci recurrence relation given by (1) to include a control input (a noise input will be added later) and nonunity (and possibly complex) coefficients. Specifically, the systems considered are the second-order discrete-time systems with control input , output , and difference equation (as will be shown below, and are the locations of the system’s poles) For different choices of and , the system’s impulse response (which results when and there are zero initial conditions, that is, will be different Fibonacci sequences. For example, if and (i.e., and , where is the golden ratio), then the impulse response is the classical, real-valued Fibonacci sequence defined above; if (where ) and , the impulse response is one of the complex Fibonacci sequences whose number-theoretic properties were studied in [8]. A subclass of systems introduced here and which is of particular interest to this work is one whose poles are at and . These complex (not conjugate) poles are at the same distances from the origin as the poles and of the classical, real-valued Fibonacci system, but lie along the ray that makes an angle with the real axis. The classical system is the special case of . As shown below, this pair of complex poles gives rise to a subclass of systems whose impulse responses are complex Fibonacci sequences with magnitudes equal to the real-valued Fibonacci numbers but which rotate in the complex plane.
The system transfer function for the general complex Fibonacci system can be obtained by z-transforming (3) to yield Expanding in partial fractions yields and since the transfer function is the -transform of the impulse response, taking the inverse -transform of each term in (5) gives If the poles and for the real-valued Fibonacci system are substituted into this expression, becomes Binet’s formula given by (2) for calculating the real-valued Fibonacci numbers. If the generalized poles and for the complex-valued Fibonacci system mentioned above are substituted, becomes that is, the impulse response of this system is a complex-valued Fibonacci sequence whose magnitudes are the real-valued Fibonacci numbers and whose angles are .
It is assumed in what follows that one pole of the system given by (3) is outside the unit circle and one is inside. The impulse responses of different systems in the class of Fibonacci systems defined by (3) are different versions of Fibonacci sequences, sequences which grow in magnitude, so one of the poles must be unstable. It will be assumed here that .
If the states for the system in (3) are defined to be and , and if the system is now further modified to include an excitation noise input and an output equation with measurement noise , that is, if (3) becomesthen a state description for the system is where
3. The Kalman Filter for the Zero-Input-Noise System
The Kalman filter design considered in this section is for the case where the complex Fibonacci system given by (9a), (9b), and (9c) has measurement noise —assumed to be white with variance —but no excitation noise (i.e., ). Letting , and represent, in turn, the Kalman gain, the a priori state-estimate error covariance (the error covariance at time before the measurement is used), and the a posteriori state-estimate error covariance (the error covariance after the measurement at is used), then (following [10, equations (5.19) and (5.26)], [11, Appendix 2.C], and noting that because the matrix can now be complex, its conjugate also appears in the Riccati equation) the a priori error covariance and filter gain are given by (where is the transpose of the complex conjugate of ) If the unknown initial state at is assumed to have error covariance , then the initial a priori error covariance at (to be used in (10a)) is given by . It is assumed here that the first measurement is obtained at .
The steady-state error covariance is the solution to (10a) with both and replaced by , and the steady-state gain follows from (10b). Using the Hamiltonian matrix approach as described in [12, page 203] to solve (10a) results in (which can be verified most easily by direct substitution into (10a) and (10b)). That these results depend only on the unstable pole is not surprising since as time goes on the system’s response is dominated by that pole. If attention is restricted to any Fibonacci system—real or complex—of the form (3) whose unstable pole is anywhere in the -plane at a distance from the origin, that is, , and is inside the unit circle, then the elements of the steady-state gain and error covariance will be functions of the golden ratio. Specifically, in this case (11) becomes This result is a more general, complex version of a similar result found in [3] for the classical, real-valued Fibonacci system (which results here if ), and it makes it clear that for the real-valued system too the matrices depend only on the location of the system’s unstable pole at .
Finally, if attention is further restricted to the subclass of systems whose poles are complex generalizations of the real pole locations of the real-valued system, and whose impulse responses are given by (7), that is, systems whose poles are and , then it can be shown that under certain restrictions on the initial a posteriori error covariance and the measurement noise variance, the elements of the a priori error covariance are functions of the real-valued Fibonacci numbers. Specifically, if the initial a posteriori error covariance is selected to be the identity matrix, and the measurement variance , then successive iterations of (10a) show that where , and , given by (1), is the th real-valued Fibonacci number. Since , it is seen that either expression for leads to in the limit. Again, this result is a more general, complex version of a result found in [3] for the real-valued Fibonacci system. This class of complex Fibonacci systems—whose complex poles are at the same distances from the origin as the poles of the real-valued Fibonacci system—is studied in more depth in the next section, where an expression for is derived for the case of nonzero excitation noise.
4. The Kalman Filter for the , Subclass of Systems with Excitation and Measurement Noise
4.1. System Characteristics
The difference equation for this subclass of systems is, from (8a), where the matrices in its state description are, from (9a), (9b), and (9c), and where its impulse response—denoted here as—is, from (7), The magnitudes of this complex-valued Fibonacci sequence have the same values as the real-valued sequence, but as increases they rotate in the complex plane (counterclockwise for , clockwise for ). This rotating sequence can be looked upon as the result of sampling a continuous spiral that is readily obtained from (16) by replacing the integer variable by a real variable as was done in [13–15] to generalize the real-valued sequence to a generalized Fibonacci function. Specifically, if Binet’s formula from (2) is substituted into (16) for , and then is replaced by , the resulting generalized Fibonacci function is given as While the function produces the real-valued Fibonacci numbers when takes on positive integer values, is, in general, a complex-valued function. A polar plot of the function for and positive is shown in Figure 1. The magnitudes of the function at the circled points are Fibonacci numbers and correspond to integer values for ranging from 0 to 7. The function exhibits a spiral shape for both positive and negative values of and for all values of except (for the complex sequence reduces to the real-valued sequence , and the generalized function for that case has a sinus-like shape [13–15]). For large positive , approaches a logarithmic spiral, that is, a spiral whose magnitude at any angle is given by , where and are arbitrary constants. Specifically, since for large the generalized function , the angle for a particular is given as and the magnitude is given as . Using the two expressions for , one can follow [14] and write , which in turn leads to and . From (17) it can be seen that the speed with which approaches the logarithmic spiral is the speed with which approaches zero. Finally, it is interesting to note that when , the approximating logarithmic spiral is a golden spiral. For this value of , the complex Fibonacci numbers —which are samples of the spiral at integer values of —are located alternately on the real and imaginary axes, that is, the spiral goes through a quarter turn from number to number in the sequence. Examination of (16) shows that the ratio of the magnitudes of successive values of is equal to the ratio of successive values of , and so it approaches for large . Thus, for this case where , the magnitude of the spiral increases through a quarter turn by a factor of approximately , where an increase by a factor of exactly in a quarter turn is a characteristic of a golden spiral. The curve shown in Figure 1 is for this case of .
4.2. The Kalman Filter
In developing the Kalman filter for this subclass of complex Fibonacci systems—the systems given by (15)—it is assumed that the measurement noise and the excitation noise are uncorrelated, zero-mean, white noise processes with variances and , respectively. The Riccati equation for the a priori error covariance is the same as (10a) except for an additional term to account for the presence of the excitation noise . Specifically, now satisfies and expression (10b) for the gain is unaltered. The steady-state error covariance is the solution to (18) with both and replaced by , and the steady-state gain follows from (10b). Again using the Hamiltonian matrix approach [12, page 203] to solve (18) results in where and the function . That (19) is a solution to (18) is verified most easily by direct substitution and use of the fact that . The function is real, and it is not difficult to show that . When there is measurement noise but no excitation noise, , , and (19) reduces to (12), as expected. Also, if the system’s matrix and error covariance matrix are written as and , respectively (to explicitly show their dependence on the locations of the poles), and if, further, it is noted from (15) that , where is the complex conjugate of , then it follows from conjugating both sides of (18) that *. Finally, noting that the , if the excitation noise variance is large relative to the measurement noise variance, then the a priori error covariance on the first state approaches from below, and the error covariance on the second state approaches , also from below.
5. Summary and Future Work
This paper has determined analytic expressions for the Kalman filter’s steady-state gain and a priori error covariance matrices for a broad class of suitably defined complex Fibonacci systems. The impulse responses of these systems are complex Fibonacci sequences, and of particular interest to this work was the subclass of systems investigated in Section 4 whose pole locations are complex generalizations of the real pole locations of the classical, real-valued Fibonacci system. The impulse responses of this subclass of systems are Fibonacci sequences that rotate in the complex plane, and which can be viewed as being samples of continuous spirals that are approximately logarithmic. Although the goal in this paper was to study the characteristics of Kalman filters for suitably defined complex Fibonacci systems, and not to study the number-theoretic properties of the Fibonacci sequences which arose as impulse responses of those systems, it would nevertheless be interesting to investigate the properties of the spiral sequence given by (16). It may also be useful to investigate—as was done in [3]—the characteristics of controllers such as LQR and deadbeat controllers for complex Fibonacci systems.