Sierpinski-垫片(matalb)

functionsierpinski_ifs(n,w1,w2,w3)%sierpinski_ifs(10000,1/3,1/3,1/3)%w1,w2,w3出现频率n=10000000;w1=1/3;w2=1/3;w3=1/3;M1=[0.50000.50];M2=[0.500.500.50];M3=[0.500.2500.50.5];x=0;y=0;%r为[0,1]区间内产生的n维随机数组r=rand(1,n);B=zeros(2,n);k=1;%当0r(i)1/3时,进行M1对应的压缩映射;%当1/3=r(i)2/3时,进行M2对应的压缩映射;%当2/3=r(i)1时,进行M3对应的压缩映射;fori=1:nifr(i)w1a=M1(1);b=M1(2);e=M1(3);c=M1(4);d=M1(5);f=M1(6);elseifr(i)w1+w2a=M2(1);b=M2(2);e=M2(3);c=M2(4);d=M2(5);f=M2(6);elseifr(i)w1+w2+w3a=M3(1);b=M3(2);e=M3(3);c=M3(4);d=M3(5);f=M3(6);endendendx=a*x+b*y+e;y=c*x+d*y+f;B(1,k)=x;B(2,k)=y;k=k+1;endplot(B(1,:),B(2,:),'.','markersize',0.1)%title({['kappa=',num2str(kappa)],['起始点为',strcat('z1=[',num2str(z1),']')],['起始点为',strcat('z2=[',num2str(z2),']')]});xlabel('实轴');ylabel('虚轴');gtext('Sierpinski垫片');%将标注函数放在相应曲线上