MP3解碼的IMDCT硬件加速器方案

MPEG-l／2 Audio Layer 3(簡(jiǎn)稱(chēng)MP3)，是專(zhuān)門(mén)針對音樂(lè )和語(yǔ)音數據進(jìn)行設計的有損壓縮算法。在市場(chǎng)需求的推動(dòng)下，會(huì )有越來(lái)越多的嵌入式系統支持MP3應用。因此，如何在嵌入式系統中利用有限的計算能力實(shí)現MP3的實(shí)時(shí)解碼，成為值得關(guān)注的問(wèn)題。MP3解碼算法流程主要包括：幀同步和邊帶信息解碼、Huffman解壓縮、反量化、立體聲解碼、反鋸齒、IMDCT和子帶合成運算等。解碼軟件效率評估得到的數據表明，IMDCT過(guò)程的運算量占到了整個(gè)解碼運算總量的19％，由此設想在系統級芯片內部添加一小塊專(zhuān)用電路，專(zhuān)門(mén)負責處理MP3解碼過(guò)程中IMDCT部分的運算。我們把這一小塊專(zhuān)用電路稱(chēng)為“IMDCT硬件加速器”。將這部分代碼通過(guò)硬件實(shí)現，并以硬件速度執行這些運算，可以有效地提高系統的解碼性能。
本文采用軟硬件協(xié)同設計的方法，提出了一種針對IMDCT運算的硬件加速器方案，目的在于使MP3解碼過(guò)程中的速度和成本這兩大要素得到優(yōu)化。

1 IMDCT運算
MP3解碼算法中使用反向修正離散余弦變換IMDCT(Inverse Modified Discrete Cosine Transform)將輸入數據從頻域變換到余弦域，對子帶濾波進(jìn)行補償運算，如下式：

式中：在長(cháng)窗類(lèi)型幀中n取36；在短窗類(lèi)型幀中n取12。
完成IMDCT變換后，結果xi必須再與視窗函數Wi作乘積運算。視窗函數是由邊帶信息中bllk_type位的值來(lái)確定的。長(cháng)窗類(lèi)型叉可以根據幀首部定義進(jìn)一步劃分為aormal、start、stop三種子類(lèi)型，IMDCT變換得到36個(gè)結果數據；而在短窗類(lèi)型下，MP3解碼器將執行3次IM-DCT變換產(chǎn)生12個(gè)輸出結果，然后互相疊加補零后也得到36個(gè)數據，類(lèi)似于長(cháng)窗類(lèi)型的輸出結果。將這36個(gè)數據與上一次的結果進(jìn)行疊加，得到IMDCT變換的18個(gè)最終輸出結果。音頻數據分單聲道和雙聲道，包含若干顆粒；每個(gè)顆粒有576個(gè)數據項，共包括32個(gè)數據塊，需要分別進(jìn)行IMDCT變換。一個(gè)顆粒變換完成后，經(jīng)過(guò)較為簡(jiǎn)單的頻率反轉后得到的18個(gè)子帶(每個(gè)子帶包含32個(gè)數據)即可作為子帶合成的輸入信息。我們把從18項輸入數據到18項輸出數據的整個(gè)過(guò)程(包括IMDCT變換、數據加窗運算、疊加)稱(chēng)為“IMDCT運算”，如圖1所示。本文所討論的IMDCT硬件加速器，即實(shí)現這部分功能的專(zhuān)用電路。

2 IMDCT變換算法的選擇
ISO標準解碼代碼里，IMDCT變換的算法沒(méi)有進(jìn)行優(yōu)化。當式(1)中n=12時(shí)，做一次IMDCT需要72次乘法和66次加法；n=36時(shí)，做一次IMDCT需要648次乘法和630次加法?？梢?jiàn)，IMDCT變換占用了大量的CPU時(shí)間，成為MP3解碼過(guò)程中主要的性能瓶頸之一。將IMDCT做成硬件加速模塊，也需要使用更快速的IMDCT算法，以進(jìn)一步提高速度。這里引入一種新型的IMDCT算法。運用這種算法，當式(1)中n=12時(shí)，做一次IMDCT只需13次乘法和39次加法；n=36時(shí)，做一次IMDCT需要47次乘法和165次加法(詳見(jiàn)參考文獻)。同時(shí)，采用余弦查找表來(lái)代替實(shí)際的cos()函數運算，以加快長(cháng)窗下的cos 36和短窗下的cos 12的變換速度。IMDCT算法改進(jìn)后，運算過(guò)程得到簡(jiǎn)化，乘法數量大大減少，從而提高了系統性能。

3 疊加運算的優(yōu)化
速度和成本是設計的兩大要素?？紤]到硬件成本，IMDCT硬件單元的面積應受到嚴格控制。這里提出一種疊加運算的優(yōu)化算法，利用該算法，可以節省23118個(gè)字的存儲電路單元。下面對此優(yōu)化算法進(jìn)行詳細介紹。
IMDCT運算的主要數據結構如下：

可以看到，對于雙聲道立體聲數據(stereo=2)，常見(jiàn)的算法是對前一區塊2個(gè)顆粒(共23218個(gè)數據項)進(jìn)行IMDCT運算，保存23218個(gè)字的高18項數據，再對下一區塊的23218個(gè)數據項進(jìn)行IMDCT運算，將得到的23218個(gè)字的低18項數據與之前保存的上一區塊的高18項數據相疊加，得到輸出結果。將IMDCT運算的數據疊加部分用硬件實(shí)現時(shí)，需要一塊大小為23218個(gè)字的存儲電路(prevbuf[2][32][18])來(lái)存儲高18項數據，用于下一步的疊加運算。對嵌入式SoC來(lái)說(shuō)，降低存儲電路的需求意味著(zhù)減小芯片面積。于是采用軟硬件協(xié)同設計的方式對算法進(jìn)行改進(jìn)，通過(guò)改變IMDCT運算所需數據的輸入次序，相對減少I(mǎi)MDCT變換輸出數據的量。如此一來(lái)，有效地減少了硬件加速器的存儲單元，減小了電路面積。
首先，在解碼軟件里，完成2個(gè)區塊4個(gè)顆粒數據的反鋸齒運算，并將結果存儲到一片連續的內存區域內。具體做法如下：

然后，采用交錯讀取數據的方法，先將處在區塊相同位置、需要進(jìn)行疊加的子帶進(jìn)行IMDCT運算。這樣，可以將存儲中間結果的prevbuf縮小到218個(gè)字，大大降低了對存儲單元的需求，減小了電路面積。具體實(shí)現如下：

4 硬件實(shí)現
由圖l可知，IMDCT運算主要包含3部分：IMCCT變換、數據加窗運算以及疊加運算。IMDCT變換部分主要是把由反鋸齒運算得到的數據IN與cos系數相乘和累加，并把最后的結果放入寄存器SUM中；加窗運算是把計算好的SUM與加窗系數Wi相乘；疊加部分則是將加窗運算后的數據進(jìn)行疊加?？傮w結構圖如圖2所示。

由圖2可知，通過(guò)合理的流程設計和利用多路選擇器(MUX)，整個(gè)硬件加速器只需要1個(gè)乘法器和1個(gè)加法器，大大降低了通過(guò)硬件實(shí)現IMDCT運算的成本。

結 語(yǔ)
本文通過(guò)引進(jìn)新的IMDCT變換算法，優(yōu)化IMDCT運算過(guò)程中的疊加運算，加快了整體的運算速度，降低了對存儲單元的需求，為高速度、低成本地實(shí)現IMDCT硬件加速器提供了一種方案。