(8.6). 8.5 .
8.5.
f (x) = xer (1 x) | (8.8) | |
f (x) = [1 (2 x 1)4] | (8.9) |
(8.8) r ~ 2 r ~ 2.7? , (8.9) r 0 1. , ? (8.8) () , . , (8.6), , , . r, r , x , , .
, . max , f (x) , .. df / dx = 0 = max. x = max dmf / dxm = 0 m < n, a dnf / dxn < 0, f (x) n - . , , (8.6) (8.8) . , (8.9) , .. d2f / dx2 = d3f / dx3 = 0, a d4f / dx4 < 0 = 1/2.
, , (8.6) . f (x, r) r, f (x, r) 4 (. 8.4).
. 8.4. f (x) r = 0.88. , r f (x) , 4.
g (x, r) = f (f (x, r)) r (. 8.5).
. 8.5. g (x) r = 0.88, x 0 = 0.5. g (x) , 1* = 0.373 x 2* = 0.512.
, g (x, r) 2 x 0 = 0.5 x 1* = 0.373 2* = 0.512. , g (x, r) f (x, r '), r ' r, f (x, r ') 2. . 8.6 f (x, r ') r ' = 0.8. , f (x, r = 0.8) 0.513 0.799. g (x, r = 0.88) . 8.5 f (x, r ' = 0.8) , . 8.6.
|
|
. 8.6. f (x) r = 0.8. f (x) , x 1* = 0.513 x 2* = 0.799, g (x) , . 8.5.
, g , , g , f g , . a, , g (x, r = 0.88) 0.512-0.373 = 0.139, f (x, r = 0.8) 0.799-0.513 = 0.286. , a = 0.286/0.139 = 2.06, . . 8.7, f (, r = 0.8) g (x, r = 0.88).
. 8.7. f (, r = 0.8) g (x, r = 0.88).
; g f, g () r. . , r ' = 0.8. , , = 1/2 , = 1/2. , , h (x) g (x) .., , f (x).
a . , . 8.3 8.8. , , . f (x) r = r (n) 1/2.
. 8.8. r. d * =1/2 2 n.
, 8.4 r (1) ≈ 0.809 r (2) ≈ 0.875. dn = xn * 1/2, n * , * = 1/2. dn . 7.8, d 1 ≈ 0.309 d 2 ≈ 0.117. , , = 1/2, = 1/2 . a
(8.10) |
a = 0.309/0.117 = 2.64 a = 2.5029078750958928485... ( , , !)
. , rn r, 2 n . 7.4 r 1 = 0.75, r 2 = 0.862 r 3 = 0.880. ( ), n rn rc :
|
|
rn rc = A d n. | (8.11) |
, d, a, , .. d f (x), . , f (x). (8.11) , d
(8.12) |
r 1, r 2 r 3 d ≈ 6.22; d = 4.66920160910299097....
8.6. a d
. , , , g (x) = 1/2 r = r (2) ≈ 0.875 f (x) r (1) ≈ 0.809. a f g.
. , . , h (x) r = r (3) ≈ 0.880 g (x) r = r (2). .
. rn 8 → 16 16 → 32. rn d.
. r (n), r, 2 n f (x) * = 1/2. d, (8.12) rn r (n).
. a d (8.8) (8.9).
, a d, f (). ? , (8.6) . , , , , . , , .
, , . ? RLC-, . , a d , .
, , . , . , , , . , . ( , ) , . , . , , , . , . , .
|
|
. , ? , , , . , , . , .
% logistic_5
clear all
for r=0.5:0.001:1
lambda=4*r;
x(1)=0.5;
for n=2:1000
x(n)=lambda*x(n-1)*(1-x(n-1));
if n>150
plot(lambda,x(n)),hold on
end
end
end
title(' \lambda')
xlabel('\lambda');
ylabel('x');
% logistic_5
clear all
for r=0.73:0.001:1
lambda=4*r;
x(1)=0.5;
for n=2:1000
x(n)=lambda*x(n-1)*(1-x(n-1));
if n>150
plot(lambda,x(n)),hold on
end
end
end
title(' \lambda')
xlabel('\lambda');
ylabel('x');
% logistic_5
clear all
for r=0.85:0.001:1
lambda=4*r;
x(1)=0.5;
for n=2:1000
x(n)=lambda*x(n-1)*(1-x(n-1));
if n>150
plot(lambda,x(n)),hold on
end
end
end
title(' \lambda')
xlabel('\lambda');
ylabel('x');
% logistic_5
clear all
for r=0.95:0.0001:0.9(75
lambda=4*r;
x(1)=0.5;
for n=2:1000
x(n)=lambda*x(n-1)*(1-x(n-1));
if n>150
plot(lambda,x(n)),hold on
end
end
end
title(' \lambda')
xlabel('\lambda');
ylabel('x');
(Henon)
( )
, , , λ. , , , , , , . , , , xn. , , λ, . , , , . :
:
, . , . . , :
. , ,
|
|
. .. , , , .
, .
,
, , . : ? , . , (n + 1)- n - . , , .. (n + 1)- n- , (n - 1)- . , , .
b .
yn + 1 = xn. (63)
. ,
(M. Henon), . , , , .
, 2-. .
,
,
, , , .. .
2- .
(70),
,
, (71) x 22 + x 12 = (x 2 + x 1)2 − 2 x 2 x 1,
, , , x 1 x 2
:
, ,
, 2-
. 38 b = 0,3. , 2- , , . , , . . (, , . . 33) , : , . , , .
. 38. (65) b=0,3.