Try this —
x = 210:0.5:227;
y = randn(size(x));
bands = [214 216; 221 223]; % X-Coordinates Of Band Limits: [x(1,1) x(1,2); x(2,1) x(2,2)]
figure
plot(x, y)
hold on
xp = [bands fliplr(bands)]; % X-Coordinate Band Definitions
yp = ([[1;1]*min(ylim); [1;1]*max(ylim)]*ones(1,size(bands,1))).'; % Y-Coordinate Band Definitions
for k = 1:size(bands,1) % Plot Bands
patch(xp(k,:), yp(k,:), [1 1 1]*0.25, 'FaceAlpha',0.5, 'EdgeColor',[1 1 1]*0.25)
end
hold off
This is fully adaptive to the size of the plot limits, requiring only that the ‘bands’ matrix be defined.
I doubt if there is a way to do this without the loop.
.
