该文讨论了当量化器的输出电平数固定时,量化器对信号失真最小化的问题。失真被建模为量化器的输入和输出之间的误差函数的期望值。推导了失真最小的量化器参数的方程。我们着重讨论一下均匀分布的情况。
1. 量化背景
在许多数据传输系统中,模拟输入信号首先在发射机上转换为数字形式,然后以数字形式传输,最后在接收机上重构为模拟信号。所产生的输出通常与输入信号相似,但并不完全相同,因为发射机上的量化器对位于有限数量的振幅范围内的所有输入振幅产生相同的数字。假设数字传输无误差,输入信号和输出信号之间的差值是量化误差。由于任何系统的数字传输速率都是有限的,因此必须使用一个量化器来将输入分类为有限数量N的范围。
2. 数学表达
对于给定的输出区间数量N,
x
k
x_k
xk表示为各个区间的边界值,
y
k
y_k
yk表示为重建值:
C
i
=
{
x
∈
R
:
Q
(
x
)
=
x
i
}
C_i=\{x\in R:Q(x)=x_i\}
Ci={x∈R:Q(x)=xi}
⋃
i
=
0
N
−
1
C
i
=
R
w
i
t
h
⋁
i
≠
j
:
C
i
⋂
C
j
=
∅
\bigcup \limits_{i=0}^{N-1}C_i=R \ with \bigvee_i \neq j:C_i\bigcap C_j=\emptyset
i=0⋃N−1Ci=R with⋁i=j:Ci⋂Cj=∅
下图中边界值用的
u
i
u_i
ui,重建值用的
S
i
′
S_i^{'}
Si′。
3. Minimum Distortion
通常为量化过程定义一个失真度量是合适的,这将是量化误差的一些统计量。
f
(
x
)
f(x)
f(x)为损失定义函数。
D
=
E
[
f
(
s
i
n
−
s
o
u
t
)
]
=
∑
i
=
1
N
∫
x
i
x
i
+
1
f
(
x
−
y
i
)
p
(
x
)
d
x
\begin{aligned} D&=E[f(s_{in}-s_{out})] \\ &=\sum \limits_{i=1}^{N} \int_{x_i}^{x_{i+1}}f(x-y_i)p(x)dx \end{aligned}
D=E[f(sin−sout)]=i=1∑N∫xixi+1f(x−yi)p(x)dx
对于定义域无限的x来说,
x
N
+
1
x_{N+1}
xN+1和
x
1
=
−
∞
x_{1}=-\infty
x1=−∞。
如果我们希望使固定N的D最小化,我们通过将D分别对于
x
i
x_i
xi和
y
i
y_i
yi求偏导,并将导数设为零来得到必要条件:
∂
D
∂
x
i
=
f
(
x
i
−
y
i
−
1
)
p
(
x
)
−
f
(
x
i
−
y
i
)
p
(
x
i
)
=
0
i
=
2
,
.
.
.
,
N
(
1
)
\begin{aligned}\frac{\partial D}{\partial x_i}=f(x_i -y_{i-1})p(x)-f(x_i-y_i)p(x_i)=0 \\i=2,...,N\ \ \ \ \ \ \ (1) \end{aligned}
∂xi∂D=f(xi−yi−1)p(x)−f(xi−yi)p(xi)=0i=2,...,N (1)
∂ D ∂ y i = − ∫ x i x i + 1 f ′ ( x − y i ) p ( x ) d x = 0 i = 1 , . . . , N ( 2 ) \begin{aligned}\frac{\partial D}{\partial y_i} = -\int_{x_i}^{x_{i+1}} f^{'}(x-y_i)\ p(x)\ dx=0 \\i=1,...,N\ \ \ \ \ \ \ (2)\end{aligned} ∂yi∂D=−∫xixi+1f′(x−yi) p(x) dx=0i=1,...,N (2)
(
1
)
b
e
c
o
m
e
s
(
f
o
r
p
(
x
i
)
≠
0
)
f
(
x
i
−
y
i
−
1
)
=
f
(
x
i
−
y
i
)
i
=
2
,
.
.
.
,
N
(
3
)
\begin{aligned}(1)becomes (for\ p(x_i) \neq 0) \\& f(x_i -y_{i-1})=f(x_i-y_i)& &i=2,...,N\ \ \ \ \ \ \ (3) \end{aligned}
(1)becomes(for p(xi)=0)f(xi−yi−1)=f(xi−yi)i=2,...,N (3)
(
2
)
b
e
c
o
m
e
s
∫
x
i
x
i
+
1
f
′
(
x
−
y
i
)
p
(
x
)
d
x
=
0
i
=
2
,
.
.
.
,
N
(
4
)
\begin{aligned}(2)becomes \\ &&&&&&\int_{x_i}^{x_{i+1}} f^{'}(x-y_i)\ p(x)\ dx=0& &i=2,...,N\ \ \ \ \ \ \ (4) \end{aligned}
(2)becomes∫xixi+1f′(x−yi) p(x) dx=0i=2,...,N (4)
在一般情况下最好的解是,如果所有的二阶偏导数的D对
x
i
x_i
xi和
y
i
y_i
yi的存在,那么临界点由条件(3)和(4)是一个最小如果矩阵第
i
i
i行和
j
j
j列元素满足正定:
∂
2
D
∂
p
i
∂
p
j
∣
c
r
i
t
i
c
a
l
p
o
i
n
t
\frac{\partial^{2}D}{\partial p_i \partial p_j}|_{critical \ point}
∂pi∂pj∂2D∣critical point
通常的情况下,我们采取
f
(
x
)
=
x
2
f(x)=x^2
f(x)=x2。
(
3
)
−
>
(
j
=
2
,
.
.
.
,
N
)
(3)->(j=2,...,N)
(3)−>(j=2,...,N)
x
i
=
(
y
i
+
y
i
+
1
)
/
2
o
r
y
i
=
2
x
i
−
y
i
−
1
x_i=(y_i+y_{i+1})/2\ or y_i=2x_i-y_{i-1}
xi=(yi+yi+1)/2 oryi=2xi−yi−1
(
4
)
−
>
(
i
=
1
,
.
.
.
,
N
)
(4)->(i=1,...,N)
(4)−>(i=1,...,N)
∫
x
i
x
i
+
1
(
x
−
y
i
)
p
(
x
)
d
x
=
0
\int_{x_i}^{x_{i+1}}(x-y_i)p(x)dx=0
∫xixi+1(x−yi)p(x)dx=0
y i y_i yi是在 x i , x i + 1 x_i, x_{i+1} xi,xi+1之间的质心。
3.1 例子
针对信源分布为
p
(
x
)
=
1
/
2
e
−
x
2
/
2
p(x)=1/ \sqrt{2}e^{-x^2/2}
p(x)=1/2e−x2/2,其中
x
N
/
2
+
1
=
0
f
o
r
N
e
v
e
n
o
r
y
(
N
+
1
)
/
2
=
0
f
o
r
N
o
d
d
x_{N/2+1}=0 for\ N\ even\ or\ y_{(N+1)/2}=0\ for\ N\ odd
xN/2+1=0for N even or y(N+1)/2=0 for N odd,这样会有一个对称的结果。具体为:
图中是失真与输出水平数的对数-对数图1。
当N很大时,振幅概率密度从单个输入范围的一端到另一端没有明显的变化,除了非常大的振幅,这是非常不可能的,因此它们的影响很小。因此,大多数输出级别都非常接近于对应输入范围的端点的平均值。现在,量化均匀分布的输入信号的最佳方法是均匀地划分输出电平,并将输入范围的端点位于输出电平之间,如图2
为这个分布产生具有2N输出水平的量化器的最佳方法是将原来的区间分成两半,并将新的输出水平放在这些范围的中点,如图所3
4. Uniform Quantization PCM
商用高速模数转换设备目前仅限于将等量的输入范围转换为介于输入范围两端之间的输出。在许多应用程序中,人们希望知道要使用的最佳间隔长度,即,对于给定数量的输出级别N,产生最小失真的间隔长度。这是一个比第一个问题更容易的问题,因为它只是二维的(对于
N
≥
2
N\ge 2
N≥2),也就是说,D是间隔的长度r和任何特定输出
y
k
y_k
yk的函数。如果输入具有对称分布,并且需要对称答案,则问题就变成一维问题;此时,我们的重建值为各个区间的中间值。如果p(x)是输入振幅概率密度,f(x)是函数,使得失真D是
E
[
f
(
s
o
u
t
−
s
i
n
)
]
E[f(s_{out}-s_{in})]
E[f(sout−sin)],那么,对于偶数为2N的输出:
D
=
2
∑
i
=
1
N
−
1
∫
(
i
−
1
)
r
i
r
f
(
x
−
[
2
i
−
1
2
]
)
p
(
x
)
d
x
+
2
∫
(
N
−
1
)
r
∞
f
(
x
−
[
2
N
−
1
2
]
)
p
(
x
)
d
x
\begin{aligned} D=2\sum \limits_{i=1}^{N-1}\int_{(i-1)r}^{ir}f(x-[\frac{2i-1}{2}])p(x)dx+2\int_{(N-1)r}^{\infty}f(x-[\frac{2N-1}{2}])p(x)dx \end{aligned}
D=2i=1∑N−1∫(i−1)rirf(x−[22i−1])p(x)dx+2∫(N−1)r∞f(x−[22N−1])p(x)dx
我们对r求偏导:
d
D
d
r
=
−
∑
i
=
1
N
−
1
(
2
i
−
1
)
∫
(
i
−
1
)
r
i
r
f
′
(
x
−
[
2
i
−
1
2
]
)
p
(
x
)
d
x
−
(
2
N
−
1
)
∫
(
N
−
1
)
r
∞
f
′
(
x
−
[
2
N
−
1
2
]
)
p
(
x
)
d
x
=
0
\begin{aligned} \frac{dD}{dr}=- \sum \limits_{i=1}^{N-1}(2i-1)\int_{(i-1)r}^{ir}f^{'}(x-[\frac{2i-1}{2}])p(x)dx-(2N-1)\int_{(N-1)r}^{\infty}f^{'}(x-[\frac{2N-1}{2}])p(x)dx=0 \end{aligned}
drdD=−i=1∑N−1(2i−1)∫(i−1)rirf′(x−[22i−1])p(x)dx−(2N−1)∫(N−1)r∞f′(x−[22N−1])p(x)dx=0
对于奇数的情况也存在类似的表达式。在这两种情况下,当指定了f(x)、p(x)和N时,这个问题都很容易受到机器计算的影响。
针对信源分布为
p
(
x
)
=
1
/
2
e
−
x
2
/
2
p(x)=1/ \sqrt{2}e^{-x^2/2}
p(x)=1/2e−x2/2,其中
x
N
/
2
+
1
=
0
f
o
r
N
e
v
e
n
o
r
y
(
N
+
1
)
/
2
=
0
f
o
r
N
o
d
d
x_{N/2+1}=0 for\ N\ even\ or\ y_{(N+1)/2}=0\ for\ N\ odd
xN/2+1=0for N even or y(N+1)/2=0 for N odd,N=2 to 36。
图中显示了等间距的输出水平间距与输出数的对数-对数图4:
最后,在图中可以看到最佳量化器的Distortion比率与最佳等距电平量化器的比值图5
5. 均匀量化与非均匀量化
5.1 Computation of the Compressor Function
我们根据C()的斜率得到:
d
x
(
s
)
d
s
∣
s
=
s
i
′
≈
Δ
Δ
i
\frac{dx(s)}{ds}|_{s=s_i^{'}}\approx \frac{\Delta}{\Delta_i}
dsdx(s)∣s=si′≈ΔiΔ
每区间MSE贡献相等的条件为:(考虑的是用FLC进行标量量化器的高码率近似)
f
(
s
i
′
)
3
Δ
i
≈
1
K
∫
−
∞
∞
f
(
s
i
)
3
d
s
\sqrt[3]{f(s_i^{'})} \Delta_i \approx \frac{1}{K}\int_{-\infty}^{\infty}\sqrt[3]{f(s_i)}ds
3f(si′)Δi≈K1∫−∞∞3f(si)ds
结合上述两个式子:
d
x
(
s
)
d
s
∣
s
=
s
i
′
≈
Δ
Δ
i
=
Δ
f
(
s
i
′
)
3
Δ
i
f
(
s
i
′
)
3
=
Δ
⋅
K
f
(
s
i
′
)
3
∫
−
∞
∞
f
(
x
)
3
d
x
\frac{dx(s)}{ds}|_{s=s_i^{'}}\approx \frac{\Delta}{\Delta_i}=\frac{\Delta \sqrt[3]{f(s_i^{'})}}{\Delta_i \sqrt[3]{f(s_i^{'})}}=\Delta \cdot K\frac{\sqrt[3]{f(s_i^{'})}}{\int_{-\infty}^{\infty}\sqrt[3]{f(x)}dx}
dsdx(s)∣s=si′≈ΔiΔ=Δi3f(si′)Δ3f(si′)=Δ⋅K∫−∞∞3f(x)dx3f(si′)
Compressor function is obtained by integrating over its slope:
c
(
s
)
=
K
Δ
∫
−
∞
∞
f
(
x
)
3
⋅
∫
−
∞
s
f
(
x
)
3
d
x
−
K
Δ
2
c(s)=\frac{K \Delta}{\int_{-\infty}^{\infty}\sqrt[3]{f(x)}} \cdot \int_{-\infty}^{s}\sqrt[3]{f(x)}dx-\frac{K \Delta}{2}
c(s)=∫−∞∞3f(x)KΔ⋅∫−∞s3f(x)dx−2KΔ
Reference:1.source coding
2.Quantizing for Minimum Distortion*