The leading Lyapunov vector is computed as follows:
1) Start with an arbitrary perturbation of arbitrary size
2) Evolve it from
t to using the TLM
3) Repeat 2) for the succeeding time intervals
After a sufficiently long time , the perturbation
converges to the leading Lyapunov
vector. The direction of this vector is independent of
the initial perturbation, the
length of the time interval or the choice of norm of the
perturbation,
properties not shared by singular vectors. If during the
repeated application of
the TLM the LV becomes too large it may be scaled down to
avoid computational
blow up.
Additional LVs can be obtained by the same procedure,
except that after
each time step the perturbation has to be orthogonalized
with respect to
the subspace of the previous LVs, since otherwise all the
LVs would
converge to the leading LV.