Signal observations - Revision history

Uadmin: /* modélisation du « signal » */

2007-04-14T01:22:00Z

‎modélisation du « signal »

Uadmin: /* modélisation du « signal » */

2006-01-18T13:18:01Z

‎modélisation du « signal »

Uadmin: /* Critique de la méthodologie et des paramètres choisis */

2006-01-18T13:16:18Z

‎Critique de la méthodologie et des paramètres choisis

Uadmin: /* Critique de la méthodologie et des paramètres choisis */

2006-01-18T13:16:11Z

‎Critique de la méthodologie et des paramètres choisis

Uadmin at 13:15, 18 janvier 2006

2006-01-18T13:15:31Z

Uadmin at 13:15, 18 janvier 2006

2006-01-18T13:15:24Z

Uadmin at 13:15, 18 janvier 2006

2006-01-18T13:15:16Z

Uadmin at 13:15, 18 janvier 2006

2006-01-18T13:15:09Z

Nouvelle page

=modélisation du « signal »=

Nous essayons ici de reconstruire la courbe du signal créé par les observation sur plusieurs dizaines d'années. En s'appuyant sur des harmoniques, il est possible de dégager des fréquences privilégiées.

Les graphiques présentés ci-après et utilisés pour cette étude sont la propriété de [[Larry Hatch]] (larryhatch(at)sysmatrix.net) que je remercie. Si vous êtes intéressés, vous pouvez trouver plus d'informations sur son site web:
{{Ref:Web | Titre=*U* UFO Database | Sujet= Un catalogue de données de 18370 objets, représentant plus de 20 années de recherche | URL=http://www.larryhatch.net | Date=02/2005 pour cette étude}}

[[UFO Database|A propos de la consistance de cette base de données]]

<u>Autres sources de données</u>

Nous prévoyons d'étendre le champ de cette étude à d'autres sources de données, notamment à tous les catalogues librement accessibles, par exemple celui de [[Luc Chastan]].
En s'appuyant sur celui-ci nous proposons un modèle partagé d'édition et de vérification des ressources et des connaissances. Ce qui correspond aussi à la mise au point aussi d'un modèle de saisie s'appuyant sur la base de Luc Chastan et accessible ici : [[UFO:Cases:2A]]. Nous souhaitons mettre, évidemment, l'accent sur ''l'environnement'' proche pour chaque cas étudié (sources d'effets néguentropiques pour le biotope à proximité ou dans la direction empruntée par l'engin).

Qui savait par exemple qu'il y avait une [[Radar HF de Valensole|station radar HF]] non loin de valensole à l'époque ou feu [[Maurice Masse]] avait été le témoin d'un atterissage ? Connaissant l'impact possible sur l'environnement/la biosphère que peut avoir un radar HF, le témoignage d'un être penché sur un pied de lavande, peut prendre une perspective radicalement nouvelle.

<u>Méthodologie d'analyse et d'évaluation de l'information</u>

<b>[[Information-Coherence-Substantialite|Sur l'information sa cohérence, et la substantialité.]]</b>

====Observations, sur 50 ans====

{{Ref:Database| Titre=Nombre d'observations mensuelles| URL=http://www.larryhatch.net/50YEARS.html | Date=01/1946 au 01/1996}}
[[image:50YEARSG.gif]]<BR>
[http://www.larryhatch.net/50YEARS.html 50CLIP.gif © 1985-2005 Larry Hatch]
<i>(1 pixel = 1 mois)</i>

====Spectre fréquentiel des observations sur 50 ans====

Un graphique réalisé par Jacques Vallée au cours des années 1970 propose visualiser les observations sous formes de fréquences. Le graphique ci-après est également issu du site de Larry Hatch.

{{Ref:Database| Titre=Observations en fréquence| URL=http://www.larryhatch.net/FREQSPC.html | Date=01/1946 au 01/1996}}
[[image:FREQSPC.gif]]<BR>
[http://www.larryhatch.net/FREQSPC.html FREQSPC.gif © 1985-2005 Larry Hatch]

Voici la légende qui accompagne ce graphique sur le site de Larry Hatch :

<i>"Every UFO sighting in these data is compared to all subsequent listings,
regardless of location. A special array is incremented for each interval
in weeks elapsed between sightings.

Counts for each interval show as white dots. The violet curve is a running
average of these counts. As expected, time intervals between major waves show
as peaks in this curve, and are marked with a green carat ^ sign. [...]"</i>

Ce graphique montre les vagues les plus caractéristiques avec leurs fréquences propres. La dispersion des observations liée à chaque vague.

====Modélisation des fréquences remarquables====

Le graphique en fréquence a été repris afin d'essayer pour extraire un ou plusieurs modèles basés sur une description "fonctionnelle". Notons immédiatement qu'il n'y a pas qu'un unique modèle qui saura convenir. Néanmoins et pour restreindre l'espace de recherche, nous essayons de privilégier les modèles les plus "économiques" possibles en termes de paramètres et afin de modéliser le signal. Un principe de parcimonie que l'on retrouve dans les réseaux neuronaux et qui peut aussi se perçevoir en termes de "compression d'information"... (la structure sous-jacente contenue dans l'information compressée plus le programme permettant de la lire, mais je m'égare...).

Les dispersions statistiques sont normalement caractérisées par des <b>[[loi normale|lois normales]]</b> (ou [http://fr.wikipedia.org/wiki/Loi_normale cf wikipedia]). Cependant, dans ce cas, nous observons que les enveloppes semblent difficiles à approximer, notamment pour les "pics d'observation" qui nécessitent un nombre d'individu considérable (n). Ce qui ne semble pas très réaliste aussi, toujours dans le sens d'une certaine "parcimonie", nous avons préféré employer des <b>fonctions linéaires par parties</b>.

Nous suppposons qu'un certain niveau de « bruit » (créé par les canulars, erreurs, etc.) est supportable: la stratégie est de travailler sur un grand nombre de cas (plusieurs milliers) et non seulement quelques uns: ces cas gênants devraient alors être répartis plus ou moins uniformément.

Ci-après une fonction d'approximation linéaire par partie pour l'ensemble des fréquences:

[[image:FREQSPC_lin.gif]]

[[Image:File-excel.gif]] http://www.u-sphere.com/data/FREQXLS.zip

<u>Regroupement des fréquences par classes</u>

<i>Ce regroupement est subjectif dans la mesure où les critères choisis ne s'appuient que sur une ou deux caractéristiques partagées. Ces critères sont indiqués. (Globalement, les regroupements en classes se font par ordre d'intensité (n)).</i>

<b>A chaque classe correspond une "typologie" de phénomènes, s'organisant dans la durée et en intensité de façon précise.</b>

Ci après, nous écrivons:

[[image:mu_.gif]] position (x en semaines) à laquelle le pic est maximal,
[[image:sigma.gif]] ratio correspondant à la vitesse de croissance/décroissance (pente) de l'événénement,
<i>n</i> amplitude maximale,
y = MAX(-[[image:sigma.gif]]*|x-[[image:mu_.gif]]|+<i>n</i>,0)
CC : lorsqu'indiqué, le coefficient de corrélation calculé pour la portion de courbe.

<i>Les pics sont représentés par des lettres en allant de gauche à droite</i>

=====<div style="background: #ffff00;">[[image:SJaune.gif]] Classe I : cycles d'événements importants mais de courte durée</div>=====

60 < n < 100
[[image:sigma.gif]] = 1/23*x (0,044*x)

<table border=0>
<tr><td>A:  </td><td width=200>[[image:mu_.gif]] = 115 </td><td><i>n</i> = 85</td></tr>
<tr><td>B:  </td><td width=200>[[image:mu_.gif]] = 159 </td><td><i>n</i> = 60</td></tr>
<tr><td>C:  </td><td width=200>[[image:mu_.gif]] = 264 </td><td><i>n</i> = 100</td></tr>
<tr><td>D:  </td><td width=200>[[image:mu_.gif]] = 378 (A + C)</td><td><i>n</i> = 100</td></tr>
<tr><td>E:  </td><td width=200>[[image:mu_.gif]] = 537 (B + D)</td><td><i>n</i> = 80</td></tr>
</table>
<br>
=====<div style="background: #ff9900;">[[image:SOrange.gif]] Classe I.1 : cycles d'événements associés à I, plus persistants dans la durée</div>=====

40 < n < 60
[[image:mu_.gif]] associé à la classe I
[[image:sigma.gif]]+ = 2,5
(vérifier si <i>n</i> est également associé à la classe I)

<table border=0>
<tr><td>A:  </td><td width=200>[[image:mu_.gif]] = 117 (I.A+2)</td><td width=200>[[image:sigma.gif]]- = 1,3</td><td><i>n</i> = 63</td></tr>
<tr><td>B:  </td><td width=200>[[image:mu_.gif]] = 266 (I.B+2)</td><td width=200>[[image:sigma.gif]]- = 2,8</td><td><i>n</i> = 66</td></tr>
<tr><td>C:  </td><td width=200>[[image:mu_.gif]] = 539 (I.E+2)</td><td width=200>[[image:sigma.gif]]- = 2,5</td><td><i>n</i> = 40</td></tr>
</table>

<i><u>remarque 1</u>: cette classe n'est peut-être qu'un artefact lié au choix de l'utilisation de fonctions linéaires par partie et ce au lieu de lois normales (gaussiennes). En effet, les événements rattachés sont temporellement liés à ceux de la classe 1 et correspondent à des durées plus longues: la base est plus "évasée"). Cependant, trois éléments militent contre:
* [[image:mu_.gif]] ne dépend pas de x,
* la forme globale des courbes (1 + 1.1) serait convexe non concave,
* la répartition des points apparait "strictement" linéaire là où elle est visible. </i><br>
<i><u>remarque 2</u>: il est tout à fait possible qu'il y ait des événements masqués (car confondus) avec 1.B et 1.D</i>

=====<div style="background: #ffAAAA;">[[image:SRose.gif]] Classe II</div>=====

18 < n < 50 ; (fonction linéairement croissante de x, apparement)
[[image:sigma.gif]]+ = 1/200*x (0,0052*x)
[[image:sigma.gif]]- = 1/50*x (0,0208*x)

<table border=0>
<tr><td>A:  </td><td width=200>[[image:mu_.gif]] = 48</td><td><i>n</i> = 18</td></tr>
<tr><td>B:  </td><td width=200>[[image:mu_.gif]] = 140</td><td><i>n</i> = 35</td></tr>
<tr><td>C:  </td><td width=200>[[image:mu_.gif]] = 248</td><td><i>n</i> = 40</td></tr>
<tr><td>D:  </td><td width=200>[[image:mu_.gif]] = 597</td><td><i>n</i> = 50</td></tr>
</table>

=====<div style="background: #99ff99;">[[image:SVert.gif]] Classe III : pics de faible intensité, durée des événements sensible à la fréquence</div>=====

<i>n</i> = 27
[[image:sigma.gif]]+ = 9/10 (indépendant de x)
[[image:sigma.gif]]- = 6/5 (indépendant de x)

<table border=0>
<tr><td>A:  </td><td width=200>[[image:mu_.gif]] = 209</td><td> </td></tr>
<tr><td>B:  </td><td width=200>[[image:mu_.gif]] = 280</td><td> </td></tr>
<tr><td>C:  </td><td width=200>[[image:mu_.gif]] = 365</td><td> </td></tr>
<tr><td>D:  </td><td width=200>[[image:mu_.gif]] = 390</td><td> </td></tr>
<tr><td>E:  </td><td width=200>[[image:mu_.gif]] = 422</td><td> </td></tr>
<tr><td>F:  </td><td width=200>[[image:mu_.gif]] = 473</td><td> </td></tr>
<tr><td>G:  </td><td width=200>[[image:mu_.gif]] = 579</td><td><i>n</i> = 39</td></tr>
</table>

<i>le dernier événement (G) n'appartient peut-être pas à cette classe. </i>

=====<div style="background: #9999ff;">[[image:SBleu.gif]] Classe IV.1 : Evénements suspects 1</div>=====

<i>Deux événements qui ne mériteraient pas la création d'une classe à eux seuls...Peut-être est-il possible de les rattacher à la classe I.1 si ce n'est qu'ils ne sont pas en relation avec un événement de I</i>

30 < <i>n</i> < 40
[[image:sigma.gif]]+ = 2 (indépendant de x)
[[image:sigma.gif]]- = 2,8 (indépendant de x)

<table border=0>
<tr><td>A:  </td><td width=200>[[image:mu_.gif]] = 31</td><td> </td></tr>
<tr><td>B:  </td><td width=200>[[image:mu_.gif]] = 39</td><td> </td></tr>
</table>

=====<div style="background: #ffffff;">[[image:SGris.gif]] Classe IV.2 : événements suspects 2 (en est-ce vraiment ?)</div>=====

Nous regroupons ici les cas à faible niveau d'intensité, c'est à dire, proche du niveau de bruit de fond (événements ponctuels) évalué à n~8.

<i>n</i> < 30

<table border=0>
<tr><td>A:  </td><td width=200>[[image:mu_.gif]] = 307</td><td width=200>[[image:sigma.gif]]+ = 0,3</td><td width=50>[[image:sigma.gif]]- = 2</td><td><i>n</i> = 20</td></tr>
<tr><td>B:  </td><td width=200>[[image:mu_.gif]] = 335</td><td width=200>[[image:sigma.gif]]+ = 0,5</td><td width=50>[[image:sigma.gif]]- = 2</td><td><i>n</i> = 20</td></tr>
<tr><td>C:  </td><td width=200>[[image:mu_.gif]] = 522</td><td width=200>[[image:sigma.gif]]+ = 0,9</td><td width=50>[[image:sigma.gif]]- = 0,5</td><td><i>n</i> = 27</td></tr>
</table>

<u>Bruit de fond</u> :
<table>
<tr><td>X:  </td><td width=200>[[image:mu_.gif]] = 0</td><td width=200>[[image:sigma.gif]] = 0</td><td><i>n</i> = 8</td></tr>
</table>

==Reconstruction du signal==

A l’instar d’une décomposition en séries de fourrier, ayant modélisé les « harmoniques » principales nous pouvons reconstruire le signal. Ci-après, la présentation n'est plus en fréquence mais en durée:

[[image:50YEARSG_Simulation_params.gif]]<br><br>
[[image:50YEARSG_Simulation.gif]]

- Axe des x: mois écoulés partant du 01/1946.
- Courbe violette: la courbe des observations telles qu’issues de la BD de Larry Hatch
- Courbe verte: la courbe « corrélée ».
- Coefficient de Pearson : ~0.76. Ce qui n'est pas très bon (!),
et ce probablement compte tenu du faible nombre de variables pris en compte.

<i><u>remarque 1</u>:ici ce "pré-travail" ne correspond pas directement à ce qui a été présenté plus haut en termes de paramètres mesurés et dans la mesure où seuls ici 4 paramètres principaux sont retenus.
<br><span style="background:#ffff00">Un travail exhaustif reprenant l'ensemble des fonctions linéaires approximées précédemment sera réalisé dès que je pourrais</span>.

<u>remarque 2</u>: la courbe de corrélée s’appuie essentiellement sur la période du soleil (T) et deux autres paramètres A et B (avec B~T/2). Il s’agit de somme de sinusoïdes et de lois normales de répartition s’exprimant sur ces sinusoïdes sur des fréquences précises (moyennant l’application de certains offsets en x & y). </i>

==Critique de la méthodologie et des paramètres choisis==

Nous soulevons des biais expérimentaux, liés à :

* la qualité de la base de données:
** la validité des informations collectées,
** la réalité non réductible simplement des phénomènes observés,
** l'exhaustivité des informations collectées.

Tous ces biais sont directements conditionnés par la qualité du matériel expérimental utilisé. D'où, l'extrême importance qu'il faut y accorder.

==Modélisation par apprentissage: réseau neuronal ==

De façon amusante, une modélisation semble tout à fait possible par l'utilisation d'un [[réseau neuronal]].

Nous proposons d'approximer le signal ci-dessus, à partir d'un réseau neuronal.

Il y a un côté arbitraire dans le choix du RN (couches, neurones, fonctions de transfert), qu'il doit être ultérieurement possible de corréler à des variables environnementales.

← Version précédente		Version du 01:22, 14 avril 2007
Ligne 1:		Ligne 1:
−	=modélisation du « signal »=	+	=Jeu de modélisation du « signal »=

	Nous essayons ici de reconstruire la courbe du signal créé par les observation sur plusieurs dizaines d'années. En s'appuyant sur des harmoniques, il est possible de dégager des fréquences privilégiées.		Nous essayons ici de reconstruire la courbe du signal créé par les observation sur plusieurs dizaines d'années. En s'appuyant sur des harmoniques, il est possible de dégager des fréquences privilégiées.

@@ Ligne 200: / Ligne 200: @@
 Ces biais sont conditionnés par la qualité du matériel expérimental utilisé. D'où, l'extrême importance qu'il faut y accorder.
-[autre biais :...]
+[autres biais :...]
 ==Modélisation par apprentissage: réseau neuronal ==

← Version précédente		Version du 13:16, 18 janvier 2006
Ligne 199:		Ligne 199:

	Ces biais sont conditionnés par la qualité du matériel expérimental utilisé. D'où, l'extrême importance qu'il faut y accorder.		Ces biais sont conditionnés par la qualité du matériel expérimental utilisé. D'où, l'extrême importance qu'il faut y accorder.
		+
		+	[autre biais :...]

	==Modélisation par apprentissage: réseau neuronal ==		==Modélisation par apprentissage: réseau neuronal ==

@@ Ligne 198: / Ligne 198: @@
 ** l'exhaustivité des informations collectées.
-Ces biais sont directement conditionnés par la qualité du matériel expérimental utilisé. D'où, l'extrême importance qu'il faut y accorder.
+Ces biais sont conditionnés par la qualité du matériel expérimental utilisé. D'où, l'extrême importance qu'il faut y accorder.
 ==Modélisation par apprentissage: réseau neuronal ==

@@ Ligne 198: / Ligne 198: @@
 ** l'exhaustivité des informations collectées.
-Tous ces biais sont directement conditionnés par la qualité du matériel expérimental utilisé. D'où, l'extrême importance qu'il faut y accorder.
+Ces biais sont directement conditionnés par la qualité du matériel expérimental utilisé. D'où, l'extrême importance qu'il faut y accorder.
 ==Modélisation par apprentissage: réseau neuronal ==

@@ Ligne 5: / Ligne 5: @@
 Les graphiques présentés ci-après et utilisés pour cette étude sont la propriété de [[Larry Hatch]] (larryhatch(at)sysmatrix.net) que je remercie. Si vous êtes intéressés, vous pouvez trouver plus d'informations sur son site web:
 {{Ref:Web | Titre=*U* UFO Database | Sujet= Un catalogue de données de 18370 objets, représentant plus de 20 années de recherche | URL=http://www.larryhatch.net | Date=02/2005 pour cette étude}}
 [[UFO Database|A propos de la consistance de cette base de données]]