How to calculate 8 days mean from ten years data?

Kshitij Singh answered . 2024-12-21 16:59:40

%% read data file
[lines, ~] = readlog('Test2.csv'); % fex : https://fr.mathworks.com/matlabcentral/fileexchange/96364-readlog?s_tid=srchtitle
data_matrix = split(lines,',');
[m,n] = size(data_matrix);
ind_start_year = find(startsWith(data_matrix(:,1),'1-Jan-')); % find row corresponding to first day of the year
ind_stop_year = find(startsWith(data_matrix(:,1),'31-Dec-')); % find row corresponding to last day of the year
toc

%% main loop
windowSize = 8;  % define buffer length for averaging
cols = 3:n; % define cols to be processed
dataAverages_total = [];

% if the data file stops before a 31 Dec , we must add this : 
if numel(ind_stop_year) < numel(ind_start_year)
    ind_stop_year = [ind_stop_year; m];
end

for ck = 1:numel(ind_start_year)% loop over year number
    ind_start_thisyear = ind_start_year(ck);
    ind_stop_thisyear = ind_stop_year(ck);
    start = ind_start_thisyear:windowSize:ind_stop_thisyear; 
    stop = [start(2:end)-1, ind_stop_thisyear]; 
    dur = stop-start+1;
    nbins = numel(start); 
    dataAverages = nan(nbins, numel(cols)); 
    for ci = 1:nbins
        rows = start(ci) : stop(ci);
        tmp = data_matrix(rows, :);
        dataAverages(ci,:) = mean(str2double(tmp(:,3:n)), 1, 'omitnan'); % results of one year
    end
    sum(dur)
    dataAverages_total = [dataAverages_total ; dataAverages]; % concatenate every year's result (dataAverages) 
end

Not satisfied with the answer ?? ASK NOW