Accueil > > > COURBE DE GUMOWSKI-MIRA
 COURBE DE GUMOWSKI-MIRA
 Information sur la source
 Description
En 1980, deux physiciens I. Gumowski et C. Mira, du centre de recherche CERN de Genève en Suisse utilisèrent le système d'équations xn+1 = Byn + f(xn) yn+1 = -xn + f(xn+1) où : f(x) = Ax + (2*(1-A)x² )/(1+x²) pour simuler la trajectoire de particules se déplaçant à très haute vitesse dans un accélérateur de particules de la forme d'une mince boîte cylindrique de plusieurs mètres de long. Ils donnèrent à la constante A des valeurs entre -1 et 1, à la constante B (beaucoup plus sensible) des valeurs très près de 1 et à xo et yo des valeurs entre -20 et 20. Ils découvrirent à leur grande surprise que les trajectoires issues de ce système et portées sur un plan cartésien produisent des images surprenantes nous rappelant certaines formes de vie marine. Le programme présent ci-dessous permet de sauvegarder les points x et y dans un fichier. J'ai utilisé GnuPlot pour tracer les points sauvegardés dans le fichier output.dat
 Source
- /* Valeurs à essayer éventuellement:
- A=-0.31; B=1; x0=3.0; y0=1.0
- A=-0.23; B=1; x0=0.8; y0=0.4
- A=-0.05; B=1; x0=15; y0=0.0
- A=-0.6; B=0.99; x0=0.09; y0=-2.76
- A=0.31; B=1; x0=3; y0=1
- A=-0.77; B=0.95; x0=3; y=1 */
-
- #include <stdio.h>
- #include <math.h>
- #define A -0.48
- #define B 0.93
-
- int main(int argc, char *argv[])
- {
-
- int i,n;
- float x,y,xp,yp;
- FILE *fichier;
- fichier = fopen("output.dat","w");
-
- x = 4.0;
- y = 4.0;
- n = 30000;
- for(i=1;i<=200;i++) {
- xp = B*y+A*x+2*(1-A)*x*x/(1+xp*xp);
- y = -x+A*xp+2*(1-A)*xp*xp/(1+xp*xp);
- x = xp;
-
- }
- for(i=1;i<=n;i++) {
- xp = B*y+A*x+2*(1-A)*x*x/(1+x*x);
- y = -x+A*xp+2*(1-A)*xp*xp/(1+xp*xp);
- x = xp;
- fprintf(fichier,"%f %f \n",x,y);
- }
- fclose(fichier);
- return 0;
- }
/* Valeurs à essayer éventuellement:
A=-0.31; B=1; x0=3.0; y0=1.0
A=-0.23; B=1; x0=0.8; y0=0.4
A=-0.05; B=1; x0=15; y0=0.0
A=-0.6; B=0.99; x0=0.09; y0=-2.76
A=0.31; B=1; x0=3; y0=1
A=-0.77; B=0.95; x0=3; y=1 */
#include <stdio.h>
#include <math.h>
#define A -0.48
#define B 0.93
int main(int argc, char *argv[])
{
int i,n;
float x,y,xp,yp;
FILE *fichier;
fichier = fopen("output.dat","w");
x = 4.0;
y = 4.0;
n = 30000;
for(i=1;i<=200;i++) {
xp = B*y+A*x+2*(1-A)*x*x/(1+xp*xp);
y = -x+A*xp+2*(1-A)*xp*xp/(1+xp*xp);
x = xp;
}
for(i=1;i<=n;i++) {
xp = B*y+A*x+2*(1-A)*x*x/(1+x*x);
y = -x+A*xp+2*(1-A)*xp*xp/(1+xp*xp);
x = xp;
fprintf(fichier,"%f %f \n",x,y);
}
fclose(fichier);
return 0;
}
 Sources du même auteur
 Sources de la même categorie
 Commentaires et avis
 Discussions en rapport avec ce code source dans le forum
Problème avec le copier coller [ par Unknown ]
J'ai un projet où j'utilise un composant image. Il me permet l'affichage d'une courbe a partir d'un oscilloscope numérique. Mon problème est que je do
faire une courbe toute conne [ par Axool ]
pouvez-vous m'aider svp: j'aimerais savoir faire une courbe toute conne, vraiment le + simple possible...mercips: si au passage vous savez comment cré
desiner une courbe a partir des valeurs d'un tableau [ par aymentri ]
je veut afficher un graphe en utilisant des valeurs stokées dans un tableau
Sismographe [ par nahs ]
Bonjours,Je cherche un moyen de tracer une courbe de la meme fason qu'un sismographe c'est à dire: sur une fenetre on verrai la feuille qui bouge de d
courbe C++ [ par tguinel ]
salut,je veut afficher un graphe en utilisant des valeurs stokées dans une base données(en C++, visualC++ 6.0).merci de m'aider.
Rafraichissement [ par gus2647 ]
Bonjour,J utilise VC++ 6 et les MFC. lorsque je trace une courbe, elle s affiche correctement, mais lorsque je cree par dessus par exemple une CListCt
Comment faire un Defilemant de courbe en API [ par nahs ]
Bonjours,J'aurai voulu saloire comment faire un effet defilent en API; les instructions qui permet de faire cela.Merci
utilisation du Tchart??(builder) [ par kach23 ]
bonjour,je fais une acquisition de mesure(carte NI) avec builder 6. en utilisant le timer de builder je pe realiser une courbe (Tchart) qui evolue au
Line To rapide [ par BarthOlivier ]
Salut,Je fait plusieurs occurence de MoveTo + LineTo sur un CDC.Est-il possible d'optimiser LineTo avec une autre fonction ?J'avais essayé en passant
tracer courbe en excel à partir de MFC [ par mabrouka ]
bonjour,comment tracer une courbe en excel à partir de visual c++ MFCmerci
|
Derniers Blogs
 CSDL FUNCTIONCSDL FUNCTION par Matthieu MEZIL
Dans mon post précédent , j'ai utilisé une CSDL Function afin de générer une requête SQL avec un DateDiff utilisant la date courante sur la BD à partir d'une requête LINQ. Dans le cadre de ce post , vous avez probablement remarqué que dans le cadre de plu...
Cliquez pour lire la suite de l'article par Matthieu MEZIL LINQ TO ENTITIESLINQ TO ENTITIES par Matthieu MEZIL
Cette semaine je suis à Montréal en tant que speaker sur Entity Framework pour l'évènement confoo . J'en profite pour remercier les organisateurs de cet évènement de m'avoir fait confiance et Access-IT de m'avoir permis d'y participer. En parallèle, j'ai ...
Cliquez pour lire la suite de l'article par Matthieu MEZIL FAIRE APPARAITRE L'ONGLET 'DéVELOPPEUR' DANS OFFICE 2010FAIRE APPARAITRE L'ONGLET 'DéVELOPPEUR' DANS OFFICE 2010 par neodante
La nouvelle interface d'Office 2010 à amener quelques modifications par rapport à celle de 2007. Certes mineures, ces modifications ont fait disparaître la case à cocher de l'onglet 'Développeur' en première page du panneau du 'bouton Office' (dans Office...
Cliquez pour lire la suite de l'article par neodante [ASTUCE] PATCH POUR MICROSOFT FORUMS NNTP BRIDGE V1[ASTUCE] PATCH POUR MICROSOFT FORUMS NNTP BRIDGE V1 par pierre
Si vous avez téléchargé comme moi Microsoft Forums NNTP Bridge V1 avant le 11 mars 2010 (voir [Astuce] Disponibilité de Microsoft Forum NNTP Bridge Version 1.0), un problème de date localisée pour les non anglais était présent. Un patch est disponibl...
Cliquez pour lire la suite de l'article par pierre PB LORS DE L'INSTALLATION SHAREPOINT 2010.PB LORS DE L'INSTALLATION SHAREPOINT 2010. par Patrick Guimonet
Lors de l'installation de SharePoint 2010, j'ai rencontré un problème de plantage à l'étape 5 du configuration Wizard. Ca se termine sur cet écran : Et en analysant le fichier de journalisation, on remarque vers la fin des 15000 et quelques lign...
Cliquez pour lire la suite de l'article par Patrick Guimonet
Forum
 RE : HELP PLEEEASRE : HELP PLEEEAS par PhilLopez
Cliquez pour lire la suite par PhilLopez
Logiciels
 Xilisoft Convertisseur Vidéo Ultimate (5.1.39.0305)XILISOFT CONVERTISSEUR VIDéO ULTIMATE (5.1.39.0305)Xilisoft Convertisseur Vidéo Ultimate est un outil puissant de conversion vidéo, facile à utilise... Cliquez pour télécharger Xilisoft Convertisseur Vidéo Ultimate Xilisoft DVD Ripper Ultimate (5.0.64.0304)XILISOFT DVD RIPPER ULTIMATE (5.0.64.0304)Xilisoft DVD Ripper Ultimate est un logiciel excellent pour copier et convertir DVD vers presque ... Cliquez pour télécharger Xilisoft DVD Ripper Ultimate Rigs of Rods (63.3)RIGS OF RODS (63.3)c'est un jeu de multi-simulation camions,autobus voitures, avions, bateaux, hélicoptère avec défo... Cliquez pour télécharger Rigs of Rods Konvertor (4.00)KONVERTOR (4.00)Le logiciel est un gestionnaire multimedia affichant, jouant et convertissant plus de 2000 format... Cliquez pour télécharger Konvertor
|