Try something like this —
% clc,clearvars,
% close all
% Example X and Y data
x = linspace(0, 50, 100); % Common X-axis for all datasets
y1 = exp(-0.1*(x-10).^2);
y2 = exp(-0.1*(x-20).^2);
y3 = exp(-0.1*(x-30).^2);
% % Plotting
figure;
hold on
patch([x flip(x)], [zeros(1,size(x,2)*2)], [zeros(size(x)) flip(y1)], 'r', EdgeColor='r', DisplayName='y_1')
patch([x flip(x)], [zeros(1,size(x,2)*2)]+0.3, [zeros(size(x)) flip(y2)], 'g', EdgeColor='g', DisplayName='y_2')
patch([x flip(x)], [zeros(1,size(x,2)*2)]+0.6, [zeros(size(x)) flip(y3)], 'b', EdgeColor='b', DisplayName='y_3')
hold off
view([-45,25])
grid on
zlim([-0.5 2.5])
xlabel('Particle size');
ylabel('Particle density');
zlabel('De');
legend(Location='best')
Here, the patch funciton will produce 3D plots if given three arguments. The first of these is the common ‘x’ row vector, the second is an spppropriate-length zeros vector with an added offset for the ‘y’ coordinate, and the ‘y#’ vectors are the ‘z’ coordinates. The patch funciton fills a closed curve, and the closures are created by using flip function to duplicate the original vectors with their orders reversed.
Make appropriate changes to get the result you want.
.
