How do I make just one plot stay in a figure, while other plots are removed/updated as I iterate through my loop

Question

John Klint 2017-3-2

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/327805-how-do-i-make-just-one-plot-stay-in-a-figure-while-other-plots-are-removed-updated-as-i-iterate-thr

评论： John Klint 2017-3-6

I'm trying to make the plot of a potential visible in my figure while I plot the time evolution of the wave function of a particle in the potential.

 V=diag(xvec.^2);  %V is potential and xvec is a vector
 figure(1)
 hold on
 plot(xvec,V);

And after this point I want "hold on" command to be turned off for the animation to make sense, but the plot of the potential V to still appear in my figure, while the loop below iterates through the time evolution.

 for t=0:dt:100*dt
     v=exp(-i*E(1:N)*t/hbar);
     Psi=v'*ev(:,1:N)'/sqrt(N);
     plot(xvec,abs(Psi).^2/dx)
     axis([-0.1*a 1.1*a 0 4])
     text(2.5*a,0.45,num2str(t))
     pause(.2)
 end

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

请先登录，再进行评论。

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

Answer 1

Walter Roberson 2017-3-2

1
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/327805-how-do-i-make-just-one-plot-stay-in-a-figure-while-other-plots-are-removed-updated-as-i-iterate-thr#answer_257033

在 MATLAB Online 中打开

   V=diag(xvec.^2);  %V is potential and xvec is a vector
   fig = figure(1);
   ax = axes('Parent', fig);
   hold(ax, 'on')
   plot(ax, xvec, V);
   for t=0:dt:100*dt
       v = exp(-i*E(1:N)*t/hbar);
       Psi = v'*ev(:,1:N)'/sqrt(N);
       y = abs(Psi).^2/dx;
       if t == 0
          ph = plot(ax, xvec, y)
          axis(ax, [-0.1*a 1.1*a 0 4])
          th = text(ax, 2.5*a, 0.45, num2str(t));
       else
         set(ph, 'YData', y);
         set(th, 'String, num2str(t));
       end
       pause(.2)
   end
   hold(ax, 'off')