コード例 #1
0
ファイル: funcoesBasicas.py プロジェクト: clebsonpy/HidroComp
 def __init__(self, nome_db):
     self.nome_db = nome_db
     self.ler = l.acessarDados(self.nome_db)
     self.ins = inserir.inserirDados(self.nome_db)
     self.estat = est.estatistica()
コード例 #2
0
ファイル: plotTR.py プロジェクト: clebsonpy/HidroComp
import Hidro2016.Estatistica as est
import Ana.Series as se
import Ana.LimiteParcial as lp
import scipy.stats as sc
import numpy as n

#1,4
dadoNa = da.DadosVazao('F:\\Clebson\\Programas\\HidroComp\\Ana\\BancoHidro', 6, 1995, 2012).Dados()
dadAnNa = se.Series(dadoNa).serieMaxAnual()
limiar = lp.LimiteParcial(dadoNa).AchaLimite(3)
dadPaNa = se.Series(dadoNa).serieMaxParcial(limiar)
dataNa, vazaoNa = se.Series(dadoNa).separaDados()
dataAnNa, anualNa = se.Series(dadAnNa).separaDados()
dataPaNa, parcNa = se.Series(dadPaNa).separaDados()
print(limiar)
estaNaAn = est.estatistica()
print(anualNa)
paraAnual = estaNaAn.estimaParametros(anualNa, 'MML', 'GEV')

'''
magAnual = estaNaAn.EstimaMagnitudes(paraAnual)
estaNaPar = est.estatistica(parcNa, 'Parcial')
paraParc = estaNaPar.CalculaParametros()
magParc = estaNaPar.EstimaMagnitudes(paraParc)
proNaAn = estaNaAn.EstimaProbabilidade(magAnual, paraAnual)
proNaPar = estaNaPar.EstimaProbabilidade(magParc, paraParc)
'''
print(paraAnual)

'''
plt.plot(magAnual, proNaAn, 'bo', label = 'Anual')
コード例 #3
0
ファイル: funcoesBasicas.py プロジェクト: clebsonpy/HidroComp

#========================================
c = controle('Banco_Hidro2')
prob = [0.001, 0.05, 0.1, 0.3, 0.5, 0.7, 0.9, 0.95, 0.999]
metodo = ['MML','MOM','MVS']
tamanhoAmostra = 100

aux = []
for m in metodo:
    aux.append(c.prepGra(metodo=m, tamanhoAmostra=tamanhoAmostra, probabilidade=prob))

dadosExt = []
for i in prob:
    dadosExt.append(gev.ppf(i, -0.168462, 6286.926278, 1819.961392))
e = est.estatistica()
accu = []
for me in range(3):
    ac = []
    for d in range(9):
        ac.append(dadosExt[d] - e.calculoAccu(dadosAmostra=aux[me][d], dadoSintetico=dadosExt[d]))
    accu.append(ac)
print(accu)

ax1 = plt.subplot(221)
ax2 = plt.subplot(223)
ax3 = plt.subplot(122)

c.plotGraf2(axes=ax1, accu=accu[1], dadosExt=dadosExt, yd=prob, quan=aux[1], metodo='MOM', tamanhoAmostra=tamanhoAmostra, prob=prob)
c.plotGraf2(axes=ax2, accu=accu[2], dadosExt=dadosExt, yd=prob, quan=aux[2], metodo='MVS', tamanhoAmostra=tamanhoAmostra, prob=prob)
c.plotGraf2(axes=ax3, accu=accu[0], dadosExt=dadosExt, yd=prob, quan=aux[0], metodo='MML', tamanhoAmostra=tamanhoAmostra, prob=prob)