Slice of a surface plot created by 2D data

33 visualizaciones (últimos 30 días)
Mai Sakuragi
Mai Sakuragi el 9 de Feb. de 2020
Respondida: Rajani Mishra el 13 de Mzo. de 2020
Hi, I have been trying to create a slice though my surface plot which is created with stacked 2D data.
Fs = 3000;
tspan = 0 : 1/Fs :15;
topt = 3
j = topt-1.0 : 0.05 : topt+4.0
---------------------------------------------------------------
delta_r1 = zeros(length(tspan), length(j));
delta_r2 = zeros(length(tspan), length(j));
for ii = 1:numel(j)
t_THz = j(ii);
delta_r1(:,ii) = Fourier1P - fourier1P(:,ii);
delta_r2(:,ii) = Fourier2P - fourier2P(:,ii);
end
delay = zeros(1,length(j));
for ii = 1:numel(j)
t_THz = j(ii);
delay(ii) = t_THz - topt;
end
yMat = repmat(delay, numel(tspan),1);
xMat = repmat(F1.',1,length(j));
plot3(xMat,yMat,delta_r1)
xlim([358 362.5])
zlim([-0.4e-4 2e-4])
xlabel('Frequency (THz)','Rotation', -15)
ylabel('Delay (ps)','Rotation', 23)
title('x_1(f)')
view([19.245913 34.150412])
figure
surf(xMat,yMat,delta_r1)
xlim([358.5 361.7])
zlim([-0.4e-4 2e-4])
xlabel('Frequency (THz)','Rotation',-15)
ylabel('Delay (ps)','Rotation',23)
view([40.018393 23.670967])
colormap(jet)
shading interp
where
delta_r1
and
delta_r2
are 2D arrays of size 45001×101.
I want to create a slice though x = 360.9 so that I can get the side view of the oscillation.
What I have now is this:
x = linspace(358,362.5,length(tspan))
delay
[X,Y] = meshgrid(x,delay)
yq = linspace(-1,4,101)
xq = linspace(360.9,360.9,101)
zq = interp2(X,Y,delta_r1',xq,yq);
dq = sqrt((xq-360.9).^2 + (yq+1).^2);
plot(dq,zq)
axis([min(dq),max(dq),-0.4e-4,2e-4])
I get this plot
It's a straight line rather than an oscillation.
Does anyone know how to fix it?
  1 comentario
Thiago Henrique Gomes Lobato
Thiago Henrique Gomes Lobato el 9 de Feb. de 2020
The problem may be in your interpolation function since you're attributing different position values to your calculated results. You may try something like this:
zq = interp2(xMat,yMat,delta_r1,xq,yq);
dq = sqrt((xq-360.9).^2 + (yq+1).^2);
plot(dq,zq)

Iniciar sesión para comentar.

Respuestas (1)

Rajani Mishra
Rajani Mishra el 13 de Mzo. de 2020
You can try using function “slice()”, refer to this link for more information : https://www.mathworks.com/help/matlab/ref/slice.html#d118e1325646. Also I found a similar question here : https://www.mathworks.com/matlabcentral/answers/488992-how-to-take-a-horizontal-slice-of-a-3d-surface-plot

Categorías

Más información sobre Interpolation en Help Center y File Exchange.

Productos


Versión

R2018b

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by