多分辨率动态模式分解

%% 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） && ...