Scribd
Upload
8 views3 pages

FERNwithDtAlgorithm

Uploaded by

Ashish Jakhar
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as RTF, PDF, TXT or read online on Scribd
8 views3 pages

FERNwithDtAlgorithm

Uploaded by

Ashish Jakhar
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as RTF, PDF, TXT or read online on Scribd
You are on page 1/ 3

% Set the number of iterations

num_iterations = 13;

% Call the main function to generate and plot the Barnsley Fern

barnsley_fern(num_iterations);

% Function to generate and plot the Barnsley Fern

function barnsley_fern(num_iterations)

% Initial figure: a closed equilateral triangle

S = [0, 0; 4, 0; 2, 4 * sqrt(3)/2; 0, 0]; % Triangle with vertices (0,0), (4,0), (2, height) and closing the loop
to the first point

% Define the transformation matrices and translation vectors

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];

% Initialize the figure

figure;

axis equal;

hold on;

axis off;

% Iterate to generate the Barnsley Fern

for i = 1:num_iterations
% Apply transformations deterministically

S = apply_transformations(S, A1, t1, A2, t2, A3, t3, A4, t4);

end

% Plot the final iteration

plot(S(:,1), S(:,2), 'o-', 'MarkerSize', 2, 'Color', 'b');

hold off;

end

% Function to apply transformations deterministically

function new_S = apply_transformations(S, A1, t1, A2, t2, A3, t3, A4, t4)

% Initialize new set

new_S = zeros(size(S, 1) * 4, 2); % Four times the number of original points

% Apply each transformation to the entire figure

new_S(1:size(S, 1), :) = (A1 * S')' + repmat(t1', size(S, 1), 1);

new_S(size(S, 1)+1:2*size(S, 1), :) = (A2 * S')' + repmat(t2', size(S, 1), 1);

new_S(2*size(S, 1)+1:3*size(S, 1), :) = (A3 * S')' + repmat(t3', size(S, 1), 1);

new_S(3*size(S, 1)+1:end, :) = (A4 * S')' + repmat(t4', size(S, 1), 1);

end

OUTPUT WITH 13 iteration BY starting with triangle

You might also like