Para un registro de datos hidrometeorologicos medios se emplea la dist.

Normal
Dada la serie historica de caudales medios anuales, en m3/s, correspondiente a un registro de
40 años, los mismos que muestran en la siguiente tabla.

Datos: 200 480 430 880 470 350 630 765

555 370 695 690 730 340 800 420
540 400 250 420 750 780 360 530
580 520 550 610 590 536 690 548
530 570 450 650 320 680 610 290

Realizar la prueba de bondad de ajuste de Smirnov Kolmogorov, para verificar si los datos de la
serie se ajustan a la distribucion Normal para α=0.05. (α: nivel de significancia)
En caso que se ajuste la distribucion teorica mencionada, calcular.
a. La probabilidad que el caudal anual, sea menor o igual que 850m 3/s.
b. La probabilidad que el caudal anual, sea mayor o igual que 520m 3/s
c. La probabilidad que el caudal este comprendido entre 350 m 3/s y 800 m3/s.
d. El caudal anual que corresponde a un periodo de retorno de 50 años.
e. El periodo de retorno correspondiente a un caudal medio anual de 780m 3/s.

1 Selección de distribución teórica

Se elige la distribución Normal o Gaussiana, porque la seria de datos son valores medios

2 Estimación de parámetros
Los parámetros en la distribución Normal son el parametro de posición que es igual a la
media aritmética, y el parámetro de escala que numericamente es igual a la desviación

Descripción Fórmula Resultado Unidades

μ: Parametro de posición
σ𝑛𝑖=1 𝑥𝑖
𝑥ҧ = 538.975 m3/s
(media aritmetica). 𝑛

α: Parametro de escala, σ𝑛𝑖=1 𝑥𝑖 − 𝑥ҧ 2

(desviación estandar). 𝑆= 161.8256 m3/s
𝑛−1
3 Prueba de bondad de ajuste
Se emplea la prueba de bondad de ajuste de Smirnov Kolmogorov, (Probabilidad empírica
Weibull).

3.1 Calculo de la diferencia abs(P(x)-F(x))

Orden Caudal 𝑚 𝑥 − 𝑥ҧ ∆
𝑃 𝑥 = 𝑧= 𝐹
m m3/s 𝑛+1 𝑆 𝑧 𝑃 𝑥 −𝐹 𝑧

1 200 0.024 -2.095 0.018 0.006

2 250 0.049 -1.786 0.037 0.012
3 290 0.073 -1.539 0.062 0.011
4 320 0.098 -1.353 0.088 0.010
5 340 0.122 -1.230 0.109 0.013
6 350 0.146 -1.168 0.121 0.025
7 360 0.171 -1.106 0.134 0.036
8 370 0.195 -1.044 0.148 0.047
9 400 0.220 -0.859 0.195 0.024
10 420 0.244 -0.735 0.231 0.013
11 420 0.268 -0.735 0.231 0.037
12 430 0.293 -0.673 0.250 0.042
13 450 0.317 -0.550 0.291 0.026
14 470 0.341 -0.426 0.335 0.006
15 480 0.366 -0.364 0.358 0.008
16 520 0.390 -0.117 0.453 0.063
17 530 0.415 -0.055 0.478 0.063
18 530 0.439 -0.055 0.478 0.039
19 536 0.463 -0.018 0.493 0.029
20 540 0.488 0.006 0.503 0.015
21 548 0.512 0.056 0.522 0.010
22 550 0.537 0.068 0.527 0.009
23 555 0.561 0.099 0.539 0.022
24 570 0.585 0.192 0.576 0.009
25 580 0.610 0.254 0.600 0.010
26 590 0.634 0.315 0.624 0.010
27 610 0.659 0.439 0.670 0.011
28 610 0.683 0.439 0.670 0.013
29 630 0.707 0.562 0.713 0.006
30 650 0.732 0.686 0.754 0.022
31 680 0.756 0.871 0.808 0.052
32 690 0.780 0.933 0.825 0.044
33 690 0.805 0.933 0.825 0.020
34 695 0.829 0.964 0.833 0.003
35 730 0.854 1.180 0.881 0.027
36 750 0.878 1.304 0.904 0.026
37 765 0.902 1.397 0.919 0.016
38 780 0.927 1.489 0.932 0.005
39 800 0.951 1.613 0.947 0.005
40 880 0.976 2.107 0.982 0.007
 3.2 Determinar la máxima diferencia ∆

∆= 𝑚á𝑥 𝑃 𝑥 −𝐹 𝑥

∆= 0.063

3.3 Calcular el valor crítico del estadístico ∆o para un nivel de significancia(𝜶), número de datos(n).
Los valores de ∆o, se obtines a partir de la tabla a partir de nivel de significancia y número de
datos.
Datos: α: 0.05 Nivel de significancia
N: 40 Número de datos

Tamaño Nivel de significancia α

muestral N 0.2 0.1 0.05 0.01
5 0.45 0.51 0.56 0.67
10 0.32 0.37 0.41 0.49
15 0.27 0.3 0.34 0.4
20 0.23 0.26 0.29 0.36
25 0.21 0.24 0.27 0.32
30 0.19 0.22 0.24 0.29
35 0.18 0.2 0.23 0.27
40 0.17 0.19 0.21 0.25
45 0.16 0.18 0.2 0.24
50 0.15 0.17 0.19 0.23
1.07 1.22 1.36 1.63
𝑁 > 50
𝑁 𝑁 𝑁 𝑁

∆o= 0.21 Valor obtenido al intepolar los datos

3.4 Comparar el valor del estadístico ∆, con el valor crítico ∆o.

Comparando los valores de ∆, con el valor crítico ∆o (tabla) en base a los siguientes críterios
de decisión.

∆< ∆𝑜 → 𝐸𝑙 𝑎𝑗𝑢𝑠𝑡𝑒 𝑒𝑠 𝑏𝑢𝑒𝑛𝑜, 𝑎𝑙𝑎 𝑛𝑖𝑣𝑒𝑙 𝑑𝑒 𝑠𝑖𝑔𝑛𝑖𝑓𝑖𝑐𝑎𝑛𝑐𝑖𝑎 𝑠𝑒ñ𝑎𝑙𝑎𝑑𝑜

∆≥ ∆𝑜 → 𝐸𝑙 𝑎𝑗𝑢𝑠𝑡𝑒 𝑛𝑜 𝑒𝑠 𝑏𝑢𝑒𝑛𝑜, 𝑎𝑙𝑎 𝑛𝑖𝑣𝑒𝑙 𝑑𝑒 𝑠𝑖𝑔𝑛𝑖𝑓𝑖𝑐𝑎𝑛𝑐𝑖𝑎 𝑠𝑒ñ𝑎𝑙𝑎𝑑𝑜

0.06325 ˂ 0.21 Ajuste es bueno

Se concluye que los datos de caudales se ajustan a la distribución normal, con un nivel de
significancia del 5% o una probabilidad del 95%
Como el ajuste el bueno se puede utlilizar la distribución NORMAL.
 4 Resolución de los enunciados propuestos
Variable reducida, función densidad, función de distribución acumulada, función densidad
estandarizada y función de distribución.

𝑥 − 𝑥ҧ Variable 1
𝑧= reducida
𝑆
1 1 𝑥−𝑥ҧ 2 1 1 2 Función densidad
−2 𝑆
𝑓𝑥 = 𝑒 Función densidad 𝑓𝑧 = 𝑒 −2𝑧 estandarizada
𝑥 2𝜋𝑆 𝑆 2𝜋
𝑥 2 𝑧
1 1 𝑥−𝑥ҧ
−2 𝑆 Función de distribución 1 1 2 Función de distribución
𝐹 𝑥 = න 𝑒 𝑑𝑥 acum. 𝐹 𝑍 = න 𝑒 −2𝑧 𝑑𝑧 acum. Est.
2𝜋𝑆 −∞ 2𝜋 −∞

a. La probabilidad de que un caudal anual, sea menor o igual que 850m 3/s.

𝑃 𝑄 ≤ 𝑞 = 𝑃 𝑄 ≤ 850 = 𝐹 𝑧

Calculo de Z para q= 850 m3/s

𝑞 − 𝑄ത
𝑧=
𝑆𝑞
z= 1.92198
𝑃 𝑄 ≤ 850 = 0.9727

b. La probabilidad que el caudal anual, sea mayor o igual que 520m 3/s

𝑃 𝑄 ≥ 𝑞 = 𝑃 𝑄 ≥ 520 = 1 − 𝐹 𝑧

Calculo de Z para q= 520 m3/s

𝑞 − 𝑄ത
𝑧= z= -0.1173
𝑆𝑞
𝑃 𝑄 ≤ 520 = 0.4533

𝑃 𝑄 ≥ 520 = 1 − 𝐹 𝑧 0.5467

c. La probabilidad que el caudal este comprendido entre 350 m 3/s y 800 m3/s.

𝑃 𝑞1 ≤ 𝑄 ≤ 𝑞2 = 𝑃 350 ≤ 𝑄 ≤ 800 = 𝐹 𝑍2 −𝐹 𝑍1 q1= 350 m3/s

q2= 800 m3/s
𝑞 − 𝑄ത
𝑧=
𝑆𝑞
z1= -1.1678
z2= 1.6130

𝐹 𝑍2 𝑃 𝑄 ≤ 800 = 0.946628 𝐹 𝑍1 𝑃 𝑄 ≤ 350 = 0.1214

𝑃 350 ≤ 𝑄 ≤ 800 = 𝐹 𝑍2 −𝐹 𝑍1 = 0.8252

 d. El caudal anual que corresponde a un periodo de retorno de 50 años.

1 1
𝑃 𝑄≥𝑞 = 𝑃 𝑄 ≤𝑞 =1−
𝑇 𝑇
T= 50 años

𝑃 𝑄≤𝑞 = 0.98

mediante tablas calculo de z.

z= 2.0537

𝑞 − 𝑄ത q= 3
871.3241 m /s
𝑧= → 𝑞 = 𝑧𝑆𝑞 + 𝑄
𝑆𝑞

e. El periodo de retorno correspondiente a un caudal medio anual de 780m 3/s.

1 1
𝑃 𝑄≥𝑞 = 𝑇=
𝑇 1−𝑃 𝑄 ≤ 𝑞

q= 780 m3/s

Cálculo de z:

𝑞 − 𝑄ത
𝑧=
𝑆𝑞

z= 1.48941

1
𝑃 𝑄 ≤ 180 = 0.931811 remplazando en 𝑇= T= 14.665 años
1−𝑃 𝑄 ≤ 𝑞
T= 15 años
o ∆o para un nivel de significancia(𝜶), número de datos(n).