MPEG聲音編碼的單片DSP實現(xiàn)

MPEG聲音編碼的單片DSP實現(xiàn)

1　MPEG聲音編碼原理 

    MPEG聲音編碼是一種基于人耳聽覺特性的子帶聲音編碼算法，它屬于一種感覺聲音編碼方法.感覺聲音編碼算法的基本結(jié)構(gòu)如圖1所示.根據(jù)編碼器著重于頻率分辨率還是時間分辨率，可分為子帶編碼器和變換編碼器.MPEG聲音第2層編碼算法在頻域上把聲音信號劃分為32個子帶，屬于一種子帶編碼器.在圖1中，時頻映射也稱濾波器組，用于把輸入的聲音信號映射成亞抽樣的頻率分量.根據(jù)使用的濾波器組的性質(zhì)，即濾波器組在頻域的分辨率的大小，這些頻率分量又可叫做子帶樣值或頻率線.

     (a)  　　                                                  (b) 　　圖1　感覺聲音***結(jié)構(gòu)框圖

    濾波器的輸出或者與濾波器組并行的時頻變換的輸出，提供給心理聲學(xué)模型以估計時間相關(guān)的聲音掩蔽門限.心理聲學(xué)模型使用了人們所知的同時掩蔽效應(yīng)，包括有調(diào)音的掩蔽特性和無調(diào)音的掩蔽特性.如果使用聲音的前后掩蔽效應(yīng)，還可進(jìn)一步提高掩蔽門限估計的準(zhǔn)確性.子帶樣值或頻率線按照盡量保證量化噪聲的頻譜處于掩蔽門限以下的準(zhǔn)則進(jìn)行量化和編碼，這樣能保證被人耳感知的量化引入的噪聲最小.根據(jù)對復(fù)雜度的要求，可以使用塊壓擴(kuò)或熵編碼的分析合成方法.

    幀打包把量化編碼的輸出和相關(guān)邊信息按照規(guī)定的格式組合起來，以便供***使用.

    2　編碼質(zhì)量和DSP速度

    單片ADSP-2181實現(xiàn)MPEG聲音編碼關(guān)鍵需要解決兩個問題：一是如何保證聲音編碼質(zhì)量；其次是如何充分利用DSP的運算速度.而這兩個問題往往又是一對矛盾，需要找到其最佳結(jié)合點.

    一般而言，決定MPEG聲音編碼器的優(yōu)劣主要是聲學(xué)模型的好壞.但是，對于使用單片16bit定點DSP的應(yīng)用而言，這個結(jié)論就不再適用了.分析表明，此時有限字長效應(yīng)對編碼質(zhì)量的影響成了主要矛盾.特別是分析濾波器組，截尾效應(yīng)竟帶來了33倍于16bitAD轉(zhuǎn)換量化誤差的噪聲，而窗系數(shù)的有限長度表示則使本來高達(dá)96dB旁瓣衰減的濾波器響應(yīng)降低到不到70dB.因此，要保證聲音編碼質(zhì)量，分析濾波器組算法必須進(jìn)行精度擴(kuò)展.

    關(guān)于速度問題，首先想到的是使用快速算法，我們也嘗試了在子帶濾波中使用快速算法［4］.
但是，實踐證明，這些快速算法使用在DSP上效果并不理想，其原因有以下3條：（1）只考慮了加法和乘法的次數(shù)，而對附值、尋址等操作毫不關(guān)心，但對所有指令都是單周期的DSP而言，乘法和加法的次數(shù)相對其他操作并不顯得特別重要；（2）沒有考慮DSP的硬件特點，其算法不能充分發(fā)揮DSP的乘累加器（MAC）并行處理的能力；（3）ADSP-2181是為16位算法操作優(yōu)化的，在需要精度擴(kuò)展的情況下，運算量將以數(shù)量級的速度急劇增加.

    基于以上質(zhì)量和速度要求的分析，我們選用了適合DSP乘累加指令的多相結(jié)構(gòu)濾波器組實現(xiàn)方式，且采用基于MAC結(jié)構(gòu)的精度擴(kuò)展方法，較好地解決了編碼質(zhì)量和DSP速度之間的矛盾.另外，對抽樣數(shù)據(jù)的輸入方式、心理聲學(xué)模型、比例因子編碼都進(jìn)行了適于ADSP-2181的改進(jìn)，減少了運算量，保證了實時性.

    3　算法的軟件設(shè)計

    軟件設(shè)計是MPEG聲音編碼的單片DSP實現(xiàn)的核心，編碼質(zhì)量和速度的要求都需要通過精心設(shè)計DSP軟件才能實現(xiàn).
    (1)基于MAC結(jié)構(gòu)的精度擴(kuò)展　MPEG聲音編碼的分析濾波器組可以有許多種實現(xiàn)方式，多相結(jié)構(gòu)是MPEG標(biāo)準(zhǔn)推薦的一種，其數(shù)學(xué)表示為

    (1) 
    (2) 

    分析表明，對Yk進(jìn)行雙字?jǐn)U展可將截尾效應(yīng)帶來的噪聲降低33倍.但是，考慮到ADSP-2181只支持16bit的乘累加運算，需要對式（1）進(jìn)行轉(zhuǎn)化，即

     (3)　 Yk=HYk+2-16LYk　 

    這樣，就可利用DSP的乘累加器結(jié)構(gòu)，運算量只增加約1倍，存儲量只增加64個字.
                  　　

    (2)輸入數(shù)據(jù)的組織　輸入數(shù)據(jù)的組織不但要考慮方便地從數(shù)模轉(zhuǎn)換器取得聲音原始數(shù)據(jù)，還要考慮輸入數(shù)據(jù)在片內(nèi)數(shù)據(jù)RAM的存儲適合作為多相濾波器組和聲學(xué)模型的FFT運算的輸入.多相濾波器組每次移入32個新的聲音數(shù)據(jù)，移出32個舊的樣值，操作如下：

          Xi=Xi-32，i=511,510,…,32
          Xi=next-input-audio-sample，i=31,30,…,0

    然而ADSP-2181并不適于實現(xiàn)數(shù)據(jù)的移動，每個賦值運算需要兩個指令才能完成，每次分析濾波操作需要1024個指令周期.如果利用ADSP-2181的多通道自動緩沖串口及間接尋址能力，適當(dāng)?shù)亟M織輸入聲音數(shù)據(jù)，就可利用滑動窗的方法實現(xiàn)數(shù)據(jù)的移入和移出，如圖2所示. 

                               圖2　多相濾波的滑動窗技術(shù)

    為了保證幀邊界處理的連續(xù)性，輸入數(shù)據(jù)緩存應(yīng)該設(shè)計成圓緩沖的形式，其長度應(yīng)能存儲兩幀聲音輸入數(shù)據(jù).當(dāng)DSP在處理一幀數(shù)據(jù)時，輸入數(shù)據(jù)可以緩沖到另一幀.這樣，數(shù)據(jù)移動的開銷就節(jié)約了.同時，輸入數(shù)據(jù)的組織還要利于聲學(xué)模型的FFT運算，F(xiàn)FT需要利用ADSP-2181的地址反轉(zhuǎn)尋址模式.由于FFT計算和輸入數(shù)據(jù)的緩存是同時進(jìn)行的，所以FFT計算的指針需要地址反轉(zhuǎn)，而輸入緩沖的指針卻不能地址反轉(zhuǎn)，否則會導(dǎo)致輸入聲音數(shù)據(jù)排列混亂.ADSP-2181提供這種能力，它的第一地址指針組I0，I1，I2，I3有地址反轉(zhuǎn)能力，而第二地址指針組I4，I5，I6，I7卻不受地址反轉(zhuǎn)模式的影響.所以從第二地址指針組中選擇指針進(jìn)行輸入緩沖，從第一地址指針組中選擇指針進(jìn)行FFT計算.

    (3)聲學(xué)模型的改進(jìn)　用DSP實現(xiàn)心理聲學(xué)模型的一個難題是其中有大量的對數(shù)運算，雖然可以用多項式逼近求得其近似值，但是其巨大的運算量說明這不是一個明智的選擇.在改進(jìn)的心理聲學(xué)模型中，F(xiàn)FT運算后并不立即換算到對數(shù)域，而是用分段折線逼近線性域的掩蔽效應(yīng)曲線.為簡單起見，使用與標(biāo)準(zhǔn)一致的分段方法.逼近采用取指數(shù)的多項式展開的一次項的方法，這種方法雖然比較粗糙，但正如前面分析的那樣，聲學(xué)模型在16bit定點實現(xiàn)時不是主要矛盾，因而還是可以接受的.

    得到掩蔽門限以后，為計算信掩比供比特分配使用，還是需要從線性域轉(zhuǎn)換到對數(shù)域.這時，我們采用一種利用ADSP-2181移位器的近似計算方法.通過EXP指令，可以提取2進(jìn)制補(bǔ)碼小數(shù)的指數(shù)，對能量而言又有1bit約3dB.因而指數(shù)值乘3就近似得到該補(bǔ)碼小數(shù)的dB值，尾數(shù)部分的影響忽略不計.

    (4)比例因子的編碼　MPEG聲音編碼標(biāo)準(zhǔn)中一共給出了63個比例因子，但是并不是所有這些比例因子都可以用16bit的2進(jìn)制數(shù)表示.如果用雙字進(jìn)行精度擴(kuò)展，在量化時又將面臨雙字除法的巨大開銷，因此，只使用其中可以用16bit的2進(jìn)制補(bǔ)碼小數(shù)精確表示的子集，即序號為3的倍數(shù)且小于等于45的比例因子.
采用比例因子子集后，比例因子編碼就可以不再通過比較的方法得到，而可以直接通過計算子帶最大幅度的指數(shù)獲得，簡化了比例因子的編碼.

    (5)軟件仿真結(jié)果　結(jié)合上述各項算法改進(jìn)，根據(jù)ADSP-2181的特點和MPEG標(biāo)準(zhǔn)，用AD公司的開發(fā)軟件進(jìn)行了軟件仿真.表1列出了仿真得到的各個模塊對運算量和存儲量要求進(jìn)行的估算結(jié)果.仿真在抽樣率為48kHz，編碼模式為立體聲，輸入信號為頻率為1kHz的正弦波，輸出碼率為 192kbit/s的情況下進(jìn)行. 由表1可知，ADSP-2181的性能得到了較充分的利用.仿真結(jié)果表明，在以上的條件下，解碼輸出的信噪比可達(dá)80dB左右.可見，所作的算法改進(jìn)是比較有效的.

    表1　各模塊的運算量和存儲量要求

    運算量/(106指令/s) 程序存儲量/103字?jǐn)?shù)據(jù)存儲量/103字 子帶濾波 18 3.0 6.5 聲學(xué)模型 103.5 1.5 比特分配和量化 2 2.0 — 格式化比特流 1 0.5 1.0

    4　硬件設(shè)計

    硬件結(jié)構(gòu)框圖如圖3所示.各模塊的基本功能如下：

      
                  　　
    DSP核：除完成所有編碼算法以外，還要完成對模數(shù)轉(zhuǎn)換電路的初始化配置；通過輔助控制電路選擇抽樣時鐘，通過接口電路接受主機(jī)的編碼參數(shù).

    輔助控制電路：由FPGA及附屬電路實現(xiàn)，完成時鐘的產(chǎn)生、FIFO狀態(tài)的監(jiān)測、地址譯碼等功能. 

    輸出緩沖：編碼碼流的暫時存儲區(qū)，同時提供完全異步的輸出接口方式.在需要實現(xiàn)圖象聲音唇形同步的應(yīng)用中特別有用.

    外部存儲器：包括BDMA空間、I/O空間.

    模數(shù)轉(zhuǎn)換電路：完成聲音的數(shù)字化，直接與DSP的串口0連接.抽樣頻率由外部提供的256倍抽樣時鐘的頻率決定，正常工作前需要進(jìn)行初始化.

    接口電路：接口電路分為兩部分，一部分是編碼輸出接口，另一部分是與主機(jī)連接的接口.主機(jī)接口使用RS232接口芯片完成DSP串口1與主機(jī)串口的連接，DSP使用中斷和內(nèi)部計時器實現(xiàn)異步串型通信.

    上述方案已經(jīng)在“九五”科技攻關(guān)項目中實現(xiàn)，實時編解碼的聲音通過了主觀測試.
     *國家“九五”重點科技攻關(guān)資助項目
     作者單位：林勝　門愛東　北京郵電大學(xué)電信工程學(xué)院，北京100876；第一作者25歲，男，博士生 
參考文獻(xiàn)
[1]　ISO/IEC 11172-3-1993　Coding of moving pictures and associated audio for digital storage media at up to about 1.5 Mbit/s——part 3: audio
[2]　Brandenburg K, Dehry Y F, Johnston J D, et al. ISO-MPEG-1 audi a generic standard for coding of high-quality digital audio. J Audio Eng Soc, 1994, 42 (10): 780～791
[3]　王建昕，董在望，印日方強(qiáng)MPEG音頻編碼算法的研究與實時實現(xiàn)清華大學(xué)學(xué)報， 1997， 37(10)： 45～48
[4]　Konstantinides K, Fast subband filtering in MPEG audio coding. IEEE Signal Processing Letters,　1994。