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 ################################################################################## ################################################################################## # =================================================================================