Gegeven is een reeks van
We wensen een functie te fitten aan deze punten, en daartoe stellen we in elk interval
De functie
De coëfficiënten
waarbij de eerste 3 vergelijkingen gelden voor
Schrijf de functie spline3()
met als argumenten de 1D NumPy- rijen x
en y
, elk bestaande uit N
+1 reële getallen. Het resultaat van de functie is de 1D NumPy-rij, met [b0, c0, d0, b1, c2, d1, ...]
.
Het is hierbij de bedoeling dit stelsel volledig te construeren en vervolgens op te lossen. Gebruik zo veel mogelijk de geziene ingebouwde NumPy-functionaliteit.
x = np.array([-7.0, -6.0, -5.0, -4.0, -3.0, -2.0]) y = np.array([5.089685242514852, 2.4208420859231627, 10.82540893401563, 6.503962993840951, 2.1015572909559777, 9.552280583187677]) s = spline3(x, y) print(s) #[ -6.49096844 0. 3.82212528 4.97540741 11.46637585 # -8.03721641 3.79650988 -12.64527338 4.52731756 -7.9120842 # 0.93667929 2.57299921 1.68027201 8.65567692 -2.88522564]