PYTHON   24
battements py
Guest on 14th March 2023 01:10:06 PM


  1. import numpy as np
  2. import matplotlib.pyplot as plt
  3. import matplotlib.animation as animation
  4. import time
  5.  
  6. fig = plt.figure(figsize=(12,4))
  7. ax = fig.add_subplot()
  8. #tmax = 0.0002
  9. tmax = 0.0018
  10. c = 350
  11.  
  12. ech = 10000
  13. t = np.linspace (0, tmax, ech)
  14. E1 = np.linspace(0, 0, ech)
  15. E2 = np.linspace(0, 0, ech)
  16. f1 = 40000.
  17. f2 = 42000.
  18. w1 = 2. * np.pi * f1
  19. w2 = 2. * np.pi * f2
  20. for e in range(ech):
  21.     E1[e]=np.cos(w1*t[e])
  22.     E2[e]=np.cos(w2*t[e])
  23.  
  24. #plt.plot(t, E1)
  25. #plt.plot(t, E2)
  26. plt.plot(t, E1 + E2)
  27.  
  28. #print("f1 = ", f1, "(s-1)")
  29. #print("f2 = ", f2, "(s-1)")
  30.  
  31. plt.show()
  32.  
  33.  
  34. #Calcul de la vitesse de l'onde
  35. d = 0.0425
  36. n = 5
  37. la = d/n
  38. c = la * f1
  39. s_exp = c * np.sqrt(np.square(0.001/d))
  40. print("vitesse de l'onde : ", c, "(m.s-1)")
  41. print("incertitude sur la vitesse de l'onde : ", s_exp, "(m.s-1)")
  42.  
  43.  
  44. #calcul vitesse du train
  45. sT = 0.00024
  46. fm = f2 - f1
  47. vitesse_emetteur = c*fm/(fm+f1)
  48. print("vitesse de l'Ć©metteur = " ", vitesse_emetteur, " (m.s-1)")
  49. sv_exp = vitesse_emetteur*np.sqrt(np.square(sT*fm)+ np.square(s_exp/c))
  50. print("incertitude sur la vitesse de l'Ć©metteur : ", sv_exp, "(m.s-

Raw Paste

Login or Register to edit or fork this paste. It's free.