Пропускная способность набора фильтров CWT 3 дБ
bw = powerbw(fb)bw = powerbw(fb)fb. bw является Ns-by-4 MATLAB®
table, где Ns является количеством частот полосы пропускания вейвлета (равный количеству шкал). Для каждого фильтра в fb таблица содержит соответствующую полосовую частоту, пропускную способность на 3 дБ, и более низкую частоту и верхние пределы частоты пропускной способности на 3 дБ.
Пропускная способность на 3 дБ ограничивает метку, где степень фильтра является половиной ее пикового значения. Частотная характеристика значения в пределах равна 1 / √ 2 раза пиковое значение. Поскольку полосы пропускания в fb нормированы с пиковыми значениями, приблизительно равняются 2, частотная характеристика значения в каждом пределе приблизительно равна 2 / √ 2. Пропускная способность на 3 дБ также известна как пропускную способность на уровне половинной мощности потому что 20 log101/√2 ≅ -3.
Создайте набор фильтров CWT.
fb = cwtfilterbank;
Получите пропускную способность (на уровне половинной мощности) на 3 дБ набора фильтров. Получите частотные характеристики вейвлетов.
bw = powerbw(fb); [psidft,f] = freqz(fb);
Осмотрите первые несколько строк в таблице bw.
bw(1:10,:)
ans=10×4 table
Frequencies HalfPowerBandwidth LowFrequencyBorder HighFrequencyBorder
___________ __________________ __________________ ___________________
0.43412 0.09323 0.38751 0.48074
0.40505 0.086988 0.36156 0.44855
0.37793 0.081162 0.33735 0.41851
0.35262 0.075726 0.31476 0.39048
0.329 0.070655 0.29368 0.36434
0.30697 0.065924 0.27401 0.33994
0.28642 0.061509 0.25566 0.31717
0.26723 0.05739 0.23854 0.29593
0.24934 0.053547 0.22257 0.27611
0.23264 0.049961 0.20766 0.25762
Выберите полосовой фильтр вейвлета из набора фильтров. Выписка из таблицы bw пределы на 3 дБ полосового фильтра.
wv = 5; frq = bw.Frequencies(wv); lfb = bw.LowFrequencyBorder(wv); hfb = bw.HighFrequencyBorder(wv);
Постройте частотную характеристику и пределы на 3 дБ. Поскольку частотная характеристика масштабируется, чтобы иметь максимальное значение, равное 2, осмотреть график подтвердить, что более низкий и верхний крест границ частоты пересекает частотную характеристику в sqrt(2).
plot(f,psidft(wv,:)) grid on hold on plot([lfb lfb],[0 2],'r') plot([hfb hfb],[0 2],'r') xlabel('Normalized Frequency (cycles/sample)') ylabel('Magnitude') title(['Bandpass Frequency: ' num2str(frq) ' cycles/sample'])
