Tema: Problemas con parámetros Weibull

Hola compañeros. Mi problema es que tengo un monton de datos de velocidades de viento en excel, pero no se como determinar los parámetros de la distribución de Weibull. Se que se puede hacer mediante una regresión lineal, que te da una recta a partir de la cual obtener el parámetro de forma y de escala, así como la velocidad media. ¿Alguien sabe como hacerlo o donde puedo encontrar dicha información?

Gracias por anticipado.

Si tienes todos los datos, con fórmulas de excel puedes obtener el parametro forma y escala.

Si cuelgas la hoja con unos pocos datos lo vemos.

Un saludo.

A lo mejor este hilo te sirve:

http://www.solarweb.net/forosolar/en...r-weibull.html

Adjunto la hoja de excel (en realidad es openoffice) con unos cuantos datos, gracias por vuestras respuestas.

En cuanto a lo que me dices Sei, ya lo mire, pero los enlaces a la hojas de calculo están rotos.

La hoja que adjuntaste venía bastante mal, dos columnas que supongo serán de velocidades, no se observan direcciones.

Te puse las fórmulillas ahora la tienes que completar para que coja todos los datos, si tienes alguna duda comentamos.

Un saludo.

Hola wenner, el la hoja de cálculo aparecian de izquierda a derecha: fecha, velocidad tomada cada 20 min y velocidad maxima de racha en 20 min, deberian aparecer tres columnas.

Lo he echado un vistazo y deduzco que lo que se hace es linearizar la distribución de weibull mediante Ln y obtener la recta. La pendiente es el factor de forma, pero el factor de escala no se muy bien como se obtiene.
¿podrias ponerme la ecuación de la recta mas detalladamente, para saber donde se encuentran los parámetros (ya sabes, pendiente, ordenada en el origen, corte en abcisas, etc)

Simplemente tienes que hacer una gráfica con:
Eje X: ln(v)
Eje Y: ln(-ln(1-p(v))

ya, pero lo que realmente quiero saber es como se obtiene esa ecuación a partir de la distribución de weibull, vamos ... que de donde sale?

Gracias por tu rapidez en contestar.

Tienes que hacer el ajuste por mínimos cuadrados. Con excel sacas la ecuación con la línea de tendencia.

Si, y obtengo una recta con una pendiente y una ordenada en el origen, del tipo:

y = 4 x + 67

según lo que tu me dices y=ln(-ln(1-p(v)); x=ln v. Lo que no entiendo es de donde sale que x e y sean esos neperianos.

Perdón por meterme, me interesa este tema. "p" y "v" qué son?

Iniciado por Breitling

Perdón por meterme, me interesa este tema. "p" y "v" qué son?

A lo que llamo "p" es la probabilidad o frecuencia de esa velocidad "v".

Iniciado por wenner

A lo que llamo "p" es la probabilidad o frecuencia de esa velocidad "v".

Algo se me escapa, me echas una mano? Qué falla en esta gráfica con unas velocidades y frecuencias? Es Excel normal.

Si excel es normal, no te lo he explicado correctamente es la probabilidad o frecuencia acumulada, yo lo hago con la función frecuencia de excel y lo divido entre el número de lecturas del anemómetro, prueba ahora y ya nos contarás.

Iniciado por wenner

Si excel es normal, no te lo he explicado correctamente es la probabilidad o frecuencia acumulada, yo lo hago con la función frecuencia de excel y lo divido entre el número de lecturas del anemómetro, prueba ahora y ya nos contarás.

Si no entiendo mal eso es lo que ya tengo: número de horas totales en las que el viento ha soplado a una determinada velocidad. O lo que es lo mismo: número total de minutos (muestreo el anemo cada minuto) en quela velocidad medida ha sido "v". Eso es frecuencia acumulada, no?

Perdón, perdón.... Sumo las frecuencias de todas las velocidades, y a cada una de ellas la divido por esa suma. Ahora sí, no?

Pero frecuencia relativa, por eso lo divido entre el número de mediciones.

Cuando llegues a la frecuencia para la velocidad máxima, el probabilidad acumulada será igual a 1 ya que recoge la suma de todas las frecuencias relativas acumuladas hasta esa velocidad, como es la máxima recogerá todas las que se han medido.

P.D. te has adelantado

Creo que ahora sí...

Yo no cojo la velocidad 0 y la probabilidad solo hasta que sea igual a 1 (quiere decir que están todos los registros de viento).

Iniciado por wenner

Yo no cojo la velocidad 0 y la probabilidad solo hasta que sea igual a 1 (quiere decir que están todos los registros de viento).

Lógico, lo otro era falsear datos. Muchas gracias !!

Hola de nuevo, perdona si me pongo pesado, pero necesito saber porque haces la gráfica de Ln v Vs Ln-(Ln(1-p(v)), el resto de cálculos en la hoja los entiendo perfectamente.

Partiendo de la función de distribución de Weibull:
F(v)=1-e^-(v/c)^k

Para determinar c y k, hacemos una aproximación por mínimos cuadrados, tomamos logaritmos por ambos lados:
-Ln(1-F(v)) = (v/c)^k

Volvemos a emplear logaritmo por ambos lados.
Ln(-Ln(1-F(v)) = kLnv - kLnc

y = ax + b
y=primer término
x=Lnv
b=-kLnc

Donde a = k y c=e^-b/k.

Ok wenner, ahora sí.
Casí me vuelvo loco intentando linearizar otra expresión diferente.
Muchas gracias por todo.

wenner

Muy buena la explicación

Llegué al foro buscando un método para determinar los parámetros de la curva weibull y ahora tengo una forma simple y buena.

Pero a pesar de todo tengo un problema

¿Que pasa si los datos no se ajustan de buena forma a la curva de weibull?
¿Hay otra distribución probabilista que se aplique en estos casos? binomial, t de student, etc

La verdad que la curva de weibull no me queda muy bien con los datos.

Despues se me ocurrio hacer la curva pero considerando la direccion del viento por octante, tipo norte, noreste, etc. a ver si asi mejoraba.

Esto se hace, o se toma la distribucion de viento independiente la direccion de este.

Me parece que si se trata de un molino que enfrenta siempre la misma direccion de viento, la distribucion debería hacerse solo con los vientos que vienen desde esa direccion.

Ojalá haya quedado clara mi duda.

Saludos y gracias

Como no contesta nadie, te doy mi opinión.

Iniciado por finflaman

¿Que pasa si los datos no se ajustan de buena forma a la curva de weibull?
¿Hay otra distribución probabilista que se aplique en estos casos? binomial, t de student, etc
No conozco otra aplicable.

La verdad que la curva de weibull no me queda muy bien con los datos.
Utilizas lo registros de una año completo?. Cuantos más datos tengas el ajuste puede ser mejor.

Puedes probar con la distribución de Rayleigh (no sé si lo he escrito bien), es una caso particular de Weibull, tomando como factor de forma k = 2.

Despues se me ocurrio hacer la curva pero considerando la direccion del viento por octante, tipo norte, noreste, etc. a ver si asi mejoraba.
Y como te ha quedado?.

Esto se hace, o se toma la distribucion de viento independiente la direccion de este.

Me parece que si se trata de un molino que enfrenta siempre la misma direccion de viento, la distribucion debería hacerse solo con los vientos que vienen desde esa direccion.
No lo tengo muy claro.

Si cuelgas lo que tengas le echamos un vistazo.

Un saludo.