www.pudn.com > zhichixiangliangjijiemian.rar > CUBWGT.M

function [w,dwdx,dwdy]=cubwgt(dif,t,v,dmax,dm) 
%  CUBIC  SPLINE  WEIGHT  FUNCTION 
l=length(v); 
for i=1:l 
drdx=sign(dif(1,i))/dm(1,v(i)); 
drdy=sign(dif(2,i))/dm(2,v(i)); 
rx=abs(dif(1,i))/dm(1,v(i)); 
ry=abs(dif(2,i))/dm(2,v(i)); 
if rx>1 
wx=0; 
dwx=0; 
elseif rx<=1 
wx=1-6*rx*rx+8*rx^3-3*rx^4; 
dwx=(-12*rx+24*rx^2-12*rx^3)*drdx; 
end 
if ry>1 
wy=0; 
dwy=0; 
elseif ry<=1 
wy=1-6*ry*ry+8*ry^3-3*ry^4; 
dwy=(-12*ry+24*ry^2-12*ry^3)*drdy; 
end 
w(i)=wx*wy; 
dwdx(i)=wy*dwx; 
dwdy(i)=wx*dwy; 
end