Différences

Ci-dessous, les différences entre deux révisions de la page.

--- fractale_de_lyapunov [2019/11/15 04:00] – créée Mushussu
+++ fractale_de_lyapunov [2020/04/07 16:53] (Version actuelle) – Mushussu
@@ Ligne 1: / Ligne 1: @@
 ====== Fractale de Lyapunov ======
+===== C++ =====
+<code c++>
+// Compilation : g++ Lyapunov.cpp -o Lyapunov
+// Execution   : ./Lyapunov "AB" 1.8 2.2 3.6 4.0 200 200 2000 -1
+#include <fstream>
+#include <iostream>
+#include <cstdlib>
+#include <cmath>
+#include <ctime>
+#include <string>
+using namespace std;
+int largeur, hauteur;
+int iteration;
+string sequence;
+float bornes[4];
+float limite;
+float proportion(float valeur, float valeurMin, float valeurMax, float valeurOutMin, float valeurOutMax) {
+    float valeurOut = (valeur - valeurMin) / (valeurMax - valeurMin);
+    valeurOut = valeurOutMin + valeurOut * (valeurOutMax - valeurOutMin);
+    return valeurOut;
+}
+int main(int argc, char *argv[]) {
+    if (argc != 10) {
+        cout << "Format : ./Lyapunov \"AB\" 1.8 2.2 3.6 4.0 200 200 2000 -1" << endl;
+        return 0;
+    }
+    // Initialisation des constantes
+    sequence = argv[1];
+    for (int i = 0; i < 4; i++) {
+        bornes[i] = stof(argv[i + 2]);
+    }
+    largeur = stoi(argv[6]);
+    hauteur = stoi(argv[7]);
+    iteration = stoi(argv[8]);
+    limite = stof(argv[9]);
+    time_t debut, fin;
+    time(&debut);
+    float  tableau[iteration];
+    tableau[0]= 0.5;
+    float a, b;
+    float* sequenceTab[sequence.size()];
+    for (int i = 0; i < sequence.size(); i++) {
+        sequenceTab[i] = (sequence[i] == 'A') ? &a : &b;
+    }
+    float aValeur[largeur];
+    float bValeur[hauteur];
+    for (int i = 0; i < largeur; i++) {
+        aValeur[i] = proportion(i, 0, largeur, bornes[0], bornes[1]);
+    }
+    for (int i = 0; i < hauteur; i++) {
+        bValeur[i] = proportion(i, 0, hauteur, bornes[2], bornes[3]);
+    }
+    // Formation du nom du fichier
+    string nom  =  "Lyapunov_" + sequence;
+    for (int i = 0; i < 4; i++) {
+        string s = "_" + to_string(bornes[i]);
+        int position = s.find('.');
+        if (position != -1) {
+            s[position] = ',';
+        }
+        nom +=  s;
+    }
+    nom += ".pgm";
+    // Création de l'image
+    ofstream image(nom);
+    image << "P2" << endl;
+    image << largeur << " " << hauteur << endl;
+    image << "255" << endl;
+    for (int y = 0; y < hauteur; y++) {
+        b = bValeur[y];
+        for (int x = 0; x < largeur; x++) {
+            a = aValeur[x];
+            for (int k = 1; k < iteration; k++) {
+                tableau[k] = *sequenceTab[k % sequence.size()] * tableau[k - 1] * (1 - tableau[k - 1]);
+            }
+            float exposant = 0;
+            for (int k = 1; k < iteration; k++) {
+              exposant += log(abs(*sequenceTab[k % sequence.size()] * (1 - 2 * tableau[k])));
+            }
+            exposant = exposant / iteration;
+            int c;
+            if ((exposant > 0) || (exposant < limite)) {
+              c = 0;
+            } else {
+              c = floor(proportion(exposant, 0, limite, 255, 0));
+            }
+            image << c << endl;;
+        }
+    }
+    time(&fin);
+    cout << "Durée = " << difftime(fin, debut) << endl;
+    const string fonction = "open " + nom;
+    system(fonction.c_str());
+    return 0;
+}
+</code>
+====== Processing =====
+<code java>
+int iteration = 20000;
+int valMin = 0;
+void setup() {
+size(200, 200);
+}
+void draw() {
+println(hour() + ":" + minute() + ":" + second());
+loadPixels();
+for (int i = 0; i < width; i++) {
+ for (int j = 0; j < height; j++) {
+   float a = map(i, 0, width, 1.8, 2.2);
+   float b = map(j, 0, height, 3.6,4);
+   float seq[] = {a, a, a, a, b, b, b, b, b,b ,b , b,b , b, b, b, b, b, b, b}; // AAABBBBBAB
+   float[]  A = new float[iteration];
+   A[0]= 0.5;
+   for (int k = 1; k < iteration; k++) {
+     A[k] = seq[k % seq.length] * A[k - 1] * (1 - A[k - 1]);
+   }
+   float exposant = 0;
+   for (int k = 1; k < iteration; k++) {
+     exposant += log(abs(seq[k % seq.length] * (1 - 2 * A[k])));
+   }
+   exposant = exposant / iteration;
+   int c;
+   if ((exposant > 0) || (exposant < -1)) {
+     c = 0;
+   } else {
+     c = floor(map(exposant, 0, -1, 255, 0));
+   }
+   pixels[i + j * width] = color(c);
+ }
+}
+updatePixels();
+saveFrame("Lyapunov" + "-" + year() + "-" + month() + "-" + day() + "-" + hour() + "-" + minute() + "-" + second() + ".png");
+println(hour() + ":" + minute() + ":" + second());
+exit();
+}
+</code>
 ===== Liens =====
-https://solarianprogrammer.com/2018/11/19/cpp-reading-writing-bmp-images/
+[[https://solarianprogrammer.com/2018/11/19/cpp-reading-writing-bmp-images|C++ lire et écrire des images BMP]]
+[[https://danielbeard.wordpress.com/2011/06/06/image-saving-code-c|C++ Ecrire une image TGA]]
+=== PPM ===
+[[https://www.youtube.com/watch?v=fbH005SzEMc|C++ Tuto vidéo écrire une image PPM]]
+[[https://www.scratchapixel.com/lessons/digital-imaging/simple-image-manipulations|Manipulation d'images simples]]
+{{tag>c++ Processing sylvain}}