function esrefmul
% Dominic Grenier octobre 2009

%N = nombre de lignes

RgRc = [50 50]; %charges aux extrémités (1,2)
Z0   = [50 20]; %impédances caractéristiques des lignes [Ohms] (1,N)
L    = [1200 600]; %longueurs des lignes [m] (1,N)
Vp   = [100 66.7]; %vitesse sur les lignes [% de c] (1,N)
Sinit= [60 0]; %conditions initiales sur les lignes (1,N)
J    = [0 1/200 0]; %jonctions à 3 résistances en T - R_s1 G_p R_s2; J=[] si N=1 (N-1,3)
zdep = 2; %jonction de départ 0:à gauche N:à droite
sdep = [80 0]; %signaux de départ vers à la gauche et vers la droite (1,2)
type = 'i'; %type de signaux 'v':tension 'i':courant 'p':puissance
rmin = 0.01; %optionnel - rapport minimum de s/smax; défaut rmin=0.001

signaux = refmul(RgRc,[Z0' L' Vp'],J,zdep,type,sdep,rmin);
traceref(signaux,[Z0' L' Vp'],Sinit,{type 'porte' 0.4e-6},'z','vst',[(4/5)*1800]);

%diagrammes possibles de traceref
%'z' : diagramme en Z
%'vst',[z_obs K t_fin] : évolution à z=z_obs, de t=0 à t=t_fin sur K points
%'vsz',[t_obs K] : distribution à t=t_obs, de z=0 à z=somme(l_k) sur K points
%     défaut : t_fin=t_max, K=200
%{type g(t) T} : type de signaux, g(t)='échelon' 'impulsion' 'porte' de largeur T (optionnel)
%     défaut : 'échelon', T=t_max/10
