program_2c.m

《Dynamical Systems with Applications using MATLAB》源代码

% Chapter 2 - Nonlinear Discrete Dynamical Systems.
% Program 2c - Computing a Lyapunov exponent for the logistic map.
% Copyright Birkhauser 2004. Stephen Lynch.

% Lyapunov exponent when mu=4 (Table 2.1).
clear;
mu=4;
x=0.1;xo=x;
itermax=49999;
    for n=1:itermax
        xn=mu*xo*(1-xo);
        x=[x xn];
        xo=xn;
    end

Liap_exp=vpa(sum(log(abs(mu*(1-2*x))))/itermax,6)

% End of Program 2c.