# Create B-Spline curve using MATLAB

Praba89
I need a MATLAB expert to guide me on how to create a b-spline curve using MATLAB Software. I understand the B-spline basis function calculations for zeroth and first degree but I have no idea on how to calculate for the 2nd degree. I need a favor on that part. I am currently working on my final year project titled INTERPOLATION OF PLANAR CURVE WITH DIFFERENT PARAMETERISATION. Thanks in Advance.

Last edited:

## Answers and Replies

timthereaper
Just like you linearly interpolate the 0th degree basis functions to get the 1st degree basis functions, you have to linearly interpolate the 1st degree basis functions to get the 2nd degree. The recursive B-spline basis formula tells you how to do that.

Here's some code snippets for MATLAB does B-spline curve evaluation for a degree-p B-spline curve with knot vector U at parameter value u with control point array C:

Code:
function S = bcurve_eval(u,p,U,C)
uspan = findKnotSpan(u,p,U);
Nu = getBasisFuncs(u,p,U,uspan);

S = zeros(1,numel(C(1,:)));
for i=0:p
index = uspan-p+i;
S = S + C(index+1,:) * Nu(i+1);
end
end

Code:
function i = findKnotSpan(u,p,U)
k = numel(U)-1;
n = k-p-1;
if u == U(n+1 +1)
i = n;
return
end
l = p +1;
h = n+1 +1;
i = floor((l+h)/2);
while u < U(i) || u >= U(i+1)
if u < U(i)
h = i;
else
l = i;
end;
i = floor((l+h)/2);
end
i = i-1;
end

Code:
function N = getBasisFuncs(u,p,U,i)
left = zeros(1,p+1);
right = zeros(1,p+1);
N = zeros(1,p+1);
N(1) = 1;
for j=1:p
left(j+1) = u - U(i+1-j +1);
right(j+1) = U(i+j +1) - u;
s = 0;
for r=0:j-1
temp = N(r +1) / (right(r+1 + 1) + left(j-r +1));
N(r +1) = s + right(r+1 +1) * temp;
s = left(j-r +1) * temp;
end
N(j +1) = s;
end
end

Hope that helps. Let me know if you have questions.

• Yaaqob Saad and (deleted member)