第四章正交变换1介绍.pptVIP

第四章 正交变换 * 进行归一化 * 进行归一化 * * * * * * * * * * * * * * * * * * * * 进行归一化 * 进行归一化 * 进行归一化 * 进行归一化 * 进行归一化 * 进行归一化 * 进行归一化 * 进行归一化 * 进行归一化 * 进行归一化 * 进行归一化 * 进行归一化 * 进行归一化 * 进行归一化 * 进行归一化 * 进行归一化 * 进行归一化 * 进行归一化 * 进行归一化 * 进行归一化 * 进行归一化 * 进行归一化 * 进行归一化 * 进行归一化 * 进行归一化 * 进行归一化 * 进行归一化 * 进行归一化 30 傅里叶变换（离散） 31 傅里叶变换（离散） 利用计算机实现快速傅里叶变换（FFT） 利用FFT蝶式流程图算法在计算机上实现快速傅里叶变换 必须解决如下问题： 1 、迭代次数r的确定； 2 、对偶节点的计算； 3 、加权系数 的计算； 4 、重新排序问题 第四章 正交变换 32 傅里叶变换（离散） 1 迭代次数r的确定 由蝶式流程图可见，迭代次数r与N有关。值可由下式确定 式中 N 是变换序列的长度。对于前述基数2的蝶式流程图是2的整数次幂。例如，序列长度为8则要三次迭代，序列长度为16时就要4次迭代等等。 第四章 正交变换 33 傅里叶变换（离散） 2 对偶节点的计算 在流程图中把标有 的点称为节点。其中下标 为列数，也就是第几次迭代，例如， 则说明它是第一次迭代的结果。 代表流程图中的行数，也就是序列的序号数。其中每一节点的值均是用前一节点对计算得来的。例如， 和 均是x 0 和x 4 计算得来的。在蝶式流程图中，把具有相同来源的一对节点叫做对偶节点。 第四章 正交变换 34 傅里叶变换（离散） 2 对偶节点的计算 如: 和 就是一对对偶节点，因为它们均来源于x 0 和x 4 。对偶节点的计算也就是求出在每次迭代中对偶节点的间隔或者节距。由流程图可见，第一次迭代的节距为 ，第二次迭代的节距为 ，第三次迭代的节距为 等等。由以上分析可得到如下对偶节点的计算方法。 如果某一节点为 ，那么，它的对偶节点为 l是表明第几次迭代的数字，k是序列的序号数，N是序列长度。 第四章 正交变换 35 傅里叶变换（离散） 例：如果序列长度N 8，求 的对偶节点。 其正确性不难由流程图来验证。 利用上式可得： 第四章 正交变换 36 傅里叶变换（离散） 3 加权系数 的计算 的计算主要是确定p值。p值可用下述方法求得。 （１）把k值写成r位的二进制数（k是序列的序号数，r是迭代次数）； （２）把这个二进制数右移r-l位，并把左边的空位补零（结果仍为r位）； （３）把这个右移后的二进制数进行比特倒转； （４）把这比特倒转后的二进制数翻成十进制数就得到p值。 第四章 正交变换 37 傅里叶变换（离散） 例子：求 的加权系数 由 和 可知k 2， 2，n 8，则 （１）因为k 2，所以写成二进制数为010； （２）r- 3-2 1，把010右移一位得到001； （３）把001做位序颠倒，即做比特倒转，得到100； （４）把100译成十进制数，得到4，所以p 4， 的加权值为 。 第四章 正交变换 38 傅里叶变换（离散） 结合对偶节点的计算，可以看出 具有下述规律：如果某一节点上的加权系数为 ，则其对偶节点的加权系数必然是 而且 ，所以一对对偶节点可用下式计算 第四章 正交变换 39 傅里叶变换（离散） 4 重新排序 由蝶式流程图可见，如果序列 是按顺序排列的，经过蝶式运算后，其变换序列 是非顺序排列的，即乱序的； 具体排序方法如下： （１）将最后一次迭代结果 中的序号数k写成二进制数，即 （２）将r位的二进制数比特倒转，即： 也就是 （３）求出倒置后的二进制数代表的十进制数，就可以得到与 相对应的 的序号数 第四章 正交变换 40 傅里叶变换（离散） 例子： N 8 的最后迭代结果： x3 0 →000→倒置→000→十进制（0） x3 1 →001→倒置→100→十进制（4） x3 2 →010→倒置→010→十进制（2） x3 3 →011→倒置→110→十进制（6） x3 4 →100→倒置→001→十进制（1） x3 5 →101→倒置→101→十进制（5） x3 6 →110→倒置→011→十进制（3） x3 7 →111→倒置→111→十