Advertisement
Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- function E = E_field_LP_Mode(l,gamma_c,gamma_g, r_c, theta,rho)
- % Cette fonction calcule la distribution du champ électrique d'un mode
- % guidé LP_lm.
- % Variables d'entrée :
- % l : valeur du paramèttre l (L) du mode;
- % gamma_c et gamma_g : valeur de ces paramètres définit en cours pour
- % un mode donné
- % r_c : rayon du coeur
- % theta,rho : coordonnées polaires des points où le champ doit être
- % calculé.
- %-------------------------------------------------------------------------
- f = zeros(size(rho)); %il est essentiel de prédéfinir les dimensions de f
- Rho1 = rho <= r_c; % Dans coeur, fonction Bessel_J
- Rho2 = rho > r_c; % Dans gaine, fonction Bessel_K
- g = cos(l*theta);
- B = besselj(l,gammma_c*r_c)/besselk(l,gamma_g*r_c);
- f(Rho1) = besselj(l, gamma_c*rho(Rho1)); %calcule à l’intérieur du coeur
- f(Rho2) = B*besselk(l,gamma_g*rho(Rho2)); %calcule à l’extérieur du cœur (i.e. dans gaine)
- E = g.*f;
- end
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement