Skip to main content
replaced http://pt.stackoverflow.com/ with https://pt.stackoverflow.com/
Fonte Link

Utilizando o mesmo algoritmo da resposta do @Marcos Banikda resposta do @Marcos Banik, escrevi um pequeno algoritmo em C baseado no método dos mínimos quadrados:

void lms(double *x, double *y, int n, double *m, double *b)
{
    int i;
    double sumYX = 0.;
    double sumX = 0.;
    double sumY = 0.;
    double sumX2 = 0.;
    double sum2X = 0.;
    for(i = 0; i < n; i++) {
        sumYX += x[i] * y[i];
        sumX += x[i];
        sumY += y[i];
        sumX2 += x[i] * x[i];
    }
    sum2X = sumX * sumX;
    
    *m = (sumYX - (sumX * sumY) / (double)n) / (sumX2 - sum2X / (double)n);
    *b = sumY / (double)n - *m * sumX / (double)n;
}

O algoritmo é exatamente o mesmo utilizado na resposta dele. O laço for inicial é feito para calcular todos as somatórias e posteriormente seus resultados são utilizados para achar m e b. Acredito que isso seja facilmente traduzível para sua linguagem de preferência. :)

Com m e b você tem a equação da reta. Para calcular qual o valor de y, basta fazer uma função:

double trend(double m, double b, double x)
{
    return m*x + b;
}

Para chamar as funções:

int main(void) 
{
    double x[] = {1., 2., 3., 4.};
    double y[] = {10., 20., 30., 40.};
    double m, b;
    double nx = 5.;
    double ny;

    lms(x, y, 4, &m, &b);
    ny = trend(m, b, nx);
    
    printf("m: %lf \nb: %lf \nnx: %lf \nny: %lf\n", m, b, nx, ny);

    return 0;
}

Utilizando o mesmo algoritmo da resposta do @Marcos Banik, escrevi um pequeno algoritmo em C baseado no método dos mínimos quadrados:

void lms(double *x, double *y, int n, double *m, double *b)
{
    int i;
    double sumYX = 0.;
    double sumX = 0.;
    double sumY = 0.;
    double sumX2 = 0.;
    double sum2X = 0.;
    for(i = 0; i < n; i++) {
        sumYX += x[i] * y[i];
        sumX += x[i];
        sumY += y[i];
        sumX2 += x[i] * x[i];
    }
    sum2X = sumX * sumX;
    
    *m = (sumYX - (sumX * sumY) / (double)n) / (sumX2 - sum2X / (double)n);
    *b = sumY / (double)n - *m * sumX / (double)n;
}

O algoritmo é exatamente o mesmo utilizado na resposta dele. O laço for inicial é feito para calcular todos as somatórias e posteriormente seus resultados são utilizados para achar m e b. Acredito que isso seja facilmente traduzível para sua linguagem de preferência. :)

Com m e b você tem a equação da reta. Para calcular qual o valor de y, basta fazer uma função:

double trend(double m, double b, double x)
{
    return m*x + b;
}

Para chamar as funções:

int main(void) 
{
    double x[] = {1., 2., 3., 4.};
    double y[] = {10., 20., 30., 40.};
    double m, b;
    double nx = 5.;
    double ny;

    lms(x, y, 4, &m, &b);
    ny = trend(m, b, nx);
    
    printf("m: %lf \nb: %lf \nnx: %lf \nny: %lf\n", m, b, nx, ny);

    return 0;
}

Utilizando o mesmo algoritmo da resposta do @Marcos Banik, escrevi um pequeno algoritmo em C baseado no método dos mínimos quadrados:

void lms(double *x, double *y, int n, double *m, double *b)
{
    int i;
    double sumYX = 0.;
    double sumX = 0.;
    double sumY = 0.;
    double sumX2 = 0.;
    double sum2X = 0.;
    for(i = 0; i < n; i++) {
        sumYX += x[i] * y[i];
        sumX += x[i];
        sumY += y[i];
        sumX2 += x[i] * x[i];
    }
    sum2X = sumX * sumX;
    
    *m = (sumYX - (sumX * sumY) / (double)n) / (sumX2 - sum2X / (double)n);
    *b = sumY / (double)n - *m * sumX / (double)n;
}

O algoritmo é exatamente o mesmo utilizado na resposta dele. O laço for inicial é feito para calcular todos as somatórias e posteriormente seus resultados são utilizados para achar m e b. Acredito que isso seja facilmente traduzível para sua linguagem de preferência. :)

Com m e b você tem a equação da reta. Para calcular qual o valor de y, basta fazer uma função:

double trend(double m, double b, double x)
{
    return m*x + b;
}

Para chamar as funções:

int main(void) 
{
    double x[] = {1., 2., 3., 4.};
    double y[] = {10., 20., 30., 40.};
    double m, b;
    double nx = 5.;
    double ny;

    lms(x, y, 4, &m, &b);
    ny = trend(m, b, nx);
    
    printf("m: %lf \nb: %lf \nnx: %lf \nny: %lf\n", m, b, nx, ny);

    return 0;
}
added 290 characters in body
Fonte Link

Utilizando o mesmo algoritmo da resposta do @Marcos Banik, escrevi um pequeno algoritmo em C baseado no métodosmétodo dos mínimos quadrados:

void lms(double *x, double *y, int n, double *m, double *b)
{
    int i;
    double sumYX = 0.;
    double sumX = 0.;
    double sumY = 0.;
    double sumX2 = 0.;
    double sum2X = 0.;
    for(i = 0; i < n; i++) {
        sumYX += x[i] * y[i];
        sumX += x[i];
        sumY += y[i];
        sumX2 += x[i] * x[i];
    }
    sum2X = sumX * sumX;
    
    *m = (sumYX - (sumX * sumY) / (double)n) / (sumX2 - sum2X / (double)n);
    *b = sumY / (double)n - *m * sumX / (double)n;
}

O algoritmo é exatamente o mesmo utilizado na resposta dele. O laço for inicial é feito para calcular todos as somatórias e posteriormente seus resultados são utilizados para achar m e b. Acredito que isso seja facilmente traduzível para sua linguagem de preferência. :)

Com m e b você tem a equação da reta. Para chamarcalcular qual o métodovalor de y, basta fazer uma função:

double trend(double m, double b, double x)
{
    return m*x + b;
}

Para chamar as funções:

int main(void) 
{
    double x[] = {1., 2., 3., 4.};
    double y[] = {10., 20., 30., 40.};
    double m, b;
    double nx = 5.;
    double ny;

    lms(x, y, 4, &m, &b);
    ny = trend(m, b, nx);
    
    printf("m: %f%lf \n\nb: b%lf \nnx: %f"%lf \nny: %lf\n", m, b, nx, ny);

    return 0;
}

Utilizando o mesmo algoritmo da resposta do @Marcos Banik, escrevi um pequeno algoritmo em C baseado no métodos dos mínimos quadrados:

void lms(double *x, double *y, int n, double *m, double *b)
{
    int i;
    double sumYX = 0.;
    double sumX = 0.;
    double sumY = 0.;
    double sumX2 = 0.;
    double sum2X = 0.;
    for(i = 0; i < n; i++) {
        sumYX += x[i] * y[i];
        sumX += x[i];
        sumY += y[i];
        sumX2 += x[i] * x[i];
    }
    sum2X = sumX * sumX;
    
    *m = (sumYX - (sumX * sumY) / (double)n) / (sumX2 - sum2X / (double)n);
    *b = sumY / (double)n - *m * sumX / (double)n;
}

O algoritmo é exatamente o mesmo utilizado na resposta dele. O laço for inicial é feito para calcular todos as somatórias e posteriormente seus resultados são utilizados para achar m e b. Acredito que isso seja facilmente traduzível para sua linguagem de preferência. :)

Para chamar o método:

int main(void) 
{
    double x[] = {1., 2., 3., 4.};
    double y[] = {10., 20., 30., 40.};
    double m, b;

    lms(x, y, 4, &m, &b);
    
    printf("m: %f \n b: %f", m, b);

return 0;
}

Utilizando o mesmo algoritmo da resposta do @Marcos Banik, escrevi um pequeno algoritmo em C baseado no método dos mínimos quadrados:

void lms(double *x, double *y, int n, double *m, double *b)
{
    int i;
    double sumYX = 0.;
    double sumX = 0.;
    double sumY = 0.;
    double sumX2 = 0.;
    double sum2X = 0.;
    for(i = 0; i < n; i++) {
        sumYX += x[i] * y[i];
        sumX += x[i];
        sumY += y[i];
        sumX2 += x[i] * x[i];
    }
    sum2X = sumX * sumX;
    
    *m = (sumYX - (sumX * sumY) / (double)n) / (sumX2 - sum2X / (double)n);
    *b = sumY / (double)n - *m * sumX / (double)n;
}

O algoritmo é exatamente o mesmo utilizado na resposta dele. O laço for inicial é feito para calcular todos as somatórias e posteriormente seus resultados são utilizados para achar m e b. Acredito que isso seja facilmente traduzível para sua linguagem de preferência. :)

Com m e b você tem a equação da reta. Para calcular qual o valor de y, basta fazer uma função:

double trend(double m, double b, double x)
{
    return m*x + b;
}

Para chamar as funções:

int main(void) 
{
    double x[] = {1., 2., 3., 4.};
    double y[] = {10., 20., 30., 40.};
    double m, b;
    double nx = 5.;
    double ny;

    lms(x, y, 4, &m, &b);
    ny = trend(m, b, nx);
    
    printf("m: %lf \nb: %lf \nnx: %lf \nny: %lf\n", m, b, nx, ny);

    return 0;
}
added 16 characters in body
Fonte Link

Utilizando o mesmo algoritmo da resposta do @Marcos Banik, escrevi um pequeno algoritmo em C baseado no métodos dos mínimos quadrados:

void lms(floatdouble *x, floatdouble *y, int n, floatdouble *m, floatdouble *b)
{
    int i;
    floatdouble sumYX = 0.;
    floatdouble sumX = 0.;
    floatdouble sumY = 0.;
    floatdouble sumX2 = 0.;
    floatdouble sum2X = 0.;
    for(i = 0; i < n; i++) {
        sumYX += x[i] * y[i];
        sumX += x[i];
        sumY += y[i];
        sumX2 += x[i] * x[i];
    }
    sum2X = sumX * sumX;
    
    *m = (sumYX - (sumX * sumY) / (floatdouble)n) / (sumX2 - sum2X / (floatdouble)n);
    *b = sumY / (floatdouble)n - *m * sumX / (floatdouble)n;
}

O algoritmo é exatamente o mesmo utilizado na resposta dele. O laço for inicial é feito para calcular todos as somatórias e posteriormente seus resultados são utilizados para achar m e b. Acredito que isso seja facilmente traduzível para sua linguagem de preferência. :)

Para chamar o método:

int main(void) 
{
    floatdouble x[] = {1., 2., 3., 4.};
    floatdouble y[] = {10., 20., 30., 40.};
    floatdouble m, b;

    lms(x, y, 4, &m, &b);
    
    printf("m: %f \n b: %f", m, b);

return 0;
}

Utilizando o mesmo algoritmo da resposta do @Marcos Banik, escrevi um pequeno algoritmo em C baseado no métodos dos mínimos quadrados:

void lms(float *x, float *y, int n, float *m, float *b)
{
    int i;
    float sumYX = 0.;
    float sumX = 0.;
    float sumY = 0.;
    float sumX2 = 0.;
    float sum2X = 0.;
    for(i = 0; i < n; i++) {
        sumYX += x[i] * y[i];
        sumX += x[i];
        sumY += y[i];
        sumX2 += x[i] * x[i];
    }
    sum2X = sumX * sumX;
    
    *m = (sumYX - (sumX * sumY) / (float)n) / (sumX2 - sum2X / (float)n);
    *b = sumY / (float)n - *m * sumX / (float)n;
}

O algoritmo é exatamente o mesmo utilizado na resposta dele. O laço for inicial é feito para calcular todos as somatórias e posteriormente seus resultados são utilizados para achar m e b. Acredito que isso seja facilmente traduzível para sua linguagem de preferência. :)

Para chamar o método:

int main(void) 
{
    float x[] = {1., 2., 3., 4.};
    float y[] = {10., 20., 30., 40.};
    float m, b;

    lms(x, y, 4, &m, &b);
    
    printf("m: %f \n b: %f", m, b);

return 0;
}

Utilizando o mesmo algoritmo da resposta do @Marcos Banik, escrevi um pequeno algoritmo em C baseado no métodos dos mínimos quadrados:

void lms(double *x, double *y, int n, double *m, double *b)
{
    int i;
    double sumYX = 0.;
    double sumX = 0.;
    double sumY = 0.;
    double sumX2 = 0.;
    double sum2X = 0.;
    for(i = 0; i < n; i++) {
        sumYX += x[i] * y[i];
        sumX += x[i];
        sumY += y[i];
        sumX2 += x[i] * x[i];
    }
    sum2X = sumX * sumX;
    
    *m = (sumYX - (sumX * sumY) / (double)n) / (sumX2 - sum2X / (double)n);
    *b = sumY / (double)n - *m * sumX / (double)n;
}

O algoritmo é exatamente o mesmo utilizado na resposta dele. O laço for inicial é feito para calcular todos as somatórias e posteriormente seus resultados são utilizados para achar m e b. Acredito que isso seja facilmente traduzível para sua linguagem de preferência. :)

Para chamar o método:

int main(void) 
{
    double x[] = {1., 2., 3., 4.};
    double y[] = {10., 20., 30., 40.};
    double m, b;

    lms(x, y, 4, &m, &b);
    
    printf("m: %f \n b: %f", m, b);

return 0;
}
Fonte Link
Carregando