poission.m

来自「自编遗传算法程序」· M 代码 · 共 23 行

clc
lambda=[2 2 3 1 0.8 0.5 1];

M=size(lambda);

p=ones(M);
n=zeros(M);

while(1)
    u=rand(M);
    p=p.*u;

    if all(p<exp(-lambda)),break;end
    n(p>=exp(-lambda))=n(p>=exp(-lambda))+1;
    
    %if p<exp(-lambda), break; end
    %n=n+1;
    
    %for i=1:length(lambda(:))
    %    if p(i)<exp(-lambda(i)),   end
    %end
end