Шкалы банка фильтров CWT
Создайте банк фильтров CWT с периодом дискретизации 0,001 секунды.
fb = cwtfilterbank('SamplingPeriod',seconds(0.001));Получите необработанные и преобразованные шкалы, используемые при создании вейвлет-полосно-пропускающих фильтров.
[rs,cs] = scales(fb);
Получите периоды полосы пропускания банка фильтров.
P = centerPeriods(fb);
Сравните самую мелкую преобразованную шкалу с наименьшим периодом центра полосы пропускания, нормированным периодом дискретизации.
min(cs)
ans = 2.3035
min(P)/seconds(0.001)
ans = 2.3035
Шкалы должны увеличиться примерно в один раз , где - количество голосов на октаву. Значение по умолчанию равен 10. Постройте график коэффициентов последовательных шкал и сравните с .
len = length(rs); plot(rs(2:len)./rs(1:len-1),'rx-') hold on plot(1:len-1,2^(1/10)*ones(1,len-1),'b') title('Successive Scale Ratios') legend('Scale Ratio','Scale Factor')
