www.gusucode.com > signal 案例源码程序 matlab代码 > signal/DynamicTimeWarpingOfRealChirpAndSinusoidExample.m

    %% Dynamic Time Warping of Chirp and Sinusoid
% Generate two real signals: a chirp and a sinusoid.

%%

x = cos(2*pi*(3*(1:1000)/1000).^2);
y = cos(2*pi*9*(1:399)/400);

%%
% Use dynamic time warping to align the signals such that the sum of the
% Euclidean distances between their points is smallest. Display the aligned
% signals and the distance.

dtw(x,y);

%%
% Change the sinusoid frequency to twice its initial value. Repeat the
% computation.

y = cos(2*pi*18*(1:399)/400);

dtw(x,y);

%%
% Add an imaginary part to each signal. Restore the initial sinusoid
% frequency. Use dynamic time warping to align the signals by minimizing
% the sum of squared Euclidean distances.

x = exp(2i*pi*(3*(1:1000)/1000).^2);
y = exp(2i*pi*9*(1:399)/400);

dtw(x,y,'squared');