我知道这种类型的问题之前已经被问过很多次了,但我无法理解我的代码的问题。我也是 Octave 的绝对初学者。
function Z=nat_spline(n, T, Y)
for i=0:n-1
H(i) = T(i+1) - T(i);
B(i) = 6*(Y(i+1) - Y(i))/H(i);
end
U(1) = 2*(H(0)+H(1));
V(1) = B(1) - B(0);
for i=2:n-1
U(i) = 2 * (H(i) + H(i-1)) - ( (H(i-1))^2 / U(i-1) );
V(i) = B(i) - B(i-1) - H(i-1)*V(i-1)/U(i-1);
end
Z(n) = 0
for i=n-1:1
Z(i) = (V(i)-H(i)*Z(i+1))/U(i);
end
Z(0) = 0;
end
它的代码很短,所以我想很容易发现任何错误。非常感谢您的帮助。