%數據擬合
%板寬B＝1000
clc;
clear;
load d:\matlab\work\model\nihedata.txt;
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%   SK,DIR,DWR,SHFI,SHFW,B  %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%



%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%數據擬合a1
a1=zeros(5120,151);
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
for i=1:5120
    a1(i,1)=nihedata(i,1);
    a1(i,2)=nihedata(i,1)^2;
    a1(i,3)=nihedata(i,1)^3;
    a1(i,4)=nihedata(i,6);
    a1(i,5)=nihedata(i,6)^2;
    a1(i,6)=nihedata(i,6)^3;
    a1(i,7)=nihedata(i,6)^4;
    a1(i,8)=nihedata(i,1)*nihedata(i,6);
    a1(i,9)=nihedata(i,1)*nihedata(i,6)^2;
    a1(i,10)=nihedata(i,1)*nihedata(i,6)^3;
    a1(i,11)=nihedata(i,1)*nihedata(i,6)^4;
    a1(i,12)=nihedata(i,1)^2*nihedata(i,6);
    a1(i,13)=nihedata(i,1)^2*nihedata(i,6)^2;
    a1(i,14)=nihedata(i,1)^2*nihedata(i,6)^3;
    a1(i,15)=nihedata(i,1)^2*nihedata(i,6)^4;
    a1(i,16)=nihedata(i,1)^3*nihedata(i,6);
    a1(i,17)=nihedata(i,1)^3*nihedata(i,6)^2;
    a1(i,18)=nihedata(i,1)^3*nihedata(i,6)^3;
    a1(i,19)=nihedata(i,1)^3*nihedata(i,6)^4;
    
    a1(i,20)=nihedata(i,5)*nihedata(i,1);
    a1(i,21)=nihedata(i,5)*nihedata(i,1)^2;
    a1(i,22)=nihedata(i,5)*nihedata(i,1)^3;
    a1(i,23)=nihedata(i,5)*nihedata(i,6);
    a1(i,24)=nihedata(i,5)*nihedata(i,6)^2;
    a1(i,25)=nihedata(i,5)*nihedata(i,6)^3;
    a1(i,26)=nihedata(i,5)*nihedata(i,6)^4;
    a1(i,27)=nihedata(i,5)*nihedata(i,1)*nihedata(i,6);
    a1(i,28)=nihedata(i,5)*nihedata(i,1)*nihedata(i,6)^2;
    a1(i,29)=nihedata(i,5)*nihedata(i,1)*nihedata(i,6)^3;
    a1(i,30)=nihedata(i,5)*nihedata(i,1)*nihedata(i,6)^4;
    a1(i,31)=nihedata(i,5)*nihedata(i,1)^2*nihedata(i,6);
    a1(i,32)=nihedata(i,5)*nihedata(i,1)^2*nihedata(i,6)^2;
    a1(i,33)=nihedata(i,5)*nihedata(i,1)^2*nihedata(i,6)^3;
    a1(i,34)=nihedata(i,5)*nihedata(i,1)^2*nihedata(i,6)^4;
    a1(i,35)=nihedata(i,5)*nihedata(i,1)^3*nihedata(i,6);
    a1(i,36)=nihedata(i,5)*nihedata(i,1)^3*nihedata(i,6)^2;
    a1(i,37)=nihedata(i,5)*nihedata(i,1)^3*nihedata(i,6)^3;
    a1(i,38)=nihedata(i,5)*nihedata(i,1)^3*nihedata(i,6)^4;
    
    a1(i,39)=nihedata(i,5)^2*nihedata(i,1);
    a1(i,40)=nihedata(i,5)^2*nihedata(i,1)^2;
    a1(i,41)=nihedata(i,5)^2*nihedata(i,1)^3;
    a1(i,42)=nihedata(i,5)^2*nihedata(i,6);
    a1(i,43)=nihedata(i,5)^2*nihedata(i,6)^2;
    a1(i,44)=nihedata(i,5)^2*nihedata(i,6)^3;
    a1(i,45)=nihedata(i,5)^2*nihedata(i,6)^4;
    a1(i,46)=nihedata(i,5)^2*nihedata(i,1)*nihedata(i,6);
    a1(i,47)=nihedata(i,5)^2*nihedata(i,1)*nihedata(i,6)^2;
    a1(i,48)=nihedata(i,5)^2*nihedata(i,1)*nihedata(i,6)^3;
    a1(i,49)=nihedata(i,5)^2*nihedata(i,1)*nihedata(i,6)^4;
    a1(i,50)=nihedata(i,5)^2*nihedata(i,1)^2*nihedata(i,6);
    a1(i,51)=nihedata(i,5)^2*nihedata(i,1)^2*nihedata(i,6)^2;
    a1(i,52)=nihedata(i,5)^2*nihedata(i,1)^2*nihedata(i,6)^3;
    a1(i,53)=nihedata(i,5)^2*nihedata(i,1)^2*nihedata(i,6)^4;
    a1(i,54)=nihedata(i,5)^2*nihedata(i,1)^3*nihedata(i,6);
    a1(i,55)=nihedata(i,5)^2*nihedata(i,1)^3*nihedata(i,6)^2;
    a1(i,56)=nihedata(i,5)^2*nihedata(i,1)^3*nihedata(i,6)^3;
    a1(i,57)=nihedata(i,5)^2*nihedata(i,1)^3*nihedata(i,6)^4;
    
    a1(i,58)=nihedata(i,5)*nihedata(i,4)*nihedata(i,1);
    a1(i,59)=nihedata(i,5)*nihedata(i,4)*nihedata(i,1)^2;
    a1(i,60)=nihedata(i,5)*nihedata(i,4)*nihedata(i,1)^3;
    a1(i,61)=nihedata(i,5)*nihedata(i,4)*nihedata(i,6);
    a1(i,62)=nihedata(i,5)*nihedata(i,4)*nihedata(i,6)^2;
    a1(i,63)=nihedata(i,5)*nihedata(i,4)*nihedata(i,6)^3;
    a1(i,64)=nihedata(i,5)*nihedata(i,4)*nihedata(i,6)^4;
    a1(i,65)=nihedata(i,5)*nihedata(i,4)*nihedata(i,1)*nihedata(i,6);
    a1(i,66)=nihedata(i,5)*nihedata(i,4)*nihedata(i,1)*nihedata(i,6)^2;
    a1(i,67)=nihedata(i,5)*nihedata(i,4)*nihedata(i,1)*nihedata(i,6)^3;
    a1(i,68)=nihedata(i,5)*nihedata(i,4)*nihedata(i,1)*nihedata(i,6)^4;
    a1(i,69)=nihedata(i,5)*nihedata(i,4)*nihedata(i,1)^2*nihedata(i,6);
    a1(i,70)=nihedata(i,5)*nihedata(i,4)*nihedata(i,1)^2*nihedata(i,6)^2;
    a1(i,71)=nihedata(i,5)*nihedata(i,4)*nihedata(i,1)^2*nihedata(i,6)^3;
    a1(i,72)=nihedata(i,5)*nihedata(i,4)*nihedata(i,1)^2*nihedata(i,6)^4;
    a1(i,73)=nihedata(i,5)*nihedata(i,4)*nihedata(i,1)^3*nihedata(i,6);
    a1(i,74)=nihedata(i,5)*nihedata(i,4)*nihedata(i,1)^3*nihedata(i,6)^2;
    a1(i,75)=nihedata(i,5)*nihedata(i,4)*nihedata(i,1)^3*nihedata(i,6)^3;
    a1(i,76)=nihedata(i,5)*nihedata(i,4)*nihedata(i,1)^3*nihedata(i,6)^4;
    
    a1(i,77)=nihedata(i,5)*nihedata(i,4)^2*nihedata(i,1);
    a1(i,78)=nihedata(i,5)*nihedata(i,4)^2*nihedata(i,1)^2;
    a1(i,79)=nihedata(i,5)*nihedata(i,4)^2*nihedata(i,1)^3;
    a1(i,80)=nihedata(i,5)*nihedata(i,4)^2*nihedata(i,6);
    a1(i,81)=nihedata(i,5)*nihedata(i,4)^2*nihedata(i,6)^2;
    a1(i,82)=nihedata(i,5)*nihedata(i,4)^2*nihedata(i,6)^3;
    a1(i,83)=nihedata(i,5)*nihedata(i,4)^2*nihedata(i,6)^4;
    a1(i,84)=nihedata(i,5)*nihedata(i,4)^2*nihedata(i,1)*nihedata(i,6);
    a1(i,85)=nihedata(i,5)*nihedata(i,4)^2*nihedata(i,1)*nihedata(i,6)^2;
    a1(i,86)=nihedata(i,5)*nihedata(i,4)^2*nihedata(i,1)*nihedata(i,6)^3;
    a1(i,87)=nihedata(i,5)*nihedata(i,4)^2*nihedata(i,1)*nihedata(i,6)^4;
    a1(i,88)=nihedata(i,5)*nihedata(i,4)^2*nihedata(i,1)^2*nihedata(i,6);
    a1(i,89)=nihedata(i,5)*nihedata(i,4)^2*nihedata(i,1)^2*nihedata(i,6)^2;
    a1(i,90)=nihedata(i,5)*nihedata(i,4)^2*nihedata(i,1)^2*nihedata(i,6)^3;
    a1(i,91)=nihedata(i,5)*nihedata(i,4)^2*nihedata(i,1)^2*nihedata(i,6)^4;
    a1(i,92)=nihedata(i,5)*nihedata(i,4)^2*nihedata(i,1)^3*nihedata(i,6);
    a1(i,93)=nihedata(i,5)*nihedata(i,4)^2*nihedata(i,1)^3*nihedata(i,6)^2;
    a1(i,94)=nihedata(i,5)*nihedata(i,4)^2*nihedata(i,1)^3*nihedata(i,6)^3;
    a1(i,95)=nihedata(i,5)*nihedata(i,4)^2*nihedata(i,1)^3*nihedata(i,6)^4;
    
    a1(i,96)=nihedata(i,5)^2*nihedata(i,4)*nihedata(i,1);
    a1(i,97)=nihedata(i,5)^2*nihedata(i,4)*nihedata(i,1)^2;
    a1(i,98)=nihedata(i,5)^2*nihedata(i,4)*nihedata(i,1)^3;
    a1(i,99)=nihedata(i,5)^2*nihedata(i,4)*nihedata(i,6);
    a1(i,100)=nihedata(i,5)^2*nihedata(i,4)*nihedata(i,6)^2;
    a1(i,101)=nihedata(i,5)^2*nihedata(i,4)*nihedata(i,6)^3;
    a1(i,102)=nihedata(i,5)^2*nihedata(i,4)*nihedata(i,6)^4;
    a1(i,103)=nihedata(i,5)^2*nihedata(i,4)*nihedata(i,1)*nihedata(i,6);
    a1(i,104)=nihedata(i,5)^2*nihedata(i,4)*nihedata(i,1)*nihedata(i,6)^2;
    a1(i,105)=nihedata(i,5)^2*nihedata(i,4)*nihedata(i,1)*nihedata(i,6)^3;
    a1(i,106)=nihedata(i,5)^2*nihedata(i,4)*nihedata(i,1)*nihedata(i,6)^4;
    a1(i,107)=nihedata(i,5)^2*nihedata(i,4)*nihedata(i,1)^2*nihedata(i,6);
    a1(i,108)=nihedata(i,5)^2*nihedata(i,4)*nihedata(i,1)^2*nihedata(i,6)^2;
    a1(i,109)=nihedata(i,5)^2*nihedata(i,4)*nihedata(i,1)^2*nihedata(i,6)^3;
    a1(i,110)=nihedata(i,5)^2*nihedata(i,4)*nihedata(i,1)^2*nihedata(i,6)^4;
    a1(i,111)=nihedata(i,5)^2*nihedata(i,4)*nihedata(i,1)^3*nihedata(i,6);
    a1(i,112)=nihedata(i,5)^2*nihedata(i,4)*nihedata(i,1)^3*nihedata(i,6)^2;
    a1(i,113)=nihedata(i,5)^2*nihedata(i,4)*nihedata(i,1)^3*nihedata(i,6)^3;
    a1(i,114)=nihedata(i,5)^2*nihedata(i,4)*nihedata(i,1)^3*nihedata(i,6)^4;
    
    a1(i,115)=nihedata(i,5)^2*nihedata(i,4)^2*nihedata(i,1);
    a1(i,116)=nihedata(i,5)^2*nihedata(i,4)^2*nihedata(i,1)^2;
    a1(i,117)=nihedata(i,5)^2*nihedata(i,4)^2*nihedata(i,1)^3;
    a1(i,118)=nihedata(i,5)^2*nihedata(i,4)^2*nihedata(i,6);
    a1(i,119)=nihedata(i,5)^2*nihedata(i,4)^2*nihedata(i,6)^2;
    a1(i,120)=nihedata(i,5)^2*nihedata(i,4)^2*nihedata(i,6)^3;
    a1(i,121)=nihedata(i,5)^2*nihedata(i,4)^2*nihedata(i,6)^4;
    a1(i,122)=nihedata(i,5)^2*nihedata(i,4)^2*nihedata(i,1)*nihedata(i,6);
    a1(i,123)=nihedata(i,5)^2*nihedata(i,4)^2*nihedata(i,1)*nihedata(i,6)^2;
    a1(i,124)=nihedata(i,5)^2*nihedata(i,4)^2*nihedata(i,1)*nihedata(i,6)^3;
    a1(i,125)=nihedata(i,5)^2*nihedata(i,4)^2*nihedata(i,1)*nihedata(i,6)^4;
    a1(i,126)=nihedata(i,5)^2*nihedata(i,4)^2*nihedata(i,1)^2*nihedata(i,6);
    a1(i,127)=nihedata(i,5)^2*nihedata(i,4)^2*nihedata(i,1)^2*nihedata(i,6)^2;
    a1(i,128)=nihedata(i,5)^2*nihedata(i,4)^2*nihedata(i,1)^2*nihedata(i,6)^3;
    a1(i,129)=nihedata(i,5)^2*nihedata(i,4)^2*nihedata(i,1)^2*nihedata(i,6)^4;
    a1(i,130)=nihedata(i,5)^2*nihedata(i,4)^2*nihedata(i,1)^3*nihedata(i,6);
    a1(i,131)=nihedata(i,5)^2*nihedata(i,4)^2*nihedata(i,1)^3*nihedata(i,6)^2;
    a1(i,132)=nihedata(i,5)^2*nihedata(i,4)^2*nihedata(i,1)^3*nihedata(i,6)^3;
    a1(i,133)=nihedata(i,5)^2*nihedata(i,4)^2*nihedata(i,1)^3*nihedata(i,6)^4;
    
    
    
    
   %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 
    a1(i,134)=nihedata(i,3);
    a1(i,135)=nihedata(i,2);
    a1(i,136)=nihedata(i,3)*nihedata(i,2);
    a1(i,137)=nihedata(i,3)*nihedata(i,2)^2;
    a1(i,138)=nihedata(i,3)^2*nihedata(i,2);
    a1(i,139)=nihedata(i,3)^2*nihedata(i,2)^2;    
    
    
    a1(i,140)=nihedata(i,6)*nihedata(i,3);
    a1(i,141)=nihedata(i,6)*nihedata(i,2);
    a1(i,142)=nihedata(i,6)*nihedata(i,3)*nihedata(i,2);
    a1(i,143)=nihedata(i,6)*nihedata(i,3)*nihedata(i,2)^2;
    a1(i,144)=nihedata(i,6)*nihedata(i,3)^2*nihedata(i,2);
    a1(i,145)=nihedata(i,6)*nihedata(i,3)^2*nihedata(i,2)^2;    
    
    
    a1(i,146)=nihedata(i,6)^2*nihedata(i,3);
    a1(i,147)=nihedata(i,6)^2*nihedata(i,2);
    a1(i,148)=nihedata(i,6)^2*nihedata(i,3)*nihedata(i,2);
    a1(i,149)=nihedata(i,6)^2*nihedata(i,3)*nihedata(i,2)^2;
    a1(i,150)=nihedata(i,6)^2*nihedata(i,3)^2*nihedata(i,2);
    a1(i,151)=nihedata(i,6)^2*nihedata(i,3)^2*nihedata(i,2)^2;    
   %%%%%%4%%%%%%%%%%%%%%%%%%3%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %%%
   %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %%%
   
   %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %%%


   
end
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%  
a2=zeros(5120,151);
a3=zeros(5120,151);
a4=zeros(5120,151);
a5=zeros(5120,151);
a6=zeros(5120,151);
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%  
load d:\matlab\work\model\data.txt;
for i=1:5120
    a2(i,:)=a1(i,:)*data(i,5);
    a3(i,:)=a1(i,:)*data(i,4);
    a4(i,:)=a1(i,:)*data(i,3);
    a5(i,:)=a1(i,:)*data(i,2);
    a6(i,:)=a1(i,:)*data(i,1);
end

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%  
xxx=[a1,a2,a3,a4,a5,a6];

load d:\matlab\work\model\y1.txt;


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


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 
k1=kkk(1:151,1)';
k2=kkk(152:302,1)';
k3=kkk(303:453,1)';
k4=kkk(454:604,1)';
k5=kkk(605:755,1)';
k6=kkk(756:906,1)';

yy1=zeros(5120,1);
y0=zeros(1,151);
for j=1:5120

y0=k1.*a1(j,:)+k2.*a2(j,:)+k3.*a3(j,:)+k4.*a4(j,:)+k5.*a5(j,:)+k6.*a6(j,:);
t=0;
for i=1:151
    t=t+y0(1,i);
    yy1(j,1)=t;
end
y0=zeros(1,151);
end
save D:\matlab\work\model\K_result1\Kc.txt k1 -ASCII;
save D:\matlab\work\model\K_result1\Krf.txt k2 -ASCII;
save D:\matlab\work\model\K_result1\Kbfw.txt k3 -ASCII;
save D:\matlab\work\model\K_result1\Kbfi.txt k4 -ASCII;
save D:\matlab\work\model\K_result1\Kcw2.txt k5 -ASCII;
save D:\matlab\work\model\K_result1\Kcw4.txt k6 -ASCII;