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
 TECHDAYS PARIS 2012 : COMMENT SHAREPOINT A SAUVé MES TECHDAYSTECHDAYS PARIS 2012 : COMMENT SHAREPOINT A SAUVé MES TECHDAYS par ROMELARD Fabrice
Speakers : Lionel Limozin et Alain Marty La session commence par une découverte de SharePoint à travers la mise en place d'un environnement SharePoint pour la gestion des Sessions animées par BeWise. Le besoin est très ba...
Cliquez pour lire la suite de l'article par ROMELARD Fabrice PERSPECTIVE 3.0 POUR SILVERLIGHT 5.0PERSPECTIVE 3.0 POUR SILVERLIGHT 5.0 par odewit
Je viens de publier la version 3.0 de Perspective pour Silverlight, qui regroupe un portage sous Silverlight 5.0 des fonctionnalités de Perspective 2.0, le framework 3D de haut-niveau introduit récemment et de nouveaux exemples de code. En voici la li...
Cliquez pour lire la suite de l'article par odewit TECHDAYS PARIS 2012 : TOP 10 DES BEST PRACTICES POUR SQL SERVERTECHDAYS PARIS 2012 : TOP 10 DES BEST PRACTICES POUR SQL SERVER par ROMELARD Fabrice
Speaker : Nadia Ben El Kadi Configuration machine La session commence par la toute première question à se poser lors de la mise en place d'environnement SQL Server, la configuration des machines : Type de mac...
Cliquez pour lire la suite de l'article par ROMELARD Fabrice TECHDAYS PARIS 2012 : KINECT + OFFICE 365 UN BON GESTE POUR VOTRE SITECHDAYS PARIS 2012 : KINECT + OFFICE 365 UN BON GESTE POUR VOTRE SI par ROMELARD Fabrice
Speakers : Fabrice Barbin, Samuel Blanchard, Julien Lo Presti Titre Prometteur et attractif invitant à voir comment lier le composant ludique Kinect dans le cadre d'une structure IT classique, notamment au travers de la plat...
Cliquez pour lire la suite de l'article par ROMELARD Fabrice TECHDAYS PARIS 2012 : PLEINIèRE DU PREMIER JOURTECHDAYS PARIS 2012 : PLEINIèRE DU PREMIER JOUR par ROMELARD Fabrice
KeyNotes du premier jour pour les développeurs. La session est principalement axée sur une des principales directions prise par Microsoft à travers tous ses nouveaux produits : Cloud privé ou public (Solution Azure) ...
Cliquez pour lire la suite de l'article par ROMELARD Fabrice
Forum
 C++ C++ par yesoun1
Cliquez pour lire la suite par yesoun1 OPNETOPNET par hth21
Cliquez pour lire la suite par hth21 RE : ARBRE BINAIRERE : ARBRE BINAIRE par pacotheking
Cliquez pour lire la suite par pacotheking
Logiciels
 Academy System (17.2.1.0)ACADEMY SYSTEM (17.2.1.0)Logiciel de gestion des établissements.
- élèves/étudiants (inscription, dossier, absence...)
-... Cliquez pour télécharger Academy System Easy-Planning (1.0.0.1)EASY-PLANNING (1.0.0.1)Basé sur les mêmes principes que MyPlanning, Easy-Planning permet de créer des plannings sous la ... Cliquez pour télécharger Easy-Planning COLLECTOR PLUS (3.00B)COLLECTOR PLUS (3.00B)COLLECTOR PLUS version 3.00B est un logiciel utilisant une base de données alimentée par :
- L... Cliquez pour télécharger COLLECTOR PLUS PONAMEDIA PREMIUM - HELLLOOO FLASH DEMO (V7.4)PONAMEDIA PREMIUM - HELLLOOO FLASH DEMO (V7.4)PONAMEDIA TV DEVIENS HELLLOOO FLASH
LA TV SUR VOTRE ORDINATEUR.
Toute une plateforme Multi... Cliquez pour télécharger PONAMEDIA PREMIUM - HELLLOOO FLASH DEMO LettresFaciles 2011 (8.0.0.1)LETTRESFACILES 2011 (8.0.0.1)LettresFaciles est un logiciel facilitant la création et la rédaction de lettres types.
Son inte... Cliquez pour télécharger LettresFaciles 2011
|