Непрерывный 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-ая строка 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.