PythonLabFisica/LeiDeHooke.py at master · DanielSenaDev/PythonLabFisica · GitHub

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
from cmath import sqrt
import math
from xml.etree.ElementTree import tostring


# Funcoes

# Calcular o valor médio
def calculo_valorMedio (lista_x, n):

    soma_dos_x = 0
    for x in lista_x:
        soma_dos_x = soma_dos_x+x
    valor_medio = soma_dos_x/n

    print("Soma dos x: ",soma_dos_x)
    print("Valor médio: ",valor_medio)
    print("")
    return valor_medio

# Calcular Desvio Padrão da Medida - DPM
def calculo_dpm (lista_x, valor_medio):
    somatorio_dpm = 0
    for x in lista_x:
        valor_elemento = (x-valor_medio)**2
        somatorio_dpm = somatorio_dpm+valor_elemento

    divisao_dpm = somatorio_dpm/(n-1)
    raiz_dpm = (divisao_dpm)**(1/2)
    dpm = raiz_dpm

    print("Somatorio DPM: ",somatorio_dpm)
    print("Divisão DPM: ",divisao_dpm)
    print("Raiz DPM: ",raiz_dpm)
    print("DPM: ",dpm)
    print("")
    return dpm

#Calcular o sigmaA = DPM/raiz(n)
def calculo_sigmaA (dpm, n):
    sigmaA = dpm/(n**(1/2))

    print("SigmaA: ",sigmaA)
    print("")
    return sigmaA


# Calcular Incerteza Total: sigmaT
def calculo_sigmaT (sigmaA, sigmaB):
    print("Calculo da Incerteza Total: sigmaT")

    sigmaT = (sigmaA**2+sigmaB**2)**(1/2)

    print("Sigma T: ", sigmaT)
    print("")
    return sigmaT

# Calcular Volume do Cilindro
def calculo_VolumeMedio (dMedio, LMedio):
    volumeMedio = (math.pi * (dMedio**2) * LMedio)/4
    print("Volume Médio: ",volumeMedio)
    return volumeMedio


##################### Valores constantes:


# Valor da Incerteza Instrumental - sigmaB
sigmaB = 0.0005


# Valor da Massa:
massa = 0.01

# Quantidade de medidas
#n = int(input("Quantidade de medições (n): "))
n = 3

# Valores medidos de D:
lista_a = [0.260, 0.260, 0.260]

# Valores medidis de L:
#lista_l = [9.25,9.21,9.22,9.26,9.30,9.30,9.25,9.25,9.25,9.25]


##########################
# Funcionamento do Programa

print("\nSigmaB do Instrumento: ", sigmaB)
print("Valor da massa: ", massa)
print("Quantidade colhida: ",n)
print("Valores de A: ",lista_a)
#print("Valores de L: ",lista_l)
print("")


print("\n====\nAgora os outros Calculos\n")

valorMedio = calculo_valorMedio(lista_a, n)
print(valorMedio)