Longcycle Code

clear

current_density=input('Please input the current density(A/g):');

filename_all=dir('*.txt'); %得到当前文件夹的所有TXT文件信息

dircell=struct2cell(filename_all); %对读取信息的格式转换,struct转cell,注意cell的数据类型为char

filename=dircell(1,:); %提取文件名

sortfilename=sortnat(dircell(1,:));

file_size=size(sortfilename);file_num=file_size(1,2);%得到文件的个数

for n=1:file_num

filename_n=char(sortfilename(1,n));%将字符型数据转化为字符串string

A=csvread(filename_n,20,0);%读取TXT文件

sizeA=size(A);row_num=sizeA(1,1);

for l=2:(row_num-1)

if A(l,2)>A(l+1,2)&&A(l,2)>A(l-1,2)

t_maxV=A(l,1);V_max=A(l,2);

else

t_minV=A(row_num,1);V_min=A(row_num,2);

end

end

specific_capacitance(n,1)=current_density*(t_minV-t_maxV)/(V_max-V_min);

end

specific_capacitance