Lissajous曲线动画演示(Matlab实现)

在别人的博文里看到 Lissajous曲线的动态图, 感觉还挺帅,于是自己也想玩一下。

最终效果

Lissajous曲线

在数学上,Lissajous曲线(又称Lissajous图形或Bowditch曲线)是由如下参数方程构成的曲线族的图形 [1]

![参数方程][Equation]
[Equation]: http://latex.codecogs.com/svg.latex?\begin{cases}x=A\sin(at+\varphi)\y=B\sin(bt)\end{cases}

纳撒尼尔·鲍迪奇(Nathaniel Bowditch)在1815年首先研究了这一曲线族,之后在1857年由朱尔斯·安东尼·利萨如(Jules Antoine Lissajous)进行了更详细的研究。

在物理上,Lissajous曲线可以看作是一个质点同时在X轴和Y轴方向上做简谐运动形成的运动轨迹(以我仅剩的一点物理知识来说,这大概是叫"运动的合成")。
当这两个相互垂直方向上的简谐振动的频率为任意值时,那么它们合成的运动可能会比较复杂,轨迹是不稳定的,示意图如下。(关于稳定,我的理解大概就是呈周期性吧)

一个示意图

而如果这两个振动的频率成简单的整数比,就能合成一个规则的、稳定的闭合曲线,这就是Lissajous曲线(示意图如下)。

另一个示意图

具体实现

  • 环境: Matlab R2017a

Animated Line 实现动画效果

首先用animatedline创建初始动画线条对象 [2] ,接着在循环中动态地向线条中添加点,并使用 drawnow 在屏幕上显示该新添加的点,然后用getframe捕获的当前图像。

%%
% Create a new figure window,
% Background color is white
% MenuBar is hidden
% Image size is 480*480
figure('Color','w','MenuBar','none','InnerPosition',[0 0 480 480]);
co = [65/255, 131/255, 196/255]; % specify the color
% Create animated line object
headMarker = animatedline('LineStyle','none','Marker','o','Color',co,'MarkerFaceColor',co);
body = animatedline('LineStyle','-','LineWidth',1.3,'Color',co); % trajectory
axis([-1,1,-1,1]); axis off; % hide axis

%% 
% Parameters
a = 13;
b = 18;
phi = pi * 0;

% Parametric Equation of Lissajous Curve
t = linspace(0,2*pi,1000); % generate 1000 linearly equally spaced points
x = sin(a*t + phi);
y = sin(b*t);

%%
% Display
for idx = 1:length(t)
    addpoints(headMarker,x(idx),y(idx));
    addpoints(body,x(idx),y(idx));
    drawnow 
    % Capture the current plot
    frame = getframe; % snapshot of the current axes
    im{idx} = frame2im(frame); % write image data to `im`
    % Remove previous head marker except for the last one
    if idx~=length(t)
        clearpoints(headMarker);
    end
end

写入GIF动图

getframe捕获的每一个图像,写入 GIF 动画文件 [3]

%%
% Save images into a GIF file. 
% Because three-dimensional data is not supported for GIF files, 
% call rgb2ind to convert the RGB data in the image to an indexed image A with a colormap map. 
filename = 'Lissajous.gif'; % specify the output file name
for idx = 1:length(t)
    [A,map] = rgb2ind(im{idx},256);
    if idx == 1
        % Set DelayTime = 0 to display next images as fast as your hardware allows
        imwrite(A,map,filename,'gif','LoopCount',Inf,'DelayTime',0);
    else
        % To append multiple images to the first image
        imwrite(A,map,filename,'gif','WriteMode','append','DelayTime',0);
    end
end

完整代码链接

参考

  1. http://www.matrix67.com/blog/archives/6947
  2. https://codepen.io/handsomeone/full/Kgmbqg
  3. https://cn.mathworks.com/matlabcentral/fileexchange/28987-lissajous-curve
  4. https://cn.mathworks.com/examples/matlab/community/22655-lissajous-curves
  5. http://www.baike.com/wiki/%E5%88%A9%E8%90%A8%E5%A6%82%E5%9B%BE%E5%BD%A2

推荐阅读更多精彩内容

  • 天还很亮 夜晚却快来临 我害怕孤寂 也害怕星星的眼睛 那些眼睛都是你 让我无法躲避 在你的眼睛里 我失去了我自己 ...
    潜水的鸟阅读 22评论 0 3
  • 长大,是个名词,它是骨骼的伸张,是面容的成熟,更是处事的不惊,待人对事的的从容和执着。 随着年龄的增长,我们其实丢...
    er新心弍阅读 38评论 0 1
  • 这是“心灵对话写作”小组的第四篇文章 今天到现在提笔想写点什么,除了雾霾似乎想不到其它。户外的活动,让我感到眼睛和...
    墨芷轩阅读 27评论 3 2
  • 首先仰卧到垫子上,双臂体侧平举,双手掌心朝下,放到地垫上,继续找站立山式的感觉,把脚后跟蹬出去
    随笔s阅读 575评论 0 0