29 lines
681 B
Python
29 lines
681 B
Python
import numpy as np
|
|
import matplotlib.pyplot as plt
|
|
|
|
# Parametry ruchu harmonicznego
|
|
A_v = 0.6 # Amplituda prędkości (m/s)
|
|
T = 0.4 # Okres (s)
|
|
f = 1 / T # Częstotliwość (Hz)
|
|
omega = 2 * np.pi * f # Pulsacja (rad/s)
|
|
C = 0 # Przyjmujemy C = 0 dla uproszczenia
|
|
|
|
# Funkcja x(t)
|
|
def x_t(t):
|
|
return (-A_v / omega) * np.cos(omega * t) + C
|
|
|
|
# Przedział czasu
|
|
t = np.linspace(0, 1, 1000) # 1 sekunda ruchu
|
|
x = x_t(t)
|
|
|
|
# Wykres funkcji x(t)
|
|
plt.figure(figsize=(8, 5))
|
|
plt.plot(t, x, label="$x(t) = -0.0382 \cos(5\pi t) + C$", color='b')
|
|
plt.xlabel("Czas [s]")
|
|
plt.ylabel("Położenie x(t) [m]")
|
|
plt.title("Wykres funkcji położenia x(t)")
|
|
plt.legend()
|
|
plt.grid()
|
|
plt.show()
|
|
|