TW583575B - Digital audio sampling scheme - Google Patents

Digital audio sampling scheme Download PDF

Info

Publication number
TW583575B
TW583575B TW91135532A TW91135532A TW583575B TW 583575 B TW583575 B TW 583575B TW 91135532 A TW91135532 A TW 91135532A TW 91135532 A TW91135532 A TW 91135532A TW 583575 B TW583575 B TW 583575B
Authority
TW
Taiwan
Prior art keywords
digital audio
filter
mean square
root mean
patent application
Prior art date
Application number
TW91135532A
Other languages
English (en)
Other versions
TW200410126A (en
Inventor
Chuan Liu
Chih-Hsien Tsou
Yu-Chih Chin
Original Assignee
Silicon Integrated Sys Corp
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Silicon Integrated Sys Corp filed Critical Silicon Integrated Sys Corp
Priority to TW91135532A priority Critical patent/TW583575B/zh
Application granted granted Critical
Publication of TW583575B publication Critical patent/TW583575B/zh
Publication of TW200410126A publication Critical patent/TW200410126A/zh

Links

Landscapes

  • Complex Calculations (AREA)
  • Tone Control, Compression And Expansion, Limiting Amplitude (AREA)

Description

五、發明說明(1) 曼懇所屬之枯^ 本發明係有關於_ ^ -ψ a 加權最小均方根演算、去==音訊取樣的方,其使用 所需的係數。 展開要配置-SRC的脈衝響應時 先前技術 在音訊應用例之中,如第i圖所示,一數位 15:見在可將來自例如一音效卡u、一音源機12、一光了 放杰13及一卡匣播放器14的許多:播 混合,#蕤ha ^ 07汗夕曰訊,以不同取樣速率的 此口並猎此透過一後端電路16及左右各一揚聲器1?、 ’以享受更多彩多姿的數位式音m。轉換這些數位音 料至所需的取樣速率是混音器的主要功能。第2圖是一M、 赫茲取樣速率的波形圖。第3圖是一藉由不連續傅利葉轉 換(DFT)所產生的第i圖波形的頻率響應圖。第4圖是取樣 頻率為512赫茲時所產生的第1圖的頻率響應圖。第5圖是 一用以改變取樣速率(例如,從第2圖的32赫茲至第3圖的 512赫兹)的取樣速率轉換器(SRC)的方塊圖。在第5圖中, 貫務上’為了簡化計算起見,取樣速率轉換是在時間領域 上執行。如第5圖所示,在時間T時,輸入資料χκ經過一用 以在每一對樣本間插入零的内插裝置41、一用以執行迴旋 计箅(D F Τ)的低通滤波器及一用以產生上述迴旋計算所需 的時間刻度的定標器43,因而產生一輸出值γκ。為了在一 輸位式混音系統中安裝這樣的SRC功能,典型的方式是使 用一有限脈衝響應低通濾波器(此後,稱之為F IR濾波器) 。使用上述内插計算以取得在取樣點間的信號振幅資訊的
0702-7568TWF(nl);90P144;Sue.ptd 第5頁 583575 五、發明說明(2) ------- 超取樣法係為F I r濾波器設計領域中的主要興趣,尤其是 ^ ^ f於超取樣法中的脈衝響應。傳統上係利用Remez互 、演算法來達成。然而,在以最少計算次方來產生最佳化 IR濾波器的設計中,這個演算法既複雜又不容易安裝於 軟體計算中。
JiJ月内岑 有鑑於此,本發明之一目的係提供一種數位式音訊取 樣的方法’其使用加權最小均方根演算法以展開要配置一 SRC的脈衝響應時所需的係數。 根據本發明,本數位式音訊取樣的方法包含安裝一軟 體私式於一電腦内’用以藉由内置的加權最小均方根演算 法來計算一取樣速率轉換(SRC)濾波器的脈衝響應。藉此 ’整體發明可效地極小化軟體配置時所需的計算次方。 實施方式 全文中’類似元件功能以相同元件編號代表之。 第6圖顯示一本發明數位式混音器方塊圖。在第6圖中 ,所有輸入的數位式音訊資料1)1-1)11係在經過各自對應的 取樣速率轉換器(SRC)有限脈衝響應(FIR)低通濾波器61 (今後簡稱為FIR濾波器)後,由加法器62混音,藉以產生 一輸出波形。如第6圖所示,上述Jr I r濾波器β 1及加法器6 2 形成一混音器6 5。混音器6 5的配置係相同於習知的架構, 但疋在關鍵性F I R濾波器的配置方法上是使用加權最小均 方根演算法。使用加權最小均方根演算法所產生的本發明 F IR濾波器的頻率響應可以公式表示於下:
0702-7568TWF(η1);90P144;Sue.ptd 第 6 頁 583575 五、發明說明(3)
N εΐζ)= ς ρ-ζ^η *«〇 (1 ) 其中’ P(Z)是傅利葉轉換函式(此後稱之為z函式),係數 Pn是FIR渡波器的脈衝響應,其中n是FIR濾波器的長度(階 .. Λ 度)°令沒(Ζ)為想要的FIR濾波器頻率響應而頻率響應誤 差函式E為: ^(z) (2) E£z)= 2 P«z 方私式(2) 了 ·平畺於一線性分佈從ω = 〇到冗的頻率密度格 點上。對於一具有長度N的FIR濾波器而言,4Ν頻率格點就 足以表示濾波器的功效。若頻帶邊緣未落在所設的頻率格 點上時,就須要增加額外的格點以對應上述遺漏的頻帶邊 緣。此時,可以下列向量方程式來表示: E = Ua Λ Η 其中 (3) Ε [Ε(ζ»),Ε(ζζ)> ]' •1«,·.·,〆 (4) (5) (6) (7) a [p〇.…·ρ«Γ [舍(ZI),為(Z2),·· ]· 其中 Zi + 1 > Zi·
0702-7568TWF(n1);90P144;Sue.p t d 第7頁 583575 五、發明說明(4) 在本發明加權最小均方根方法的設計中, ( ) }被極小化,其中,rn是最小均方根加權值 (least square weighting vaiue)。最佳化的解決方案 使用下列方程式來表示: a = (UTRUrJU7R^ (8) 可表示成如下的一對角、線^車,只_角線 數值q、r2、···: q 〇 …〇 '」 (9) =線性相位低通濾波器的範例,其使用—指數函 入上述方程式(1)-(9)的線性相位項中,並與 罄ΪΪ較結果示於第7圖。上述比較基準為取樣頻率之‘ 帶邊緣在0.15、濾波器長度為51且所有對角線頻 件下。從第7圖中可看到,使用本發明最小均根方々條 配置的遽波器具有遠小於習知方式的漣波振盪音很量方式所 (ripple magnei tude)。又,使用本發明最小 所配置的濾波器在靠近頻帶邊緣處的效率,▲ 方式 緣處使用一較大的rn值時,相較其它地方可減=爭,帶邊 率損耗進而顯著地增加此處執行的效率。 夕夕的效 最小均方根技術的重點所在。 疋發明採用 如上述,本發明數位式音訊取樣的方法勺人… 體程式於一電腦内,用以藉由内置的 安裝〜軟 權最小均方根演算 0702-7568TWF(nl);90P144;Sue.ptd 第8頁 583575 五、發明說明(5) 法來計算一取樣速率轉換(SRC)濾波器的脈衝響應。該加 權最小均方根演算法也可配置在一數位信號處理器或一特 定積體電路内,藉以執行,並不限定於只經由上述電腦來 執行。如此,本發明可效地極小化軟體配置時所需的計算 次方。 雖然本發明已以較佳實施例揭露如上,然其並非用以 限定本發明,任何熟知此技術之人士,在不脫離本發明之 精神及範圍内,當可做更動與潤飾,因此本發明之保護範 圍當視後附之申請專利範圍所界定者為準。
0702-7568TWF(nl);90P144;Sue.ptd 第9頁 583575 圖式簡單說明 為讓本發明之上述及其它目的、特徵、與優點能更顯 而易見’下文特舉一較佳實施例,並配合所附圖式,作詳 細說明如下: 第1圖顯示一典型數位式音訊系統的方塊圖; 第2圖是一 3 2赫茲取樣速率的波形圖; 、第3圖是一藉由不連續傅利葉轉換(DFT)所產生的 圖波形的頻率響應圖; 響應ί4圖疋取樣頻率為512赫兹時所產生的第1圖的頻率 弟5圖是一用 的方塊圖; 以改變取樣速率的取樣速率 轉換器(SRC) 第6圖顯示一本發明數位式混音 繁Mr - σ方龙圖,及 第7圖顯不一使用習知方式與本發 方式進行最佳化所產生的濾波器 權最小均方根 [#㈣則 ❹羊響應的比較圖。 11 音效卡 12 音源機 13 光碟播放器 14 卡匣播放器 15 數位式混音器 16 後端電路 17 左揚聲器 18 右揚聲器 41 内插裝置
0702-7568TWF(nl);90P144;Sue.ptd 第10頁 583575 圖式簡單說明 42 低 通 濾 波 器 43 定 標 器 61 有 限 脈 衝 響應低通濾波器(F I R濾波器) 62 加 法 器 65 混 音 器 ϋϋ·ΙΙ 0702-7568TWF(nl);90P144;Sue.ptd 第 11 頁

Claims (1)

  1. 583575
    數位式音訊的方法,包含安裝-使用加權最 i鄕@ ’貝鼻法來計算一取樣速率轉換(SRC)濾波器的脈 衝a μ的軟體程式於一電腦内,其中,上述脈衝響應的相 關頻譜在一取樣頻率的倍數處具有波谷(n〇tches),可供 該SRC濾波器使用以進一步產生想要的頻率響應。 2 ·如申請專利範圍第1項之數位式音訊的方法,其中 該加權最小均方根演算法配置於一數位信號處理器。 3·如申請專利範圍第1項之數位式音訊的方法,其中 ’該加權最小均方根演算法配置於一專用的積體電路。 4·如申請專利範圍第1項之數位式音訊的方法,其中 ,該脈衝響應係以方程式a= (LTRLiruTRA 表示,其中,A 是想要的頻率響應,U是濾波函式,R是一其對角線是想要 的最小均方根加權值的對角矩陣,及T是反轉運算 (reverse operation) 〇 5 ·如申請專利範圍第4項之數位式音訊的方法,其 中 上述U、R及έ可以下列方程式表示: (7) 丨),為(ζ2),…Γ
    583575 六、申請專利範圍 — [uKW u = 1« …,z, L.....」 (5 ) f中W…為該最小均方根加權函式;Z為傅利 茶轉換函式,用以轉換時間領域上的外部輸入音訊传號 為頻率領域上所要的頻率響應,以及當i=1到 限脈衝響應(FIR)低通濾波器的長度時,Zi+1>zi。 6.—種具有一數位式音訊取樣方法的1音器,該數位 式音訊取樣方法包含一使用加權最小均方根演算法來計算 脈衝響應的裝置(means),該混音器包括: # _複數個並行的SRC濾波器,每個濾波器具有一或更多 脈衝響應,且每個濾波器連接至一外部不同音訊來源,以 接收該外部不同音訊來源的樣本並將該一或更多脈衝響應 與該樣本進行迴旋運算,以產生想要的輸出係數來形成想 要的頻率響應;及 一加法器,連接至該負數個並行的SRC濾波器,以結 合各輪出係數來產生一音訊輸出。 、、 如申請專利範圍第6項之具有一數位式音訊取樣方 法的'心音器’其中,該裝置(means)是一安裝著該加權最 小均方根演算法的數位信號處理器。 8 ·如申請專利範圍第6項之具有一數位式音訊取樣方 法的混音器,其中,該裝置(means)是一安裝著該加權最 小均方根演算法的專用積體電路。
TW91135532A 2002-12-09 2002-12-09 Digital audio sampling scheme TW583575B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
TW91135532A TW583575B (en) 2002-12-09 2002-12-09 Digital audio sampling scheme

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
TW91135532A TW583575B (en) 2002-12-09 2002-12-09 Digital audio sampling scheme

Publications (2)

Publication Number Publication Date
TW583575B true TW583575B (en) 2004-04-11
TW200410126A TW200410126A (en) 2004-06-16

Family

ID=34058043

Family Applications (1)

Application Number Title Priority Date Filing Date
TW91135532A TW583575B (en) 2002-12-09 2002-12-09 Digital audio sampling scheme

Country Status (1)

Country Link
TW (1) TW583575B (zh)

Also Published As

Publication number Publication date
TW200410126A (en) 2004-06-16

Similar Documents

Publication Publication Date Title
Pirkle Designing audio effect plugins in C++: for AAX, AU, and VST3 with DSP theory
JP3636361B2 (ja) 高精度及び高効率を有するディジタルフィルタ
EP2063413B1 (en) Reverberation effect adding device
CN104217714B (zh) 使用音频信号的固定分区大小卷积的声音合成
JP2001356788A (ja) 周波数補間装置、周波数補間方法及び記録媒体
Kahles et al. Oversampling for nonlinear waveshaping: Choosing the right filters
Bank Perceptually motivated audio equalization using fixed-pole parallel second-order filters
US9112479B2 (en) Filtering discrete time signals using a notch filter
JP2004294712A (ja) 残響音生成装置およびプログラム
EP2128853B1 (en) Reverberation imparting apparatus
TW583575B (en) Digital audio sampling scheme
Freed Real-time inverse transform additive synthesis for additive and pitch synchronous noise and sound spatialization
JP2779983B2 (ja) 電子楽器
JPS6217759B2 (zh)
JP3971193B2 (ja) 雑音発生装置、雑音発生方法およびプログラム
TWI254276B (en) Sound synthesizer and recording medium to record the program
CN1245683C (zh) 数字音讯取样的方法
JP5092902B2 (ja) Firフィルタ係数算出装置、firフィルタ装置、および、firフィルタ係数算出プログラム
Norilo Exploring the Vectored Time Variant Comb Filter.
JP3783552B2 (ja) 楽音信号合成方法、楽音信号合成装置および記録媒体
Yadav et al. Design of Modified RNS-PPA Based FIR Filter for High-Speed Application
JP2002271886A (ja) イコライザ及びそのフィルタ係数決定方法
JP2509009B2 (ja) 残響付加装置の係数設定方法
JPH04136992A (ja) 電子楽器の効果装置
Chassaing A senior project course in digital signal processing with the TMS320

Legal Events

Date Code Title Description
MM4A Annulment or lapse of patent due to non-payment of fees