Calculation of diurnal cycle

Hi everyone!

I have some .csv files (attached you can see one of them) containing hourly values of temperature and humidity from 1/7/2019 to 30/9/2019. Frome these values I want to extract the diurnal cycle of temperature and humidity, which means 1 mean value for all dates at 00:00, one for 01:00, till 23:00.

What is the most efficient way to do so?

Thank you in advance!

[darova]

What about mean function? Did you try it?

[Daphne PARLIARI]

Sure but I want to take the mean considering the hours, that's what I am not sure how to achieve...

Answer 1

[Star Strider]

Try this:

D = readtable('Airport.csv');
Hrs = hour(D.Date);                                                                 % Hours
[UHr,~,ic] = unique(Hrs);                                                           % Hours Reference Index Vector
TRH = accumarray(ic, (1:numel(Hrs)).', [], @(x){mean([D.Temp(x),D.RelHum(x)])});    % Means By Hour

TRHmtx = cell2mat(TRH);                                                             % Matrix From Cell Array: [Temp RelHum]

figure
yyaxis left
plot(UHr, TRHmtx(:,1))
ylabel('Temperature (°C)')
yyaxis right
plot(UHr, TRHmtx(:,2))
ylabel('Relative Humidity (%)')
grid
xlabel('Time (Hours)')
set(gca, 'Xtick',(0:23))
xlim([0 23])

producing:

[Daphne PARLIARI]

Thank you!!!

[Star Strider]

As always, my pleasure!