Loop animation with quiver plot

Question

Emil Brandt 2016-5-1

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/281945-loop-animation-with-quiver-plot

回答： Alessandro Masullo 2016-5-3

Hey guys. I am fairly new with matlab. Currently i am trying to animate the change of a velocity field in a fluid in time. The velocity field is dependent on depth (z) and distance (x). I would like to present the velocity field in a x/z coordinate system and by an animation show how it changes with time. I can plot the vector field by constant x with quiver(t,z,u,w). Is it possible to use a for-loop for the animation? I tried a lot of different solutions:

 x=0
 z=linspace(0, -h, 10);
 t=linspace(0, 16*pi, 10);
for i=1:length(t)
    for j=1:length(z)
        u(j,i)=(pi*H)/T*(cosh(k)*(z(j)+h))/(sinh(k*h))*cos(omega*t(i)-k*x); %Particle velocity x direction
        w(j,i)=-(pi*H)/T*(sinh(k)*(z(j)+h))/(sinh(k*h))*sin(omega*t(i)-k*x) %Particle velocity z direction
    end
 end
 figure (4)
 for t=linspace(0, 16*pi, 10);
   for z=linspace(0, -h, 10);
    quiver(x,z,u,w) 
    title('Velocity')
    xlabel('Distance [m]') % x-axis label
    ylabel('Depth [m]') % y-axis label
    ylim([-h 10]) 
    pause(1/60)
   end
 end

5 个评论
显示 3更早的评论隐藏 3更早的评论

KSSV 2016-5-2

many parameters are not defined.....give the values...it is easy to put it in a loop..

Emil Brandt 2016-5-3

All values are defined in the document. I can easily plot the vector field components (u an w) in a time vs. depth diagram for a constant x (x=0) - as done in the attached file. The loop thing just seems not to be able to run the quiver plot in a x vs. depth for somehow.I want matlab to update the attached diagram by the time t, so that the vector field changes by time.

请先登录，再进行评论。

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

Answer 1

Alessandro Masullo 2016-5-3

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/281945-loop-animation-with-quiver-plot#answer_220557

Every time that you want to update the figure, you need to call the drawnow function. You won't probably be able to get a framerate of 60fps, but you may speed up the code by getting rid of the xlabel title functions.

If you want to generate a video, you may have a look at getframe.