Exemple #1
0
        public static double[] Calcular(int n, string fx, double e, string x_ini)
        {
            double[] xk   = x_ini.SplitToDoubles();
            double[] grad = Gradiente.CalculaGradiente(n, fx, x_ini.SplitToDoubles()).ToArray();
            double[] gk   = (double[])grad.Clone();
            double[] gk_1;
            double[] dk = gk.Negativo();
            double   beta;

            while (Interpretadores.NormaVetor(grad) > e)
            {
                for (int i = 0; i < n; i++)
                {
                    string xk_s = Interpretadores.GeraVetorY(xk, Interpretadores.LambdaVDirec(dk));
                    double lamb = Newton.Calcular(0, ((Interpretadores.GeraFy(fx, xk_s)).Replace(',', '.')).Replace("lamb", "x[1]"), 0.1);

                    xk   = Interpretadores.SubsLambda(lamb, xk_s.Replace(',', '.'));
                    gk_1 = Gradiente.CalculaGradiente(n, fx, xk).ToArray();
                    if (i >= n - 1)
                    {
                        break;
                    }
                    beta = GradConjugado.MultEscalar(gk_1, gk_1) / GradConjugado.MultEscalar(gk, gk);
                    dk   = Interpretadores.SomaVetor(gk_1.Negativo(), dk.MultConstante(beta));
                    gk   = (double[])gk_1.Clone();
                }
                grad = Gradiente.CalculaGradiente(n, fx, xk).ToArray();
            }
            return(xk);
        }
        public static double[] Calcular(int n, string fx, double e, string x_ini)
        {
            double[] x0 = x_ini.SplitToDoubles();
            double[] xk;
            double[] gk = (double[])Gradiente.CalculaGradiente(n, fx, x0).ToArray().Clone();
            double[] gk_1;
            double[] dk;
            double[] qk;
            double[] pk;
            double[] xi;
            double   lambdak;

            double[,] Sk = MatrizSimetricaPositivaDefinida(n);
            int k, i;

            k = i = 0;
volta:
            xk = (double[])x0.Clone();
            while (Interpretadores.NormaVetor(gk) > e)
            {
                dk = GradConjugado.MultMatriz(Sk, gk.ToMatriz()).ToVector().Negativo();

                string vy = Interpretadores.GeraVetorY(xk, Interpretadores.LambdaVDirec(dk));
                string fy = Interpretadores.GeraFy(fx, vy.Replace(',', '.')).Replace("lamb", "x[1]");
                lambdak = Newton.Calcular(3, fy, 0.1);

                xk = Interpretadores.SubsLambda(lambdak, vy.Replace(',', '.'));

                if (k < n - 1)
                {
                    gk_1 = Gradiente.CalculaGradiente(n, fx, xk).ToArray();

                    qk = Interpretadores.SubtracaoVetor(gk_1, gk);

                    pk = dk.MultConstante(lambdak);


                    var parte1 = GradConjugado.MultMatriz(pk.ToMatriz(), pk.TranspVetor()).DivideConstante(GradConjugado.MultEscalar(pk, qk));
                    var parte2 = GradConjugado.MultMatriz(GradConjugado.MultMatriz(GradConjugado.MultMatriz(Sk, qk.ToMatriz()), qk.TranspVetor()), Sk).DivideConstante(GradConjugado.MultEscalar((GradConjugado.MultMatriz(qk.TranspVetor(), Sk)).ToVector(), qk));

                    Sk = SubtracaoMatriz(SomaMatriz(Sk, parte1), parte2);

                    k += 1;
                    gk = (double[])gk_1.Clone();
                }
                else
                {
                    xi = (double[])xk.Clone();
                    i += 1;
                    x0 = (double[])xk.Clone();
                    gk = Gradiente.CalculaGradiente(n, fx, x0).ToArray();
                    k  = 0;
                    goto volta;
                }
            }
            return(xk);
        }