StationaryWaveletPacketTransform[data]
gives the stationary wavelet packet transform (SWPT) of an array of data.
StationaryWaveletPacketTransform[data,wave]
gives the stationary wavelet packet transform using the wavelet wave.
StationaryWaveletPacketTransform[data,wave,r]
gives the stationary wavelet packet transform using r levels of refinement.
StationaryWaveletPacketTransform
StationaryWaveletPacketTransform[data]
gives the stationary wavelet packet transform (SWPT) of an array of data.
StationaryWaveletPacketTransform[data,wave]
gives the stationary wavelet packet transform using the wavelet wave.
StationaryWaveletPacketTransform[data,wave,r]
gives the stationary wavelet packet transform using r levels of refinement.
Details and Options
- StationaryWaveletPacketTransform is a generalization of StationaryWaveletTransform where the full tree of wavelet coefficients is computed.
- StationaryWaveletPacketTransform gives a DiscreteWaveletData object.
- Properties of the DiscreteWaveletData dwd can be found using dwd["prop"], and a list of available properties can be found using dwd["Properties"].
- The resulting wavelet coefficients are arrays of the same depth and dimensions as the input data.
- The data can be any of the following:
-
list arbitrary-rank numerical array image arbitrary Image object audio an Audio or sampled Sound object - The possible wavelets wave include:
-
BattleLemarieWavelet[…] Battle–Lemarié wavelets based on B-spline BiorthogonalSplineWavelet[…] B-spline-based wavelet CoifletWavelet[…] symmetric variant of Daubechies wavelets DaubechiesWavelet[…] the Daubechies wavelets HaarWavelet[…] classic Haar wavelet MeyerWavelet[…] wavelet defined in the frequency domain ReverseBiorthogonalSplineWavelet[…] B-spline-based wavelet (reverse dual and primal) ShannonWavelet[…] sinc function-based wavelet SymletWavelet[…] least asymmetric orthogonal wavelet - The default wave is HaarWavelet[].
- With higher settings for the refinement level r, larger scale features are resolved.
- The default refinement level r is given by
, where
is the minimum dimension of data. - With refinement level Full, r is given by
. - The tree of wavelet coefficients at level
consists of coarse coefficients
and detail coefficients
, with
representing the input data. - The forward transform is given by
,
,
, and
, where
is the filter length for the corresponding wspec and
is the length of input data. - The inverse transform is given by
. - The
are lowpass filter coefficients and
are highpass filter coefficients that are defined for each wavelet family. - The following options can be given:
-
Method Automatic method to use WorkingPrecision MachinePrecision precision to use in internal computations - StationaryWaveletPacketTransform uses periodic padding of data.
- InverseWaveletTransform gives the inverse transform.
- By default, InverseWaveletTransform uses coefficients represented by dwd["BasisIndex"] for reconstruction. Use WaveletBestBasis to compute and set an optimal basis.
Examples
open all close allBasic Examples (3)
Compute a stationary wavelet packet transform:
dwd = StationaryWaveletPacketTransform[{0, 0, 1, 0, 0}]The resulting DiscreteWaveletData represents a full tree of wavelet coefficients:
dwd["TreeView"]The inverse transform reconstructs the input:
InverseWaveletTransform[dwd]a = Audio["ExampleData/rule30.wav"]dwd = StationaryWaveletPacketTransform[a, Automatic, 2]Use dwd[…,"Audio"] to extract coefficient signals:
dwd[All, "Audio"]AudioLength[#[[2]]]& /@ %Compute the inverse transform:
InverseWaveletTransform[dwd]Transform an Image object:
dwd = StationaryWaveletPacketTransform[[image], Automatic, 1]Use dwd[…,"Image"] to extract coefficient images:
dwd[All, "Image"]Compute the inverse transform:
InverseWaveletTransform[dwd]Scope (33)
Basic Uses (4)
Useful properties can be extracted from the DiscreteWaveletData object:
dwd = StationaryWaveletPacketTransform[RandomReal[1, {16}], DaubechiesWavelet[4], 2]Get a full list of properties:
dwd["Properties"]Get data and coefficient dimensions:
dwd["DataDimensions"]dwd["Dimensions"]Use Normal to get all wavelet coefficients explicitly:
dwd = StationaryWaveletPacketTransform[Range[5]];Normal[dwd]Also use All as an argument to get all coefficients:
dwd[All]Use Automatic to get only the coefficients used in the inverse transform:
dwd[Automatic]Use the "TreeView" or "WaveletIndex" to find out what wavelet coefficients are available:
dwd = StationaryWaveletPacketTransform[Range[5]];dwd["TreeView"]dwd["WaveletIndex"]Extract specific coefficient arrays:
dwd[{0}]dwd[{0, 0}]Extract several wavelet coefficients corresponding to a list of wavelet index specifications:
dwd[{{0}, {0, 1}}]Extract all coefficients whose wavelet indexes match a pattern:
dwd[{_}]dwd[{{_, 0}, {_, 1}}]Use a higher refinement level to increase the frequency resolution:
data = Table[Sin[x ^ 2] + RandomReal[{-0.2, 0.2}], {x, 0, 10, 0.02}];ListLinePlot[data]With a smaller refinement level, more of the signal energy is left in {0,0}:
dwd1 = StationaryWaveletPacketTransform[data, DaubechiesWavelet[4], 2];WaveletListPlot[dwd1, PlotLayout -> "CommonYAxis", Ticks -> Full]With further refinement, {0,0} is resolved into further components:
dwd2 = StationaryWaveletPacketTransform[data, DaubechiesWavelet[4], 3];WaveletListPlot[dwd2, PlotLayout -> "CommonYAxis", Ticks -> Full]Wavelet Families (10)
Compute the wavelet packet transform using different wavelet families:
data = Table[Sin[x ^ 2], {x, 0, 10, 0.02}];dwd1 = StationaryWaveletPacketTransform[data, DaubechiesWavelet[4], 3];
dwd2 = StationaryWaveletPacketTransform[data, SymletWavelet[4], 3];{WaveletListPlot[dwd1, PlotLayout -> "CommonXAxis"], WaveletListPlot[dwd2, PlotLayout -> "CommonXAxis"]}Use different families of wavelets to capture different features:
data = Table[4 Exp[-100 (t - 0.5) ^ 2] + Sin[5Pi t], {t, 0, 1, 0.01}];ListLinePlot[data]HaarWavelet (default):
dwd = StationaryWaveletPacketTransform[data, HaarWavelet[], 3];WaveletListPlot[dwd, All, PlotLayout -> "CommonYAxis", Filling -> Axis]data = Table[4 Exp[-100 (t - 0.5) ^ 2] + Sin[5Pi t], {t, 0, 1, 0.01}];dwd = StationaryWaveletPacketTransform[data, DaubechiesWavelet[2], 3];WaveletListPlot[dwd, All, PlotLayout -> "CommonYAxis", Filling -> Axis]data = Table[4 Exp[-100 (t - 0.5) ^ 2] + Sin[5Pi t], {t, 0, 1, 0.01}];dwd = StationaryWaveletPacketTransform[data, BattleLemarieWavelet[3], 3];WaveletListPlot[dwd, All, PlotLayout -> "CommonYAxis", Filling -> Axis]data = Table[4 Exp[-100 (t - 0.5) ^ 2] + Sin[5Pi t], {t, 0, 1, 0.01}];dwd = StationaryWaveletPacketTransform[data, BiorthogonalSplineWavelet[4, 2], 3];WaveletListPlot[dwd, All, PlotLayout -> "CommonYAxis", Filling -> Axis]data = Table[4 Exp[-100 (t - 0.5) ^ 2] + Sin[5Pi t], {t, 0, 1, 0.01}];dwd = StationaryWaveletPacketTransform[data, CoifletWavelet[2], 3];WaveletListPlot[dwd, All, PlotLayout -> "CommonYAxis", Filling -> Axis]data = Table[4 Exp[-100 (t - 0.5) ^ 2] + Sin[5Pi t], {t, 0, 1, 0.01}];dwd = StationaryWaveletPacketTransform[data, MeyerWavelet[3], 3];WaveletListPlot[dwd, All, PlotLayout -> "CommonYAxis", Filling -> Axis]ReverseBiorthogonalSplineWavelet:
data = Table[4 Exp[-100 (t - 0.5) ^ 2] + Sin[5Pi t], {t, 0, 1, 0.01}];dwd = StationaryWaveletPacketTransform[data, ReverseBiorthogonalSplineWavelet[4, 2], 3];WaveletListPlot[dwd, All, PlotLayout -> "CommonYAxis", Filling -> Axis]data = Table[4 Exp[-100 (t - 0.5) ^ 2] + Sin[5Pi t], {t, 0, 1, 0.01}];dwd = StationaryWaveletPacketTransform[data, ShannonWavelet[8], 3];WaveletListPlot[dwd, All, PlotLayout -> "CommonYAxis", Filling -> Axis]data = Table[4 Exp[-100 (t - 0.5) ^ 2] + Sin[5Pi t], {t, 0, 1, 0.01}];dwd = StationaryWaveletPacketTransform[data, SymletWavelet[3], 3];WaveletListPlot[dwd, All, PlotLayout -> "CommonYAxis", Filling -> Axis]1-Dimensional Data (6)
Plot the coefficients over a common horizontal axis using WaveletListPlot:
dwd = StationaryWaveletPacketTransform[Table[Sin[x ^ 2], {x, 0, 10, 0.02}], Automatic, 3];WaveletListPlot[dwd]Plot against a common vertical axis:
WaveletListPlot[dwd, PlotLayout -> "CommonYAxis"]Visualize coefficients as a function of time and refinement level using WaveletScalogram:
dwd = StationaryWaveletPacketTransform[Table[Sin[20x] + Sin[10x], {x, 1, 10, 0.01}]];The coefficient indexes appear as tooltips when the mouse pointer is moved over a coefficient:
WaveletScalogram[dwd, All]data = Table[1, {x, 0, 10, 0.02}];ListLinePlot[data]All coefficients are small except coarse coefficients {0,0,…}:
dwd = StationaryWaveletPacketTransform[data, DaubechiesWavelet[4], 3];WaveletListPlot[dwd, All, PlotLayout -> "CommonYAxis", ImageSize -> Medium, Filling -> Axis, Ticks -> Full]Data oscillating at the highest resolvable frequency (Nyquist frequency):
data = Table[(-1)^n, {n, 80}];ListLinePlot[data]Only the first detail coefficient {1} and its coarse child coefficients {1,0,0,…} are not small:
dwd = StationaryWaveletPacketTransform[data, Automatic, 3];WaveletListPlot[dwd, All, PlotLayout -> "CommonYAxis", Ticks -> Full]Data with large discontinuities:
data = Table[Piecewise[{{1, 1 < x < 3.1}}, 0], {x, 0, 4, 0.02}];ListLinePlot[data]Coarse coefficients {0,…} have the same large-scale structure as the data:
dwd = StationaryWaveletPacketTransform[data, Automatic, 3];WaveletListPlot[dwd, {0..}, FrameTicks -> Full]Detail coefficients are sensitive to discontinuities:
WaveletListPlot[dwd, Except[{0..}]]Data with both spatial and frequency structure:
data = Table[Piecewise[{{5, n <= 30}, {2 + (-1)^n, Inequality[30, Less, n, LessEqual, 60]}, {5, 60 < n}}], {n, 90}];ListLinePlot[data, PlotRange -> {0, 5}]Coarse coefficients {0,…} track the local mean of the data:
dwd = StationaryWaveletPacketTransform[data, Automatic, 3];WaveletListPlot[dwd, {0..}, Ticks -> Full]First detail coefficient {1} and its coarse child coefficients {1,0,…} represent the oscillations:
WaveletListPlot[dwd, {1, 0...}, Ticks -> Full]All coefficients on a common vertical axis:
WaveletListPlot[dwd, All, PlotLayout -> "CommonYAxis", Ticks -> Full]2-Dimensional Data (5)
Compute a two-dimensional stationary wavelet packet transform:
dwd = StationaryWaveletPacketTransform[(| | | | | |
| - | - | - | - | - |
| 0 | 0 | 1 | 0 | 0 |
| 0 | 1 | 1 | 1 | 0 |
| 1 | 1 | 1 | 1 | 1 |
| 0 | 1 | 1 | 1 | 0 |
| 0 | 0 | 1 | 0 | 0 |)]View the tree of wavelet coefficients:
dwd[{"TreeView", Center}]Inverse transform to get back the original signal:
InverseWaveletTransform[dwd]//Chop//MatrixFormUse dwd[…,"MatrixPlot"] to visualize each coefficient as a MatrixPlot:
dwd = StationaryWaveletPacketTransform[DiamondMatrix[32], Automatic, 1]dwd[Automatic, "MatrixPlot"]Visualize diagonal detail coefficient {3} and its child coefficients {3,__}:
dwd2 = StationaryWaveletPacketTransform[DiamondMatrix[32], Automatic, 2];dwd2[{{3}, {3, __}}, "MatrixPlot"]In two dimensions, the vector of filtering operations in each direction can be computed:
Tuples[{0, 1}, 2] /. {0 -> "lowpass", 1 -> "highpass"}Interpreting these vectors as binary digit expansions, you get wavelet index numbers:
FromDigits[#, 2]& /@ Tuples[{0, 1}, 2]Get the lowpass and highpass filters for a Haar wavelet:
{lp, hp} = Map[Last, WaveletFilterCoefficients[HaarWavelet[], {"PrimalLowpass", "PrimalHighpass"}], {2}]The resulting 2D filters are outer products of filters in the two directions:
Apply[KroneckerProduct, Tuples[{lp, hp}, 2], {1}]Map[MatrixPlot, %]Wavelet transform of step data:
step[{a_, b_, c_, d_}] :=
ArrayFlatten[{{ConstantArray[a, {3, 3}], ConstantArray[b, {3, 5}]}, {ConstantArray[c, {5, 3}], ConstantArray[d, {5, 5}]}}];Data with a vertical discontinuity:
MatrixPlot[step[{0, 1, 0, 1}], FrameTicks -> False]All horizontal and diagonal detail coefficients, wavelet index {___,2|3,___}, are zero:
StationaryWaveletPacketTransform[step[{0, 1, 0, 1}], Automatic, 2][Automatic, {"MatrixPlot", ImageSize -> 40}]Data with horizontal discontinuity:
MatrixPlot[step[{0, 0, 1, 1}], FrameTicks -> False]All vertical and diagonal detail coefficients, wavelet index {___,1|3,___}, are zero:
StationaryWaveletPacketTransform[step[{0, 0, 1, 1}], Automatic, 2][Automatic, {"MatrixPlot", ImageSize -> 40}]Higher-Dimensional Data (2)
Compute a three-dimensional wavelet packet transform:
data = RandomReal[1, {16, 16, 16}];dwd = StationaryWaveletPacketTransform[data, Automatic, 2]List all computed wavelet coefficients:
dwd["WaveletIndex"]Inverse transform to get back the original signal:
Norm[Flatten[data - InverseWaveletTransform[dwd]]]Wavelet transform of a three-dimensional cross array:
data = CrossMatrix[All, {8, 8, 8}];Graphics3D[{Red, Cuboid /@ Position[data, 1]}]dwd = StationaryWaveletPacketTransform[data, Automatic, 2]Visualize lowpass wavelet coefficients {___,0}:
Table[i -> Graphics3D[{If[Positive[Extract[dwd[i][[1, 2]], #]], Red, Green], Cuboid[#]}& /@ Position[dwd[i][[1, 2]], u_ /; Abs[u] > 0], PlotRange -> {{1, 9}, {1, 9}, {1, 9}}], {i, Cases[dwd["WaveletIndex"], {___, 0}]}]Energy of the original data is conserved within the transformed coefficients:
Total[Flatten[data]^2] == Total[Flatten[dwd[Automatic][[All, 2]]]^2]Audio Data (2)
Transform an Audio object:
a = ExampleData[{"Audio", "Apollo11ReturnSafely"}];dwd = StationaryWaveletPacketTransform[a]The inverse transform yields a reconstructed audio:
InverseWaveletTransform[dwd]By default, coefficients are given as lists of data for each sound channel:
dwd = StationaryWaveletPacketTransform[ExampleData[{"Audio", "PianoScale"}]];Dimensions[{1, 1} /. dwd[{1, 1}]]Get the {1,1} coefficient as an Audio object:
dwd[{1, 1}, "Audio"]Inverse transform of {1,1} coefficient as an Audio object:
dwd[{1, 1}, {"Audio", "Inverse"}]Sound Data (2)
Transform a Sound object:
snd = ExampleData[{"Sound", "Apollo11ReturnSafely"}]dwd = StationaryWaveletPacketTransform[snd]The inverse transform yields a reconstructed audio object:
InverseWaveletTransform[dwd]Browse all coefficients using a MenuView:
dwd = StationaryWaveletPacketTransform[ExampleData[{"Sound", "Clarinet"}]];MenuView[dwd[All, "Sound"]]Image Data (2)
Transform an Image object:
img = Image[DiamondMatrix[All, {64, 64}]]dwd = StationaryWaveletPacketTransform[img, HaarWavelet[], 3]The inverse transform yields a reconstructed Image object:
InverseWaveletTransform[dwd]Wavelet coefficients are normally given as lists of data for each image channel:
dwd = StationaryWaveletPacketTransform[[image], HaarWavelet[], 2];Dimensions[{1, 1} /. dwd[{1, 1}]]Get all coefficients as Image objects instead:
dwd[All, "Image"]//ShortGet raw Image objects with no rescaling of color levels:
dwd[All, {"Image", "ImageFunction" -> Identity}]//ShortGet the inverse transform of the {0,1} coefficient as an Image object:
dwd[{0, 1}, {"Image", "Inverse"}]Generalizations & Extensions (3)
StationaryWaveletPacketTransform works on arrays of symbolic quantities:
dwd = StationaryWaveletPacketTransform[{a, b, c, d}, WorkingPrecision -> ∞];dwd[Automatic]//SimplifyInverse transform recovers the input exactly:
InverseWaveletTransform[dwd]//SimplifySpecify any internal working precision:
dwd = StationaryWaveletPacketTransform[{1, 2, 5, 2}, WorkingPrecision -> 20];dwd[Automatic]data = Exp[I RandomReal[1, 4]];dwd = StationaryWaveletPacketTransform[data]The wavelets coefficients are complex:
dwd[Automatic]Options (3)
WorkingPrecision (3)
By default, WorkingPrecision->MachinePrecision is used:
data = RandomInteger[1, {10}];dwd1 = StationaryWaveletPacketTransform[data]dwd2 = StationaryWaveletPacketTransform[data, WorkingPrecision -> MachinePrecision]dwd1 == dwd2Use higher-precision computation:
data = {0, 0, 1, 1};dwd = Normal@StationaryWaveletPacketTransform[data, Automatic, 1, WorkingPrecision -> 20]{Precision[dwd], Accuracy[dwd]}Use WorkingPrecision->∞ for exact computation:
data = RandomInteger[10, {4}];Normal@StationaryWaveletPacketTransform[data, WorkingPrecision -> ∞]//SimplifyProperties & Relations (10)
StationaryWaveletPacketTransform computes the full tree of wavelet coefficients:
swpt = StationaryWaveletPacketTransform[{1, 1, 3, 1, 1}];swpt[{"TreeView", Left}]StationaryWaveletTransform computes a subset of the full tree of coefficients:
swt = StationaryWaveletTransform[{1, 1, 3, 1, 1}];swt[{"TreeView", Left}]DiscreteWaveletPacketTransform coefficients halve in length with each level of refinement:
Normal[DiscreteWaveletPacketTransform[{1, 2, 3, 4}]]Rotated data gives different coefficients:
Normal[DiscreteWaveletPacketTransform[{2, 3, 4, 1}]]StationaryWaveletPacketTransform coefficients have the same length as the data:
Normal[StationaryWaveletPacketTransform[{1, 2, 3, 4}]]Rotated data gives rotated coefficients:
Normal[StationaryWaveletPacketTransform[{2, 3, 4, 1}]]The default refinement is given by Min[Round[Log2[Min[Dimensions[data]]]],4]:
data = RandomReal[1, {100}];r = Min[Round[Log2[Min[Dimensions[data]]]], 4]StationaryWaveletPacketTransform[data] == StationaryWaveletPacketTransform[data, Automatic, r]data = RandomReal[1, {100, 10, 10}];r = Min[Round[Log2[Min[Dimensions[data]]]], 4]StationaryWaveletPacketTransform[data] == StationaryWaveletPacketTransform[data, Automatic, r]The energy norm is conserved for orthogonal wavelet families:
data = RandomReal[1, {100}];swt = StationaryWaveletPacketTransform[data];Norm[data] == Norm[Flatten[Last /@ swt[Automatic]]]The energy norm is approximately conserved for biorthogonal wavelet families:
data = RandomReal[1, {100}];swt = StationaryWaveletPacketTransform[data, BiorthogonalSplineWavelet[2, 4]];Norm[data]Norm[Flatten[Last /@ swt[Automatic]]]The mean of the data is captured at the maximum refinement level of the transform:
data = RandomReal[1, 16];swt = StationaryWaveletPacketTransform[data, HaarWavelet[]];Extract the coefficient for the maximum refinement level:
r = swt["Refinement"]swt[ConstantArray[0, {r}]]Mean[data]The sum of inverse transforms from individual coefficient arrays gives the original data:
data = Table[DiscreteDelta[n], {n, -2, 2}]dwd = StationaryWaveletPacketTransform[data];dwd["TreeView"]Individually inverse transform each wavelet coefficient array:
data1 = InverseWaveletTransform[dwd, Automatic, {0, 0}]data2 = InverseWaveletTransform[dwd, Automatic, {0, 1}]data3 = InverseWaveletTransform[dwd, Automatic, {1, 0}]data4 = InverseWaveletTransform[dwd, Automatic, {1, 1}]The sum gives the original data:
data1 + data2 + data3 + data4Compute a Haar stationary wavelet packet transform in one dimension:
low[v_] := ListCorrelate[{1 / 2, 1 / 2}, v, -1]
high[v_] := ListCorrelate[{-1 / 2, 1 / 2}, v, -1]HaarWaveletTransform[v_] := {{0} -> low[v], {1} -> high[v]};Compute {0} and {1} wavelet coefficients:
HaarWaveletTransform[{a, b, c, d}]Compare with StationaryWaveletPacketTransform:
StationaryWaveletPacketTransform[{a, b, c, d}, HaarWavelet[], 1, WorkingPrecision -> ∞][All]In two dimensions, a separate filter is applied in each dimension:
f2d[fx_, fy_] := Composition[Map[fy, #]&, Map[fx, #]&]Lowpass and highpass filters for Haar wavelet:
low[v_] := ListCorrelate[{1 / 2, 1 / 2}, v, -1]
high[v_] := ListCorrelate[{-1 / 2, 1 / 2}, v, -1]Haar wavelet transform of matrix data:
data = Table[Sin[x y], {x, -2, 2, 4 / 63}, {y, -2, 2, 4 / 63}];Table[MatrixPlot[f2d[fx, fy][data], PlotLabel -> {fx, fy}, FrameTicks -> None], {fx, {low, high}}, {fy, {low, high}}]//FlattenCompare with StationaryWaveletPacketTransform using HaarWavelet:
dwd = StationaryWaveletPacketTransform[data, HaarWavelet[], 1];Table[MatrixPlot[Last[p], PlotLabel -> First[p], FrameTicks -> None], {p, Normal[dwd]}]Image channels are transformed individually:
dwds = Map[StationaryWaveletPacketTransform, ColorSeparate[[image]]]Combine {0} coefficients of separately transformed image channels:
w = Image[Table[First[{0} /. Normal[t]], {t, dwds}], Interleaving -> False]Compare with {0} coefficient of StationaryWaveletPacketTransform of original image:
dwd = StationaryWaveletPacketTransform[[image]];{0} /. dwd[All, {"Image", "ImageFunction" -> Identity}]ImageSubtract[w, %]Related Guides
Text
Wolfram Research (2010), StationaryWaveletPacketTransform, Wolfram Language function, https://reference.wolfram.com/language/ref/StationaryWaveletPacketTransform.html (updated 2017).
CMS
Wolfram Language. 2010. "StationaryWaveletPacketTransform." Wolfram Language & System Documentation Center. Wolfram Research. Last Modified 2017. https://reference.wolfram.com/language/ref/StationaryWaveletPacketTransform.html.
APA
Wolfram Language. (2010). StationaryWaveletPacketTransform. Wolfram Language & System Documentation Center. Retrieved from https://reference.wolfram.com/language/ref/StationaryWaveletPacketTransform.html
BibTeX
@misc{reference.wolfram_2026_stationarywaveletpackettransform, author="Wolfram Research", title="{StationaryWaveletPacketTransform}", year="2017", howpublished="\url{https://reference.wolfram.com/language/ref/StationaryWaveletPacketTransform.html}", note=[Accessed: 13-June-2026]}
BibLaTeX
@online{reference.wolfram_2026_stationarywaveletpackettransform, organization={Wolfram Research}, title={StationaryWaveletPacketTransform}, year={2017}, url={https://reference.wolfram.com/language/ref/StationaryWaveletPacketTransform.html}, note=[Accessed: 13-June-2026]}