专利名称：：一种基于先进音频编码器的量化器码率失真控制方法
技术领域：
：本发明涉及一种基于先进音频编码(AdvancedAudioCoding,AAC)的量化器码率失真控制方法。
背景技术：
：AAC是变换域编码的一种，它的目标是在一定的比特约束之下将失真的某种度量最优化。失真是由量化过程引起的，在量化过程中，码率失真控制为量化编码器的参数赋值，以获得准确的比特消耗数，从而能够将量化失真控制在感知容忍的范围内，以便保证编码的质量和效率。AAC编码使用了非均匀的量化器和变长编码，这使得量化算法的设计变得非常困难。首先，非均匀量化器无法根据感知容忍噪声来制定编码器参数，并且，因为变长编码的原因，比特消耗数只能通过量化结果计算获得。为控制量化失真和比特消耗，需要恰当地选择量化器的量化阶(scalefactor)等参数以保证码率失真控制有效i也控制量化失真和码率。双循环嵌套搜索(TLS)算法是目前使用得比较广泛的R-D控制算法，它是AAC标准推荐的算法，但这种算法在进行强制搜索最优解时会带来算法的不可实现性，它的复杂度为0((6(Tu广)(60种量化阶的可能性，11个码书，49个独立的子带)。为了降低R-D算法的复杂度，数种修正方法被提出，其中有设置最大循环次数，縮小量化阶搜索范围，使用统计线性模型预测比特消耗，使用上一帧的参数作为当前帧的初始值等等，但均未取得良好效果。
发明内容有鉴于此，本发明所要解决的技术问题在于，提供一种基于先进音频编码的量化器码率失真控制方法，该方法可以降低5马率失真控制的复杂度，提高编码器的编码质量和效率。为了实现上述目的，本发明采用的技术方案如下一种基于先进音频编码器的量化器码率失真控制方法，包含如下处理过程A、由待编码的码流获得子带可分配比特数的预测值，根据子带可分配比特数的预测值计算子带平均采样可分配比特数；B、由上述子带平均采样可分配比特数通过预定原则确定谷带量化阶选取系数；由待编码的码流的改进离散余弦变换谱和谷带量化阶选取系数计算得到谷带量化阶；C、以上述谷带量化阶为参数进行码率失真控制，得到最终参考比特消耗数。所述步骤B中的由待编码的码流的改进离散余弦变换谱和谷带量化阶选取系数计算得到谷带量化阶通过下式实现假设量化公式用如下函数表述c=G(x(vv)，D)，x(w)为待编码的码流的改进离散余弦变换谱，D为量阶，c为量化值。令c=v"/fe^^e，G"()为量化函数e()的逆函数，则求得的量阶D就是谷带量化阶，表达式为va〃e_y=g-1(x(w),va〃ey6a化)，特别的，如果在AAC中，表达式为max(lx(w)卜)其中，v。/—为谷带量化阶，x(vv)为待编码的码流的改进离散余弦变换谱，^//^^化为谷带量化阶选取系数。所述步骤B中确定谷带量化阶选取系数的预定原则为若to尸enSW，/ep<1.5，v"//^yZase=1;若1.5StoPenS薩j7/<2.5，=2;若2.5S6zYPenS謹p/ep<3.5，v"〃^y6a化=4;若3.5Sto尸enS扁//^<4.7，va/一6"化=7;若4.7S6"尸en&，/ep<6，va/fej^ose=12;若6^6"尸eriSfl，/，v"〃^Zos^=15;其中，w/一^e为谷带量化阶选取系数，toPeWamp为子带平均采样可分配比特数。所述子带平均采样可分配比特数通过下式计算得到6腐—p=sfbAvailableW+。.l5[l。g2(bW)-2],其中，toi^&邵、为子带平均采样可分配比特数，sfbAvailable("为当前子带可分配比特数，bw""为当前子带带宽。所述当前子带可分配比特数通过下式计算得到sfbAvailable("=sfbAvailable(A:-1)-sfbCost(A:—1)+0.7sfbBits(A0sfbAvailable(A:5=max(0.95sfb6its(i),sfbAvailable〖A:》，其中，sfbAvailable(fc)为当前子带可分配比特数，sfbBits("为所述步骤A中得到的子带可分配比特数的预测值，sfbCost(ll)为上一子带的参考比特消耗数，所述参考比特消耗数通过谷带量化阶按预定方式计算得到。所述预定方式是利用谷带量化阶对子带进行量化，得到整数量化谱，并以整数量化谱为索引检索霍夫曼码书。所述的量化器码率失真控制方法，步骤C中还包括对谷带量化阶的修正处理71)通过下式计算得到量化阶修正幅度因子y^ctorC7a"ge=int+0.5其中>T^cfo厂c^a"gg为:、0.19,化阶修正幅度因子，int为取整运算符，to&nfomp为子带平均采样可分配比特数，toPen&附^/、,为子带平均采样参考比特消耗数，由下式计算得到=￡^^，其中，sfbbit("为当前子带参考比特消耗数，bw彿("为当前子带带宽；72)利用量化阶修正幅度因子按预定原则对谷带量化阶进行修正。所述步骤72)中谷带量化阶修正的预定原则为-若量化阶修正幅度因子^cC^"ge小于-4，且谷带量化阶选取系数为1，放弃当前子带信息，不进行编码；若量化阶修正幅度因子/"cC^"ge大于等于-4，小于0，且谷带量化阶选取系数为l，子带带宽大于4，将谷带量化阶增加l内c^C^"ge卜若量化阶修正幅度因子/""orCA朋ge大于等于0，小于等于1，或者量化阶修正幅度因子力ctoK^朋ge小于0且子带带宽为4，谷带量化阶不变；若量化阶修正幅度因子/actorC/^"ge大于1，或者/actorC/朋ge小于0且子带带宽大于4，进行窗口搜索，选择子带内加权量化误差最小的谷带量化阶作为当前子带的谷带量化阶。当量化阶修正采用窗口搜索方式时，包括如下步骤91)确定窗口范围以[sear/^啤,化arc緣e/+1,…,5earc緣e/+化arcWVw附]作为量化阶搜索的窗口范围，其中，=max(化arc/ziS啤，va〃^y—附ox5te/+1)<formula>formulaseeoriginaldocumentpage9</formula>，式中，iW^y为谷带量化阶，x(w)为待编码的码流的改进离散余弦变换谱，w/^6c^为谷带量化阶选取系数，/actorOz朋ge为量化阶修正幅度因子，bw("为子带k的带宽；92)计算子带内加权量化误差-sfbhigh(A)Error—(&)=艺rate(w)(0.06416aO)4+0.5948a(3+1.775a(w)2》其中，Error^W为子带内加权量化误差，w是谱序号，sfblow(^，sfbhigh^)分别为子带A的上界和下界，ot(w)为相对取整误差，由下式计算得到-mt16'saclefactor(A)—x(w)卩x216'saclefactor(&)x(w)4x2其中lx(w)l:x2—f,是量化后未取整的谱，saclefactorW为子带k进行窗口搜索时的候选量化阶；rate(w)为加权系数，由下式计算得到x(w)2x(w)4.rate(vv)=0.1sfbhigh(/t)w=sfblow(A;)+0,993)选择子带内加权量化误差最小的谷带量化阶作为当前子带的谷带量化阶。所述量化器码率失真控制方法，步骤C中的子带参考比特消耗数可以直接用谷带量化阶，或者更精确地用通过对谷带量化阶调整得到的量化阶对编码子带进行量化(如果没有调整，无需重新量化)，得到整数量化谱，并以整数量化谱为索引检索霍夫曼码书得到的。所述的量化器码率失真控制方法，还包括利用幅度畸变因子对量化阶的再修正处理若i<0.8，saclefactor=saclefactor(A:)+1;若0.8S"1.25，不做修正；若i>1.25，saclefactor=saclefactor(A:)-1;其中，saclefactor为再修正后的量化阶，saclefactor(A:)为再修正前的量化阶，R为幅度畸变因子，通过下式计算得到<formula>formulaseeoriginaldocumentpage10</formula>其中，sfblow("，sfbhigh("分别为子带A:的上界和下界，x(w)为待编码的码流的改进离散余弦变换谱，vv为谱序号，m"g/c为量化跳跃点，取值0.4054。本发明的有益效果在于，改进了量化阶的确定方法，以此方法确定的量化阶作为码率失真控制方法的参数，有效地降低了码率失真控制的复杂度，并因此降低了量化编码时的运算量，从而提高了编码器的编码质量和效率。图1是运用本发明所述方法的MegalAAC编码器的结构框图；图2是本发明具体实施方式的流程图3是图1所示MegalAAC编码器与另两种现有的编码器的ODG指标比较示意图4是图1所示MegalAAC编码器与另两种现有的编码器的丽R指标比较示意图5是图1所示MegalAAC编码器与另两种现有的编码器的0DG分布示意图6是图1所示MegalAAC编码器与另两种现有的编码器的丽R分布示意图。具体实施例方式下面对照附图结合具体实施方式对本发明作详细说明。图1是运用本发明所述方法的MegalAAC编码器的结构框图，主要包括心理声学模型处理模块、谱处理模块、量化编码模块和码流合成模块，其中，量化编码模块中采用了本发明所述的量化器码率失真控制方法，以便降低量化编码时的运算量。本发明具体实施方式的所述方法的流程参看图2，包含如下步骤1.获取谷带量化阶进行量化编码处理，需要先确定量化阶；本发明确定的量化阶，称为谷带量化阶，其计算步骤如下la)计算子带平均采样可分配比特数对于待编码的码流的编码处理是一个连续的过程，通常以一帧为单位，并将每一帧数据分为若干子带进行处理。设子带可分配比特数的预测值为sfbBits("，该预测值可以通过AAC标准提供的传统心理声学模型算法由待编码的码流计算得出。一般公式为子带可分配比特数的预测值=子带感知熵x当前帧可用的比特数/所有子带感知熵和。其中，如编码是定比特率编码(CBR)，那么当前帧可用的比特数是一个定值，等于比特率X1024/采样率；如果是随着使用情况变化，那么就是可变比特率编码(VBR)，这种情况下的当前帧可用的比特数一般由帧间比特控制算法提供。将该预测值减去边带信息所占据的一部分比特数比例，该比例通常在0.05到0.5之间，本实施例取0.3。若上一子带的参考比特消耗数为sfbCost(A-l)，上一子带可分配比特数的预测值为sfbAvailable(;t-l)，那么当前子带可分配比特数为<formula>formulaseeoriginaldocumentpage11</formula>式(1)中，0.95是比特控制的一个系数，可取0.5到1之间的数，本实施例取(下式例取0.95;再利用子带带宽bw^修正获得子带平均采样可分配比特数，如<formula>formulaseeoriginaldocumentpage11</formula>lb)谷带量化阶选取系数的确定谷带量化阶选取系数按如下原则确定<table>tableseeoriginaldocumentpage12</column></row><table>其中，将^//^6。^=1时的谷带量化阶称为第一谷带量化阶，等于2时的谷带量化阶称为第二谷带量化阶，等于4时的谷带量化阶称为第四谷带量化阶，余者类推。做//^6"化为谷带量化阶选取系数。在本实施例应用于AAC中，根据AAC的量化公式，以谷带量化阶选取系数作为量化的结果，由待编码的码流的改进离散余弦变换(MDCT)谱反推计算出谷带量化阶，其具体通过下式得到max(lx(w)卜)(3)其中，v"/—为谷带量化阶，x(w)为待编码的码流的改进离散余弦变换(MDCT)谱，当然，.在本发明应用于其他的音频编码器时，由于量化公式的不同，式(3)也相应有所区别。2.利用谷带量化阶进行量化，得到参考比特消耗数2a)使用谷带量化阶进行量化，获得参考量化谱AAC中的参考量化谱是由式(4)产生，其中scalefactor[/t]为子带A;的量化阶，m"g/c为量化跳跃点，AAC的标准推荐为0.4054，x(w)为MDCT谱，w为谱序号，满足条件sfbLow("《w《sfbHigh("，sfbLow(A)和sfbfflgh("分别是子带it的下界和上界。|2■--scalefactor[A:]X("VV)4X216+m喂'c(4)以式(3)获得的谷带量化阶^//^作为当前子带的量化阶，利用式(4)获得当前子带的参考量化谱。2b)子带参考比特消耗数的计算子带的参考比特消耗数需要使用参考量化谱查找码书来获得。本实施例中的码书，使用霍夫曼码书，AAC中的霍夫曼码书一共有13个，其规格如表i，其中o号码书用来表示量化谱数据全为0，1到11号码书用来对量化谱进行编码，縮放因子码书用来对縮放因子进行编码。表lAAC中的霍夫曼码书<table>tableseeoriginaldocumentpage13</column></row><table>计算量化阶的消耗码长，具体过程为获得上一子带的量化阶scalefactor[Ll]和当前子带的量化阶scalefactor[it]，量化阶使用差分编码，HUf_SfLen是縮放因子的码长表，由下式确定量化阶的消耗码长=scalefactor[A:_1]_scalefactor[A:]⑤sfLen("=Huf—sfLe争/+60]量化谱的消耗码长的计算过程为首先求得子带中整数量化谱的最大绝对值、^依次将子带中相邻的四个整数量化谱组成四维矢量["Z2^,4]作为索引去查找，累加求得不同码书编码所需码长的总数，其中，有符号码书的索引要加上最大绝对值作为偏移，而无符号码书的索引则加上1。检索方法为corfe丄e"l=HufJLenl[d+1]|>2+1]|>3+1]["+1]cocfe丄e"2=Huf—Len2[fl+1]|>2+1][。+l][f4+1]coffe丄ew5=Huf—Len5[d+4]|>2+4]+Huf—Len5[f3+4][,4+4]a^e丄e"6=Huf_Len6|>l+4][f2+4]+Huf—Len6[,3+4][f4+4]n=|,l|,,2,|,f3=|,3|,M=|M|Ze"=!W+!!f2+!W+!!M^e￡e4=Huf—Len4[fl]02]03][M]+WgwZewcocfe丄ew7=HufJLen7[f1]02]+Huf—Len7[f3]04]+w'gw丄ewcwfe丄ew8=Huf—Len8[>l]|>2]+Huf—Len8[。][M]+s/gw丄ewcocfe丄e"9=Huf—Len9|>l][/2]+Huf—Len9[f3][f4]+cocfeZew10=Huf_LenlO[Yl][>2]+Huf—LenlO[。][M]+cocfe丄e"l1=Huf—Lenll[>l][f2]+Huf—Lenll[f3][74]+w'g"丄e"其中符号！！的运算原则为若r0，！!Z检索码书的流程根据/_判断如下1，,ax=1，检索1，2，3，4，5,6，7，8，9，10，11号码书"ax=2，检索3，4，5，6，7，8，9，10，11号码书2</maxC]4,检索5，6，7，8，9，10，11号码书4<、ax□，检索7，8，9，10，11号码书7<,max012，检索9，10，11号码书I2<、ax□"，检索11号码书16￡L，检索ll号码书，启动溢出保护机制(6)7=0，！!,=0。其中，溢出保护机制，在U16时先以16替代大于16的整数谱d乍为索引进行检索，然后在这个霍夫曼码字后面，为等于16的索引依次添加逃避序列。逃避序列(escapesequence)由以下三部分组成[escape_prefix][escape—separator][escape—word]其中[escape—prefix]是；v个'1，，[escape—separator]是一个'0，，[escape_word]是一个w+4比特的无符号整型数，,满足f=+[escape—word],增力卩比年寺数为2int[log2(f+1)+l]-5。检索码长码书，得到第《子带整数量化谱使用m号码书的参考比特消耗值，将其记录在参考比特消耗矩阵reffiits(m,"中，没有计算到的则在矩阵相应位置记为65535。当前子带k的参考比特消耗数为sfbbit("=sfLen(A:)+min(reffiits(m，A:))，m=1,2...,11(7)为了进一步提高码率失真控制的效果，本发明具体实施方式还包括对量化阶的修正处理步骤，在此，将前述步骤中出现的没有经过修正处理的谷带量化阶称为初始谷带量化阶，而将此后出现的经过修正处理的量化阶称为最终谷带量化阶。2c)量化阶修正幅度因子计算计算子带平均采样参考比特消耗数sfbbit("计算量化阶修正幅度因子(8)_/^ctorC%a"ge=int0.19+0.5(9)int为取整运算符，to/^&^^、为子带平均采样可分配比特数。2d)谷带量化阶修正判断1)若/"c^C/^"ge小于-4，并且初始谷带量化阶是第一谷带量化阶，那么当前子带信息被放弃，不进行编码；2)若力c^C/w"ge大于等于-4，小于O，并且初始谷带量化阶是第一谷带量化阶且子带带宽大于4，那么，在初始谷带量化阶的基础上增加|/aC^O^ge|，以此作为最终谷带量化阶；3)若/actorCA朋ge大于等于0，小于等于l，或者/flCO^"ge小于0且带宽为4，那么以初始谷带量化阶作为最终谷带量化阶；4)若/actorC/w"ge大于1，或者/fl"orOw"ge小于0且带宽大于4，那么进行窗口搜索，修正判断决定进行窗口搜索时，首先确定窗口范围，过程如下searc/ziS啤=ra/——y^ctorC7a"ge+int+0.(10)化flrc/ziS啤=max—"rcM啤，va〃ey—附ox5V^p+1)其中mw&印起到量化谱最大值限制的作用，由下式得到:p16,,xS鄉=ylog22435max'(lx(w)卩)化arc/L/Vw附=min24bw()t)十mm乂乂6，afo(j^加/^awg^)cos(ii)A2-2001450,(12)以[searc/^e/,sewc/^e;+1,，saan^iSVe;+sea/rAM/附]作为量<七卩介，叟索的齒口范围，计算带内加权量化误差，过程为计算相对取整余差inta(w)=—l|一--saclefactor(t)||_---(13)i,2--saclefactor("Ix(w)卜x216其中|X(w)|*x2-》她—"是量化后未取整的谱，saclefactor("为子带k搜索时的候选量化阶。加权系数为—=0.lstbhigy2+，4(14)带内加权误差为sfbhigh(A)Error^(A:)=艺rateO)(0.06416otO)4+0.5948a(w)3+1.775a(2)其中sfblow(/t)，sfbhigh(/t)分别为子带it的上界和下界，x(w)为MDCT谱，w为谱序号。对于需要进行窗口范围量化阶搜索的子带，进行sewcM^n+1次带内加权量化误差计算后，选择其中带内加权量化误差最小的量化阶作为当前子带的最终谷带量化阶。2e)经过量化阶修正处理后，由初始谷带量化阶得到最终谷带量化阶，使用最终谷带量化阶利用式(4)进行最终量化(如果没有修正，无需重新量化)，计算得到最终的子带参考比特消耗数。2f)对最终谷带量化阶进行再修正为了解码器更好的反量化，还可以对最终谷带量化阶进行再修正，再修正后的量阶可用于解码器的反量化。再修正处理首先对当前子带的量化取整谱计算其幅度畸变因子<formula>formulaseeoriginaldocumentpage16</formula>若0.8S"1.25，不做修正；若i〉1.25，另卩么saclefactor=sacIefactor(A)—1。其中saclefactor'(it)为再修正后的量化阶，saclefactor("为再修正前的量化阶，sfblow("，sfbhigh("分别为子带A的上界和下界，x(w)为MDCT谱，w为谱序号。确定了初始谷带量化阶后，即可以该量化阶为参数进行码率失真控制。一方面，由初始谷带量化阶获得子带参考比特消耗数，可直接作为后续编码步骤的输入，完成熵编码和码流合成，最终得到编码码流。另一方面，更精确的，子带参考比特消耗数可作为反馈提供给量化阶确定过程而对初始谷带量化阶进行修正，以得到更好的码率失真控制效果，从而最终提升整个编码器的整体编码效率与质量。以上公式和数据均是通过实验结果或者是综合待求参数的相关影响因素拟合所给出的较优值，本发明并不限于仅使用上述具体数值。有关运算复杂度评估通过比较未使用本发明的FreeAdvancedAudioCoder(FAAC)和使用了本发明的MegalAACEncoder(Megal)，在对44100Hz采样，16位量化的立体声音频进行平均比特率编码的情况下进行，参考指标为每秒百万操作数，其中FAAC是德国FraunhoferIIS开发的AAC编码器，FraunhoferIIS是Mpeg标准的主要制定者，其FAAC编码器是AAC标准的验证编码器。表2编码器的运算复杂度<table>tableseeoriginaldocumentpage17</column></row><table>算消耗比特数时，仅需调用码书4.35次,相比较FAAC需要61.37次。Megal的量化编码算法的计算量为12.8M0PS，其总计算量为21.2M0PS，相比较FAAC分别为35M0PS和90.8M0PS，分别降低63.4%和76.7%。编码器的音质评估使用EAQUAL1.3，这是使用感知音频客观评价标准PEAQ的客观评价程序，PEAQ提供的音质评价指标名称及其描述见表3。表3eaqual软件输出的评价指标及意义<table>tableseeoriginaldocumentpage18</column></row><table>这里选择总体指标(ODG)和两个重要的单项指标(BandwidthTest和丽R)作为主要参考指标。音质评估使用三种编码器参照进行，对它们进行横向评估，这三个编码器分别为本发明的Megal，现有技术的NCTUAACEncoder(以下简称NCTU)和FAAC。其中NCTU是台湾交通大学感知音频小组开发的AAC编码器，音质评估源使用美国惠威公司提供的音响试音碟的第一张和第二张，剔除了重复的曲目，选取了其中37支音乐剪辑，这些剪辑包含了乐曲的基本类型，它们的名称以及描述见表4。表4测试曲目<table>tableseeoriginaldocumentpage18</column></row><table><table>tableseeoriginaldocumentpage19</column></row><table>测试结果见表5。表5测试结果<table>tableseeoriginaldocumentpage20</column></row><table><table>tableseeoriginaldocumentpage21</column></row><table>从以上测试数据可以得到，NCTU的平均0DG相对FAAC提高0.163，而使用本发明的Megal的平均0DG相对NCTU又提高了0.187;NCTU的平均丽R相对FAAC降低了1.06dB，而使用本发明的Megal的平均丽R相对NCTU又降低了L08dB。类似的结论可以在关于MegalAAC编码器与另两种编码器的ODG与丽R指标比较和分布的图3、图4、图5和图6中得到。计算量评估与音质客观评价都说明，本发明可以使AAC编码器运算量大大降低，而音质得到大幅提高，整个编码器的质量和效率具有显著的提升。权利要求1.一种基于先进音频编码器的量化器码率失真控制方法，其特征在于，包含如下处理过程A、由待编码的码流获得子带可分配比特数的预测值，根据子带可分配比特数的预测值计算子带平均采样可分配比特数；B、由上述子带平均采样可分配比特数通过预定原则确定谷带量化阶选取系数；由待编码的码流的改进离散余弦变换谱和谷带量化阶选取系数计算得到谷带量化阶；C、以上述谷带量化阶为参数进行码率失真控制，得到最终参考比特消耗数。2.如权利要求1所述的量化器码率失真控制方法，其特征在于，所述步骤B中的由待编码的码流的改进离散余弦变换谱和谷带量化阶选取系数计算得到谷带量化阶通过下式实现max(lx(w)卜)其中，^//^为谷带量化阶，x(w)为待编码的码流的改进离散余弦变换谱，w/^y^e为谷带量化阶选取系数。3.如权利要求1或2所述的量化器码率失真控制方法，其特征在于，所述步骤B中确定谷带量化阶选取系数的预定原则为若6"尸enSa，/<1.5，va/一6a化=1，若1.5S6"尸^S画—p<2.5，va/一Z)^e二2若2.5S6"尸e/;Sfl，/ep<3.5，=4若3.5S<4.7，va//ej^xye=7若4.7S6"尸en&，/ep<6，va/一6(xse=12;若6SZzY尸enSWmp/ep，va//ej^as^=15;其中，^//^Z^e为谷带量化阶选取系数，^尸^&，罕为子带平均采样可分配比特数。4.如权利要求3所述的量化器码率失真控制方法，其特征在于，所述子带平均采样可分配比特数通过下式计算得到<formula>formulaseeoriginaldocumentpage2</formula>其中，toP^Sa，/为子带平均采样可分配比特数，sfbAvailable("为当前子带可分配比特数，b("为当前子带带宽。5.如权利要求4所述的量化器码率失真控制方法，其特征在于，所述当前子带可分配比特数通过下式计算得到sfbAvailable("=sfbAvailable(A:-1)-sfbCost(A-1)+0.7sfbBits(A:)sfbAvailable(A^=max(0.95sfbBits(A:),sfbAvailable&))，其中，sfbAvailable("为当前子带可分配比特数，sfbBits("为所述步骤A中得到的子带可分配比特数的预测值，sfbCost(H)为上一子带的参考比特消耗数，所述参考比特消耗数通过谷带量化阶按预定方式计算得到。6.如权利要求5所述的量化器码率失真控制方法，其特征在于，所述预定方式是利用谷带量化阶对子带进行量化，得到整数量化谱，并以整数量化谱为索弓I检索霍夫曼码书。7.如权利要求6所述的量化器码率失真控制方法，其特征在于，所述步骤C中还包括对谷带量化阶的修正处理71)通过下式计算得到量化阶修正幅度因子<formula>formulaseeoriginaldocumentpage3</formula>其中，々ctorC7j朋ge为:化阶修正幅度因子，int为取整运算符，to/^&zm;^p为子带平均采样可分配比特数，Wi^&m/^^为子带平均采样参考比特消耗数，由下式计算得到to尸enfom;^,=!￡^1，其中，sfbb啡)为当前子带参考比特消耗数，b("为当前子带带宽；72)利用量化阶修正幅度因子按预定原则对谷带量化阶进行修正。8.如权利要求7所述的量化器码率失真控制方法，其特征在于，所述步骤72)中谷带量化阶修正的预定原则为若量化阶修正幅度因子^c^CA朋^小于-4，且谷带量化阶选取系数为放弃当前子带信息，不进行编码；若量化阶修正幅度因子/"c^C^wge大于等于-4，小于0，且谷带量化阶选取系数为l，子带带宽大于4，将谷带量化阶增加l/actorC^"gel;若量化阶修正幅度因子/"c^C/^ge大于等于0，小于等于1，或者量化阶修正幅度因子/"ctorC/^zge小于0且子带带宽为4，谷带量化阶不变；若量化阶修正幅度因子/"ctorC7^wge大于1，或者/actorOw"ge小于0且子带带宽大于4，进行窗口搜索，选择子带内加权量化误差最小的谷带量化阶作为当前子带的谷带量化阶。9.如权利要求8所述的量化器码率失真控制方法，其特征在于，当量化阶修正采用窗口搜索方式时，包括如下步骤91)确定窗口范围以[舰rc緣e;,s匿c雄e/+1,…，s麼c緣e;+,rc研臓]作为量化阶搜索的窗口范围，其中，<formula>formulaseeoriginaldocumentpage4</formula>，式中，^//^为谷带量化阶，x(w)为待编码的码流的改进离散余弦变换谱，ra〃e_yZ^e为谷带量化阶选取系数，/actorC/2朋ge为量化阶修正幅度因子，bw(A:)为子带k的带宽；92)计算子带内加权量化误差<formula>formulaseeoriginaldocumentpage4</formula>其中，Error^("为子带内加权量化误差，w是谱序号，sfblowW,sfbhigh(。分别为子带A的上界和下界，(x(vv)为相对取整误差，由下式计算得到<formula>formulaseeoriginaldocumentpage4</formula>其中lx(v^x2—)是量化后未取整的谱，saclefactor("为子带k进行窗口搜索时的候选量化阶；rate(w)为加权系数，由下式计算得到<formula>formulaseeoriginaldocumentpage4</formula>93)选择子带内加权量化误差最小的谷带量化阶作为当前子带的谷带量化阶。10.如权利要求9所述的量化器码率失真控制方法，其特征在于，还包括利用幅度畸变因子对量化阶的再修正处理若/<0.8，saclefactor=saclefactor("+1;若0.8^^1.25，不估嫌正；若>1.25，saclefactor=saclefactor("-1;其中，saclefactor'(/t)为再修正后的量化阶，saclefactor(/t)为再修正前的量化阶，R为幅度畸变因子，通过下式计算得到<formula>formulaseeoriginaldocumentpage5</formula><formula>formulaseeoriginaldocumentpage5</formula><formula>formulaseeoriginaldocumentpage5</formula>其中，sfbkm("，sfbhigh(yt)分别为子带/t的上界和下界，x(w)为待编码的码流的改进离散余弦变换谱，w为谱序号，m"g,c为量化跳跃点，取值0.4054。全文摘要本发明公开了一种基于先进音频编码器的量化器码率失真控制方法，包含如下处理过程A.由待编码的码流获得子带可分配比特数的预测值，根据子带可分配比特数的预测值计算子带平均采样可分配比特数；B.由上述子带平均采样可分配比特数通过预定原则确定谷带量化阶选取系数；由待编码的码流的改进离散余弦变换谱和谷带量化阶选取系数计算得到谷带量化阶；C.以上述谷带量化阶为参数进行码率失真控制，得到最终参考比特消耗数。本发明有效地降低了码率失真控制的复杂度，并因此降低了量化编码时的运算量，从而提高了编码器的编码质量和效率。文档编号H03M7/40GK101308661SQ20071014845公开日2008年11月19日申请日期2007年8月21日优先权日2007年5月16日发明者晟吴,邱小军,黎家力申请人:中兴通讯股份有限公司