Ejemplo n.º 1
0
        private void GEO_Click(object sender, EventArgs e)
        {
            //zmienne od użytkownika
            double a     = Convert.ToDouble(textBoxA.Text);
            double b     = Convert.ToDouble(textBoxB.Text);
            double d     = Convert.ToDouble(textBoxD.Text);
            string dHelp = textBoxD.Text.Split(',').Last();
            double l     = Math.Ceiling(Math.Log(((1 / d) * (b - a) + 1), 2)); // obliczanie l czyli na ilu bitach się liczba zmiejsci
            int    T     = Convert.ToInt32(textBoxT.Text);
            double tau   = Convert.ToDouble(textBoxTau.Text);


            List <double> xreal       = new List <double>();
            List <double> xrealChange = new List <double>();
            List <double> fx          = new List <double>();
            List <int>    xint        = new List <int>();
            //List<string> xbin = new List<string>();
            string xbin;
            int    i       = 0;
            int    counter = 0;
            //najlepsze wyniki
            double xrealbest = 0;
            double fxbest    = 0;
            string xbinbest  = "0";

            string xbinb;
            //wyniki po zmianie bitu
            List <string> xbinChange = new List <string>();
            //rangi
            List <int> rangs = new List <int>();
            //Losowe liczby
            List <double> randomNumbers = new List <double>();

            List <double> px = new List <double>(); //Prawdopodobieństwo wystąpienia zmiany
            List <int>    r  = new List <int>();    //nadanie rang


            ConvertFromAndToDecimal convertFromAndToDecimal = new ConvertFromAndToDecimal();
            GEO geo = new GEO();

            //czyszczenie serii dla wykresu, zeby dla kazdego nowego start generował od zera
            chart2.Titles.Clear();
            chart2.Series["Fx"].Points.Clear();
            chart2.Series["FxBest"].Points.Clear();
            chart2.Titles.Add("Historia");
            counter = 0;

            xreal = InitPopulation(a, b, 1);

            do
            {
                fx   = CalculateFunction(xreal);
                xint = FromXrealToXint(xreal, a, b, l);
                xbin = convertFromAndToDecimal.intToBinarty(xint[0], l);

                counter++;
                //przypisanie najlepszego wyniku
                if (fxbest < fx[0])
                {
                    xrealbest = xreal[0];
                    fxbest    = fx[0];
                    xbinbest  = xbin;
                    //generowanie pliku
                    geo.generateFileGEO("plikGEO.txt", xrealbest, xbinbest, fxbest, counter);
                }

                chart2.Series["Fx"].Points.AddXY((counter), fx[0]);
                chart2.Series["FxBest"].Points.AddXY((counter), fxbest);

                //zmiana bitów
                xbinChange.Clear();
                for (i = 0; i < xbin.Length; i++)
                {
                    xbinChange.Add(geo.changeBit(xbin, i));
                }

                //zmina z xbin -> xint aby wyznaczyć wartości funkcji
                xint.Clear();
                foreach (var item in xbinChange)    // wygenerowanie xint (xbin -> xint)
                {
                    xint.Add(convertFromAndToDecimal.convertFrom(item, 2));
                }

                xrealChange = FromXintToXreal(xint, a, b, l); //wygenerowanie z xint do xreal

                //obliczanie funkji
                fx = CalculateFunction(xrealChange);

                //nadawanie rangi
                rangs = geo.doRang(fx);

                //obliczanie prawdopodobieństwa
                px.Clear();
                for (i = 0; i < fx.Count(); i++)
                {
                    px.Add(geo.px(i + 1, tau));
                }
                //losowanie l liczb z zakresu 0-1
                randomNumbers = RandomNumber(l);

                xbinb = "";
                for (i = 0; i < xbin.Length; i++)
                {
                    xbinb = xbinb + geo.doMutationGEO(xbin[i], randomNumbers[i], px[i]);
                }

                xint.Clear();
                xint.Add(convertFromAndToDecimal.convertFrom(xbinb, 2)); //xbin --> xint)

                xreal = FromXintToXreal(xint, a, b, l);                  //wygenerowanie z xint do xreal
                T--;
            } while (T > 0);


            richTextBox.Text = "GEO Najlepszy wynik: \nxreal: " + xrealbest + " \nxbin: " + xbinbest + " \nfx: " + fxbest;
        }
Ejemplo n.º 2
0
        //lab5
        private void AlgorytmGenetyczny_Click(object sender, EventArgs e)
        {
            //do pomiaru czasu
            Stopwatch stopWatch = new Stopwatch();
            TimeSpan  ts;
            string    elapsedTime;

            dataGridView1.Rows.Clear();
            double a     = Convert.ToDouble(textBoxA.Text);
            double b     = Convert.ToDouble(textBoxB.Text);
            double d     = Convert.ToDouble(textBoxD.Text);
            string dHelp = textBoxD.Text.Split(',').Last();

            double n = Convert.ToDouble(textBoxN.Text);
            double l = Math.Ceiling(Math.Log(((1 / d) * (b - a) + 1), 2)); // obliczanie l czyli na ilu bitach się liczba zmiejsci

            double Pk = Convert.ToDouble(textBoxPk.Text);
            double Pm = Convert.ToDouble(textBoxPm.Text);
            int    T  = Convert.ToInt32(textBoxT.Text);

            Boolean eliteChecked = checkBoxEilte.Checked;

            List <double> xreal  = new List <double>();
            List <double> fx     = new List <double>();
            List <double> gx     = new List <double>();
            List <double> px     = new List <double>();
            List <double> qx     = new List <double>();
            List <double> r      = new List <double>();
            List <double> xreal2 = new List <double>();


            List <int>    xreal2int    = new List <int>();
            List <string> xreal2bin    = new List <string>();
            List <string> xParents     = new List <string>();
            List <int>    Pc           = new List <int>();
            List <string> children     = new List <string>();
            List <string> xbinChildren = new List <string>();

            List <List <int> > pointsMutation = new List <List <int> >();
            List <string>      xbinMutation   = new List <string>();

            List <int>    xreal3int = new List <int>();
            List <double> xreal3    = new List <double>();

            List <double> fxmax = new List <double>();
            List <double> fxavg = new List <double>();
            List <double> fxmin = new List <double>();

            double elitefx    = 0;
            double elitexreal = 0;
            //List<double> elitefxlist = new List<double>();
            // List<double> elitexreallist = new List<double>();

            Dictionary <double, int> dict = new Dictionary <double, int>(); //do obliczania wystapien danego elementu

            //do wyników ostatecznych
            List <double> resultXreal      = new List <double>();
            List <double> resultFx         = new List <double>();
            List <double> resultOccurrence = new List <double>();

            GenerateFileModul8(a, b, d, n, T, Pk, Pm);//Generowanie pliku do modulu 8

            int j          = 0;
            int el         = 0;
            int elitePoint = 0;

            Lab3 lab3 = new Lab3(dHelp.Length);

            ConvertFromAndToDecimal convertFromAndToDecimal = new ConvertFromAndToDecimal();

            Lab4 lab4 = new Lab4();

            stopWatch.Start();

            xreal = InitPopulation(a, b, n); //wygenerowanie populacji xReal


            for (int z = 0; z < T; z++)
            {
                fx = CalculateFunction(xreal); // wyliczenie funkcji f(x)

                if (eliteChecked)
                {
                    if (elitefx > fx.Max())
                    {
                        //wstawianie elity
                        elitePoint = random.Next(0, (int)n);
                        fx[el]     = elitefx;
                        xreal3[el] = elitexreal;
                    }
                    else
                    {
                        //generowanie elity
                        el      = 0;
                        elitefx = fx.Max();
                        foreach (var elite in fx)
                        {
                            if (fx.Max() == elite)
                            {
                                elitexreal = xreal[el];
                                break;
                            }
                            el++;
                        }
                    }
                }
                WriteGenerationForModul8(xreal, fx, elitexreal, z); //zapisywanie do pliku pokolenia

                fxmin.Add(fx.Min());
                fxmax.Add(fx.Max());
                fxavg.Add(fx.Average());

                gx     = lab3.searchGx(fx, d);
                px     = lab3.countPx(gx);
                qx     = lab3.countQx(px);
                r      = RandomNumber(n);
                xreal2 = lab3.pickXreal(xreal, qx, r, n); // populacja po selekcji

                //4laby

                xreal2int = FromXrealToXint(xreal2, a, b, l); //wygenerowanie xreal2int (xreal2 -> xint)
                xreal2bin.Clear();
                foreach (var item in xreal2int)               // wygenerowanie ciagu binarnego (xreal2int -> xreal2bin)
                {
                    //xreal2bin.Add(convertFromAndToDecimal.convertTo(item, 2, l));
                    xreal2bin.Add(convertFromAndToDecimal.intToBinarty(item, l));
                }

                r = RandomNumber(n);

                xParents = lab4.serachParets(xreal2bin, r, Pk);  //generowania pokolenia rodziców

                Pc = lab4.searchSection(xParents, l);

                children     = lab4.doChildren(xParents, Pc, n);
                xbinChildren = lab4.doXbinChildren(children, xreal2bin);

                pointsMutation.Clear();
                foreach (var item in xbinChildren)
                {
                    r = RandomNumber(l);
                    pointsMutation.Add(lab4.searchMutation(r, Pm));
                }

                j = 0;
                xbinMutation.Clear();
                foreach (var item in xbinChildren)
                {
                    xbinMutation.Add(lab4.doMutation(item, pointsMutation[j]));
                    j++;
                }
                xreal3int.Clear();
                foreach (var item in xbinMutation)    // wygenerowanie xint (xbin -> xint)
                {
                    xreal3int.Add(convertFromAndToDecimal.convertFrom(item, 2));
                }

                xreal3 = FromXintToXreal(xreal3int, a, b, l); //wygenerowanie z xint1 do xreal1

                fx = CalculateFunction(xreal3);               // wyliczenie funkcji f(x)

                // elitefxlist.Add(elitefx);
                // elitexreallist.Add(elitexreal);

                xreal = xreal3;
            }

            //czyszczenie serii dla wykresu, zeby dla kazdego nowego start generował od zera
            chart1.Titles.Clear();
            chart1.Series["FxMAX"].Points.Clear();
            chart1.Series["FxAVG"].Points.Clear();
            chart1.Series["FxMIN"].Points.Clear();

            Chart_Draw(fxmax, fxavg, fxmin);

            dict = countOccurrence(xreal3);

            foreach (var item in dict)
            {
                resultXreal.Add(item.Key);

                resultOccurrence.Add((((double)item.Value / (double)xreal.Count()) * 100));
            }

            resultFx = CalculateFunction(resultXreal);

            for (int i = 0; i < resultXreal.Count; i++)
            {
                genTable.Add(new GenerateTable(i + 1, resultXreal[i], resultFx[i], resultOccurrence[i]));
            }
            dataGridView1.DataSource = genTable;
            dataGridView1.DataSource = generateTableBindingSource;

            ((CurrencyManager)BindingContext[genTable]).Refresh();

            stopWatch.Stop();
            ts = stopWatch.Elapsed;

            // Format and display the TimeSpan value.
            elapsedTime = String.Format("{0:00}:{1:00}:{2:00}.{3:00}",
                                        ts.Hours, ts.Minutes, ts.Seconds,
                                        ts.Milliseconds / 10);


            //Generowanie pliku do modulu 9
            GenerateFileModul9(a, b, d, n, T, Pk, Pm, fxmax, fxavg, fxmin, elapsedTime);
        }
Ejemplo n.º 3
0
        private void ANW_Click(object sender, EventArgs e)
        {
            //zmienne od użytkownika
            double a     = Convert.ToDouble(textBoxA.Text);
            double b     = Convert.ToDouble(textBoxB.Text);
            double d     = Convert.ToDouble(textBoxD.Text);
            string dHelp = textBoxD.Text.Split(',').Last();
            double l     = Math.Ceiling(Math.Log(((1 / d) * (b - a) + 1), 2)); // obliczanie l czyli na ilu bitach się liczba zmiejsci
            int    T     = Convert.ToInt32(textBoxT.Text);

            List <double> xreal       = new List <double>();
            List <double> xrealChange = new List <double>();
            double        xrealc      = 0;
            List <double> fx          = new List <double>();
            List <double> fxc         = new List <double>();
            List <double> fxlocal     = new List <double>();

            double     fxcc = 0;
            List <int> xint = new List <int>();
            //wyniki po zmianie bitu
            List <string> xbinChange = new List <string>();
            string        xbin       = "";
            string        xbinc      = "";

            //zmienne do zapisania najlepszych wyników
            double xrealbest = 0;
            double fxbest    = 0;
            string xbinbest  = "";

            //zmienne pętle
            int    t       = 0;
            int    i       = 0;
            double w       = 0.0;
            int    licznik = 0;
            double measure = 0.0;
            bool   local   = false;

            ConvertFromAndToDecimal convertFromAndToDecimal = new ConvertFromAndToDecimal();
            GEO geo = new GEO();

            chart3.Series.Clear();
            chart3.Titles.Clear();
            chart3.Series.Add("FxBest");
            chart3.Series["FxBest"].ChartType  = System.Windows.Forms.DataVisualization.Charting.SeriesChartType.Line;
            chart3.Series["FxBest"].XValueType = System.Windows.Forms.DataVisualization.Charting.ChartValueType.Int32;
            chart3.Titles.Add("Historia");

            for (t = 0; t < T; t++)
            {
                xreal   = InitPopulation(a, b, 1);
                local   = false;
                licznik = 0;
                while (!local)
                {
                    fxc  = CalculateFunction(xreal);
                    xint = FromXrealToXint(xreal, a, b, l);
                    xbin = convertFromAndToDecimal.intToBinarty(xint[0], l);

                    //zmiana bitów
                    xbinChange.Clear();
                    for (i = 0; i < xbin.Length; i++)
                    {
                        xbinChange.Add(geo.changeBit(xbin, i));
                    }

                    //zmina z xbin -> xint aby wyznaczyć wartości funkcji
                    xint.Clear();
                    foreach (var item in xbinChange)    // wygenerowanie xint (xbin -> xint)
                    {
                        xint.Add(convertFromAndToDecimal.convertFrom(item, 2));
                    }

                    xrealChange = FromXintToXreal(xint, a, b, l); //wygenerowanie z xint do xreal

                    //obliczanie funkji
                    fx = CalculateFunction(xrealChange);

                    i = 0;
                    //obliczanie na którym elemencie jest najlepsze fx
                    foreach (var os in fx)
                    {
                        if (os == fx.Max())
                        {
                            break;
                        }
                        i++;
                    }

                    //sprawdzanie warunku stopu
                    if (fxc[0] < fx.Max())
                    {
                        xreal[0] = xrealChange[i];
                        fxc[0]   = fx.Max();


                        xrealc = xrealChange[i];
                        fxcc   = fx.Max();
                        xbinc  = xbinChange[i];
                    }
                    else
                    {
                        local = true;
                    }
                    fxlocal.Add(fxcc);
                    licznik++;

                    geo.generateFileANW("plikANW.txt", xrealc, xbinc, fxcc, licznik, t);
                }

                measure = (double)1 / (double)licznik;
                i       = 0;
                chart3.Series.Add("local" + t);
                chart3.Series["local" + t].ChartType = System.Windows.Forms.DataVisualization.Charting.SeriesChartType.Line;
                for (w = measure; w <= 1; w = w + measure)
                {
                    //dodawanie lokalnych wartości do wykresu
                    chart3.Series["local" + t].Points.AddXY((t + w), fxlocal[i]);
                    if (fxbest < fxlocal[i])
                    {
                        chart3.Series["FxBest"].Points.AddXY((t + w), fxlocal[i]);
                    }

                    i++;
                }
                fxlocal.Clear();

                if (fxbest < fxcc)
                {
                    xrealbest = xrealc;
                    fxbest    = fxcc;
                    xbinbest  = xbinc;
                }

                chart3.Series["FxBest"].Points.AddXY(t + 1, fxbest);
            }


            //wypisywanie na okienko
            richTextBox.Text = "ANW Najlepszy wynik: \nxreal: " + xrealbest + " \nxbin: " + xbinbest + " \nfx: " + fxbest;
        }