音响IC资料-EAD-DSP系列之算法.docVIP

細說EAD-DSP系列DSP之演算法 文/何志誠 在剛出現Over Sampling技術的那幾年，專業級的數位/類比轉換器的價位是那麼地高不可攀；最近幾年隨著IC製程技術的進步，如果一部幾千元的CD player連最簡單的四倍Over Sampling都沒有，似乎有點說不過去，除了CD/LD的解碼器之外，Over Sampling也是THX、AC-3、MPEG2的基本技術之一。老朽會拿EAD-DSP系列(為美國Enlightened Audio Designs公司所生產)來舉例說明乃是因為它們採用了日本NPC(SEIKO)公司所生產的DSP ASIC(Application Specific IC) SM5813(這一顆IC與Burr-Brown公司之DF-1700規格相同，目前有多家產品採用這一顆IC搭配Burr-Brown之PCM63PK之D/A converter。例如國產的雅瑟音響之DA20ST即採用了與EAD-DSP7000幾乎完全相同之關鍵零組件)，這一顆IC可以說是目前市面上所能看到之Over Sampling ASIC之代表作，也就是說它們都採用了相同的數學(演算法)架構。在這裡順便提一下，在準備本文之前，蒲總編曾經向我提及Crystal公司的CS-4303、CS-4328……一系列的IC，希望我能寫一系列的文章來談這一系列用於Multichannels Audio System與廣泛用途之Over Sampling與D/A converter IC。其實Crystal這一系列IC的演算法架構與NPC-SM5813並沒有甚麼不同，只是基於成本的考量與音質的需求不太高的情況下，才將OverSampling與D/A converter做在一起。當然，它們內部所用的數位濾波器長度也沒那麼長(這意味著計算結果的準確度亦較差)，而且其D/A converter也採用了較廉價的△-Σ(唸Delta-Sigma)架構，也就是俗稱的one-bit或bit-stream。至於NPC-SM5813或DF-1700沒有將D/A converter合併成一顆IC是為了Hi-End的考量，使其能搭配Burr-Brown的PCM63PK的緣故。 在談OverSampling的演算法之前，我們先假定看倌老爺您對Discrete Fourier Transform(DFT)或Fast Fourier Transform(FFT)、Convolution、Low pass filter與加權平均(Weighted Average method這個方法，老朽在高傳真視聽第252期「淺談超取樣」一文中曾提及)有個基本概念；至於詳細內容我們將在後面或另闢文章討論。首先，我們先來看一個例子。圖1a為函數：y = 1.5e-2t．sin(6πt)之連續波形。圖1b為圖1a之數位取樣的圖形，其取樣頻率為32Hz，也就是說圖1a波形之頻寬不高於16Hz。圖1c為圖1b經過富利葉轉換(Fourier Transform，其實是利用Discrete Fourier Transform算出來的)後之頻譜圖形，假如我們將圖1b中每兩個Sample間都「塞」入15個零，也就是說原來有33個Sample，現在則變成513個Sample，再將該513個取富利葉轉換，則可得到如圖1d之頻譜。我們試著比較圖1c與1d，我們可以發現圖1c之頻譜重複出現了16次(這時取樣頻率為512Hz)，相信這時候看倌老爺您已經了解，只要使用截止頻率為16Hz的低通濾波器就可以將圖1b變成圖1a的樣子；用相同的方法，我們也可以將圖2a“Up Sample”成圖2b。在這個例子裡我們可以用圖3a來描述“Up Sample”的過程，然而，在實際產品應用的情況(Real Time Case)卻不是這樣「硬幹」的，我們必須回到Time Domain來處理，除了方便製作，更可以降低成本，其過程如圖3b所示。接下來我們將深入探討如何實現該演算法，並將其簡化以使其更適合硬體製作與降低製作成本。 我們試著比較圖3a、3b與圖1a~1d的關係，很明顯的，在圖3a與3b中，Z=15，也就是“Up Sample”16倍，有一點我們必須注意，當我們在frequency domain(頻域)乘上一個低通濾波器，這就相當於在Time Domain(時域)與一個低通濾波器作Convolution的運算，先來看一下Discrete Convolution的定義。 (1) 習慣上，我們將上式表示為 (2) 我們舉一個例子來說明，圖4a與圖4b分別代表X(KT)與h(KT)，根據Convolution的定義(即方程式(1))知道要計算Convolution的結果y(KT)需要h[(0-i)T], h[(1-i)T], h[(