%%?make?the?example?movie
%?parameters
nx?=?80;?%?horizontal?pixels
ny?=?80;?%?vertical?pixels
n??=?nx?*?ny;
T?=?10;?%?sec
dt?=?0.01;
t?=?dt:dt:T;
sig?=?0.1;?%?magnitude?of?gaussian?noise?added

%?3?underlying?modes
[Xgrid?Ygrid]?=?meshgrid（1:nx?1:ny）;

mode1.u?=?exp（-（（Xgrid-40）.^2/250+（Ygrid-40）.^2/250））;
mode1.f?=?5.55;?%?cycles/sec
mode1.A?=?1;
mode1.lambda?=?exp（1i?*?mode1.f*2*pi*dt）;
mode1.range?=?[0?5];

mode2.u?=?zeros（nx?ny）;?
mode2.u（nx-40:nx-10?ny-40:ny-10）?=?1;
mode2.f?=?0.9;?%?cycles/sec?
mode2.A?=?1;
mode2.lambda?=?exp（1i?*?mode2.f*2*pi*dt）;
mode2.range?=?[3?7];

mode3.u?=?zeros（nx?ny）;?
mode3.u（1:nx-20?1:ny-20）?=?1;
mode3.f?=?0.15;?%?cycles/sec?
mode3.A?=?0.5;
mode3.lambda?=?exp（1i?*?mode3.f*2*pi*dt）;
mode3.range?=?[0?T];

modes（1）?=?mode1;
modes（2）?=?mode2;
modes（3）?=?mode3;

%?make?the?movie
Xclean?=?zeros（n?numel（T））;
for?ti?=?1:numel（t）
????Snap?=?zeros（nx?ny）;
????for?mi?=?1:numel（modes）
????????mymode?=?modes（mi）;
????????if?ti?>?round（mymode.range（1）/dt）?&&?...
???