Eksempel 15.1 ballfart.py

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

# Leser t- og s-verdier inn i lister og lager en tom v-liste
data = np.loadtxt('ballslipp.txt') # laster inn tekstfil
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]

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

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