This isn't exactly the same as your plot, but I hope it helps. If you want a patterned fill on your bars, you can check out Hatchfill2 (see some examples here). You may have to convert it to a histogram or play around with the bars if you don't want space between the bars.
% Generate data
y = zeros(5,4);
for i = 1:5
y(i,:) = (6-i)*10 + randi(10,[1,4]);
end
% Plot bar graph
figure(1);clf
b = bar(1:5,y,'FaceAlpha',0.55);
b(1).FaceColor = [0.9290 0.6940 0.1250];
b(2).FaceColor = [0.4660 0.6740 0.1880];
b(3).FaceColor = [0.4940 0.1840 0.5560];
b(4).FaceColor = [1 1 0];
% Axis ticks
xticklabels({'R','200','400','600','800'})
a = gca;
a.TickLabelInterpreter = 'latex';
a.FontSize = 12;
% Labels and legend
xlabel('Temperature in $^\circ$C','Interpreter','latex','FontSize',13)
ylabel('Compressive Strength (MPa)','Interpreter','latex','FontSize',13)
legend({'GPC (Control)','10\% GP-A','20\% GP-A','30\% GP-A'}, ...
'Interpreter','latex','FontSize',12)
% Errorbars
hold on
sd = randi(3,5,4);
errorbar([b.XEndPoints],[b.YEndPoints],sd(:), ...
'Color','k','LineStyle','none','HandleVisibility','off')
hold off
