Στα natural splines στο βιβλίο του κ.Γραββάνη , δημιουργείται ένα σύστημα Ax = b ,όπου
Για να επιλύσει το παραπάνω σύστημα χρησιμοποιεί τον εξής κώδικα στο matlab:
- Κώδικας: Επιλογή όλων
clear;
clc;
x = [...] ;
a = [...];
x0 = ...;
n = length(x) - 1 ;
for i = 0 : (n-1)
h(i+1) = x(i+2) - x(i+1) ;
end
b= zeros( n+1 , 1 ) ;
for i =2: n
b(i,1) = 3 *(a(i+1)-a(i))/h(i) - 3/h(i-1)*(a(i) - a(i-1) ) ;
end
%epilisi grammikou sistimatos
l(1) =0 ; m(1) = 0 ; z(1) = 0 ;
for i =1:(n-1)
l(i+1) = 2*( x(i+2) - x(i) ) - h(i)* m(i) ;
m(i+1) = h(i+1)/l(i+1);
z(i+1) = ( s(i+1) - h(i)*z(i) ) / l ( i+1) ;
end
l(n+1) =1;
z(n+1) = 0 ;
c(n+1) = 0 ;
for j = ( n-1) : (-1) : 0
c(j+1) = z(j+1) - m(j+1)*c(j+2) ;
end
Δεν μπορώ να καταλάβω όμως ποία μέθοδο χρησιμοποιεί από τις άμεσες μεθόδους επίλυσης γραμμικών συστημάτων....'Εχει κανείς καμία ιδέα;
) ας στείλει ενα pm.