MATLAB Answers

0

Plot Circle In For Loop

Asked by Kelly McGuire on 27 Apr 2018
Latest activity Commented on by Kelly McGuire on 30 Apr 2018
I have this code that is set up to animate a circle plot based on the circle's diameter. I am struggling with how to plot the circle based on the diameter in the for loop. I have one column vector of 1,667 columns. In each of those columns is the diameter (i.e. 1,667 diameters from data). I want to animate the circle changing based on those diameters. So, I need the plot in the for loop to go through, iterate through the columns, and update the plot, and but show that as an animation. I've attached my code. Please help with the plot portion of the code, thanks!

  5 Comments

Ok, I tried the code in that link for a circle, and that code works. Now, the trouble I'm having is with my matrix dimensions. I'm not too good with the iteration coding yet, but here is my updated code. I'm getting the error: Matrix dimensions must agree Error in circle lin 12 Error in AnimatedCirclePlot line 11. I've attached both .m files.
Looks like I forgot all of the (i,j) in the xunit and yunit functions. That fixed the matrix dimensions problem. Now, a new error, one that really don't know how to solve. Index in position 2 exceeds array bounds (must not exceed 1). I know this is talking about the cos(th(i,j)) and sin(th(i,j)) in circle.m, but I'm not sure what the problem is really saying...

Sign in to comment.

1 Answer

Ameer Hamza
Answer by Ameer Hamza
on 27 Apr 2018
Edited by Ameer Hamza
on 27 Apr 2018
 Accepted Answer

@Kelly, As I answered on your previous question, that you need to fill the plot for circle yourself. If you are still not able to make that work then here is the complete code.
B = 10:10:1000; % test B vector. You will use your own B
animationWriter = VideoWriter('HSPtoHSPCircle');
open(animationWriter);
l = line();
l.XData = [];
l.YData = [];
xlim([-max(B) max(B)]);
ylim([-max(B) max(B)]);
theta = 0:pi/50:2*pi;
for i = 1:length(B)
% create data for circle plot or any other
% processing you want to do
xData = B(i)*cos(theta);
yData = B(i)*sin(theta);
l.XData = xData;
l.YData = yData;
frame = getframe(gcf);
writeVideo(animationWriter, frame);
end
close(animationWriter);
I created this animation using given code.

  6 Comments

Hey Ameer, wasn't sure how to direct this question to you in the Matlab Answers section, but it is related to the previous question and your code worked perfectly. Another plot I am trying to make is the same idea as the animated circle changing by diameter, but this time a line plot changing based on the number of water molecules at a given amino acid inside an ion channel. I am trying to modify the code you gave me to do this. Here is what I have so far, but I'm not sure what to put in the xlim or xData, because the x-data is just the name of the residue, so the xdata doesn't change, but a line needs a y-coordinate and an x-coordinate. An example of the plot and the matrix is attached, which is a plot from just the first row of # of water molecules in the 1667 x 4 matrix, as well as the modified code so far. At each of those 4 residues, the # number of waters is changing each frame, but I didn't want to plot just moving points, but an animated line of the # of waters changing. Maybe this same code can't be used for this though.
Ameer Hamza
on 30 Apr 2018
As far as changing the labels on X-axis is concerned, the following code can do it,
ax.XLim = [1 4]
ax.XTick = 1:4
ax.XTickLabel = {'residue 1', 'residue 2', 'residue 3', 'residue 4'}
As far as making an animation with moving line is concerned. It is quite different to this question. Here is an outline of how can you do it.
  • Make a linear interpolation of edges between two rows of your data matrix.
  • Instead of just drawing the rows of your matrix, also plot all the interpolated points. This way, the animation will look smooth.
If you still face a problem, you might want to start a new question, because the new problem is quite out of the scope for this question.
Sounds good, I'll post a new question. Was hoping the same code you gave me for the circle plot could do this.

Sign in to comment.