递归最小二乘法参数辨识

% 最小二乘的递推算法 
% U（k）-U_ocv（k）=a1*I（k）+a2*（U（k-1）-U_ocv（k-1））+a3*I（k-1）;
%clear  
tic

%----------------------------读入数据--------------------------------
data= xlsread（‘1110Fe1udds25.xls‘）;%1电流 2电压 3soc 4soc带格式
for i=1:5200 
   U_ocv（i）=Uocv（ data（i+304） ） ;%为了避免soc=1从第三十个开始读取5200个数
   z（i）=data（i+302）-U_ocv（i）;
   M（i）=-data（i+301）;
end
%----------------------------递推求解--------------------------------
P=10*eye（3）;   % 估计方差
Pstore=zeros（35001）; 
Pstore（:1）=[P（11）;P（22）;P（33）]; % P矩阵
Theta=zeros（35002）; %参数的估计值，存放中间过程估值 
Theta（:1）=[1;1;1]; 
% K=zeros