查看“︁量化 (信号处理)”︁的源代码

{{NoteTA
|G1=IT
|G2=Communication
|G3=Math
}}
[[File:Quantized.signal.svg|right|thumb|量化信号]]
[[File:Digital.signal.svg|right|thumb|数字信号]]
'''量化'''在[[数字信号处理]]领域是指将信号的连续取值（或者大量可能的离散取值）近似为有限多个（或较少的）离散值的过程。量化主要应用于从[[连续信号]]到[[数字信号]]的转换中。连续信号经过[[采样]]成为[[离散信号]]，离散信号经过量化即成为数字信号。注意离散信号并不需要经过采样的过程。信号的采样和量化通常都是由ADC实现的。

例如[[CD]]音频信号就是按照44100[[Hz]]的频率采样，按16[[位元]]量化为有着65536（=<math>2^{16}</math>）个可能取值的数字信号。

量化就是将模拟声音的波形转换为数字，表示采样值的二进制位数决定了量化的精度。量化的过程是先将整个幅度划分成有限个小幅度（量化阶距）的集合，把落入某个阶距内的样值归为一类，并赋予相同的量化值。

== 数学描述 ==
最简单最易懂的量化是[[标量 (数学)|标量]]（有别于多维[[向量|矢量]]）量化，开始标量量化之前先要给出输入数据。
通常，一个标量量化操作可以给出下面的描述

:<math>Q(x) = g(\lfloor f(x) \rfloor)</math>

其中
* <math>x</math>是实数，
* <math>\lfloor x \rfloor</math>是[[取整函数|下取整函数]]，生成整数<math>i = \lfloor f(x) \rfloor</math>
* <math>f(x)</math>和<math>g(i)</math>是任意的实值函数。

整数<math>i</math>是表示的数值，它通常被存储或者传输，然后在后来需要解释的时候使用<math>g(i)</math>进行最终的解释重建。整数<math>i</math>有时也称作''量化指数''。

在计算机或者其它应用，一个已知的量化方法''均匀量化''。在均匀量化方法里共有两个变量，叫''mid-rise''和''mid-tread''。

如果<math>x</math>是一个－1到1之间的数，一个mid-rise uniform量化操作，可以用"M"bit来表示量化的精度。

:<math>Q(x) = \frac{\left\lfloor 2^{M-1}x \right\rfloor+0.5}{2^{M-1}}</math>.

在这个例子中<math>f(x)</math>和<math>g(i)</math>运算符都是乘以比例因子（其中一个是另外一个的逆），并且在''g''（''i''）中带有一个偏移量以使得每个量化表示都位于输入区域的中间位置。<math>2^{-(M-1)}</math>经常称为''量化步长''。按照这个量化定律，假定在整个量化步长上量化噪声大致是[[连续型均匀分布|均匀分布]]的，并且假定量化的输入信号<math>x</math>在整个-1到1的区间大致均匀分布，量化的[[信噪比]]（SNR）可以用下面的公式计算，

:<math>
\frac{S}{N_q} \approx 20 \log_{10}(2^M)
=
6.0206 M \ \operatorname{dB}</math>.

根据这个等式，人们常说SNR大约是每[[位]]6 [[decibel|dB]]。

在mid-tread一致量化中，偏移0.5将加在下取整函数内部而不是外部。

有时候，mid-rise量化使用时不加偏移0.5。这将信号与噪声比减小了大约6.02 dB，但是当步距小的时候为了简化这是可接受的。

在数字[[电话]]系统中，两个流行的量化机制是'[[A-law algorithm|A-law]]'（在[[欧洲]]占据主导地位）和'[[Mu-law algorithm|μ-law]]'（在[[北美]]和[[日本]]占据主导地位）。这些机制将离散的模拟数值映射到8位尺度，在小值的时候近似线性随着幅度增长按照对数增加。由于人耳对于[[音量]]的感知近似对数曲线，这就使用一定的位数在可听见的声音强度范围提供了更高的信噪比。

== 忽略熵约束：Lloyd–Max量化 ==

在上面的陈述中，若令 <math>\lambda</math> 等于 0，从而忽略掉比特率约束，或等价地假设要用定长码（FLC）而非用{{le|变长码|variable-length code}}（或其他[[熵編碼法]]，如[[算术编码]]在率失真上就比定长码好）来表示量化数据，这个最优化问题就简化为了只需最小化失真 <math>D</math> 的问题了。

<math>M</math> 级量化器产生的索引可以用 <math> R = \lceil \log_2 M \rceil </math> 比特/符号的定长码。例如当 <math>M=</math>256 阶时，定长码的比特率 <math>R</math> 为 8 比特/符号。由于这个原因，这样的量化器有时称作8比特量化器。不过使用定长码消除了压缩改进，但可以通过更好的熵编码来改善。

假设 <math>M</math> 阶定长码，率失真最小化问题可以简化为失真最小化问题。
简化的问题可以陈述为：给定一个[[機率密度函數]]为 <math>f(x)</math> 的信源 <math>X</math> ，并约束量化器必须仅使用 <math>M</math> 个分类区域，求得决策边界 <math>\{b_k\}_{k=1}^{M-1} </math> 与重建层级 <math>\{y_k\}_{k=1}^M</math> 来最小化得到的失真
:<math> D=E[(x-Q(x))^2] = \int_{-\infty}^{\infty} (x-Q(x))^2f(x)dx = \sum_{k=1}^{M} \int_{b_{k-1}}^{b_k} (x-y_k)^2 f(x)dx =\sum_{k=1}^{M} d_k </math>.

对上述问题求最优解得到的量化器有时叫做MMSQE（最小均方量化误差）解，而得到的概率密度函数最优化的（非均匀）量化器叫做''Lloyd–Max''量化器，是用独立发现迭代方法<ref name=GrayNeuhoff>Robert M. Gray and David L. Neuhoff, "Quantization", ''IEEE Transactions on Information Theory'', Vol. IT-44, No. 6, pp. 2325–2383, Oct. 1998. {{doi|10.1109/18.720541}}</ref><ref>Stuart P. Lloyd, "Least Squares Quantization in PCM", ''{{tsl|en|IEEE Transactions on Information Theory||IEEE Transactions on Information Theory}}'', Vol. IT-28, pp. 129–137, No. 2, March 1982 {{doi|10.1109/TIT.1982.1056489}} (work documented in a manuscript circulated for comments at [[贝尔实验室|Bell Laboratories]] with a department log date of 31 July 1957 and also presented at the 1957 meeting of the [[Institute of Mathematical Statistics]], although not formally published until 1982).</ref><ref>Joel Max, "Quantizing for Minimum Distortion", ''{{tsl|en|IEEE Transactions on Information Theory||IRE Transactions on Information Theory}}'', Vol. IT-6, pp. 7–12, March 1960. {{doi|10.1109/TIT.1960.1057548}}</ref>从 <math> {\partial D / \partial b_k} = 0 </math> 和 <math>{\partial D/ \partial y_k} = 0 </math> 求解两组联立方程的两个人来命名的，如下：
:<math> {\partial D \over\partial b_k} = 0 \Rightarrow b_k = {y_k + y_{k+1} \over 2} </math>,
会将阈值置于每对重建值的中点，而
:<math> {\partial D \over\partial y_k} = 0 \Rightarrow y_k = { \int_{b_{k-1}}^{b_k} x f(x) dx \over \int_{b_{k-1}}^{b_k} f(x)dx } = \frac1{p_k} \int_{b_{k-1}}^{b_k} x f(x) dx </math>
会让重建值位于其相关分类区间的质心（条件期望值）。

{{le|Lloyd算法|Lloyd's algorithm|Lloyd方法I算法}}，最初于1957提出，并可以直接推广到用于[[向量量化|向量]]数据。这个推广会得到{{le|Linde–Buzo–Gray算法|Linde–Buzo–Gray algorithm|Linde–Buzo–Gray（LBG）}}或[[K-平均算法|K-平均]]分类器最优化方法。此外，此方法还可以进一步推广到对向量数据包含一个熵约束。<ref name=ChouLookabaughGray>Philip A. Chou, Tom Lookabaugh, and {{tsl|en|Robert M. Gray||Robert M. Gray}}, "Entropy-Constrained Vector Quantization", ''IEEE Transactions on Acoustics, Speech, and Signal Processing'', Vol. ASSP-37, No. 1, Jan. 1989. {{doi|10.1109/29.17498}}</ref>

== 量化与数据压缩 ==

量化在[[有损数据压缩]]中起着相当重要的作用。很多情况下，量化可以被当作将有损数据压缩同[[无损数据压缩]]相区别的标志之一。量化的目的通常是为了减少数据量。一些压缩算法，例如[[MP3]]和[[Vorbis]]，以有选择地丢弃部分数据作为压缩的一种方法，这种手段可以被认为是量化的过程也可以被看作是一种有损压缩的形式。

[[JPEG]]是一种利用了量化的图像有损压缩。JPEG的编码过程对原始的图像数据作[[离散余弦变换]]，然后对变换结果进行量化并作[[熵编码]]。通过量化可以降低变换值的精度，从而减少图像的数据量。当然，精度的损失意味着图像质量的下降。然而图像的质量可以通过量化位数的选择加以控制。例如，JPEG在每像素3[[位元|比特]]的精度下得到的图像质量还让人可以接受的，相对于[[脈衝編碼調變|PCM]]抽样得到的每个像素24比特的原始图像来说，数据量大大下降了。

现代压缩技术通常以量化输出的[[信息熵]]，而不是输出值集合的大小度量信息量的多少。

== 自然界中的量子化 ==
从最基本的意义上来说，所有的[[物理量]]都是量子化的，这是[[量子力学]]的结论。为了数学上的明晰性，在宏观的尺度上可以将量子的性质忽略，因此信号可以表示为连续的形式。

在实际应用中，这种内在的量子或量化的性质并不需要考虑。首先，量子效应会被信号的噪声淹没，因为任何观察对象在实际系统中总会伴随有其他物理现象。其次，测量仪器不可能绝对精确，被测的信号必然会被测量噪声污染。

== 量化误差 ==
量化误差是指在量化过程引起的误差，表现为量化结果和被量化模拟量之间存在差值。这种差值在输出端体现为引入了[[量化噪声]]。

== 相關條目 ==
* [[模-数转换器]]，[[数-模转换器]]
* [[量化误差]], [[量化噪声]]
* [[离散信号]]，[[数字信号]]
* [[抖动]]
* [[信息论]]
* [[率失真]]
* [[矢量量化]]

== 参考文献 ==
{{reflist}}

== 外部链接 ==
* [http://www.math.ucdavis.edu/~saito/courses/ACHA/44it06-gray.pdf Paper on mathematical theory and analysis of quantization] {{Wayback|url=http://www.math.ucdavis.edu/~saito/courses/ACHA/44it06-gray.pdf |date=20050331081934 }}
* [http://www.dsprelated.com/comp.dsp/keyword/Quantization.php Quantization threads in Comp.DSP] {{Wayback|url=http://www.dsprelated.com/comp.dsp/keyword/Quantization.php |date=20201024150249 }}

{{DSP}}

[[Category:信号处理]]