Plot legend labels using loop

I am plotting the power spectrum from 16 channels. PSD was calculated using the pwelch function and its output has 3 dimensions: (Channels) x (Frequency index) x (Trial). Is there a way to produce the legend (data 1, data 2,...) with the channel numbers using a for loop, in order to avoid writing down all 16 labels?

thank you!

Here is what i get by the way:

figure(1), hold on
plot(f, mean(pwrsignal(:,:,:), 3), 'LineWidth', 2)
xlim([0 200])
legend

Answer 1

[Cris LaPierre]

data = rand(100,16);
plot(data)
legend("data " + (1:size(data,2)))

[Panos Kerezoudis]

great, thank you so much!