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

%----------------------------讀入數(shù)據(jù)--------------------------------
data=?xlsread（‘1110Fe1udds25.xls‘）;%1電流?2電壓?3soc?4soc帶格式
for?i=1:5200?
???U_ocv（i）=Uocv（?data（i+304）?）?;%為了避免soc=1從第三十個(gè)開(kāi)始讀取5200個(gè)數(shù)
???z（i）=data（i+302）-U_ocv（i）;
???M（i）=-data（i+301）;
end
%----------------------------遞推求解--------------------------------
P=10*eye（3）;???%?估計(jì)方差
Pstore=zeros（35001）;?
Pstore（:1）=[P（11）;P（22）;P（33）];?%?P矩陣
Theta=zeros（35002）;?%參數(shù)的估計(jì)值，存放中間過(guò)程估值?
Theta（:1）=[1;1;1];?
%?K=zeros