public static double[] Calcular(int n, string fx, double e, string x_ini) { double[] xk = Interpretadores.SplitToDoubles(x_ini); double[] x = new double[n]; for (int i = 0; i < x.Length; i++) { x[i] = double.MaxValue; } double[] Grad = Gradiente.CalculaGradiente(n, fx, xk).ToArray(); double[] W; double[,] Hess; while (Math.Abs(Interpretadores.NormaVetor(Grad)) >= e) { Hess = Hessiano(n, fx, xk); for (int i = 0; i < Grad.Length; i++) { Grad[i] = -Grad[i]; } W = MétodosSistemasLineares.SistemaLinearTriangularSuperior(n, Hess, Grad); for (int i = 0; i < xk.Length; i++) { x[i] = xk[i]; } xk = Interpretadores.SomaVetor(xk, W); if (Interpretadores.NormaVetor(Interpretadores.SubtracaoVetor(xk, x)) <= e) { break; } 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); }
public static double[] Calcular(string func, double[] xk, double erro, int n) { int k = 0; double lambda, norma = 1000; double[] d, yk, vetNorma, x; vetNorma = new double[n]; d = new double[n]; yk = new double[n]; string lambdaD, vetYj, Fy; do { k++; x = (double[])xk.Clone(); yk = (double[])xk.Clone(); for (int j = 0; j < n; j++) { /*Monta o vetor d*/ for (int b = 0; b < n; b++) { d[b] = 0; } d[j] = 1; /**/ //multiplica lambda por d lambdaD = Interpretadores.LambdaVDirec(d); //vetor lambda*d + vetor yk vetYj = Interpretadores.GeraVetorY(yk, lambdaD); //substitui lambda por x nos lugares corretos -> Fy só com lambda de variável Fy = Interpretadores.GeraFy(func, vetYj.Replace(',', '.')); //Newton para calcular o valor de lambda lambda = Newton.Calcular(0, (Fy.Replace("lamb", "x[1]")).Replace(',', '.'), 0.1); //substitui o valor de lambda na função Fy -> gera vetor xk novo yk = Interpretadores.SubsLambda(lambda, vetYj.Replace(',', '.')); } xk = (double[])yk.Clone(); vetNorma = Interpretadores.SubtracaoVetor(xk, x); //resultado da norma norma = Interpretadores.NormaVetor(vetNorma); } while (norma > erro); return(xk); }