computeCubicSplineCoefficients.m
041b8992f8b7master

computeCubicSplineCoefficients.m
View Options

			function [ c ] = computeCubicSplineCoefficients(s, halfBoundary)
			% computeCubicSplineCoefficients(s, halfBoundary) applies the causal and
			% anticausal filters in order to obtain the coefficients c[k].

			n = length(s);
			c = zeros(1, n);
			h = 1 / (length(s) - 1 - 2 * halfBoundary);
			% c0 = 6 / (sqrt(3) + 2); %%%%% WHAT THE FUCK????
			c0 = 6.0;
			a = (sqrt(3.0) - 2.0);
			tempo = zeros(1, n);

			initialValueCausal = computeInitialValueCausal(s, a);
			tempo(1) = initialValueCausal;
			for i=2:n
			tempo(i) = s(i) + a * tempo(i-1);
			end

			initialValueAntiCausal = computeInitialValueAntiCausal(tempo, a);
			% c(n) = tempo(n) + a * initialValueAntiCausal;
			%
			% for i=n-1:-1:1
			% c(i) = tempo(i) + a * c(i+1);
			% end
			c(n) = initialValueAntiCausal;
			for i=n-1:-1:1
			c(i)=a*(c(i+1)-tempo(i));
			end
			c = c0 * c;
			end