curfit3.m

生成系数矩阵

function hh=curfit3(t1,t2)
p=t1(:,1);
p=p';
pp=t2(:,1);
pp=pp';
h=(t1(:,2)-t1(:,3))';
[pn,minp,maxp,hn,minh,maxh]=premnmx(p,h);%归一化处理
%P，T分别为原始输入和输出数据，minp和maxp分别为P中的最小值和最大值。mint和maxt分别为T的最小值和最大值。
%premnmx函数用于将网络的输入数据或输出数据进行归一化，归一化后的数据将分布在[-1,1]区间内。
% b=son(pn',1);
% b1=b'*b;
% b2=b'*hn';
% a=inv(b1)*b2;
a=polyfit(pn,hn,3);
pp1=tramnmx(pp,minp,maxp);%归一化处理

%用plot做二维图象!-----------------------------
figure(2)
hold off
xx1=[-1.0:0.05:1.0];
yy1=a(4)+a(3)*xx1+a(2)*xx1.^2+a(1)*xx1.^3;
xx=postmnmx(xx1,minp,maxp);
yy=postmnmx(yy1,minh,maxh);
plot(p,h,'*')
hold on
plot(xx,yy,'-r')
%---------------------------------

hh=a(4)+a(3)*pp1+a(2)*pp1.^2+a(1)*pp1.^3;
hh=postmnmx(hh,minh,maxh)';