Oppgave 15.01 ballfart1.py

# Program ballfart1.py
import numpy as np
import matplotlib.pyplot as plt

# Leser t- og s-verdier inn i lister og lager v-liste
data = np.loadtxt("ballslipp.txt") # tekstfil med t og s
t = data[:,0] # innleste t-verdier [s] 
s = data[:,1] # innleste s-verdier [m] 
n = len(s) # lengde av t-liste
v = np.zeros(n) # v-liste [m/s]
t1 = np.zeros(n) # t1-liste [s]

# Stegvis beregning av v
for i in range(1,n): 
  v[i] = (s[i]-s[i-1])/(t[i]-t[i-1]) # farten beregnes
  t1[i] = (t[i]+t[i-1])/2

# Skriver ut figur
plt.figure(1)
plt.plot(t1,v)
plt.grid()
plt.title('Fart')
plt.xlabel('$t$ / s')
plt.ylabel('$v$ / m/s')
plt.show()