Изменение коэффициентов мультисигнала 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-канальных данных ЭЭГ 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. Берлин, Гейдельберг: Springer Berlin Heidelberg, 2005. https://doi.org/10.1007/11578079_96 .