Измените коэффициенты разложения мультисигнала 1-D
DEC = chgwdeccfs(DEC,'ca',COEFS)
DEC = chgwdeccfs(DEC,'cd',COEFS,LEV)
DEC = chgwdeccfs(DEC,'all',CA,CD)
DEC
= chgwdeccfs(DEC,'all',V)
DEC = chgwdeccfs(...,IDXSIG)
DEC = chgwdeccfs(DEC,'ca',COEFS) заменяет коэффициенты приближения на уровне DEC.level с содержавшимися в матричном COEFS. Если COEFS одно значение V, все коэффициенты заменяются V.
DEC = chgwdeccfs(DEC,'cd',COEFS,LEV) заменяет коэффициенты детали на уровне LEV с содержавшимися в матричном COEFS. Если COEFS одно значение V, затем LEV может быть вектор уровней и всех коэффициентов, которые принадлежат этим уровням, заменяются V. LEV должно быть таково что 1 ≤ LEV ≤ DEC.level
DEC = chgwdeccfs(DEC,'all',CA,CD) замены все приближение и коэффициенты детали. CA должна быть матрица и CD должен быть массив ячеек длины DEC.level.
Если COEFS (или CA или CD) один номер, затем он заменяет все связанные коэффициенты. В противном случае, COEFS (или CA, или CD) должна быть матрица соответствующего размера.
Для действительного значения V, DEC
= chgwdeccfs(DEC,'all',V) замены все коэффициенты V.
DEC = chgwdeccfs(...,IDXSIG) заменяет коэффициенты для сигналов, индексы которых даны векторным IDXSIG. Если исходные данные хранятся построчные или постолбцовые в матричном X, затем IDXSIG содержит строку или индексы столбца, соответственно, данных.
Загрузите 23 канала данные EEG Espiga3 [1]. Каналы располагаются по столбцам. Данные производятся на уровне 200 Гц.
load Espiga3Выполните разложение на уровне 2 с помощью db2 вейвлет.
dec = mdwtdec('c',Espiga3,2,'db2')
dec = struct with fields:
dirDec: 'c'
level: 2
wname: 'db2'
dwtFilters: [1x1 struct]
dwtEXTM: 'sym'
dwtShift: 0
dataSize: [995 23]
ca: [251x23 double]
cd: {[499x23 double] [251x23 double]}
Измените коэффициенты деталей на уровне 1. Замените все значения 0.
decBis = chgwdeccfs(dec,'cd',0,1);Измените коэффициенты деталей на уровне 1 и уровне 2 для сигналов 11 - 15. Замените все значения 0.
decTer = chgwdeccfs(dec,'cd',0,1:2,11:15);Сравните исходные и новые коэффициенты для деталей на уровне 1 для сигналов 11 - 15.
plot(dec.cd{1}(:,11:15),'b')
hold on
plot(decTer.cd{1}(:,11:15),'r')
legend('Original','Changed')
[1] Мезаструктура, Гектор. “Адаптированные Вейвлеты для Обнаружения Шаблона”. Происходящий в Распознавании образов, Анализе изображения и Приложениях, отредактированных Альберто Санфелиу и Мануелем Лазо Кортесом, 3773:933–44. Берлин, Гейдельберг: Спрингер Берлин Гейдельберг, 2005. https://doi.org/10.1007/11578079_96.