Загрузите файл данных, содержащий запись свистка train, выбранного с частотой 8192 Гц. Найдите 10 точек, в которых среднеквадратичный уровень сигнала изменяется наиболее значительно.

load train

findchangepts(y,'MaxNumChanges',10,'Statistic','rms')

Вычислите кратковременную спектральную плотность степени сигнала. Разделите сигнал на сегменты с 128 выборками и окончите каждый сегмент с Окном Хэмминга. Задайте 120 выборки перекрытия между смежными сегментами и 128 точками ДПФ. Найдите 10 точек, в которых среднее значение спектральной плотности степени изменяется наиболее значительно.

[s,f,t,pxx] = spectrogram(y,128,120,128,Fs);

findchangepts(pow2db(pxx),'MaxNumChanges',10)

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

rng('default')

lr = 20;

mns = [0 1 4 -5 2 0 1];
nm = length(mns);

vrs = [1 4 6 1 3];
nv = length(vrs);

v = randn(1,lr*nm*nv)/2;

f = reshape(repmat(mns,lr*nv,1),1,lr*nm*nv);
y = reshape(repmat(vrs,lr*nm,1),1,lr*nm*nv);

t = v.*y+f;

Постройте график сигнала, выделяя шаги его конструкции.

subplot(2,2,1)
plot(v)
title('Original')
xlim([0 700])

subplot(2,2,2)
plot([f;v+f]')
title('Means')
xlim([0 700])

subplot(2,2,3)
plot([y;v.*y]')
title('Variances')
xlim([0 700])

subplot(2,2,4)
plot(t)
title('Final')
xlim([0 700])

Найдите пять точек, где среднее значение сигнала изменяется наиболее значительно.

figure
findchangepts(t,'MaxNumChanges',5)

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

findchangepts(t,'MaxNumChanges',5,'Statistic','rms')

Найдите точку, где среднее и стандартное отклонение сигнала изменяются больше всего.

findchangepts(t,'Statistic','std')

Загрузите речевой сигнал, дискретизированный $$F_s=7418Hz$$. Файл содержит запись женского голоса, говорящего слово «MATLAB ®».

load mtlb

Разберите гласные и согласные в слове, найдя точки, в которых отклонение сигнала значительно изменяется. Ограничьте количество точек изменения пятью.

numc = 5;

[q,r] = findchangepts(mtlb,'Statistic','rms','MaxNumChanges',numc)

q = 5×1

         132
         778
        1646
        2500
        3454

r = -4.4055e+03

Постройте график сигнала и отобразите точки изменения.

findchangepts(mtlb,'Statistic','rms','MaxNumChanges',numc)

Чтобы воспроизвести звук с паузой после каждого из сегментов, раскомментируйте следующие линии.

% soundsc(1:q(1),Fs)
% for k = 1:length(q)-1
%     soundsc(mtlb(q(k):q(k+1)),Fs)
%     pause(1)
% end
% soundsc(q(end):length(mtlb),Fs)

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

vc = sin(2*pi*(0:201)/17).*sin(2*pi*(0:201)/19).* ...
    [sqrt(0:0.01:1) (1:-0.01:0).^2]+(0:201)/401;

Найдите точки, где средний сигнал изменяется наиболее значительно. The 'Statistic' в этом случае пара "имя-значение" является необязательной. Задайте минимальное улучшение остаточной ошибки 1.

findchangepts(vc,'Statistic','mean','MinThreshold',1)

Найдите точки, где среднеквадратичный уровень сигнала изменяется больше всего. Задайте минимальное улучшение остаточной ошибки 6.

findchangepts(vc,'Statistic','rms','MinThreshold',6)

Найдите точки, где стандартное отклонение сигнала изменяется наиболее значительно. Задайте минимальное улучшение остаточной ошибки 10.

findchangepts(vc,'Statistic','std','MinThreshold',10)

Найдите точки, где среднее значение и наклон сигнала изменяются наиболее резко. Задайте минимальное улучшение остаточной ошибки 0,6.

findchangepts(vc,'Statistic','linear','MinThreshold',0.6)

Сгенерируйте двумерную кривую Безье с 1000 образцами с 20 случайными контрольными точками. Кривая Безье определяется:

$$C(t)={\displaystyle \sum _{k=0}^m\left(\genfrac{}{}{0ex}{}{m}{k}\right)t^k(1-t{)}^{m-k}{P}_k}$$,

где ${\mathit{P}}_{\mathit{k}}$ является $\mathit{k}$th из $\mathit{m}$ контрольные точки, $\mathit{t}$ находится в диапазоне от 0 до 1, и $$\left(\genfrac{}{}{0ex}{}{m}{k}\right)$$ является биномиальным коэффициентом. Постройте график кривой и контрольных точек.

m = 20;
P = randn(m,2);
t = linspace(0,1,1000)';

pol = t.^(0:m-1).*(1-t).^(m-1:-1:0);
bin = gamma(m)./gamma(1:m)./gamma(m:-1:1);
crv = bin.*pol*P;

plot(crv(:,1),crv(:,2),P(:,1),P(:,2),'o:')

Разделите кривую на три сегмента, так что точки в каждом сегменте находятся на минимальном расстоянии от среднего сегмента.

findchangepts(crv','MaxNumChanges',3)

Разделите кривую на 20 сегментов, которые лучше подгонка прямым линиям.

findchangepts(crv','Statistic','linear','MaxNumChanges',19)

Сгенерируйте и постройте график трехмерной кривой Безье с 20 случайными контрольными точками.

P = rand(m,3);
crv = bin.*pol*P;

plot3(crv(:,1),crv(:,2),crv(:,3),P(:,1),P(:,2),P(:,3),'o:')
xlabel('x')
ylabel('y')

Визуализируйте кривую сверху.

view([0 0 1])

Разделите кривую на три сегмента, так что точки в каждом сегменте находятся на минимальном расстоянии от среднего сегмента.

findchangepts(crv','MaxNumChanges',3)

Разделите кривую на 20 сегментов, которые лучше подгонка прямым линиям.

findchangepts(crv','Statistic','linear','MaxNumChanges',19)