<jupyter_start><jupyter_text><jupyter_code>from google.colab import drive

drive.mount('/content/drive')

#Instalar biblioteca lmfit
!pip install lmfit

import matplotlib.pyplot as plt
import pandas as pd
import numpy as np
#import lmfit<jupyter_output>Mounted at /content/drive
Collecting lmfit
[?25l  Downloading https://files.pythonhosted.org/packages/d1/2e/9e4cdcebf93603ce57bc7577707ad56ed1282ae306f611e525612b9cfeea/lmfit-1.0.1.tar.gz (258kB)
[K     |████████████████████████████████| 266kB 8.3MB/s 
[?25hCollecting asteval>=0.9.16
[?25l  Downloading https://files.pythonhosted.org/packages/4f/d5/ce490afa6f6ab1bb97ffd66a80bfb032d22407fc215d60b21ca174e25c51/asteval-0.9.21.tar.gz (54kB)
[K     |████████████████████████████████| 61kB 7.5MB/s 
[?25hRequirement already satisfied: numpy>=1.16 in /usr/local/lib/python3.6/dist-packages (from lmfit) (1.19.4)
Requirement already satisfied: scipy>=1.2 in /usr/local/lib/python3.6/dist-packages (from lmfit) (1.4.1)
Collecting uncertainties>=3.0.1
[?25l  Downloading https://files.pythonhosted.org/packages/45/41/fc7e7b73b603e7c2c9e040b7aa8caf4a88d74b6faa567601ed82b6f0d8e1/uncertainties-3.1.5-py2.py3-none-any.whl (246kB)
[K     |████████████████████████████████| 256kB 12.7MB/s 
[?25hRequirem[...]<jupyter_text>Importar data frames<jupyter_code>#garrafa térmica aberta
sheet_id_garrafa_aberta='1sN2kTrNoocruN402abYPbnojdXnaoUhNwtNgtOJrALY'

garrafa_aberta=pd.read_csv(f'https://docs.google.com/spreadsheets/d/{sheet_id_garrafa_aberta}/export?format=csv')

#calorímetro aberto
cal_aberto=pd.read_excel('/content/drive/MyDrive/F 229/Resfriamento/calorímetro aberto.xlsx')
t_cal_aberto=cal_aberto['tempo médio']
T_cal_aberto=cal_aberto['Kelvin']

#calorímetro fechado
cal_fechado=pd.read_excel('/content/drive/MyDrive/F 229/Resfriamento/calfechado.xlsx')
t_cal_fechado=cal_fechado['t(s)']
T_cal_fechado=cal_fechado['T(K)']


#G<NAME>
garrafa_fechada=pd.read_excel('/content/drive/MyDrive/F 229/Resfriamento/Garrafa Térmica Fechada.xlsx')
t_garrafa_fechada=garrafa_fechada['tempo']
T_garrafa_fechada=garrafa_fechada['Temperatura']

#Garrafa de Vidro
garrafa_vidro=pd.read_excel('/content/drive/MyDrive/F 229/Resfriamento/Garrafa de vidro.xlsx')
t_vidro=garrafa_vidro['média']
T_vidro=garrafa_vidro['temperat']

garrafa_vidro<jupyter_output><empty_output><jupyter_text>Gráficos dos Dados Experimentais preliminares<jupyter_code>plt.scatter(garrafa_aberta['Tempo (s)'],garrafa_aberta['Temperatura (K)'],label='Pontos Experimentais')
plt.scatter(garrafa_aberta['Tempo2 (s)'][0:6],garrafa_aberta['Temperatura2 (K)'][0:6],label='Ponto Médio')
plt.errorbar(garrafa_aberta['Tempo2 (s)'][0:6],garrafa_aberta['Temperatura2 (K)'][0:6],xerr=np.array([3.75277675,15.29978213,13.85640646,16.16580754,10.96965511,2.020725942]),yerr=1,fmt='none',label='Incerteza',ecolor='black',elinewidth=2)
plt.xlabel('Tempo (s)',fontsize=12)
plt.ylabel('Temperatura (K)',fontsize=12)
plt.grid()
plt.legend(fontsize='small')



fig,ax1=plt.subplots(1,1)

ax1.scatter(garrafa_aberta['Tempo2 (s)'][0:6],garrafa_aberta['Temperatura2 (K)'][0:6],label='Média dos pontos')
ax1.errorbar(garrafa_aberta['Tempo2 (s)'][0:6],garrafa_aberta['Temperatura2 (K)'][0:6],xerr=np.array([0.5773502692,15.29978213,13.85640646,16.16580754,10.96965511,2.020725942]),yerr=1,fmt='none',label='Incerteza',ecolor='black',elinewidth=2)

ax1.set_ylabel('Temperatura (K)',fontsize=12)
ax1.set_xlabel('Tempo (s)',fontsize=12)

ax1.legend(fontsize='small')
ax1.grid()

plt.show()<jupyter_output><empty_output><jupyter_text>Linearização através da primeira equação<jupyter_code>fig,ax1 = plt.subplots(1,1)

from scipy.stats import linregress

t=garrafa_aberta['Tempo2 (s)'][1:6]
T=garrafa_aberta['Temperatura2 (K)'][1:6]

t3= np.log(t)
T3 = np.log(T)
ax1.plot(t, T3,'o',label='<NAME>')
ax1.errorbar(t,T3,xerr=np.array([15.29978213,13.85640646,16.16580754,10.96965511,2.020725942]),yerr=0.002,fmt='none',label='Incerteza',ecolor='black',elinewidth=2)

a,b,_,_,_,= linregress(t,T3)
t_r = np.linspace(3,250)
ax1.plot(t_r,a*t_r+b, label='Reta de Regressão Linear')

plt.xlabel('(Tempo (s)')
plt.ylabel('$\ln$ (Temperatura (K))')
plt.legend()
ax1.grid()
# plt.savefig(pathfig+'fig_garrafa_termica_fechada.png')
plt.show()

a<jupyter_output><empty_output><jupyter_text>Teste da Equação 1 Linearizada<jupyter_code>T_infinito=26.5+273.15
T_0=335.15
gamma=-a

fig,ax1 = plt.subplots(1,1)



t1=np.linspace(0,300)
T1= (T_0 - T_infinito)*np.exp(-gamma*t1) + T_infinito

ax1.scatter(t,T,label='Garra Térmica Aberta')
ax1.errorbar(t,T,xerr=np.array([15.29978213,13.85640646,16.16580754,10.96965511,2.020725942]),yerr=1,fmt='none',label='Incerteza',ecolor='black',elinewidth=2)
ax1.plot(t1,T1,'r',label='Curva Teórica Modelo 1')

plt.xlabel('Tempo (s)')
plt.ylabel('Temperatura (K)')
plt.legend()
plt.grid()

plt.show()

a<jupyter_output><empty_output><jupyter_text>Definição das variações de temperatura<jupyter_code>T_infinito=26.5+273.15

#garrafa térmica aberta
T_0=335.15
Delta_T3=1/(T-T_infinito)**(1/3)
Delta_T4=1/(T-T_infinito)**(1/4)
Delta_T6=1/(T-T_infinito)**(1/6)

#calorímetro aberto
T_0_cal_aberto=340.15
Delta_T3_cal_aberto=1/(T_cal_aberto - T_infinito)**(1/3)
Delta_T4_cal_aberto=1/(T_cal_aberto - T_infinito)**(1/4)
Delta_T6_cal_aberto=1/(T_cal_aberto - T_infinito)**(1/6)

#calorímetro fechado
T_0_cal_aberto=343
Delta_T3_cal_fechado=1/(T_cal_fechado - T_infinito)**(1/3)
Delta_T4_cal_fechado=1/(T_cal_fechado - T_infinito)**(1/4)
Delta_T6_cal_fechado=1/(T_cal_fechado - T_infinito)**(1/6)

#<NAME>
T_0_cal_aberto=341.15
Delta_T3_gar_fechada=1/(T_garrafa_fechada - T_infinito)**(1/3)
Delta_T4_gar_fechada=1/(T_garrafa_fechada - T_infinito)**(1/4)
Delta_T6_gar_fechada=1/(T_garrafa_fechada - T_infinito)**(1/6)

#<NAME>
T_0_cal_aberto=340.15
Delta_T3_vidro=1/(T_vidro - T_infinito)**(1/3)
Delta_T4_vidro=1/(T_vidro - T_infinito)**(1/4)
Delta_T6_vidro=1/(T_vidro - T_infinito)**(1/6)<jupyter_output><empty_output><jupyter_text>Modelo de Equação 2  n=1/3<jupyter_code>fig,ax1 = plt.subplots(1,1)

#<NAME>
ax1.scatter(Delta_T3,t,label='Garra Térmica Aberta n=1/3')
a_3,b_3,_,_,_,= linregress(Delta_T3,t)
t_r = np.linspace(0.28,0.32)
ax1.plot(t_r,a_3*t_r+b_3)

#Calorímetro Aberto
ax1.scatter(Delta_T3_cal_aberto,t_cal_aberto,label='Calorímetro Aberto n=1/3')
a_3_cal_aberto,b_3_cal_aberto,_,_,_,= linregress(Delta_T3_cal_aberto,t_cal_aberto)
t_r = np.linspace(0.28,0.32)
ax1.plot(t_r,a_3_cal_aberto*t_r+b_3_cal_aberto)

#Calorímetro Fechado
ax1.scatter(Delta_T3_cal_fechado,t_cal_fechado,label='Calorímetro fechado n=1/3')
a_3_cal_fechado,b_3_cal_fechado,_,_,_,= linregress(Delta_T3_cal_fechado,t_cal_fechado)
t_r = np.linspace(0.27,0.32)
ax1.plot(t_r,a_3_cal_fechado*t_r+b_3_cal_fechado)

#Garrafa Térmica Fechada
ax1.scatter(Delta_T3_gar_fechada,t_garrafa_fechada,label='Garrafa Térmica Fechada n=1/3')
a_3_gar_fechada,b_3_gar_fechada,_,_,_,= linregress(Delta_T3_gar_fechada,t_garrafa_fechada)
t_r = np.linspace(0.27,0.35)
ax1.plot(t_r,a_3_gar_fechada*t_r+b_3_gar_fechada)

#Garrafa de Vidro
ax1.scatter(Delta_T3_vidro,t_vidro,label='Garrafa de Vidro n=1/3')
a_3_vidro,b_3_vidro,_,_,_,= linregress(Delta_T3_vidro,t_vidro)
t_r = np.linspace(0.27,0.32)
ax1.plot(t_r,a_3_vidro*t_r+b_3_vidro)

plt.ylabel('Tempo (s)')
plt.xlabel('Temperatura (K)')
plt.legend()
plt.grid()

fig,ax2 = plt.subplots(1,1)

#Garrafa Térmica Aberta
T_3=(a_3**3)/(-b_3+t)**3 +T_infinito
ax2.plot(t,T_3,label="Garrafa Térmica Aberta")
#ax2.scatter(t,T)

#Calorímetro Aberto
#T_3_cal_aberto=(a_3_cal_aberto**3)/(-b_3_cal_aberto+t_cal_aberto)**3 +T_infinito
#ax2.plot(t_cal_aberto,T_3_cal_aberto,label="Calorímetro Aberto")
#ax2.scatter(t_cal_aberto,T_cal_aberto)

#Calorímetro Fechado
T_3_cal_fechado=(a_3_cal_fechado**3)/(-b_3_cal_fechado+t_cal_fechado)**3 +T_infinito
ax2.plot(t_cal_fechado,T_3_cal_fechado,label="Calorímetro Fechado")
#ax2.scatter(t_cal_fechado,T_cal_fechado)

#Garrafa Térmica Fechada
T_3_gar_fechada=(a_3_gar_fechada**3)/(-b_3_gar_fechada+t_garrafa_fechada)**3 +T_infinito
ax2.plot(t_garrafa_fechada,T_3_gar_fechada,label="<NAME>")
#ax2.scatter(t_garrafa_fechada,T_garrafa_fechada)

#G<NAME> Vidro
T_3_vidro=(a_3_vidro**3)/(-b_3_vidro+t_vidro)**3 +T_infinito
ax2.plot(t_vidro,T_3_vidro,label="<NAME>")
#ax2.scatter(t_vidro,T_vidro)

plt.xlabel('Tempo (s)')
plt.ylabel('Temperatura (K)')
plt.legend(fontsize=9, loc='upper rigth')
plt.grid()

plt.show()
Delta_T3_cal_aberto<jupyter_output><empty_output><jupyter_text>n=1/4<jupyter_code>fig,ax1 = plt.subplots(1,1)

#<NAME>
ax1.scatter(Delta_T4,t,label='Garra Térmica Aberta n=1/4')
a_4,b_4,_,_,_,= linregress(Delta_T4,t)
t_r = np.linspace(0.38,0.43)
ax1.plot(t_r,a_4*t_r+b_4,'g')

#Calorímetro Aberto
ax1.scatter(Delta_T4_cal_aberto,t_cal_aberto,label='Calorímetro Aberto n=1/4')
a_4_cal_aberto,b_4_cal_aberto,_,_,_,= linregress(Delta_T4_cal_aberto,t_cal_aberto)
t_r = np.linspace(0.38,0.43)
ax1.plot(t_r,a_4_cal_aberto*t_r+b_4_cal_aberto)

#Calorímetro Fechado
ax1.scatter(Delta_T4_cal_fechado,t_cal_fechado,label='Calorímetro fechado n=1/4')
a_4_cal_fechado,b_4_cal_fechado,_,_,_,= linregress(Delta_T4_cal_fechado,t_cal_fechado)
t_r = np.linspace(0.38,0.42)
ax1.plot(t_r,a_4_cal_fechado*t_r+b_4_cal_fechado)

#<NAME>
ax1.scatter(Delta_T4_gar_fechada,t_garrafa_fechada,label='<NAME> n=1/3')
a_4_gar_fechada,b_4_gar_fechada,_,_,_,= linregress(Delta_T4_gar_fechada,t_garrafa_fechada)
t_r = np.linspace(0.27,0.35)
ax1.plot(t_r,a_4_gar_fechada*t_r+b_4_gar_fechada)

#<NAME>
ax1.scatter(Delta_T4_vidro,t_vidro,label='<NAME> n=1/4')
a_4_vidro,b_4_vidro,_,_,_,= linregress(Delta_T4_vidro,t_vidro)
t_r = np.linspace(0.27,0.32)
ax1.plot(t_r,a_4_vidro*t_r+b_4_vidro)


plt.ylabel('Tempo (s)')
plt.xlabel('Temperatura (K)')
plt.legend()
plt.grid()

print(a_4)
print(b_4)

fig,ax2 = plt.subplots(1,1)

#<NAME>
T_4=(a_4**4)/(-b_4+t)**4 +T_infinito
ax2.plot(t,T_4,label="Garrafa Térmica Aberta")
#ax2.scatter(t,T)

#Calorímetro Aberto
#T_4_cal_aberto=(a_4_cal_aberto**4)/(-b_4_cal_aberto+t_cal_aberto)**4 +T_infinito
#ax2.plot(t_cal_aberto,T_4_cal_aberto,label="Calorímetro Aberto")
#ax2.scatter(t_cal_aberto,T_cal_aberto)

#Calorímetro Fechado
T_4_cal_fechado=(a_4_cal_fechado**4)/(-b_4_cal_fechado+t_cal_fechado)**4 +T_infinito
ax2.plot(t_cal_fechado,T_4_cal_fechado,label="Calorímetro Fechado")
#ax2.scatter(t_cal_fechado,T_cal_fechado)

#Garrafa Térmica Fechada
T_4_gar_fechada=(a_4_gar_fechada**4)/(-b_4_gar_fechada+t_garrafa_fechada)**4 +T_infinito
ax2.plot(t_garrafa_fechada,T_4_gar_fechada,label="Garrafa Térmica Fechada")
#ax2.scatter(t_garrafa_fechada,T_garrafa_fechada)

#Garrafa de Vidro
T_4_vidro=(a_4_vidro**4)/(-b_4_vidro+t_vidro)**4 +T_infinito
ax2.plot(t_vidro,T_4_vidro,label="Garrafa de Vidro")
#ax2.scatter(t_vidro,T_vidro)


plt.xlabel('Tempo (s)')
plt.ylabel('Temperatura (K)')
plt.legend(fontsize=9, loc='upper rigth')
plt.grid()

plt.show()<jupyter_output>16423.328925798905
-6743.385340579813<jupyter_text>n=1/6<jupyter_code>fig,ax1 = plt.subplots(1,1)

#<NAME>
ax1.scatter(Delta_T6,t,label='G<NAME> Aberta n=1/6')
a_6,b_6,_,_,_,= linregress(Delta_T6,t)
t_r = np.linspace(0.53,0.57)
ax1.plot(t_r,a_6*t_r+b_6,'y')

#Calorímetro Aberto
ax1.scatter(Delta_T6_cal_aberto,t_cal_aberto,label='Calorímetro Aberto n=1/4')
a_6_cal_aberto,b_6_cal_aberto,_,_,_,= linregress(Delta_T6_cal_aberto,t_cal_aberto)
t_r = np.linspace(0.53,0.56)
ax1.plot(t_r,a_6_cal_aberto*t_r+b_6_cal_aberto)

#Calorímetro Fechado
ax1.scatter(Delta_T6_cal_fechado,t_cal_fechado,label='Calorímetro Fechado n=1/6')
a_6_cal_fechado,b_6_cal_fechado,_,_,_,= linregress(Delta_T6_cal_fechado,t_cal_fechado)
t_r = np.linspace(0.53,0.56)
ax1.plot(t_r,a_6_cal_fechado*t_r+b_6_cal_fechado)

#<NAME>
ax1.scatter(Delta_T6_gar_fechada,t_garrafa_fechada,label='Garrafa Térmica Fechada n=1/6')
a_6_gar_fechada,b_6_gar_fechada,_,_,_,= linregress(Delta_T6_gar_fechada,t_garrafa_fechada)
t_r = np.linspace(0.27,0.35)
ax1.plot(t_r,a_6_gar_fechada*t_r+b_6_gar_fechada)

#<NAME>
ax1.scatter(Delta_T4_vidro,t_vidro,label='<NAME> n=1/6')
a_6_vidro,b_6_vidro,_,_,_,= linregress(Delta_T6_vidro,t_vidro)
t_r = np.linspace(0.27,0.32)
ax1.plot(t_r,a_6_vidro*t_r+b_6_vidro)

plt.ylabel('Tempo (s)')
plt.xlabel('Temperatura (K)')
plt.legend()
plt.grid()

fig,ax2 = plt.subplots(1,1)

#<NAME>
T_6=(a_6**6)/(-b_6+t)**6 +T_infinito
ax2.plot(t,T_6,label="<NAME>")
ax2.scatter(t,T)
plt.errorbar(garrafa_aberta['Tempo2 (s)'][1:6],garrafa_aberta['Temperatura2 (K)'][1:6],xerr=np.array([15.29978213,13.85640646,16.16580754,10.96965511,2.020725942]),yerr=1,fmt='none',label='Incerteza',ecolor='black',elinewidth=2)

#Calorímetro Aberto
#T_6_cal_aberto=(a_6_cal_aberto**6)/(-b_6_cal_aberto+t_cal_aberto)**6 +T_infinito
#ax2.plot(t_cal_aberto,T_6_cal_aberto,label="Calorímetro Aberto")
#ax2.scatter(t_cal_aberto,T_cal_aberto)

#Calorímetro Fechado
#T_6_cal_fechado=(a_6_cal_fechado**6)/(-b_6_cal_fechado+t_cal_fechado)**6 +T_infinito
#ax2.plot(t_cal_fechado,T_6_cal_fechado,label="Calorímetro Fechado")
#ax2.scatter(t_cal_fechado,T_cal_fechado)

#Garrafa Térmica Fechada
#T_6_gar_fechada=(a_6_gar_fechada**6)/(-b_6_gar_fechada+t_garrafa_fechada)**6 +T_infinito
#ax2.plot(t_garrafa_fechada,T_6_gar_fechada,label="Garrafa Térmica Fechada")
#ax2.scatter(t_garrafa_fechada,T_garrafa_fechada)

#Garrafa de Vidro
#T_6_vidro=(a_6_vidro**6)/(-b_6_vidro+t_vidro)**6 +T_infinito
#ax2.plot(t_vidro,T_6_vidro,label="Garrafa de Vidro")
#ax2.scatter(t_vidro,T_vidro)


plt.xlabel('Tempo (s)')
plt.ylabel('Temperatura (K)')
plt.legend(fontsize=8, loc='upper right')
plt.grid()
#plt.xticks(range(0,250,20))
#plt.yticks(np.arange(330,343,0.9))

plt.show()
#garrafa térmica aberta
print(a_6)
print(b_6)

#calorímetro Fechado
print(a_6_cal_fechado)
print(b_6_cal_fechado)
#<NAME>
print(a_6_gar_fechada)
print(b_6_gar_fechada)
#<NAME>
print(a_6_vidro)
print(b_6_vidro)<jupyter_output><empty_output>
{'fraction_non_alphanumeric': 0.16706833002551744, 'fraction_numerical': 0.07960022682166147, 'mean_word_length': 8.084996780424984, 'pattern_counts': {'":': 0, '<': 58, '<?xml version=': 0, '>': 62, 'https://': 4, 'lorem ipsum': 0, 'www.': 0, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 3, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
{'chain_length': '9', 'dir': 'jupyter-structured-clean-dedup', 'id': '65524', 'license': 'MIT', 'n_tokens_mistral': 7933, 'n_tokens_neox': 7387, 'n_words': 600, 'path': 'Experimento_de_Resfriamento.ipynb', 'repo_name': 'willian-gabriel-unicamp/F-229', 'size': '14419'}