poffit.m

生成系数矩阵

function hh=poffit(t1,t2)
k=0;%圆滑因子!

p=t1(:,1:2);
p=p';
h=(t1(:,3)-t1(:,4))';
pp=t2(:,1:2);
pp=pp';
[pn,minp,maxp,hn,minh,maxh]=premnmx(p,h);%归一化处理
%P，T分别为原始输入和输出数据，minp和maxp分别为P中的最小值和最大值。mint和maxt分别为T的最小值和最大值。
%premnmx函数用于将网络的输入数据或输出数据进行归一化，归一化后的数据将分布在[-1,1]区间内。
n=length(pn);
m=n;% % % m:结点个数

b=zeros(n,m);
for i=1:n
    for j =1:m
        b(i,j)=sqrt((pn(1,i)-pn(1,j))^2+(pn(2,i)-pn(2,j))^2+k^2);
    end
end
x=inv(b)*hn';

pp1=tramnmx(pp,minp,maxp);%归一化处理
nn=length(pp1);
p=zeros(nn,m);
for i=1:nn
    for j =1:m
        p(i,j)=sqrt((pp1(1,i)-pp1(1,j))^2+(pp1(2,j)-pp1(2,j))^2+k^2);
    end
end
hh=p*x;
hh=postmnmx(hh,minh,maxh);