How to create stacked quiver plots?

Question

doruk isik 2022-9-1

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1793140-how-to-create-stacked-quiver-plots

评论： Bjorn Gustavsson 2022-9-2

采纳的回答： Bjorn Gustavsson

Hello,

I have a dataset with 6 columns corresponding to slices taken at different x-locations (see below)

The dataset is created as (x_coordinate,y_coordinate,z_coordinate,Vel_x, Vel_y, Vel_z) where Vel_x, Vel_y and Vel_z are x,y and z velocities, respectively.

Based on the x-coordinates, I created individual arrays at each x-location named slice1, slice2, slice3 etc. where each array has a similar structure with the first column of all being constant values (corresponding fixed x-locations).

I would like to create a quiver plot at each x-location to see the velocity arrows in Y-Z plane (refer to image).

I am able to create quiver plots seperately using

velRad_mag=sqrt(slice2(:,5).^2+slice2(:,6).^2);

h1=quiver(slice2(:,2),slice2(:,3),slice2(:,5)./velRad_mag,slice2(:,6)./velRad_mag,'k');

xlabel('y')

ylabel('z')

axis equal

which produces

I was wondering how I could add a third axis which will be the x-axis and have quiver plots in one single figure to see the change of velocity vectors in Y-Z plane along the x-direction.

How can I place my quiver plots with an offset, is there a setting or an object property that controls the let's say x-offset?

Thanks!

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

Bjorn Gustavsson 2022-9-2

1
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1793140-how-to-create-stacked-quiver-plots#answer_1040450

在 MATLAB Online 中打开

This sounds like a straight use of quiver3? Either

h1=quiver3(slice2(:,1),...
           slice2(:,2),...
           slice2(:,3),...
           slice2(:,5)/velRad_mag,...
           slice2(:,6)./velRad_mag,...
           0*slice2(:,6)./velRad_mag,...
           'k');

Or just put the calls into a loop:

dx_offset = 1;
n_slices = 5;
slices_all = cat(3,slice1,slice2,slice3,slice4,slice5); % Since I just realised the "2"
for ix = n_slices:-1:1
  h1(ix) = quiver3(dx_offset*(ix-1)+0*slices_all(:,2,ix),...
                   slices_all(:,2,ix),...
                   slice_all(:,3,ix),...
                   slice_all(:,5,ix)./velRad_mag,...
                   slices_all(:,6,ix)./velRad_mag,...
                   0*slices_all(:,6,ix)./velRad_mag,...
                   'k'); 
end

HTH

2 个评论
显示无隐藏无

doruk isik 2022-9-2

在 MATLAB Online 中打开

I wasn't sure how to use quiver3 as the structure looks like quiver3(X,Y,Z,U,V,W) and I only am intrested in Y,Z and V,W arguments. But never thought of manipulating X and U components as you suggested.

I changed lines 0*slices_all(:,6,ix),.. and slices_all(:,6,ix),... to get what I initially wanted.

hold on
dx_offset = 1;
n_slices = 5;
slices_all = cat(3,slice1,slice2,slice3,slice4,slice5); % Since I just realised the "2"
for ix = n_slices:-1:1
  h1(ix) = quiver3(dx_offset*(ix-1)+0*slices_all(:,2,ix),...
                   slices_all(:,2,ix),...
                   slices_all(:,3,ix),...
                   0*slices_all(:,6,ix),...
                   slices_all(:,5,ix),...
                   slices_all(:,6,ix),...
                   'k'); 
end
view(3)