public string Polinomial(double[] xi, double[] yi, int xl)
{
int n = xi.Length;
double[,] matriz = new double[xl + 1, xl + 2];
double sx = 0, sy = 0, st = 0, sr = 0;
for (int i = 0; i < n; i++)
{
sx += xi[i];
sy += yi[i];
}
double promy = sy / n;
for (int i = 0; i < n; i++)
{
for (int j = 0; j < xl + 1; j++)
{
for (int k = 0; k < xl + 1; k++)
{
matriz[j, k] += Math.Pow(xi[i], j + k);
}
matriz[j, xl + 1] += yi[i] * Math.Pow(xi[i], j);
}
}
var a = new Gauss();
double[] sol = a.GaussJordan(matriz, xl + 1);//Cambio xl por xl+1
for (int i = 0; i < n; i++)
{
st += Math.Pow(promy - yi[i], 2); //st += Math.Pow(promy / yi[i], 2)
double s = 0;
for (int j = 0; j < xl + 1; j++)
{
s += (sol[j] * Math.Pow(xi[i], j)); //s += (Math.Pow(sol[j], i) * Math.Pow(xi[i], j))
}
sr += Math.Pow((s - yi[i]), 2);
}
//sr += yi[0];
double r = Math.Sqrt((st - sr) / st) * 100;
Coeficiente = r;
string SolF = null;
for (int i = xl + 1; i > 0; i--) //for (int i = 0; i < xl+1; i++)
{
if (i == 1)
{
SolF = SolF + Convert.ToString(Math.Round(sol[i - 1], 4)) + (i - 1) + " ";
}
else
{
SolF = SolF + Convert.ToString(Math.Round(sol[i - 1], 4)) + "X" + (i - 1) + " + ";
}
}
return(SolF);
}
public string[] RegresionPolinomial(int puntos, double[] xi, double[] yi) //NO ANDA
{
string[] final;
double[] rdo = null;
double CoeCo = 0;
int n = 2;
while (CoeCo < 80) //Mientras el coeficiente sea menor que 80, vuelvue a resolver aumentando la potencia.
{
n = n + 1; //Incrementa la potencia.
double[,] matriz = new double[n, n + 1]; //armando la matriz
double x;
double y = 0;
for (int k = 0; k < puntos; k++)
{
y = y + yi[k]; //Sumatoria de Yi
}
for (int i = 0; i < n; i++)
{
for (int j = 0; j <= n; j++) //empieza a recorrerla matriz de nxn+1 elementos, es decir grado n-1.
{
if (i == 0 && j == 0) //Pregunta si es el primero.
{
matriz[0, 0] = puntos; //El primer valor es el unico diferente.
}
else //No es el primero, sigue operando normalmente.
{
double valf = 0;
for (int k = 0; k < puntos; k++)
{
if (j < n) //Pregunta si el valor no es el ultimo de cada fila.
{
x = Math.Pow(xi[k], i + j);
valf = x; //Calcula el valor de esa pos.
}
else
{
x = Math.Pow(xi[k], i); // Calcula el valor de la ultima pos de cada fila.
valf = y * x;
}
}
matriz[i, j] = valf; //Arma la matriz con los valores.
}
}
}
var a = new Gauss();
rdo = a.GaussJordan(matriz, n); //Resuelve la matriz con el metodo Gauss-Jordan.
double recP = y / puntos; //Recta promedio entre los puntos
double sr = 0, st = 0;
double sri = 0;
double xres = 0;
for (int k = 0; k < puntos; k++)
{
for (int i = 1; i < n; i++)
{
xres = xres - rdo[i] * xi[k];
}
sri = sri + Math.Pow(yi[k] - n - xres, 2); //Valor de la funcion en ese punto (xi).
st = st + Math.Abs(yi[k] - recP); //Dif de los puntos con la recta promedio.
sr = sr + Math.Abs(yi[k] - sri); //Dif de los puntos con la funcion calculada.
}
CoeCo = Math.Sqrt((st - sr) / st) * 100; //Coeficiente de correlacion.
}
final = new string[n];
for (int i = 0; i < n; i++)
{
final[i] = Convert.ToString(rdo[i]);
}
return(final);
}