Генерируют два реальных сигнала: чирп и синусоида.

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

Используйте динамическую деформацию по времени для выравнивания сигналов таким образом, чтобы сумма евклидовых расстояний между их точками была наименьшей. Отображение выровненных сигналов и расстояния.

dtw(x,y);

Измените частоту синусоиды в два раза больше ее начального значения. Повторите вычисления.

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

dtw(x,y);

Добавьте мнимую часть к каждому сигналу. Восстановите начальную частоту синусоиды. Используйте динамическую временную деформацию для выравнивания сигналов путем минимизации суммы возведенных в квадрат евклидовых расстояний.

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

dtw(x,y,'squared');

Разработайте шрифт, напоминающий выходные данные ранних компьютеров. Используйте его для записи слова MATLAB ®.

chr = @(x)dec2bin(x')-48;

M = chr([34 34 54 42 34 34 34]);
A = chr([08 20 34 34 62 34 34]);
T = chr([62 08 08 08 08 08 08]);
L = chr([32 32 32 32 32 32 62]);
B = chr([60 34 34 60 34 34 60]);

MATLAB = [M A T L A B];

Повреждайте слово, повторяя случайные столбцы букв и изменяя интервал. Отображение исходного слова и трех поврежденных версий. Сбросьте генератор случайных чисел для воспроизводимых результатов.

rng('default')

c = @(x)x(:,sort([1:6 randi(6,1,3)]));

subplot(4,1,1,'XLim',[0 60])
spy(MATLAB)
xlabel('')
ylabel('Original')

for kj = 2:4
    subplot(4,1,kj,'XLim',[0 60])
    spy([c(M) c(A) c(T) c(L) c(A) c(B)])
    xlabel('')
    ylabel('Corrupted')
end

Создайте еще две поврежденные версии слова. Выровняйте их с помощью динамического искажения по времени.

one = [c(M) c(A) c(T) c(L) c(A) c(B)];
two = [c(M) c(A) c(T) c(L) c(A) c(B)];

[ds,ix,iy] = dtw(one,two);

onewarp = one(:,ix);
twowarp = two(:,iy);

Отображение неориентированных и выровненных слов.

figure

subplot(4,1,1)
spy(one)
xlabel('')
ylabel('one')

subplot(4,1,2)
spy(two,'r')
xlabel('')
ylabel('two')

subplot(4,1,3)
spy(onewarp)
xlabel('')
ylabel('onewarp')

subplot(4,1,4)
spy(twowarp,'r')
xlabel('')
ylabel('twowarp')

Повторите вычисления, используя встроенные функциональные возможности dtw.

dtw(one,two);

Генерировать два сигнала, состоящих из двух различных пиков, разделенных впадинами разной длины. Постройте график сигналов.

x1 = [0 1 0 0 0 0 0 0 0 0 0 1 0]*.95;
x2 = [0 1 0 1 0]*.95;

subplot(2,1,1)
plot(x1)
xl = xlim;
subplot(2,1,2)
plot(x2)
xlim(xl)

Выровняйте сигналы без ограничений на пути деформации. Чтобы обеспечить идеальное выравнивание, функция должна повторять только одну выборку более короткого сигнала.

figure
dtw(x1,x2);

Постройте график пути деформации и прямой посадки между двумя сигналами. Для достижения выравнивания функция щедро расширяет желоб между пиками.

[d,i1,i2] = dtw(x1,x2);

figure
plot(i1,i2,'o-',[i1(1) i1(end)],[i2(1) i2(end)])

Повторите вычисления, но теперь ограничьте траекторию деформации отклонением максимум трех элементов от прямолинейной посадки. Постройте график растянутых сигналов и траектории деформации.

[dc,i1c,i2c] = dtw(x1,x2,3);

subplot(2,1,1)
plot([x1(i1c);x2(i2c)]','.-')
title(['Distance: ' num2str(dc)])
subplot(2,1,2)
plot(i1c,i2c,'o-',[i1(1) i1(end)],[i2(1) i2(end)])

Ограничение исключает чрезмерную концентрацию деформации на небольшом подмножестве образцов за счет качества выравнивания. Повторите расчет с ограничением в один образец.

dtw(x1,x2,1);

Загрузите речевой сигнал, дискретизированный при $${}_{\mathrm{Fs}}\mathrm{=}\mathrm{}\mathrm{}\mathrm{}$$7418 Гц. Файл содержит запись женского голоса, говорящего слово «MATLAB ®».

load mtlb

% To hear, type soundsc(mtlb,Fs)

Извлеките два сегмента, которые соответствуют двум экземплярам//фонемы. Первый происходит приблизительно между 150 мс и 250 мс, а второй - между 370 мс и 450 мс.

a1 = mtlb(round(0.15*Fs):round(0.25*Fs));
a2 = mtlb(round(0.37*Fs):round(0.45*Fs));

subplot(2,1,1)
plot((0:numel(a1)-1)/Fs+0.15,a1)
title('a_1')
subplot(2,1,2)
plot((0:numel(a2)-1)/Fs+0.37,a2)
title('a_2')
xlabel('Time (seconds)')

% To hear, type soundsc(a1,Fs), pause(1), soundsc(a2,Fs)

Деформируйте временные оси так, чтобы евклидово расстояние между сигналами было сведено к минимуму. Вычислите общую «длительность» искаженных сигналов и постройте их график.

[d,i1,i2] = dtw(a1,a2);

a1w = a1(i1);
a2w = a2(i2);

t = (0:numel(i1)-1)/Fs;
duration = t(end)

duration = 0.1297

subplot(2,1,1)
plot(t,a1w)
title('a_1, Warped')
subplot(2,1,2)
plot(t,a2w)
title('a_2, Warped')
xlabel('Time (seconds)')

% To hear, type soundsc(a1w,Fs), pause(1), sound(a2w,Fs)

Повторите эксперимент с полным словом. Загрузите файл, содержащий слово, «сильное», произнесенное женщиной и мужчиной. Сигналы выбраны на уровне 8 кГц.

load('strong.mat')

% To hear, type soundsc(her,fs), pause(2), soundsc(him,fs)

Деформировать временные оси так, чтобы абсолютное расстояние между сигналами было сведено к минимуму. Постройте график исходных и преобразованных сигналов. Вычислите общую деформированную «длительность».

dtw(her,him,'absolute');
legend('her','him')

[d,iher,ihim] = dtw(her,him,'absolute');
duration = numel(iher)/fs

duration = 0.8394

% To hear, type soundsc(her(iher),fs), pause(2), soundsc(him(ihim),fs)

Файлы MATLAB1.gif и MATLAB2.gif содержат два рукописных образца слова «MATLAB ®». Загрузите файлы и выровняйте их по оси X с помощью динамического временного искажения.

samp1 = 'MATLAB1.gif';
samp2 = 'MATLAB2.gif';

x = double(imread(samp1));
y = double(imread(samp2));

dtw(x,y);