本申請(qǐng)是申請(qǐng)?zhí)枮?01380037816.5、申請(qǐng)日為2013年7月16日、發(fā)明名稱為“用于渲染音頻聲場(chǎng)表示以供音頻回放的方法和設(shè)備”的發(fā)明專利申請(qǐng)的分案申請(qǐng)。

本發(fā)明涉及一種用于渲染(render)音頻聲場(chǎng)表示(具體地，涉及高保真度立體聲響復(fù)制格式的音頻表示)以供音頻回放的方法和設(shè)備。

背景技術(shù)：

精確定位是任何空間音頻再現(xiàn)系統(tǒng)的關(guān)鍵目標(biāo)。這種再現(xiàn)系統(tǒng)高度可應(yīng)用于會(huì)議系統(tǒng)、游戲或受益于3d聲音的其他虛擬環(huán)境。3d中的聲音場(chǎng)景可被合成或捕獲為自然聲場(chǎng)。聲場(chǎng)信號(hào)(例如，高保真度立體聲響復(fù)制(ambisonics))攜帶期望聲場(chǎng)的表示。高保真度立體聲響復(fù)制格式基于聲場(chǎng)的球面諧波分解。雖然基本高保真度立體聲響復(fù)制格式或b格式使用階數(shù)0和1的球面諧波，所謂的高階高保真度立體聲響復(fù)制(hoa)還使用至少2階的其他球面諧波。需要解碼或渲染過(guò)程來(lái)從這種高保真度立體聲響復(fù)制格式的信號(hào)中獲得單獨(dú)的擴(kuò)音器信號(hào)。在本文中將擴(kuò)音器的空間布置稱為擴(kuò)音器設(shè)置。然而，雖然已知的渲染方案僅適于常規(guī)的擴(kuò)音器設(shè)置，任意的擴(kuò)音器設(shè)置更加普遍。如果將這種渲染方案應(yīng)用于任何擴(kuò)音器設(shè)置，聲音方向性受到損傷。

技術(shù)實(shí)現(xiàn)要素：

本發(fā)明描述了一種用于渲染/解碼針對(duì)于常規(guī)和非常規(guī)空間擴(kuò)音器分布的音頻聲場(chǎng)表示的方法，其中，該渲染/解碼提供了高度改進(jìn)的定位特性并節(jié)省了能量。具體地，本發(fā)明提供了用于獲得(例如hoa格式的)聲場(chǎng)數(shù)據(jù)的解碼矩陣的新方式。因?yàn)閔oa格式描述了不直接涉及擴(kuò)音器位置的聲場(chǎng)，且因?yàn)橐@得的擴(kuò)音器信號(hào)必然處于基于通道的音頻格式，hoa信號(hào)的解碼與音頻信號(hào)的渲染始終緊密相關(guān)。因此，本發(fā)明涉及解碼和渲染聲場(chǎng)相關(guān)的音頻格式。

本發(fā)明的一個(gè)優(yōu)點(diǎn)是實(shí)現(xiàn)節(jié)能解碼以及非常好的定向特性。術(shù)語(yǔ)“節(jié)能”指的是在解碼之后保留hoa定向信號(hào)中的能量，使得例如將以恒定的響度感知恒定幅度定向空間掃描。術(shù)語(yǔ)“好的定向特性”指的是通過(guò)定向的主波瓣和較小的副波瓣來(lái)表征的揚(yáng)聲器方向性，其中，與常規(guī)的渲染/解碼相比，提高了方向性。

本發(fā)明公開(kāi)了針對(duì)任意的擴(kuò)音器設(shè)置渲染聲場(chǎng)信號(hào)(例如，高階高保真度立體聲響復(fù)制(hoa))，其中，該渲染導(dǎo)致高度改進(jìn)的定位特性，且是節(jié)能的。這是通過(guò)用于聲場(chǎng)數(shù)據(jù)的新類(lèi)型的解碼矩陣以及用于獲得解碼矩陣的新方式來(lái)獲得的。在針對(duì)任意的空間擴(kuò)音器設(shè)置渲染音頻聲場(chǎng)表示的方法中，通過(guò)以下步驟獲得針對(duì)目標(biāo)擴(kuò)音器的給定布置進(jìn)行渲染的解碼矩陣：獲得目標(biāo)揚(yáng)聲器的數(shù)目及它們的位置、球面建模網(wǎng)格的位置以及hoa階數(shù)，根據(jù)建模網(wǎng)格的位置和揚(yáng)聲器的位置來(lái)生成混合矩陣，根據(jù)球面建模網(wǎng)格的位置和hoa階數(shù)來(lái)生成模式矩陣，根據(jù)混合矩陣和模式矩陣計(jì)算第一解碼矩陣，以及使用平滑和縮放系數(shù)來(lái)平滑和縮放第一解碼矩陣，以獲得節(jié)能的解碼矩陣。

在一個(gè)實(shí)施例中，如權(quán)利要求1中所述，本發(fā)明涉及一種用于解碼和/或渲染音頻聲場(chǎng)表示以供音頻回放的方法。在另一實(shí)施例中，如權(quán)利要求9中所述，本發(fā)明涉及一種用于解碼和/或渲染音頻聲場(chǎng)表示以供音頻回放的設(shè)備。在又一實(shí)施例中，如權(quán)利要求15中所述，本發(fā)明涉及其上存儲(chǔ)可執(zhí)行指令的計(jì)算機(jī)可讀介質(zhì)，所述可執(zhí)行指令使計(jì)算機(jī)執(zhí)行一種用于解碼和/或渲染音頻聲場(chǎng)表示以供音頻回放的方法。

一般而言，本發(fā)明使用以下方案。首先，導(dǎo)出取決于用于回放的擴(kuò)音器設(shè)置的平移(panning)函數(shù)。其次，針對(duì)擴(kuò)音器設(shè)置中的所有擴(kuò)音器，根據(jù)這些平移函數(shù)(或由平移函數(shù)所獲得的混合矩陣)計(jì)算解碼矩陣(例如，高保真度立體聲響復(fù)制解碼矩陣)。在第三步驟中，生成解碼矩陣并將其處理為是節(jié)能的。最后，對(duì)解碼矩陣濾波，以平滑擴(kuò)音器平移主波瓣并抑制副波瓣。針對(duì)給定的擴(kuò)音器設(shè)置，使用已濾波的解碼矩陣來(lái)渲染音頻信號(hào)。副波瓣是渲染的副作用，并在不想要的方向上提供音頻信號(hào)。因?yàn)獒槍?duì)給定的擴(kuò)音器設(shè)置來(lái)優(yōu)化渲染，副波瓣是令人煩惱的。本發(fā)明的優(yōu)點(diǎn)之一是最小化副波瓣，使得擴(kuò)音器信號(hào)的方向性得到改進(jìn)。

根據(jù)本發(fā)明一個(gè)實(shí)施例，一種用于解碼和/或渲染音頻聲場(chǎng)表示以供音頻回放的方法包括以下步驟：緩沖接收到的hoa時(shí)間采樣b(t)，其中，形成m個(gè)采樣的塊以及時(shí)間索引μ，對(duì)系數(shù)b(μ)進(jìn)行濾波以獲得已頻率濾波的系數(shù)使用解碼矩陣(d)將所述已頻率濾波的系數(shù)渲染(33)到空間域，其中，獲得空間信號(hào)w(μ)。在一個(gè)實(shí)施例中，其他步驟包括：在延遲線中針對(duì)所述l個(gè)通道中的每個(gè)通道單獨(dú)延遲所述時(shí)間采樣w(t)，其中，獲得l個(gè)數(shù)字信號(hào)，以及對(duì)所述l個(gè)數(shù)字信號(hào)進(jìn)行數(shù)模(d/a)變換和放大，其中，獲得l個(gè)模擬擴(kuò)音器信號(hào)。

用于渲染步驟(即，用來(lái)針對(duì)目標(biāo)揚(yáng)聲器的給定布置進(jìn)行渲染)的解碼矩陣d是通過(guò)以下步驟獲得的：獲得目標(biāo)揚(yáng)聲器的數(shù)目和揚(yáng)聲器的位置，確定球面建模網(wǎng)格的位置和hoa階數(shù)，根據(jù)球面建模網(wǎng)格的位置和揚(yáng)聲器的位置生成混合矩陣，根據(jù)球面建模網(wǎng)格的位置和hoa階數(shù)生成模式矩陣，根據(jù)混合矩陣g和模式矩陣計(jì)算第一解碼矩陣，以及利用平滑和縮放系數(shù)對(duì)第一解碼矩陣進(jìn)行平滑和縮放，其中，獲得解碼矩陣。

根據(jù)另一方面，一種用于解碼和/或渲染音頻聲場(chǎng)表示以供音頻回放的設(shè)備包括具有用于獲得解碼矩陣d的解碼矩陣計(jì)算單元的渲染處理單元，解碼矩陣計(jì)算單元包括：用于獲得目標(biāo)揚(yáng)聲器的數(shù)目l的裝置和用于獲得揚(yáng)聲器的位置的裝置；用于確定球面建模網(wǎng)格的位置的裝置和用于獲得hoa階數(shù)n的裝置；以及用于根據(jù)球面建模網(wǎng)格的位置和揚(yáng)聲器的位置生成混合矩陣g的第一處理單元；用于根據(jù)球面建模網(wǎng)格和hoa階數(shù)n生成模式矩陣的第二處理單元；用于根據(jù)執(zhí)行對(duì)模式矩陣與厄米特轉(zhuǎn)置混合矩陣g的乘積的緊致奇異值分解的第三處理單元(其中，u、v是根據(jù)酉矩陣導(dǎo)出的，以及s是具有奇異值項(xiàng)的對(duì)角矩陣)；用于根據(jù)來(lái)由u、v矩陣計(jì)算第一解碼矩陣的計(jì)算裝置，其中是單位矩陣或?qū)蔷仃嚕鰧?duì)角矩陣是根據(jù)所述具有奇異值項(xiàng)的對(duì)角矩陣導(dǎo)出的；以及用于利用平滑系數(shù)對(duì)第一解碼矩陣進(jìn)行平滑和縮放的平滑及縮放單元，其中，獲得解碼矩陣d。

根據(jù)又一方面，一種計(jì)算機(jī)可讀介質(zhì)在其上存儲(chǔ)有可執(zhí)行指令，當(dāng)在計(jì)算機(jī)上執(zhí)行該可執(zhí)行指令時(shí)，該可執(zhí)行指令使計(jì)算機(jī)執(zhí)行上述用于解碼音頻聲場(chǎng)表示以供音頻回放的方法。

當(dāng)與附圖結(jié)合時(shí)，根據(jù)對(duì)以下描述和所附權(quán)利要求的考慮，本發(fā)明的其他目標(biāo)、特征和優(yōu)點(diǎn)將變得顯而易見(jiàn)。

附圖說(shuō)明

參考附圖對(duì)本發(fā)明的示例性實(shí)施方式進(jìn)行描述，在附圖中：

圖1是根據(jù)本發(fā)明的一個(gè)實(shí)施例的方法的流程圖；

圖2是用于構(gòu)建混合矩陣g的方法的流程圖；

圖3是渲染器的框圖；

圖4是解碼矩陣生成過(guò)程的示意性步驟的流程圖；

圖5是解碼矩陣生成單元的框圖；

圖6是示例性的16個(gè)揚(yáng)聲器的設(shè)置，其中，揚(yáng)聲器被示出為連接的節(jié)點(diǎn)；

圖7是自然視角下的示例性的16個(gè)揚(yáng)聲器的設(shè)置，其中，節(jié)點(diǎn)被示出為揚(yáng)聲器；

圖8是示出比率的能量圖，該比率針對(duì)于利用現(xiàn)有技術(shù)[14]獲得的解碼矩陣的完美節(jié)能特征是恒定的，其中，n＝3；

圖9是針對(duì)根據(jù)現(xiàn)有技術(shù)[14](n＝3)設(shè)計(jì)的解碼矩陣的聲壓圖，其中，中心揚(yáng)聲器的平移(panning)波束具有強(qiáng)的副波瓣；

圖10是示出比率的能量圖，該比率的波動(dòng)比利用現(xiàn)有技術(shù)[2]獲得的解碼矩陣的4db大，其中，n＝3；

圖11是針對(duì)根據(jù)現(xiàn)有技術(shù)[2](n＝3)設(shè)計(jì)的解碼矩陣的聲壓圖，其中，中心揚(yáng)聲器的平移波束具有較小的副波瓣；

圖12是示出比率的能量圖，該比率的波動(dòng)比通過(guò)根據(jù)本發(fā)明的方法或裝置獲得的1db小，其中，利用相等的響度感知到具有恒定幅度的空間平移；

圖13是針對(duì)利用根據(jù)本發(fā)明的方法設(shè)計(jì)的解碼矩陣的聲壓圖，其中，中心揚(yáng)聲器具有帶有較小副波瓣的平移波束。

具體實(shí)施方式

大體上，本發(fā)明涉及向擴(kuò)音器渲染(即，解碼)聲場(chǎng)格式的音頻信號(hào)(例如，高階高保真度立體聲響復(fù)制(hoa)音頻信號(hào))，其中，擴(kuò)音器位于對(duì)稱或不對(duì)稱、常規(guī)或非常規(guī)位置。音頻信號(hào)可適于饋送比可用擴(kuò)音器更多的擴(kuò)音器，例如，hoa系數(shù)的數(shù)目可以大于擴(kuò)音器的數(shù)目。本發(fā)明以非常好的定向特性為解碼器提供了節(jié)能解碼矩陣，即，與利用常規(guī)解碼矩陣所獲得的揚(yáng)聲器方向性波瓣相比，揚(yáng)聲器方向性波瓣大體上包括更強(qiáng)的定向主波瓣和較小的副波瓣。節(jié)能指的是在解碼之后保留hoa方向性信號(hào)中的能量，使得例如以恒定的響度感知恒定幅度定向空間掃描。

圖1輸出了根據(jù)本發(fā)明的一個(gè)實(shí)施例的方法的流程圖。在該實(shí)施例中，用于渲染(即，解碼)hoa音頻聲場(chǎng)表示以供音頻回放的方法使用如下生成的解碼矩陣：首先，確定11目標(biāo)擴(kuò)音器的數(shù)目l、擴(kuò)音器的位置球面建模網(wǎng)格和階數(shù)n(例如，hoa階數(shù))。根據(jù)揚(yáng)聲器的位置和球面建模網(wǎng)格生成12混合矩陣g，以及根據(jù)球面建模網(wǎng)格和hoa階數(shù)n，生成13模式矩陣根據(jù)混合矩陣g和模式矩陣計(jì)算14第一解碼矩陣利用平滑系數(shù)平滑15第一解碼矩陣其中，獲得已平滑解碼矩陣以及利用根據(jù)已平滑解碼矩陣d獲得的縮放因子來(lái)縮放16已平滑解碼矩陣其中，獲得解碼矩陣d。在一個(gè)實(shí)施例中，平滑15和縮放16在單個(gè)步驟中執(zhí)行。

在一個(gè)實(shí)施例中，通過(guò)兩種不同方法中的一種來(lái)獲得平滑系數(shù)這取決于擴(kuò)音器的數(shù)目l和hoa系數(shù)通道的數(shù)目o3d＝(n+1)²。如果擴(kuò)音器的數(shù)目l低于hoa系數(shù)通道的數(shù)目o3d，使用用于獲得平滑系數(shù)的新方法。

在一個(gè)實(shí)施例中，生成對(duì)應(yīng)于多個(gè)不同擴(kuò)音器布置的多個(gè)解碼矩陣，并將其存儲(chǔ)以供后續(xù)使用。不同的擴(kuò)音器布置可以在以下至少一方面不同：擴(kuò)音器的數(shù)目、一個(gè)或多個(gè)擴(kuò)音器的位置、以及輸入音頻信號(hào)的階數(shù)n。因此，在初始化渲染系統(tǒng)時(shí)，確定匹配解碼矩陣，根據(jù)當(dāng)前需要從存儲(chǔ)器獲取該匹配解碼矩陣，并將其用于解碼。

在一個(gè)實(shí)施例中，通過(guò)根據(jù)執(zhí)行對(duì)模式矩陣與厄米特轉(zhuǎn)置混合矩陣g^h的乘積的緊致奇異值分解(compactsingularvaluedecomposition)，并根據(jù)由矩陣u、v計(jì)算第一解碼矩陣獲得解碼矩陣d。u、v是根據(jù)酉矩陣導(dǎo)出的，以及s是具有模式矩陣與厄米特轉(zhuǎn)置混合矩陣g^h的乘積的緊致奇異值分解的奇異值元素的對(duì)角矩陣。根據(jù)本實(shí)施例獲得的解碼矩陣通常在數(shù)值上比利用下述備選實(shí)施例獲得的解碼矩陣更加穩(wěn)定。矩陣的厄米特轉(zhuǎn)置是該矩陣的復(fù)共軛轉(zhuǎn)置。

在備選實(shí)施例中，通過(guò)根據(jù)執(zhí)行對(duì)厄米特轉(zhuǎn)置模式矩陣與混合矩陣g的乘積的緊致奇異值分解，獲得解碼矩陣d，其中，通過(guò)導(dǎo)出第一解碼矩陣。

在一個(gè)實(shí)施例中，根據(jù)對(duì)模式矩陣和混合矩陣g執(zhí)行緊致奇異值分解，其中，通過(guò)導(dǎo)出第一解碼矩陣，其中，是通過(guò)將大于等于閾值thr的所有奇異值替換為1并將小于該閾值thr的元素替換為0，根據(jù)奇異值分解矩陣s導(dǎo)出的截?cái)嗑o致奇異值分解矩陣。閾值thr取決于奇異值分解矩陣的實(shí)際值，且示例性地，可以在0.06*s1(s的最大元素)的量級(jí)上。

在一個(gè)實(shí)施例中，根據(jù)對(duì)模式矩陣和混合矩陣g執(zhí)行緊致奇異值分解，其中，通過(guò)導(dǎo)出第一解碼矩陣。和閾值thr如上針對(duì)之前實(shí)施例所述。閾值thr通常是根據(jù)最大奇異值導(dǎo)出的。

在一個(gè)實(shí)施例中，根據(jù)hoa階數(shù)n和目標(biāo)揚(yáng)聲器的數(shù)目l，使用兩種不同的方法來(lái)計(jì)算平滑系數(shù)：如果存在比hoa通道少的目標(biāo)揚(yáng)聲器，即，如果o3d＝(n²+1)＞l，則平滑和縮放系數(shù)對(duì)應(yīng)于常規(guī)的maxre系數(shù)集合，常規(guī)的maxre系數(shù)集合是根據(jù)階數(shù)n+1的勒讓德多項(xiàng)式的零導(dǎo)出的；否則如果存在足夠的目標(biāo)揚(yáng)聲器，即，如果o3d＝(n²+1)≤l，則根據(jù)由長(zhǎng)度等于(2n+1)且寬帶等于2n的凱撒窗的元素來(lái)構(gòu)建系數(shù)其中，縮放因子為cf。所使用的凱撒窗的元素開(kāi)始于僅被使用一次的第(n+1)個(gè)元素，并以被重復(fù)使用的后續(xù)元素繼續(xù)：第(n+2)個(gè)元素被使用3次，等等。

在一個(gè)實(shí)施例中，根據(jù)已平滑的解碼矩陣來(lái)獲得縮放因子。具體地，在一個(gè)實(shí)施例中，根據(jù)以下等式獲得縮放因子

下面描述完整的渲染系統(tǒng)。本發(fā)明的主要關(guān)注點(diǎn)是渲染器的初始化階段，在該階段如上所述地生成解碼矩陣d。在此，主要的關(guān)注點(diǎn)是用于導(dǎo)出一個(gè)或多個(gè)解碼矩陣(例如，用于碼本)的技術(shù)。為了生成解碼矩陣，有多少目標(biāo)擴(kuò)音器可用以及它們位于何處(即，它們的位置)是已知的。

圖2示出了根據(jù)本發(fā)明的一個(gè)實(shí)施例的用于構(gòu)建混合矩陣g的方法的流程圖。在該實(shí)施例中，創(chuàng)建21僅具有零的初始混合矩陣，以及針對(duì)每個(gè)具有角方向ωs＝[θs，φs]^t和半徑rs的虛擬源，執(zhí)行下面的步驟。首先，確定22環(huán)繞位置的三個(gè)擴(kuò)音器l1、l2、l3，其中，采用單位半徑，且構(gòu)建23矩陣其中根據(jù)lt＝spherical_to_cartesian(r)，將矩陣r變換24為笛卡爾坐標(biāo)。然后，根據(jù)s＝(sinθscosφs，sinθssinφs，cosθs)^t構(gòu)建25虛擬源位置，且根據(jù)g＝lt^-1s計(jì)算26增益g，其中，根據(jù)g＝g/||g||2來(lái)歸一化27增益，且將g的對(duì)應(yīng)元素gl，s替換為歸一化的增益：

下面一節(jié)給出對(duì)高階高保真度立體聲響復(fù)制(hoa)的簡(jiǎn)要介紹，并定義要處理(即，針對(duì)擴(kuò)音器渲染)的信號(hào)。

高階高保真度立體聲響復(fù)制(hoa)基于對(duì)感興趣的緊致區(qū)域內(nèi)的聲場(chǎng)的描述，該感興趣的緊致區(qū)域被假定為與聲源無(wú)關(guān)。在該情況下，在感興趣的區(qū)域內(nèi)，在時(shí)間t以及位置x＝[r，θ，φ]^t(球面坐標(biāo)：半徑r，傾角θ，方位角φ)處的聲壓p(t,x)的時(shí)空行為在物理上完全通過(guò)均勻波(homogeneouswave)等式來(lái)決定。可以示出，可根據(jù)[13]將聲壓相對(duì)于時(shí)間的傅里葉變換(即，其中，ω表示角頻率，且對(duì)應(yīng)于)擴(kuò)展到球面諧波(sh)序列中：

在等式(2)中，cs表示聲音的速度，以及是角波數(shù)(angularwavenumber)。此外，jn(·)指示第一類(lèi)的且階數(shù)為n的球面貝塞爾函數(shù)，以及表示階數(shù)為n且度數(shù)為m的球面諧波(sh)。與聲場(chǎng)有關(guān)的完整信息實(shí)際上包含在聲場(chǎng)系數(shù)內(nèi)。

應(yīng)該注意到，sh一般是復(fù)數(shù)取值的函數(shù)。然而，通過(guò)它們的適當(dāng)線性組合，可以獲得取實(shí)數(shù)值的函數(shù)，并關(guān)于這些函數(shù)執(zhí)行擴(kuò)展。

關(guān)于等式(2)中的壓力聲場(chǎng)，可以將源場(chǎng)定義為：

其中，源場(chǎng)或幅度密度[12]d(kcs，ω)取決于角波數(shù)和角度方向ω＝[θ，φ]^t。源場(chǎng)可以由遠(yuǎn)場(chǎng)/近場(chǎng)的離散/連續(xù)源組成[1]。通過(guò)以下等式，源場(chǎng)系數(shù)與聲場(chǎng)系數(shù)有關(guān)[1]：

其中，是第二類(lèi)的球面漢克爾函數(shù)，以及rs是相對(duì)于原點(diǎn)的源距離。

可以在頻域或時(shí)域中將hoa域中的信號(hào)表示為源場(chǎng)或聲場(chǎng)系數(shù)的傅里葉逆變換。下面的描述將假設(shè)使用有限數(shù)目的源場(chǎng)系數(shù)的時(shí)域表示：

：在n＝n處截?cái)嗟仁?3)中的無(wú)限序列。截?cái)鄬?duì)應(yīng)于空間帶寬限制。

系數(shù)(或hoa通道)的數(shù)目被給出為：

o3d＝(n+1)²對(duì)于3d(6)

或?qū)τ趦H2d的描述，給出為o2d＝2n+1。系數(shù)包括用于后續(xù)擴(kuò)音器的再現(xiàn)的一個(gè)時(shí)間采樣t處的音頻信息。它們可被存儲(chǔ)或發(fā)送，并因此經(jīng)受到數(shù)據(jù)速率壓縮?？赏ㄟ^(guò)具有o3d個(gè)元素的矢量b(t)來(lái)表示系數(shù)的單個(gè)時(shí)間采樣t：

以及通過(guò)矩陣來(lái)表示m個(gè)時(shí)間采樣的塊

b：＝[b(tstart+1)，b(tstart+2)，..，b(tstart+m)](8)

可通過(guò)利用環(huán)形諧波(circularharmonic)的擴(kuò)展來(lái)導(dǎo)出聲場(chǎng)的二維表示。這是上述的一般性描述的特殊情況，該特殊情況使用固定傾角不同的系數(shù)加權(quán)以及精簡(jiǎn)至o2d個(gè)系數(shù)的集合(m＝±n)。因此，所有下面的考慮也適用于2d表示；術(shù)語(yǔ)“球面”從而需要替換為術(shù)語(yǔ)“環(huán)形”。

在一個(gè)實(shí)施例中，將元數(shù)據(jù)與系數(shù)數(shù)據(jù)一起發(fā)送，允許明確地識(shí)別系數(shù)數(shù)據(jù)。通過(guò)所發(fā)送的元數(shù)據(jù)或因?yàn)榻o定的上下文，給出用于導(dǎo)出時(shí)間采樣系數(shù)矢量b(t)的所有必要信息。此外，要注意到的是，hoa階數(shù)n或o3d中的至少一個(gè)，且在一個(gè)實(shí)施例中還包括特殊標(biāo)記以及用于指示近場(chǎng)記錄的rs，在解碼器處是已知的。接下來(lái)描述向擴(kuò)音器渲染hoa信號(hào)。該節(jié)示出了解碼以及某些數(shù)學(xué)特性的基本原理。

基本解碼假設(shè)：第一，平面波擴(kuò)音器信號(hào)，以及第二，可以忽略從揚(yáng)聲器到原點(diǎn)的距離。可以將針對(duì)位于球面方向(l＝1，...，l)處的l個(gè)擴(kuò)音器所渲染的hoa系數(shù)b的時(shí)間采樣描述為[10]：

w＝db(9)

其中，表示l個(gè)揚(yáng)聲器信號(hào)的時(shí)間采樣，以及解碼矩陣可通過(guò)以下等式導(dǎo)出解碼矩陣

d＝ψ⁺(10)

其中，ψ⁺是模式矩陣ψ的偽取逆。模式矩陣ψ被定義為

ψ＝[y1，…yl](11)

其中，以及由揚(yáng)聲器方向的球面諧波組成，其中，h表示復(fù)共軛轉(zhuǎn)置(也稱為厄米特)。

接下來(lái)，描述由奇異值分解(svd)進(jìn)行的矩陣的偽取逆。導(dǎo)出偽取逆的一種通用方式是首先計(jì)算緊致svd：

ψ＝usv^h(12)

其中，是根據(jù)旋轉(zhuǎn)矩陣導(dǎo)出的，以及是降序排列的奇異值s1≥s2≥…≥sk的對(duì)角矩陣，其中，k＞0以及k≤min(o3d，l)。

通過(guò)以下等式確定偽取逆：

其中，對(duì)于sk具有非常小的值的不好的條件矩陣，將對(duì)應(yīng)的逆值替換為0。這被稱為截?cái)嗥娈愔捣纸狻Ｍǔ?，選擇關(guān)于最大奇異值s1的檢測(cè)閾值來(lái)識(shí)別要被替換為0的對(duì)應(yīng)逆值。

下面描述節(jié)能特性。通過(guò)以下等式給出hoa域中的信號(hào)能量：

e＝b^hb(14)

以及通過(guò)以下等式給出空間域中的對(duì)應(yīng)能量：

節(jié)能解碼器矩陣的比率是(基本)恒定的。這僅在d^hd＝ci時(shí)可實(shí)現(xiàn)，其中，單位矩陣為i，以及常數(shù)這要求d的norm-2(范2)條件數(shù)cond(d)＝1。再次地，這要求d的svd(奇異值分解)產(chǎn)生相同的奇異值：d＝usv^h，其中，s＝diag(sk，...，sk)。

一般而言，節(jié)能渲染器設(shè)計(jì)是本領(lǐng)域已知的。在[14]中通過(guò)以下等式提出了針對(duì)l≥o3d的節(jié)能解碼器矩陣設(shè)計(jì)：

d＝vu^h(16)

其中，將來(lái)自等式(13)的強(qiáng)制為并因此可以在等式(16)中丟棄乘積d^hd＝uv^hvu^h＝i，以及比率變?yōu)?。該設(shè)計(jì)方法的好處是能量節(jié)省，該能量節(jié)省保證了同質(zhì)空間聲印象的，其中，空間平移在感知到的響度上沒(méi)有波動(dòng)。該設(shè)計(jì)的缺陷是：針對(duì)不對(duì)稱的、非常規(guī)的擴(kuò)音器位置(參見(jiàn)圖8-9)，方向精確度的丟失以及較強(qiáng)的擴(kuò)音器波束側(cè)波瓣。本發(fā)明可以克服該缺陷。

針對(duì)非常規(guī)位置的揚(yáng)聲器的渲染器設(shè)計(jì)也是本領(lǐng)域已知的。在[2]中描述了針對(duì)l≥o3d和l＜o(jì)3d的解碼器設(shè)計(jì)方法，該解碼器設(shè)計(jì)方法允許在再現(xiàn)方向性上以較高精確度進(jìn)行渲染。該設(shè)計(jì)方法的缺陷是所導(dǎo)出的渲染器不節(jié)能(參見(jiàn)圖10-11)。

可以將球面卷積用于空間平滑。這是空間濾波過(guò)程，或是系數(shù)域中的窗口化(windowing)(卷積)。其目的是最小化副波瓣，稱為平移波瓣。通過(guò)原始hoa系數(shù)與帶狀系數(shù)(zonalcoefficient)的加權(quán)乘積來(lái)給出新的系數(shù)[5]：

這等效于在空間域?qū)²的左卷積[5]。在[5]中，將這方便地用于在通過(guò)對(duì)hoa系數(shù)b加權(quán)來(lái)進(jìn)行渲染/解碼之前對(duì)擴(kuò)音器信號(hào)的定向特性進(jìn)行平滑，該加權(quán)通過(guò)以下等式進(jìn)行：

其中，矢量通常包含取實(shí)數(shù)值的加權(quán)系數(shù)和常數(shù)因子df。平滑的概念是隨著增加的階數(shù)索引n對(duì)hoa系數(shù)進(jìn)行衰減。平滑加權(quán)系數(shù)的已知示例是所謂的maxrv、和maxre以及同相系數(shù)[4]。第一項(xiàng)提供缺省的幅度波束(不重要的(trivial)，長(zhǎng)度為o3d的全1矢量)，第二項(xiàng)提供均勻分布的角功率以及同相特征全副波瓣抑制。

下面描述所公開(kāi)解決方案的其他細(xì)節(jié)和實(shí)施例。首先，在初始化、啟動(dòng)行為和處理方面對(duì)渲染器架構(gòu)進(jìn)行描述。

每次擴(kuò)音器設(shè)置(即，擴(kuò)音器的數(shù)目以及任何擴(kuò)音器相對(duì)于收聽(tīng)位置的位置發(fā)生改變)，渲染器需要執(zhí)行初始化過(guò)程，以確定針對(duì)所支持hoa輸入信號(hào)具有的任何hoa階數(shù)的解碼矩陣組。同樣地，根據(jù)揚(yáng)聲器與收聽(tīng)位置之間的距離來(lái)確定延遲線的單獨(dú)揚(yáng)聲器延遲dl和揚(yáng)聲器增益下面描述該過(guò)程。在一個(gè)實(shí)施例中，所導(dǎo)出的解碼矩陣存儲(chǔ)在碼本內(nèi)。每次hoa音頻輸入特征改變，渲染器控制單元確定當(dāng)前有效的特征，并從碼本選擇匹配的解碼矩陣。碼本關(guān)鍵字可以是hoa階數(shù)n，或等效地，o3d(參見(jiàn)等式(6))。

參考圖3解釋用于渲染的數(shù)據(jù)處理的示意性步驟，圖3示出了渲染器的處理框的框圖。它們是第一緩沖器31、頻域?yàn)V波單元32、渲染處理單元33、第二緩沖器34、用于l個(gè)通道的延遲單元35、以及數(shù)模變換器及放大器36。

首先在第一緩沖器31中存儲(chǔ)具有時(shí)間索引t和o3dhoa系數(shù)通道的hoa時(shí)間采樣b(t)，以形成塊索引為μ的m個(gè)采樣的塊。在頻域?yàn)V波單元32中對(duì)的系數(shù)b(μ)進(jìn)行頻率濾波，以獲得已頻率濾波的塊該技術(shù)已知(參見(jiàn)[3])用于補(bǔ)償球面擴(kuò)音器源的距離，并用于使得可處理近場(chǎng)記錄。在渲染處理單元33中通過(guò)以下等式來(lái)向空間域渲染已頻率濾波的塊

其中，表示具有m個(gè)時(shí)間采樣的塊的l個(gè)通道中的空間信號(hào)。該信號(hào)在第二緩沖器34中緩沖，并被串行化，以形成在l個(gè)通道中具有時(shí)間索引t的單個(gè)時(shí)間采樣，在圖3中稱之為w(t)。這是饋送到延遲單元35中的l個(gè)數(shù)字延遲線的串行信號(hào)。延遲線補(bǔ)償收聽(tīng)位置到延遲為dl個(gè)采樣的單獨(dú)的揚(yáng)聲器l之間的不同距離。理論上，每個(gè)延遲線是fifo(先進(jìn)先出存儲(chǔ)器)。然后，在數(shù)模變換器及放大器36中對(duì)已延遲補(bǔ)償?shù)男盘?hào)355進(jìn)行d/a變換并放大，數(shù)模變換器及放大器36提供可饋送到l個(gè)擴(kuò)音器的信號(hào)365?？梢栽赿/a變換之前或通過(guò)在模擬域中采用揚(yáng)聲器通道放大來(lái)考慮揚(yáng)聲器增益補(bǔ)償

渲染器初始化如下進(jìn)行。

首先，需要知道揚(yáng)聲器的數(shù)目和位置。初始化的第一步驟是使新的揚(yáng)聲器數(shù)目l及相關(guān)位置可用，其中，其中，rl是從收聽(tīng)位置到揚(yáng)聲器l的距離，以及和是相關(guān)球面角?？蓱?yīng)用各種方法，例如，揚(yáng)聲器位置的手動(dòng)輸入，或使用測(cè)試信號(hào)的自動(dòng)初始化?？墒褂眠m當(dāng)?shù)慕涌?如，已連接的移動(dòng)設(shè)備或集成于設(shè)備的用于選擇預(yù)定義位置集合的用戶界面)來(lái)進(jìn)行揚(yáng)聲器位置的手動(dòng)輸入。可利用評(píng)估單元，使用麥克風(fēng)陣列和專用的揚(yáng)聲器測(cè)試信號(hào)來(lái)進(jìn)行自動(dòng)初始化，以導(dǎo)出通過(guò)rmax＝max(r1，...，rl)確定最大距離rmax，通過(guò)rmin＝min(r1，...，rl)確定最小距離rmin。

將l個(gè)距離rl和rmax輸入到延遲線和增益補(bǔ)償35。通過(guò)以下等式來(lái)確定針對(duì)每個(gè)揚(yáng)聲器通道dl的延遲采樣的數(shù)目：

其中，采樣速率為fs，聲音速度為c(溫度為20攝氏度時(shí)，)，以及指示向下一個(gè)整數(shù)取整。為了補(bǔ)償針對(duì)不同rl的揚(yáng)聲器增益，通過(guò)確定擴(kuò)音器增益或使用聲學(xué)測(cè)量來(lái)導(dǎo)出擴(kuò)音器增益

如下進(jìn)行(例如，針對(duì)于碼本的)解碼矩陣的計(jì)算。圖4示出了在一個(gè)實(shí)施例中用于生成解碼矩陣的方法的示意性步驟。圖5示出了在一個(gè)實(shí)施例中用于生成解碼矩陣的對(duì)應(yīng)設(shè)備的處理框。輸入是揚(yáng)聲器方向球面建模網(wǎng)格和hoa階數(shù)n。

可以將揚(yáng)聲器方向表達(dá)為球面角以及通過(guò)球面角ωs＝[θs，φs]^t來(lái)表達(dá)球面建模網(wǎng)格將方向的數(shù)目選擇為大于揚(yáng)聲器的數(shù)目(s＞l)，且大于hoa系數(shù)的數(shù)目(s＞o3d)。網(wǎng)格的方向應(yīng)該通過(guò)非常規(guī)則的方式來(lái)采樣單位球面。在[6]、[9]中討論了適合的網(wǎng)格，且可在[7]、[8]中找到適合的網(wǎng)格。一次性選擇網(wǎng)格作為示例，根據(jù)[6],s＝324個(gè)網(wǎng)格足以用于解碼最多hoa階數(shù)n＝9的矩陣。針對(duì)不同的hoa階數(shù)，可以使用其他網(wǎng)格。遞增地選擇hoa階數(shù)n，以根據(jù)n＝1，...，nmax填充碼本，其中，nmax是所支持的hoa輸入內(nèi)容的最大hoa階數(shù)。

將揚(yáng)聲器方向和球面建模網(wǎng)格輸入到構(gòu)建混合矩陣框41，構(gòu)建混合矩陣框41生成其混合矩陣g。將球面建模網(wǎng)格和hoa階數(shù)n輸入到構(gòu)建模式矩陣框42，構(gòu)建模式矩陣框42生成其模式矩陣將混合矩陣g和模式矩陣輸入到構(gòu)建解碼矩陣框43，構(gòu)建解碼矩陣框43生成其解碼矩陣將該解碼矩陣輸入到平滑解碼矩陣框44，平滑解碼矩陣框44平滑并縮放解碼矩陣。下面提供其他細(xì)節(jié)。平滑解碼矩陣框44的輸出是解碼矩陣d，利用相關(guān)的關(guān)鍵字n(或備選地o3d)將解碼矩陣d存儲(chǔ)在碼本中。在構(gòu)建模式矩陣框42中，球面建模網(wǎng)格被用于構(gòu)建類(lèi)似于等式(11)的模式矩陣：其中，要注意到的是，在[2]中將模式矩陣稱為ξ。

在構(gòu)建混合矩陣框41中，利用來(lái)創(chuàng)建混合矩陣g。要注意到的是，在[2]中將混合矩陣g稱為w?；旌暇仃噂的第l行由從方向到揚(yáng)聲器l的針對(duì)混合s虛擬源的混合增益組成。在一個(gè)實(shí)施例中，矢量基幅度平移(vbap)[11]被用于導(dǎo)出這些混合增益，[2]中也是如此。用于導(dǎo)出g的算法總結(jié)如下：

1利用0值來(lái)創(chuàng)建g(即，初始化g)

2針對(duì)每個(gè)s＝1...s

3{

4找到圍繞位置的3個(gè)揚(yáng)聲器l1，l2，l3，假設(shè)單位半徑以及構(gòu)建矩陣其中，

5在笛卡爾坐標(biāo)中計(jì)算lt＝spherical_to_cartesian(r)。

6構(gòu)建虛擬源位置s＝(sinθscosφs，sinθssinφs，cosθs)^t。

7計(jì)算g＝lt^-1s，其中

8歸一化增益：g＝g/||g||2

9利用g的元素來(lái)填充g的相關(guān)元素gl，s：

10}

在構(gòu)建解碼矩陣框43中，計(jì)算對(duì)模式矩陣與轉(zhuǎn)置混合矩陣的矩陣乘積的緊致奇異值分解。這是本發(fā)明的一個(gè)重要方面，可通過(guò)各種方式來(lái)執(zhí)行。在一個(gè)實(shí)施例中，根據(jù)以下等式來(lái)計(jì)算模式矩陣與轉(zhuǎn)置混合矩陣g^t的矩陣乘積的緊致奇異值分解s：

在備選實(shí)施例中，根據(jù)以下等式來(lái)計(jì)算模式矩陣與偽逆混合矩陣g⁺的矩陣乘積的緊致奇異值分解s：

其中，g⁺是混合矩陣g的偽取逆。

在一個(gè)實(shí)施例中，創(chuàng)建對(duì)角矩陣，在該對(duì)角矩陣中，其中，第一對(duì)角元素是s的逆對(duì)角元素：以及接下來(lái)的對(duì)角元素被設(shè)置為值1(如果其中，a是閾值)，或被設(shè)置為值0(如果)。

發(fā)現(xiàn)適合的閾值大約是0.06。在例如±0.01的范圍內(nèi)或在±10％的范圍內(nèi)的較小偏差是可接受的。然后，如下計(jì)算解碼矩陣：

在平滑解碼矩陣框44中，平滑解碼矩陣。替代現(xiàn)有技術(shù)中已知的在解碼之前向hoa系數(shù)應(yīng)用平滑系數(shù)，可以將其與解碼矩陣相組合。這節(jié)省了一個(gè)處理步驟或相應(yīng)節(jié)省了處理框。

為了針對(duì)用于hoa內(nèi)容比擴(kuò)音器具有更多系數(shù)(即，o3d＞l)的解碼器也獲得良好的節(jié)能特性，根據(jù)hoa階數(shù)n(o3d＝(n+1)²)來(lái)選擇所應(yīng)用的平滑系數(shù)

與在[4]中一樣，針對(duì)l≥o3d，對(duì)應(yīng)于根據(jù)階數(shù)n+1的勒讓德多項(xiàng)式的零導(dǎo)出的maxre系數(shù)。

針對(duì)l＜o(jì)3d，根據(jù)凱撒窗構(gòu)建的的系數(shù)如下所示：

其中，len＝2n+1，width＝2n，其中，是具有2n+1個(gè)取實(shí)數(shù)值的元素的矢量。該元素是通過(guò)凱撒窗公式創(chuàng)建的

其中，io()表示第一類(lèi)的零階修正貝塞爾函數(shù)。矢量是根據(jù)以下項(xiàng)構(gòu)建的：

其中，針對(duì)hoa階數(shù)索引n＝o..n，每個(gè)元素具有2n+1個(gè)重復(fù)，以及cf是用于在不同的hoa階數(shù)節(jié)目(program)之間保持相等的響度的恒定縮放因子。亦即，所使用的凱撒窗的元素開(kāi)始于僅被使用一次的第(n+1)個(gè)元素，并以被重復(fù)使用的后續(xù)元素繼續(xù)：第(n+2)個(gè)元素被使用3次，等等。

在一個(gè)實(shí)施例中，對(duì)已平滑的解碼矩陣進(jìn)行縮放。在一個(gè)實(shí)施例中，在圖4a)中示出的平滑解碼矩陣框44中執(zhí)行縮放。在不同的實(shí)施例中，在圖4b)中示出的縮放矩陣框45中將縮放作為單獨(dú)步驟執(zhí)行。

在一個(gè)實(shí)施例中，根據(jù)解碼矩陣來(lái)獲得恒定縮放因子。具體地，其可根據(jù)解碼矩陣的所謂弗羅比舍范數(shù)來(lái)獲得：

其中，是(已平滑的)矩陣的第l行第q列的矩陣元素。已歸一化的矩陣是

圖5示出了根據(jù)本發(fā)明的一方面的一種用于解碼音頻聲場(chǎng)表示以供音頻回放的設(shè)備。該設(shè)備包括具有用于獲得解碼矩陣d的解碼矩陣計(jì)算單元140的渲染處理單元33，解碼矩陣計(jì)算單元140包括用于獲得目標(biāo)揚(yáng)聲器的數(shù)目l的裝置1x和用于獲得揚(yáng)聲器的位置的裝置，用于確定球面建模網(wǎng)格的位置的裝置1y和用于獲得hoa階數(shù)n的裝置1z，以及用于根據(jù)球面建模網(wǎng)格的位置和揚(yáng)聲器的位置生成混合矩陣g的第一處理單元141，用于根據(jù)球面建模網(wǎng)格和hoa階數(shù)n生成模式矩陣的第二處理單元142，用于根據(jù)執(zhí)行對(duì)模式矩陣與厄米特轉(zhuǎn)置混合矩陣g的乘積的緊致奇異值分解的第三處理單元143(其中，u、v是根據(jù)酉矩陣導(dǎo)出的，以及s是具有奇異值元素的對(duì)角矩陣)，用于根據(jù)來(lái)由矩陣u、v計(jì)算第一解碼矩陣的計(jì)算裝置144，以及用于利用平滑系數(shù)對(duì)第一解碼矩陣進(jìn)行平滑和縮放的平滑及縮放單元145(其中，獲得解碼矩陣d)。在一個(gè)實(shí)施例中，平滑及縮放單元145例如是用于平滑第一解碼矩陣的平滑單元1451(其中獲得已平滑的解碼矩陣)和用于對(duì)已平滑的解碼矩陣進(jìn)行縮放的縮放單元1452(其中，獲得解碼矩陣d)。

圖6示出了節(jié)點(diǎn)示意圖中的示例性的16個(gè)揚(yáng)聲器的設(shè)置中的揚(yáng)聲器位置，其中，將揚(yáng)聲器示出為已連接節(jié)點(diǎn)。前臺(tái)連接示出為實(shí)線，后臺(tái)連接示出為虛線。圖7通過(guò)用透視法縮小繪制的視圖的形式示出了具有16個(gè)揚(yáng)聲器的相同設(shè)置。

下面描述利用圖5和圖6中的揚(yáng)聲器設(shè)置獲得的示例結(jié)果。在2個(gè)球面(所有的測(cè)試方向)上以db示出聲音信號(hào)的能量分布，以及具體地，比率的分布。示出中心揚(yáng)聲器波束(圖6中的揚(yáng)聲器7)作為擴(kuò)音器平移波束的示例。例如，在[14]中設(shè)計(jì)的解碼器矩陣(n＝3)產(chǎn)生圖8中示出的比率其提供了幾乎完美的節(jié)能特性，因?yàn)楸嚷?imgfile="bda0001244681640000192.gif"wi="154"he="79"img-content="drawing"img-format="gif"orientation="portrait"inline="no"/>幾乎是恒定的：暗區(qū)(對(duì)應(yīng)于較低音量)與亮區(qū)(對(duì)應(yīng)于較高音量)之間的差異小于0.01db。然而，如圖9中示出的，中心揚(yáng)聲器的對(duì)應(yīng)平移波束具有較強(qiáng)的副波瓣。特別是對(duì)于偏離中心的收聽(tīng)者而言，這妨礙了空間感知。

另一方面，在[2]中設(shè)計(jì)的解碼器矩陣(n＝3)產(chǎn)生圖9中示出的比率在圖10中使用的刻度中，暗區(qū)對(duì)應(yīng)于下降到-2db的較低音量，以及亮區(qū)對(duì)應(yīng)于上升到+2db的較高音量。因此，比率示出了大于4db的波動(dòng)，這是不利的，因?yàn)橐韵嗤捻懚炔荒芨兄胶愣ǚ鹊睦鐝捻敳康街行膿P(yáng)聲器位置的空間平移。然而，如圖11中示出的，中心揚(yáng)聲器的對(duì)應(yīng)平移波束具有非常小的副波瓣，這對(duì)于偏離中心的收聽(tīng)位置而言是有益的。

圖12示出了利用根據(jù)本發(fā)明的解碼器矩陣獲得的聲音信號(hào)的能量分布，為了易于比較，示例性地針對(duì)于n＝3。(在圖12的右側(cè)示出的)比率的刻度范圍從3.15到3.45db。因此，該比率中的波動(dòng)小于0.31db，以及聲場(chǎng)中的能量分布是非常均勻的。因此，以相同的響度感知到了具有恒定幅度的任何空間平移。如圖13中示出的，中心揚(yáng)聲器的平移波束具有非常小的副波瓣。這對(duì)于偏離中心的收聽(tīng)位置而言是有益的，在該收聽(tīng)位置處，副波瓣可能是可聽(tīng)到的，并因此將會(huì)是令人煩惱的。因此，本發(fā)明提供了利用[14]和[2]中的現(xiàn)有技術(shù)可獲得的組合優(yōu)點(diǎn)，而無(wú)需忍受其相應(yīng)缺點(diǎn)。

要注意到的是，在本文中只要提及揚(yáng)聲器，表示的是聲音發(fā)射設(shè)備，例如擴(kuò)音器。

附圖中的流程圖和/或框圖示出了根據(jù)本發(fā)明的各種實(shí)施例的系統(tǒng)、方法和計(jì)算機(jī)程序產(chǎn)品的可能實(shí)現(xiàn)的配置、操作和功能。在這一點(diǎn)上，流程圖或框圖中的每個(gè)框可以表示代碼的模塊、片段或部分，該代碼包括用于實(shí)現(xiàn)所指定的邏輯功能的一個(gè)或多個(gè)可執(zhí)行指令。

還應(yīng)該注意到，在一些備選實(shí)現(xiàn)中，框中提到的功能可以通過(guò)附圖中所提到的順序之外的順序來(lái)發(fā)生。例如，取決于所涉及的功能，連續(xù)示出的兩個(gè)框事實(shí)上可以實(shí)質(zhì)上同時(shí)執(zhí)行，或有時(shí)以相反的順序來(lái)執(zhí)行該框，或者可以通過(guò)備選順序執(zhí)行該框。還要注意到的是，框圖和/或流程圖示意圖的每個(gè)框，以及框圖和/或流程圖示意圖中的框組合可通過(guò)基于專用硬件的系統(tǒng)或者專用硬件和計(jì)算機(jī)指令的組合來(lái)實(shí)現(xiàn)，該基于專用硬件的系統(tǒng)執(zhí)行特定功能或動(dòng)作。雖然沒(méi)有明確描述，可以在任何組合或子組合中使用當(dāng)前實(shí)施例。

此外，本領(lǐng)域技術(shù)人員應(yīng)該意識(shí)到的是，本原理的各方面可以體現(xiàn)為系統(tǒng)、方法或計(jì)算機(jī)可讀介質(zhì)。因此，本原理的各方面可以采用以下形式：完全硬件實(shí)施例、完全軟件實(shí)施例(包括固件、駐機(jī)軟件、微代碼等)、或組合了軟件和硬件方面(在本文中可全部統(tǒng)稱為“電路”、“模塊”或“系統(tǒng)”)的實(shí)施例。此外，本原理的各方面可以采用計(jì)算機(jī)可讀存儲(chǔ)介質(zhì)的形式?？梢岳靡粋€(gè)或多個(gè)計(jì)算機(jī)可讀存儲(chǔ)介質(zhì)的任何組合。本文中使用的計(jì)算機(jī)可讀存儲(chǔ)介質(zhì)被認(rèn)為是非暫時(shí)存儲(chǔ)介質(zhì)，該非暫時(shí)存儲(chǔ)介質(zhì)被給定了在其中存儲(chǔ)信息的固有能力以及從其提供對(duì)信息的獲取的固有能力。

此外，本領(lǐng)域技術(shù)人員應(yīng)該意識(shí)到，這里所呈現(xiàn)的框圖表示了體現(xiàn)本發(fā)明原理的示意性系統(tǒng)組件和/或電路的概念視圖。類(lèi)似地，應(yīng)該意識(shí)到，任何流程圖、流程圖表、狀態(tài)轉(zhuǎn)移圖、偽代碼等表示各種過(guò)程，其中可以在計(jì)算機(jī)可讀存儲(chǔ)介質(zhì)中表示這些過(guò)程并且由計(jì)算機(jī)或處理器執(zhí)行這些過(guò)程，不管是否明確示出了這樣的計(jì)算機(jī)或處理器。

所引用的參考文獻(xiàn)

[1]t.d.abhayapala.generalizedframeworkforsphericalmicrophonearrays：spatialandfrequencydecomposition.inproc.ieeeinternationalconferenceonacoustics，speech，andsignalprocessing(icassp)，(accepted)vol.x，pp.，april2008，lasvegas，usa.

[2]johann-markusbatke，floriankeiler，andjohannesboehm.methodanddevicefordecodinganaudiosoundfieldrepresentationforaudioplayback.internationalpatentapplicationwo20111117399(pd100011).

[3]daniel，rozennnicol，andsébastienmoreau.furtherinvestigationsofhighorderambisonicsandwavefieldsynthesisforholophonicsoundimaging.inaesconventionpaper5788presentedatthe114thconvention，march2003.paper4795presentedatthe114thconvention.

[4]daniel.représentationdechampsacoustiques，applicationalatransmissionetalareproductiondescenessonorescomplexesdansuncontextemultimedia.phdthesis，universiteparis6，2001.

[5]jamesr.driscollanddennism.healyjr.computingfouriertransformsandconvolutionsonthe2-sphere.advancesinappliedmathematics，15：202-250，1994.

[6]fliege.integrationnodesforthesphere.

http：//www.personal.soton.ac.uk/jf1w07/nodes/nodes.html，online，accessed2012-06-01.

[7]fliegeandulrikemaier.atwo-stageapproachforcomputingcubatureformulaeforthesphere.technicalreport，fachbereichmathematik，dortmund，1999.

[8]r.h.hardinandn.j.a.sloane.webpage：sphericaldesigns，sphericalt-designs.http：//www2.research.att.com/～njas/sphdesigns/.

[9]r.h.hardinandn.j.a.sloane.mclaren′simprovedsnubcubeandothernewsphericaldesignsinthreedimensions.discreteandcomputationalgeometry，15：429-441，1996.

[10]m.a.poletfi.three-dimensionalsurroundsoundsystemsbasedonsphericalharmonics.j.audioeng.soc.，53(11)：1004-1025，november2005.

[11]villepulkki.spatialsoundgenerationandperceptionbyamplitudepanningtechniques.phdthesis，helsinkiuniversityoftechnology，2001.

[12]boazrafaely.plane-wavedecompositionofthesoundfieldonaspherebysphericalconvolution.j.acoust.soc.am.，4(116)：2149-2157，october2004.

[13]earlg.williams.fourieracoustics，volume93ofappliedmathematicalsciences.academicpress，1999.

[14]f.zotter.h.pomberger，andm.noisternig.energy-preservingambisonicdecoding.actaacusticaunitedwithacustica，98(1)：37-47.january/february2012.