% Calcul de la propagation d'une onde electromagnetique
% dans un milieu peu dense.
%
% Licence de Physique, UL5A
% Universite Paul Sabatier
% Licence de Physique UL5, TD numeriques, 2003-2004
% Version 08/02/2004 LK
% -----------------------------------------------

clear all; 

% ---------
initiales=' ??'

c0 = 1

w1 = ??
w1_2 = w1^2;

W = ??
W_2 = W^2;

% Initialisation courbe plus (a completer)
% -----------------------------------------

wi =  W+0.001;            % w initial pour le graphe   (de W a 10)
wf = ??;                  % w final pour le graphe
Nw = 1000;                % Nombre de points   
dw = (wf-wi)/(Nw-1);      % pour avancer ds la propagation
vmax = 5;                 % Vitesse de phase maxi

% Etude de k en fonction de w (a completer)
% -----------------------------------------

for lw=1:Nw
  w_p(lw) = wi + (lw-1) * dw;
  w_2     = w_p(lw)^2;

  vph_p(lw) = ??;                    % Vitesse de phase : c
  k_p (lw)  = ??;                    % Nb d'onde :   k = w/c
  vg_p (lw) = ??;      % Vitesse de groupe

  if (vph_p(lw) > vmax)
    vph_p(lw) = vmax;
  end
end
kmax = k_p(Nw);

% Initialisation courbe moins (a completer)
% -----------------------------------------
wi = ??;                       % w initial
wf = w1-0.001;                % w final
dw = (wf-wi)/(Nw-1);          

for lw=1:Nw
   w_m(lw)=wi+(lw-1)*dw;
   w_2=w_m(lw)^2;

   vph_m(lw) = ??;                   % Vitesse de phase
   k_m (lw)  = ??;                   % Nb d'onde :   k = w/c
   vg_m (lw) = ??;                   % Vitesse de groupe

   if (vph_m(lw) > vmax)
       vph_m(lw) = vmax;
    end
   if (k_m(lw) > kmax)
      k_m(lw) = kmax;
   end
end

% Partie graphique
% ----------------
plot(k_p,w_p,k_m,w_m)
xlabel('k')
ylabel('w')
title(['Relation de dispersion', initiales])
pause

plot(w_p,vg_p,w_m,vg_m)
xlabel('w')
ylabel('Vitesse de groupe')
title(['Vitesse de groupe', initiales])
pause

plot(w_p,vph_p,w_m,vph_m)
xlabel('w')
ylabel('Vitesse de phase')
title(['Relation de phase', initiales])


