Создание анимации эволюции диаграммы Вороного

@Elena_Vasilieva · Регистрация: 30.05.2012

Author24 — интернет-сервис помощи студентам

помогите разобраться, что не так.
Задача: изначально есть 10 файлов с координатами точек [x y] в каждой строчке; каждый файл на конкретный момент времени t
Мне нужно получить анимацию эволюции некоторого параметра psi соответствующего ячейке Вороного конкретной точки в каждый момент времени.

я вначале пишу функцию которая рисует цветную диаграмму Вороного в зависимости от psi и на выходе дает дескриптор текущего графика

Matlab M

function[h]= psi6_voronoi(name)
 
[a, psi]= psi6_local(name);        % функция выдает исходные координаты в массиве a и параметр psi       соответствующий данной точке
 
[V,C]=voronoin([a(:,1),a(:,2)]);
 
 
 
axis([0.1 1 0.1 1]);
for k=1:length(C)    
    if all(C{k}~=1)             % конечные многоугольники рисуем полигональными объектами
       VertCell = V(C{k},:);
    
        if psi(k,2)==6
        patch('XData',VertCell(:,1),'YData',VertCell(:,2),'FaceColor','g','FaceAlpha',1)
        
        elseif psi(k,2)<6
        patch('XData',VertCell(:,1),'YData',VertCell(:,2),'FaceColor','r','FaceAlpha',1)
        
         elseif psi(k,2)>6
        patch('XData',VertCell(:,1),'YData',VertCell(:,2),'FaceColor','b','FaceAlpha',1)
        end
    end
end
   h=gcf; 
    
end

а дальше начинаются трудности. Когда я пытаюсь все картинки свести в анимацию:
clear all

for j=1:10
name='j'
[h]= psi6_voronoi(name);

F(j) = getframe(h);

end

movie(F)

У меня во первых на анимации все кадры накладываются друг на дружку, и появляются две системы координат сдвинутые друг относительно друга после команды movie.
Подскажите как мне очищать окно графика прежде чем там будет нарисована диаграмма Вороного в следующий момент времени и как сделать чтобы оси не дублировались?

Зосима · 02.06.2012, 18:58

Есть у меня несколько догадок, но они останутся догадками, пока сам не попробую поиграться. Можете скинуть мне на мыло zosima0@gmail.com файлы с данными и m-файл с программой? Я во вторник доберусь до матлаба и гляну

P.S.: диапазон значений на графике у вас задается ф-цией axis(...), однако для того, чтобы оно работало для plot, (собственно и grid on) нужно поместить после ф-ции, в данном случае, patch (правда я не уверен, т.к. с ней не работал) тогда все графики будут строиться в одинаковых систеиах координат и прыгать не должно.
PPS: чтобы графики не налазили один поверх другого при анимации можно попробовать поставить после F(j)=getframe(h); clear h; или всиавлять пустые фреймы между графиками:

Matlab M

for j = 1:20
    name = num2str( (j+1)/2 );
    if mod(j, 2) == 0 % проверка на четность
         h = [];
    else 
        h = psi_voron...
    end
    F(j) = getframe(j);
    clear h;
end
movie(F)

Точнее смогу сказать только во вторник.

@Elena_Vasilieva · 02.06.2012, 20:00 **[ТС]**

спасибо) на самом деле я пошла немного другим путем и все решилось. Просто не стала выводить все графики в одном окне, а в цикле открывала и закрывала полученную картинку с запоминанием дескриптора на нее. Тем не менее спасибо за помощь)
А если самому интересно, то можно например задать рандомные координаты - суть от этого не изменится.

@Elena_Vasilieva 1 / 1 / 0 Регистрация: 30.05.2012 Сообщений: 22
	02.06.2012, 20:00 [ТС]	3
	спасибо) на самом деле я пошла немного другим путем и все решилось. Просто не стала выводить все графики в одном окне, а в цикле открывала и закрывала полученную картинку с запоминанием дескриптора на нее. Тем не менее спасибо за помощь) А если самому интересно, то можно например задать рандомные координаты - суть от этого не изменится. 0