我正在使用MATLAB的movie()
函数来制作包含大量时间序列数据的电影。使用当前实例化,大约需要14个小时才能完成。我该怎么做才能更好地对其进行优化?
我想要做的最大的事情是抑制每次将绘图绘制到屏幕上,同时仍然为getframe()函数更新它。
我的代码的缩写版本是:
t = 0:1000000; % I have about 10^6 data points
x = sin(t); % Let's pretend the data I'm plotting is the sine function
y = cos(t); % I have multiple data series being plotted in 'snapshots'
num_frames = length(t);
movie_length = 100; % seconds
fps = 60;
for k = 1:num_frames
clf; hold on;
plot(1, x(k),'bo')
plot(2, y(k),'bo')
ax = gca;
ax.XLim = [0 1];
ax.YLim = [-1,1];
f(k) = getframe;
end
movie(f,1,fps)
答案 0 :(得分:1)
这里是一个版本,可以在我的计算机上将其速度提高大约两倍。我已将点数减少到10^3
。
clear f g;
t = 0:10^3; % I have about 10^6 data points
x = sin(t); % Let's pretend the data I'm plotting is the sine function
y = cos(t); % I have multiple data series being plotted in 'snapshots'
num_frames = length(t);
tic;
for k = 1:num_frames
clf; hold on;
plot(1, x(k),'bo')
plot(2, y(k),'bo')
ax = gca;
ax.XLim = [.8 2.2];
ax.YLim = [-1,1];
f(k) = getframe();
end
toc
% This is faster
h_fig = figure;
g(num_frames) = struct('cdata', [], 'colormap', []);
tic;
ax = axes(h_fig);
ax.XLim = [.8 2.2];
ax.YLim = [-1,1];
hold on;
p1 = plot(1,x(1), 'bo');
p2 = plot(2,y(1), 'bo');
drawnow;
g(1) = getframe(ax);
for k = 2:num_frames
p1.YData = x(k);
p2.YData = y(k);
g(k) = getframe(ax);
end
toc
请注意,如果没有getframe
,则第二个版本的速度大约要快100倍。因此,如果您知道如何计算单个帧的cdata
,它可能比绘制数据并使用getframe
要快得多。
一句话:我无法正确地复制,粘贴和运行您的代码。即使记起来很快,也可以在发布前对其进行错误测试。