niherf.m

轧制力拟合程序 基本思想是线性方程加BP补偿

%数据拟合
%板宽B＝1000
clc;
clear;
load d:\matlab\work\RF\nihedata.txt;
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%   SK,DIR,DWR,SHFI,SHFW,B  %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
a1=zeros(640,72);
for i=1:640
    a1(i,1)=nihedata(i,2);
    a1(i,2)=nihedata(i,3);
    a1(i,3)=nihedata(i,5);
    a1(i,4)=nihedata(i,6);
    a1(i,5)=nihedata(i,2)^2;
    a1(i,6)=nihedata(i,3)^2;
    a1(i,7)=nihedata(i,5)^2;
    a1(i,8)=nihedata(i,6)^2;

    a1(i,9)=nihedata(i,1)*nihedata(i,2);
    a1(i,10)=nihedata(i,1)*nihedata(i,3);
    a1(i,11)=nihedata(i,1)*nihedata(i,5);
    a1(i,12)=nihedata(i,1)*nihedata(i,6);
    a1(i,13)=nihedata(i,1)*nihedata(i,2)^2;
    a1(i,14)=nihedata(i,1)*nihedata(i,3)^2;
    a1(i,15)=nihedata(i,1)*nihedata(i,5)^2;
    a1(i,16)=nihedata(i,1)*nihedata(i,6)^2;

    a1(i,17)=nihedata(i,4)*nihedata(i,2);
    a1(i,18)=nihedata(i,4)*nihedata(i,3);
    a1(i,19)=nihedata(i,4)*nihedata(i,5);
    a1(i,20)=nihedata(i,4)*nihedata(i,6);
    a1(i,21)=nihedata(i,4)*nihedata(i,2)^2;
    a1(i,22)=nihedata(i,4)*nihedata(i,3)^2;
    a1(i,23)=nihedata(i,4)*nihedata(i,5)^2;
    a1(i,24)=nihedata(i,4)*nihedata(i,6)^2;

    a1(i,25)=nihedata(i,7)*nihedata(i,2);
    a1(i,26)=nihedata(i,7)*nihedata(i,3);
    a1(i,27)=nihedata(i,7)*nihedata(i,5);
    a1(i,28)=nihedata(i,7)*nihedata(i,6);
    a1(i,29)=nihedata(i,7)*nihedata(i,2)^2;
    a1(i,30)=nihedata(i,7)*nihedata(i,3)^2;
    a1(i,31)=nihedata(i,7)*nihedata(i,5)^2;
    a1(i,32)=nihedata(i,7)*nihedata(i,6)^2;

    a1(i,33)=nihedata(i,8)*nihedata(i,2);
    a1(i,34)=nihedata(i,8)*nihedata(i,3);
    a1(i,35)=nihedata(i,8)*nihedata(i,5);
    a1(i,36)=nihedata(i,8)*nihedata(i,6);
    a1(i,37)=nihedata(i,8)*nihedata(i,2)^2;
    a1(i,38)=nihedata(i,8)*nihedata(i,3)^2;
    a1(i,39)=nihedata(i,8)*nihedata(i,5)^2;
    a1(i,40)=nihedata(i,8)*nihedata(i,6)^2;
    
    a1(i,41)=nihedata(i,1)^2*nihedata(i,2);
    a1(i,42)=nihedata(i,1)^2*nihedata(i,3);
    a1(i,43)=nihedata(i,1)^2*nihedata(i,5);
    a1(i,44)=nihedata(i,1)^2*nihedata(i,6);
    a1(i,45)=nihedata(i,1)^2*nihedata(i,2)^2;
    a1(i,46)=nihedata(i,1)^2*nihedata(i,3)^2;
    a1(i,47)=nihedata(i,1)^2*nihedata(i,5)^2;
    a1(i,48)=nihedata(i,1)^2*nihedata(i,6)^2;

    a1(i,49)=nihedata(i,4)^2*nihedata(i,2);
    a1(i,50)=nihedata(i,4)^2*nihedata(i,3);
    a1(i,51)=nihedata(i,4)^2*nihedata(i,5);
    a1(i,52)=nihedata(i,4)^2*nihedata(i,6);
    a1(i,53)=nihedata(i,4)^2*nihedata(i,2)^2;
    a1(i,54)=nihedata(i,4)^2*nihedata(i,3)^2;
    a1(i,55)=nihedata(i,4)^2*nihedata(i,5)^2;
    a1(i,56)=nihedata(i,4)^2*nihedata(i,6)^2;

    a1(i,57)=nihedata(i,7)^2*nihedata(i,2);
    a1(i,58)=nihedata(i,7)^2*nihedata(i,3);
    a1(i,59)=nihedata(i,7)^2*nihedata(i,5);
    a1(i,60)=nihedata(i,7)^2*nihedata(i,6);
    a1(i,61)=nihedata(i,7)^2*nihedata(i,2)^2;
    a1(i,62)=nihedata(i,7)^2*nihedata(i,3)^2;
    a1(i,63)=nihedata(i,7)^2*nihedata(i,5)^2;
    a1(i,64)=nihedata(i,7)^2*nihedata(i,6)^2;

    a1(i,65)=nihedata(i,8)^2*nihedata(i,2);
    a1(i,66)=nihedata(i,8)^2*nihedata(i,3);
    a1(i,67)=nihedata(i,8)^2*nihedata(i,5);
    a1(i,68)=nihedata(i,8)^2*nihedata(i,6);
    a1(i,69)=nihedata(i,8)^2*nihedata(i,2)^2;
    a1(i,70)=nihedata(i,8)^2*nihedata(i,3)^2;
    a1(i,71)=nihedata(i,8)^2*nihedata(i,5)^2;
    a1(i,72)=nihedata(i,8)^2*nihedata(i,6)^2;


end

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
xxx=a1;


load d:\matlab\work\RF\RF.txt;


%coeff=xxx\y;
kkk=pinv(xxx)*RF;



save D:\matlab\work\RF\K_result\KKK.txt kkk -ASCII;


t=0;
yy=zeros(640,1);
for i=1:640
    for j=1:72
        t=xxx(i,j)*kkk(j,1);
        yy(i,1)=yy(i,1)+t;
    end
end