3462403612

Jacek Złydach (JW2)

Metody Numeryczne - Zadanie Dodatkowe III

Interpolacja krzywymi sklejanymi stopnia drugiego (SPLINE-2)

Implementacja praktyczna

Wstęp

Poniższa praktyczna implementacja stanowi uzupełnienie teoretycznych rozważań na temat interpolacji za pomocą krzywych sklejanych drugiego stopnia. Wykorzystuje ona kod omówiony dokładnie w III sprawozdaniu z Metod Numerycznych (metody numerycznego rozwiązywania układów równań liniowych).

Funkcje

Poniżej umieściłem garść funkcji pomocniczych dokładnie omówionych w trzecim sprawozdaniu (między innymi faktoryzacja LU).

sklej_macierze_poziomo( lhs, rhs) :=

lhs

return (—1) if rows(lhs) * rows(rhs) for ie 0..cols(rhs) - 1 (cols(lhs)+i) jj^i)

return lhs

wydziel_macierz_lewa(lhs, leftCols) :=

for i e 0.. leftCols - 1

leftMatrix^<' <— lhs^ return leftMatrix

wydziel_macierz_prawa(lhs, leftCols) :=

podstawienia_U(A,b) :=

last(b) a

last(b),last(b) for i e last(b) - 1.. 0 last(b)

S (^Ai.j ■ ’ „ , j=i+l

for i € leftCols.. cols(lhs) - 1

• u.n/i . • (i-leftCols) ,    (i)

nghtMatrix    <— lhs

return rightMatrix ^last(b)

str. 1/19