for loop which changes color depending on the month

I am trying to write a for loop where each month on my graph is a different colour. I am very unfamiliar with for loop and im unsure on how i work it into the code. I want the first subplot to have a rainbow essentially. Here is what I currently have.

c=hsv(12);

plotIndex = month(windTurbineData.disc_TimeStamp);

meancapacity = windTurbineData.mean_Power_kW/7000

colormap(c)

for i = 1:12

subplot (2,2,[1,2]);

plot(windTurbineData.disc_TimeStamp,meancapacity)

title('Power Output 2019')

ylabel('Power output/rated power [-]')

ylim([-0.2 1.2])

grid on

end

subplot (2,2,3);

title('Power Output vs Wind Speed 2019')

xlabel('Wind Speed m/s')

ylabel('Power output/rated power [-]')

colorbar

subplot(2,2,4);

title('Standard Deviation in Power Output 2019')

xlabel('Wind Speed m/s')

ylabel('Std Power output/rated power [-]')

colorbar('Ticks', linspace(0, 1, 12), 'TickLabels',{'Jan','Feb','March','April','May','June','July','August','Sept','Oct','Nov','Dec'})

Answer 1

[DGM]

I don't have your data, so I just have to assume some placeholder data instead.

It's probably better to use a legend for identifying line plots like this. If you want to actually use a colorbar, you'll need to correctly configure it so that the labels line up with the discrete color segments, otherwise it's ambiguous.

Depending on what your actual data looks like, overlaying that many line plots might end up being pretty hard to read anyway.

monthnames = {'Jan','Feb','March','April','May','June','July','August','Sept','Oct','Nov','Dec'};

% placeholder data
t = datenum([2000 2000],[01 12],[01 31]);
t = t(1):t(2);
P = rand(1,numel(t));
S = P + 0.05*randn(1,numel(t));

% normalized power per month
plotIndex = month(t);
meancapacity = P;

subplot (2,2,[1,2]);
hold on
for moy = 1:12
    mask = plotIndex==moy;
    plot(day(t(mask)),meancapacity(mask))
end
title('Power Output 2019')
ylabel('Power output/rated power [-]')
xlabel('Day of Month')
ylim([-0.2 1.2])
grid on
set(gca,'colororder',hsv(12));
legend(get(gca,'children'),monthnames,'location','eastoutside');


% power vs speed
[Ssorted idx] = sort(S);
Psorted = P(idx);

subplot (2,2,3);
plot(Ssorted,Psorted);
title('Power Output vs Wind Speed 2019')
xlabel('Wind Speed m/s')
ylabel('Power output/rated power [-]')


% std per month
Pstd = zeros(1,12);
for moy = 1:12
    mask = plotIndex==moy;
    Pstd(moy) = std(P(mask));
end

subplot(2,2,4);
plot(1:12,Pstd)
title('Standard Deviation in Power Output 2019')
ylabel('Std Power output/rated power [-]')
set(gca,'xtick',1:12,'xticklabels',monthnames,'xticklabelrotation',90)

The titles and labels all collide here on the web display, but should be better if the figure is allowed to be larger.

EDIT:

https://www.mathworks.com/matlabcentral/answers/1600915-creating-a-plot-that-changes-colour-as-the-month-changes

Related question that should be closer to the goals.