cwt
Непрерывный 1D вейвлет преобразовывает
Эта версия cwt больше не рекомендуется. Используйте обновленный cwt вместо этого.
Синтаксис
coefs = cwt(x,scales,'wname')
coefs = cwt(x,scales,'wname','plot')
coefs = cwt(x,scales,'wname','coloration')
coefs = cwt(x,scales,'wname','coloration',xlim)
[coefs,sgram] = cwt(x,scales,'wname','scal')
[coefs,sgram] = cwt(x,scales,'wname','scalCNT')
[coefs,frequencies] = cwt(x,scales,wname,
samplingperiod)
[coefs,sgram,frequencies] = cwt(x,scales,wname,
samplingperiod,'scal')
Описание
coefs = cwt(x,scales,' возвращает непрерывный вейвлет преобразовывает (CWT) wname')x сигнала с действительным знаком. Вейвлет преобразовывает, вычисляется для заданных шкал с помощью вейвлета анализа wname. scales является 1D вектором с положительными элементами. Скаляр вектора символов или строки wname обозначает вейвлет, распознанный wavemngr. coefs является матрицей с количеством строк, равных длине scales и количеству столбцов, равных длине входного сигнала. k-th строка coefs соответствует коэффициентам CWT для k-й элемента в векторе шкал.
coefs = cwt(x,scales,' графики непрерывный вейвлет преобразовывают коэффициенты, с помощью wname','plot')coloration по умолчанию 'absglb'.
coefs = cwt(x,scales,' использует заданную окраску. См. Определения для опций окраски.wname','coloration')
coefs = cwt(x,scales,' окрашивает коэффициенты с помощью окраски и wname','coloration',xlim)xlim, где xlim является вектором, [x1 x2], с 1 ≤ x1 <x2 ≤ length(x).
[coefs,sgram] = cwt(x,scales,' возвращает и строит scalogram. wname','scal')'scal' производит график изображений scalogram.
[coefs,sgram] = cwt(x,scales,' отображает представление контура scalogram. wname','scalCNT')
[coefs,frequencies] = cwt(x,scales, возвращает частоты в циклах в единицу времени, соответствуя шкалам и вейвлету анализа wname,
samplingperiod)wname. samplingperiod является положительным скаляром с действительным знаком. Если модули samplingperiod являются секундами, частоты находятся в герц.
[coefs,sgram,frequencies] = cwt(x,scales, возвращает scalogram и частоты, соответствующие шкалам и вейвлету анализа. Если у вас есть по крайней мере два элемента в wname,
samplingperiod,'scal')scales, можно также использовать флаг 'scalCNT' , чтобы вывести scalogram. samplingperiod используется только в преобразовании шкал к частотам. Определение samplingperiod не влияет на внешний вид графиков, сгенерированных cwt.
Примеры
График непрерывный вейвлет преобразовывает и scalogram, использующий sym2 вейвлет во всех целочисленных шкалах от 1 до 32, с помощью фрактального сигнала, как введено:
load vonkoch
vonkoch=vonkoch(1:510);
len = length(vonkoch);
cw1 = cwt(vonkoch,1:32,'sym2','plot');
title('Continuous Transform, absolute coefficients.')
ylabel('Scale')
[cw1,sc] = cwt(vonkoch,1:32,'sym2','scal');
title('Scalogram')
ylabel('Scale')
Сравните дискретный и непрерывный вейвлет преобразовывает, с помощью фрактального сигнала, как введено:
load vonkoch
vonkoch=vonkoch(1:510);
len=length(vonkoch);
[c,l]=wavedec(vonkoch,5,'sym2');
% Compute and reshape DWT to compare with CWT.
cfd=zeros(5,len);
for k=1:5
d=detcoef(c,l,k);
d=d(ones(1,2^k),:);
cfd(k,:)=wkeep(d(:)',len);
end
cfd=cfd(:);
I=find(abs(cfd) <sqrt(eps));
cfd(I)=zeros(size(I));
cfd=reshape(cfd,5,len);
% Plot DWT.
subplot(311); plot(vonkoch); title('Analyzed signal.');
set(gca,'xlim',[0 510]);
subplot(312);
image(flipud(wcodemat(cfd,255,'row')));
colormap(pink(255));
set(gca,'yticklabel',[]);
title('Discrete Transform,absolute coefficients');
ylabel('Level');
% Compute CWT and compare with DWT
subplot(313);
ccfs=cwt(vonkoch,1:32,'sym2','plot');
title('Continuous Transform, absolute coefficients');
set(gca,'yticklabel',[]);
ylabel('Scale');
Больше о
Ссылки
Daubechies, я. Десять лекций по вейвлетам, Филадельфии, PA: общество промышленной и прикладной математики (SIAM), 1992.
Mallat, S. Тур вейвлета по обработке сигналов, Сан-Диего, CA: Academic Press, 1998.