步骤 5：关于苹果和桔子
如果你一直跟着我，我们几乎完成了通向傅里叶变换的旅程我们学了需要多少正弦波，它的数量依 赖于我们查看的采样的数量有一个频率上下限界，并且不知道怎么确定单个分量的幅度以完成我们的配 方我们一直不清楚究竟如何从我们的采样来确定实际的配方直观上我们可以断定能找到正弦波的幅度， 设法把一个已知频率正弦波和采样作对比，我们测量找出它们有多么接近如果它们精确地相等，我们知 道该正弦波存在着相同的幅度，如果我们发现我们的信号与参考正弦波一点也不匹配，我们将认为这个不 存在尽管如此，我们如何高效地把一个已知的正弦波同采样信号进行比较？幸运地是，数字信号处理工 作者早已解决了如何作这些事实上，这象加法和乘法一样容易－我们取一个已知频率的单位正弦波（这 意味着它的振幅是 1，可从我们的计算器或计算机中精确地获得）和我们的信号采样相乘累加乘积之后， 我们将得到我们正在观测的这个频率上正弦波分量的幅度这是个举例，一个简单的完成这些工作的 C 代 码片段：
Listing 1.1: The direct realization of the Discrete Sine Transform (DST):
#define M_PI 3.14159265358979323846
long bin,k; double arg;
for (bin = 0; bin < transformLength; bin++) {
transformData[bin] = 0.;
for (k = 0; k < transformLength; k++) {
arg = (float)bin * M_PI *(float)k / (float)transformLength; transformData[bin] += inputData[k] * sin(arg);
}
}
这段代码变换存储在 inputData[0...transformLength-1]中我们测量的采样点成为一个正弦波分量的幅 度队列 transformData[0...transformLength-1]根据通用术语，我们称参考正弦波的频率步长为盒（bin）， 这意味着它们被认为象是一个我们放置我们估计的任意分量波的幅度的容器离散正弦变换(DST)是一个
普通程序，它假设我们无法想象我们的信号看起来象什么样，否则我们能使用一个更加高效率的方法来确 定正弦波分量的幅度（例如，我们予先知道, 我们的信号是一个已知频率的正弦波我们能直接地找出它 的高度而不用计算正弦波的整个范围实现这个有效的逼近是基于傅里叶原理，它能在文献的戈策尔(Goe rtzel)算法条目下找到）
这些就是你坚持想要的我们为什么用这样的方法计算正弦变换的一个解释：对我们用一个已知频率正弦 波的乘积来作一种非常直观逼近的理由，可以设想，这大致相当于一个固有频率的“共振”在系统内发生时 物理世界发生的事情sin(arg)项本质上是一个获得由输入信号波形激励的谐振器如果输入（信号）有在 我们正观测的频率上的分量，它的输出将是参考正弦波谐振的幅度因为我们的参考波是单位幅度的，输 出是一个在那个频率上的分量的实际幅度的一个直接测量因为谐振器只是简单的滤波器，变换（不可否 认是在稍微宽松条件下）被认为有极窄的带通滤波器组的特征，它位于我们估值的频率中心的周围这有 助于解释一个事实，为什么傅立叶变换提供了对信号进行过滤的一个高效工具
只是为了完备性: 当然，上述程序是可逆的，当我们知道它的正弦波分量时，我们的信号(在数字精确度 极限内)能完全被重建，通过简单地把正弦波加起来这留下给读者做为一个练习同样程序能改变使用余 弦波做为基本函数工作-我们只需简单地改变 sin(arg)条件到 cos(arg)来获得离散余弦变换的直接实现(DC T)
现在，就象在这篇文章较前面的段落中我们讨论过的那样，我们在实践中没有办法区分一个被测量的正 弦类函数象是正弦波还是余弦波做为代替我们总是测量正弦信号，且正弦和余弦变换在实践中没有太大 的用途，除了一些特殊情况（象图象压缩的地方，即毎块图象具有能用一个基本的余弦或正弦函数较好模 拟特性，例如能用余弦基本函数较好表现的相同颜色的大区域）正弦信号是一个比正弦或余弦波更一般 的片断，因为它可以开始在在它的周期中的一个任意位置我们记得，当余弦波开始于 1 时，正弦波总开 始在 0 当我们采取正弦波作为参考，余弦波开始在它的周期的最后 1/4 之处一般用度或弧度测量它们 的偏移量，这是两个一般与三角函数相关的单位一个完整的周期等于 360°（代表 \度\或 2π 个弧度(代 表\π\，\π\发音象“pie”π 是希腊字表示数 3.14159265358979323846... 在三角学方面有重要意义) 余弦波因而有一个 90°或 π/2 的偏移这偏移叫正弦信号的相位，因此余弦波相对正弦信号有 90°或 π/2 相 位
相位的事情就有这些内容因为我们一直不能限定信号在 0°或 90°相位开始（因为我们正观测一个我们 可能无法控制的信号），它对同时直接唯一的描述信号的频率振幅相位至关重要以正弦或余弦做变
换相位限制在 0°或 90°，一个具有任意相位的正弦信号将引起相邻频率出现假峰（因为它们试图“帮助”分 析，强制给被测信号加上一个 0°或 90°的相位作用）它有些象用一圆石头去填满一个方孔：你需要小一 些的圆石头去填充剩余的空间，并且更小的石头填好依然留出空的空间，等等我们所需要的是能处理一 般信号的变换，它能处理任意相位正弦波构成的信号
步骤 6：离散傅叶变换
从正弦变换到傅里叶变换的步骤是简单的，只需用更一般的方法在正弦变换中对每个频率上的测度 使用正弦波，在傅里叶变换中正弦余弦波二者都使用就是说，对任意的当前频率，我们以同一频率的 正弦和余弦波来“比较”（或“共振”）被测信号如果我们的信号看起来很象正弦波，变换的正弦部份将有 一个大的幅值如果它看起来象余弦波，变换的余弦部份将有一个大的幅值如果看起来象反相的正弦波， 也就是说，它开始于 0 但下降至－1 取代上升至 1，它的正弦部份将有一个大的负幅值这表明用＋－ 符号和正弦余弦相位能表示任意给定频率的正弦信号[2]
Listing 1.2: The direct realization of the Discrete Fourier Transform[3]:
#define M_PI 3.14159265358979323846
long bin, k;
double arg, sign = -1.; /* sign = -1 -> FFT, 1 -> iFFT */ for (bin = 0; bin <= transformLength/2; bin++) {
cosPart[bin] = (sinPart[bin] = 0.); for (k = 0; k < transformLength; k++) {
arg = 2.*(float)bin*M_PI*(float)k/(float)transformLength; sinPart[bin] += inputData[k] * sign * sin(arg); cosPart[bin] += inputData[k] * cos(arg);
}
}
我们仍遗留着一个问题，就是如何获得傅里叶变换所缺乏的那些有用的东西我说过傳里叶变换的优赿性 超过正弦和余弦变换是因为用正弦信号工作但至今我们还未看到任何正弦信号，仍只有正弦和余弦好， 这需要一点附加处理步聚：
#define M_PI 3.14159265358979323846
long bin;
for (bin = 0; bin <= transformLength/2; bin++) { /* frequency */
frequency[bin] = (float)bin * sampleRate / (float)transformLength;
/* magnitude */
magnitude[bin] = 20. * log10( 2. * sqrt(sinPart[bin]*sinPart[bin] + osPart[bin]) / (float)transformLength);
cosPart[bin]*c
/* phase */
phase[bin] = 180.*atan2(sinPart[bin], cosPart[bin]) / M_PI - 90.;
}
在运行清单 1.3 所示的关于 DFT 输出的代码段之后，我们结束被看作以正弦信号波的和的输入信号表示 K 序正弦信号是用 frequency[k],magnitude[k]和 phase[k]来描述的单位是 Hz(Hertz,周/秒)，dB(Decibel), 和°(Degree)请注意在经过清单 1.3 的后加工（处理）即把正弦和余弦函数部份转换成一个单一的正弦信 号之后，我们命名 K 序正弦信号的振幅－DFT 存贮为幅度，且它总是取相对值我们可以说一个-1.0 的振 幅对应于 1.0 的幅度，对应于相位+或-180°在文献中，做做傅里叶变换的场合，队列 magnitude[] 被称 作被测信号的幅度谱，队列 phase[]被称作被测信号的相位谱
如用分贝测量存贮幅度的参考，输入波也期望有一个在[-1.0,1.0]之间的采样值，相对于 0dB 幅度满刻 度数字做为一个 DFT 的有趣应用，比如清单 1.3 就可被用于写一个基于离散傅里叶变换的谱分析
结论
象我们已知那样，傅里叶变换和其系列的离散正弦和余弦变换，提供了把一个信号分解成一束分波的 便利工具结果有正弦或余弦之一，或正弦信号（用正弦和余弦波的组合来描述）在傅里叶变换中同时 使用正弦和余弦波的好处是我们因而能引入相位的概念，它使变换更一般化，因而我们能用它有效清楚地 分析既不是纯正弦也不是纯余弦的正弦信号，当然其它信号也一样
傅里叶变换与被考察信号无关，因而无论我们正分析的信号是一个正弦信号或是一些其它的更复杂的， 变换需要相同的操作数这就是为什么傅里叶变换被称做无参数变换的原因，这意味着它对需要的信号“智 能的”分析没有直接的帮助（在考察一个我们已知是一个信号是正弦曲线的情况下，我们更喜欢精确地获 得关于相位，频率，幅度的信息以代替一串在一些预定频率上的正弦和余弦波）
现在我们也知道了我们是在求输入信号在一组固定频率栅格上的值，输入信号实际存在的频率组在这 组栅格上可能不起作用我们在分析中利用的栅格是人为的，因为我们几乎按照关于它们的频率的尝试来 选择参考正弦余弦波说到了这些，马上清楚了一个将要很容易遇到的要点，即被测信号的频率位于变 换栅格的频率之间因此，有一个频率发生在位于两个频率栅格之间的正弦曲线，在变换中将不好被描述 包围着与输入信号频率最接近的栅格的相邻的栅格将试图?改正?频率的背离因而，输入信号的能量将拖 尾至数个相邻的栅格这也是傅里叶变换不能迅速地分析声音返回它的基波和谐波（并且，这也是为什么 我们称正弦和余弦波为分波而不谐波和泛音）
简单的说，没有进一步的快速处理，DFT 和一个狭窄的坝一样，细小并行的带通滤波器组（“通道”） 和每个通道带有附加的相位信息这对分析信号做滤波器和运用其它的技巧是有益的（改变一个信号的 音调而不改变它的速度是它们其中之一，说明在 DSPdimension.com 上另一篇不同的文章中），但它需 要对少量普通任务附加快速处理同样，它能被认为是使用除了正弦和余弦波基本函数的变换系列的一个 特例在这个方向上展开概念超出了这篇文章的范围
最后，重要的是要提及一个更高效的 DFT 工具，也就是一个被称做快速傅里叶变换的算法它最初 是由库利和图克在 1969 年构思的（它的根源仍然要追溯到高斯和其它人的工作）FFT 只是一个高效的 算法，它比上面给出的以直接逼近计算 DFT 所化的时间少，它是结果完全相同的其它方法无论如何，F FT 是以库利/图克算法实施的，它需要变换长度是 2 的幂在实践中，对大多数应用来说这是一个可以接 受的限制有大量的以不同方法实施 FFT 的可利用的文献，因而，可以说有足够多不同的 FFT 实现，其 中一些
一天征服傅里叶变换 - 图文
