Estimación del Hidrograma Unitario en Hidrología
Estimación del Hidrograma Unitario en Hidrología
Abstract: Resumen:
Hydrology is an earth science with a very broad La hidrología es una ciencia que se encarga de
field of action because it is responsible for analizar el recurso más importante de la tierra que
analyzing the most important resource of the earth, es el agua. Por lo tanto, para muchos académicos
which is water. Therefore, for many academics it is es importante conocer diferentes factores de los
important to know different factors of water bodies cuerpos hídricos como lagos, lagunas, ríos, mares
such as lakes, lagoons, rivers, seas, etc., that allow etc., que permitan un mayor entendimiento de los
a better understanding of them. mismos.
Therefore, one of many studies carried out is Por consiguiente, uno de tantos estudios realizados
surface hydrology with the main objective of es la hidrología superficial con objetivo principal
calculating the runoff of rainwater through the de calcular el escurrimiento del agua de lluvia por
drainage network until it reaches the river network, la red de drenaje hasta alcanzar la red fluvial,
known more as runoff, which is produced by a conocido mas como escorrentía, que se produce
given precipitation. por una precipitación determinada.
And, all this is possible thanks to different Y, todo esto es posible gracias a diferentes métodos
analytical methods developed to analyze the analíticos desarrollados para analizar los
behaviors present in the runoff such as, for comportamientos presentes en las escorrentías
example, the unit hydrograph that also has several como lo es, por ejemplo, el hidrograma unitario
ways of being developed according to different que asimismo tiene varias formas de ser
authors who managed through convolution and desarrollado según diferentes autores que lograron
deconvolution processes to provide data that help a través de procesos de convolución y
the design of hydraulic structures or even for deconvolución brindar datos que ayudan al diseño
prevention of natural phenomena such as floods or de estructuras hidráulicas o, inclusive, para
overflows of rivers due to high rainfall. prevención de fenómenos naturales como crecidas
o desbordamientos de ríos por precipitaciones
Keywords:
elevadas.
Flow, hydrograph, runoff, linear regression, model.
Palabras claves: Caudal, hidrograma, escorrentía,
regresión lineal, modelo.
1
INTRODUCCIÓN OBJETIVOS
A través de los años, en la hidrología se ha estudiado el • Realizar un modelo con valores supuestos de
comportamiento de los diferentes cuerpos hídricos que caudales para poder evidenciar el desarrollo de la
existen, intentando aplicar modelos matemáticos que convolución y la desconvolución.
asemejen resultados similares a lo que ocurre en campo • Analizar las graficas generadas con los valores
pero siempre teniendo presente que esto es muy poco que representan el cambio de los caudales con
probable que pase, ya que en la mayoría de los casos de respecto el tiempo.
estudio de temáticas relacionadas con la naturaleza y los
recursos no se cuenta con situaciones precisas sino
suposiciones de lo que podría ocurrir debido a que MARCO TEÓRICO
siempre existirán agentes externos que no se podrán
calcular con exactitud. HIDROGRAMA UNITARIO.
Siendo así, que para la hidrología superficial es un reto
completo entender como procederán a actuar los ríos,
Es una herramienta utilizada en la hidrología para
cauces o demás cuerpos hídricos, a través del tiempo.
Ya que, por las precipitaciones que se lleven en la zona
realizar análisis, previsiones y estimaciones de
pueden tener consecuencias negativas para su entorno, inundaciones, todo esto a partir de hidrogramas de
ya sea que por la baja cantidad de lluvia se presenten escorrentía directa realizados en base a una
episodios de sequía, o en cambio, existan largos precipitación efectiva uniforme sobre una cuenca
periodos de lluvia generando desbordamiento de ríos en cierta cantidad de tiempo. Siendo así, que es uno
que afecten zonas rurales o urbanas cercanas. de los métodos más utilizados para crear
Debido a estas situaciones los científicos han tenido que predicciones de avenidas futuras generadas por
incursionar en encontrar diferentes métodos para aguaceros, investigaciones de la gestión de los
calcular la escorrentía que habrá en una cuenca recursos hidráulicos y/o efectos de las obras
especifica y esto es posible de realizar si se hacen
hidráulicas a largo plazo. [1]
estudios relaciones donde se realice un hidrograma de
Fue propuesto por Sherman en 1932 quien declaro
cierto periodo de tiempo y un hietograma de un exceso
de lluvia, porque con estos dos datos a través de que el método de hidrograma unitario se usaría
métodos matriciales, estadísticos y aritméticos se únicamente para la escorrentía superficial sin tener
pueden hacer un calculo aproximado de un hidrograma en cuenta la del agua subterránea, y tampoco
unitario, donde se asemeja una lluvia que aunque no se funciona para escorrentía generada por nieve o
repetirá bajo la ley de singularidad, se obtienen datos hielo. Además, aunque el modelo fue desarrollado
cercanos a lo que podría presentarse con cierto nivel de para cuencas grandes se ha descubierto que
precipitación en el área estudiada. funciona para cuencas pequeñas con al menos 0,5
Por último, es importante entender que aunque estos hectáreas hasta 25 km2.
métodos han sido estudiados a fondo para intentar ser Además, esta teoría conforma un modelo donde se
lo más veraces con relación a los verdaderos resultados
plantea que la cuenca es un sistema lineal e
que tienen las precipitaciones, no se puede asegurar
invariable a través del tiempo donde la entrada es
100% una predicción correcta, ya que como se dijo
previamente existen demasiados factores externos que la lluvia y la salida es el caudal de escorrentía
pueden cambiar la interpretación de los datos ya sea por directa. No obstante, esta teoría al ser aplicada
el clima, la zona geológica, los sedimentos en la cuenca, presenta situaciones en donde los comportamientos
presencia de plantas o árboles, entre otros. de las cuencas no son lineales ya que son afectados
por agentes externos, pero es utilizada por la
2
simplicidad y el conocimiento de los métodos
numéricos necesarios. [1]
PROCEDIMIENTO
3
Siendo así, que esto es empleado en la siguiente
ecuación:
𝑈 = (𝑁 − 𝑀 + 1)
En donde
𝑈 = 𝐻𝑖𝑑𝑟𝑜𝑔𝑟𝑎𝑚𝑎𝑠 𝑢𝑛𝑖𝑡𝑎𝑟𝑖𝑜𝑠
𝑀 = 𝑁° 𝑑𝑒 𝑝𝑢𝑙𝑠𝑜𝑠 𝑑𝑒 𝑒𝑥𝑐𝑒𝑠𝑜 𝑑𝑒 𝑙𝑙𝑢𝑣𝑖𝑎
𝑁 = 𝑁° 𝑑𝑒 𝑝𝑢𝑙𝑠𝑜𝑠 𝑑𝑒 𝑒𝑠𝑐𝑜𝑟𝑟𝑒𝑛𝑡í𝑎
MÉTODO MATRICIAL.
Ilustración 3. Procedimiento para realizar una regresión lineal
múltiple en Excel. Fuente: propia
Para este método, se trabajó en Excel en base a los
valores de la escorrentía directa con 15 pulsos y el
MÉTODO GAUSS JORDAN. hietograma con 4 pulsos, en donde cumpliendo la
formula:
Para este método se tiene el hidrograma de
escorrentía directa con unos datos de hietograma 𝑈 = (𝑁 − 𝑀 + 1)
de exceso de lluvia. Ahora, se quieren estructuras
las ecuaciones, en donde se tiene la relación de En donde
caudales con pulsos y los hidrogramas unitarios. Se
despeja U de cada ecuación para hacer el proceso 𝑈 = 𝐻𝑖𝑑𝑟𝑜𝑔𝑟𝑎𝑚𝑎𝑠 𝑢𝑛𝑖𝑡𝑎𝑟𝑖𝑜𝑠
de deconvolución. En la ilustración 4, se puede 𝑀 = 𝑁° 𝑑𝑒 𝑝𝑢𝑙𝑠𝑜𝑠 𝑑𝑒 𝑒𝑥𝑐𝑒𝑠𝑜 𝑑𝑒 𝑙𝑙𝑢𝑣𝑖𝑎
identificar las ecuaciones para un caso donde se 𝑁 = 𝑁° 𝑑𝑒 𝑝𝑢𝑙𝑠𝑜𝑠 𝑑𝑒 𝑒𝑠𝑐𝑜𝑟𝑟𝑒𝑛𝑡í𝑎
tuvieran cinco hidrogramas unitarios, cuatro pulsos
de hietograma de exceso de lluvia y cuatro pulsos 12 = (15 − 4 + 1)
de escorrentía directa.
4
Se obtuvieron 12 valores de hidrograma unitario. diferentes datos, pero todos en base a los valores
Posteriormente, se procede a realizar el siguiente propuestos del hidrograma unitario a través del
procedimiento. intervalo del tiempo de la hora 0 hasta la hora 11 y
del hietograma de exceso de lluvia en un rango de
4 horas. Con lo cual, se obtienen los caudales de
escorrentía directa y ya sea que se realice el
proceso de convolución o deconvolución se pueden
desarrollar.
Tabla 1. HU
t (horas) q (m3/s/cm)
0 0
Ilustración 5. Procedimiento para realizar el cálculo del 1 10
hidrograma unitario por el método matricial. Fuente: propia 2 29
3 46
4 54
5 72
Siendo así, que se realizan estos pasos para llegar 6 85
7 101
a la ecuación: 8 84
9 65
10 43
[𝑈] = [𝑍]−1 [𝑃]𝑇 [𝑄] 11 21
Ilustración 6. Valores propuestos de HU.
En donde
Tabla 2. Hietograma - ERH (Exceso de lluvia)
𝑈 = 𝐻𝑖𝑑𝑟𝑜𝑔𝑟𝑎𝑚𝑎 𝑢𝑛𝑖𝑡𝑎𝑟𝑖𝑜. t (horas) P neta (cm)
𝑍 = 𝑃𝑟𝑜𝑑𝑢𝑐𝑡𝑜 𝑑𝑒 𝑙𝑎 𝑡𝑟𝑎𝑛𝑠𝑝𝑢𝑒𝑠𝑡𝑎 𝑑𝑒 𝑃 𝑦 𝑃. 1 5
𝑄 = 𝐶𝑎𝑢𝑑𝑎𝑙𝑒𝑠 𝑑𝑒 𝑒𝑠𝑐𝑜𝑟𝑟𝑒𝑛𝑡í𝑎 𝑑𝑖𝑟𝑒𝑐𝑡𝑎. 2 4
𝑃 = 𝑀𝑎𝑡𝑟𝑖𝑧 𝑐𝑜𝑛 𝑑𝑎𝑡𝑜𝑠 𝑑𝑒𝑙 ℎ𝑖𝑒𝑡𝑜𝑔𝑟𝑎𝑚𝑎. 3 3
4 1
Ilustración 7. Valores propuestos para el hietograma.
Cabe mencionar, que esta fórmula se desarrolló a
través de la deconvolución analizando un
hidrograma multiplico complejo, pero con Ahora, se realiza el primer proceso de convolución
optimizaciones y ajustes de mínimos cuadrados. con los datos del hidrograma unitario y del
hietograma de exceso de lluvia, generando los
RESULTADOS valores de la escorrentía directa o caudal de salida
en unidades de metro cúbico sobre segundo. Esto
Ahora, al realizar el procedimiento de los tres es posible al realizar el producto del valor del
métodos utilizados previamente se obtuvieron hietograma respectivo por el valor del hidrograma
5
unitario que corresponda. Después, se realiza una En esta última gráfica, se puede visualizar el
sumatoria horizontal de los valores acumulados en comportamiento que se predice tendrá la cuenca
el transcurso de las cuatro horas. según los cálculos, evidenciando que el valor mas
alto de caudal calculado será 1151 m3/s que
Caudal de salida
t H.U. 5 4 3 1 Q resultante
ocurrirá aproximadamente entre un rango de la
(horas) q (m3/s/cm) (cm) (cm) (cm) (cm) (m3/s)
0 0,0 0,0 0,00
hora 7 y la 8. Asimismo, se sabe que los máximos
1 10,0 50,0 0,0 50,00
2 29,0 145,0 40,0 0,0 185,00 en cada hora estarán centrados en esa zona y
3 46,0 230,0 116,0 30,0 0,0 376,00
4 54,0 270,0 184,0 87,0 10,0 551,00 simetría negativa, pero en un rango mas abierto de
5 72,0 360,0 216,0 138,0 29,0 743,00
6 85,0 425,0 288,0 162,0 46,0 921,00 la hora 6 hasta la hora 10.
7 101,0 505,0 340,0 216,0 54,0 1115,00
8 84,0 420,0 404,0 255,0 72,0 1151,00 Ahora, para comprobar que el proceso se realizo
9 65,0 325,0 336,0 303,0 85,0 1049,00
10 43,0 215,0 260,0 252,0 101,0 828,00 debidamente se implementara el primer proceso de
11 21,0 105,0 172,0 195,0 84,0 556,00
12 84,0 129,0 65,0 278,00 convolución con la represión lineal múltiple.
13 63,0 43,0 106,00
14 21,0 21,00 Primero, se crea la matriz con los valores del
Ilustración 8. Caudal de salida (escorrentía directa)
hietograma y luego se sigue el procedimiento de
análisis de datos de Excel para la regresión.
A continuación, se visualizan las graficas
Q U1 U2 U3 U4 U5 U6 U7 U8 U9 U10 U11 U12
6
Hietograma de exceso de lluvia e hidrograma de escorrentía directa
METODO GAUSS JORDAN Hietograma
Escorrentía
Tiempo (Exceso de H.U.
directa n
lluvia)
3
(horas) (cm) (m3/s/cm)
Para este método, se tomaron los valores de 1 5
(m /s)
0,00 U1 0,0
escorrentía directa generados en la solución 2 4 50,00 U2 10,0
3 3 185,00 U3 29,0
previamente realizada por el método de la 4 1 376,00 U4 46,0
5 551,00 U5 54,0
regresión lineal, en donde junto al hietograma se 6 743,00 U6 72,0
7 921,00 U7 85,0
realiza el procedimiento para la deconvolución, 8 1115,00 U8 101,0
9 1151,00 U9 84,0
que al realizar el proceso inverso a la convolución 10 1049,00 U10 65,0
11 828,00 U11 43,0
conllevara a los valores del hidrograma unitario 12 556,00 U12 21,0
1 2 3 4 5 6 7 8 9 10 11 12
1 5,00 0,00 0,00 0,00 0,00 0,00 0,00 0,00 0,00 0,00 0,00 0,00
2 4,00 5,00 0,00 0,00 0,00 0,00 0,00 0,00 0,00 0,00 0,00 0,00
3 3,00 4,00 5,00 0,00 0,00 0,00 0,00 0,00 0,00 0,00 0,00 0,00
4 1,00 3,00 4,00 5,00 0,00 0,00 0,00 0,00 0,00 0,00 0,00 0,00
5 0,00 1,00 3,00 4,00 5,00 0,00 0,00 0,00 0,00 0,00 0,00 0,00
6 0,00 0,00 1,00 3,00 4,00 5,00 0,00 0,00 0,00 0,00 0,00 0,00
7 0,00 0,00 0,00 1,00 3,00 4,00 5,00 0,00 0,00 0,00 0,00 0,00
[P] 8 0,00 0,00 0,00 0,00 1,00 3,00 4,00 5,00 0,00 0,00 0,00 0,00
9 0,00 0,00 0,00 0,00 0,00 1,00 3,00 4,00 5,00 0,00 0,00 0,00
10 0,00 0,00 0,00 0,00 0,00 0,00 1,00 3,00 4,00 5,00 0,00 0,00
11 0,00 0,00 0,00 0,00 0,00 0,00 0,00 1,00 3,00 4,00 5,00 0,00
12 0,00 0,00 0,00 0,00 0,00 0,00 0,00 0,00 1,00 3,00 4,00 5,00
13 0,00 0,00 0,00 0,00 0,00 0,00 0,00 0,00 0,00 1,00 3,00 4,00
14 0,00 0,00 0,00 0,00 0,00 0,00 0,00 0,00 0,00 0,00 1,00 3,00
15 0,00 0,00 0,00 0,00 0,00 0,00 0,00 0,00 0,00 0,00 0,00 1,00
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15
7
Paso 4: Multiplicación Matriz [P] T * Matriz [P ]:
[P] T[P]=[Z]
1
2
3
4
1
51,0
35,0
19,0
5,0
2
35,0
51,0
35,0
19,0
3
19,0
35,0
51,0
35,0
4
5,0
19,0
35,0
51,0
5
0,0
5,0
19,0
35,0
6
0,0
0,0
5,0
19,0
7
0,0
0,0
0,0
5,0
8
0,0
0,0
0,0
0,0
9
0
0
0
0
10
0
0
0
0
11
0
0
0
0
12
0
0
0
0
generan los mismos valores en cada calculo, siendo
5 0,0 5,0 19,0 35,0 51,0 35,0 19,0 5,0 0 0 0 0
6 0,0 0,0 5,0 19,0 35,0 51,0 35,0 19,0 5 0 0 0
7
8
9
10
11
0,0
0,0
0,0
0,0
0,0
0,0
0,0
0,0
0,0
0,0
0,0
0,0
0,0
0,0
0,0
5,0
0,0
0,0
0,0
0,0
19,0
5,0
0,0
0,0
0,0
35,0
19,0
5,0
0,0
0,0
51,0
35,0
19,0
5,0
0,0
35,0
51,0
35,0
19,0
5,0
19
35
51
35
19
5
19
35
51
35
0
5
19
35
51
0
0
5
19
35
así que permiten una exactitud similar en cada
12 0,0 0,0 0,0 0,0 0,0 0,0 0,0 0,0 5 19 35 51
[Z] -1
2
3
4
5
6
-0,031999
0,001602
0,009917
-0,002497
-0,004268
0,065596
-0,033276
-0,006331
0,011904
0,000922
-0,033276
0,065653
-0,032882
-0,006410
0,011723
-0,006331
-0,032882
0,068111
-0,033494
-0,007472
0,011904
-0,006410
-0,033494
0,068213
-0,033202
0,000922
0,011723
-0,007472
-0,033202
0,068656
-0,006590
0,000998
0,012436
-0,007548
-0,033566
0,002337
-0,006502
0,001192
0,012237
-0,007548
0,001850125
0,002313567
-0,006852516
0,001191884
0,012436415
-0,001528062
0,001717533
0,002313567
-0,006501944
0,000998475
-0,000256869
-0,001528062
0,001850125
0,002336793
-0,006590239
0,000564176
0,000181985
-0,001459795
0,00070176
0,002931476
un ejercicio académico que permita el
7 0,002931 -0,006590 0,000998 0,012436 -0,007548 -0,033566 0,068656 -0,033202 -0,007471552 0,011722937 0,000922447 -0,004267915
8 0,000702 0,002337 -0,006502 0,001192 0,012237 -0,007548 -0,033202 0,068213 -0,033494291 -0,006410194 0,011903721 -0,002497355
9
10
11
12
-0,001460
0,000182
0,000564
-0,000312
0,001850
-0,001528
-0,000257
0,000564
0,002314
0,001718
-0,001528
0,000182
-0,006853
0,002314
0,001850
-0,001460
0,001192
-0,006502
0,002337
0,000702
0,012436
0,000998
-0,006590
0,002931
-0,007472
0,011723
0,000922
-0,004268
-0,033494
-0,006410
0,011904
-0,002497
0,068110545
-0,032881751
-0,00633079
0,00991722
-0,032881751
0,065652548
-0,033276268
0,001601563
-0,00633079
-0,033276268
0,06559638
-0,031999417
0,00991722
0,001601563
-0,031999417
0,03999929
entendimiento de la relación entre los hidrogramas,
Paso 6: Multiplicación Matriz [Z] -1 * Matriz [P ] T:
1
2
3
1
0,2000
-0,1600
0,0080
2
0,0000
0,2000
-0,1600
3
0,0000
0,0000
0,2000
4
0,0000
0,0000
0,0000
5
0,0000
0,0000
0,0001
6
0,0000
0,0000
0,0000
7
0,0000
0,0001
-0,0002
8
-0,0001
0,0000
0,0002
9
0,0000
-0,0003
0,0003
10
0,0001
0,0002
-0,0007
11
-0,0001
0,0005
-0,0003
12
-0,000216441
-0,000940658
0,00226384
13
0,00062707
-4,19657E-05
-0,002138715
14
-0,000371406
0,001435659
-0,000982108
15
-0,000311861
0,000564176
0,000181985
los hietogramas y el calculo del hidrograma
4 0,0496 0,0080 -0,1600 0,2000 0,0000 0,0001 -0,0001 -0,0002 0,0003 0,0002 -0,0009 0,000189708 0,00202476 -0,002529261 -0,001459795
[Z] -1[P] T=[L]
5 -0,0125 0,0495 0,0081 -0,1599 0,1998 -0,0001 0,0006 -0,0004 -0,0009 0,0014 0,0015 -0,005457976 0,003315475 0,004442073 0,00070176
6
7
8
9
10
11
-0,0213
0,0147
0,0035
-0,0073
0,0009
0,0028
-0,0125
-0,0212
0,0145
0,0034
-0,0069
0,0010
0,0495
-0,0126
-0,0211
0,0146
0,0030
-0,0070
0,0080
0,0493
-0,0123
-0,0209
0,0140
0,0029
-0,1598
0,0084
0,0488
-0,0127
-0,0196
0,0142
0,1998
-0,1597
0,0083
0,0487
-0,0124
-0,0196
-0,0003
0,1988
-0,1583
0,0091
0,0454
-0,0129
0,0008
0,0004
0,1979
-0,1589
0,0113
0,0458
-0,0001
0,0022
-0,0018
0,1966
-0,1539
0,0121
-0,0015
-0,0022
0,0054
0,0001
0,1892
-0,1551
0,0008
-0,0041
0,0016
0,0077
-0,0088
0,1878
0,003728263
0,01004747
-0,017596763
-0,006271771
0,038978637
-0,003771158
-0,007046339
-0,002581384
0,01931155
-0,012205243
-0,027770003
0,035515204
0,002204189
-0,011881299
0,004411656
0,023420869
-0,028471579
-0,03040187
0,002931476
-0,004267915
-0,002497355
0,00991722
0,001601563
-0,031999417
unitario.
12 -0,0016 0,0016 0,0022 -0,0052 -0,0012 0,0133 -0,0090 -0,0201 0,0297 0,0359 -0,1263 0,086720692 0,065600473 0,087998453 0,03999929
Ilustración 16. Procedimiento para el método matricial de Sin embargo, es importante tener en cuenta que
convolución. estos métodos para ser utilizados en diseños
hidráulicos reales tienen que tener estudios veraces
Por consiguiente, se evidencio la realización de y con mucha precisión, ya que los valores
todos los cálculos mencionados en el generados podrían prevenir desastres o aprovechar
procedimiento para al final obtener los valores del situaciones que se puedan presentar.
hidrograma unitario acordes con los previamente Además, es importante mencionar que la
calculados y propuestos. hidrología y toda la rama de la mecánica de fluidos
abarcan temas muy importantes y complejos que
Paso 7: Multiplicación Matriz [L] * [Q]:
HU
implican estudios detallados y experimentos
1 0,00 q (m3/s/cm)
2
3
50,00
185,00 1
1
0
específicos para situaciones modeladas que
4 376,00 2 10
5
6
551,00
743,00
3
4
29
46
permitan comprobar la veracidad de los
7 921,00 5 54
[Q] 8
9
1115,00
1151,00
[U]=[L]*[Q]
6
7
72
85 procedimientos realizados en campo.
10 1049,00 8 101
11 828,00 9 84
12 556,00 10 65
13 278,00 11 43
14
15
106,00
21,00
12 21 REFERENCIAS
Paso 8: Deconvolución:
Ilustración 17. Valores calculados por el método matricial.
Caudal de salida
t
(horas)
H.U.
(m3/s)
5
(cm)
4
(cm)
3
(cm)
1
(cm)
Q resultante
(m3/s)
[1] A. P. José Luis Ayuso, «ESTIMACION DEL
0
1
0,0
10,0
0,0
50,0 0,0
0
50
HIDROGRAMA UNITARIO. ESTUDIO
2 29,0 145,0 40,0 0,0 185
Y una vez mas para comprobar que los datos están
3 46,0 230,0 116,0 30,0 0,0 376 COMPARATIVO DE CUATRO METROS
4 54,0 270,0 184,0 87,0 10,0 551
5 72,0 360,0 216,0 138,0 29,0 743 LINEALES,» E.T.S Ingenieros Agronomos y de
bien calculados, se realiza una vez mas el proceso
6
7
85,0
101,0
425,0
505,0
288,0
340,0
162,0
216,0
46,0
54,0
921
1115
Montes, vol. 1, nº 10, p. 12, 2018.
8 84,0 420,0 404,0 255,0 72,0 1151
CONCLUSIONES
Para finalizar, se puede concluir que los métodos
de regresión lineal, Gauss Jordan y el matricial