Види та порядок проведення вейвлет-аналізу
1. Обчислення багатовіконне перетворення
Нехай x(t) – сигнал, який необхідно аналізувати. Вибирається материнський вейвлет, що буде прототипом для всіх функцій (вікон), які можна отримати з нього шляхом стиснення (розширення). Існує кілька функцій, що застосовуються як материнські вейвлети. Двома прикладами є вейвлети Морле та Мексиканський капелюх, які й використовуються у прикладах розділу.
Після вибору материнської функції обчислення починаються з масштабу s=1. БВП обчислюється для всіх значень s, менших і більших 1. Однак повне перетворення зазвичай не потрібно, оскільки реальні сигнали обмежені за смугою. Тому число масштабів може бути обмежено. У прикладах цього розділу ми також використовуємо обмежену кількість масштабів.
Процедура аналізу стартує з масштабу s=1 і триває при значеннях s, що збільшуються, тобто аналіз починається з високих частот і продовжується убік низьких частот. Перше значення s відповідає найбільш стислому вейвлету. При збільшенні значення s вейвлет розширюється. Вейвлет розміщується в початку сигналу, у точці, що відповідає часу =0. Вейвлет-функція масштабу «1» помножується з сигналом та інтегрується на всьому часовому інтервалі. Вейвлет масштабу s=1 потім зміщується вправо на до точки t=, і процедура повторюється. Одержуємо ще одне значення, яке відповідає t=, s=1 на частотно-часовому плані. Ця процедура повторюється доти, поки вейвлет не досягне кінця сигналу. У такий спосіб одержуємо ряд точок на масштабно-часовому плані для масштабу s=1. Тепер збільшимо s на деяке значення. Точно кажучи, оскільки перетворення безперервне, то і s мають змінюватися безперервно. Під час виконання перетворення в комп'ютері ми обчислюємо апроксимацію, збільшуючи обидва параметри на деяке мале значення. Тим самим здійснюємо дискретизацію масштабно-часової площини.
Наведена вище процедура повторюється для кожного значення s. При цьому рядок за рядком заповнюється масштабно-часова площина. Так обчислюється БВП. Нижче наведені рисунки ілюструють процес перетворення крок за кроком.
На рис. 1 показаний сигнал і вейвлет-функція для чотирьох різних значень . Значення масштабу, 1 , відповідає найменшому значенню, або найбільшій частоті. Відзначте, яким компактним є носій. Він має бути таким само вузьким, як і час життя найвищої частоти сигналу. На рисунку показані чотири різних позиції вейвлет-функції в точках to=2 , to=40, to=90 та to=140. У кожній позиції вона перемножується з сигналом. Добуток буде ненульовим лише у випадку, коли сигнал перетинається з носієм вейвлета, і нульовим – в інших випадках. Зміщенням вейвлета у часі відповідає часова локалізація сигналу, а зміщенням у масштабі – масштабна (частотна) локалізація.
Рисунок 1 – Сигнал і вейвлет-функція для чотирьох різних значень для масштабу
Якщо в сигналі присутні спектральні компоненти, що відповідають поточному значенню s (яке в цьому випадку 1), то добуток вейвлета з сигналом в інтервалі, де цей спектральний компонент присутній, дає відносно велике значення. У протилежному випадку – добуток малий або дорівнює нулю. Сигнал, показаний на рис. 1, має спектральні компоненти, порівняні із шириною вікна при s=1 на інтервалі біля t=100мс. БВП сигналу, показаного на рис. 1, дає більші значення для низьких масштабів біля часу t=100мс і малі – в інших інтервалах часу. ............
