Matlab模拟旋转圆的动画

原创

matlab科研助手 2022-02-23 12:29:17 ©著作权

文章标签 ide matlab代码参考文献 文章分类 Python 后端开发

©著作权归作者所有：来自51CTO博客作者matlab科研助手的原创作品，请联系作者获取转载授权，否则将追究法律责任

1 简介

Matlab模拟旋转圆的动画

2 完整代码

%% Circle Illusion - Moving dots to rotating circle
% Animation of a rotating circle that emerges from a set of dots moving
% harmonically in straight lines with different orientations.
%
%%
clear ; close all ; clc
%% Parameters
% Video
tf  = 30;                       % Final time                    [s]
fR  = 30;                       % Frame rate                    [fps]
dt  = 1/fR;                     % Time resolution               [s]
t   = linspace(0,tf,tf*fR);     % Time                          [s]
% Dots
T   = 10;                       % Oscillation period            [s]
f   = 1/T;                      % Oscillation frequency         [Hz]
A   = 1;                        % Amplitude                     [m]
N   = 8;                        % Number of dots
ph  = pi/N;                     % Phase                         [rad]
th  = pi/N;                     % Orientation                   [rad]
%% Defining motion
% Preallocating
px = zeros(length(t),N);
py = zeros(length(t),N);
% Position matrices
for i = 1:N
    px(:,i) = A*sin(2*pi*f*t + i*ph) * cos(i*th);
    py(:,i) = A*sin(2*pi*f*t + i*ph) * sin(i*th);
end
%% Animation
figure
% set(gcf,'Position',[50 50 1280 720])  % YouTube: 720p
% set(gcf,'Position',[50 50 854 480])   % YouTube: 480p
set(gcf,'Position',[50 50 640 640])     % Social
hold on ; grid on ; box on ; axis equal
set(gca,'xlim',[-1.1 1.1],'ylim',[-1.1 1.1])
set(gca,'XTick',[],'YTick',[])
set(gca,'FontName','Verdana','FontSize',22)
title('Circle Illusion')
% Create and open video writer object
v = VideoWriter('circle_illusion.mp4','MPEG-4');
v.Quality   = 100;
v.FrameRate = fR;
open(v);
% Generating frames
for i=1:length(t)
    cla
    for j=1:N
        plot([-A*cos(j*th) A*cos(j*th)],[-A*sin(j*th) A*sin(j*th)],':','Color',[100 100 100]/255,'LineWidth',1.5)
    end
    plot(px(i,:),py(i,:),'ro','MarkerFaceColor','r','MarkerEdgeColor','k','MarkerSize',15)
    frame = getframe(gcf);
    writeVideo(v,frame);
end
close(v);