uy = [1] * 41
uxtempo = [1] * 41
uxtemperatura = [1] * 41
#Execução do MT_PEU
Estime = EstimacaoNaoLinear(Modelo,
simbolos_x=['t', 'T'],
simbolos_y=['y'],
simbolos_param=['ko', 'E'],
Folder='Exemplo1')
Estime.setDados(0, (tempo, uxtempo), (temperatura, uxtemperatura))
Estime.setDados(1, (y, uy))
Estime.setConjunto(tipo='estimacao')
Estime.optimize(initial_estimative=[0.5, 25000], algoritmo='ipopt')
Estime.incertezaParametros(metodoIncerteza='Geral' '')
Estime.predicao()
# Valores originais
hessian = array([[5.46625676e+00, -7.50116238e-03],
[-7.50116238e-03, 1.02960714e-05]])
sensibilidade = array([[-1.06335451e-01, 1.52826132e-04],
[-5.59353653e-02, 8.03907392e-05],
[-1.28132438e-01, 1.80542088e-04],
[-2.26100338e-01, 3.18581523e-04],
[-2.26100338e-01, 3.18581523e-04],
[-1.28132438e-01, 1.80542088e-04],
[-2.07847997e-01, 2.89084593e-04],
[-1.16257751e-01, 1.61696650e-04],
[-6.14815711e-02, 8.55114088e-05],
[-2.07847997e-01, 2.89084593e-04],
# Otimização
Estime.otimiza(limite_superior=sup,
limite_inferior=inf,
algoritmo='PSOFamily',
itmax=500,
Num_particulas=30,
metodo={
'busca': 'Otimo',
'algoritmo': 'PSO',
'inercia': 'TVIW-linear'
},
args=[tipo],
printit=True)
#Estime.SETparametro([3.0784923406298477, 0.10246446218717528, 5.0419496833229021, 0.43170222912546136])
Estime.incertezaParametros(delta=1e-5, metodo='SensibilidadeModelo')
Estime.predicao()
Estime.analiseResiduos()
etapas = [
'otimizacao', 'grandezas-entrada', 'predicao', 'grandezas-calculadas',
'analiseResiduos', 'regiaoAbrangencia'
]
Estime.graficos(etapas)
Estime.relatorio(export_y=True, export_cov_y=True, quebra='\r\n')
##################################################################################
##################################################################################
# EXEMPLOS PARA MODELOS LINEARES
##################################################################################
##################################################################################
# =================================================================================