container : Array, predefined Array for holding results

inplace   : truefalse, specifies that output should overwrite input

The DCT(A) command computes the discrete cosine transform (DCT) of the Array A and returns the result in an Array with datatype float[8].

The InverseDCT(A) command computes the inverse discrete cosine transform of the Array A and returns the result in an Array with datatype float[8].

Before the code performing the computation runs, A is converted to datatype float[8] if it does not have that datatype already. For this reason, it is most efficient if A has this datatype beforehand.

The discrete cosine transform B of a sample of $N$ elements with initial index 1 is defined by the formula

where $C_k$ is given by

Samples may be of arbitrary length, but when the length is a power of $2$, a faster algorithm is used,

If the container=C option is provided, then the results are put into C and C is returned. With this option, no additional memory is allocated to store the result. The container must be an Array having datatype float[8] and size equal to those of A.

If the inplace or inplace=true option is provided, then A is overwritten with the results. In this case, the container option is ignored. Furthermore, A must have datatype float[8].

The SignalProcessing[DCT] and SignalProcessing[InverseDCT] commands are thread-safe as of Maple 17.

For more information on thread safety, see index/threadsafe.

$\mathrm{with}\left(\mathrm{SignalProcessing}\right)\:$

$\mathrm{with}\left(\mathrm{plots}\right)\:$

$N≔10\:$

$A≔\mathrm{GenerateTone}\left(N\,1\,0.1\,0\right)\+\mathrm{GenerateTone}\left(N\,3\,0.4\,0.2\right)\:$

$B≔\mathrm{DCT}\left(A\right)$

$\left[\begin{array}{cccccccccc}-5.617333549722722{10}^{-16}& 1.8408150749374697& 2.1266270208800986& 1.07340522632521& -5.998696806194468{10}^{-16}& 1.8674493435445694& -1.233922005235723{10}^{-15}& 4.5320177608978245& 3.29911475027503& -2.454780728041238\end{array}\right]$

$\mathrm{display}\left(⟨⟨⟨\mathrm{listplot}\left(A\,\'\mathrm{title}\'\=''Signal''\right)⟩\,⟨\mathrm{listplot}\left(B\,\'\mathrm{title}\'\=''Cosine\; Transform''\right)⟩⟩⟩\right)$

$C≔\mathrm{InverseDCT}\left(B\right)$

$\left[\begin{array}{cccccccccc}3.940199733523725& -1.9199792650748926& 1.78442596381239& 0.03271740521865585& -2.8373638374798786& 1.9401997335237255& -3.5380132538247864& 1.1663919750624943& 0.6507513939685514& -1.2193298487299862\end{array}\right]$

$\mathrm{map}\left(\mathrm{fnormal}\,A-C\right)$

$\left[\begin{array}{cccccccccc}0.& -0.& 0.& 0.& 0.& 0.& -0.& -0.& 0.& 0.\end{array}\right]$

$A≔\mathrm{GenerateJaehne}\left(10\,3\right)$

$\left[\begin{array}{cccccccccc}0.0& 0.46930339512069263& 1.7633557568774194& 2.9630650217854133& 1.7633557568774196& -2.1213203435596424& -1.76335575687742& 2.963065021785413& -1.7633557568774183& 0.46930339512069374\end{array}\right]$

$C≔\mathrm{Array}\left(1..10\,\'\mathrm{datatype}\'\={\'\mathrm{float}\'}_8\right)\:$

$\mathrm{DCT}\left(A\,\'\mathrm{container}\'\=C\right)$

$\left[\begin{array}{cccccccccc}1.5000000000000004& 1.534351019701625& -0.3036659479510252& -2.5278309009474733& -2.0603989756276833& 2.2480137296073375& 1.2781204298267193& -2.6781162095883118& 2.163239928271285& -1.0329588850112363\end{array}\right]$

$C$

$\left[\begin{array}{cccccccccc}1.5000000000000004& 1.534351019701625& -0.3036659479510252& -2.5278309009474733& -2.0603989756276833& 2.2480137296073375& 1.2781204298267193& -2.6781162095883118& 2.163239928271285& -1.0329588850112363\end{array}\right]$

$\mathrm{InverseDCT}\left(C\,\'\mathrm{inplace}\'\right)$

$\left[\begin{array}{cccccccccc}-1.3877787807814457{10}^{-16}& 0.46930339512069275& 1.7633557568774196& 2.963065021785413& 1.7633557568774196& -2.1213203435596424& -1.76335575687742& 2.9630650217854124& -1.7633557568774183& 0.46930339512069374\end{array}\right]$

$C$

$\left[\begin{array}{cccccccccc}-1.3877787807814457{10}^{-16}& 0.46930339512069275& 1.7633557568774196& 2.963065021785413& 1.7633557568774196& -2.1213203435596424& -1.76335575687742& 2.9630650217854124& -1.7633557568774183& 0.46930339512069374\end{array}\right]$

$\mathrm{map}\left(\mathrm{fnormal}\,A-C\right)$

$\left[\begin{array}{cccccccccc}0.& -0.& -0.& 0.& 0.& 0.& 0.& 0.& 0.& 0.\end{array}\right]$

The SignalProcessing[DCT] and SignalProcessing[InverseDCT] commands were introduced in Maple 17.

For more information on Maple 17 changes, see Updates in Maple 17.