离散傅里叶变换 - 快速计算方法及C实现 - 第三篇

最新推荐文章于 2023-03-14 23:11:21 发布
VIP文章 最新推荐文章于 2023-03-14 23:11:21 发布
阅读量913 收藏 6
点赞数 1
分类专栏： 编程 文章标签： Radix-2 Radix-3 DFT FFT C
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/YU_Xianguo/article/details/99641213
版权

DFT – Fast algorithms and C implementations - Part3

Radix-3 DFT

我在第一篇博客里分析了Radix-2 DFT的计算原理，并据此写出了Radix-2 DFT的代码，并验证了其有效性。现在我们有信心认为，只要有了公式，就能写出正确的程序，产生正确的结果！

下面我们推导一下Radix-3 DFT的公式，其实非常简单：

设N长序列\lbrace x_j \rbrace_{j<N}的DFT为\lbrace X_k \rbrace_{k<N}，且N能被3整除，则：

\begin{align} X_k &=\sum_{j<N}{x_jW_N^{jk}} \nonumber \\ &=\sum_{j<N/3}{x_{3j}W_N^{3jk}}+\sum_{j<N/3}{x_{3j+1}W_N^{(3j+1)k}} +\sum_{j<N/3}{x_{3j+2}W_N^{(3j+2)k}}\nonumber \\ &= \sum{x_{3j}W_{N/3}^{jk}} + W_N^k \sum{x_{3j+1}W_{N/3}^{jk}} + W_N^{2k} \sum{x_{3j+2}W_{N/3}^{jk}} \nonumber \end{align}

\begin{align} X_{k+N/3} &=\sum_{j<N/3}{x_{3j}W_N^{3j(k+N/3)}}+\sum_{j<N/3}{x_{3j+1}W_N^{(3j+1)(k+N/3)}} +\sum_{j<N/3}{x_{3j+2}W_N^{(3j+2)(k+N/3)}}\nonumber \\ &= \sum{x_{3j}W_{N/3}^{jk}} + W_3^1W_N^k \sum{x_{3j+1}W_{N/3}^{jk}} + W_3^2 W_N^{2k} \sum{x_{3j+2}W_{N/3}^{jk}} \nonumber \end{align}

最低0.47元/天 解锁文章
yu_xianguo
关注
  • 1
    点赞
  • 踩
  • 6
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
专栏目录
快速离散傅里叶变换c语言版
08-27
利用c语言实现的快速离散傅里叶变换，代码都是自己写的，绝无copy!
C++实现一维离散傅里叶变换
ShadowN1ght的博客
08-22 9606
本文介绍如何用C++实现一维离散傅里叶变换（Discrete Fourier Transform）。 1. 一维傅里叶变换公式简介 如果读者对傅里叶变换的认识处于一知半解的程度，那么建议先阅读相关文章，例如网络上广泛流传的《如果看了此文你还不懂傅里叶变换，那就过来掐死我吧》，以加深理解。 一维连续傅里叶变换的公式表达如下： 经过离散化后，上述公式变为： 其中M为x的个
快速傅里叶变换（FFT）的原理、实现及代码解析（附C#源码）
09-09 1万+
一、离散傅里叶变换回顾与FFT的引出 对于长度为N点的数字信号序列 ，定义其离散傅里叶变换为：  我们知道，利用系数 的性质可以大大减少DFT的计算量，这种算法就是快速离散傅里叶变换FFT。需要说明的是，FFT不是一种新的变换，而是一种求DFT的快速计算机算法。 对
离散傅里叶变换(DFT) C代码
07-31
是《数字信号处理C语言程序集》的DFT的C语言源代码，已测试！
C语言直接DFT和FFTW两种方式实现离散傅里叶变换
BloomW的博客
05-16 4681
/*  * dft.c  *  *  Created on: 2018年5月16日  *      Author: ChiWang  */#include &lt;stdio.h&gt; #include &lt;math.h&gt;#define N   16              //采样个数 #define PI  3.1415926535 typedef struct         ...
关于离散傅里叶计算的使用心得
qq_30058483的博客
06-08 129
离散傅里叶使用注意
浅谈快速离散傅里叶变换的实现
dechebiao6765的博客
09-29 351
在运用之前我们需要知道他是什么？是怎么来的？怎么去应用。 傅立叶变换是一种分析信号的方法，它可分析信号的组成成分，也可用这些成分合成信号。许多波形可作为信号的成分，比如正弦波、方波、锯齿波等，傅立叶变换用正弦波作为信号的组成成分，在时域他们是相互重叠在一起的，我们需要运用傅里叶变换把他们分开并在频域显示出来。 连续傅里叶变换（Fourier Transform）如下： ...
第三章 离散傅里叶变换及其快速算法习题答案参考.doc
05-12
第三章 离散傅里叶变换及其快速算法习题答案参考.doc
《数字图像处理》冈萨雷斯 第三版part1(共1-3部分）
11-12
　4.1 二维离散傅里叶变换 　4.2 在MATLAB中计算并可视化二维DFT 　4.3 频域滤波 　4.3.1 基本概念 　4.3.2 DFT滤波的基本步骤 　4.3.3 用于频域滤波的M函数 　4.4 从空间滤波器获得频域滤波器 　4.5 在频域中直接...
并行计算导论（原书第2版）.[美]Ananth Grama(带详细书签).pdf
04-24
本书系统介绍涉及并行计算的体系结构、编程范例、... 13.2.3 并行快速傅里叶变换中的额外计算 13.3 转置算法 13.3.1 二维转置算法 13.3.2 转置算法的推广 13.4 书目评注 习题 附录A 函数的复杂度与阶次分析 索引
IOI国家集训队论文集1999-2019
09-20
任 恺 -《图论的基本思想及方法》 杨 俊 -《二分策略在信息学竞赛中的应用》 张伟达 -《用改进算法的思想解决规模维数增大的问题》 黄 刚 -《数据结构的联合》 杨 弋 -《从“小H的小屋”的解法谈算法的优化》 ...
MATLAB教学视频数模资料及源程序代码下载及说明-MATLAB教学视频，数模资料及源程序代码下载及说明.pdf
08-13
本期视频时长约120分钟，通过三个具体的数学建模案例，详细地讲解了MATLAB中一维插值和二维插值的应用和实现方法。另外，还通过自编程的方式，实现拉格朗日插值方法。视频的最后，还对多维插值做了基本的介绍。 ...
使用fftw3-3进行快速傅里叶变换
libzhouxi的专栏
08-25 815
使用fftw3-3进行快速傅里叶变换
数字信号处理3: 快速傅里叶变换（FFT）（含代码）
zuguorui的专栏
02-09 3855
文章目录1. FFT推导2. FFT为什么快?3. 一些加速措施3.1 查表法计算三角函数3.2 奇偶分解 在之前的文章《傅里叶变换》中，我们已经推导了连续傅里叶变换和离散傅里叶变换。由于计算机的发展，离散傅里叶变换（DFT）可谓是信号处理的杀手锏。但是离散傅里叶变换计算量巨大，通常在实时信号处理时是无法使用的，直到快速傅里叶变换（FFT）算法被发现。 与DFT不同，FFT是一种算法而非理论，因此...
快速DFT计算——DIT基2算法
weixin_42385536的博客
08-19 1566
不同Radix实现方式的快速傅里叶变换复杂度matlab仿真分析,对比基2,基4以及分裂基
fpga/matlab/simulink算法仿真工程
03-14 556
这种分解既有基2的部分，又有基4的部分，因此称为分裂基分解。上面的N/2N/2点DFT 又可分解为一个N/4N/4点的DFT 和两个N/8N/8点的DFT, 而两个N/4N/4点的DFT也分别可以分解为一个N/8N/8点的DFT和两个N/16N/16点的DFT。这样一个N点变换就分解为了两个N/2点变换，这里F 1 ( k ) F_1(k)F 1 (k)和F 2 ( k ) F_2(k)F 2(k)分别是序列x中的奇数号和偶数号序列的N / 2 N/2N/2点DFT变换。
用C语言实现DFT算法
热门推荐
姚明明的博客
01-22 1万+
一. 简介 离散傅里叶变换（Discrete Fourier Transform，缩写为DFT），是傅里叶变换在时域和频域上都呈离散的形式，将信号的时域采样变换为其DTFT的频域采样。在形式上，变换两端（时域和频域上）的序列是有限长的，而实际上这两组序列都应当被认为是离散周期信号的主值序列。即使对有限长的离散信号作DFT，也应当将其看作其周期延拓的变换。 X（k）=∑n=0N−1x(n)exp(−...
课设复习之DFT的快速实现
dreamer5z的博客
05-24 1324
代码核心思想参考数字信号处理(第三版) 高西全 丁美玉 编著 第四章 快速傅里叶变换 FFT主要流程图如上图所示，主要包括三层循环，最外层为级循环，即计算基2FFT每一级的蝶形运算，中间为旋转因子的循环，即该级下所有旋转因子种类，内部为该旋转因子下计算的蝶形运算。书中所讲步骤较为详细，此处不...
证件图像校正
YU_Xianguo的博客
03-20 6565
引子： 你是否遇到过这种情况：用手机对着身份证拍张照片，然后想打印出扫描件的效果？应用商城里有个热门app好像叫“扫描全能王”，大概提供了这个功能，但是它是收费的，普通用户用起来有些许不便。最近我也遇到了这个问题，别人用手机拍了个身份证照片传给我，如上面所示，让我把它打印出来：我最先想用PS，但是我发现无法用一个规则的矩形把它框得很好，我那没入门的PS技术无法解决这个问题，于是我就花了一晚上...
已知一个数组形式的LSF(filtered_lsf[i]，长度为filtered_lsf_size),如何不使用第三方库的情况下用C++代码实现对LSF归一化离散傅里叶变换得到sfr曲线
最新发布
05-26
LSF（Line Spectral Frequencies）是一种表示语音信号的方法，它是由线性预测分析（Linear Predictive Analysis，LPA）求得的。归一化离散傅里叶变换（Normalized Discrete Fourier Transform，NDFT）是一种将时域信号转换为频域信号的方法。在语音信号处理中，基于LSF的NDFT可以用来计算声道频率响应（Spectral Frequency Response，SFR）曲线，用于声学特征提取和语音识别。 下面是一个基于LSF的NDFT计算SFR曲线的C++代码实现： ```c++ #include <iostream> #include <cmath> #define PI 3.14159265358979323846 int main() { const int filtered_lsf_size = 10; // LSF长度 const double filtered_lsf[filtered_lsf_size] = {0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, 1.0}; // LSF数组 const int sfr_size = 512; // SFR长度 double sfr[sfr_size]; // SFR数组 // 归一化LSF double normalized_lsf[filtered_lsf_size]; for (int i = 0; i < filtered_lsf_size; i++) { normalized_lsf[i] = filtered_lsf[i] * PI; } // NDFT计算SFR for (int i = 0; i < sfr_size; i++) { double omega = i * 2 * PI / sfr_size; double re = 0.0, im = 0.0; for (int j = 0; j < filtered_lsf_size; j++) { re += cos(omega * normalized_lsf[j]); im -= sin(omega * normalized_lsf[j]); } sfr[i] = sqrt(re * re + im * im); } // 输出SFR for (int i = 0; i < sfr_size; i++) { std::cout << sfr[i] << std::endl; } return 0; } ``` 在上面的代码中，首先定义了LSF数组和长度，以及SFR数组和长度。为了归一化LSF，将其每个元素都乘以PI得到normalized_lsf数组。然后，对于每个频率点i，计算其对应的角频率omega，并通过循环计算NDFT中的实部和虚部。最后，将实部和虚部的平方和开方得到SFR曲线的幅度值。 需要注意的是，SFR曲线的长度应当与输入信号的长度相同，通常为512或1024。此外，由于NDFT计算复杂度较高，需要进行优化或使用快速傅里叶变换（Fast Fourier Transform，FFT）算法。

“相关推荐”对你有帮助么？

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值