代码21-10.txt

matlab宝典书籍源码 全部代码 希望大家喜欢

N=50000;                       % 初始值
v=rand(N,1);                   % 随机数
ABC=[-1, sqrt(2)*i, 1];
xy=[i; zeros(N-1,1)];
for kk=2:N
if v(kk)<1/3, xy(kk)=0.5*xy(kk-1)+0.5*ABC(1);
elseif v(kk)<2/3, xy(kk)=0.5*xy(kk-1)+0.5*ABC(2);
else xy(kk)=0.5*xy(kk-1)+0.5*ABC(3);
end
end
plot(xy, '.' ,'MarkerSize',4)  % 画图
axis equal                     % 设置显示比例