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