matlab 基2时域抽取法,MATLAB辅助DSP实现基2时域抽取法FFT

微机保护

第２１卷第５期 制造业信息化２００８年９月

机电产品开发与创新

Ｄｅｖｅｌｏｐｍｅｎｔ＆ＩｎｎｏｖａｔｉｏｎｏｆＭａｃｈｉｎｅｒｙ＆ＥｌｅｃｔｒｉｃａｌＰｒｏｄｕｃｔｓ

Ｖｏｌ．２１，Ｎｏ．５

Ｓｅｐ．，２００８

ＭＡＴＬＡＢ辅助ＤＳＰ实现基２时域抽取法ＦＦＴ

郭小莉，黄钉劲，阮照军

(西安工业大学光电工程学院，陕西西安７１００３２)

摘要：介绍了基２时域抽取法ＦＦＴ的原理和算法，并在ＭＡＴＬＡＢ仿真软件的辅助下、在数字信号处理

ＤＳＰ上实现。

关键词：时域抽取；ＤＳＰ；ＦＦＴ；仿真；ＭＡＴＬＡＢ中图分类号：ＴＰ３９１．４

文献标识码：Ａ

文章编号：１００２－６６７３(２００８)０５－１３４－０２

０引言

对于确定时域信号的频谱来说，通常采用离散傅立叶变换ＤＦＴ的方式来实现。随着微电子技术和计算机技术的发展，数字信号处理ＤＳＰ得以广泛应用。实际中，在ＤＳＰ的软件中实现信号从时域到频域的转换并未应用ＤＴＦ而是快速傅立叶变换ＦＦＴ，这是因为ＦＦＴ的输出与ＤＦＴ的输出相同，但运算量要小得多。因此采用ＤＳＰ实现ＦＦＴ，可进行信号高速的实时处理。在具有强大数值计算和图形、图像显示能力的ＭＡＴＬＡＢ软件的辅助下实现ＦＦＴ，缩短设计周期，验证结果，大大提高设计效率。

Ｘ(ｋ)＝!ｘ(ｎ)ＷＮ＝!ｘ(２ｍ)ＷＮ

ｎ＝０

ｍ＝０

Ｎ－１

ｋｎ

(Ｎ／２)－１

２ｍｋ

＋!ｘ(２ｍ＋１)ＷＮ

ｍ＝０

(Ｎ／２)－１

(２ｍ＋１)ｋ

(Ｋ＝０，１，．．．，Ｎ－１)(４)

因为：

ＷＮ＝ｅ

２ｍｋ

－ｊ２!２ｍｋ＝ｅ

－ｊ２!ｍｋ＝ＷＮ／２

ｍｋ

(５)

所以式(４)可写成：

Ｘ(ｋ)＝!ｘ１(ｍ)ＷＮ／２＋ＷＮ!ｘ２(ｍ)ＷＮ／２(ｋ＝０，１，．．．Ｎ－１)

ｍ＝０

ｍ＝０

(Ｎ／２)－１

ｍｋｋ

(Ｎ／２)－１

ｍｋ

(６)

式(６)便是基２时域抽取法ＦＦＴ的结果，它是把长度的序列一分为二，将点ＤＦＴ表示成为两个Ｎ／２点

１基２时域抽取法ＦＦＴ的原理及设计方法

１．１基２时域抽取法ＦＦＴ的原理

设输入时域信号为ｘ(ｎ)，其转换后的频域信号为Ｘ(ｎ)，则离散傅立叶变换ＤＦＴ为：

Ｎ－１

－ｊ２!ｋｎ

ＤＦＴ的线性组合。然后再把Ｎ／２点ＤＦＴ一分为二，表示

为两个Ｎ／４点的ＤＦＴ。如此重复下去，直至分解成两点

ＤＦＴ的运算，两点ＤＦＴ实际上只是加减运算。这就是基２时域抽取ＦＦＴ算法的原理

［１］

。

１．２ＭＡＴＬＡＢ设计仿真

现就使用ＭＡＴＬＡＢ软件来实现如下所示的时域信

(１)

号ｘ(ｔ)的频谱。设此时域信号为ｘ(ｔ)＝０．０７５ｃｏｓ(６００!)，对此信号以３０ＫＨｚ采样频率进行采集得到５１２个采样点，则ｘ(ｔ)经采样后为ｘ(ｎ)＝０．０７５ｃｏｓ(ｎ２!)。按照这一要

Ｘ(ｎ)＝!ｘ(ｎ)ｅ

ｎ＝０

(ｋ＝０，１，．．．，Ｎ－１)

在基２时域抽取法中，将此时域信号为ｘ(ｎ)根据采样号ｎ分解为偶采样点和奇采样点，则偶数序列为：

ｘ１(ｍ)＝ｘ(２ｍ)，ｍ＝０，１，．．．，(Ｎ／２)－１

奇数序列为：

(２)(３)

求，通过ＭＡＴＬＡＢ进行设计

［２］

，其程序如下：

ｘ２(ｍ)＝ｘ(２ｍ＋１)，ｍ＝０，１，．．．，(Ｎ／２)－１ＤＦＴ。即：

收稿日期：２００８－０８－０７

Ｎ＝５１２；

ｘ＝０．０７５＊ｃｏｓ(２＊ｐｉ＊(０：Ｎ－１)／１００)；ｓｕｂｐｌｏｔ(２，１，１)；ｐｌｏｔ(ｘ(１：Ｎ))；ｙｌａｂｅｌ(＇原始信号＇)Ｎ＝０：Ｎ－１；Ｘ＝ｆｆｔ(ｘ)ｓｕｂｐｌｏｔ(２，１，２)；

对于式(１)表示的ＤＦＴ，便可分成两个长度为Ｎ／２的

作者简介：郭小莉(１９８３－)，女，在读硕士研究生。主要研究方向：精密仪器及机械；黄钉劲(１９６５－)，男，副教授，硕士生导师。主要研究方向：测控技术。

１３４