program zasadniczo dział dobrze, ale okazało się, że muszę teraz założyć że po iteracjach mam osiągnięty( wyznaczony) punkt, do którego teraz atraktor musi zbiegac, musze zrobić jakąś pętle zewnętrzną...
już teraz totalnie nie wiem o co chodzi? czy ktoś by pomógł?Dziś jest ostateczny dzień na ten program stąd poszukuję pomocy

Kopiuj
% Dana jest funkcja F: (x,z) --> (xx,zz) zdefiniowana układem równań
%           xx = x + a - b*(lambda*z + (1 - lambda)*sin(2*pi*x)
%           yy = lambda*z - (1 - lambda)*sin(2*pi*x)
% W porównaniu z oryginalną treścią użyłam:
%           x   - zamiast beta
%           xx  - zamiast beta'
%           zz  - zamiast z'.
% Ponadto w treści zadania jest zastrzeżenie, że x = x(mod 1),
% tym samym dziedziną odwzorowania F jest zbiór [0,1)xR.
% Rozumiem, że zadanie polega na wykreśleniu obrazu pasa |z| <= 1,
% więc F([0,1)x[-1,1])

clc
clear
% Parametry układu równań a, b, lambda:
 
a = 0.8;
b = 1.4;
lambda = 0.5;
% Powyższe parametry można zmieniać,
% obserwując jaki ma to wpływ na kształt atraktora
%
% W programie można narysować skończoną ilość punktów.
% Wybieram tu siatkę punktów o gęstości dx i dy
dx = 0.01;
dz = 0.01;
%
n = 0;  % indeks kolejnych, rozpatrywanych punktów
%
% Kolejno wybieram pary argumentów:
%       x - (dotyczy beta) celowo z przedziału [0,3) co 0.01,
%            by zademonstrować sklejanie mod 1, tym samym np. 2.13 jest
%            utożsamiane z 0.13 dzięki zastosowaniu funkcji mod(x,1),
%       z - z przedziału [-1,1]
n=1;

    x=0;
    z=0;
  
for i = 0:200
    for j = 0:200
%         x = i * dx;         % kolejny x
%         z = - 1 + j * dz;   % kolejny z
        % w wektorze xx oraz zz zapiszemy "nowe" x i "nowe" z
        n = n + 1;  % zgodnie z wybraną siatką n = 1,2...(301*201)
                    % tyle punktów atraktora wykreślam
        % Zgodnie z układem równań obliczam (xx,zz) = F(x,z)
        xx(n) = mod(x, 1) + a - b * (lambda * (z) + (1 - lambda) * sin(2 * pi * (mod(x, 1))));
        xx(n) = mod(xx(n), 1);  % wartości x sklejają się, więc mod(xx,1)
        zz(n) = lambda * z + (1 - lambda) * sin(2 * pi * (mod(x, 1)));
x=xx(n);
z=zz(n);

    end
end
%
% Przygotowałam 60051 punktów atraktora:
%       w wektorze xx odcięte tych punktów,
%       w wektorze zz rzędne tych punktów.
%
plot(xx, zz, 'r.','MarkerSize',1);     % wykreślam punkty atraktora kolorem czerwonym (r) markerem kropką (.)
% 
% W osi x (dotyczy beta) wybrałam zakres [0,3] (i=0:300), by zademonstrować
% sklejanie się punktów x mod 1.
% Taki sam efekt da ograniczenie się do przedziału [0,1) (i=0:299).