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 

皮尔逊三型曲线