:Ȩϵжάֵ(Weight Function)

ʽ:
InterpolationWF2D(x0,y0,f0,x,y,n,a,k,t)
InterpolationWF2D(x0,y0,f0,x,y,n,a,k)
InterpolationWF2D(x0,y0,f0,x,y,n,a)
InterpolationWF2D(x0,y0,f0,x,y,n)
InterpolationWF2D(x0,y0,f0,x,y)

x0:洢ɢx
y0:洢ɢy,˲x0άһ
f0:洢ɢ(x,y)¶Ӧĺֵ,,˲x0y0άһ
x:洢Ҫֵx
y:洢Ҫֵy,˲xάһ
n:ѡȡֵ֪,ĬΪ4
a:˥ϵ,һĬΪ1tѡ3ʱ,ֵֻܽ(0,1)֮
k:һ1*2ľ,Ԫطֱʾx0y0Բֵf0Ȩϵ,ĬϾһ,ȡ1ע,ȨϵС0.
t:Ȩϵ,ɲο˵.ֵĬΪ1


˵:

1ɹxyͬάĲֵֵ,ɳԸӦƽȲ,ʹú{SmoothGrid<\SmoothGrid>}Խмӹ

2tѡȡʱӦ
t=1, f(r)=exp(-a*r)/R
t=2, f(r)=(1-r)^a       //ѡģʱ,齫nȡ8
t=3, r<a,f(r)= 1+(0.1-1)*r/a ;; ra,f(r)=0.1*(1-r)/(1-a)  //עѡʱ,aֵ(0,1)֮
t=4, f(r)=4/(3-exp(1-r^a))
t=5, f(r)=sin((1-r)^a)
t=6, f(r)=sin(1-r)^a

ע:ֵıֵڸֵֵ֪֮䡣

ԭ:{ɢݲֵԤ<http://blog.163.com/shikang999@126/blog/static/17262489620150193586584/>}

:

//Сľ{ά<http://emuch.net/bbs/viewthread.php?tid=5999724&fpage=1>}Ϊ,data,λسִиɿ


x=data(,1);//data1
y=data(,2);
z=data(,3);
a=max(x);
b=min(x);
x0=MeshGrid2(b,a,50);//x
a=max(y);
b=min(y);
y0=MeshGrid2(b,a,50);//y
y0=y0';//עҪתһ
z0=InterpolationWF2D(x,y,z,x0,y0);//ñ
cloud(x0,y0,z0)//ԻƳͼ1Ч

z0=InterpolationWF2D(x,y,z,x0,y0,8,1,2);cloud(x0,y0,z0)//ԻƳͼ2Ч
//ͼ1
{<http://img0.ph.126.net/ewg8NYE2lex6DU5b2BqrxA==/6630917132605487004.png>}
//ͼ2
{<http://img2.ph.126.net/sBCJi-meVHwo0rL4fPwINw==/6630149673489293162.png>}