**线谱对**（LSP）或**线谱频率**（LSF）用于表示在信道上传输的线性预测系数（LPC）。LSP具有一些特性，如对量化噪声的敏感度较小，优于LPC的直接量化。因此，LSP在语音编码中非常有用。

LSP表示法是日本电信电话的板倉文忠于1975年发明的。1975年到1981年，他研究了基于LSP的语音分析与合成问题。1980年，他的团队开发出了基于LSP的语音合成芯片。LSP是语音合成和编码的一项重要技术，20世纪90年代的几乎所有国际语音编码标准都将其作为重要组成，为提高全球移动信道和互联网上的数字语音通信水平做出了很大贡献。1985年，Bishnu S. Atal、Manfred R. Schroeder基于LSP开发了CELP算法。

数学原理
----

线谱多项式{\displaystyle A(z)=1-\sum \_{k=1}^{p}a\_{k}z^{-k}}可写作{\displaystyle A(z)=0.5[P(z)+Q(z)]}，其中：

* {\displaystyle P(z)=A(z)+z^{-(p+1)}A(z^{-1})}
* {\displaystyle Q(z)=A(z)-z^{-(p+1)}A(z^{-1})}

根据构造，*P*是**回文多项式**，*Q*是**反回文多项式**；物理上，*P*(*z*)对应声门关闭时的声道，*Q*(*z*)则对应声门打开时的声道。可证明：

* *P*、*Q*零点位于复平面中的单位圆。
* 绕圆运动时，*P*与*Q*的根交替出现。
* 由于*P*、*Q*系数都是实数，因此根以共轭对的形式出现。

LP多项式的线谱对表示简单地包含了P、Q根的位置（即使{\displaystyle z=e^{i\omega },P(z)=0}的{\displaystyle \omega }）。由于根成对出现，因此只需传输一半的根（一般{\displaystyle \in (0,\ \pi )}）。因此，*P*与*Q*的系数总数等于原LP系数数*p*（不计{\displaystyle a\_{0}=1}）。
确定系数的常用算法是在单位圆上间隔较近的点串上求多项式值，观察结果何时变号；变号时，根必定位于测试点之间。由于*P*与*Q*的根穿插在一起，因此只要一次就能找到两个多项式的根。

LSP分析
-----

要转回LPC，就要计算{\displaystyle P(z),\ Q(z)}的根。下面只考虑线性预测多项式{\displaystyle A(z)}阶数为偶数{\displaystyle N}的情形，这时LSP多项式的{\displaystyle P(z),\ Q(z)}为{\displaystyle N+1}多项式。

LSP多项式的{\displaystyle P(z),\ Q(z)}可分别被{\displaystyle (1+z^{-1})}、{\displaystyle (1-z^{-1})}整除，剩余多项式用{\displaystyle (z+z^{-1})/2}除，在单位圆上可表为{\displaystyle (z+z^{-1})/2=\cos \omega }。即，{\displaystyle P(z),\ Q(z)}可进行如下因式分解：

{\displaystyle P(z)=\left(1+z^{-1}\right)\prod \_{i=1,3,...N-1}\left(1-2\cos \omega \_{i}z^{-1}+z^{-2}\right)}
{\displaystyle Q(z)=\left(1-z^{-1}\right)\prod \_{i=2,4,...N}\left(1-2\cos \omega \_{i}z^{-1}+z^{-2}\right)}
求出该式的根，便能计算线谱对{\displaystyle \omega \_{i}}。更具体地，如下</ref>
：

(1) 由线性预测系数{\displaystyle a\_{i}}计算{\displaystyle P(z),\ Q(z)}各系数

由{\displaystyle P(z),\ Q(z)}的定义，用下式计算。多项式系数{\displaystyle p\_{i},q\_{i}}，
{\displaystyle p\_{0}=p\_{N+1}=1\,}
{\displaystyle p\_{i}=p\_{N-i+1}=a\_{i}+a\_{N-i+1}\,}
{\displaystyle q\_{0}=-q\_{N+1}=1\,}
{\displaystyle q\_{i}=-q\_{N-i+1}=a\_{i}-a\_{N-i+1}\,}
(2) {\displaystyle P(z),\ Q(z)}分别除以{\displaystyle (1+z^{-1})}、{\displaystyle (1-z^{-1})}

相当于从单位圆上的根上除去实根。
此多项式除法可通过系数加减来计算。将商式系数记作{\displaystyle p',q'}，
{\displaystyle p'\_{0}=1\,}
{\displaystyle p'\_{i}=p\_{i}-p'\_{i-1}\,}
{\displaystyle q'\_{0}=1\,}
{\displaystyle q'\_{i}=q\_{i}+q'\_{i-1}\,}
(3) 商式{\displaystyle P'(z),\ Q'(z)}用{\displaystyle x=(z+z^{-1})/2}置换变量

相当于剩余复共轭根在实轴上的投影。置换后的式子可用切比雪夫多项式表示。
{\displaystyle P'(z),\ Q'(z)}是关于{\displaystyle x}的{\displaystyle N/2}次多项式，系数可从{\displaystyle p',\ q'}机械计算。
(4) 用牛顿-拉弗森法解{\displaystyle x}的两个方程

在区间<math(-1,\ 1)</math>内，根{\displaystyle x\_{i}}交替存在，则可交替求解两个方程。
(5) 由求得的根计算线谱{\displaystyle \omega \_{i}}

由求得的N个根{\displaystyle x\_{i}}求下式中的{\displaystyle \omega \_{i}}：
{\displaystyle \omega \_{i}=\arccos(x\_{i})\,}
将线谱对变换为线性预测系数时更简单，与上述相反，从线谱对{\displaystyle \omega \_{i}}求{\displaystyle P(z),\ Q(z)}各系数即可：

{\displaystyle A(z)={\frac {1}{2}}\left(P(z)+Q(z)\right)}
{\displaystyle P(z),\ Q(z)}各系数为{\displaystyle (1-2\cos \omega \_{i}z^{-1}+z^{-2})}形式的二次多项式的积，进而可作为乘以{\displaystyle (1\pm z^{-1})}的式子的系数，可以机械计算。

{\displaystyle P(z),\ Q(z)}的系数有对称性，因此能从{\displaystyle N/2}次系数通过以下公式变换为线性预测系数：

{\displaystyle {\begin{array}{lcl}a\_{i}&=&{\frac {1}{2}}\left(p\_{i}+q\_{i}\right)\\a\_{N-i+1}&=&{\frac {1}{2}}\left(p\_{i}-q\_{i}\right)\end{array}}}

性质
--

线谱对有几个有趣而有用的性质。*P*(*z*)、*Q*(*z*)的根交错排列时，只有根单调递减，滤波器的稳定性才有保证。另外，两个根越近，滤波器在相应频率上的谐振就越大。由于LSP对量化噪声不过分敏感，因此被广泛用于量化LPC滤波器。线谱频率可以内插。

另见
--

* 对数面积比

资料
--

* Speex manual （页面存档备份，存于） and source code (lsp.c)
* "The Computation of Line Spectral Frequencies Using Chebyshev Polynomials" （页面存档备份，存于）/ P. Kabal and R. P. Ramachandran. IEEE Trans. Acoustics, Speech, Signal Processing, vol. 34, no. 6, pp. 1419–1426, Dec. 1986.

Includes an overview in relation to LPC.

* "Line Spectral Pairs" chapter （页面存档备份，存于） as an online excerpt (pdf) / "Digital Signal Processing - A Computer Science Perspective" (ISBN 0-471-29546-9) Jonathan Stein.

参考文献
----

1. ↑  Sahidullah, Md.; Chakroborty, Sandipan; Saha, Goutam.. International Journal of Biometrics. 2010-01, **2** (4): 358–378. doi:10.1504/ijbm.2010.035450.
2. ↑  Zheng, F.; Song, Z.; Li, L.; Yu, W.  (PDF). Proceedings of the 5th International Conference on Spoken Language Processing (ICSLP'98). 1998, (3): 1123–6  [2023-10-21]. （原始内容存档 (PDF)于2022-10-09）.
3. 1 2 . IEEE.  [2019-07-15]. （原始内容存档于2015-02-10）.
4. ↑ . IEEE Global History Network. 2009-05-20  [2009-07-21]. （原始内容存档于2015-04-03）.
5. ↑   （页面存档备份，存于） Tony Robinson: Speech Analysis
6. ↑  e.g. lsf.c in  （页面存档备份，存于）
7. ↑  嵯峨山茂樹. *応用音響学: 音声分析(5) LSP分析*.(pdf) 東京大学 応用音響学 講義資料.
8. 1 2  Peter Kabal, Ravi P. Ramachandran. *The Computation of Line Spectral Frequencies Using Chebyshev Polynomials （页面存档备份，存于）*.(pdf) IEEE Trans. Acoustics, Speech, Signal Processing, vol.34, no. 6, pp.1419-1426, Dec. 1986.
9. 1 2  Wai C. Chu. *Speech Coding Algorithms: Foundation and Evolution of Standardized Coders*. pp.239-250, 2003.

|  |  |  |  |  |  |  |  |  |  |  |  |  |  |  |  |  |  |  |  |  |  |  |  |  |  |  |  |  |  |  |  |  |  |  |  |  |  |  |  |  |  |  |  |  |  |  |  |  |  |
| --- | --- | --- | --- | --- | --- | --- | --- | --- | --- | --- | --- | --- | --- | --- | --- | --- | --- | --- | --- | --- | --- | --- | --- | --- | --- | --- | --- | --- | --- | --- | --- | --- | --- | --- | --- | --- | --- | --- | --- | --- | --- | --- | --- | --- | --- | --- | --- | --- | --- |
| 

| * 
* 论
* 编
方法 |
| --- |
|  |
| 理论 | 
* 熵
* 语法（英语：）
* 信息论
* 柯氏复杂度
* 量化
* 率失真理论
* 

 |
|  |
|  | 

|  |  |
| --- | --- |
| 熵編碼 | 
* 算术编码
* 非对称数字系统编码（英语：）
* 格倫布編碼
* 霍夫曼编码
	+ 适应性
	+ 范型
	+ 改进
* 區間編碼
* 香农编码
* 香农-范诺编码
* 香農-范諾-埃利西斯編碼
* Tunstall编码（英语：）
* 一元编码（英语：）
* 通用編碼（英语：）
	+ 指數格倫布編碼
	+ 斐波那契編碼
	+ 以利亞伽瑪編碼
	+ 萊文斯坦編碼（英语：）

 |
|  |
| 字典編碼（英语：） | 
* 字节对编码
* LZ77 / LZ78
	+ 842（英语：）
	+ Brotli
	+ Deflate
	+ LZ4
	+ LZFSE
	+ LZJB（英语：）
	+ LZMA
	+ LZO
	+ LZRW（英语：）
	+ LZS（英语：）
	+ LZSS
	+ LZW
	+ LZWL（英语：）
	+ LZX
	+ Snappy
	+ Zstandard

 |
|  |
| 其他 | 
* BWT
* CTW（英语：）
* 差分编码
* DMC
* 差分脉冲编码调制
* 无损离散余弦变换
* MTF
* PAQ（英语：）
* PPM（英语：）
* RLE

 |

 |
|  |
|  | 

|  |  |
| --- | --- |
| 变换编码 | 
* DCT
	+ MDCT
* DST
* FFT
* 小波变换
	+ 多贝西
	+ DWT
	+ SPIHT（英语：）

 |
|  |
| 预测编码 | 
* DPCM
	+ ADPCM（英语：）
* LPC
	+ ACELP（英语：）
	+ CELP
	+ LAR（英语：）
	+ LSP
	+ WLPC（英语：）
* 运动
	+ 运动补偿
	+ 运动预测
	+ 运动矢量
* 心理声学

 |

 |
|  |
|  | 

|  |  |
| --- | --- |
| 概念 | 
* 码率
	+ 平均码率
	+ 恒定码率
	+ 可变码率
* 压扩（英语：）
* 卷积
* 动态范围
* 延迟（英语：）
* 采样
* 采样定理
* 音质
* 语音编码
* 子带编码

 |
|  |
| 编解码组件 | 
* A-law（英语：）
* μ-law（英语：）
* DPCM
	+ ADPCM（英语：）
	+ DM
* FT
	+ FFT
* LPC
	+ ACELP（英语：）
	+ CELP
	+ LAR（英语：）
	+ LSP
	+ WLPC（英语：）
	+ CELP
	+ MDCT
* 心理聲學模型

 |

 |
|  |
|  | 

|  |  |
| --- | --- |
| 概念 | 
* 色度抽样
* 编码树单元
* 色彩空間
* 壓縮失真
* 图像分辨率
* 宏块
* 像素
* 峰值信噪比
* 量化
* 标准测试图像

 |
|  |
| 方法 | 
* DCT
* Deflate
* 分形压缩
* K-L变换
* LP（英语：）
* RLE
* 小波变换
	+ 多贝西
	+ DWT
	+ SPIHT（英语：）

 |

 |
|  |
|  | 

|  |  |
| --- | --- |
| 概念 | 
* 码率
	+ 平均码率
	+ 恒定码率
	+ 可变码率
* 显示分辨率
* 帧（英语：）
* 帧率
* 帧类型
* 隔行扫描
* 視訊特性
* 視訊質量

 |
|  |
| 编解码组件 | 
* DCT
* DPCM
* 去区块滤波器
* 重叠变换（英语：）
* 运动
	+ 运动补偿
	+ 运动预测
	+ 运动矢量
* 量化

 |

 |
|  |
| 另见压缩格式和数据压缩 |

 |