function peano_curve(n)
peano_old=[0,0;0,1;0.5,1;0.5,0;1,0;1,1];
x=peano_old(:,1);
y=peano_old(:,2);
figure
set(gcf,'position',[0,0,1920,1080]);
plot(x,y)
axis equal
for i=1:n-1
p1=[peano_old(:,1),2+1/(3^i-1)-peano_old(:,2)];
p1=p1(length(p1):-1:1,:);
p2=[p1(:,1),4+3/(3^i-1)-p1(:,2)];
p2=p2(length(p2):-1:1,:);
peano_new=[peano_old;p1;p2];
p1=[2+1/(3^i-1)-peano_new(:,1),peano_new(:,2)];
p1=p1(length(p1):-1:1,:);
p2=[4+3/(3^i-1)-p1(:,1),p1(:,2)];
p2=p2(length(p2):-1:1,:);
peano_new=[peano_new;p1;p2];
peano_old=peano_new/(3+2/(3^i-1));
x=peano_old(:,1);
y=peano_old(:,2);
plot(x,y)
axis equal
pause(1)
end
end
皮尔逊三型曲线