rbf.m

本代码为基于matlab的求解偏微分方程的利用径向基函数的无网格方法的源码

function y =  rbf (rr)            % standard rbfs
global RBFscale;
global RBFpar;
global RBFtype;   % 1=MQ, 2=G, 3=IMQ, 4=W2, 5=TPS, 6=W4
r=rr/RBFscale;
if RBFtype==1
    y=-(1+r.*r).^(1/2);             % Multiquadric
    return;
end
if RBFtype==2
    y=exp(-r.*r);                  % Gaussian
    return;
end
if RBFtype==3
    y=(1+r.*r).^(-1/2);            % inverse Multiquadric
    return;
end
if RBFtype==4
  y=(1+4*r).*max(0,(1-r)).^4;      % Wendland 2
   return;
end
if RBFtype==5
    y=r.*r.*log(r+eps);                % TPS
    return;
end
if RBFtype==6
    y=(3+18*r+35*r.^2).*max(0,(1-r)).^6;      % Wendland 4
    return;
end
if RBFtype==7
    y=log(sqrt(r.^2+RBFpar.^2));      % regularized log = AFS of Laplace in R^2
    return;
end
error('RBF type not implemented')