Beispiel #1
0
        private Bitmap MakeLevels(List <double> fx, List <double> x1, List <double> x2, double botBorder, double upBorder, double offset, int width, int height, int size, double deltaX)
        {
            var bitmap   = new Bitmap(width, height);
            var graphics = Graphics.FromImage(bitmap);

            var painter = new PicPainter(graphics, width, height, picOffset);

            painter.ClearGraphics();

            //First step
            painter.DrawCoordinateLines(botBorder, upBorder, offset, size, 0.5d);

            var list = new List <double>();
            //Second step
            var res = new Dictionary <double, List <int> >();

            for (var i = 0; i < fx.Count; i++)
            {
                if (!res.Keys.Contains(fx[i]))
                {
                    res.Add(fx[i], new List <int> {
                        i
                    });
                }
                else
                {
                    res[fx[i]].Add(i);
                }
            }


            foreach (var fun in res)
            {
                if (fun.Value.Count > 1)
                {
                    var alreadyDrawn = new List <int>();
                    for (var i = 0; i < fun.Value.Count - 1; i++)
                    {
                        var minLength = double.PositiveInfinity;
                        var index     = -1;
                        for (var j = 0; j < fun.Value.Count; j++)
                        {
                            if (i == j)
                            {
                                continue;
                            }
                            var _i         = fun.Value[i];
                            var _j         = fun.Value[j];
                            var firstPart  = x1[_i] - x1[_j];
                            var secondPart = x2[_i] - x2[_j];
                            var length     = Math.Sqrt(firstPart * firstPart + secondPart * secondPart);

                            if (minLength > length & !alreadyDrawn.Contains(j))
                            {
                                index     = j;
                                minLength = length;
                            }
                        }
                        if (index == -1)
                        {
                            break;
                        }
                        var iCoord = fun.Value[i];
                        var jCoord = fun.Value[index];
                        painter.DrawLine(x1[iCoord], x2[iCoord], x1[jCoord], x2[jCoord], botBorder, upBorder, offset, Pens.Red);

                        alreadyDrawn.Add(index);
                    }
                }
            }

            return(bitmap);
        }
Beispiel #2
0
        private void HookeJeevesMethodWithPenalty(string expression)
        {
            var image      = new Bitmap(levels);
            var picPainter = new PicPainter(Graphics.FromImage(image), pictureBox.Width, pictureBox.Height, picOffset);

            if (!Double.TryParse(addInfoTextBox.Text, out var deltaX) ||
                !Double.TryParse(botBorderTextBox.Text, out var botBorder) ||
                !Double.TryParse(upBorderTextBox.Text, out var upBorder) ||
                !Double.TryParse(hjX1TextBox.Text, out var x1) ||
                !Double.TryParse(hjX2TextBox.Text, out var x2) ||
                !Double.TryParse(rTextBox.Text, out var r))
            {
                MessageBox.Show("oops");
                return;
            }
            var horizontalVector = 0;
            var verticalVector   = 0;
            var size             = CalcSize(botBorder, upBorder, deltaX);

            size++;
            var offset = CalcOffset(botBorder, upBorder) + picOffset;

            object[,] limits;
            try
            {
                limits = CreateLimits();
            }
            catch (Exception ex)
            {
                MessageBox.Show(ex.Message);
                return;
            }
            var fx = CalcWithPenalty(x1, x2, r, expression, limits);

            while (deltaX > 10E-4)
            {
                double newFx = double.MinValue;
                #region Research
                horizontalVector = 0;
                verticalVector   = 0;
                //Horizontal
                double minusHorizontalFx = CalcWithPenalty(x1, x2 - deltaX, r, expression, limits);
                double plusHorizontalFx  = CalcWithPenalty(x1, x2 + deltaX, r, expression, limits);
                if (minusHorizontalFx < fx || plusHorizontalFx < fx)
                {
                    horizontalVector = minusHorizontalFx <= plusHorizontalFx ? -1 : 1;
                }
                fx = CalcWithPenalty(x1, x2 + deltaX * horizontalVector, r, expression, limits);
                picPainter.DrawLine(x1, x2, x1, x2 + deltaX * horizontalVector, botBorder, upBorder, offset, Pens.Black);
                x2 += deltaX * horizontalVector;

                //Vertical
                double minusVerticalFx = CalcWithPenalty(x1 - deltaX, x2, r, expression, limits);
                double plusVerticalFx  = CalcWithPenalty(x1 + deltaX, x2, r, expression, limits);
                if (minusVerticalFx < fx || plusVerticalFx < fx)
                {
                    verticalVector = minusVerticalFx <= plusVerticalFx ? -1 : 1;
                }
                fx = CalcFunction(expression, x1 + deltaX * verticalVector, x2);
                picPainter.DrawLine(x1, x2, x1 + deltaX * verticalVector, x2, botBorder, upBorder, offset, Pens.Black);
                x1 += deltaX * verticalVector;
                #endregion
                //
                if (horizontalVector == 0 && verticalVector == 0)
                {
                    deltaX /= 10;
                    continue;
                }
                #region Pattern maching
                while (true)
                {
                    newFx = CalcWithPenalty(x1 + deltaX * verticalVector, x2 + deltaX * horizontalVector, r, expression, limits);
                    if (newFx < fx)
                    {
                        picPainter.DrawLine(x1, x2, x1 + deltaX * verticalVector, x2 + deltaX * horizontalVector, botBorder, upBorder, offset, Pens.Black);
                        x2 += deltaX * horizontalVector;
                        x1 += deltaX * verticalVector;
                        fx  = newFx;
                    }
                    else
                    {
                        break;
                    }
                }
                #endregion
            }
            pictureBox.Image = image;
            fxTextBox.Text   = fx.ToString();
            x1TextBox.Text   = x1.ToString();
            x2TextBox.Text   = x2.ToString();
        }
Beispiel #3
0
        private void BruteMethod(string expression)
        {
            var image      = new Bitmap(levels);
            var picPainter = new PicPainter(Graphics.FromImage(image), pictureBox.Width, pictureBox.Height, picOffset);
            var minFx      = Double.MaxValue;
            var minX1      = -1d;
            var minX2      = -1d;

            if (!Double.TryParse(addInfoTextBox.Text, out var deltaX) ||
                !Double.TryParse(botBorderTextBox.Text, out var botBorder) ||
                !Double.TryParse(upBorderTextBox.Text, out var upBorder))
            {
                MessageBox.Show("oops");
                return;
            }
            var size = CalcSize(botBorder, upBorder, deltaX);

            size++;
            object[,] limits;
            try
            {
                limits = CreateLimits();
            }
            catch (Exception ex)
            {
                MessageBox.Show(ex.Message);
                return;
            }
            resultGridView.ColumnCount = size + 1;
            resultGridView.RowCount    = size + 1;
            for (var i = 0; i < size + 1; i++)
            {
                resultGridView.Columns[i].Width = 40;
            }
            CreateNumbering(size, deltaX, botBorder);
            var    offset = CalcOffset(botBorder, upBorder) + picOffset;
            bool   isOk;
            double fx;
            var    x1    = botBorder;
            var    x2    = botBorder;
            var    allFx = new List <double>();
            var    allX1 = new List <double>();
            var    allX2 = new List <double>();
            var    minI  = -1;
            var    minJ  = -1;

            for (var i = 0; i < size; i++)
            {
                for (var j = 0; j < size; j++)
                {
                    try
                    {
                        isOk = CalcLimits(limits, x1, x2);
                        fx   = Math.Round(CalcFunction(expression, x1, x2), 2);
                    }
                    catch (Exception ex)
                    {
                        MessageBox.Show(ex.Message + " " + ex.StackTrace + " Ошибка при вычислении", "Ошибка", MessageBoxButtons.OK, MessageBoxIcon.Error);
                        return;
                    }
                    if (isOk)
                    {
                        picPainter.DrawPoint(x1, x2, botBorder, upBorder, offset);
                        allFx.Add(fx);
                        allX1.Add(x1);
                        allX2.Add(x2);
                        if (minFx > fx)
                        {
                            minFx = fx;
                            minX1 = x1;
                            minX2 = x2;
                            minI  = i;
                            minJ  = j;
                        }
                        resultGridView.Rows[size - i - 1].Cells[j + 1].Style.BackColor = Color.White;
                        resultGridView.Rows[size - i - 1].Cells[j + 1].Value           = fx;
                    }
                    else
                    {
                        resultGridView.Rows[size - i - 1].Cells[j + 1].Value           = fx;
                        resultGridView.Rows[size - i - 1].Cells[j + 1].Style.BackColor = Color.Red;
                    }
                    x2 += deltaX;
                }
                x1 += deltaX;
                x2  = botBorder;
            }
            resultGridView.Rows[size - minI - 1].Cells[minJ + 1].Style.BackColor = Color.Green;
            picPainter.DrawPoint(minX1, minX2, botBorder, upBorder, offset, true);
            fxTextBox.Text = Math.Round(minFx, 2).ToString();
            x1TextBox.Text = Math.Round(minX1, 2).ToString();
            x2TextBox.Text = Math.Round(minX2, 2).ToString();
            for (var i = 0; i < allFx.Count; i++)
            {
                allFx[i] = Math.Round(allFx[i], 2);
            }
            pictureBox.Image = image;
        }
Beispiel #4
0
        private void MonteCarloMethod(string expression)
        {
            var image      = new Bitmap(levels);
            var picPainter = new PicPainter(Graphics.FromImage(image), pictureBox.Width, pictureBox.Height, picOffset);
            var fxList     = new List <double>();
            var fx         = Double.MaxValue;
            var x1         = -1d;
            var x2         = -1d;
            var allFx      = new List <double>();
            var allX1      = new List <double>();
            var allX2      = new List <double>();

            if (!int.TryParse(addInfoTextBox.Text, out var stepsNum) ||
                !Double.TryParse(botBorderTextBox.Text, out var botBorder) ||
                !Double.TryParse(upBorderTextBox.Text, out var upBorder))
            {
                MessageBox.Show("oops");
                return;
            }
            var offset = CalcOffset(botBorder, upBorder) + picOffset;

            object[,] limits;
            try
            {
                limits = CreateLimits();
            }
            catch (Exception ex)
            {
                MessageBox.Show(ex.Message);
                return;
            }
            var rnd        = new Random();
            var pointCount = 0;

            for (var i = 0; i < stepsNum; i++)
            {
                var _x1 = rnd.NextDouble() * (upBorder - botBorder) + botBorder;
                var _x2 = rnd.NextDouble() * (upBorder - botBorder) + botBorder;
                try
                {
                    var isOk = CalcLimits(limits, _x1, _x2);
                    if (isOk)
                    {
                        var tempFx = CalcFunction(expression, _x1, _x2);
                        picPainter.DrawPoint(_x1, _x2, botBorder, upBorder, offset);
                        allFx.Add(tempFx);
                        allX1.Add(_x1);
                        allX2.Add(_x2);
                        if (tempFx < fx)
                        {
                            fx = tempFx;
                            x1 = _x1;
                            x2 = _x2;
                        }
                        pointCount++;
                    }
                }
                catch
                {
                    MessageBox.Show("Ошибка при вычислении", "Ошибка", MessageBoxButtons.OK, MessageBoxIcon.Error);
                    return;
                }
            }
            picPainter.DrawPoint(x1, x2, botBorder, upBorder, offset, true);
            for (var i = 0; i < allFx.Count; i++)
            {
                allFx[i] = Math.Round(allFx[i], 1);
            }
            fxTextBox.Text        = fx.ToString();
            x1TextBox.Text        = x1.ToString();
            x2TextBox.Text        = x2.ToString();
            addResultTextBox.Text = pointCount.ToString();
            pictureBox.Image      = image;
        }