快速直接DCT实现
发布时间:2008/12/18 0:00:00 访问次数:724
dct的对称属性已经被byeong lee[141]用来构造类似fft的dct算法。由于其与radix-2 cooley-tukey fft的相似性,所以最终的算法称为快速dot或简称fct。换句话说,就是快速dct算法可以用矩阵结构开发[142]。由于dct是正交变换,所以可以通过转置逆dct(idct)得到dct。idct-ⅱ型有:
重复这一过程就可以进一步分解dct。给出的(6.62)与radix-2 fft旋转因子之间的比较表明,除法对fct似乎是必要的。所以旋转因子1/(2cn,kn)就应该预先被计算出来并储存在表中。这样的制表方法对于cooley-tukey fft也是适合的,因为在线计算二角函数一般是非常耗时间的。接下来用一个示例来说明fct。
例 8点 fct
对于8点fct,等式(6.60)至(6.65)式就变成:
这样,重构就变成:
等式(6.66)和(6.67)构成了图1中流程图的第一级,而(6.70)式和(6.71)式构成了流程图的最后一级。
图1 采用速记符号c[p]=1/(2cos(pπ16))的8点快速dot流程图
在图1中,输入序列x[k]是位逆序的。输出序列x[n]的顺序按下面的方式生成:由集合(0,1)开始通过增加—个前缀0和1形成新的集合。前缀是1时,前面格式中所有的位都是颠倒的。例如:从序列10得到两个子序列010和110=101。图2给出了这种帧格式的图解。
图2 8点快速dct的输入输出的置换
欢迎转载,信息来源维库电子市场网(www.dzsc.com)
dct的对称属性已经被byeong lee[141]用来构造类似fft的dct算法。由于其与radix-2 cooley-tukey fft的相似性,所以最终的算法称为快速dot或简称fct。换句话说,就是快速dct算法可以用矩阵结构开发[142]。由于dct是正交变换,所以可以通过转置逆dct(idct)得到dct。idct-ⅱ型有:
重复这一过程就可以进一步分解dct。给出的(6.62)与radix-2 fft旋转因子之间的比较表明,除法对fct似乎是必要的。所以旋转因子1/(2cn,kn)就应该预先被计算出来并储存在表中。这样的制表方法对于cooley-tukey fft也是适合的,因为在线计算二角函数一般是非常耗时间的。接下来用一个示例来说明fct。
例 8点 fct
对于8点fct,等式(6.60)至(6.65)式就变成:
这样,重构就变成:
等式(6.66)和(6.67)构成了图1中流程图的第一级,而(6.70)式和(6.71)式构成了流程图的最后一级。
图1 采用速记符号c[p]=1/(2cos(pπ16))的8点快速dot流程图
在图1中,输入序列x[k]是位逆序的。输出序列x[n]的顺序按下面的方式生成:由集合(0,1)开始通过增加—个前缀0和1形成新的集合。前缀是1时,前面格式中所有的位都是颠倒的。例如:从序列10得到两个子序列010和110=101。图2给出了这种帧格式的图解。
图2 8点快速dct的输入输出的置换
欢迎转载,信息来源维库电子市场网(www.dzsc.com)
上一篇:利用DFT计算DCT
热门点击
- D/A转换器的基本原理
- AD转换器的选择
- 语音信号的μ/A律压缩
- 并行A/D转换器AD574
- Bluestein Chirp-z变换
- 语音信号的采集和播放
- 语音信号模数/数模转换
- Cooley-Tukey FFT算法
- DFT和FFT算法的比较
- DFT的属性
推荐技术资料
- DS2202型示波器试用
- 说起数字示波器,普源算是国内的老牌子了,FQP8N60... [详细]