function [I,Q,R,omega_t,theta_t] = simulation_sans(dist,nb_ondes)
%La distance est en mètre

%Constantes
lambda=2.9*10^8/(900000000);
Beta=0;
C=10*10^-6;

%On avance à 10m/s, On échantillonne à 1kHz
n=1000*dist/10;

%On crée les vecteurs, remplis avec des zéros
%Le temps en X, le numéro de l'onde en Y
phi_t=zeros(n,nb_ondes);
alpha_t=zeros(n,nb_ondes);
omega_t=zeros(n,nb_ondes);
theta_t=zeros(n,nb_ondes);

%On génère maitenant les deux vecteurs phi_t et alpha_t complets 
for i = 1:nb_ondes,
    k=1;
    while k~=n,
        phi=rand(1,1)*2*pi;
        alpha=rand(1,1)*2*pi;
        k=k+1;
        phi_t(k,i)=phi;
        alpha_t(k,i)=alpha;
    end
end


%On construit omega_t et theta_t
omega_t=(2*pi*10/lambda)*cos(alpha_t);
theta_t=phi_t;

%On calcule maintenant I et Q
for i = 1:n,
   I(i)=0;
   Q(i)=0;
	for j = 1:nb_ondes,
   	  I(i)=I(i)+C*cos(omega_t(i,j)*i/1000+theta_t(i,j));
      Q(i)=Q(i)+C*sin(omega_t(i,j)*i/1000+theta_t(i,j));
   end
   R(i)=sqrt((I(i))^2+(Q(i))^2);
end

