Анализ с помощью вейвлетов представляется в виде дискретного преобразования с помощью умножения матриц. Поэтому неудивительно, я полагаю, что при построении, например, D4, пакет R wmtsa возвращает график:
require(wmtsa)
filters <- wavDaubechies("d4")
plot(filters)
Вопрос в том, как перейти от этого дискретного сюжета к сюжету в записи Википедии :
Обратите внимание, что я не заинтересован в генерации этих кривых именно с wmtsa. Подойдет любой другой пакет - у меня нет Matlab или Mathematica. Но мне интересно, стоит ли начинать с перевода этого фрагмента кода Mathematica из этой статьи в R вместо использования встроенных функций:
Wave1etTransform.m
c[k-1 := c[k] = Daubechies[4][[k+l]];
phi[l] = (l+Sqrt[3])/2 // N;
phi[2] = (l-Sqrt[3])/2 // N;
phi[xJ; xc=0 II x>=3] : = 0
phi[x-?NumberQ] := phi[x] =
N[Sqrt[2]] Sum[c[k] phi[2x-k],{k,0,3}];