% Simulation d'un Processus de Poisson Simple 
% à partir de ses temps de sauts exponentiels.
clf
clear
lambda=.25;  % Intensité des sauts
n=30;        % Nombre de sauts
lastx=0.
lasty=0.
for i=0:n-1,
 X(1+2*i)=lastx;
 X(2+2*i)=lastx+rexpweib(1,lambda);
 Y(1+2*i)=lasty+1;
 Y(2+2*i)=Y(1+2*i);
 lasty=Y(2+2*i);
 lastx=X(2+2*i);
end
%
figure(1)
plot(X,Y,'b-')
hold on
xlabel('t')
ylabel('Nt')
legend('Trajectoire d''un Processus de Poisson Simple',4)
%
n=length(X);
Z=Y(2:n)./X(2:n);
figure(2)
clf
plot(X(2:n),Z,'b-')
hold on
plot(X(2:n),lambda*ones(n-1),'r-.')
xlabel('t')
ylabel('Nt/t')
legend('Nt/t','Intensite lambda')
