WO2012094952A1

WO2012094952A1 - 一种信号处理方法和装置

Info

Publication number: WO2012094952A1
Application number: PCT/CN2011/085197
Authority: WO
Inventors: 张德明; 李海婷; 阿里斯⋅塔勒布; 许剑峰
Original assignee: Huawei Technologies Co Ltd
Current assignee: Huawei Technologies Co Ltd
Priority date: 2011-01-10
Filing date: 2011-12-31
Publication date: 2012-07-19
Anticipated expiration: 2013-07-10
Also published as: EP2664995A4; EP2664995A1; KR20130116904A; US9792257B2; EP3518121C0; PL3518121T3; JP5783395B2; EP4550325A3; US9996503B2; US9519619B2; US20170075860A1; KR101627900B1; US20170364479A1; US20130304784A1; ES3026524T3; EP3518121B1; EP3518121A1; EP4550325A2; JP2014503093A; CN102592601A

Description

一一一种信号处理方法和装置

本申请要求于 2011 年 01 月 10 日提交中国专利局、申请号为 201110004032.5、发明名称为 "一种信号处理方法和装置" 的中国专利申请的优先权，其全部内容通过引用结合在本申请中。技术领域

本发明涉及数字信号处理技术领域，尤其是一种信号处理方法和装置。背景技术

在数字通信领域，语音、图像、音频、视频的传输有着非常广泛的应用需求，如手机通话、音视频会议、广播电视、多媒体娱乐等。为了降低音视频信号存储或者传输过程中占用的资源，音视频压缩编码技术应运而生。音视频压缩编码技术的发展中涌现出了很多不同的技术分支，其中将信号由时域变换到频域后再进行编码处理的技术，又称为变换域编码技术由于具有很好的压缩特性，得到了非常广泛的应用。

变换域编码技术中将信号由时域变换到频域的方法有很多种，其中傅立叶变换 ( Discrete Fourier transform, DFT )、离散余弦变换 ( Discrete Cosine Transform, DCT )、离散正弦变换（ Discrete sine transform, DST )以及修正离散余弦变换（ Modified Discrete Cosine Transform, MDCT )等时频变换有着广泛的应用，在频谱分析、图像编码和语音编码等领域的应用尤其广泛。信号经过时频变换后可以使用量化技术进行压缩编码，也可以使用其他一些参数音频编码的方法进行编码，从而达到数据压缩的目的。

然而，发明人发现，直接根据变换公式进行 DCT-IV或 MDCT正变换与逆变换会导致很高的计算复杂度和存储量，提供一种低存储量的时域频域变换方法成为一种迫切的需要。发明内容

本发明实施例的目的在于提供一种数据处理方法，减少音视频编码过程中时域频域变换处理的存储量。 - - 根据本发明的一实施例，一种数据处理方法，包括：

对输入数据进行重组，得到重组后的数据；

使用一个对称的旋转因子对重组后的数据进行预旋转，旋转因子为 ^a-^ , P = 0, ..., L I 2 - l ^ _a为常数；

对予贞旋转处理后的数据进行 2个点的†夬速傅里 ( Fast Fourier Transform, FFT ) 变换， L为输入数据的长度；

使用一个对称的旋转因子对 FFT变换后的数据进行后旋转处理，旋转因子为

H , _b为常数；获得输出数据。根据本发明的另一实施例，一种时域到频域的信号处理方法，包括：对时域数据进行预处理，得到预处理后数据；

使用旋转因子 ^a ' ^ 对所述预处理后的数据进行预旋转处理；

对所述预旋转处理后的数据进行 N/4点的快速傅里叶变换；

使用旋转因子 ^b ' ^^k+°'⁵对所述离散傅里叶变换后的数据进行后旋转处理，获得频域数据；

其中，在获得频域数据之前还包括：利用固定旋转补偿因子进定旋转补偿的步骤；所述 a、 b为常数，所述 N为时域数据长度，所述 ^{= e N}。根据本发明的另一实施例，一种频域到时域的信号处理方法，包括：对频域数据进行进行预处理，得到预处理后的数据；

使用旋转因子 ^e ' ^W^⁺⁰'⁵对所述预处理后的数据进行预旋转处理；

对所述预旋转处理后的数据进行 N/4点的快速傅里叶变换；

使用旋转因子 ^d · 。·⁵对所述快速傅里叶变换后的数据进行后旋转处理；对所述后旋转处理后的数据进行后处理，获得时域数据；

其中，在获得时域数据之前还包括：利用固定旋转补偿因子进行固定旋转补偿的步骤；所述 c、 d为常数，所述 N为两倍的频域数据长度，所述 ^{] N}。根据本发明的另一实施例，一种信号处理装置，包括： - - 重组单元，用于对输入数据进行重组，得到重组后的数据；

预旋转单元，用于使用一个对称的旋转因子对重组后的数据进行预旋转，旋转因子为 a ^¹, Ρ = 0，···，^/2_1, _a为常数；

变换单元，用于对预旋转处理后的数据进行个点的快速傅里（Fast Fourier Transform, FFT) 变换， L为输入数据的长度；

后旋转单元，用于使用一个对称的旋转因子对 FFT变换后的数据进行后旋转处理，旋转因子为 b'^^ , ? = 0，…， /2-1, b为常数；

输出单元，用于获得输出数据。根据本发明的另一实施例，一种时域到频域的信号处理装置，包括：预处理单元，用于对时域数据进行预处理，得到预处理后的数据；预旋转单元，用于使用旋转因子对所述预处理后的数据进行预旋转处理；

变换单元，用于对所述预旋转处理后的数据进行 N/4点的快速傅里叶变换；后旋转单元，用于使用旋转因子 b'^^+o⁵对所述离散傅里叶变换后的数据进行后旋转处理，获得频域数据；其中，还包括：

固定补偿单元，用于利用固定旋转补偿因子进₂行固定旋转补偿处理；所述 a、 b为常数，所述 N为时域数据长度，所述 = 。根据本发明的另一实施例，一种频域到时域的信号处理装置，包括：预处理单元，用于对频域数据进行预处理，得到预处理后的数据；预旋转单元，用于使用旋转因子 ^e · ^W^⁺⁰'⁵对所述预处理后的数据进行预旋转处理；

变换单元，对所述预旋转处理后的数据进行 N/4点的快速傅里叶变换；后旋转单元，使用旋转因子 d'^w^⁵对所述快速傅里叶变换后的数据进行后旋转处理；其中，还包括：

固定补偿单元，用于利用固定旋转补偿因子进行固定旋转补偿处理；所述 c、 d为常数，所述 N为两倍的频域数据长度，所述 =_e— ^。一一本发明实施例，预旋转和后旋转步骤使用的旋转因子具有对称性，可以降低数据存储量。同时，使用 FFT可以加快变换的速度，降低其计算复杂度。

附图说明

为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提下，还可以根据这些附图获得其他的附图。

图 1为本发明提供的时域到频域的 DCT-IV变换方法一个实施例的流程示意图；

图 2为本发明提供的时域到频域的 DCT-IV变换方法另一个实施例的流程示意图；

图 3为本发明提供的时域到频域的 DCT-IV变换方法另一个实施例的流程示意图；

图 4为本发明提供的时域到频域的 MDCT变换方法一个实施例的流程示意图；

图 5为本发明提供的频域到时域的 MDCT变换方法一个实施例的流程示意图；

图 6为本发明提供的时域到频域的 MDCT变换方法另一个实施例的流程示意图；

图 7为本发明提供的频域到时域的 MDCT变换方法另一个实施例的流程示意图；

图 8为本发明提供的信号处理装置一个实施例的结构示意图；

图 9为本发明提供的时域到频域的信号处理装置一个实施例的结构示意图；

图 10为本发明提供的频域到时域的信号处理装置一个实施例的结构示意图。具体实施方式

下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。

数字信号处理领域，音频编解码器、视频编解码器广泛应用于各种电子设备中，例如：移动电话，无线装置，个人数据助理（PDA ), 手持式或便携式计算机， GPS接收机 /导航器，照相机，音频 /视频播放器，摄像机，录像机，监控设备等。通常，这类电子设备中包括音频编码器或音频解码器，音频编码器或者解码器可以直接由数字电路或芯片例如 DSP ( digital signal processor )实现，或者由软件代码驱动处理器执行软件代码中的流程而实现。

例如一种音频编码器，首先对输入信号进行分帧处理，得到 20ms—帧的时域数据；然后对时域数据进行加窗处理，得到加窗后的信号；对加窗后的时域信号进行频域变换，例如 MDCT变换或者 DCT-IV型变换，将信号由时域变换到频域；再对频域信号进行分带处理，得到分带以后的频域信号；然后对每一个子带信号的能量进行计算，并对子带能量进行量化编码，传输到解码端；接下来根据量化后的子带能量值进行基于听觉掩蔽效应的自适应比特分配，得到各个子带用于量化编码的比特数；最后对各个子带内的频点进行归一化处理，根据所分配的编码比特数，使用矢量量化技术对子带内归一化处理后的频点进行矢量量化，得到矢量量化的码书索引，编码后传输到解码端。解码端接收到由编码端传输出来的压缩码流以后，才艮据相应的解码步骤从码流中查找出各子带信号能量的码书索引，并得到各子带信号能量的量化值；根据这些量化值釆用和编码端一致的比特分配技术，得到各个子带所分配的比特数；根据各个子带所分配的比特数，以及从码流中获取的各个子带的矢量量化的码书索引，得到各个子带量化后的归一化频域系数；根据各子带信号能量的量化值将各个子带量化后的归一化频域系数进行去归一化处理，得到完整的频域信号；对解码的到的频域信号进行编码端所用变换相对应的逆变换，将信号由频域变换到时域，对时域信号进行后处理后得到合成信号，即为解码端输出信号。时域到频域的信号处理方法又可称作正变换，频域到时域的信号处理方法又可称一一作逆变换。

DCT作为一种空间变换，其最大特点是具有能量紧致性，这就使得以 DCT 为基础的编码系统会具有很好的压缩性能。第四型 DCT(DCT-I V)经常用于音频和视频数据压缩。 DCT-IV变换的公式为：

HU-¹其中， k为 0到 L-l间的整数。可见，直接根据变换公式进行 DCT-IV 正变换与逆变换会导致很高的计算复杂度和存储量。由于 DCT-IV变换在实时通信领域特别是音频编码中的广泛应用，减少 DCT-V变换方法的存储量成为一种迫切的需要。

参考图 1, 本发明一个实施例提供的一种信号处理方法，用于实现编码过程中的时域到频域的 DCT-IV变换，减少变换过程的存储量，该方法包括如下步骤：

S101: 对时域数据进行重组，得到重组后的数据。

假设为需要进行 DCT-IV型变换的数据，该数据可以是经过加窗等预处理步骤后的数据。对数据进行重组，得到重组后的数据 ^z (^：

P) = ^P) +j-x(L-\-2pXp = 0X2,- ₂-\

S102: 使用一个对称的旋转因子对重组后的数据进行预旋转，旋转因子为 - = 0,l,2,…， /2— 1, _a为常数。将重组数据 ^Z(P)进行预旋转处理，旋转因子为 a'^ , P = 0，1，2，—， /2_1。 , _a为常数。

由于旋转因子中的项可以写成以下形式：

Wl = cos¾tll—j sin^l), p = 0,1,2, ···, - 1 满足 _COS^¾^ = _COS^¾^ = _Sin^½ = 0,l,2,...,¾- 1 和 s^m

= 0,1,2,… - 1的条件，因此具体的实现过程中只需要存储一个％点的余弦数据表格 a' ^' ^^⁰，¹，²…^— ¹或者正 — — 弦数据表格 ^a · ^， = 0，1，2…， ½ _ 1即可。

S103: 对预旋转处理后的数据进行 /²个点的快速傅里叶变换（Fast Fourier Transform, FFT )。

S104: 使用一个对称的旋转因子对 FFT变换后的数据进行后旋转处理，旋转因子为 b ₄r , q = Q L/2-\, b为常数；对 FFT变换后的数据进行后旋转处理，旋转因子为 b'^ , ? = 0，1，2，〜，J/2_1 即 _q 为 0 到 L/2-1 间的整数。其中 r4L ~^C0S" 4L - sm^^ = 0,l,2,.,¾-l_{? b}为常数。由于旋因子中的项同样可以写成以下形式：

W

q = 0,1,2,…， - 1 因此具体的实现过程中只需要存储一个％点的余弦数据表格 b-oos^, = 0Λ2,···^- i或者正弦数据表格 b'sm^S^ = 0,1,2,···, ½_1即可。在正变换中两个旋转因子的常数 a、 b乘积与逆变换中两个旋转因子的常数 c、 d乘积等于 ²/ 时，以正变换的输出数据直接作为逆变换的输入数据，逆变换的结果可以完成完全重构，即恢复获得进行正变换前的数据，通常在实际操作中并不一定需要实现完全重构。关于常数 a和 b取值的选择，为了实现完全重构，只要满足正变换 a与 b的乘积乘上反变换 c与 d的乘积为 ²/ 即可。在一个实

施例中， a和 b的乘积等于

，这样预旋转与后旋

— cos^2r(2p+1), π = 012 ··· ^L/-\

转只需要存储一个^点的余弦数据表格 ^4L " ' ' ² 即可。

S105: 获得频域数据。

后旋转之后的数据的实部记为 ²?), 是频域数据的奇数频率；后旋转之后的数据的虚部的相反数记为 (^-1_²），是频域数据的偶数频率₍

2-1

为后旋转之后的数据 y(k) =

由于原始 DCT-IV 变换的公式

Α· = 0,1,2,…， -1 ：

其中

L/2-1

W_L p,q = Q,...,LI2— \

z(p) = x(2 p) + jx(L

可以重新写为：

其中，为了简化计算，可以取^ 近似为 1, 近似为 1,

2-l

= x(2p) + jx(L-l- 2p), p = 0,1,2,… , J / 2— 1

： 2g(2pfl)

cos. _ _/· sin^^ , p = 0,1,2,…^― 1

旋转因子具备对称性特征，即满足

2,(2(_½：₁-_¾)+1) = _cos2,(L (_2q+1J) = ^ = 0,1,2,···,¾-1 s.^ 2,(2(_½-₁-_¾)₊₁) = _s^ 2,₍L-₄(2_q+1)) = cosi^iii , _q = o,i,2,… — 1

同样旋转因子也满足这样的对称性。 - - 本实施例中，由于预旋转和后旋转步骤使用的旋转因子具有对称性，实现的时候，只需要为存储一个 /²个点的余弦表或一个 /²个点的正弦表，从而可以降低数据存储量。同时，使用 FFT可以加快 DCT-IV变换的速度，降低其计算复杂度。进一步地，省去固定旋转步骤可以在变换满足重构特性情况下进一步降低计算复杂度。

在另一个实施例中，在获得频域数据之前还包括：利用固定旋转补偿因子进行固定旋转补偿的步骤。由于变换公式中固定旋转补偿与其他部分的操作是相乘的关系，根据乘法的交换特性，可以进行一次固定旋转补偿也可以进行多次固定旋转补偿，固定旋转补偿的执行顺序可以在获得频域数据步骤之前的任意顺序。

在一个实施例中，釆用 ^ϊ³进行一次固定旋转补偿；该固定旋转补偿的步骤可以是在预旋转之前或者在预旋转之后进行，还可以在后旋转之前进行或者在后旋转之后进行。执行一次补偿时，固定旋转补偿因子可以是 ^ϊ³。为了进一步降低计算复杂度，可以用一些近似值来代替的固定旋转补偿。

= e— Ί = ^，因此可以泰勒级数展开进行近似，例如用一阶泰勒级数展开结果作为 ^ϊ³的近似值，其中， ^{1 +}· '( )表示实部为 1虚部为 fe)的复数。

在另一个实施例中，执行两次固定旋转补偿时，固定旋转补偿因子可以是 ^ϊ¹和 ^」，也可以是它们的近似值。第一次固定旋转补偿的补偿因子是 ^―和 ^w^中的任何一个，第二次固定旋转补偿的补偿因子是和 ^―」中的另一个。第一次固定旋转补偿可以是在预旋转之前或者在预旋转之后进行，第二次固定旋转补偿可以在后旋转之前进行或者在后旋转之后进行。为了进一步降低计算复杂度可以用一些近似值如泰勒级数展开来代替或进行固定旋转补偿。例如用一阶泰勒级数展开结果 ^{1 +} 作为 ^―的近似值，用一阶泰勒级数展开结果作为」的近似值。

显然的，由于乘法的交换特性，可以进行三次及以上的固定旋转补偿，各补偿因子的乘积为 ^ϊ³或者至少一个补偿因子为乘积为的至少一个因子的近似值。固定旋转补偿的因子也可以为 1。

本实施例中，由于增加了固定旋转补偿的步骤，可以保证变换后得到数据与原始 DCT-IV型变换后得到的数据一致，提高 DCT-IV变换的准确度。 - - 另外， DCT-IV的反变换与正变换的步骤是基本一致的，不同的是，逆变换中首先是对频域数据进行重组，而最终变换以后得到的是时域数据，旋转因子中的常数 a和 b变为常数 c和 d。

参考图 2, 本发明一个实施例提供的一种信号处理方法，用于实现编码过程中的时域到频域的 DCT-IV变换，减少变换过程的存储量，该方法包括如下步骤：

S201: 对时域数据进行重组，得到重组后的数据。

z(p) = x(2p)+j-x(L-\-2p\p = 0X --;^L/₂-\

S202: 使用一个对称的旋转因子对重组后的数据进行预旋转，旋转因子为 - = 0,l,2,…， /2— 1 , _a为常数。将重组数据 ^Z(P)进行预旋转处理，旋转因子为 a'^ , P = 0，1，2，—， /2_1。

a为常数。

由于旋转因子中的项可以写成以下形式：

Wl = cos¾tll—j sin^l), p = 0,1,2, ···, - 1 满足

和 s'^¾^ = s''"^ ^ = = 0，1，2，···，¾ _ 1的条件，因此具体的实现过程中只需要存储一个％点的余弦数据表格 ^a'^C°^sl£^， :⁰'¹'²-^-¹或者正弦数据表格 ^a'^sinl£^' = 0'1'2·"，¾_1即可。

S203: 进行第一次固定旋转补偿。

对预旋转后的数据进行固定旋转补偿，固定旋转补偿因子为。为了进一步降低计算复杂度可以用一些近似值如泰勒级数展开来代替进行固定旋转补偿。例如用一阶泰勒级数展开结果 ^{1 +} 作为的近似值进行固定旋转 — — 补偿。

S204: 对固定旋转补偿后的数据进行个点的 FFT变换。

S205: 进行第二次固定旋转补偿。

将 FFT变换后的数据乘以进行固定旋转补偿，或者将 FFT变换后的数据乘以的近似值进行固定旋转补偿，近似值可以使用的泰勒级数展开获得，如使用一阶泰勒级数展开结果 "⁷'^)作为的近似值。

S206:对使用一个对称的旋转因子对固定旋转补偿后的数据进行后旋转处理，旋转因子为 b^^¹ , qH .，L/2-l 。其中 Wl

₌ ο,ΐ, 2, · · ·, ¾ - 1 _{? b}为常数。由于旋转因子中的 ^ ¹项同样可以写成以下形式：

Wl = cos^¾±ll—j sin^¾±ll, q = 0,1,2, ···, _ 1 因此具体的实现过程中只需要存储一个％点的余弦数据表格

0,1,2,··· - 1

即可。

关于常数 a和 b取值的选择，为了实现完全重构，只要满足正变换两个旋转因子的常数 a与 b的乘积乘上反变换两个旋转因子的常数 c与 d的乘积为 ²/ 即可。一个实施例中， a和 b的乘积等于 ,c和 d的乘积也等于。另一个实施例中 a和 b的乘积等于 l,c和 d的乘积等于 ²/ 。另一个实施例中选取

这样预旋转与后旋转只需要存储一个％点的余弦数据表格

即可。

S207: 获得变换后的频域数据。

后旋转之后的数据的实部记为 ²?), 是频域数据的奇数频率；后旋转之 - - 后的数据的虚部的相反数记为 -ι-¾ , 是频域数据的偶数频率。

y(2q) = Re{Z(q)}

r— Λ -> Q- 0，1，2，…， /2— 1

y(L-\-2q) = -lm{Z(q)} 0".., /2— 1

，或者 p,q = 0,...,L/2-l

，为后旋转之后的数据。值得注意的是，实施例中的使用进行固定旋转补偿的步骤不仅可以在预旋转之后执行，也可以在预旋转之前进行；使用进行固定旋转补偿的步骤不仅可以在后旋转之前执行，也可以在后旋转之后进行。另外变换公式中两个固定旋转补偿与其他部分的操作是相乘的关系，可以适用乘法的交换特性，因此使用进行固定旋转补偿的步骤可以与使用进行固定旋转补偿的步骤互换。本实施例中，釆用两次固定旋转补偿步骤，可以保证 FFT变换的输入数据与原 DCT-IV变换中 FFT的输入数据一致，同时可以保证变换后得到数据与原始 DCT-IV型变换后得到的数据一致，提高 DCT-IV变换的准确度。另外， DCT-IV的反变换与正变换的步骤是基本一致的，不同的是，逆变换中首先是对频域数据进行重组，而最终变换以后得到的是时域数据，旋转因子中的常数 a和 b变为常数 c和 d 参考图 3, 本发明一个实施例提供的一种信号处理方法，用于实现编码过程中的时域到频域的 DCT-IV变换，减少变换过程的存储量，该方法包括如下步骤：

S301: 对时域数据进行重组，得到重组后的数据。 — — 假设为需要进行 DCT-IV型变换的数据，该数据可以是经过加窗等预处理步骤后的数据。对数据进行重组，得到重组后的数据

ζθ) = ¾：2/7) + j' (Z _ 1— 2/?) , = 0, 1,2,…^― 1

S302: 使用一个对称的旋转因子对重组后的数据进行预旋转，旋转因子为 - ^W = 0, …， /2_1, _a为常数。将重组数据 ^Z(W进行预旋转处理，旋转因子为 a'^ , P = 0 D /2-l。其中 ffT=_C。^-/_Si_n^ ^ = 0，l，2，'''， -l， _a为常数。由于旋转因子中的项可以写成以下形式：

Wl = cos¾tll—j sin^l), p = 0,1,2, ·■·, - 1 满足 cos ^ _^cos ^ _^sm ,/₄ i 和 £ =。Ο5^ ? = 0,1,2 ··,¾— 1的条件,因此具体的实现过程中只需要存储一个点的余弦数据表格 ^a'^COs ^' ^'¹'²… 或者正弦数据表格 a'sin^H ¹'²"'^— ¹即可。

S303: 对预旋转处理后的数据进行个点的快速傅里叶变换（Fast Fourier Transform, FFT )。

S304: 进行固定旋转补偿。将 FFT变换后的数据乘以进行固定旋转补偿，或者将 FFT变换后的数据乘以的近似值进行固定旋转补偿，近似值可以使用的泰勒级数展开获得，如使用一阶泰勒级数展开结果 ^{1 +} 作为的近似值。 S305: 使用一个对称的旋转因子对固定旋转补偿后的数据进行后旋转处理，旋转因子为 b¾]⁺¹ , ? = 0,1,2,…， /2-1 。其中 WT = cos^l -_jsn^q = 0,1,2, ···, - 1 , _b为常数。 — ~ 由于旋转因子中的项同样可以写成以下形式：

0,1,2,··· - 1

即可。关于常数 a和 b取值的选择，为了实现完全重构，只要满足正变换两个旋转因子的常数 a与 b的乘积乘上反变换两个旋转因子的常数 c与 d的乘积为 ²/ 即可。一个实施例中， a和 b的乘积等于

, 这样预旋转与后旋转只需要存储一个％点的余弦数据表格

即可。

S306: 获得频域数据。后旋转之后的数据的实部记为 ²?), 是频域数据奇数频率；后旋转之后的数据的虚部的相反数记为 -ι-¾ ，是频域数据偶数频率。

或

「，

= ∑ ( + 7¾) * PY C' p,q = 0,...,L/2-l

或。 L 」，为后旋转之后的数据' - - 实施例中的进行固定旋转补偿的步骤不仅可以在后旋转之前执行，也可以在预旋转之前或者 FFT变换之前或者后旋转之后进行，因为变换公式中固定旋转补偿与其他部分的操作是相乘的关系，可以适用乘法的交换特性。

本实施例中，由于增加了固定旋转补偿的步骤，可以保证变换后得到数据与原始 DCT-IV型变换后得到的数据一致，提高 DCT-IV变换的准确度。

DCT-IV的反变换与正变换的步骤是类似的，不同的是反变换的输入为频和时域数据时，即构成 DCT-IV的反变换的实施例。另外， DCT-IV的反变换与正变换在进行固定旋转补偿的顺序上可以不同。例如，可以正变换时在后旋转之后进行固定旋转补偿，而逆变换时在预旋转之前进行固定旋转补偿。语音 /音频编码中， MDCT变换由于利用了时域混叠消除（ TDAC )技术减

"边界效应 " 也得到了广泛的应用。 MDCT变换的公式为：

2π

A = 0 1… i— l

N " ，其中， A是归一化因子， A 为常数。可见，直接根据变换公式进行 MDCT正变换与逆变换会导致很高的计算复杂度和存储量，特别是对于点数较大的 MDCT变换。由于 MDCT变换在实时通信领域特别是音频编码中的广泛应用，提供一种快速的 MDCT变换方法也同样成为一种迫切的需要。

参考图 4, 本发明一个实施例提供的一种信号处理方法，用于实现编码过程中的时域到频域的 MDCT变换，减少变换过程的存储量，该方法包括如下步骤：

S401 : 对时域数据进行预处理，得到预处理后的数据。 - 假设 ^ "为需要进行 MDCT变换的数据，该数据可以是经过加窗等步骤处理后的数据。对数据 ^ "进行重组，得到重组后的数据 ^M "：

= 0₂„ + jz_N/i― ₂ n = 0,1,2, ···, _ 1

或者表示成：

= ~y ¾— « = 0,1,2,···, ¾-l

= y_n-% _ y « = ¾,¾ + !,···,¾-!

S402: 使用一

为 a-C" = 0,1,2, ··· /— 1。将重组数据 ^M "进行预旋转处理，旋转因子为 a'^ « = 0,1,2,···, ¾-l JL 中 ^{=e =cos(}^^⁾" ^sm(^^⁾ , a为常数。

由于旋转因子中的 C项可以写成以下形式： y ^ ⁼

满足 _COSe (" + 0.5)) = _Sin ( (w/— 1— " + 0.5)) 和 sin(^(« + 0.5)) = ( (Υ + 0.5))的条件，因此具体的实现过程中只需要存储一个 ^Ν 点的余弦数据表格 ^a'^c°s( (" ^{+ 5}))， "^Ί'²···^-¹或者正弦数据表格 a-sin( (w + 0.5))， « = 0,1,2,· -·,¾-1即可

S403: 进行第一次固定旋转补偿。对预旋转后的数据进行固定旋转补偿，固定旋转补偿因子为 ^^—° ³⁷⁵。为了进一步降低计算复杂度可以用一些近似值如泰勒级数展开来代替^^—° ³⁷⁵进行固定旋转补偿。例如用一阶泰勒级数展开结果 ^{1 +} (^作为的近似值进行固定旋转补偿。 S404: 对固定旋转补偿后的数据进行 Ν/⁴个点的 FFT变换。

S405: 进行第二次固定旋转补偿。

将 FFT变换后的数据乘以 ^{W 5}进行固定旋转补偿，或者将 FFT变换后的数据乘以 ^{W 5}的近似值进行固定旋转补偿，近似值可以使用 ^{W 5}的泰勒级数展开获得，如使用一阶泰勒级数展开结果 ^{1 +} 作为 ^{W 5}的近似值。

S406: 釆用一个对称的旋转因子对固定旋转补偿后的数据进行后旋转处理，旋转因子为 ^{b 0}.⁵ , k =

。其中

^{Wn = e = cos(}^^⁾" ^sm(^^⁾ , b为常数。

由于旋转因子中的^ 项同样可以写成以下形式： yy

因此具体的实现过程中只需要存储一个 N/⁴点的余弦数据表格 a-cos(^(^ ₊ 0.5) k = 0X2,..,%-\ 或者正弦数据表格 a · sin( + 0.5))， A: = 0,1,2, ···,%- 1即可。

S407: 获得频域数据。

后旋转之后的数据的实部记为是频域数据的奇数频率；后旋转之后的数据的虚部的相反数记为 ^λ-^，是频域数据的偶数频率。

频域数据也即最终的谱为 ^ ' ^k = °^²-^-ί , 可表述为：

X_2k = = 0,1,2, ···,¾-! u_n · W-^{0 75} · a k = 0,1,2, ···,¾-!

Έ^Μ»·(^{1 +} · ^a · n 0,1,2,… , — 1

- -

X%-,-_k = - Im (! + 》 · ^b ·

· + · ^a ⁵ ζ t = ο,ι,2,···, ¾-i

值得注意的是，实施例中的使用进行第一次固定旋转补偿的步骤不仅可以在预旋转之后执行，也可以在预旋转之前进行；使用 ° ³⁷⁵进行第二次固定旋转补偿的步骤不仅可以在后旋转之前执行，也可以在后旋转之后进行。

由于变换公式中固定旋转补偿与其他部分的操作是相乘的关系，可以适用乘法的交换特性，可以进行一次固定旋转补偿也可以进行多次固定旋转补偿，固定旋转补偿的执行顺序可以在获得频域数据步骤之前的任意顺序。各补偿因子的乘积为 ^W 或者乘积为 ^W 的至少一个因子的近似值。

本实施例中，釆用具有对称性的旋转因子可以减少存储量，改进前的方法的存储量为 N/²个点，改进后的方法存储量为 NM个点。同时增加固定旋转补偿的步骤，提高 MDCT变换的准确度，保证变换后得到的数据与原始 MDCT变换后得到的数据一致。原始的基于 N/4点 FFT的 MDCT快速变换的公式为：

可以重新写为：

其中，

W;⁺⁰-⁵ =e'^{j (}"^+ft5^co ^{2g( +ft5}^- _Sm(^2g("_N ^+ft5)) = 0,1,2···,%— 1

很容易证明这个改写后的旋转因子具备对称性特征，即 C满足： — ~ cos(^ (« + 0.5)) = sin(^ (¾— 1— w + 0.5)) .

sin( (w + 0.5)) = cos(^ ¾— 1— w + 0.5)). 同样^ ⁵也满足这样的对称性。

利用这个对称性特征可以减少存储量。实现的时候，只需要为^" 存储一个 /⁴个点的余弦表或一个 /⁴个点的正弦表，并在做 /⁴个点的 FFT变换之前进行^^— ° ³⁷⁵的固定旋转补偿，在做 N"个点的 FFT变换之后进行 ° ³⁷⁵的固定旋转补偿。可以证明这个变换也是完全重构的。

为了进一步降低计算复杂度可以用一些近似值如泰勒级数展开来代替

^{W 5}进行固定旋转补偿。例如用一阶泰勒级数展开结果 ^{1 +} (^作为 ^{W 5}的近似值。

参考图 5, 本发明一个实施例提供的一种信号处理方法，用于实现编码过程中的频域到时域的 MDCT变换，减少变换过程的存储量，该方法包括如下步骤：

S501: 对频域数据进行重组，得到重组数据。

数据重组后获取的中间变量为（ ^{m +jX}^-^ , k = 0,1,2,···,%— 1。

S502: 使用一个对称的旋转因子对重组数据进行预旋转处理。

将重组数据（ ^+^ ^)进行预旋转处理，旋转因子为 ^{c +}。.⁵ , k = 0,l,2,---,¾-l_o 其中^ ^^)— (^^) , _c为常数。

S503: 进行第一次固定旋转补偿。

对预旋转后的数据进行固定旋转补偿，固定旋转补偿因子为 ^^—° ³⁷⁵。为了进一步降低计算复杂度可以用一些近似值如泰勒级数展开来代替^^—° ³⁷⁵进行固定旋转补偿。例如用一阶泰勒级数展开结果 ^{1 +} (^作为 ^{w 5}的近似值进行固 ~ ~ 定旋转补偿。

S504: 对固定旋转补偿后的数据进行 N/⁴个点的 FFT变换。 S505: 进行第二次固定旋转补偿。

将 FFT变换后的数据乘以 ^{W S}进行固定旋转补偿，或者将 FFT变换后的数据乘以 ^{W 5}的近似值进行固定旋转补偿，近似值可以使用 ^{W 5}的泰勒级数展开获得，如使用一阶泰勒级数展开结果 ¹ 作为 ^{W 5}的近似值。

S506: 使用一个对称的旋转因子对固定旋转后的数据进行后旋转处理。对固定旋转补偿后的数据进行后旋转处理，旋转因子为 ^d'c，！

// = 0,1,2,···,¾-1_ο 其中^ ， _d为常数。 S507: 获得时域数据。

得到时域数据 '" = GH — ¹。

¾-ι = ¾₊₂„ = Im{w_n } « = 0,1,2,···, ¾-l

其中

，值得注意的是，关于常数、 b、 c和 d取值的选择，为了实现完全重构，只要满足正变换 a与 b的乘积乘上反变换 c与 d的乘积为 ⁴/ 即可。本实施例中选取 _a=h=c=d=^l

a" ~^C~ , 这只需要存储一个 Ν/⁴ 点的余弦数据表格

即可。

同样的，实施例中的使用 ^^—° ³⁷⁵进行第一次固定旋转补偿的步骤不仅可以在预旋转之后执行，也可以在预旋转之前进行；使用^— ^ΰ37ί进行第二次固定旋 - . 转补偿的步骤不仅可以在后旋转之前执行，也可以在后旋转之后进行。由于乘法的交换特性，也可以进行三次及以上的固定旋转补偿，各补偿因子的乘积为或者乘积为 ^{W 75}的至少一个因子的近似值。

本实施例中，釆用两次固定旋转补偿的步骤，提高 MDCT变换准确度，保证变换后得到的数据与原始 MDCT变换后得到的数据一致。参考图 6, 本发明一个实施例提供的一种信号处理方法，用于实现编码过程中的时域到频域的 MDCT变换，减少变换过程的存储量，该方法包括如下步

S601: 对时域数据进行预处理，得到预处理后的数据 _t 假设 ^ "为需要进行 MDCT变换的数据，该数据可以是经过加窗等步骤处理后的数据。对数据 ^ "进行重组，得到重组后的数据 ^M "：

= 0₂„ + ¾-ι-2„) n = 0,1,2, ···, _ 1 其中：

或者表示成：

S602: 使用一

为 a ^+U5， " = 0,1,2,···,/— 1。将重组数据 ^M "进行预旋转处理，旋转因子为 ^¾'^^+⁵, " = 0，1，²，···，％-1。其中 ^{=e =cos(}^^⁾" ^sm(^^⁾ , a为常数。

由于旋转因子中的项可以写成以下形式: - -

满足 _COSe (" + 0.5)) = _Sin ( (w/— 1— " + 0.5)) 和

0.5))的条件，因此具体的实现过程中只需要存储一个 ^Ν 点的余弦数据表格 ^a'^c°s( (" ^{+ 5}))， "^Ί'²···^-¹或者正弦数据表格 a-sin( (w + 0.5))， « = 0,1,2,· -·,¾-1即可

S603: 进行固定旋转补偿。

对预旋转后的数据进行固定旋转补偿，固定旋转补偿因子为 ^—°⁷⁵。为了进一步降低计算复杂度可以用一些近似值如泰勒级数展开来代替^^—°⁷⁵进行固定旋转补偿。例如用一阶泰勒级数展开结果 ^{1 +} 作为 ^W 的近似值进行固定旋转补偿。

S604: 对固定旋转补偿后的数据进行 N/⁴个点的 FFT变换。

S605：釆用一个对称的旋转因子对 FFT变换后的数据进行后旋转处理。对 FFT变换后的数据进行后旋转处理，旋转因子为

Λ = 0,1,2,···,¾-1_ο 其中^ ^^)— (^^) , _b为常数。

由于旋转因子中的项同样可以写成以下形式：因此具体的实现过程中只需要存储一个 N/⁴点的余弦数据表格 a-cos(^(^ ₊ 0.5) h = 0X2,..,%-\ 或者正弦数据表格 a · sin( + 0.5))， A: = 0,1,2, ···,%- 1即可。

S606: 获得频域数据。

后旋转之后的数据的实部记为是频域数据的奇数频率；后旋转之后的数据的虚部的相反数记为 ^Xn - , 是频域数据的偶数频率。 — ~ 频域数据也即最终的谱为 ^ ' ^k = °^²-^-¹ , 可表述为： 0,… ,/— l 0,… , — l

或者

0,… ,/— l k = 0,-,y₄-l

值得注意的是，实施例中的进行固定旋转补偿的步骤可以在预旋转之后或者预旋转之前进行，也可以在后旋转之前或者后旋转之后进行，因为变换公式中固定旋转补偿与其他部分的操作是相乘的关系，可以适用乘法的交换特性。参考图 7, 本发明一个实施例提供的一种信号处理方法，用于实现编码过程中的频域到时域的 MDCT变换，减少变换过程的存储量，该方法包括如下步骤：

S701: 对频域数据进行重组，得到重组数据。

数据重组后获取的中间变量为 ^{∑lk +jX}"A-^ , k = 0,l,2,-,¾-l_o

S702: 使用一个对称的旋转因子对重组数据进行预旋转处理。

将重组数据（ ^ +^ ^)进行预旋转处理，旋转因子为 ^{c +}。.⁵ , k = 0,l,2,---,¾-l_o 其中^ ^^)— (^^) , _c为常数。

S703: 进行固定旋转补偿。

对预旋转后的数据进行固定旋转补偿，固定旋转补偿因子为 ^—°⁷⁵。为了 - - 进一步降低计算复杂度可以用一些近似值如泰勒级数展开来代替^^—°⁷⁵进行固定旋转补偿。例如用一阶泰勒级数展开结果 ^{1 +} 作为 W的近似值进行固定旋转补偿。

S704: 对固定旋转补偿后的数据进行 N/⁴个点的 FFT变换。

S705: 釆用一个对称性的旋转因子对 FFT变换后的数据进行后旋转处理。对 FFT变换后的数据进行后旋转处理，旋转因子为 ^d'W，

+0.5 2π +0.5) 2π +0.5)

« = 0,1,2,···,¾-1_ο 其中^ — (^^) , _d为常数。

S706: 获得时域数据。得到时域数据^'"^'¹'²'…'^—¹。

¾-ι-2„ = ¾₊2„ = ^lm } " = Ο,Ι,², — 1

= ¾_¾-i-2„ = ^lm } " = 1, ···,%- 1

_ ^/8 « = 0,1,2,· -·,¾-1

其中 ί^Μ"+% = « ,

值得注意的是，关于常数&、 b、 c和 d取值的选择，为了实现完全重构，只要满足正变换 a与 b的乘积乘上反变换 c与 d的乘积为 ⁴/ 即可。本实施例中选取

^{¾ C} ^ , 这样只需要存储一个 ΝΜ 点的余弦数据表格

^=cos(^(A + 0.5)), Λ = 0,1,2,···,¾-1

yJN 即可。同样的，实施例中的进行固定旋转补偿的步骤可以在预旋转之后或者预旋转之前进行，也可以在后旋转之前或者后旋转之后进行，因为变换公式中固定旋转补偿与其他部分的操作是相乘的关系，可以适用乘法的交换特性。 - - 本领域普通技术人员可以理解实现上述实施例方法中的全部或部分流程，是可以通过计算机程序来指令相关的硬件来完成，所述的程序可存储于一计算机可读取存储介质中，该程序在执行时，可包括如上述各方法的实施例的流程。其中，所述的存储介质可为磁碟、光盘、只读存储记忆体（Read-Only Memory, ROM)或随机存储记忆体（Random Access Memory, RAM)等。参考图 8与上述方法实施例相关联，本发明信号处理装置的一个实施例包括：

重组单元 801, 用于对输入数据进行重组，得到重组后的数据；

预旋转单元 802, 用于使用一个对称的旋转因子对重组后的数据进行预旋转，旋转因子为 a ^¹, Ρ = υ - , _a为常数；

Fourier Transform, FFT ) 变换， L为输入数据的长度；

后旋转单元 804, 用于使用一个对称的旋转因子对 FFT变换后的数据进行后旋转处理，旋转因子为 b'^^ , ? = 0，…， /2-1, b为常数；

输出单元 805, 用于获得输出数据。

上述信号处理装置可以用于实现编解码过程中的时域到频域或频域到时域的 DCT-IV变换，正变换时输入数据为时域数据，输出数据为频域数据；逆变换时输入数据为频域数据，输出数据为时域数据。

在另一个实施例中，信号处理装置还包括：

固定旋转补偿单元，用于利用固定旋转补偿因子进行固定旋转补偿处理。一个实施例中，固定旋转补偿单元用于进行至少一次固定旋转补偿，至少 - - 一次固定旋转补偿的旋转补偿因子的乘积为^^。

另一个实施例中，固定旋转补偿单元用于进行至少一次固定旋转补偿，所述至少一次固定旋转补偿的旋转补偿因子为乘积的至少一个因子的一阶泰勒级数展开值。为了满足重构， a和 b的乘积可以等于

。本发明提供的时域到频域的信号处理装置的一个实施例，用于实现编码过程中的时域到频域的 MDCT变换，减少变换过程的存储量，参考图 9, 该信号处理装置包括：

预处理单元 901 , 用于对时域数据进行预处理，得到预处理后的数据；预旋转单元 902，用于使用旋转因子 ^a · C对所述预处理后的时域数据进行预旋转处理；

变换单元 903 ,用于对所述预旋转处理后的数据进行 N/4点的快速傅里叶变换；

后旋转单元 904, 用于使用旋转因子 b ' ^^+o⁵对所述离散傅里叶变换后的数据进行后旋转处理，获得频域数据；其中，还包括：

固定补偿单元 905 , 用于利用固定旋转补偿因子进行固定旋转补偿处理；所述 a、 b为常数，所述 N为时域数据长度，所述 ^{= e N}。

一个实施例中，固定旋转补偿单元用于进行至少一次固定旋转补偿，该至少一次固定旋转补偿的旋转补偿因子的乘积为 W。

另一个实施例中，固定旋转补偿单元用于进行至少一次固定旋转补偿，该 - - 至少一次固定旋转补偿的旋转补偿因子为乘积为的至少一个因子的一阶泰勒级数展开值。

本发明提供的频域到时域的信号处理装置的一个实施例，用于实现编码过程中的频域到时域的 MDCT变换，减少变换过程的存储量，参考图 10, 该信号处理装置包括：

重组单元 1001 , 用于对频域数据进行重组，得到重组后的数据；预旋转单元 1002，用于使用旋转因子 ^e · ^W^⁺⁰'⁵对所述重组后的频域数据进行预旋转处理；

变换单元 1003 ,对所述预旋转处理后的数据进行 N/4点的快速傅里叶变换；后旋转单元 1004，使用旋转因子 ^d · C对所述快速傅里叶变换后的数据进行后旋转处理；其中，还包括：

固定补偿单元 1005, 用于利用固定旋转补偿因子进行固定旋转补偿处理；所述 c、 d为常数，所述 N为两倍的频域数据长度，所述 ^{= e N}。

另一个实施例中，固定旋转补偿单元用于进行至少一次固定旋转补偿，该至少一次固定旋转补偿的旋转补偿因子为乘积为 W的至少一个因子的一阶泰勒级数展开值。与本文所揭示实施例关联描述的中例示性逻辑块、模块及电路均可由下列装置构建或实施:通用处理器、数字信号处理器 (DSP)、专用集成电路 (ASIC)、场可编程门阵列 (FPGA)或其它可编程逻辑装置、离散门或晶体管逻辑、离散 _ _ 硬件组件、或其设计用于实施上文所述功能的任何组合。通用处理器可为微处理器，但另一选择为，处理器也可为任何常规处理器、控制器、微控制器或状态机。处理器也可构建为计算装置的组合，例如， DSP与微处理器的组合、多个微处理器的组合、一个或一个以上微处理器与 DSP核心的联合，或任一其它此类配置。以上所述仅为本发明的几个实施例，本领域的技术人员依据申请文件公开

+

Claims

权利要求

1、一种数据处理方法，其特征在于，包括：

对输入数据进行重组，得到重组后的数据；

使用一个对称的旋转因子对重组后的数据进行预旋转，旋转因子为 ^a-^ , P = , ..., L/2 - l ^ _a为常数；

对予贞旋转处理后的数据进行 2个点的'！ "夬速傅里 ( Fast Fourier Transform, FFT ) 变换， L为输入数据的长度；

使用一个对称的旋转因子对 FFT 变换后的数据进行后旋转处理，旋转因子为 b ₄r , ? = 0，···， 2_1 , _b为常数；

获得输出数据。

2、根据权利要求 1所述的方法，其特征在于，在获得输出数据之前还包括：利用固定旋转补偿因子进行固定旋转补偿的步骤。

3、根据权利要求 2所述的方法，其特征在于，利用固定旋转补偿因子进行固定旋转补偿的步骤包括：

进行至少一次固定旋转补偿，所述至少一次固定旋转补偿的旋转补偿因子的乘积为 ^ϊ³。

4、根据权利要求 2所述的方法，其特征在于，利用固定旋转补偿因子进行固定旋转补偿的步骤包括：

进行至少一次固定旋转补偿，所述至少一次固定旋转补偿的旋转补偿因子为乘积为的至少一个因子的一阶泰勒级数展开值。

5、根据权利要求 1所述的方法，其特征在于，所述 a和 b的乘积等于 I。

6、一种时域到频域的信号处理方法，其特征在于，包括：

对时域数据进行预处理，得到预处理后数据；

对所述预旋转处理后的数据进行 N/4点的快速傅里叶变换；

其中，在获得频域数据之前还包括：利用固定旋转补偿因子进行固定旋转补偿的步骤；所述 a、 b为常数，所述 N为时域数据长度，所述 ^{= e} 。

7、根据权利要求 6所述的方法，其特征在于，利用固定旋转补偿因子进行固定旋转补偿的步骤包括：

进行至少一次固定旋转补偿，所述至少一次固定旋转补偿的旋转补偿因子的乘积为^ 。

8、根据权利要求 6所述的方法，其特征在于，利用固定旋转补偿因子进行固定旋转补偿的步骤包括：

进行至少一次固定旋转补偿，所述至少一次固定旋转补偿的旋转补偿因子为乘积为 ^{w 75}的至少一个因子的一阶泰勒级数展开值。

2 9、根据权利要求 1所述的方法，其特征在于，所述 a和 b的乘积等于。

10、一种频域到时域的信号处理方法，其特征在于，包括：

对频域数据进行进行重组，得到重组后的数据；

使用旋转因子 ^e ' ^^k+°'⁵对所述重组后的数据进行预旋转处理；

对所述预旋转处理后的数据进行 N/4点的快速傅里叶变换；

使用旋转因子 ^d · ^W 。.⁵对所述快速傅里叶变换后的数据进行后旋转处理；对所述后旋转处理后的数据进行后处理，获得时域数据；

其中，在获得时域数据之前还包括：利用固定旋转补偿因子进行固定旋转补偿的步骤；所述 c、 d为常数，所述 N为两倍的频域数据长度，所述 = ^e ' ^Ν。

11、根据权利要求 10所述的方法，其特征在于，利用固定旋转补偿因子进行固定旋转补偿的步骤包括：

12、根据权利要求 10所述的方法，其特征在于，利用固定旋转补偿因子进行固定旋转补偿的步骤包括：

13、一种信号处理装置，其特征在于，包括：重组单元，用于对输入数据进行重组，得到重组后的数据；

预旋转单元，用于使用一个对称的旋转因子对重组后的数据进行预旋转，旋转因子为 a ^¹ , Ρ = 0，· · ·，^/2 _1 , _a为常数；

变换单元，用于对预旋转处理后的数据进行个点的快速傅里（Fast Fourier Transform, FFT ) 变换， L为输入数据的长度；

后旋转单元，用于使用一个对称的旋转因子对 FFT 变换后的数据进行后旋转处理，旋转因子为 b '^^ , ? = 0，…， /2 -1 , b为常数；

输出单元，用于获得输出数据。

14、根据权利要求 13所述的装置，其特征在于，还包括：

固定旋转补偿单元，用于利用固定旋转补偿因子进行固定旋转补偿处理。

15、根据权利要求 14所述的装置，其特征在于，所述固定旋转补偿单元用于进行至少一次固定旋转补偿，所述至少一次固定旋转补偿的旋转补偿因子的乘积为 ^ϊ³。

16、根据权利要求 14所述的装置，其特征在于，所述固定旋转补偿单元用于进行至少一次固定旋转补偿，所述至少一次固定旋转补偿的旋转补偿因子为乘积为的至少一个因子的一阶泰勒级数展开值。

17、根据权利要求 14所述的装置，其特征在于，所述 a和 b的乘积等于。

18、一种时域到频域的信号处理装置，其特征在于，包括：

预处理单元，用于对时域数据进行预处理，得到预处理后的数据；预旋转单元，用于使用旋转因子对所述预处理后的数据进行预旋转处理；

变换单元，用于对所述预旋转处理后的数据进行 N/4 点的快速傅里叶变换；

后旋转单元 , 用于使用旋转因子 ^b ' 对所述离散傅里叶变换后的数据进行后旋转处理，获得频域数据；其中，还包括：

固定补偿单元，用于利用固定旋转补偿因子进行固定旋转补偿处理；所述 a、 b为常数，所述 N为时域数据长度，所述 = ^e—^{j r}。

19、根据权利要求 18所述的装置，其特征在于，所述固定旋转补偿单元用于进行至少一次固定旋转补偿，所述至少一次固定旋转补偿的旋转补偿因子的乘积为^ 。

20、根据权利要求 18所述的装置，其特征在于，所述固定旋转补偿单元用于进行至少一次固定旋转补偿，所述至少一次固定旋转补偿的旋转补偿因子为乘积为 ^{W 75}的至少一个因子的一阶泰勒级数展开值。

21、一种频域到时域的信号处理装置，其特征在于，包括：

重组单元，用于对频域数据进行重组，得到重组后的数据；

预旋转单元，用于使用旋转因子 ^e ' ^^k+°'⁵对所述重组后的数据进行预旋转处理；

变换单元，对所述预旋转处理后的数据进行 N/4点的快速傅里叶变换；后旋转单元，使用旋转因子 ^d · W对所述快速傅里叶变换后的数据进行后旋转处理；其中，还包括：

固定补偿单元，用于利用固定旋转补偿因子进行固定旋转补偿处理；所述 c、 d为常数，所述 N为两倍的频域数据长度，所述 ^{= e N}。

22、根据权利要求 21所述的装置，其特征在于，所述固定旋转补偿单元用于进行至少一次固定旋转补偿，所述至少一次固定旋转补偿的旋转补偿因子的乘积为^ 。

23、根据权利要求 21所述的装置，其特征在于，所述固定旋转补偿单元用于进行至少一次固定旋转补偿，所述至少一次固定旋转补偿的旋转补偿因子为乘积为 ^W 的至少一个因子的一阶泰勒级数展开值。