Matlab для циклической анимации - PullRequest
Новая
       40

Matlab для циклической анимации

1 голос
/ 18 сентября 2011

Я пытаюсь анимировать график функции, но не могу заставить программу отобразить правильные точки.Я хочу построить точки между 0 и 10 и анимировать этот график.Как получить график как функцию времени? 

clear;
w = 2*pi; 
t = 0:.01:10;
y = sin(w*t); 
x = cos(w*t);


for  j=1:10
  plot(x(6*j),y(6*j),'*');
  axis square;
  grid on;
  F(j) = getframe;

 end 
 movie(F,1,1);

Я уточнил код так: 

clear;
 w = 2*pi; 



 for  j=2:11
 t=j-1;
 y = sin(w*t); 
 x = cos(w*t);
  plot(x(t),y(t),'*');
  axis square;
  grid on;
  F(j) = getframe;

 end 
 movie(F);

Это должно делать то, что я пытаюсь, однако теперь я получаю«Индекс превышает размерность матрицы».Как я могу решить это?

Ответы [ 2 ]

1 голос
/ 19 сентября 2011

Вот пример, который показывает анимированную точку вдоль круговой траектории во время записи фильма AVI.

Чтобы узнать больше о создании анимации и записи фильмов в MATLAB, ознакомьтесь с этим гидом . 

%# some parameters
ERASEMODE = 'normal';        %# normal,xor
RENDERER = 'painters';       %# painters,zbuffer,opengl

%# data
t = linspace(0,2*pi,100)';   %'# adjust number of points here
D = [cos(t) -sin(t)];

%# plot circluar path
figure('DoubleBuffer','on', 'Renderer',RENDERER)
plot(D(:,1), D(:,2), 'Color','b', 'LineWidth',2)
grid on, axis([-1.5 1.5 -1.5 1.5]), axis square
xlabel('x'), ylabel('y'), title('Circle Animation')
%#set(gca, 'DrawMode','fast')

%# moving point
hPoint = line('XData',D(1,1), 'YData',D(1,2), 'EraseMode',ERASEMODE,  ...
        'Color','r', 'Marker','.', 'MarkerSize',30);

%# moving coordinates text
hTxtCoords = text(D(1,1), D(1,2), sprintf('(%.2f,%.2f)',D(1,:)), ...
    'Color',[0.2 0.2 0.2], 'FontSize',8, 'EraseMode',ERASEMODE, ...
    'HorizontalAlignment','left', 'VerticalAlignment','top');

%# angle text
hTxtAngle = text(0, 0, num2str(t(1),'%.02f'), ...
    'FontSize',15, 'EraseMode',ERASEMODE, ...
    'HorizontalAlignment','center', 'VerticalAlignment','middle');

%# prepare video output
useVideoWriter = ~verLessThan('matlab','7.11');
if useVideoWriter
    vid = VideoWriter('vid.avi');
    vidObj.Quality = 100;
    vid.FrameRate = 30;
    open(vid);
else
    vid = avifile('vid.avi', 'fps',30, 'quality',100);
end

%# loop
for i=1:numel(t)
    set(hPoint, 'XData',D(i,1), 'YData',D(i,2))    %# update point location
    set(hTxtAngle, 'String',num2str(t(i),'%.02f')) %# update angle text
    set(hTxtCoords, 'Position',D(i,:), ...         %# update angle text
        'String',sprintf('(%.3f,%.3f)',D(i,:)))
    drawnow                                        %# force refresh
    if ~ishandle(hPoint), break; end               %# if you close the figure

    %# capture frame
    if useVideoWriter
        writeVideo(vid,getframe);
    else
        vid = addframe(vid, getframe(gcf));
    end 
end

%# close and save video output
if useVideoWriter
    close(vid);
else
    vid = close(vid);
end

%# open AVI file using system default player
winopen('vid.avi')

screenshot_animation

0 голосов
/ 18 сентября 2011

Он делает именно то, что вы просите.Вы субсэмплируете x и y, поэтому это выглядит довольно забавно.

Попробуйте 

plot(x,y);
axis square;
ax = axis;  

for jx = 1 : length(t),
    plot(x(ix), y(ix), '*');
    axis(ax);  grid on;
    F(ix) = getframe;
end

movie(F, 1, 1/(t(2)-t(1)))

Я бы также использовал t = 1 : 0.1 : 10;, чтобы он отображал скорость 10 FPSвместо 100. Понижение частоты, скажем, до w = pi; также будет более плавным.

В конце концов, Matlab просто не является отличным анимационным решением.

Ответ на уточненный код вопроса

Вам нужно будет использовать plot(x,y);, но это покажет еще одну ошибку - ваш индекс кадра не начинается с 1. Он засорит F(j) впервая итерация, где j = 2.Почему бы просто не зациклить t = 1 : 10?

...