N_pendulum_multi_mv

%% 绘图
fig=figure;fig.ToolBar='none';
% fig.Renderer='painters';% 提升2D渲染性能
dpi=1.25;scr=[1920 1080];res=[1600 900];
fig.Color='k';fig.Position=[(scr-res)/2,res]/dpi;
pos1=[0 0 1 1];
subplot('Position',pos1)
plot(-N,-N,N,0.25*N);hold on

axis equal;axis tight;axis off
ax=gca;ax.Clipping='off';
n=1;
for ii=1:M
    p{ii}=plot(x(n,2:end,ii),y(n,2:end,ii),'.',MarkerSize=30,Color='y');
    s{ii}=plot(x(n,:,ii),y(n,:,ii),'-',LineWidth=1,Color='w');
end
%% 动画
mv = 1;% 是否录制视频，1录制，其它不录制
if mv==1
v=VideoWriter('D:\Programs\MATLAB\三摆100.mp4','MPEG-4');
v.Quality=100;v.FrameRate=60;open(v);
end
for n=1:Frame
    for k=1:M
        p{k}.XData=x(n,2:end,k);
        p{k}.YData=y(n,2:end,k);
        s{k}.XData=x(n,:,k);
        s{k}.YData=y(n,:,k);
    end
    axis([-N N -N 0.25*N])
    drawnow limitrate
    if mv==1;F=getframe(gcf);writeVideo(v,F);end
    clc;fprintf(['当前进度', ...
        num2str(100*n/Frame,'%.1f'),'%%\n']);
end
if mv==1;close(v);end