專利名稱：：加窗信號mdct域的能量及相位調整方法及其裝置的制作方法
技術領域：
：本發(fā)明屬于空間音頻
技術領域：
，特別是涉及加窗信號在MDCT域的能量及相位調整方法及其裝置。
背景技術：
：音頻信號屬于短時平穩(wěn)過程，實際信號處理中，通過對信號采樣分幀，將具有平穩(wěn)統(tǒng)計特性的每幀信號作為處理對象。信號分幀相當于采用普通矩形窗對信號進行截短處理，由于矩形窗有陡峭的到零的下降沿，造成分幀后幀間的頻譜泄漏，使得重建后的聲音信號在幀間產生抖動。因此，有必要在變換前對時域信號加窗，避免因分幀產生的塊效應，保證幀間過渡的平滑。多聲道音頻編碼技術追求高音質和低碼率，為了重現(xiàn)原始聲音效果，在實際信號處理中通常以加窗信號作為研究對象。本發(fā)明所稱多聲道包括雙聲道及更多聲道的情況。對信號進行加窗處理時，通常應用的窗函數主要有矩形窗、三角窗、漢寧窗、海明窗和凱撒窗等類型。其中矩形窗主瓣較集中，旁瓣較高并有負旁瓣，易使變換產生高頻干擾和頻譜泄漏。與矩形窗相比，漢寧窗可緩解頻譜泄漏，但頻率分辨力下降。海明窗與漢寧窗均為余弦窗，僅加權系數不同。每種窗型都有自己的特點，實際信號處理時一般根據輸入信號的特性選擇窗型?？臻g音頻編碼技術作為多聲道音頻技術的典型代表，將加窗后的多聲道信號所包含的音頻信息和空間信息相分離并分別編碼?，F(xiàn)有的傳統(tǒng)空間音頻編解碼系統(tǒng)如圖l所示，編碼端輸入的多聲道信號加窗后，一方面通過下混成為一個8聲道進行傳統(tǒng)單聲道音頻編碼(MDCT時頻分析，包括MDCT變換和量化編碼)，另一方面經聲源分離(HQMF/CFB/STDFT時頻分析)后進行空間參數提取和空間參數編碼；解碼端通過碼流分析分離出下混聲道和空間參數信息，相應地分別對下混聲道進行傳統(tǒng)單聲道音頻解碼、上混，對空間參數信息進行空間參數解碼，最終通過合成濾波得到輸出多聲道信號。在空間音頻編碼系統(tǒng)中，空間參數反映多聲道音頻信號(包括雙聲道)的空間信息，編碼端空間參數的有效提取是解碼端音質還原的關鍵。空間參數主要包括用以定位聲源水平方位的聲道間強度差(InterauralLevelDifference,簡稱ILD)和聲道間時間差(InterauralTimeDifference,簡稱ITD)，以及用以確定聲源聲場寬度的聲道間相關度(InterauralCorrelation,簡稱IC)。其中，ILD記錄了聲道間信號的能量比，ITD記錄了聲道間信號具有最大互相關的時間偏移，IC記錄了此時的歸一化相關度。因此，ILD和ITD分別與時域加窗信號的能量及延時特性相關，IC同時受信號能量及延時特性的影響。實際信號處理中，通常在編碼端對加窗信號進行時頻域變換，完成信號的聲源分離。在MP3、AAC等通用音頻編碼系統(tǒng)中，傳統(tǒng)的單聲道音頻編碼采用修正離散余弦變換(ModifiedDiscreteCosineTransform,簡稱MDCT)作為最常見的時頻分析工具，且MDCT具有時域抗混疊特性；而編碼端在提取聲場空間信息時則采用模擬人耳聽覺特性的時頻分析工具，包括人工耳蝸濾波器(CochlearFilterBank,簡稱CFB)、短時離散傅立葉變換(Short-TimeDFT，簡稱STDFT)以及混合正交鏡像濾波器組(HybridQuadratureMirrorFilterbank，簡稱HQMF)。現(xiàn)有技術表明，若傳統(tǒng)的單聲道編碼與空間參數提取使用相同的時頻分析工具，可以有效降低編解碼的復雜度和算法延時。針對空間參數提取的時頻分析工具CFB、STDFT和HQMF由于復雜度、臨界采樣特性以及與心理聲學的匹配性等原因不適合在傳統(tǒng)的單聲道編碼器中使用。因此，為了降低復雜度和延時，當采用基于MDCT的感知音頻編碼器作為傳統(tǒng)單聲道編碼器時，研究如何提取加窗信號MDCT域上的空間參數成為空間音頻編碼領域的研究熱點。然而，MDCT變換屬于實變換，與同時包含信號相位信息和幅度信息的復變換不同，在MDCT域上提取空間參數存在困難首先MDCT譜的能量抖動使MDCT譜的聲道間強度差ILD與原加窗信號ILD不一致，導致ILD參數提取存在差異；其次MDCT譜中不包含直接的相位信息，難以反映原時域加窗信號的延時特性，導致直接提取聲道間時間差參數ITD存在困難。此外，空間音頻編解碼系統(tǒng)通常以加窗信號作為研究對象，由于加窗處理會對時域及變換域上的信號特性產生影響，在實際信號處理中通常采用與未加窗信號不同的處理方法，因此研究加窗信號在MDCT域上的能量及相位調整方法是有別于未加窗信號的，研究本身具有不可忽視的價值。
發(fā)明內容本發(fā)明目的在于針對現(xiàn)有技術的不足，提出了加窗信號在MDCT域的能量及相位調整方法及其裝置，以統(tǒng)一空間音頻編碼系統(tǒng)中的時頻分析工具。本發(fā)明提供的加窗信號MDCT域的能量及相位調整方法，技術方案包含以下步驟l，對輸入信號進行MDCT加窗得到加窗信號，以MDCT加窗變換所用的窗函數為共軛窗函數構造MDST變換的窗函數，具體構造方法如下，將MDCT變換的窗函數^(w)分解為兩個W點窗函數w。(")和Wl(")的銜<formula>formulaseeoriginaldocumentpage11</formula>取MDST變換的窗函數w(")為MDCT變換的窗函數的共軛窗函數，則v^(")滿足如下條件<formula>formulaseeoriginaldocumentpage11</formula>步驟2，根據MDCT變換的窗函數得到加窗信號MDCT變換的形式化表達，根據該形式化表達對加窗信號進行MDCT變換，得到加窗信號的MDCT譜；由MDST變換的窗函數得到加窗信號MDST變換的形式化表達；所述加窗信號MDCT變換的形式化表達如下<formula>formulaseeoriginaldocumentpage11</formula>所述加窗信號MDST變換的形式化表達如下=Z("K(")sin[^("+++)].*=0,1，…，〃一1其中，》("、r("分別為第z'幀加窗信號的MDCT譜和MDST譜，x'(")為輸入信號x(w)的第/幀信號，其中"表示第w個采樣點；〈,〉表示向量內積，q(")、^(w)分別為MDCT變換基和MDST變換基，w為MDCT變換的窗函數，w,(w)為MDST變換的窗函數；步驟3，由加窗信號的MDCT譜和加窗信號MDST變換的形式化表達構造加窗信號的MDST譜，包括以下四個步驟，步驟3.1,針對當前幀，從加窗信號的MDCT譜內提取并保留前一幀及后一幀譜線》"及》—、并做譜線運算，得到的運算后譜線分別記為il和《，具體譜線運算方式如下f丄=(—-X")/2，》丄=(—+f-1)/2;步驟3.2，根據加窗信號MDCT變換和加窗信號MDST變換的形式化表達，提取MDCT譜至MDST譜的轉換矩陣，包括和矩陣a及差矩陣b;提取具體方式如下，首先，對于2W維的MDCT變換基q和MDST變換基&，分別分解為W維列向:，即從而得到以下四個^x^矩陣形式的子向量C。、CpS。、S1:<formula>formulaseeoriginaldocumentpage12</formula>和矩陣a及差矩陣b根據以上子向量得到，a-s^;+s:q<formula>formulaseeoriginaldocumentpage12</formula>其中，A=0,1，...,7V-1,上標"T"表示轉置，上標"0"和"1"分別用于標示變換基前半部分和后半部分的子向量；步驟3.3，稀疏化近似步驟3.2所得和矩陣A，得到近似和矩陣a"所述稀疏化具體方式為，保留和矩陣a中絕對值最大的2m個元素值，其它元素均置為0值；其中m取小于WxA72的任意自然數；步驟3.4，根據步驟3.3所得近似和矩陣Am和步驟3.2所得差矩陣B對步驟3.1所得運算后譜線》1和》二進行轉換，實現(xiàn)MDCT譜到MDST譜的轉換，得到加窗信號的MDST譜；轉換公式如下TVm+TV—其中》二和fl分別為步驟3.1中所得運算后譜線，B為步驟3.2中所得差矩陣，Am為步驟3.3中所得近似和矩陣；步驟4，依據步驟2中所得加窗信號的MDCT譜和步驟3中所得加窗信號的MDST譜，構造加窗信號的CMCLT譜；所述CMCLT譜是共軛調制復重疊變換譜，是以變換時所用窗函數為共軛窗函數的MDCT譜和MDST譜分別作為實部和虛部，從而定義MDCT的一種復數化擴展形式，構造的CMCLT譜如下式所示其中，f'("、r("分別為第z'幀加窗信號的MDCT譜和MDST譜，x'為輸入信號的第/幀信號，q、^分別為MDCT變換基和MDST變換基，W。和W,分別為MDCT變換和MDST變換的窗函數矩陣表示，即其中W^&ag。(""，W,w。(")和w,(")分別為步驟1中窗函數^(w)分解后的兩個7V點窗函數；步驟5，根據加窗信號的CMCLT譜，提取加窗信號在CMCLT域的能量信息和相位信息；步驟6，根據步驟5所得能量信息及相位信息，在MDCT域上提取空間參數。而且，所述相位信息的提取過程包括以下三個步驟，2'("=》("-=〈x、WA〉-WA〉(l)將MDCT變換的窗函數^(")表示為離散正弦變換的四型變換基^(")的線性組合，確定展開系數，即其中，A為窗函數v^(")以離散正弦變換的四型變換基^(")展開的系數，A:取O，l,...,W-l，且四型變換基A(")表示如下■S/t(w)=sin工(w+丄)(A:+丄)iV22,w=0，l,..,,iV—l，(2)對第/幀加窗信號的CMCLT譜的實部》("和虛部r("做線性運算，構造f二("及乏U"如下(3)進行CMCLT譜的線性組合，將CMCLT譜轉換為DFT譜的形式，按照所述線性組合方式形成的新譜線f(A:)如下式所示2("=—Z(-1)、-2/-1)++2/))/=0/V/2-1—S(乏W+2/+1)+-2/一2))上式表現(xiàn)的關系展示了一個CMCLT譜的線性組合，等價于輸入信號的DFT譜線附加一個與輸入信號無關的相位偏移，且線性組合的形式依賴于四型變換基A(")展開的系數。本發(fā)明還提供了實現(xiàn)上述MDCT域信號能量與相位補償方法的相應裝置，包含以下部分MDCT加窗模塊，用于對輸入信號進行MDCT加窗得到加窗信號；共軛窗構造模塊，用于以MDCT加窗變換所用的窗函數為共軛窗函數構造MDST變換的窗函數；加窗MDST表示模塊，用于MDST變換的窗函數得到加窗信號MDST變換的形式化表達；加窗MDCT變換模塊，用于根據MDCT變換的窗函數得到加窗信號MDCT變換的形式化表達，根據該形式化表達對加窗信號進行MDCT變換，得到加窗信號的MDCT譜；加窗MDST構造模塊，用于由加窗信號的MDCT譜和加窗信號MDST變換的形式化表達構造加窗信號的MDST譜；加窗CMCLT構造模塊，用于依據加窗信號的MDCT譜和加窗信號的MDST譜，構造加窗信號的CMCLT譜；能量提取模塊，用于根據加窗信號的CMCLT譜，提取加窗信號在CMCLT域的能量信息；相位提取模塊，用于根據加窗信號的CMCLT譜，提取加窗信號在CMCLT域的相位信息；空間參數提取模塊，用于根據加窗信號在CMCLT域的能量信息和相位信息，在在MDCT域上提取空間參數；其中，輸入信號進入MDCT加窗模塊，MDCT加窗模塊分別連接加窗MDCT變換模塊和共軛窗構造模塊，MDCT加窗模塊通過連接輸出加窗信號進入加窗MDCT變換模塊，輸出MDCT加窗變換所用的窗函數到共軛窗構造模塊；共軛窗構造模塊與加窗MDST表示模塊相連，共軛窗構造模塊通過連接輸出MDST變換的窗函數到加窗MDST表示模塊；加窗MDST表示模塊和加窗MDCT變換模塊的輸出分別接入加窗MDST構造模塊，加窗MDST表示模塊通過連接輸出加窗信號MDST變換的形式化表達到加窗MDST構造模塊，加窗MDCT變換模塊輸出加窗信號的MDCT譜到加窗MDST構造模塊；加窗MDST構造模塊和加窗MDCT變換模塊的輸出分別接入加窗CMCLT構造模塊，加窗MDST構造模塊通過連接輸出加窗信號的MDST譜到加窗CMCLT構造模塊，加窗MDCT變換模塊通過連接輸出加窗信號的MDCT譜到加窗CMCLT構造模塊；加窗CMCLT構造模塊的輸出分別接入能量提取模塊和相位提取模塊，能量提取模塊所得能量信息和相位提取模塊所得相位信息最終進入空間參數提取模塊，空間參數提取模塊輸出在MDCT域上提取的空間參數。本發(fā)明提供了加窗信號在MDCT域的能量及相位調整方法及其裝置，可統(tǒng)一空間音頻編碼系統(tǒng)中的時頻分析工具，有效提取多聲道信號的空間參數。統(tǒng)一時頻變換工具后，整個編碼端只需進行一次時頻變換，較大幅度降低了空間音頻編解碼系統(tǒng)的復雜度，從而優(yōu)化編解碼工作。圖1為現(xiàn)有空間音頻編解碼系統(tǒng)結構框圖。圖2為本發(fā)明的裝置結構框圖。圖3為本發(fā)明實施例的加窗MDST構造模塊實現(xiàn)框圖。圖4為本發(fā)明實施例的相位提取模塊的實現(xiàn)框圖。圖5為本發(fā)明在空間音頻編解碼系統(tǒng)中的應用示意圖。具體實施例方式本發(fā)明提供的技術方案流程如下，具體實施時可以采用計算機實現(xiàn)以下流程的自動運行步驟l，對輸入信號進行MDCT加窗得到加窗信號，以MDCT加窗變換所用的窗函數為共軛窗函數構造MDST變換的窗函數。在空間音頻
技術領域：
，MDST為常用變換手段，意為修正離散正弦變換(ModifiedDiscreteSineTransform,簡稱MDST)，即MDCT的共軛變換。以MDCT加窗變換所用的窗函數為共軛窗函數構造MDST變換的窗函數，具體構造方法如下將MDCT變換的窗函數we(")分解為兩個W點窗函數w。(")和Wl(w)的銜接取MDST變換的窗函數w,(")為MDCT變換的窗函數的共軛窗函數，則A(")滿足如下條件本發(fā)明實施例對輸入信號x(")加窗采用正弦窗型，具體實施時若采用其它窗型，都可參考該實施例實現(xiàn)。因為采用正弦窗型，信號MDCT變換的窗函數可表示為<formula>formulaseeoriginaldocumentpage17</formula>(式D分析知的共軛窗函數為余弦窗，信號MDST變換的窗函數可表示為-<formula>formulaseeoriginaldocumentpage17</formula>(式2)步驟2，根據MDCT變換的窗函數得到加窗信號MDCT變換的形式化表達，根據該形式化表達對加窗信號進行MDCT變換，得到加窗信號的MDCT譜；由MDST變換的窗函數得到加窗信號MDST變換的形式化表達。所述加窗信號MDCT變換的形式化表達如下<formula>formulaseeoriginaldocumentpage17</formula>地"X'(")，然("》2JV-17T1AA1所述加窗信號MDST變換的形式化表達如下2AMT1A/"1=Sx'(")—)sin[:("+^+,+)].&=0,1，…，〃-1其中，f'("、r("分別為第/幀加窗信號的MDCT譜和MDST譜，jc'(")為輸入信號x(")的第/幀信號，其中"表示第"個采樣點；〈，〉表示向量內積，q(")、^(")分別為MDCT變換基和MDST變換基，w為MDCT變換的窗函數，w為MDST變換的窗函數；本發(fā)明實施例根據式1和式2，取得(A)加窗信號MDCT變換的形式化表達如下=，x'(w)c。s[丄(w+丄)]sin[-("+丄+—X*+丄)].A:=O，l,...，iV-1(式3)(B)加窗信號MDST變換的形式化表達如下=《x'(")Sin[^("+會)]cos[^:("+會++1)].*=。,1"."W—1(式4)步驟3，由加窗信號的MDCT譜和加窗信號MDST變換的形式化表達構造加窗信號的MDST譜，包括以下四個步驟步驟3.1，針對當前幀，從加窗信號的MDCT譜內提取并保留前一幀及后一幀譜線f'"及X"，并做譜線運算，得到的運算后譜線分別記為^和《，具體譜線運算方式如下》(=(-f—，1)/2,《=(-X'+'+，1)/2。步驟3.2，根據加窗信號MDCT變換和加窗信號MDST變換的形式化表達，提取MDCT譜至MDST譜的轉換矩陣，包括和矩陣A及差矩陣B;提取具體方式如下，首先，對于2W維的MDCT變換基q和MDST變換基&，分別分解為7V維列向量，即")T,)T(。"")T=("0)TMDCT變換基q和MDST變換基^各有N個，因此得到以下四個WxW矩陣形式的子向量C。、C,、S。、S1:Ccs<s,和矩陣A及差矩陣B根據以上子向量得到，A-S:C。+S:q，B=S^Co-SjC，。其中，&=0，1，...,7V-1，上標"T"表示轉置，上標"0"和"1"分別用于標示變換基前半部分和后半部分的子向量。步驟3.3，稀疏化近似步驟3.2所得和矩陣A，得到近似和矩陣Am;所述稀疏化具體方式為，保留和矩陣A中絕對值最大的2m個元素值，其它元素均置為0值；其中m取小于(WxA0/2的任意自然數。在步驟3.3中，差矩陣B是一個符號單位陣的轉置，和矩陣A是一個稀疏對角陣。因此可以針對和矩陣，將和矩陣進行稀疏化近似，即保留和矩陣A中絕對值最大的2m個元素值，其它元素均置為O值，并將稀疏化近似后的和矩陣A記為近似和矩陣Am。步驟3.4，根據步驟3.3所得近似和矩陣Am和步驟3.2所得差矩陣B對步驟3.1所得運算后譜線f二和i:進行轉換，實現(xiàn)MDCT譜到MDST譜的轉換，得到加窗信號的MDST譜；轉換公式如下TVm+AT-其中i(和f二分別為步驟3.1中所得運算后譜線，B為步驟3.2中所得差矩陣，Am為步驟3.3中所得近似和矩陣。步驟4，依據步驟2中所得加窗信號的MDCT譜和步驟3中所得加窗信號的MDST譜，構造加窗信號的CMCLT譜；所述CMCLT譜是共軛調制復重疊變換譜，是以變換時所用窗函數為共軛窗函數的MDCT譜和MDST譜分別作為實部和虛部，從而定義MDCT的一種復數化擴展形式，構造的CMCLT譜如下式所示=》詢-("=〈？,WA〉—W,WA〉其中，》("、r(A:)分別為第/幀加窗信號的MDCT譜和MDST譜，^為輸入信號的第/幀信號，q、^分別為MDCT變換基和MDST變換基。值得注意的是，步驟2中x'(")、q(")、^(w)加(")是計算需要，"指第/幀信號的第"個采樣點，此次不加《表示直接針對第z'幀信號整體考慮，運算包括該幀信號內所有采樣點。Wj卩W,分別為MDCT變換和MDST變換的窗函數矩陣表示，即'w00、「w，0、,ws=s卜-wj其中\(zhòng)￥。=^^{>^(")}，W「^^(w("》，w。(")和w,(")分別為步驟1中窗函數m^(")分解后的兩個W點窗函數。Aag(w。("))表示對角矩陣W。的對角線上元素為w。(")，"/ag(w"""表示對角矩陣W,的對角線上元素為。本發(fā)明實施例中，就是用步驟2所得加正弦窗的MDCT譜》("和步驟3中構造的MDST譜r(A:)分別作為實部和虛部，構造加窗信號的CMCLT譜2'("。步驟5，根據加窗信號的CMCLT譜，提取加窗信號在CMCLT域的能量信息和相位f曰息。本發(fā)明對于加窗信號在CMCLT域的能量信息進行了一番考察將加窗信號在CMCLT域的能量記為Em,未加窗的時域信號能量記為En，令CMCLT變換矩陣R-(W。CWSS)，則RRT=(WCCWSS,CTW,、STW。=WCCTW+\V。SSTW'20八0+—2=—M20—w,o乂o)rw。oin+jnJ〔ow'oYw'ol乂0In—乂0—w,o乂2n其中Wj卩W,分別為MDCT變換和MDST變換的窗函數矩陣表示，i2n為27Vx2W的單位陣、^為7VxiV的單位陣、^是A^xiV的反對角單位陣。C和S分別為MDCT變換基q和MDST變換基&組成的矩陣C=(c。q…Cw一)S=(s0^...《n因此，有下式成立:A=0w—l4=0=(("(n7廣》f信號其中，》(A)為第z'幀加窗信號的CMCLT譜，》'(yt)和f("分別為第/幀力的MDCT譜和MDST譜。上式中加窗信號在CMCLT域的能量Em與未加窗的時域信號能量En僅相差常數2倍，該關系式即為加窗信號在CMCLT域的能量信息。2本發(fā)明對于加窗信號在CMCLT域的相位信息提供了具體提取方案(1)MDCT變換的窗函數^(")表示為離散正弦變換的四型變換基L(w)的線性組合，確定展開系數，即wc(w)="0>yo+"一十…+"一+...+一1"^-1其中，A為窗函數w(n)以基函數L(")展開的展開系數，且基函數L(")表示如下&=sin工("+丄)("丄)iV、22,w-0，l，…，iV—l，本發(fā)明實施例中，將輸入信號加窗并進行MDCT變換所采用的正弦窗w表示為離散正弦變換的四型變換基L(w)的線性組合，即;r1一wc(w)=sin[——(《+—)]="0>o27V2上式中展開系數除《。=1，其余均為o。離散正弦變化的四型變換基簡稱DST-IV變換基，MDCT變換基和MDST變換基是其變形，這三者均屬于現(xiàn)有技術。(2)對第/幀加窗信號CMCLT譜的實部》("和虛部f'(fc)做線性運算，構造之U"及乏UW如下Z'+(A:)=r(A:)+》(A:)ZU"-r("-》("(式5)(3)構造f'(&)，實現(xiàn)CMCLT譜的線性組合，將CMCLT譜轉換為DFT譜的形式，其中F(A:)如下式所示力肪-1Z("--Z(-1)、(《(A:-2/-1)++2/))W/2-l-￡(-1)、(A*+2/+1)+,(A:-2/-2))其中，e為數學常量2.71828183，/取0，1,...，^/2-1，A;取0，1,…,〃-1。上式中表達一個CMCLT譜線的線性組合，等價于輸入信號的DFT譜線附加一個與輸入信號無關的相位偏移，且線性組合的形式依賴于四型變換基L(")展開的系數。DFT是離散傅里葉變換的簡稱，由于DFT變換具有線性相位，則構造的，'(/fc)亦具有和時域延時相對應的線性相位。根據上式，將周期為2iV的第/幀加窗信.龍，_i號x'(")延時d，即《(w)二x'("-力，則延時后的CMCLT譜Zrf(A:)-Z(A:)e，因此信號的時域延時對應于CMCLT譜的線性附加相位-!汲。為了便于理解本發(fā)明技術方案，以下提供本發(fā)明實施例，'(A)的構造公式由來首先依據一般窗型CMCLT譜的線性組合方法，將CMCLT譜轉換為DFT譜的形式，其中一般窗型乏("的構造如下式所示<formula>formulaseeoriginaldocumentpage24</formula>將A("、乏UW及展開系數A代入公式6，構造2(",其中實施例對^滿足下式<formula>formulaseeoriginaldocumentpage24</formula>因此式6中僅有/=0的項有效，即<formula>formulaseeoriginaldocumentpage24</formula>將式3、式4和式5代入式8化簡得<formula>formulaseeoriginaldocumentpage24</formula>由式9知，通過對加窗信號CMCLT譜實部和虛部的線性組合，得到一個與DFT僅有相位差的變換。依據式9，將周期為2iV的第z'幀加窗信號^(w)延時d，即《(")=x'("-J)，則延時后的CMCLT譜^(fc):2'("e力、因此信號的時域延7T時對應于CMCLT譜的線性附加相位-:說。步驟6，根據步驟5所得能量信息及相位信息，在MDCT域上提取空間參數。在步驟5知道CMCLT域的能量信息和相位信息后，就可以在MDCT域上提取等同的空間參數。具體實施時，往往需要多個空間參數反映多聲道音頻信號的空間信息。在步驟6中，依據步驟5中所得加窗信號CMCLT譜的能量信息提取MDCT域上與能量相關的空間參數，依據步驟5中所得加窗信號CMCLT譜的相位信息提取MDCT域上與相位相關的空間參數。為了具體說明空間參數的提取，本步驟以雙耳線索ILD，ITD為實施例具體說明，具體實施方式如下(A)在MDCT域提取ILD信息，ILD記錄了左右聲道信號的能量比ILD在時域的值ILDt在時域上的定義公式為ILDt=101ogl。(j、2(,)A/p(0&)，但由于信號處理不是在時域上進行，而是在MDCT域上進行，因此根據CMCLT域的能量信息按下式進行調整-ILDm=201og10p-201og10其中，ILDm為在MDCT域提取的ILD的值，x,(/)和;c力)分別為左右聲道信號，《和《分別為左右聲道第f'幀加窗信號的CMCLT譜。因為存在能量的倍數關系，所以可以用左右聲道信號在MDCT域上的CMCLT譜能量之比，從而得到MDCT域上的值IIX^，從而得到等價于時域上的值ILDt的左右聲道信號的能量比。(B)在MDCT域提取ITD信息，ITD記錄了左右聲道信號最大互相關的時間偏移ITD在時域的值ITDt在時域上的定義公式為ITDt-argmaxO/Wjaz+r)*),但由于信號處理不是在時域上進行，而是在MDCT域上進行，因此根據CMCLT域的能量信息按下式進行調整ITDm=:rarg(i;(A:)《*(")欲其中ITDt為ITD在時域的值，ITDm為在MDCT域提取的ITD的值，x々)和x々)分別為左右耳入聲，A=0，1，...，W—1。r為具有最大互相關時的時間偏移。通過計算左右聲道第/幀加窗信號的CMCLT譜《("和乏二(A)相位差的群延時提取MDCT域的ITD值。具體實施時，除了采用軟件流程方式實現(xiàn)本發(fā)明提供的加窗信號MDCT域的能量及相位調整技術方案外；還可以通過模塊化技術提供加窗信號MDCT域的能量及相位調整裝置，以支持面向用戶的能量及相位調整工作。參見圖2，加窗信號MDCT域的能量及相位調整裝置可包含以下部分MDCT加窗模塊，用于對輸入信號進行MDCT加窗得到加窗信號；共軛窗構造模塊，用于以MDCT加窗變換所用的窗函數為共軛窗函數構造MDST變換的窗函數；加窗MDST表示模塊，用于MDST變換的窗函數得到加窗信號MDST變換的形式化表達；加窗MDCT變換模塊，用于根據MDCT變換的窗函數得到加窗信號MDCT變換的形式化表達，根據該形式化表達對加窗信號進行MDCT變換，得到加窗信號的MDCT譜；加窗MDST構造模塊，用于由加窗信號的MDCT譜和加窗信號MDST變換的形式化表達構造加窗信號的MDST譜；加窗CMCLT構造模塊，用于依據加窗信號的MDCT譜和加窗信號的MDST譜，構造加窗信號的CMCLT譜；能量提取模塊，用于根據加窗信號的CMCLT譜，提取加窗信號在CMCLT域的能量〈曰息；相位提取模塊，用于根據加窗信號的CMCLT譜，提取加窗信號在CMCLT域的相位信息；空間參數提取模塊，用于根據加窗信號在CMCLT域的能量信息和相位信息，在在MDCT域上提取空間參數；其中，輸入信號進入MDCT加窗模塊，MDCT加窗模塊分別連接加窗MDCT變換模塊和共軛窗構造模塊，MDCT加窗模塊通過連接輸出加窗信號進入加窗MDCT變換模塊，輸出MDCT加窗變換所用的窗函數到共軛窗構造模塊；共軛窗構造模塊與加窗MDST表示模塊相連，共軛窗構造模塊通過連接輸出MDST變換的窗函數到加窗MDST表示模塊；加窗MDST表示模塊和加窗MDCT變換模塊的輸出分別接入加窗MDST構造模塊，加窗MDST表示模塊通過連接輸出加窗信號MDST變換的形式化表達到加窗MDST構造模塊，加窗MDCT變換模塊輸出加窗信號的MDCT譜到加窗MDST構造模塊；加窗MDST構造模塊和加窗MDCT變換模塊的輸出分別接入加窗CMCLT構造模塊，加窗MDST構造模塊通過連接輸出加窗信號的MDST譜到加窗CMCLT構造模塊，加窗MDCT變換模塊通過連接輸出加窗信號的MDCT譜到加窗CMCLT構造模塊；加窗CMCLT構造模塊輸出的CMCLT譜分別接入能量提取模塊和相位提取模塊，能量提取模塊所得能量信息和相位提取模塊所得相位信息輸入空間參數提取模塊，空間參數提取模塊輸出在MDCT域上提取的空間參數。具體實施時，還可以進行更細致的劃分設計，例如加窗MDST構造模塊由譜線運算單元、加窗轉換矩陣提取單元、加窗轉換矩陣簡化單元和MDCT-MDST轉換單元構成；相位提取模塊由窗函數展開單元、CMCLT譜線運算單元和CMCLT-DFT轉換單元構成?？蓞⒁妶D3，加窗MDST構造模塊的譜線運算單元用于實現(xiàn)步驟3.1，輸出運算后譜線^和f二到加窗轉換矩陣提取單元；加窗轉換矩陣提取單元用于實現(xiàn)步驟3.2，輸出和矩陣A到加窗轉換矩陣簡化單元，輸出及差矩陣B到MDCT-MDST轉換單元；加窗轉換矩陣簡化單元用于實現(xiàn)步驟3.3，輸出近似和矩陣Am到MDCT-MDST轉換單元；MDCT-MDST轉換單元用于實現(xiàn)步驟3.4，輸出所得加窗信號的MDST譜?？蓞⒁妶D4，相位提取模塊的窗函數展開單元用于將信號MDCT變換所加窗型A(")表示為DST-IV變換基^(")的線性組合，確定展開系數^并輸出到CMCLT譜線運算單元；CMCLT譜線運算單元用于對第/幀加窗信號CMCLT譜的實部》(A:)和虛部r做線性運算，構造&④及^("并輸出到CMCLT-DFT轉換單元；CMCLT-DFT轉換單元用于構造，'("，實現(xiàn)CMCLT譜的線性組合，將CMCLT譜轉換為DFT譜的形式。將本發(fā)明用于空間音頻編解碼系統(tǒng)，基于本發(fā)明的空間音頻編解碼系統(tǒng)應用裝置如圖5所示。圖5將應用本發(fā)明前后的空間音頻編解碼系統(tǒng)作了對比。圖片左方是未應用本發(fā)明的現(xiàn)有空間音頻編解碼系統(tǒng)，其中的空間參數提取和傳統(tǒng)單聲道音頻編碼運用獨立的時頻分析工具，分別為HQMF/CFB/STDFT以及MDCT變換。圖片右方是基于本發(fā)明的空間音頻編解碼系統(tǒng)，其中的空間參數提取和傳統(tǒng)單聲道音頻編碼共同采用MDCT時頻變換。空間參數提取可通過本發(fā)明提供的裝置實現(xiàn)，在空間參數提取模塊基礎上擴展共軛窗構造模塊、加窗MDST表示模塊、加窗MDCT變換模塊、加窗MDST構造模塊、加窗CMCLT構造模塊、能量提取模塊、相位提取模塊。統(tǒng)一時頻變換工具后，整個編碼端只需進行一次時頻變換，較大幅度降低了空間音頻編碼系統(tǒng)的復雜度。相應的解碼端合成濾波采用MDCT時頻變換的逆變換(即IMDCT變換)即可。權利要求1.一種加窗信號MDCT域的能量及相位調整方法，其特征在于，包含以下步驟步驟1，對輸入信號進行MDCT加窗得到加窗信號，以MDCT加窗變換所用的窗函數為共軛窗函數構造MDST變換的窗函數，具體構造方法如下，將MDCT變換的窗函數wc(n)分解為兩個N點窗函數w0(n)和w1(n)的銜接<mathsid="math0001"num="0001"><math><![CDATA[<mrow><msub><mi>w</mi><mi>c</mi></msub><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><mo>=</mo><mfencedopen='{'close=''><mtable><mtr><mtd><msub><mi>w</mi><mn>0</mn></msub><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><mo>,</mo></mtd><mtd><mi>n</mi><mo>=</mo><mn>0,1</mn><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><mi>N</mi><mo>-</mo><mn>1</mn></mtd></mtr><mtr><mtd><msub><mi>w</mi><mn>1</mn></msub><mrow><mo>(</mo><mi>n</mi><mo>-</mo><mi>N</mi><mo>)</mo></mrow><mo>,</mo></mtd><mtd><mi>n</mi><mo>=</mo><mi>N</mi><mo>,</mo><mi>N</mi><mo>+</mo><mn>1</mn><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><mn>2</mn><mi>N</mi><mo>-</mo><mn>1</mn></mtd></mtr></mtable></mfenced><mo>,</mo></mrow>]]></math></maths>取MDST變換的窗函數ws(n)為MDCT變換的窗函數wc(n)的共軛窗函數，則ws(n)滿足如下條件<mathsid="math0002"num="0002"><math><![CDATA[<mrow><msub><mi>w</mi><mi>s</mi></msub><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><mo>=</mo><mfencedopen='{'close=''><mtable><mtr><mtd><msub><mi>w</mi><mn>1</mn></msub><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><mo>,</mo></mtd><mtd><mi>n</mi><mo>=</mo><mn>0,1</mn><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><mi>N</mi><mo>-</mo><mn>1</mn></mtd></mtr><mtr><mtd><msub><mrow><mo>-</mo><mi>w</mi></mrow><mn>0</mn></msub><mrow><mo>(</mo><mi>n</mi><mo>-</mo><mi>N</mi><mo>)</mo></mrow><mo>,</mo></mtd><mtd><mi>n</mi><mo>=</mo><mi>N</mi><mo>,</mo><mi>N</mi><mo>+</mo><mn>1</mn><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><mn>2</mn><mi>N</mi><mo>-</mo><mn>1</mn></mtd></mtr></mtable></mfenced><mo>;</mo></mrow>]]></math></maths>步驟2，根據MDCT變換的窗函數得到加窗信號MDCT變換的形式化表達，根據該形式化表達對加窗信號進行MDCT變換，得到加窗信號的MDCT譜；由MDST變換的窗函數得到加窗信號MDST變換的形式化表達；所述加窗信號MDCT變換的形式化表達如下<mathsid="math0003"num="0003"><math><![CDATA[<mrow><mo>=</mo><munderover><mi>Σ</mi><mrow><mi>n</mi><mo>=</mo><mn>0</mn></mrow><mrow><mn>2</mn><mi>N</mi><mo>-</mo><mn>1</mn></mrow></munderover><msup><mi>x</mi><mi>i</mi></msup><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><msub><mi>w</mi><mi>c</mi></msub><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><mi>cos</mi><mo>[</mo><mfrac><mi>π</mi><mi>N</mi></mfrac><mrow><mo>(</mo><mi>n</mi><mo>+</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><mo>+</mo><mfrac><mi>N</mi><mn>2</mn></mfrac><mo>)</mo></mrow><mrow><mo>(</mo><mi>k</mi><mo>+</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><mo>)</mo></mrow><mo>]</mo><mo>.</mo><mi>k</mi><mo>=</mo><mn>0,1</mn><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><mi>N</mi><mo>-</mo><mn>1</mn></mrow>]]></math></maths>所述加窗信號MDST變換的形式化表達如下<mathsid="math0004"num="0004"><math><![CDATA[<mrow><mo>=</mo><munderover><mi>Σ</mi><mrow><mi>n</mi><mo>=</mo><mn>0</mn></mrow><mrow><mn>2</mn><mi>N</mi><mo>-</mo><mn>1</mn></mrow></munderover><msup><mi>x</mi><mi>i</mi></msup><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><msub><mi>w</mi><mi>s</mi></msub><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><mi>sin</mi><mo>[</mo><mfrac><mi>π</mi><mi>N</mi></mfrac><mrow><mo>(</mo><mi>n</mi><mo>+</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><mo>+</mo><mfrac><mi>N</mi><mn>2</mn></mfrac><mo>)</mo></mrow><mrow><mo>(</mo><mi>k</mi><mo>+</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><mo>)</mo></mrow><mo>]</mo><mo>.</mo><mi>k</mi><mo>=</mo><mn>0,1</mn><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><mi>N</mi><mo>-</mo><mn>1</mn></mrow>]]></math></maths>其中，id="icf0007"file="A2009100620280002C7.tif"wi="27"he="5"top="241"left="43"img-content="drawing"img-format="tif"orientation="portrait"inline="yes"/>分別為第i幀加窗信號的MDCT譜和MDST譜，xi(n)為輸入信號x(n)的第i幀信號，其中n表示第n個采樣點；<，>表示向量內積，ck(n)、sk(n)分別為MDCT變換基和MDST變換基，wc(n)為MDCT變換的窗函數，ws(n)為MDST變換的窗函數；步驟3，由加窗信號的MDCT譜和加窗信號MDST變換的形式化表達構造加窗信號的MDST譜，包括以下四個步驟，步驟3.1，針對當前幀，從加窗信號的MDCT譜內提取并保留前一幀及后一幀譜線id="icf0008"file="A2009100620280003C1.tif"wi="7"he="4"top="91"left="38"img-content="drawing"img-format="tif"orientation="portrait"inline="yes"/>及id="icf0009"file="A2009100620280003C2.tif"wi="7"he="4"top="91"left="52"img-content="drawing"img-format="tif"orientation="portrait"inline="yes"/>，并做譜線運算，得到的運算后譜線分別記為id="icf0010"file="A2009100620280003C3.tif"wi="5"he="5"top="91"left="155"img-content="drawing"img-format="tif"orientation="portrait"inline="yes"/>和id="icf0011"file="A2009100620280003C4.tif"wi="5"he="5"top="91"left="166"img-content="drawing"img-format="tif"orientation="portrait"inline="yes"/>，具體譜線運算方式如下步驟3.2，根據加窗信號MDCT變換和加窗信號MDST變換的形式化表達，提取MDCT譜至MDST譜的轉換矩陣，包括和矩陣A及差矩陣B；提取具體方式如下，首先，對于2N維的MDCT變換基ck和MDST變換基sk，分別分解為N維列向量，即<mathsid="math0005"num="0005"><math><![CDATA[<mrow><msup><mrow><mo>(</mo><msub><mi>c</mi><mi>k</mi></msub><mo>)</mo></mrow><mi>T</mi></msup><mo>=</mo><mfencedopen='('close=')'><mtable><mtr><mtd><msup><mrow><mo>(</mo><msubsup><mi>c</mi><mi>k</mi><mn>0</mn></msubsup><mo>)</mo></mrow><mi>T</mi></msup></mtd><mtd><msup><mrow><mo>(</mo><msubsup><mi>c</mi><mi>k</mi><mn>1</mn></msubsup><mo>)</mo></mrow><mi>T</mi></msup></mtd></mtr></mtable></mfenced></mrow>]]></math></maths><mathsid="math0006"num="0006"><math><![CDATA[<mrow><msup><mrow><mo>(</mo><msub><mi>s</mi><mi>k</mi></msub><mo>)</mo></mrow><mi>T</mi></msup><mo>=</mo><mfencedopen='('close=')'><mtable><mtr><mtd><msup><mrow><mo>(</mo><msubsup><mi>s</mi><mi>k</mi><mn>0</mn></msubsup><mo>)</mo></mrow><mi>T</mi></msup></mtd><mtd><msup><mrow><mo>(</mo><msubsup><mi>s</mi><mi>k</mi><mn>1</mn></msubsup><mo>)</mo></mrow><mi>T</mi></msup></mtd></mtr></mtable></mfenced></mrow>]]></math></maths>從而得到以下四個N×N矩陣形式的子向量C0、C1、S0、S1<mathsid="math0007"num="0007"><math><![CDATA[<mrow><msub><mi>C</mi><mn>0</mn></msub><mo>=</mo><mfencedopen='('close=')'><mtable><mtr><mtd><msubsup><mi>c</mi><mn>0</mn><mn>0</mn></msubsup></mtd><mtd><msubsup><mi>c</mi><mn>1</mn><mn>0</mn></msubsup></mtd><mtd><mo>.</mo><mo>.</mo><mo>.</mo></mtd><mtd><msubsup><mi>c</mi><mrow><mi>N</mi><mo>-</mo><mn>1</mn></mrow><mn>0</mn></msubsup></mtd></mtr></mtable></mfenced></mrow>]]></math></maths><mathsid="math0008"num="0008"><math><![CDATA[<mrow><msub><mi>C</mi><mn>1</mn></msub><mo>=</mo><mfencedopen='('close=')'><mtable><mtr><mtd><msubsup><mi>c</mi><mn>0</mn><mn>1</mn></msubsup></mtd><mtd><msubsup><mi>c</mi><mn>1</mn><mn>1</mn></msubsup></mtd><mtd><mo>.</mo><mo>.</mo><mo>.</mo></mtd><mtd><msubsup><mi>c</mi><mrow><mi>N</mi><mo>-</mo><mn>1</mn></mrow><mn>1</mn></msubsup></mtd></mtr></mtable></mfenced></mrow>]]></math></maths><mathsid="math0009"num="0009"><math><![CDATA[<mrow><msub><mi>S</mi><mn>0</mn></msub><mo>=</mo><mfencedopen='('close=')'><mtable><mtr><mtd><msubsup><mi>s</mi><mn>0</mn><mn>0</mn></msubsup></mtd><mtd><msubsup><mi>s</mi><mn>1</mn><mn>0</mn></msubsup></mtd><mtd><mo>.</mo><mo>.</mo><mo>.</mo></mtd><mtd><msubsup><mi>s</mi><mrow><mi>N</mi><mo>-</mo><mn>1</mn></mrow><mn>0</mn></msubsup></mtd></mtr></mtable></mfenced></mrow>]]></math></maths><mathsid="math0010"num="0010"><math><![CDATA[<mrow><msub><mi>s</mi><mn>1</mn></msub><mo>=</mo><mfencedopen='('close=')'><mtable><mtr><mtd><msubsup><mi>s</mi><mn>0</mn><mn>1</mn></msubsup></mtd><mtd><msubsup><mi>s</mi><mn>1</mn><mn>1</mn></msubsup></mtd><mtd><mo>.</mo><mo>.</mo><mo>.</mo></mtd><mtd><msubsup><mi>s</mi><mrow><mi>N</mi><mo>-</mo><mn>1</mn></mrow><mn>1</mn></msubsup></mtd></mtr></mtable></mfenced></mrow>]]></math></maths>和矩陣A及差矩陣B根據以上子向量得到，<mathsid="math0011"num="0011"><math><![CDATA[<mrow><mi>A</mi><mo>=</mo><msubsup><mi>S</mi><mn>1</mn><mi>T</mi></msubsup><msub><mi>C</mi><mn>0</mn></msub><mo>+</mo><msubsup><mi>S</mi><mn>0</mn><mi>T</mi></msubsup><msub><mi>C</mi><mn>1</mn></msub><mo>,</mo></mrow>]]></math>id="icf0019"file="A2009100620280003C12.tif"wi="35"he="5"top="244"left="147"img-content="drawing"img-format="tif"orientation="portrait"inline="yes"/></maths><mathsid="math0012"num="0012"><math><![CDATA[<mrow><mi>B</mi><mo>=</mo><msubsup><mi>S</mi><mn>1</mn><mi>T</mi></msubsup><msub><mi>C</mi><mn>0</mn></msub><mo>-</mo><msubsup><mi>S</mi><mn>0</mn><mi>T</mi></msubsup><msub><mi>C</mi><mn>1</mn></msub><mo>,</mo></mrow>]]></math></maths>其中，k＝0，1，...，N-1，上標“T”表示轉置，上標“0”和“1”分別用于標示變換基前半部分和后半部分的子向量；步驟3.3，稀疏化近似步驟3.2所得和矩陣A，得到近似和矩陣Am；所述稀疏化具體方式為，保留和矩陣A中絕對值最大的2m個元素值，其它元素均置為0值；其中m取小于N×N/2的任意自然數；步驟3.4，根據步驟3.3所得近似和矩陣Am和步驟3.2所得差矩陣B對步驟3.1所得運算后譜線id="icf0021"file="A2009100620280004C1.tif"wi="5"he="5"top="102"left="64"img-content="drawing"img-format="tif"orientation="portrait"inline="yes"/>和id="icf0022"file="A2009100620280004C3.tif"wi="5"he="5"top="102"left="76"img-content="drawing"img-format="tif"orientation="portrait"inline="yes"/>進行轉換，實現(xiàn)MDCT譜到MDST譜的轉換，得到加窗信號的MDST譜；轉換公式如下其中id="icf0024"file="A2009100620280004C5.tif"wi="5"he="5"top="139"left="43"img-content="drawing"img-format="tif"orientation="portrait"inline="yes"/>和id="icf0025"file="A2009100620280004C6.tif"wi="5"he="5"top="139"left="54"img-content="drawing"img-format="tif"orientation="portrait"inline="yes"/>分別為步驟3.1中所得運算后譜線，B為步驟3.2中所得差矩陣，Am為步驟3.3中所得近似和矩陣；步驟4，依據步驟2中所得加窗信號的MDCT譜和步驟3中所得加窗信號的MDST譜，構造加窗信號的CMCLT譜；所述CMCLT譜是共軛調制復重疊變換譜，是以變換時所用窗函數為共軛窗函數的MDCT譜和MDST譜分別作為實部和虛部，從而定義MDCT的一種復數化擴展形式，構造的CMCLT譜如下式所示其中，id="icf0027"file="A2009100620280004C8.tif"wi="27"he="5"top="223"left="48"img-content="drawing"img-format="tif"orientation="portrait"inline="yes"/>分別為第i幀加窗信號的MDCT譜和MDST譜，xi為輸入信號的第i幀信號，ck、sk分別為MDCT變換基和MDST變換基，Wc和Ws分別為MDCT變換和MDST變換的窗函數矩陣表示，即<mathsid="math0013"num="0013"><math><![CDATA[<mrow><msub><mi>W</mi><mi>c</mi></msub><mo>=</mo><mfencedopen='('close=')'><mtable><mtr><mtd><msub><mi>W</mi><mn>0</mn></msub></mtd><mtd><mn>0</mn></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><msub><mi>W</mi><mn>1</mn></msub></mtd></mtr></mtable></mfenced><mo>,</mo><msub><mi>W</mi><mi>s</mi></msub><mo>=</mo><mfencedopen='('close=')'><mtable><mtr><mtd><msub><mi>W</mi><mn>1</mn></msub></mtd><mtd><mn>0</mn></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><msub><mrow><mo>-</mo><mi>W</mi></mrow><mn>0</mn></msub></mtd></mtr></mtable></mfenced></mrow>]]></math></maths>其中W0＝diag{w0(n)}，W1＝diag{w1(n)}，w0(n)和w1(n)分別為步驟1中窗函數wc(n)分解后的兩個N點窗函數；步驟5，根據加窗信號的CMCLT譜，提取加窗信號在CMCLT域的能量信息和相位信息；步驟6，根據步驟5所得能量信息及相位信息，在MDCT域上提取空間參數。2.如權利要求1所述的能量及相位調整方法，其特征在于所述相位信息的提取過程包括以下三個步驟，(1)將MDCT變換的窗函數表示為離散正弦變換的四型變換基^(")的線性組合，確定展開系數，即其中，a,為窗函數^(")以離散正弦變換的四型變換基L(")展開的系數，&取O，l,...,TV-l，且四型變換基^(")表示如下<formula>formulaseeoriginaldocumentpage5</formula>(2)對第/幀加窗信號的CMCLT譜的實部》'Ot)和虛部r(/t)做線性運算，構造乏U"及乏UW如下Z:(A:)-f乂A:)+上'(A:)(3)進行CMCLT譜的線性組合，將CMCLT譜轉換為DFT譜的形式，按照所述線性組合方式形成的新譜線i'^)如下式所示<formula>formulaseeoriginaldocumentpage6</formula>上式表現(xiàn)的關系展示了一個CMCLT譜的線性組合，等價于輸入信號的DFT譜線附加一個與輸入信號無關的相位偏移，且線性組合的形式依賴于四型變換基^(")展開的系數。3.用于實現(xiàn)權利要求1或2所述加窗信號MDCT域的能量及相位調整方法的裝置，其特征在于包含以下部分，MDCT加窗模塊，用于對輸入信號進行MDCT加窗得到加窗信號；共軛窗構造模塊，用于以MDCT加窗變換所用的窗函數為共軛窗函數構造MDST變換的窗函數；加窗MDST表示模塊，用于MDST變換的窗函數得到加窗信號MDST變換的形式化表達；加窗MDCT變換模塊，用于根據MDCT變換的窗函數得到加窗信號MDCT變換的形式化表達，根據該形式化表達對加窗信號進行MDCT變換，得到加窗信號的MDCT譜；加窗MDST構造模塊，用于由加窗信號的MDCT譜和加窗信號MDST變換的形式化表達構造加窗信號的MDST譜；加窗CMCLT構造模塊，用于依據加窗信號的MDCT譜和加窗信號的MDST譜，構造加窗信號的CMCLT譜；能量提取模塊，用于根據加窗信號的CMCLT譜，提取加窗信號在CMCLT域的能量信息；相位提取模塊，用于根據加窗信號的CMCLT譜，提取加窗信號在CMCLT域的相位信息；空間參數提取模塊，用于根據加窗信號在CMCLT域的能量信息和相位信息，在在MDCT域上提取空間參數；其中，輸入信號進入MDCT加窗模塊，MDCT加窗模塊分別連接加窗MDCT變換模塊和共軛窗構造模塊，MDCT加窗模塊通過連接輸出加窗信號進入加窗MDCT變換模塊，輸出MDCT加窗變換所用的窗函數到共軛窗構造模塊；共軛窗構造模塊與加窗MDST表示模塊相連，共軛窗構造模塊通過連接輸出MDST變換的窗函數到加窗MDST表示模塊；加窗MDST表示模塊和加窗MDCT變換模塊的輸出分別接入加窗MDST構造模塊，加窗MDST表示模塊通過連接輸出加窗信號MDST變換的形式化表達到加窗MDST構造模塊，加窗MDCT變換模塊輸出加窗信號的MDCT譜到加窗MDST構造模塊；加窗MDST構造模塊和加窗MDCT變換模塊的輸出分別接入加窗CMCLT構造模塊，加窗MDST構造模塊通過連接輸出加窗信號的MDST譜到加窗CMCLT構造模塊，加窗MDCT變換模塊通過連接輸出加窗信號的MDCT譜到加窗CMCLT構造模塊；加窗CMCLT構造模塊的輸出分別接入能量提取模塊和相位提取模塊，能量提取模塊所得能量信息和相位提取模塊所得相位信息最終進入空間參數提取模塊，空間參數提取模塊輸出在MDCT域上提取的空間參數。全文摘要本發(fā)明提出了加窗信號在MDCT域的能量及相位調整方法及其裝置，通過加窗信號的MDCT譜和MDST譜，構造加窗信號的CMCLT譜；根據加窗信號的CMCLT譜，提取加窗信號在CMCLT域的能量信息和相位信息，然后根據所得能量信息及相位信息，在MDCT域上提取空間參數。所述CMCLT譜是共軛調制復重疊變換譜，是以變換時所用窗函數為共軛窗函數的MDCT譜和MDST譜分別作為實部和虛部，從而定義MDCT的一種復數化擴展形式。應用本發(fā)明可統(tǒng)一空間音頻編碼系統(tǒng)中的時頻分析工具，有效提取多聲道信號的空間參數。統(tǒng)一時頻變換工具后，整個編碼端只需進行一次時頻變換，較大幅度降低了空間音頻編解碼系統(tǒng)的復雜度，從而優(yōu)化編解碼工作。文檔編號G10L19/00GK101552006SQ200910062028公開日2009年10月7日申請日期2009年5月12日優(yōu)先權日2009年5月12日發(fā)明者劉雨田,胡瑞敏,冰陳,琪陳,陳文琴,陳水仙申請人:武漢大學