%?example?script?that?recreates?the?plots?shown?in?the?README?document
close?all;
clear?all;

%?Change?this?data?file?name?below?to?plot?different?files
load（‘data/Matlab/RW25‘）

%%?First?the?constant?load?profiles?that?are?run?after?every?50?random?walk?
%?discharge?cycles?are?plotted?using

steps?=?data.step;?%?save?steps?array?to?new?variable

%?initialize?plot
figure?hold?on
xlim（[-.12.5]）
ylim（[34.25]）
xlabel（‘Time?（h）‘）;
ylabel（‘Voltage?（V）‘）;
title（‘Reference?Discharge?Profiles‘）

for?i?=?1:length（steps）?%?search?through?the?array?of?step?structures
????if?strcmp（steps（i）.comment‘reference?discharge‘）?
????????RT?=?steps（i）.relativeTime/3600;?%?relative?time?in?hours
????????V?=?steps（i）.voltage;?%?voltage
????????plot（RTV‘k‘）
????end
end


%%?We?can?benchmark?the?battery抯?capacity?by?integrating?current?over?the?
%?reference?cycles.?The?next?plot?shows?this?capacity?measurement?vs?date.

counter?=?0;?%?initialize?counters

for?i?=?1:length（steps）?%?search?through?the?array?of?step?structures
????if?strcmp（steps（i）.comment‘reference?discharge‘）?
????????counter?=?counter+1;

????????%?save?the?date?of?this?reference?discharge?step
????????date（counter）?=?datenum（steps（i）.date）;
????????
????????%?save?the?battery?capacity?measured?by?integrating?current?
????????capacity（counter）?=?trapz（steps（i）.relativeTime/3600data.step（i）.current）;
????end
end

%?plot?the?data
figure
plot（datecapacity‘o‘）
datetick（‘x‘‘mmmyyyy‘）
xlabel（‘Date‘）;
ylabel（‘Capacity?（Ah）‘）;
title（‘Degradation?of?Measured?Capacity‘）

%%?The?next?code?snippet?shows?how?the?resting?periods?after?each?reference?
%?discharge?can?be?grouped?with?the?corresponding?reference?discharge?cycle?
%?and?added?to?the?plot?above.?

%?initialize?Reference?Discharge?Struct?that?will?hold?indexes?of?reference
%?discharge?steps?and?the?rest?steps?that?come?imediately?before?and?after
RefDis.repition{1}.indexes?=?0;

%?initialize?loop?parameters
counter?=?0;
for?i?=?1:length（data.step）
????if?strcmp（steps（i）.comment‘reference?discharge‘）
????????counter?=?counter+1;
????????if?strcmp（steps（i-1）.comment‘rest?prior?reference?discharge‘）
????????????RefDis.repition{counter}.indexes?=?[i-1?i];
????????else
????????????RefDis.repition{counter}.indexes?=?i;
????????end
????end
????
????if?strcmp（steps（i）.comment‘rest?post?reference?discharge‘）
????????RefDis.repition{counter}.indexes?=?[RefDis.repition{counter}.indexes?i];
????end
end

%?now?plot?each?of?the?grouped?reference?discharge?cycles
figure?hold?on
for?i?=?1:length（RefDis.repition）
????%?stitch?together?all?of?the?substeps?identified?in?the?repitition
????%?field
????step?=?steps（RefDis.repition{i}.indexes（1））;
????
????%?if?the?first?index?is?a?rest?period?then?make?it?end?at?relativeTime?==?0?
????if（strcmp（step.comment‘rest?prior?reference?discharge‘））
????????step.relativeTime?=?step.relativeTime?-?step.relativeTime（end）;?
????end

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????目錄???????????0??2015-03-24?18:58??RW_Skewed_High_40C_DataSet_2Post\data\
?????目錄???????????0??2015-03-24?18:58??RW_Skewed_High_40C_DataSet_2Post\data\Matlab\
?????文件????11188968??2015-03-12?17:19??RW_Skewed_High_40C_DataSet_2Post\data\Matlab\RW25.mat
?????文件?????9750077??2015-03-12?18:03??RW_Skewed_High_40C_DataSet_2Post\data\Matlab\RW26.mat
?????文件????12156834??2015-03-12?19:15??RW_Skewed_High_40C_DataSet_2Post\data\Matlab\RW27.mat
?????文件????11154359??2015-03-12?19:57??RW_Skewed_High_40C_DataSet_2Post\data\Matlab\RW28.mat
?????目錄???????????0??2015-03-24?18:58??RW_Skewed_High_40C_DataSet_2Post\data\R\
?????文件?????4583572??2015-03-25?03:47??RW_Skewed_High_40C_DataSet_2Post\data\R\RW25.Rda
?????文件?????3951653??2015-03-25?03:50??RW_Skewed_High_40C_DataSet_2Post\data\R\RW26.Rda
?????文件?????4978283??2015-03-25?03:55??RW_Skewed_High_40C_DataSet_2Post\data\R\RW27.Rda
?????文件?????4607175??2015-03-25?10:08??RW_Skewed_High_40C_DataSet_2Post\data\R\RW28.Rda
?????文件???????12222??2015-03-24?19:01??RW_Skewed_High_40C_DataSet_2Post\MatlabSamplePlots.m
?????文件??????126077??2015-03-26?05:46??RW_Skewed_High_40C_DataSet_2Post\README_RW_SkewedHigh_40C.html
?????文件???????20684??2015-03-26?05:46??RW_Skewed_High_40C_DataSet_2Post\README_RW_SkewedHigh_40C.Rmd