gasket.m

这是机械工业张德丰《MATLAB与外部程序接口编程》随书源码

% gasket.m
function theImage=gasket(numPoints)
% GASKET An image of a Sierpinski Gasket.
% IM=GASKET(NUMPOINTS)
%
% Example:
% x=gasket(50000);
% imagesc(x)=;colormap([1 1 1;0 0 0]);
%axis equal tight
theImage=zeros(1000,1000);
corners=[866 1:1 500;866 1000];
startPoint=[866 1];
theRand=rand(numPoints,1);
theRand=ceil(theRand*3);
for i=1:numPoints
    startPoint =floor((corners(theRand(i),:)+startPoint)/2);
    theImage(startPoint(1),startPoint(2))=1;
end