计算物理基础（彭芳麟）

第二章课后习题3、习题4

习题3 科赫曲线

clc
clear all
close all

new = [0,1i];
subplot(2,2,1);
plot(new)
axis equal
for k = 1:3
    old = new;
    n = length(old)-1;  %线段总数比节点数少1
    diff = (old(2:n+1)-old(1:n))/3;  %将每条线段三等分
    new(1:4:4*n-3) = old(1:n);  %起点坐标
    new(2:4:4*n-2) = old(1:n)+diff;  %插入点P1坐标
    new(3:4:4*n-1) = new(2:4:4*n-2)+diff*((1+sqrt(3)*1i)/2);  %插入点P2坐标
    new(4:4:4*n) = old(1:n)+2*diff;  %插入点P3坐标
    new(4*n+1) = [1i];  %终点坐标
    subplot(2,2,k+1);
    plot(new)
    axis equal
end

 习题4 Minkowski曲线

clc
clear all
close all
new=[0,1];
subplot(2,2,1);
axis([0,1,-0.5,0.5]);
axis equal
for k = 1:4
    old = new;
    n = length(old)-1;  %线段总数比节点数少1
    diff = (old(2:n+1)-old(1:n))/4;  %将每条线段四等分
    p1 = old(1:n)+diff;
    p2 = p1+diff*(i);
    p3 = p2+diff;
    p4 = p3-diff*2*(i);
    p = p3-diff*(i);
    p5 = p4+diff;
    p6 = p5+diff*(i);
    new(1:15:15*n-14) = old(1:n);
    new(2:15:15*n-13) = p1;
    new(3:15:15*n-12) = p1;
    new(4:15:15*n-11) = p2;
    new(5:15:15*n-10) = p2;
    new(6:15:15*n-9) = p3;
    new(7:15:15*n-8) = p;
    new(8:15:15*n-7) = p;
    new(9:15:15*n-6) = p4;
    new(10:15:15*n-5) = p4;
    new(11:15:15*n-4) = p5;
    new(12:15:15*n-3) = p5;
    new(13:15:15*n-2) = p6;
    new(14:15:15*n-1) = p6;
    new(15:15:15*n) =  old(2:end);
    subplot(2,2,k)
    plot(new)
    axis([0,1,-0.5,0.5]);
end