randpos.m

来自「This program shows the distributions of 」· M 代码 · 共 19 行

function [z] = RandPOS(n)
%n=1000;
x = 1.5*rand(1,n)-1;
y = sqrt(3)*rand(1,n)-0.5*sqrt(3);

for k=1:n
    if x(k)>-1 & x(k)<-0.5 & y(k)>(1+x(k))*sqrt(3)
        x(k)=x(k)+1.5;
        y(k)=y(k)-sqrt(3)/2;
    elseif x(k)>-1 & x(k)<-0.5 & y(k)<-(1+x(k))*sqrt(3)
        x(k)=x(k)+1.5;
        y(k)=y(k)+sqrt(3)/2;
    end
    if sqrt(x(k)^2+y(k)^2)<0.1
        x(k)=0.1;
        y(k)=0;
    end
end
z=x+j*y;