lightscatpro5.htm

球形粒子Mie 散射用Matlab编写的Mie theory计算程序。在光子学中运用广泛

- T2spheresG.m : non valid&eacute; <br>
- TNspheres.m : calcule la T-matrice dans le cas de la diffusion par N sph&egrave;res</p>
<p>(Matrices de translation)<br>
  - translation_z.m : calcule la matrice de translation suivant l'axe z des fonctions sph&eacute;riques dans le cas m=1 (cas ou thata_i=0) <br>
  - translation_zc.m : calcule la matrice de translation suivant l'axe z des fonctions sph&eacute;riques dans le cas g&eacute;n&eacute;ral (m qcq) (utilis&eacute; par T2spheres.m)<br>
  - translation_fct.m : calcule la matrice de translation suivant l'axe z sur une distance de n*d &agrave; partir de la matrice de translation suivant z sur une distance d<br>
  - translation_qcq.m : calcule la matrice de translation suivant un vecteur quelconque (r,theta,phi)=(d,eta,psi)<br>
  - wigner_mat1.m : cacule de mani&egrave;re matricielle les 3J-symboles de Wigner initiaux <br>
  permettant d'initier une r&eacute;currence (utilis&eacute; par translation_z.m)<br>
  - wigner_recurrence1. m : permet le cacul des 3J-symboles de Wigner par r&eacute;currence (utilis&eacute; par translation_z.m)<br>
  - Wigner_mat2.m et wigner_recurrence2.m font la m&ecirc;me chose que wigner_mat1.m et wigner_recurrence1.m mais avec les termes sqrt((j3-m1+m2+1)!)Wigner(j1,j2,j3,m1,m2,m3)<br>
  au lieu de Wigner(j1,j2,j3,m1,m2,m3)</p>
<p><strong>Groupe 3	: exploitation du champ calcul&eacute;</strong> <br>
  - PostTraitement.m : reconstruit l'onde souhait&eacute;e dans l'espace souhait&eacute;<br>
- PostTraitement2.m :  exploite l'onde reconstruite<br>
- CoordEspace.m : Script o&ugrave; sont d&eacute;finis les coordonn&eacute;es de la zone spatiale &eacute;tudi&eacute;e <br>
- recons_champs.m : calcule les 6 composantes du champ dans l'espace &agrave; partir de leur coefficients de d&eacute;composition sur une base sph&eacute;rique&micro;<br>
- recons_champs_b.m : comme recons_champsB1.m mais en rempla&ccedil;ant les boucles par des calculs sur des matrices 3D (int&eacute;r&ecirc;t limit&eacute;)<br>
- view_champs.m : repr&eacute;sente les 6 composantes du champ dans l'espace<br>
- view_vecteur.m
: repr&eacute;sente sous forme de fl&egrave;ches dans l'espace un champ de vecteur <br>
- poynting.m : calcule et repr&eacute;sente le vecteur de poynting dans l'espace<br>
- stokes.m : Calcules les &eacute;l&eacute;ments de Stokes (polarisation) dans l'espace<br>
- axespolar.m : calcule et repr&eacute;sente les axes de la polarisation de l'onde en fonction de la direction d'observation <br>
- cross_section.m : calcule les sections efficaces de diffusion, d'exctinction et de r&eacute;trodiffusion<br>
- diffusion_diagram.m : calule et repr&eacute;sente les digramme de diffusion 1D ou 2D <br>
- stokes.m 
: calcule et repr&eacute;sente les &eacute;lements de Stokes dans l'espace. </p>
<p><strong>Autres :</strong> <br>
- Debut.m : affiche un message dans une fen&ecirc;tre au d&eacute;but du lancement d'un script <br>
- Fin.m
: idem &agrave; la fin d'ex&eacute;cution d'un script <br>
- retour.m : permet de retourner dans le r&eacute;pertoire &quot;work&quot; de matlab<br>
  - Mprint.m : permet d'enregistrer des figures matlab sous forme d'images <br>
  dont on peut choisir la taille et le format<br>
  - changerep.m : permet de passer un vecteur des rep&egrave;res sph&eacute;riques&lt;-&gt;cart&eacute;siens<br>
  - affichage.m : fonction non utilis&eacute;e permettant de mani&egrave;re simplifi&eacute;e d'afficher des r&eacute;sultats sous forme graphique<br>
  - Jlimit.m : permet de simuler le comportement de la fonction de bessel sph&eacute;rique d'un ordre quelconque en z=0<br>
  - ordre.m : calcul l'ordre conseill&eacute; de d&eacute;veloppement en fonction vectorielle <br>
  sph&eacute;rique &agrave; partir des donn&eacute;es du probl&egrave;mes. Le r&eacute;sultat est <br>
affich&eacute; dans l'interface<br>
- FichierTemp.m : fichier dans lequel vont &ecirc;tre copi&eacute; les fichiers d&eacute;crivant l'agr&eacute;gat</p>
<p>Fonctions permettant de faire des tests :<br>
  -test_multi_translation_z.m <br>
-test_translation_qcq2.m<br>
-test_translation_qcq<br>
-translation_astuce.m</p>
<p><a name="struct"></a>5) Structure du programme<br>
________________________<br>
Le programme a &eacute;t&eacute; structur&eacute;e en 3 parties diff&eacute;rentes.</p>
<table width="590" border="1">
  <tr>
    <td width="504" bgcolor="#99FF00"><p><span class="Style2">1-D&eacute;marrage</span> <strong>et acquisition des donn&eacute;es</strong> <br>
      </p>
    </td>
    <td width="70"><strong>Variables</strong></td>
  </tr>
  <tr>
    <td><p>- Lancement de l'interface graphique<br>
      ou<br>
      - Editeur de commande Matlab </p>
    </td>
    <td>data</td>
  </tr>
  <tr>
    <td bgcolor="#99FF00"><strong>2-Calcul de la T-matrice et des co&eacute;fficients de d&eacute;composition des champs </strong></td>
    <td>&nbsp;</td>
  </tr>
  <tr>
    <td><p>- Cacul des coefficients de l'onde incident <br>
      - Calcul de la T-matrice
          <br>
          - Calcul des coefficients de l'onde diffus&eacute;e <br>
          (Coefficients de l'onde
dans la bile)  </p>
    </td>
    <td><p>P_i<br>
        T<br>
      P_s<br>
      (P_o)</p>
      </td>
  </tr>
  <tr>
    <td bgcolor="#99FF00"><strong>3-Post Traitement </strong></td>
    <td>&nbsp;</td>
  </tr>
  <tr>
    <td>- D&eacute;finition de la zone spatiale consid&eacute;r&eacute;e<br>
      - Reconstruction du champ dans l'espace<br>
      - Calcul des grandeurs physiques </td>
    <td>space<br>
      Champ<br>
      Result</td>
  </tr>
</table>
<p><br>
<a name="var"></a>6) Liste de variables<br>
___________________<br>
  Pour faciliter l'&eacute;change des variables entre les programmes, celles-ci ont &eacute;t&eacute; regroup&eacute;es en structure.<br>
  Par exemple 
la variable <em>lambda</em> qui contient la valeur de la longueur d'onde est dans la structure <em>data</em>. Pour appeler cette varaible, il faut &eacute;crire <em>data.lambda</em>.</p>
<p><strong>Structure &quot;data&quot;</strong><br>
Contient les variables d&eacute;crivant l'onde incidente et l'agr&eacute;gat</p>
<table width="884" border="1">
  <tr>
<td width="190"><strong>   Name </strong></td>
<td width="532"><strong> description </strong></td>
<td width="140"><strong> format </strong></td>
  </tr>
  <tr>
<td>  lambda </td><td> wavelength in free space, in nm </td><td> number </td>
  </tr>
  <tr>
<td>  R </td><td> Radius of spheres </td><td> vector </td>
  </tr>
  <tr>
<td>  n1 </td><td> real index of the environment </td><td> number </td>
  </tr>
  <tr>
<td>  n2 </td><td> real index of spheres </td><td> vector </td>
  </tr>
  <tr>
<td>  N </td><td> complex index of spheres </td><td> vector </td>
  </tr>
  <tr>
<td>  k </td><td>  2\pi/\lambda  </td><td> number </td>
  </tr>
  <tr>
<td>  ConducInf </td><td> perfectly conductive sphere if 1 </td><td> vector </td>
  </tr>
  <tr>
<td>  sigma </td><td> conductivity of spheres </td><td> vector </td>
  </tr>
  <tr>
<td>  Nbspheres </td><td> number of spheres in the aggregate </td><td> integer </td>
  </tr>
  <tr>
<td>  d </td><td> distance between spheres </td><td> number </td>
  </tr>
  <tr>
<td></td><td> for linear regular spheres on z-axis </td><td> </td>
 </tr>
  <tr>
<td>position </td><td> radius position of each spheres </td><td> vector </td>
 </tr>
  <tr>
<td>eta </td><td> \theta position of each spheres </td><td> vector </td>
 </tr>
  <tr>
<td>psi </td><td> \varphi position of each spheres </td><td> vector</td>
 </tr>
  <tr>
<td>angledincidence </td><td> incidence angle of the propagation vector </td><td> number</td>
 </tr>
  <tr>
<td>polarisation </td><td> (1\rightarrow 6) state of polarization </td><td> integer</td>
 </tr>
  <tr>
<td>Lmin </td><td> minimum order (1 in general) </td><td> integer</td>
 </tr>
  <tr>
<td>Lmax </td><td> expansion order 0\leq l \leq \mathrm{Lmax} </td><td> integer</td>
 </tr>
  <tr>
</table>

<p><strong>Structure for the T-matrix</strong><br>
La T-matrice est d&eacute;crite par la structure <strong>T </strong>qui contient 16 variables pouvant chacune &ecirc;tre scalaire ou matricielles</p>
<table width="884" border="1">