Ejemplo n.º 1
0
        static void Main()
        {
            newtonRaphson SL = new newtonRaphson(0.0001, 0.01, 2, 2);
            GR = new GoldenRatio();
            for (SL.k = 0; ; SL.k++)
            {

                Console.WriteLine();
                Console.WriteLine("--------------------------------------");
                Console.WriteLine("\nточка:" + SL.next.x + " " + SL.next.y);
                Console.WriteLine("итерация: " + SL.k);
                Console.WriteLine("градиент: " + SL.grad.x+" "+SL.grad.y);
                Console.WriteLine("норма градиента: " + SL.normV(SL.grad));

                if (SL.normV(SL.grad) <= SL.e1) { SL.answer = SL.next; break; }

                if ((SL.m.Invert().Det()) > 0)
                {
                    SL.d = -((SL.m.Invert()) * SL.grad);

                }
                else
                {
                    SL.d.x = -SL.grad.x;
                    SL.d.y = -SL.grad.y;

                }

                SL.step = GR.Calculate(SL.next, SL.d);
                Console.WriteLine("шаг: " + SL.step);
                //if(k>1000) { break;}

                SL.next.x = SL.first.x + SL.step * SL.d.x;
                SL.next.y = SL.first.y + SL.step * SL.d.y;

                // Console.WriteLine(k);

                if (SL.normV(new Point(SL.next.x - SL.first.x, SL.next.y - SL.first.y)) < SL.e2)
                {
                    if (Math.Abs(SL.F(SL.next.x, SL.next.y) - SL.F(SL.first.x, SL.first.y)) < SL.e2)
                    { SL.answer = SL.next; break; }
                }

                SL.first.x = SL.next.x;
                SL.first.y = SL.next.y;

                SL.grad = SL.gradF(SL.next);

            }

            SL.answer.x = SL.next.x;
            SL.answer.y = SL.next.y;
            Console.WriteLine();
            Console.WriteLine("--------------------------------------");
            Console.WriteLine("--------------------------------------");
            Console.WriteLine("\nОтвет:" + SL.answer.x + " " + SL.answer.y);
            Console.WriteLine("итерация: " + SL.k);
            Console.WriteLine("норма градиента: " +SL.normV(SL.grad));
            Console.Read();
        }
Ejemplo n.º 2
0
        static void Main()
        {
            newtonRaphson SL = new newtonRaphson(0.0001, 0.01, 2, 2);

            GR = new GoldenRatio();
            for (SL.k = 0; ; SL.k++)
            {
                Console.WriteLine();
                Console.WriteLine("--------------------------------------");
                Console.WriteLine("\nточка:" + SL.next.x + " " + SL.next.y);
                Console.WriteLine("итерация: " + SL.k);
                Console.WriteLine("градиент: " + SL.grad.x + " " + SL.grad.y);
                Console.WriteLine("норма градиента: " + SL.normV(SL.grad));

                if (SL.normV(SL.grad) <= SL.e1)
                {
                    SL.answer = SL.next; break;
                }


                if ((SL.m.Invert().Det()) > 0)
                {
                    SL.d = -((SL.m.Invert()) * SL.grad);
                }
                else
                {
                    SL.d.x = -SL.grad.x;
                    SL.d.y = -SL.grad.y;
                }

                SL.step = GR.Calculate(SL.next, SL.d);
                Console.WriteLine("шаг: " + SL.step);
                //if(k>1000) { break;}

                SL.next.x = SL.first.x + SL.step * SL.d.x;
                SL.next.y = SL.first.y + SL.step * SL.d.y;

                // Console.WriteLine(k);



                if (SL.normV(new Point(SL.next.x - SL.first.x, SL.next.y - SL.first.y)) < SL.e2)
                {
                    if (Math.Abs(SL.F(SL.next.x, SL.next.y) - SL.F(SL.first.x, SL.first.y)) < SL.e2)
                    {
                        SL.answer = SL.next; break;
                    }
                }

                SL.first.x = SL.next.x;
                SL.first.y = SL.next.y;

                SL.grad = SL.gradF(SL.next);
            }

            SL.answer.x = SL.next.x;
            SL.answer.y = SL.next.y;
            Console.WriteLine();
            Console.WriteLine("--------------------------------------");
            Console.WriteLine("--------------------------------------");
            Console.WriteLine("\nОтвет:" + SL.answer.x + " " + SL.answer.y);
            Console.WriteLine("итерация: " + SL.k);
            Console.WriteLine("норма градиента: " + SL.normV(SL.grad));
            Console.Read();
        }