Exponenta Banner
exponenta event banner

pat2cwav

Создайте вейвлет из шаблона

Синтаксис

[PSI,XVAL,NC] = pat2cwav(YPAT,METHOD,POLDEGREE,REGULARITY)

Описание

[PSI,XVAL,NC] = pat2cwav(YPAT,METHOD,POLDEGREE,REGULARITY) вычисляет допустимый вейвлет для CWT (данный XVAL и PSI) адаптированный к шаблону задан векторным YPAT, и нормы равняются 1.

Базовый шаблон x-значений установлен в

xpat = linspace(0,1,length(YPAT))

Постоянный NC таково что NC\psi аппроксимирует YPAT на интервале [0,1] использованием выравнивания методом наименьших квадратов

  • полином степени POLDEGREE когда METHOD равно 'polynomial'

  • проекция на пробеле функций, ортогональных к константам, когда METHOD равно 'orthconst'

REGULARITY параметр задает граничные ограничения в точках 0 и 1. Допустимыми значениями является 'continuous', 'differentiable', и 'none'.

Когда METHOD равно 'polynomial'

  • если REGULARITY равно 'continuous', POLDEGREE должен быть больше или быть равен 3.

  • если REGULARITY равно 'differentiable', POLDEGREE должен быть больше или быть равен 5.

Примеры

Принцип для разработки нового вейвлета для CWT должен аппроксимировать данный шаблон с помощью оптимизации наименьших квадратов при ограничительном продвижении к допустимому вейвлету, которому хорошо удовлетворяют для обнаружения шаблона с помощью непрерывного вейвлета, преобразовывают (см. Misiti и др.).

load ptpssin1; 
plot(X,Y), title('Original Pattern')


[psi,xval,nc] = pat2cwav(Y, 'polynomial',6, 'continuous') ;
plot(X,Y,'-',xval,nc*psi,'--'), 
title('Original Pattern and Adapted Wavelet (dashed line)')

Можно проверять тот psi удовлетворяет определению вейвлета путем отмечания, что оно объединяется, чтобы обнулить и что его норма L2 равна 1.

dx = xval(2)-xval(1);
Mu = sum(psi*dx)
L2norm = sum(abs(psi).^2*dx)

Ссылки

Мизити, M., И. Мизити, Г. Оппенхейм, J.-M. Poggi (2003), “Приложения Les ondelettes et leurs”, Гермес.

Представлено до R2006a

Документация Wavelet Toolbox
Поддержка