test.m

来自「最邻近的RBF算法,MATLAB实现,较为有价值的资料.」· M 代码 · 共 27 行

clc
clear
p=0:0.1:5
t=sqrt(p);
r=0.2;
[c ,w ,nr ]=zjwnnc (p ,t ,r)


x=0:0.2:5
n=size(x,2);
sumd=zeros(1,n);
sumdd=zeros(1,n);
testout=zeros(1,n);
for j=1:1:n
    for i=1:1:nr
    d=(w(:,i)*exp(-(sqrt(sum((x(j)-c(:,i)).^2)))^2/r^2));
    sumd(:,j)=sumd(:,j)+d;
    dd=(exp(-(sqrt(sum((x(j)-c(:,i)).^2)))^2/r^2));
    sumdd(:,j)=sumdd(:,j)+dd;
    end
    testout(:,j)=sumd(:,j)/sumdd(:,j);
end

figure
hold on
plot(p,t,'k*')
plot(x,testout,'ro')