我想在3d中模拟粒子的轨迹,因此我使用以下代码创建了这个小模拟:
nP=100;
N=100;
z=rand(nP,N);
y=rand(nP,N); % Compute coordinates of particles for each iteration NS
z=rand(nP,N);
f = figure
view(3);
for i=N
h = animatedline('MaximumNumPoints', 1.e4,'color',rand(1,3));
for k = 1:length(x)
addpoints(h,x(i,k),y(i,k),z(i,k));
drawnow
end
end
hold on
hold off
numpoints = 500;
y2 = 3 +square(x+1);
f = figure
h = animatedline('Color','b','LineWidth',2);
h2 = animatedline('Color','r','LineWidth',2);
grid on;
%axis([0,12,-3,+6])
for k = 1:N
addpoints(h,x(k),y(k),z(k))
%addpoints(h2,x(k),y2(k))
drawnow
% Capture the plot as an image
frame = getframe(f);
im = frame2im(frame);
[imind,cm] = rgb2ind(im,256);
% Write to the GIF File
if k == 1
imwrite(imind,cm,'test.gif','gif', 'Loopcount',inf);
else
imwrite(imind,cm,'test.gif','gif','WriteMode','append');
end
end
我正在尝试将其另存为gif或mp4。一切正常,但最后保存的是2d中粒子的轨迹。关于如何进行这项工作的任何想法?
您在代码中有一个错误:
z=rand(nP,N);
y=rand(nP,N);
z=rand(nP,N);
应该是(请参见第一行):
x=rand(nP,N);
y=rand(nP,N); % Compute coordinates of particles for each iteration NS
z=rand(nP,N);
在第二部分view(3)
之前drawnow
...
for k = 1:N
addpoints(h,x(k),y(k),z(k))
%addpoints(h2,x(k),y2(k))
view(3);
drawnow
...
完整的代码对我有用:
nP=100;
N=100;
x=rand(nP,N);
y=rand(nP,N); % Compute coordinates of particles for each iteration NS
z=rand(nP,N);
f = figure
view(3);
for i=N
h = animatedline('MaximumNumPoints', 1.e4,'color',rand(1,3));
for k = 1:length(x)
addpoints(h,x(i,k),y(i,k),z(i,k));
drawnow
end
end
hold on
hold off
numpoints = 500;
y2 = 3 +square(x+1);
f = figure
h = animatedline('Color','b','LineWidth',2);
h2 = animatedline('Color','r','LineWidth',2);
grid on;
%axis([0,12,-3,+6])
for k = 1:N
addpoints(h,x(k),y(k),z(k))
%addpoints(h2,x(k),y2(k))
view(3);
drawnow
% Capture the plot as an image
frame = getframe(f);
im = frame2im(frame);
[imind,cm] = rgb2ind(im,256);
% Write to the GIF File
if k == 1
imwrite(imind,cm,'test.gif','gif', 'Loopcount',inf);
else
imwrite(imind,cm,'test.gif','gif','WriteMode','append');
end
end
本文收集自互联网,转载请注明来源。
如有侵权,请联系[email protected] 删除。
我来说两句