hm2-cheat-sheet

Splineinterpolation

Das Gleichungssystem unter Punkt 4 hat dabei folgende Form:

Dies kann in Python folgendermaßen generiert werden:

import numpy as np

def calculate_matrix(x):
    n = len(x) - 1
    h = x[1:] - x[:-1]
    A = np.zeros((n - 1, n - 1))

    for i in range(n - 1):
        A[i, i] = 2 * (h[i] + h[i + 1])

    for i in range(0, n - 2):
        A[i, i + 1] = h[i]
        A[i + 1, i] = h[i]

    return A, h

Beispiel