線譜對

線譜對（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算法。

數學原理

線譜多項式$A(z)=1-\sum _{k=1}^{p}a_{k}z^{-k}$可寫作$A(z)=0.5[P(z)+Q(z)]$，其中：

• $P(z)=A(z)+z^{-(p+1)}A(z^{-1})$
• $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根的位置（即使$z=e^{i\omega },P(z)=0$的$\omega$）。由於根成對出現，因此只需傳輸一半的根（一般$\in (0,\ \pi )$）。因此，P與Q的係數總數等於原LP係數數p（不計$a_{0}=1$）。 確定係數的常用算法是在單位圓上間隔較近的點串上求多項式值，觀察結果何時變號；變號時，根必定位於測試點之間。由於P與Q的根穿插在一起，因此只要一次就能找到兩個多項式的根。

LSP分析

要轉回LPC，就要計算$P(z),\ Q(z)$的根。下面只考慮線性預測多項式$A(z)$階數為偶數$N$的情形，這時LSP多項式的$P(z),\ Q(z)$為$N+1$多項式。

LSP多項式的$P(z),\ Q(z)$可分別被$(1+z^{-1})$、$(1-z^{-1})$整除，剩餘多項式用$(z+z^{-1})/2$除，在單位圓上可表為$(z+z^{-1})/2=\cos \omega$。即，$P(z),\ Q(z)$可進行如下因式分解：

$P(z)=\left(1+z^{-1}\right)\prod _{i=1,3,...,N-1}\left(1-2\cos \omega _{i}z^{-1}+z^{-2}\right)$

$Q(z)=\left(1-z^{-1}\right)\prod _{i=2,4,...,N}\left(1-2\cos \omega _{i}z^{-1}+z^{-2}\right)$

求出該式的根，便能計算線譜對$\omega _{i}$。更具體地，如下</ref> ：

(1) 由線性預測係數$a_{i}$計算$P(z),\ Q(z)$各係數

由$P(z),\ Q(z)$的定義，用下式計算。多項式係數$p_{i},q_{i}$，

$p_{0}=p_{N+1}=1\,$

$p_{i}=p_{N-i+1}=a_{i}+a_{N-i+1}\,$

$q_{0}=-q_{N+1}=1\,$

$q_{i}=-q_{N-i+1}=a_{i}-a_{N-i+1}\,$

(2) $P(z),\ Q(z)$分別除以$(1+z^{-1})$、$(1-z^{-1})$

相當於從單位圓上的根上除去實根。

此多項式除法可通過係數加減來計算。將商式係數記作$p',q'$，

$p'_{0}=1\,$

$p'_{i}=p_{i}-p'_{i-1}\,$

$q'_{0}=1\,$

$q'_{i}=q_{i}+q'_{i-1}\,$

(3) 商式$P'(z),\ Q'(z)$用$x=(z+z^{-1})/2$置換變量

相當於剩餘復共軛根在實軸上的投影。置換後的式子可用切比雪夫多項式表示。

$P'(z),\ Q'(z)$是關於$x$的$N/2$次多項式，係數可從$p',\ q'$機械計算。

(4) 用牛頓-拉弗森法解$x$的兩個方程

在區間<math(-1,\ 1)</math>內，根$x_{i}$交替存在，則可交替求解兩個方程。

(5) 由求得的根計算線譜$\omega _{i}$

由求得的N個根$x_{i}$求下式中的$\omega _{i}$：

$\omega _{i}=\arccos(x_{i})\,$

將線譜對變換為線性預測係數時更簡單，與上述相反，從線譜對$\omega _{i}$求$P(z),\ Q(z)$各係數即可：

$A(z)={\frac {1}{2}}\left(P(z)+Q(z)\right)$

$P(z),\ Q(z)$各係數為$(1-2\cos \omega _{i}z^{-1}+z^{-2})$形式的二次多項式的積，進而可作為乘以$(1\pm z^{-1})$的式子的係數，可以機械計算。

$P(z),\ Q(z)$的係數有對稱性，因此能從$N/2$次係數通過以下公式變換為線性預測係數：

${\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.