Plotting Sierpinski's triangle
26 views (last 30 days)
Show older comments
Sierpinski’s triangle can be implemented in MATLAB by plotting points iteratively according to one of the following three rules which are selected randomly with equal probability.
- Rule 1: x=05.*x y=0.5*y
- Rule 2: x=0.5*x+.25 y=0.5*y+ sqrt(3)/4
- Rule 3: x=0.5x+.5 y=.5*y
Write a script that calculates the x and y vectors and then plots y versus x as individual points.
Start with x1=0 and y1=0. Run the program three times, with 100, 1000, and 10000 iterations. Print all three plots.
I do not know where to start. I am thinking if and for loops but not sure.
0 Comments
Accepted Answer
Paulo Silva
on 15 Mar 2011
Edited: Sam Chak
on 1 Aug 2024
Almost forgot to post my code, hope this helps those who bother to search a little bit
clf
hold on
N=100000;
x=zeros(1,N);y=x;
for a=2:N
c=randi([0 2]);
switch c
case 0
x(a)=0.5*x(a-1);
y(a)=0.5*y(a-1);
case 1
x(a)=0.5*x(a-1)+.25;
y(a)=0.5*y(a-1)+sqrt(3)/4;
case 2
x(a)=0.5*x(a-1)+.5;
y(a)=0.5*y(a-1);
end
end
plot(x,y,'.')
title('Sierpinski’s triangle made in matlab by Paulo Silva')
legend(sprintf('N=%d Iterations',N))
More Answers (2)
Paulo Silva
on 28 Feb 2011
That is homework and nothing seems to be done so far, I won't do your homework, I will just point out what you need to know (with preallocation of memory)
- Learn how to create vectors with 1xN zeros
doc zeros
- Learn how to use the for loop
doc for
- Learn how to use the function randi
doc randi
- Learn how to use the function switch
doc switch
- Learn how to use the function plot
doc plot
Now here's the logic
#1 Decide the number of iterations you want N
#2 Allocate the space in memory for two vectors (x and y coordinates) 1xN
#3 Create a loop from 2 to N
#4 Do the randi to select what rule to apply
#5 Use the switch statement to use only one of the 3 rules
#6 After the loops ends do the plot of both vectors

%Free tip: %this updates the x value in each iteration (example of the first rule) x(IterationNumber)=0.5*x(IterationNumber-1); %Now I ask why I said to create a loop from 2 to N? can you find that out %in this tip?
Sean de Wolski
on 28 Feb 2011
You should start by asking your teacher since either they've failed at teaching you or you've failed at learning.
See Also
Categories
Find more on Annotations in Help Center and File Exchange
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!