Изменение коэффициентов многосигнального 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] Меса, Гектор. «Адаптированные вейвлеты для обнаружения шаблона». В Прогресс Pattern Recognition, Image Analysis and Applications, под редакцией Альберто Санфелиу и Мануэля Лазо Кортеса, 3773: 933-44. Берлин, Гейдельберг: Спрингер Берлин Гейдельберг, 2005. https://doi.org/10.1007/11578079_96 .