Вейвлет-преобразование
Программирование других вейвлет-преобразований
Наряду со встроенной функцией wave, профессиональные версии Mathcad 11-12 снабжены пакетом расширения для осуществления вейвлет-анализа. Пакет расширения содержит большое число дополнительных встроенных функций, имеющих отношение к вейвлет-преобразованиям. Обзор пакетов расширения выходит за рамки данной книги, поэтому ограничимся простым упоминанием об этой возможности.
Примечание
Дополнительную информацию об использовании данных встроенных функций можно найти в соответствующей.электронной книге, которую можно открыть при помощи меню Help › E-Books › Wavelet extension pack (Справка › Электронные книги › Вейвлет-анализ данных).
Помимо встроенной функции вейвлет-спектра Добеши и возможностей пакета расширения, допускается непосредственное программирование алгоритмов пользователя для расчета вейвлет-спектров. Оно сводится к аккуратному численному расчету соответствующих семейств интегралов. Один из примеров такой программы приведен в листинге 4.19. Анализу подвергается функция, составленная из суммы двух синусов, а график двухпараметрического спектра с(а,b) выведен на рис. 4.17 в виде привычных для вейвлет-анализа линий уровня на плоскости (а,b).
Примечание
Программа листинга очень проста, но исключительно далека от хорошей в смысле быстродействия. Каждый интеграл вычисляется независимо, без использования методов ускорения, типа применяемых в алгоритме БПФ. Однако простые приемы программирования вполне доступно раскрывают математический смысл вейвлет-преобразования.
Листинг 4.19. Вычисление вейвлет-спектра на основе "мексиканской шляпы":
Рис. 4.17. Вейвлет-спектр модельного сигнала на основе "мексиканской шляпы" (продолжение листинга 4.19)