Code returns, "Error using quiver3 (line 44) V and W must be the same size." How to fix this? The variables are the same size(1x25 double).

Question

Aamna Jangda 2017-10-25

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/363249-code-returns-error-using-quiver3-line-44-v-and-w-must-be-the-same-size-how-to-fix-this-the-va

回答： Mischa Kim 2017-10-26

%% Program to plot Ekman spiral

%% Input section

clear all; close all;

i=sqrt(-1) % To make sure i is the complex i

h=120; % Water depth

nut=.02; % (constant) turbulence viscosity

phi=51; % Latitude (deg.)

rho=1000; % Water density

tauwx=0; % Wind shear stress x-direction

tauwy=1; % Wind shear stress y-direction

tau=tauwx+tauwy*i; % Shear stress vector

f=2*7.27e-5*sin(phi*pi/180); % Coriolis coefficient

dek=sqrt(2*nut/f) % Ekman depth

dz=5; % Vertical step size

zek=[-h:dz:0]; % Vertical coordinate

%% Define some additional matrices for plotting purposes

xek=zeros(size(zek));

yek=zeros(size(zek));

taux=zeros(size(zek));

tauy=zeros(size(zek));

taux(end)=tauwx;

tauy(end)=tauwy;

%% Solution of Ekman spiral

s=tau*dek/rho/(1+i)/nut*exp((1+i)/dek*zek);

uek=real(s);

vek=imag(s);

%% Plot results

figure(1);

quiver3(xek,yek,zek,uek,vek,0,'linewidth',2);

hold on

quiver3(xek,yek,zek,taux,tauy,0,'r','linewidth',2);

title('Ekman spiral')

xlabel('u');ylabel('v');

zlabel('z')

print('-dpng','ekman.png')

print('-depsc','ekman.eps')

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

请先登录，再进行评论。

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

Answer 1

Mischa Kim 2017-10-26

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/363249-code-returns-error-using-quiver3-line-44-v-and-w-must-be-the-same-size-how-to-fix-this-the-va#answer_287806

在 MATLAB Online 中打开

Aamna, based on your information I am not clear exactly what you want to do. However, the input arguments for coordinates and vectors in quiver3 all must be the same size. You can use meshgrid to achieve that:

%%Solution of Ekman spiral
[XEK,YEK,ZEK] = meshgrid(xek,yek,zek);
s   = tau*dek/rho/(1+i)/nut*exp((1+i)/dek*ZEK);
uek = real(s);
vek = imag(s);
%%Plot results
figure(1);
quiver3(XEK,YEK,YEK,uek,vek,zeros(size(uek)),'linewidth',2);