Загрузите в данных о землетрясении Кобе. Данные являются измерениями сейсмографа (вертикальное ускорение в $\text{nm}/{\text{секунда}}^2$) зарегистрированный в Университете Тасмании, Хобарт, Австралия, 16 января 1995, начинаясь в 20:56:51 (GMT) и продолжаясь в течение 51 минуты во вторых интервалах.

load kobe

Откройте Signal Multiresolution Analyzer и нажмите Load Signal. Окно кажется перечисляющим все переменные рабочей области, которые может обработать приложение.

Выберите данные Кобе из диалогового окна и нажмите ОК. Появляется четырехуровневое разложение MODWTMRA сигнала. Анализируемый сигнал называют kobe1 в Анализируемой панели Сигналов. Суффиксный [modwtmra] идентифицирует разложение как базирующийся вейвлет. По умолчанию графики относительно демонстрационного индекса, и частоты находятся в циклах на выборку. Установите переключатель Sample Rate. Поскольку частота дискретизации данных составляет 1 герц, вы не должны изменять значение по умолчанию. Графики и частоты обновляются, чтобы использовать частоту дискретизации.

Графики в средней панели Разложения являются проекциями разложений вейвлета сигнала в каждой шкале на исходном подпространстве сигнала. Исходный сигнал, kobe, и реконструкция, kobe1, построены в панели Реконструкции. Панель Выбора Уровня показывает относительные энергии сигнала через шкалы, а также диапазоны частот.

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

Изменение любых настроек в панели инструментов включит кнопку Decompose.

Загрузите шумный Доплеровский сигнал. Сигнал является шумной версией Доплеровского тестового сигнала Донохо и Джонстона [1].

load noisdopp

Откройте Signal Multiresolution Analyzer и загрузите сигнал в приложение. По умолчанию приложение создает четырехуровневое разложение MODWTMRA сигнала. В Анализируемой панели Сигналов разложение вейвлета называют noisdopp1. Панель Реконструкции отображает исходные и восстановленные сигналы, построенные двумя различными цветами.

Чтобы добавить разложение EMD, нажмите Add ▼ и выберите EMD.

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

Чтобы более легко видеть различия между этими двумя реконструкциями, нажмите noisdopp в легенде графика. Текст исчезает, и график исходного сигнала скрыт. Можно использовать легенду, чтобы скрыть любой график в панели Реконструкции.

Можно изменить параметры в панели инструментов, чтобы сгенерировать различное разложение EMD. Установка курсора на параметр покажет подсказку.

Чтобы узнать больше о параметрах и алгоритме EMD, смотрите emd.

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

Загрузите данные о землетрясении Кобе в свою рабочую область. Данные являются измерениями сейсмографа (вертикальное ускорение в $\text{nm}/{\text{секунда}}^2$) зарегистрированный в Университете Тасмании, Хобарт, Австралия, 16 января 1995, начинаясь в 20:56:51 (GMT) и продолжаясь в течение 51 минуты во вторых интервалах.

load kobe

Откройте Signal Multiresolution Analyzer и загрузите данные о землетрясении в приложение. По умолчанию приложение создает четырехуровневое разложение MODWTMRA сигнала под названием kobe1 использование порядка 4 Symlet sym4. Кликните по переключателю Sample Rate так, чтобы графики были относительно времени.

Создайте новое шестиуровневое разложение с помощью порядка 4 Coiflet. Нажмите Duplicate в панели инструментов. Начиная с kobe1 в настоящее время выбранный пункт в Анализируемых Сигналах, копия первого разложения создается. Копия называется kobe1Copy. Графики в Реконструкции обновляются, чтобы включать новое разложение. За исключением цвета, копия будет идентична с первым разложением. Можно поменять имя копии путем щелчка правой кнопкой по имени в Анализируемых Сигналах.

Измените настройки в панели инструментов к следующим значениям и затем нажмите Decompose.

В Выборе Уровня отметьте, какие компоненты разложения включены в реконструкцию: приближение и детали уровня 5 и уровня 6.

Уровень 4 имеет приблизительно 60% полной энергии. Удалите уровни 5 и 6 из реконструкции, и вместо этого включайте уровень 4. Покажите только приближение и детали уровня 4 в панели Разложения. Чтобы приблизительно выровнять разложение с реконструкцией, перетащите панель Разложения ниже панели Реконструкций.

У вас есть три опции экспорта. Можно экспортировать реконструкцию или целое разложение к рабочей области, или можно создать скрипт MATLAB™. Чтобы сгенерировать скрипт, нажмите Export> Generate MATLAB Script.

Неназванный скрипт открывается в вашем редакторе следующим исполняемым кодом:

% Logical array for selecting reconstruction elements
levelForReconstruction = [false, false, false, true, false, false, true];
% Perform the decomposition using modwt
wt = modwt(kobe, 'coif4', 6);
% Construct MRA matrix using modwtmra
mra = modwtmra(wt, 'coif4');
% Sum along selected multiresolution signals
kobe1Copy = sum(mra(levelForReconstruction,:),1);

Истинно-ложные значения в levelForReconstruction соответствуйте который Include флажки устанавливаются в Выборе Уровня. Можно сохранить скрипт, как, или измените его, чтобы применить те же настройки разложения к другим сигналам. Запустите скрипт и постройте исходный сигнал и реконструкцию. За исключением возможно цветов, график будет совпадать с kobe1Copy реконструкция показана в приложении.

t = 0:numel(kobe)-1;
plot(t,kobe)
grid on
hold on
plot(t,kobe1Copy,'LineWidth',2)
xlabel('Seconds')
title('Reconstruction')
legend('Original','Reconstruction','Location','northwest')
axis tight