Exemplo n.º 1
0
        void CalculGaussCoefs(int N, out double[] LegendreRoots, out double[] GaussCoefs)
        {
            LegendrePolynom Polynom = new LegendrePolynom(N, 1e-9);

            double[,] Matrix = new double[N, N + 1];

            for (int i = 0; i < N; i++)
            {
                for (int j = 0; j < N; j++)
                {
                    if (i == 0)
                    {
                        Matrix[i, j] = 1;
                    }
                    else
                    {
                        Matrix[i, j] = Math.Pow(Polynom.Roots[j], i);
                    }
                }

                if (i % 2 == 0)
                {
                    Matrix[i, N] = 2.0 / (i + 1);
                }
                else
                {
                    Matrix[i, N] = 0;
                }
            }
            GaussCoefs    = SLE.GaussSolve(Matrix, N);
            LegendreRoots = Polynom.Roots;
        }
Exemplo n.º 2
0
 public List <double> SeekIterations(SLE m, List <double> answers, int n)
 {
     for (int i = 0; i < n; i++)
     {
         logger.Log(answers);
         answers = Iterate(m, answers);
     }
     return(answers);
 }
Exemplo n.º 3
0
 public List <double> Seek(SLE m, List <double> answers)
 {
     while (!Stop(m, answers))
     {
         logger.Log(answers);
         answers = Iterate(m, answers);
     }
     return(answers);
 }
Exemplo n.º 4
0
        protected bool Stop(SLE m, List <double> answers)
        {
            var newAnswers = new List <double>(Iterate(m, answers));

            for (int i = 0; i < m.Right.Count; i++)
            {
                if (Math.Abs(answers[i] - newAnswers[i]) >= precision)
                {
                    return(false);
                }
            }
            return(true);
        }
Exemplo n.º 5
0
        static void Input()
        {
            Console.WriteLine("Введите количество неизвестных");
            int n = int.Parse(Console.ReadLine());

            Console.WriteLine("Введите коэффиценты матрицы системы");
            double[,] matrix = new double[n, n];
            for (int i = 0; i < n; i++)
            {
                string[] str = Console.ReadLine().Split();
                for (int j = 0; j < n; j++)
                {
                    matrix[i, j] = Convert.ToDouble(str[j]);
                }
            }

            Console.WriteLine("Введите вектор свободных коэффицентов");
            double[] value = new double[n];
            string[] temp  = Console.ReadLine().Split();
            for (int i = 0; i < n; i++)
            {
                value[i] = double.Parse(temp[i]);
            }

            SLE system = new SLE(new Matrix(matrix), new Vector(value));

            Console.WriteLine();

            Console.WriteLine("Метод Крамера");
            system.Kramer();

            Console.WriteLine("Метод обратной матрицы");
            system.InvertibleMatrix();

            Console.WriteLine("Метод Гаусса");
            system.Gauss();

            Console.WriteLine("Метод прогонки");
            system.TridiagonalMatrixAlgorithm();

            Console.WriteLine("Метод квадратных корней");
            system.CholeskyDecomposition();

            Console.WriteLine("Метод простых итераций");
            system.Itera(0.001);

            Console.WriteLine("Метод Зейделя");
            system.Seidel(0.001);
        }
Exemplo n.º 6
0
        static double Gauss(double p)
        {
            double Result = 0;

            double[,] Matrix = new double[N, N + 1];

            LegendrePolynom.GetValueAt(0);
            for (int i = 0; i < N; i++)
            {
                for (int j = 0; j < N; j++)
                {
                    if (i == 0)
                    {
                        Matrix[i, j] = 1;
                    }
                    else
                    {
                        Matrix[i, j] = Math.Pow(LegendrePolynom.Roots[j], i);
                    }
                }

                if (i % 2 == 0)
                {
                    Matrix[i, N] = 2.0 / (i + 1);
                }
                else
                {
                    Matrix[i, N] = 0;
                }
            }
            double[] ResultColumn = SLE.GaussSolve(Matrix, N);

            for (int i = 0; i < N; i++)
            {
                double z = 0.5 + LegendreCoefs[i] * 0.5;
                Result += ResultColumn[i] * Integrand(p, z);
            }
            Result *= 0.5;
            return(Result);
        }
Exemplo n.º 7
0
 abstract public List <double> Iterate(SLE m, List <double> answers);