BiCubic基函数matlab代码

clear all; close all; clc;
A=-5:0.005:5;
a = -0.5;%自由参数
B=zeros(size(A));
for i = 1 : length(A)
    if (abs(A(i)) <= 1)
        B(i) = (a+2)*abs(A(i))^3 - (a+3)*abs(A(i))^2 + 1;
    elseif (abs(A(i)) < 2)
        B(i) = a*abs(A(i))^3 - 5*a*abs(A(i))^2 + 8*a*abs(A(i)) - 4*a;
    else
        B(i) = 0;
    end
end
plot(A,B);
axis([-5 5 -0.2 1]);
xlabel('x')
ylabel('w(x)')