def probabilidadPropios(vectorEstado, vectorPropio):
prob = []
bra = cM.adjMatrix(vectorEstado)[0]
for i in range(len(bra)):
bra[i] = [(bra[i])]
for j in range(len(vectorPropio)):
prob += [((compl.modulo(
cM.productMatrix(cM.transMatrix(bra), vectorPropio)[0][0]))**2, 0)]
return prob
def esta_en_mandelbrot(c):
z = [0, 0]
cont = 0
for i in range(30):
z = funcion_f(z, c)
cont += 1
#print(z)
if comp.modulo(z) > 2:
return False, cont
return True, 0
def particula_en_una_recta(n, vi):
"""
Retorna un vector de probabilidades P.
Parametros:
n: Número de puntos sobre la recta.
vi: Vector de estado inicial.
"""
P = [0 for i in range(n)]
modulo_v = c.modulo_vector(vi)
for i in range(n):
P[i] = round((c.modulo(vi[i]) / modulo_v)**2, 4)
#Graficación de resultados
return P
def probabilidadPosicion(posicion, vectorKet):
return (compl.modulo(vectorKet[posicion][0])**
2) / cM.normaVec(vectorKet)**2
def probabilidad(c):
return compl.modulo(c)**2
def test_modulo():
assert c.modulo([1, 1]) == 1.4142135623730951, 'Debe ser 1.4142'
def modulomatriz(matriz):
for i in range(len(matriz)):
for j in range(len(matriz[0])):
matriz[i][j] = [comp.modulo(matriz[i][j])**2, 0]
return matriz
def valorreal(matriz):
for i in range(len(matriz)):
for j in range(len(matriz[0])):
matriz[i][j] = [comp.modulo(matriz[i][j]), 0]
return matriz
def test_modulo(self):
pc = (4, -2)
res = 4.47213595499958
rescod = complejos.modulo(pc)
self.assertEqual(res, rescod)
