//MINIMOS CUADRADOS POLINOMIO
public ResultadoRegresionPolinomio MinimosCuadradosPolinomio(DatosParametros Datos)
{
ResultadoRegresionPolinomio Resultado = new ResultadoRegresionPolinomio();
//ResultadoRegresion resulta = new ResultadoRegresion();
ResultadosCalcular Calculo = Calcular(Datos);
double[,] Matriz = new double[Datos.Grado + 1, Datos.Grado + 2];
double ST = 0;
double SR = 0;
double R = 0;
for (int i = 0; i < Datos.NumPares; i++)
{
for (int j = 0; j < Datos.Grado + 1; j++)
{
for (int k = 0; k < Datos.Grado + 1; k++)
{
Matriz[j, k] += Math.Pow(Datos.X[i], k + j);
}
Matriz[j, Datos.Grado + 1] += (Datos.Y[i] * Math.Pow(Datos.X[i], j));
}
}
Resultado.ResultadosX = ObtenerGaussJordan(Matriz, Datos.Grado + 1);
double S = 0;
for (int di = 0; di < Datos.NumPares; di++)
{
ST = ST + Math.Pow((Calculo.SumaY / Datos.NumPares) - Datos.Y[di], 2);
S = 0;
for (int j = 0; j < Datos.Grado + 1; j++)
{
S = S + (Resultado.ResultadosX[j] * Math.Pow(Datos.X[di], j));
}
SR = SR + Math.Pow(S - Datos.Y[di], 2);
}
R = Math.Sqrt(Math.Abs(ST - SR) / ST) * 100;
Resultado.Efectividad = R;
return(Resultado);
}
//MINIMOS CUADRADOS
public ResultadoRegresion MinimosCuadrados(DatosParametros Datos)
{
ResultadoRegresion Resultados = new ResultadoRegresion();
ResultadosCalcular Calculo = Calcular(Datos);
double Nominador = (Datos.NumPares * Calculo.SumaXporY) - (Calculo.SumaX * Calculo.SumaY);
double Denominador = (Datos.NumPares * Calculo.SumaXCuadrados) - Math.Pow(Calculo.SumaX, 2);
Resultados.Pendiente = Nominador / Denominador;
Resultados.OrdenadaOrigen = Calculo.PromedioY - (Resultados.Pendiente * Calculo.PromedioX);
double[] lista = new double[5];
Resultados.Efectividad = CalcularR(Datos, Calculo, Resultados, lista);
return(Resultados);
}
public ResultadosCalcular Calcular(DatosParametros Datos)
{
ResultadosCalcular resultados = new ResultadosCalcular();
for (int i = 0; i < Datos.NumPares; i++)
{
resultados.SumaX = resultados.SumaX + Datos.X[i];
resultados.SumaY = resultados.SumaY + Datos.Y[i];
resultados.SumaXporY = resultados.SumaXporY + (Datos.X[i] * Datos.Y[i]);
resultados.SumaXCuadrados = resultados.SumaXCuadrados + (Math.Pow(Datos.X[i], 2));
}
resultados.PromedioX = resultados.SumaX / Datos.NumPares;
resultados.PromedioY = resultados.SumaY / Datos.NumPares;
return(resultados);
}
private double CalcularR(DatosParametros datos, ResultadosCalcular calculo, ResultadoRegresion regresion, double[] lista)
{
double st = 0;
double sr = 0;
double r = 0;
foreach (var item in datos.Y)
{
double aux = calculo.PromedioY - item;
st = st + Math.Pow(aux, 2);
}
for (int i = 0; i < datos.NumPares; i++)
{
double aux2 = (regresion.Pendiente * datos.X[i]) + regresion.OrdenadaOrigen - datos.Y[i];
sr = sr + Math.Pow(aux2, 2);
}
r = Math.Sqrt(Math.Abs(st - sr) / st) * 100;
return(r);
}
//LAGRANJE
public ResultadoLagranje InterpolacionPolinomioLagranje(DatosParametros Datos, double ValorX)
{
ResultadoLagranje Resultado = new ResultadoLagranje();
double[] ResultadosP = new double[Datos.NumPares];
for (int z = 0; z < Datos.NumPares; z++)
{
double Nominadores = 1;
double Denominador = 1;
for (int i = 0; i < Datos.NumPares; i++)
{
if (i != z)
{
Nominadores = Nominadores * (ValorX - Datos.X[i]);
Denominador = Denominador * (Datos.X[z] - Datos.X[i]);
}
}
Resultado.Interpolacion += Datos.Y[z] * (Nominadores / Denominador);
}
return(Resultado);
}