How do I add a colormap to a series of 3D lines?

Hello,

For a 3D plot I have:

U = 2:0.2:30;

f is a 1x20002 vector of frequencies

the Z axis is a 141x20002 matrix of values between 0 and 1.

An example figure of what I am trying to achieve is:

As you can see from the figure, there is a spectrum line for each U* and most of the frequencies are between 0-2.

Now, I can plot it using plot3 for my data just fine, please see the figure below:

However, I think the parula colormap for each U* line is quite nice. I have checked that it is not possible to do this with the plot3 function.

I have tried the surf function as well but there is something missing.

Would anyone please be able to instruct me on how could I do it?

Answer 1

[KSSV]

Read about waterfall.

[Bruno Soares]

Hi KSSV,

Thanks for replying.

Is it really possible with waterfall? From what I read it generates the coloured line as I wish but it does not fill the surface of the peaks as in the figure sent??

[KSSV]

Else try this:

x = 0:.05:2*pi;
y = sin(x);
z = zeros(size(x));  % give your z here 
col = x;  % This is the color, vary with x in this case.

surface([x;x],[y;y],[z;z],[col;col],...
        'facecol','no',...
        'edgecol','interp',...
        'linew',2);

[Bruno Soares]

Thanks again KSSV.

Quite close, but still something missing... The base now merged and became part of the surface, while it was supposed to be separate U lines. Any other advice please?

[KSSV]

It looks you have defind fine intervel in x. Reduce the number of intervels and see.

[Bruno Soares]

I divided the step by 2 and got the following:

[KSSV]

Try scatter with marker style as squares. Share your data, let me try as well. Show your code which you have tried. That should be possible with waterfall too.

[Bruno Soares]

Sure, please see the code below and data attached:

Fs=100;

spacing=1;

cut=20000; Tspan=0:0.01:600; Ured1=2:4:28;

L=length(Tspan(cut:end));

f = 2*pi*Fs*(0:(L/2))/L;

for ijj=1:length(Ured1)

p2=Y1(ijj,:);

p2=p2(:);

Y222 = fft(detrend(p2(cut:end),'constant'));

P21 = abs(Y222/L);

P11 = P21(1:L/2+1);

P11(2:end-1) = 2*P11(2:end-1);

ydft = Y222(1:L/2+1);

psd2 = (1/(Fs*L)) * abs(ydft).^2;

normalised(:,ijj) = (psd2-min(psd2))/(max(psd2)-min(psd2));

end

figure(16);clf;

set(gcf,'color','w');

MatrixUred1=ones(length(Ured1),length(f));

Matrixf=ones(length(Ured1),length(f));

for i=1:length(Ured1)

MatrixUred1(i,:)=MatrixUred1(i,:)*Ured1(i);

end

Matrixf=repmat(f,7,1);

% surf(MatrixUred1',Matrixf',normalised,...

% 'facecol','no',...

% 'edgecol','interp',...

% 'linew',2);

plot3(MatrixUred1',Matrixf',normalised,'k')

hold on

[M,I]=max(normalised);

a=I;

ylim([0 5]);xlim([2 30]);

xticks(2:4:30);yticks(0:5);zticks(0:0.1:1);