forked from csdms-contrib/slepian_alpha
-
Notifications
You must be signed in to change notification settings - Fork 0
/
blob.m
54 lines (50 loc) · 1.04 KB
/
blob.m
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
function varargout=blob(N,Nj)
% BLOB(N,nj)
% [x,y]=BLOB(N,nj)
%
% Makes (moving) picture of a random blob by superposition of random
% circles.
%
% INPUT:
%
% N Number of loops for, and if movie [default: 100]
% Nj Smoothness, roughly
%
% SEE ALSO:
%
% RANDCIRC
%
% Last modified by fjsimons-at-alum.mit.edu, 10/06/2008
defval('N',100)
defval('Nj',10);
[xold,yold]=randcirc(0,0,1,1,10);
r=linspace(0,1,Nj+1);
rm=1-r;
if nargout==0
figure(gcf)
end
for index = 1:N
[x,y]=randcirc(0,0,1,0.2,10);
col=rand(1,3);
for j=1:Nj,
% First is increasing radius, second decreasing
xx=x*r(j)+xold*rm(j);
yy=y*r(j)+yold*rm(j);
if nargout==0
% plot(xx,yy,'-','linewidth',2,'color',[1,1,1]*0);
% plot(xx,yy,'-','linewidth',2,'color',col);
axis([-2 2 -2 2])
axis off
drawnow
% pause(0.5)
fill(xx',yy',col)
axis equal ; axis off
end
end
xold=x;
yold=y;
end
varn={'xx','yy'};
for index=1:nargout
varargout{index}=eval(varn{index});
end