Métodos Numéricos -HP PRIME-UNH-Gabriel David Quispe Sanes

UNIVERSIDAD NACIONAL DE HUANCAVELICA

FACULTAD DE CIENCIAS DE INGENIERÍA ESCUELA ACADÉMICA PROFESIONAL DE INGENIERÍA CIVIL-HVCA

MANUAL DE PACK DE APLICACIONES PROGRAMADOS PROGRAMADOS EN LENGUAJE HP PPL PPL ____________________________ ______________ ____________________________ ____________________________ ___________________ _____ : Métodos Numéricos Aplicados a la Ingeniería. CATEDRÁTICO : MSc. Ing. Iván Arturo, Ayala Bizarro. : Gabriel David, Quispe Sanes. ESTUDIANTE : IV CICLO ____________________________ ______________ ____________________________ ____________________________ ____________________ ______ 







CÁTEDRA

A MI MADRE: Por haberme dado la vida que es el tesoro más preciado del mundo.

A MI PADRE: Por su invaluable apoyo, ternura y cariño que siempre me ha ofrecido.

AL MSc. Ing. IVÁN A. AYALA BIZARRO: Por su exigencia académica y sus recomendaciones motivadoras y reflexivas.

A MIS AMIGOS DE SIEMPRE: Por sus buenos deseos y alientos para seguir adelante

A MI MADRE: Por haberme dado la vida que es el tesoro más preciado del mundo.

A MI PADRE: Por su invaluable apoyo, ternura y cariño que siempre me ha ofrecido.

AL MSc. Ing. IVÁN A. AYALA BIZARRO: Por su exigencia académica y sus recomendaciones motivadoras y reflexivas.

A MIS AMIGOS DE SIEMPRE: Por sus buenos deseos y alientos para seguir adelante

INTRODUCCIÓN El presente trabajo titulado : “MANUAL DE PACK DE APLICACIONES HP PPL” , es el

PROGRAMADOS EN LENGUAJE

resultado de esfuerzo, sacrificio y empeño por parte de mi persona.

El resultado de esta iniciativa mi satisface mucho personalmente, gracias al empeño y la exigencia de los grandes maestros, me refiero especialmente al MSc. Ing. Iván Arturo, Ayala Bizarro, quien actualmente es docente de la cátedra de MÉTODOS NUMÉRICOS APLICADOS A LA INGENIERÍA en la Escuela Académica Profesional de Ingeniería Civil-Hvca, de la UNIVERSIDAD NACIONAL DE HUANCAVELICA. HUANCAVELICA. La iniciativa de este trabajo se dio desde el primer día del desarrollo de la cátedra Métodos Numéricos Aplicados a la Ingeniería, que se vino trabajando conforme se iba avanzando el silabo de dicha cátedra ya mencionada. El trabajo consistía en programar los temas que engloban dicha cátedra como por ejemplo: Solución de ecuaciones no lineales utilizando los diferentes métodos de solución aplicados a la ingeniería Civil (Cálculo de fricción en tuberías mediante mediante la ecuación de de Colebrook-White, cálculo de tirante en canales abiertos, cálculo de cuantía de aceros de refuerzo en vigas, etc.), solución de sistema de ecuaciones no lineales, cálculo de autovalores y autovectores mediante el método de Jacobi, regresiones(Lineal, logarítmica, exponencial, potencial, polinómica y múltiple lineal ) mediante solución matricial, integración numérica(Método trapezoidal, Simpson 1/3, Simpson 3/8 y la cuadratura de Gauss), etc. El lenguaje de programación programación a que mi refiero es HP-PPL, y este lenguaje es propio de la calculadora HP PRIME. Programar en este lenguaje nos facilita para desarrollar los exámenes parciales con mucha facilidad, aparte de que te ayuda a entender como los grandes software hacen grandes cálculos exactos y precisos. La base para grandes logros en programación digital en diferentes lenguajes es conocer antes que nada las diferentes funciones que cumple cada comando de programación como las condicionales (IF, ELSE, etc.), bucles(FOR, WHILE,etc), así así como las funciones de PRINT, INPUT, MSGBOX, MSGBOX, BREAK, etc. Y además además si queremos tener éxito en el campo de programación un requisito indispensable que todo programador debe manejar es la imaginación y la habilidad para dar solución a ciertos inconvenientes que se pueden presentar en el acto y finalmente un requisito que complementa para el logro eficaz en programación es el orden al momento de realizarla. Finalmente, el presente trabajo es una guía de las 6 aplicaciones que fue realizado por mi propia persona, persona , sin más que decir y agradeciendo a cada de los lectores presento las manualidades de las aplicaciones en lenguaje HP-PPL. Huancavelica 14 de febrero del 2017

G-COLEBROOK WHITE

Icono de la aplicación

AUTOR: QUISPE SANES, Gabriel David.

INTRODUCCIÓN El presente trabajo de programación en lenguaje HP PPL cuyo resultado final es una aplicación denominado “G-COLEBROOK WHITE ”, es el resultado de un gran esfuerzo y sacrifico, de mi propia persona.

Esta aplicación que presento a continuación es de mucha importancia en el campo de la Ingeniería Civil, específicamente para el cálculo de fricción en tuberías mediante la Ecuación de Manning. El cálculo de fricción en tuberías es esencial para un diseño eficaz de una red de tuberías. La ecuación de ecuación de Colebrook-White, es un tipo de ecuación no lineal, que analíticamente su solución es de mucha laboriosidad, razón por la cual, un método apropiado y eficiente es el Método Numérico, que soluciona la ecuación no lineal mediante las iteraciones y con un margen de error que debe ser controlado, de acuerdo a la precisión buscada. Hoy en día, la mayoría de las aplicaciones iterativas para la solución de ecuaciones no lineales básicamente lo desarrollan mediante el método de Newton Raphson, pero vale decir que los otros métodos tales como Método de Punto fijo, Bisección, Secante, etc., también tienen la misma precisión que el Método de Newton Raphson, a diferencia de que unos otros convergen lento o rápido. Finalmente está aplicación es aprovechable para la catedra de Métodos Numéricos, específicamente para el cálculo de fricción en tuberías mediante la solución numérica.

EL AUTOR

INSTALACION DE LA APLICACIÓN G-COLEBROOK WHITE 

Abrir la carpeta del archivo de la aplicación.



Como primer paso abrir el Kit de conectividad de HP PRIME:



Finalmente descomprimir el archivo y arrastrar encima del icono del emulador de la HP PRIME en Kit de conectividad, quedándonos de esta manera instalada la aplicación ( G-COLEBROOK WHITE).



Para ejecutar la aplicación damos click en icono de GUARDAR, para enviar la aplicación a la HP PRIME (Emulador o físico) y automáticamente se nos aparecerá en la calculadora.

MANEJO DE APLICACIÓN G-COLEBROOK WHITE 











Para el manejo de la aplicación es importante saber teóricamente sobre el cálculo de fricción en tuberías, mediante la Ecuación de Colebrook-White. Una vez ya instalada G-COLEBROOK WHITE, damos click encima del icono de la aplicación, y de pronto nos aparecerá la interfaz de la aplicación.

Damos click en [Presione aquí para continuar], para realizar los cálculos respectivos, y automáticamente nos saldrá una nueva ventana para introducir los datos necesarios para dicho cálculo. Los datos que aparecerán en el cuadro de entrada saldrá por defecto, para realizar la operación simplemente reemplazar por encima de los datos que vienen por defecto y que el programa controla tranquilamente con los nuevos datos introducidos. En la parte inferior apareceré una ventana para elegir el método deseado, esto viene por defecto con las iniciales “SALIR”, para que pueda realizar la operación es necesario elegir un método y pulsar “OK” , si en caso contrario no elige el método y pulsa “OK”, automáticamente el programa se cerrará sin realizar ningún cálculo. Antes de continuar con el ingreso de datos, recordaremos algunas fórmulas para el cálculo de fricción en tuberías. Ecuación de Colebrook White:

1  −2 ∙ቆ   2.51 ቇ 3.7∙  ∙ ඥ  ඥ  Donde:   Coeficiente de fricción.   Coeficiente de rugosidad.   Número de Reynolds.   Diámetro de la sección transversal de la tubería.

Para su cálculo respectivo se hace de la siguiente manera:

  ඥ 1  →   

  2.51∙൰   −2∙൬3.7∙    2.51∙൰ ሺሻ  2 ∙൬3.7∙  Para el Método de Newton Raphson la derivada de la función es:    3.7∙   2.51 2∙ ሺሻ  1  ሺ ∙ሻ∙ሺ10ሻ 



 

Para introducir los datos respectivos, en la parte inferior de la pantalla, una vez haya dado click en el primer casilla de entrada de datos, aparece un mensaje de que es lo debes ingresar precisamente en ese primer cuadro, y así sucesivamente para cada casilla de entrada de datos aparece el mensaje respectivo. Es importante tener en cuenta que para el Método de Punto fijo, es importante dar valores iniciales precisos, caso contrario el método no va converger.

Pondremos a prueba con el método de Punto fijo para el cálculo de fricción respectiva para los datos iniciales. Una vez ingresado lo datos respectivos y seleccionado el Método de punto fijo, damos click en “OK” y automáticamente si nos abrirá el editor de matrices (EDITMAT), donde aparece las iteraciones respectivas.









Damos click en “OK” y automáticamente nos saldrá una ventana de editor de textos, con los datos iniciales y el resultado final:

Para realizar una nueva operación ya sea utilizando un nuevo método simplemente damos clik en “Esc” de la tecla de la HP PRIME y automáticamente nos retornara al inicio. Finalmente pondremos a prueba con el Método de la Secante y lógico que tiene que llegar a los mismos resultados finales. Una vez seleccionado el Método de la Secante, damos click en “OK”, y de pronto si nos abrirá una nueva ventana con un encabezado de texto “INGRESE DATOS ADICIONALES ”, esto debido a que el método a utilizar (Secante), es un método abierto, pero calcula la raíz de la ecuación desde 2 intervalos iniciales cualquieras.







Una vez ingresado los datos adicionales (Nos va pedir para el Método de la Secante y Bisección), pulsamos “OK”, y si nos abrirá el editor de matrices, que es específicamente las iteraciones respectivas.

Para visualizar los datos iniciales y el resultado final, simplemente damos click en “OK”, y aparece la ventana de datos y resultado final. Si comparamos el resultado de ambos métodos puestos a prueba, notamos que la precisión es casi lo mismo, tiene una variación mínima para un error máximo permitido.





Y así sucesivamente puedes comprobar con cualquiera de los métodos siguiendo la misma secuencia de la explicación realizada. Para salir de la aplicación simplemente pulsamos “SALIR” y “OK”.

Próximos lanzamientos…….lenguaje de programación Python!!!

Cálculo de fricción en tuberías en lenguaje de programación en Shell de Python al 80%

GRACIAS…Seguiremos trabajando……

G-ACEREF



INTRODUCCIÓN El presente trabajo de programación en lenguaje HP PPL cuyo resultado final es una aplicación denominado “G-ACEREF”, es el resultado de perseverancia, lealtad y esfuerzo de mi persona.

Esta aplicación que presento a continuación es de mucha importancia en el campo de la Ingeniería Civil, sobre todo en las edificaciones, para calcular la cantidad de aceros de refuerzo que necesita una viga para soportar con eficiencia a la carga a la que es sometido, esto para no generar fallas estructurales. La aplicación que os presento realiza los cálculos usando la solución numérica, mediante los Métodos de Newton Raphson, Punto fijo y Secante, que son métodos para dar solución a una ecuación no lineal. En el campo de la Ingeniería Civil, es importante saber cifras exactas para un determinado proyecto u obra de construcción. Para esto, necesariamente se tiene que corroborar mediante programas, para el cálculo exacto y optimo, para un determinado condición que se requiere. La aplicación que presento, tiene justamente esa función de calcular la cantidad de aceros de refuerzo mediante el uso de Métodos Numéricos, que necesita una viga para soportar una carga, sin la ayuda de esta aplicación u otras aplicaciones que cumplen esta misma función, no se podría hacer el cálculo de la cantidad de aceros de refuerzo que necesita dicha viga. Esto es lo que nos diferencia a los Ingenieros Civiles frente a un maestro de obra, de utilizar los softwares para los cálculos exactos del material u otro insumo, mientras un maestro de obra, frente a esta situación, lo que haría es hacer un tanteo de la cantidad de aceros de refuerzo que necesita cierta viga, esto generaría un peligro de colapso con el tiempo o una inversión económica alzada. Es por tal razón, los Ingenieros Civiles sin excepción alguna debemos de por lo menos saber cómo los grandes softwares hacen los cálculos respectivos, y sin caso a alguno de nosotros nos nace esta iniciativa, si se puede lograr, con bastante esfuerzo decisión y sacrifico. Finalmente está aplicación es aprovechable para la catedra de Métodos Numéricos, específicamente para el cálculo de fricción en tuberías mediante la solución numérica.

EL AUTOR

INSTALACION DE LA APLICACIÓN G-ACEREF 









Finalmente descomprimir el archivo y arrastrar encima del icono del emulador de la HP PRIME en Kit de conectividad, quedándonos de esta manera instalada la aplicación ( G-ACEREF ).


MANEJO DE APLICACIÓN G-ACEREF 









Para el manejo de la aplicación es importante saber teóricamente sobre el cálculo de cantidad de aceros de refuerzo de una viga puesto a carga. Una vez ya instalada G-ACEREF , damos click encima del icono de la aplicación, y de pronto nos aparecerá la interfaz de la aplicación.

Damos click en [Presione aquí para continuar], para realizar los cálculos respectivos, y automáticamente nos saldrá una nueva ventana para introducir los datos necesarios para dicho cálculo. Los datos que aparecerán en el cuadro de entrada saldrá por defecto, para realizar la operación simplemente reemplazar por encima de los datos que vienen por defecto y que el programa controla tranquilamente con los nuevos datos introducidos. Antes de continuar con el ingreso de datos, recordaremos algunas fórmulas para el cálculo de cantidad de aceros de refuerzo mediante Métodos Numéricos.

Viga Longitudinal

Sección transversal de la viga

Del gráfico: Peralte de la viga. Peralte efectivo de la viga. Base de la viga. Longitud de la viga. Formulas y ecuaciones utilizadas:

ℎ    ሺሻ  0.59 ∙  ∙  −   ∙ ሺ ሻ Derivada para el Método de Newton Raphson: ሺሻ  1.18∙  ∙  −1 De la ecuación general:   ∙    ∙   ℎ −0.05    ∙ ∙   ∙    4     Donde: : Momento último.  ∶ Momento nominal.  ∶ Base de la viga.  ∶ Peralte efectivo de la viga.  ∶ Base de la viga.  ∶ Área del acero a ocupar.  ∶ Área de la sección transversal de la viga.  ∶ Cuantía de acero.  ∶ Cantidad de aceros requeridos. 



Para introducir los datos respectivos, en la parte inferior de la pantalla, una vez haya dado click en el primer casilla de entrada de datos, aparece un mensaje de que es lo debes ingresar precisamente en ese primer cuadro, y así sucesivamente para cada casilla de entrada de datos aparece el mensaje respectivo. Es importante tener en cuenta que para el Método de Punto fijo, es importante dar valores iniciales precisos, caso contrario el método no convergerá.







Pondremos a prueba con el método de Punto fijo para el cálculo de la cantidad de aceros de refuerzo con los datos iniciales. Una vez ingresado lo datos respectivos y seleccionado el Método de punto fijo, damos click en “OK” y automáticamente si nos abrirá el editor de matrices (EDITMAT), donde aparece las iteraciones respectivas.

Damos click en “OK” y automáticamente nos saldrá una ventana de editor de textos, donde aparecerá los resultados parciales y el resultado final.







Para realizar una nueva operación ya sea utilizando un nuevo método simplemente damos clik en “Esc” de la tecla de la HP PRIME y automáticamente nos retornara al inicio. Finalmente pondremos a prueba con el Método de la Secante y lógico que tiene que llegar a los mismos resultados finales. Una vez seleccionado el Método de la Secante, en esta parte debemos tener en cuenta en la casilla de entrada de dato de la parte inferior derecha último, donde aparece las iniciales que específicamente se refiere al segundo intervalo, siendo el primero , para este método.





1ሺሻ

Pulsamos “OK”, y si nos abrirá el editor de matrices, que es específicamente las iteraciones respectivas.

 

 

Para visualizar el resultado final, simplemente damos click en “OK”, y aparece la ventana de editor de textos. Si comparamos el resultado de ambos métodos puestos a prueba, notamos que el resultado son los mismos, quiere decir que ambos métodos son precisos.

Y el mismo proceso se realiza para hacer el cálculo mediante el método de punto fijo. Finalmente para salir de la aplicación simplemente pulsamos “SALIR” y “OK”.


GCANALES



INTRODUCCIÓN El presente trabajo de programación en lenguaje HP PPL cuyo resultado final es una aplicación denominado “GCANALES”, es una de las aplicaciones sencillas que fue trabajado por mi propia persona, es así que voy

trabajando por mejorar en el diseño de interfaz de la gráfica y el manejo de las matrices que es de mucha de importancia en el campo de la Ingeniería, como sabemos los grandes programas como SAP 200 utilizan matrices para solucionar. Ya pronto estaré lanzando una aplicación de Análisis matricial de estructuras mediante las rigideces, donde que me muestre sus matrices de rigidez de cada barra, paso a paso, así como factorización de Cholesky de una matriz de nxn y que ésta a la vez sea aplicado en la solución de un sistema de ecuaciones lineales y también en el tema de regresiones polinómicas y múltiples. La aplicación tiene importancia en el Diseño de Canales de tipo Rectangular, Triangular y Trapezoidal, en la cátedra de Métodos Numéricos Aplicados a la Ingeniería, en el cual es muy importante controlar el margen de error en el cálculo del tirante de los canales de secciones ya mencionadas, utilizando los 4 métodos básicos para la solución de ecuaciones no lineales en el campo de la Ingeniería. Si bien es cierto en la actualidad se cuenta con muchas aplicaciones para el diseño de canales de diferentes secciones, estas aplicaciones en su mayoría solamente utilizan un solo método (Método de Newton Raphson) para calcular el tirante del canal, pero esto no nos limita demostrar que también con los otros métodos se llega a los mismos resultados. Finalmente está aplicación es aprovechable para la catedra de Métodos Numéricos, donde se usa diferentes métodos para dar solución a una serie de problemas.

EL AUTOR

INSTALACION DE LA APLICACIÓN GCANALES 







Finalmente descomprimir el archivo y arrastrar encima del icono del emulador de la HP PRIME en Kit de conectividad, quedándonos de esta manera instalada la aplicación ( GCANALES ).




MANEJO DE APLICACIÓN GCANALES 











Para el manejo de la aplicación es importante saber teóricamente sobre diseño de canales rectangulares triangulares y trapezoidales. Una vez ya instalada GCANALES damos click encima del icono de la aplicación, y de pronto nos aparecerá la interfaz de la aplicación.

Damos click en [Presione aquí para continuar] , para realizar los cálculos respectivos y necesarios, y automáticamente nos saldrá una nueva ventana para introducir los datos respectivos. Los datos que aparecerán en el cuadro de entrada saldrá por defecto, para realizar la operación simplemente reemplazar por encima de los datos que vienen por defecto y que el programa controla tranquilamente con los nuevos datos introducidos. En la parte inferior apareceré una ventana para elegir el método deseado, esto viene por defecto con las iniciales “SALIR”, para que pueda realizar la operación es necesario elegir un método y pulsar “OK” , si en caso contrario no elige el método y pulsa “OK”, automáticamente el programa se cerrará sin realizar ningún cálculo. Antes de introducir los datos al programa os mencionaré las formulas necesarias para el cálculo del tirante de las secciones rectangular, triangular y trapezoidal.

Área del canal:

 ∙ ሺ  ሻ    2  ∙ Perímetro mojado:     ∙ቆට 1   ට 1 ቇ Ecuación de continuidad:    ∙ Ecuación de Manning:   1 ∙ ൗ  ∙ൗ  ∙ ∙  ൗ  ∙  ൗ   ቆ∙ൗ ቇ   ∙  ∙  Sustituimos:   ቀ ൗ  ቁ      ∙  ሺሻ   ∙  −  Para el Método de Newton Raphson: ሺሺሻሻ  ′ሺሻ  2∙ ∙ ∙   3 ∙ ∙  ∙  −     EN FORMA GENERAL:

A partir del diseño del proceso de cálculo de un canal trapezoidal, se minimiza las fórmulas para el cálculo de tirante de secciones rectangular y trapezoidal, haciendo de la siguiente forma:  

    0 Para la sección triangular:   0 Para la sección rectangular:

Para el cálculo del resultado final, aplicamos las siguientes formulas en forma general para las secciones mencionadas:     

   ∙ሺ+ሻ  ∙ Perímetro mojado:     ∙ቀඥ 1    ඥ 1  ቁ Radio hidráulico:   ൗ Velocidad:   ൗ Espejo de agua :    ሺ   ሻ∙ Área hidráulica:







Pondremos a prueba con el método de Punto fijo para un canal de sección trapezoidal con taludes iguales, si en algún ejercicio nos pida calcular las dimensiones para un canal triangular o rectangular simplemente pondremos 0 a algunos datos como talud para el caso de una sección rectangular o b=0 para el caso de una sección triangular. En la entrada de datos aparece los parámetros a y b, los cuales sirven si es que vamos a utilizar los métodos de BISECCION Y SECANTE , caso contrario no es necesario dar valores precisos.

Damos click en “OK” y automáticamente nos saldrá otra ventana de las iteraciones:







Para ver los resultados finales simplemente damos click en “Enter” y nos aparecerá una ventana de resultados finales.

Haciendo una comparación con el programa “HCANALES” , los resultados nos da igual y con la misma precisión.

Para realizar una nueva operación ya sea utilizando un nuevo método simplemente damos clik en tecla de la HP PRIME y automáticamente nos retornara al inicio.

“Esc” de la



Finalmente pondremos a prueba con el Método de la Secante y lógico que tiene que llegar a los mismos resultados finales.



Para salir de la aplicación simplemente pulsamos “SALIR” y “OK”.

Muy pronto nuevo lanzamiento……..

Próximo proyecto….lenguaje de programación Python..!!!

Diseño de canales de secciones rectangular, triangular y trapezoidal en lenguaje de programación en Python al 70%


G-VECTORES APLICATION


AUTOR: QUISPE SANES, Gabriel David. CONSULTAS: 

CEL: 982925212



CORREO: [email protected]



FACEBOOK: Gabriel QS

INTRODUCCIÓN El presente trabajo de programación en lenguaje HP PPL cuyo resultado final es una aplicación denominado “GVECTORES APLICATION”, es una de las aplicaciones sencillas que fue trabajado por mi propia persona, es así que voy trabajando por mejorar en el diseño de interfaz de la gráfica y el manejo de las matrices que es de mucha de importancia en el campo de la Ingeniería Civil. La aplicación tiene como fin calcular los vectores y valores propios de una matriz simétrica de NxN mediante la solución numérica por el Método de Jacobi. Más que una aplicación es una ayuda para el cálculo de vectores y valores propios de una matriz en la cátedra de Métodos Numéricos Aplicado a la Ingeniería. Calcular los vectores y valores propios de una matriz mediante la solución numérica no es imposible pero si trabajoso sin la ayuda de un software como el Excel. Esta aplicación que presento no es menos ni más que el Excel, pero con algunas ventajas para el cálculo de vectores y valores propios de una matriz llegando finalmente a la misma respuesta y precisión requerida. La aplicación es de uso muy sencillo. Para empezar simplemente introducir la matriz (Simétrica y cuadrada) y dejar que la aplicación se encargue por la solución y precisión buscada. Finalmente está aplicación es aprovechable para la catedra de Métodos Numéricos, donde se usa diferentes métodos para dar solución a una serie de problemas que se presentan en el campo de la Ingeniería.

EL AUTOR

INSTALACIÓN DE LA APLICACIÓN G-VECTORES APLICATION 

Abrir la carpeta del archivo de la aplicación:








Finalmente descomprimir el archivo y arrastrar encima del icono del emulador de la HP PRIME en Kit de conectividad, quedándonos de esta manera instalada la aplicación ( G-VECTORES APLICATION).


MANEJO DE APLICACIÓN G-VECTORES APLICATION 







Para el manejo de la aplicación es importante saber teóricamente sobre el cálculo de auto valores y auto vectores de una matriz por el método de Jacobi Una vez ya instalada G-VECTORES APLICATION damos click encima del icono de la aplicación, y de pronto nos aparecerá la interfaz de presentación de la aplicación.

Damos click en [Presione aquí para continuar] , para realizar los cálculos respectivos y necesarios, y automáticamente nos saldrá una nueva ventana para introducir los datos iniciales y únicos. En la ventana nos parecerá dos cuadros de entrada de datos renombrados “NxN” y “C” y en la parte inferior de la ventana sale un texto (“Ingrese el tamaño de la matriz simétrica ”) cuando pulsamos en una de las entradas, y es ahí donde específica a que se refiere las iniciales “NxN” y “C” y en particular en todo el trabajo que he

desarrollado programando el lenguaje HP PPL siempre especifico en la parte inferior a que se refiere específicamente las iniciales de las entradas para no tener inconvenientes a la hora de ingresar los datos.



Una vez introducido los datos correctamente, damos click en “OK”, y automáticamente nos saldrá una nueva ventana de editor de matrices ( EDITMAT), ya redimensionados (En nuestro caso 4x4) y rellenados de 0 por defecto.









En la tarde superior aparece un texto de encabezado ( “MATRIZ INICIAL” ) y es en esa parte donde especifica en que campo nos encontramos (Campo de datos, resultados).

Una vez que se nos aparezca la ventana de editor de matrices( EDITMAT), simplemente introducimos los datos de la matriz del cual queremos calcular sus valores y vectores propios y pulsamos “OK”

Una vez pulsado “OK” nos parecerá un mensaje ( “CICLO 1”), el cual nos indica que estamos iniciando el proceso iterativo desde el ciclo 1 hasta el ciclo que introducimos (en nuestro caso 3), vale decir que realizara 3 ciclos para calcular los valores y vectores propios de la matriz. Para introducir la cantidad de ciclos dependerá de la precisión buscada, es así que a más ciclos mayor precisión y exactitud. Pulsamos “OK”.

 



Una vez pulsado “OK” se nos aparecerá un nuevo mensaje de texto ( “Haciendo 0 a los coeficientes a21” ). Pulsamos “OK”.

Después de pulsar “OK” se nos parecerá la primera matriz de resultados, específicamente la matriz P_n (donde n es el contador).







Pulsamos “OK” y automáticamente si nos abrirá un nuevo editor de matrices con encabezado de “MATRIZ DE VALORES ”, que se refiere específicamente a los pre-valores propios de la matriz.

Damos click en “OK” y nos parecerá una nueva ventana de editor de matrices con encabezado de “MATRIZ DE VECTORES” , que se refiere en específico a sub-vectores propios de la matriz. Y así sucesivamente dar click en “OK” y el programa se encarga de calcular los pre-resultados y al final el resultado buscado.





Y continuación nos vamos a la parte final del “CICLO 1”, para ver cómo se va acercando al resultado en el ciclo mencionado. Esto específicamente se ve cuando el mensaje sale ( “Haciendo 0 a los coeficientes a43” ), porque es el último elemento para hacer 0 y cumplir el ciclo.

Matriz de valores propios en el CICLO 1: Como observamos las posiciones a21, a31, a32, a41, a42, a43 (simétrica), tienden a 0 en el CICLO 1.

Matriz de vectores propios en el CICLO 1 

Y finalmente nos vamos al “CICLO 3”, para ver el resultado final y para verificar la precisión a mayor ciclo que realizamos.

GINTEGRATION NUMERIC


AUTOR: QUISPE SANES, Gabriel David. CONSULTAS: 

CEL: 982925212



CORREO: [email protected]



FACEBOOK: Gabriel QS

INTRODUCCIÓN El presente trabajo de programación en lenguaje HP PPL cuyo resultado final es una aplicación denominado “GINTEGRATION NUMERIC ”, es una de las aplicaciones sencillas y a la vez de mucha importancia que fue trabajado por mi propia persona en su totalidad, habiendo tenido la necesidad de contar con un programa que calcule la integral numérica de una función de lo más simple a lo más complejo mediante la solución numérica. La aplicación tiene como finalidad calcular la integral definida de una función algebraica, logarítmica, exponencial, trigonométrica, etc. Mediante los Métodos Numéricos tales como la Regla de Trapecio, Simpson 1/3, Simpson 3/8 y finalmente la Cuadratura de Gauss, este último es uno de los métodos de gran importancia y de mucha precisión. En el campo de la Ingeniería Civil, las integrales se presentan frecuentemente, en algunos casos difíciles de operar o casi imposible desarrollar analíticamente, es por tal razón se recorre a los Métodos Numéricos para dar la solución respetiva con un margen de error controlado. Finalmente está aplicación es aprovechable para la catedra de Métodos Numéricos Aplicados a la Ingeniería, donde se usa diferentes métodos para dar solución a una serie de problemas que se presentan en el campo de la Ingeniería Civil.

EL AUTOR

INSTALACIÓN DE LA APLICACIÓN G-INTEGRATION NUMERIC 

Abrir la carpeta del archivo de la aplicación:








Finalmente descomprimir el archivo y arrastrar encima del icono del emulador de la HP PRIME en Kit de conectividad, quedándonos de esta manera instalada la aplicación ( G-INTEGRATION NUMERIC).


MANEJO DE APLICACIÓN G-INTEGRATION NUMERIC 









Para el manejo de la aplicación es importante conocer teóricamente sobre el cálculo integral analítica y numérica de una función, y en específico conocer los métodos de integración numérica tales como Método de trapecio, Simpson 1/3, Simpson 3/8 y la Cuadratura de Gauss Una vez ya instalada G-INTEGRATION NUMERIC damos click encima del icono de la aplicación, y de pronto nos aparecerá la interfaz de presentación de la aplicación.

Damos click en [Presione aquí para continuar] , para realizar los cálculos respectivos y necesarios, y automáticamente nos saldrá una nueva ventana para elegir el método respectivo con el cual deseamos calcular la integral de la función. En la ventana nos parecerá 6 sub-menús para la elección respectiva del método así como para ver el autor del programa y finalmente para salir del programa una vez realidad la operación respectiva. Seleccionamos uno de los sub-menús dando click encima del ítem deseado (en nuestro caso empezaremos con el ítem 1, Método de trapecio).









Una vez seleccionado el ítem, si nos abrirá una nueva ventana para introducir los datos respectivos tales como la función, límite superior, inferior, las particiones que deseamos realizar. Pondremos a prueba con un integral: que analíticamente hacer el cálculo es muy trabajoso. Utilizando el programa interno de la calculadora HP PRIME se calculó un resultado de: 52.0499711038

  ξ 1 

Una vez introducido los datos, pulsamos “OK”, y automáticamente si nos abrirá una ventana de editor de matrices (EDIMAT) con un encabezado “MATRIZ RESULTADO”, en esta parte aparece las variables “x” así como “P” y f(x). Para que tengan en claro que es lo que hace el programa y como lo calcula la integral de una función y aparte de eso a que se refiere las variables “x” y “P” en todos los métodos como ya verán. Realizare

un poco de teoría en esta parte para que les quede muy claro.

Nº

ℎ   −               x

f(x)

P

f(X)*P

0

a

f(a)

0.5

0.5*f(a)

1

a+h

f(a+h)

1

f(a+h)

2

a+2h

f(a+2h)

1

f(a+2h)

3

a+3h

f(a+3h)

1

f(a+3h)

4

a+4h

f(a+4h)

1

f(a+4h)

:

:

:

:

:

n

b

f(b)

0.5

0.5*f(b)

෍  ∗

    ℎ ∗ ෍ሺሻ ∗  Y específicamente el cuadro anterior nos aparece en la siguiente ventana que visualizaremos.







Una vez visualizado el cuadro de pre-resultados, damos click en “OK” y se nos parecerá un visualizador de textos con un título de “RESULTADOS” , donde aparece la sumatoria total y resultado final (Integral).

Una vez ya calculado la integral definida simplemente pulsamos “Esc” de la tecla de la HP PRIME para retornar al menú inicial y si es queremos realizar una nueva operación simplemente pulsamos en uno de los ítems. Para el caso de Simpson 1/3 y 3/8 se sigue los mismos procedimientos ya explicados a diferencia de la teoría de cada uno de los métodos, a continuación detallaré la teoría de Simpson 1/3 y 3/8 y calcularemos la misma función y compara la precisión de cada uno de los métodos.

 

El pre-requisito para aplicar Simpson 1/3 es que “n” tiene que ser mayor a 4. El pre-requisito para aplicar Simpson 3/8 es que “n” tiene que ser mayor a 6 y múltiplo de 3. SIMPSON 1/3

SIMPSON 3/8

Nº 0

x a

f(x) f(a)

P 1

f(X)*P f(a)

P 1

f(X)*P f(a)

1

a+h

f(a+h)

4

4*f(a+h)

3

3*f(a+h)

2

a+2h

f(a+2h)

2

2*f(a+2h)

3

3*f(a+2h)

3

a+3h

f(a+3h)

4

4*f(a+3h)

2

2*f(a+3h)

4

a+4h

f(a+4h)

2

2*f(a+4h)

3

3*f(a+4h)

:

:

:

:

:

:

:

n

b

f(b)

1

f(b)

1

f(b)

෍  ∗

    13  13 ෍ሺሻ∗     38  38 ෍ሺሻ∗



Procedemos a calcular con los dos métodos restantes de Simpson 1/3 y 3/8:

SIMPSON 1/3

෍  ∗ 

SIMPSON 3/8



Como vemos el método que se acerca más rápido hasta este instante es de Simpson 3/8, pero no está llegando a la precisión requerida con n=20. Finalmente comprobaremos con el Método de cuadratura de Gauss. Para ello regresamos al menú principal y damos click en el ítem 4. Una vez seleccionado el ítem, si nos aparecerá un cuadro de entrada de datos de la función, límite superior e inferior y numero de particiones(n). Si bien es cierto, este método necesita de constantes n para realizar los cálculos respectivos, el programa ya viene introducido tales constantes que varía desde n=1 hasta n=10. Se requiere calcular la integral con la mayor precisión posible simplemente para n=10. Una vez introducido los datos pulsamos “OK”.



Después de pulsar “OK” se nos parecerá la primera matriz de constantes de “n”(



Pulsamos “OK” y si nos aparecerá la matriz de resultados de las iteraciones respectivas en nuestro caso 10











, ):



En la columna 4 de esta matriz aparece “RESULT_ITER” que se refiere específicamente a la parte de sumatoria que aparece en la siguiente formula de Cuadratura de Gauss para el cálculo de integral definido.

   − ∫ ሺሻ  ൬ 2 ൰∗෍ ∗൬ −2  ∗     ൰ 2 =





Para calcular el resultado final simplemente pulsamos “OK” y se nos parecerá el resultado final, que es ya casi exacto que la precisión buscada.

Finalmente para visualizar el autor de este programa propiamente dicho pulsar el ítem 5 y si en caso desea salir del programa simplemente pulsar por encima de “SALIR ”

PROGRAMAS EN EJECUCIÓN… 

G-REGRETION NUMERIC:



SOLUCION DE ECUACIONES DIFERENCIALES(Método de Runge Kutta):



ANÁLISIS MATRICIAL DE ESTRUCTURAS (A un 40%):

Métodos Numéricos -HP PRIME-UNH-Gabriel David Quispe Sanes

Recommend Documents