J'ai un histogramme des données de vitesse du vent qui est souvent représenté en utilisant une distribution de Weibull. Je voudrais calculer la forme et les facteurs d'échelle de Weibull qui donnent le meilleur ajustement à l'histogramme.
J'ai besoin d'une solution numérique (par opposition aux solutions graphiques ) car le but est de déterminer la forme weibull par programmation.
Edit: Les échantillons sont collectés toutes les 10 minutes, la vitesse du vent est moyenne sur les 10 minutes. Les échantillons incluent également la vitesse maximale et minimale du vent enregistrée au cours de chaque intervalle, qui sont actuellement ignorées, mais je voudrais les intégrer plus tard. La largeur du bac est de 0,5 m / s
distributions
histogram
java
klonq
la source
la source
Réponses:
L'estimation du maximum de vraisemblance des paramètres de Weibull peut être une bonne idée dans votre cas. Une forme de distribution de Weibull ressemble à ceci:
Où sont des paramètres. Étant donné les observations , la fonction log-vraisemblance estθ,γ>0 X1,…,Xn
Une solution "basée sur la programmation" serait d'optimiser cette fonction en utilisant une optimisation contrainte. Résoudre pour une solution optimale:
En éliminant nous obtenons:θ
Maintenant, cela peut être résolu pour l'estimation ML . Cela peut être accompli à l'aide de procédures itératives standard qui résolvent sont utilisées pour trouver la solution de l'équation comme - Newton-Raphson ou d'autres procédures numériques.γ^
Maintenant peut être trouvé en termes de comme:θ γ^
la source
Utilisez fitdistrplus:
Besoin d'aide pour identifier une distribution par son histogramme
Voici un exemple de la façon dont la distribution Weibull est adaptée:
la source