10.08.2013 Views

Etude théorique de radars géologiques - Epublications - Université ...

Etude théorique de radars géologiques - Epublications - Université ...

Etude théorique de radars géologiques - Epublications - Université ...

SHOW MORE
SHOW LESS

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

Samuel BESSE : Étu<strong>de</strong> Théorique <strong>de</strong> Radars Géologiques ­ Analyses <strong>de</strong> sols, d'antennes et interprétation <strong>de</strong>s signaux 41<br />

Phase d'initialisation 1 ère itération : étape 1 1 ère itération : étape 2 1 ère itération : étape 3<br />

Figure 33 : Représentation <strong>de</strong> l'algorithme <strong>de</strong> la métho<strong>de</strong> dite Diamond­Square4D sur un motif élémentaire.<br />

L'application d'une itération génère huit cubes sur lesquels on peut répéter l'algorithme.<br />

Comme dans le cas <strong>de</strong>s surfaces la valeur moyenne quadratique est liée à la valeur du<br />

paramètre h et la répartition <strong>de</strong>s pondérations suit une loi gaussienne. L'expression (20) se déduit<br />

d'une étu<strong>de</strong> statistique sur un grand nombre <strong>de</strong> réalisations (une centaine). L'étu<strong>de</strong> statistique révèle<br />

également que h est indépendante du nombre <strong>de</strong> motifs (NX.NY.NZ) et converge pour un<br />

nombre d'itérations Nit supérieur ou égal à quatre.<br />

h=0,31 0,45 ⋅h−1,3 ⋅h 2 1,4 ⋅h 3<br />

Il est souvent important <strong>de</strong> connaître la pondération maximale possible notamment dans les<br />

applications <strong>de</strong> la troisième partie. Par exemple, si l'on utilise la matrice pour modéliser un sol<br />

hétérogène en permittivité avec la métho<strong>de</strong> FDTD, il faut être certain <strong>de</strong> mailler suffisamment<br />

finement pour limiter la dispersion numérique (la taille d'une maille doit être inférieure à λ milieu/10).<br />

Pour cela, l'abaque présenté figure 34 permet <strong>de</strong> déterminer la valeur <strong>de</strong> pondération maximale<br />

<strong>théorique</strong> présente dans le matériau. Il est obtenu par un calcul numérique en substituant toutes les<br />

variables aléatoires intervenant dans l'algorithme par la borne supérieure <strong>de</strong> l'intervalle dans lequel<br />

elles sont choisies. Dans la pratique, le majorant <strong>théorique</strong> a une probabilité infime d'être atteint et<br />

cette probabilité est d'autant plus petite que le nombre d'itérations est important.<br />

Le passage <strong>de</strong> la matrice normalisée M0 à la matrice <strong>de</strong> permittivité s'obtient par (21) où les<br />

paramètres 〈 r 〉 et k laissés au choix <strong>de</strong> l'utilisateur désignent respectivement la permittivité<br />

moyenne et un facteur multiplicatif.<br />

r i , j , k = 〈 r 〉k⋅M 0 i , j , k (21)<br />

(20)

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!