%?最小二乘的遞推算法?
%?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