% Initialize parameters
num_iterations = 100;
num_points = 100000; % Number of points in the initial set
% Define the initial set as a small vertical strip
S = [zeros(num_points, 1), rand(num_points, 1) * 0.1]; % Small vertical strip
% Define the transformation matrices and translation vectors for Barnsley Fern
A1 = [0 0; 0 0.16]; t1 = [0; 0];
A2 = [0.85 0.04; -0.04 0.85]; t2 = [0; 1.6];
A3 = [0.2 -0.26; 0.23 0.22]; t3 = [0; 1.6];
A4 = [-0.15 0.28; 0.26 0.24]; t4 = [0; 0.44];
% Probability weights for each transformation (Barnsley Fern)
p1 = 0.02;
p2 = 0.85;
p3 = 0.07;
p4 = 0.07;
% Iterative process
for i = 1:num_iterations
% Generate random values for selecting the transformation
rand_vals = rand(num_points, 1);
% Initialize new set
new_S = zeros(num_points, 2);
% Apply each transformation based on the random values
new_S(rand_vals < p1, :) = (A1 * S(rand_vals < p1, :)')' + repmat(t1', sum(rand_vals < p1), 1);
new_S(rand_vals >= p1 & rand_vals < p1 + p2, :) = (A2 * S(rand_vals >= p1 & rand_vals < p1 + p2, :)')' +
repmat(t2', sum(rand_vals >= p1 & rand_vals < p1 + p2), 1);
new_S(rand_vals >= p1 + p2 & rand_vals < p1 + p2 + p3, :) = (A3 * S(rand_vals >= p1 + p2 & rand_vals <
p1 + p2 + p3, :)')' + repmat(t3', sum(rand_vals >= p1 + p2 & rand_vals < p1 + p2 + p3), 1);
new_S(rand_vals >= p1 + p2 + p3, :) = (A4 * S(rand_vals >= p1 + p2 + p3, :)')' + repmat(t4',
sum(rand_vals >= p1 + p2 + p3), 1);
% Update the set
S = new_S;
end
% Create a figure for visualization
figure;
hold on;
axis equal;
% Plot the final iteration
plot(S(:,1), S(:,2), '.', 'MarkerSize', 1,'Color', 'b');
hold off;
OUTPUT with initial as vertical strip of points
