Шкалы набора фильтров 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')
