NV in Wild

Refath Bari

B_1(x,t)

B_2(x,t)

B_1(x,t)

B_2(x,t)

B_3(x,t)

B_1(x,t)

B_2(x,t)

B_3(x,t)

B_4(x,t)

def B1(t): 
  w = 4
  A = 4
  phShift = np.rand(0,1)
  return A*math.sin(w*t+phShift)

B_1(x,t)

def B1(t): 
  w = 4
  A = 4
  phShift = 1
  return A*math.sin(w*t+phShift)

def B2(t): 
  w = 4
  A = 4
  phShift = 0
  return A*math.cos(w*t)

def B3(t): 
  w = 4
  A = 6
  phShift = 0
  return A*math.sin(w*t)

def B4(t): 
  w = 4
  A = 8
  phShift = 0
  return A*(math.sin(w*t-0.4))**3

B_1(x,t)

B_2(x,t)

B_3(x,t)

B_4(x,t)

def total_B(t): 
  return B1(t)+B2(t)+B3(t)+B4(t)

f(0)=f(T)

-\int_{0}^{0.3} B(\alpha) d\alpha+\int_{0.3}^{0.6} B(\alpha) d\alpha

-\int_{0}^{0.6} B(\alpha) d\alpha+\int_{0.6}^{1.2} B(\alpha) d\alpha

-\int_{0}^{1} B(\alpha) d\alpha+\int_{1}^{2} B(\alpha) d\alpha

-\int_{0}^{1.6} B(\alpha) d\alpha+\int_{1.5}^{3.2} B(\alpha) d\alpha

-\int_{0}^{\tau} B(\alpha) d\alpha+\int_{\tau}^{2\tau} B(\alpha) d\alpha

def phi_accum(field, a, b): 
    B_int = quad(field, a, b, limit=100000)
    return B_int[0]

def phi_accum_prime(field, t_b):
  return ( -1 * phi_accum(field, 0, t_b) + phi_accum(field, t_b, 2*t_b) ) 

def HE_coherence(field, t_b):
  return math.cos(phi_accum(field, 0, t_b))

for i in x: 
  b1_graph.append(B1(i)) #REPEAT FOR OTHER 3 MAGNETIC FIELDS
  sum_graph.append(total_B(i))
  
  y1.append(phi_accum_prime(B1, i)) # REPEAT FOR OTHER 3
  y1_mod.append(math.cos(phi_accum_prime(B1, i))) # REPEAT FOR OTHER 3
  y1_free_coh.append(free_coherence(B1, i)) # REPEAT FOR OTHER 3