idct
逆离散余弦变换。
库::`工程师`
例子:
使用反离散余弦变换进行信号恢复
Details
我们将产生一个由频率为 25 Hz,以频率采样 1000 Hz期间 1 几秒钟。 将具有方差的白高斯噪声添加到正弦波中 0.01.
import EngeeDSP.Functions: randn
using Random
Random.seed!(1234)
Fs = 1000
t = 0:1/Fs:1-1/Fs
x = sin.(2π * 25 * t) + randn(length(t), 1) / 10
让我们计算序列的离散余弦变换(DCT)。 让我们来确定有多少 1000 准备系数是显着的。 让我们选择一个重要的门槛 1.
import EngeeDSP.Functions: dct
y=dct(x)
sigcoeff=abs。(y)。>= 1
howmany=sum(sigcoeff)
17
仅使用显着分量重构信号。
import EngeeDSP.Functions: idct
y[.!sigcoeff] .= 0
z = idct(y)
让我们绘制原始信号和重建信号。
p1 = plot(t, x, title = "Original")
yl = ylims(p1)
p2 = plot(t, z, title = "Reconstructed")
ylims!(p2, yl)
plot(p1, p2, layout = (2, 1), legend = false)
此外
逆离散余弦变换
Details
逆离散余弦变换从离散余弦变换(DCT)的其系数重构序列。 功能 idct 是函数的倒数 dct.
PREP有四个标准选项。 此功能仅实现DKP-2类型。 对于转换后的信号 长度 和克朗克符号 逆变换的定义如下:
-
反向准备-1:
-
反向准备-2:
-
反向准备-3:
-
反向准备-4:
这些行是从 和 而不是通常的 和 由于 AnyMath 向量从 以前 ,不是从 以前 .
所有PREP变体都是_unitar_(或者,等价地,_orthogonal):要查找直接转换,请交换 和 在每一个定义中。 DKP-1和DKP-4彼此相反。 DKP-2和DKP-3彼此相反。