最小二乘法电池参数辨识

clear all;
a=xlsread（‘C:\Users\MSI-PC\Desktop\新建文件夹\递推最小二乘法在线辨识锂离子电池参数/fe081dst‘）;
%a=xlsread（‘E:/test/15tem‘）;此处电流是DST循环工况
y=a（:2）;%votage
u=a（:3）;%current
u=u;
s=length（y）;
%y=y（5000:48000）;
%u=u（5000:48000）;
%s=length（y）;
%soc=a（:4）;
soc（1）=0.999;
for i=1:s-1
    soc（i+1）=soc（i）+u（i）/（3600*8.1）;
    if soc（i+1）<0.01
        soc（i+1）=0.01;
    end
end
%y1=y（13385:13422）;
%u1=u（13385:13422）;
%s=length（y1）;
b（:1）=[1 log（0.999） log（0.001） 0 u（1） 0]‘;%state
x=1;%forgetting factor
%chushihua
c（:1）=[0.1 0.1 0.1 0.9 0 0]‘;%estimated needed
P=10000*eye（66）;
for i=2:s
    b（:i）=[1 log（so