Непрерывное 1-D вейвлет-преобразование
Эта версия 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 - вектор 1-D с положительными элементами. Вектор символов или строковый скаляр 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,' возвращает и строит график скалограммы. wname','scal')'scal' создает график изображения скалограммы.
[coefs,sgram] = cwt(x,scales,' отображает контурное представление скалограммы. wname','scalCNT')
[coefs,frequencies] = cwt(x,scales, возвращает частоты в циклах за единицу времени, соответствующие масштабам и анализирующему вейвлету wname, samplingperiod)wname. samplingperiod - положительный действительный скаляр. Если единицы измерения samplingperiod секунды, частоты в герцах.
[coefs,sgram,frequencies] = cwt(x,scales, возвращает скалограмму и частоты, соответствующие масштабам и анализирующему вейвлету. При наличии хотя бы двух элементов в wname, samplingperiod,'scal')scales, также можно использовать флаг 'scalCNT' для вывода скалограммы. samplingperiod используется только при преобразовании шкал в частоты. Определение samplingperiod не влияет на внешний вид графиков, созданных cwt.
Постройте график непрерывного вейвлет-преобразования и скалограммы с использованием сим2-вейвлета при всех целых шкалах от 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, I. Десять лекций по вейвлетам, Филадельфия, Пенсильвания: Общество промышленной и прикладной математики (SIAM), 1992.
Маллат, S. Wavelet Tour of Signal Processing, Сан-Диего, Калифорния: Академическая пресса, 1998.