Back to Gallery
  • Remix
  • Share
  • New Entry
on 22 Nov 2023
  • 5
  • 3
  • 0
  • 0
  • 1297
drawframe(1);
Write your drawframe function below
function drawframe(f1)
[x,y] = meshgrid(-5:5,-5:5);
z = peaks(11)/1.2;
tri = delaunay(x,y);
tr = triangulation(tri,x(:),y(:),z(:));
plus_minus=1;
index=round(linspace(0,2,48));
zz=linspace(1,1.2,48);
i=index(f1);
TR_out=Fractal_Kochsnow3D_tri(tr,1,i,plus_minus);
c=jet(size(TR_out.Points,1));
% trisurf(TR_out,'FaceColor','blue')
patch('Faces',TR_out.ConnectivityList,'Vertices',TR_out.Points,'FaceColor','red','EdgeColor','y','LineWidth',0.5);
axis equal off;
view(-30+f1,30);
zoom(zz(f1))
function TR_out = Fractal_Kochsnow3D_tri(TR,m,n,plus_minus,varargin)
%% 缺省值处理。。 默认为向外扩展，可选为向内或者向外
if isempty(varargin)
inorout = 'out';
else
inorout = varargin{1};
end
%% 循环获取结果
%inorout1={'out';'in'};
if n==0
TR_out = TR;
else
for i = 1:n
P=[];
T=[];
for j = 1:size(TR.ConnectivityList,1)
TR1 = triangulation([1,2,3],TR.Points(TR.ConnectivityList(j,:),:));
Temp = KochSeed(TR1,m,plus_minus,inorout);
if j==1
P=[P;Temp.Points];
T=[T;Temp.ConnectivityList];
else
T=[T;Temp.ConnectivityList+size(P,1)];
P=[P;Temp.Points];
end
end
TR_out = triangulation(T,P);
TR=TR_out;
end
end
end
function TR_out=KochSeed(TR,m,plus_minus,varargin)
%% 缺省值处理。。 默认为向外扩展，可选为向内或者向外
if isempty(varargin)
inorout = 'out';
else
inorout = varargin{1};
end
% if strcmpi(inorout,'in')
% inorout = 1;
% else
% inorout = -1;
% end
%% 处理不同的inorout
if strcmpi(inorout,'in')
v=TR.Points([1,3,2],:);
TR = triangulation(TR.ConnectivityList,v);
end
%正三角形的边长dl
dl=norm(TR.Points(1,:)-TR.Points(2,:));
%以dl/2构建正四面体,那么它的高h
h=dl/2*sqrt(6)/3*m;
v=[TR.Points;[TR.Points+TR.Points([2,3,1],:)]/2;TR.incenter+plus_minus*TR.faceNormal*h];
f=[[1,4,6;2,5,4;3,6,5];[4,5,7;5,6,7;6,4,7]];
TR_out=triangulation(f,v);
end
end
Animation
Remix Tree