Método de Newton-Raphson con C++ Linux

Innumerables sistemas en la naturaleza pueden representarse mediante ecuaciones no lineales que vienen dadas por polinomios y la solución exacta, para muchos de ellos, está en encontrar la única raíz que tiene sentido físico para un polinomio de grado mayor que 2. En estos casos, hay que apelar a métodos numéricos y el de Newton-Raphson converge rápidamente. El código del programa es este:


#include <iostream>
#include <fstream>

using namespace std;

#include <cmath>

void evaluapoli(int n, double *f, double x, double &p); //utilizada por newton_raphson.h
void derivapoli(int n, double *f, double *fd);
void newtonraphson(int n, double *f, double *fd, double &x, int &contador);

int main(){

 ifstream label1 ("datos//datos.in");

 int n;

 cout << "El grado del polinomio y sus coeficientes están en datos/datos.in\n\n";

 //Introduce el grado del polinomio

 label1 >> n; cout << "grado = " << n << "\n\n";

 double *f, *fd; 

 f = new double [n+1];   // Arreglo para los coeficientes del polinomio

 fd = new double [n];    // Arreglo para los coeficientes de la derivada del polinomio

 // Se introducen desde datos.in y se imprimen en pantalla los coeficientes del polinomio

 cout << "Coeficientes del polinomio\n\n";

 for (int i = n; i > -1; i--) {  

  label1 >> f[i]; 

  cout << "F" << i << " = " << f[i] << "\n";

 }

 cout << "\n";

 derivapoli (n, f, fd);  // función que determina la primera derivada y
          // se imprimen sus coeficientes

 cout << "Coeficientes del polinomio primera derivada\n\n";

 for (int i = n - 1; i > -1; i--) {  

  cout << "Fd" << i << " = " << fd[i] << "\n";

 }

 cout << "\n";

 cout << "Introduzca el valor de inicio para la raiz: \n" << endl << "X = ";

 double x; 

 int contador = 0;

 // Se introduce el valor de x

 cin >> x;

 cout << endl;

 // Implementación del método de Newton-Raphson

 newtonraphson (n, f, fd, x, contador); 

 cout << "La raiz es x = " << x << " y se obtuvo en " << contador << " iteraciones\n";

 return 0;

}

void evaluapoli(int n, double *f, double x, double &p){

 int i; 

 double *a;

 a = new double [n];

 for (i = n; i > -1; i--) {  // Crea una copia de los elementos de f[i]
                             // porque el ciclo for inferior los destruye

  a[i] = f[i];                   

 } 

 for (i = n; i > 0; i--) {

  p = a[i] * x + a[i - 1];

  a[i - 1] = p;

 }

// P es el valor del polinomio en el punto x

}

void derivapoli(int n, double *f, double *fd){

 int i;

 for (i = n; i > 0; i--){

  fd[i-1] = i * f[i];

 }

}

void newtonraphson(int n, double *f, double *fd, double &x, int &contador){

 /* La función evaluapoli determina el valor del polinomio en el punto x.*/

 double p=0, y1, y2, verificador, z = 0;

 do {    evaluapoli (n, f, x, p);  // Determina el valor del polinomio en el punto

  y1 = p;

  evaluapoli (n-1, fd, x, p);  //Valor de la primera derivada en el punto

  y2 = p;

  x = x - y1/y2;

  verificador = fabs(x - z);

  z = x;

  contador += 1;  

 } while (verificador > 1e-10);

}

El programa principal llama a tres sub rutinas: una de ellas desarrolla el método de Newton-Raphson (que a su vez incluye evaluapoli). Aquí lo utilizo para encontrar las raíces de este polinomio:

P(x) = x³ – 2x² – 5x + 6

que obtuve con MONOMIO (las raíces son 1, -2, 3). El grado (3) y los coeficientes (1, -2, -5, 6) se introducen desde el archivo datos/datos.in y, por tanto, no hay que modificar el programa para incluir otro polinomio. Después de grabado el código anterior como newtonraphson.c++ en una carpeta llamada, por ejemplo, NEWTON_RAPHSON, nos ubicamos en ésta y la compilación se efectua con:

g++ newtonraphson.c++ -o newtonraphson

La ejecución del mismo (./newtonraphson [Enter]) produce la salida siguiente:

El grado del polinomio y sus coeficientes están en datos/datos.in

grado = 3

Coeficientes del polinomio

F3 = 1
F2 = -2
F1 = -5
F0 = 6

Coeficientes del polinomio primera derivada

Fd2 = 3
Fd1 = -4
Fd0 = -5

Introduzca el valor de inicio para la raiz:

X = -1000000

La raiz es x = -2 y se obtuvo en 37 iteraciones

Exageramos partiendo de un valor de prueba correspondiente a un millón (negativo) y convergió al valor de una de las tres posibles raíces en 37 iteraciones. De más está decir que mientras más cerca se está del valor verdadero la convergencia se obtiene en menos pasos.

Esta entrada fue publicada en Código C++. Guarda el enlace permanente.

9 respuestas a Método de Newton-Raphson con C++ Linux

  1. Diamar dijo:

    hola, tengo una duda respecto al archivo datos.in donde se guarda y que contiene exactamente…
    Eres tan amable y me dice por favor!

    • En el texto lo dice pero tal vez sea más ilustrativo con un ejemplo. Tiene esto:

      3
      1 -2 -5 6
      

      (también admite otras configuraciones) dentro de la carpeta datos. Todo eso lo puedes guardar, por ejemplo, en una carpeta llamada NEWTON-RAPHSON.

  2. Diamar dijo:

    Ah, se me olvido algo… el Blogg MONOMIO no aparece, es por eso que aun estoy confundida?

  3. Si usas el buscador del blog con monomio, te remitirá a varios artículos. El que buscas es éste:

    https://joseguerreroa.wordpress.com/2010/06/01/producto-de-monomios-en-c-linux/

    Ese artículo estaba en un Blog que tenía antes con Blogger y por eso es que el link apunta a “nada” (gracias a ti ya lo corregí). Por otra parte, con esa aplicación, “invente” un polinomio con base a sus raíces (1, -2, 3). Por eso sabía que una de ellas, para el método de Newton-Raphson en este artículo, era -2.

  4. jose dijo:

    hola disculpa me puedes explicar que significa esta parte
    (verificador > 1e-10)

    porfavor…!!!

    • Es la “tolerancia” o diferencia entre dos valores consecutivos (recuerda que es un método iterativo), es decir, el do…while realiza el proceso de cálculo mientras ésta sea mayor que 1.10-10. Mientras menor es la diferencia menor es el error o viceversa en la estimación de la posible raíz.

  5. araceli dijo:

    Para darle los parámetros que dices que lee desde el archivo datos/datos.in,
    ese archivo lo meto en la misma carpeta del programa y lo creo con cualquier editor?

Responder

Por favor, inicia sesión con uno de estos métodos para publicar tu comentario:

Logo de WordPress.com

Estás comentando usando tu cuenta de WordPress.com. Cerrar sesión / Cambiar )

Imagen de Twitter

Estás comentando usando tu cuenta de Twitter. Cerrar sesión / Cambiar )

Foto de Facebook

Estás comentando usando tu cuenta de Facebook. Cerrar sesión / Cambiar )

Google+ photo

Estás comentando usando tu cuenta de Google+. Cerrar sesión / Cambiar )

Conectando a %s