matlab用DFT分析三频频谱,Matlab在DFT频谱误差分析中的应用

自然科学学科研究ＮａｔｕｒａｌＳｃｉｅｎｃｅＤｉｓｃｉｐｌｉｎｅｒｅｓｅａｒｃｈ

Ｍａｔｌａｂ在ＤＦＴ频谱误差分析中的应用

吴炎波邓冠群

(华中师范大学信息技术系湖北?武汉４３００７９)

摘要研究数字信号最有效的方法通常是离散傅里叶变换(ＤＦｒ)，对信号进行ＤＦＴ时会存在各种误差．本文针对ＤＦＴ分析的主要误差来源——栅栏效应和截断效应，在理论分析的基础上，采用ＭＡＴＬＡＢ仿真，更加直观的说明了补零只能减轻栅栏效应，使频谱更加细致，只有增加有效信号的截取长度，才能提高频率分辨率，减轻截断效应．最后给出如何选择适合截断长度的结论．

关键词ＭＡＴＬＡＢＤＦＴ谱分析栅栏效应截断效应

中图分类号：ＴＰ３１文献标识码；Ａ

Ａｐｐｌｉｃａｔｉｏｎｏｆ

ＭａｔｌａｂｉｎＥｒｒｏｒＡｎａｌｙｓｉｓｏｆＤＦＴＦｒｅｑｕｅｎｃｙＳｐｅｃｔｒｕｍ

ＷＵＹａｎｂｏ，ＤＥＮＧＧｕａｎｑｕｎ

(Ｄｅｐｔ．ｏｆＩｎｆｏｒｍａｔｉｏｎＴｅｃｈｎｏｌｏｇｙ，ＣＣＮＵ，Ｗｕｈａｎ，Ｈｕｂｅｉ４３００７９)ＡｂｓｔｒａｃｔＤＦＴｉｓａｌｗａｙｓｔｈｅｍｏｓｔｅｆｆｉｃｉｅｎｔｍｅｔｈｏｄｏｆｄｉｇｉｔａｌｓｉｇｎａｌｐｒｏｃｅｓｓｉｎｇ，ｂｕｔｉｔｂｒｉｎｇｓｋｉｎｄｓｏｆｅｒｒｏｒｓ．ＡｃｃｏｒｄｉｎｇｔｏｔｈｅｍａｉｎｓｏｕｒｃｅｏｆｅｎＤｒ—一ｆ钿ｃｅｅｆｆｅｃｔａｎｄｔｒｕｎｃａｔｅｄｅｆｆｅｃｔ，ｔｈｉｓａｒｔｉｃｌｅｂａｓｅｄｏｎＭａｔｌａｂｓｉｍｕｌａｔｉｏｎ．ｍｏｒｅｉｎｔｕｉｔｉｖｅｉｌｌｕｓｔｒａｔｅｓｔｈａｔｆｉｌｌｉｎｇｚｅｆｏｃａｒｌｏｎｌｙｅａｓｅｆｅｎｃｅｅｆｆｅｃｔ，ｍａｋｅｔｈｅｓｐｅｃｔｒｕｍ１１１０１忙ｍｅｔｉｃｕｌｏｕｓｌｙ，ｏｎｌｙｉｎｅｒｅａｓｅｉｎｇｅｆｆｅｃｔｉｖｅｓｉｇｎａｌｉｎｔｅｒｃｅｐｔｌｅｎｇｔｈｃａｎｉｍ—ｐｒｏｖｅｔｈｅｆｒｅｑｕｅｎｃｙｒｅｓｏｌｕｔｉｏｎ，ｒｅｄｕｃｅｔｒｕｎｃａｔｉｏｎｅｆｆｅｃｔ．Ａｔｌａｓｔ，Ｐｒｅｓｅｎｔｔｈｅｃｏｎｃｌｕｓｉｏｎｏｆｈｏｗｔｏ∞ｌｅｃｔｓｕｉｔａｂｌｅｔｒｕｎｃａｔｉｏｎｌｅｎｇｔｈ．Ｋｅｙ

ｗｏｒｄｓＭＡＴＬＡＢ，ＤＦＴｓｐｅｃｔｒｕｍａｎａｌｙｓｉｓ，Ｆｅｎｃｅｅｆｆｅｃｔ．Ｔｒｕｎｃａｔｉｏｎｅｆｆｅｃｔ

Ｉ栅栏效应与“高密度频谱”

Ｎ点ＤＦＴ是在频率区间【Ｏ，２ｎｌ上对时域离散信号的频谱进行Ｎ点等间隔采样，而采样点之间的谱线是看不到的。这就好像从Ｎ个栅栏缝隙中观看信号的频谱情况，仅得到Ｎ个缝隙中看到的频谱样值函数值，这种现象就成为栅栏效应。如果相邻采样点之间距离比较大的话，栅栏效应有可能漏掉大的频谱分量，使得到的频谱与实际情况不相符。针对这种情况，笔者力求寻找出一种方法得到更加细致的谱线。

很容易想到，在不改变序列长度Ｎｏ的情况下，增加采样点数Ｎ，就可以使相邻“栅栏”问的距离减小。

例如：取序列ｘ(ｎ)＝ｓｉｎ(０．４６ｎ)＋ｃｏｓ(０．５０ｎ)(取Ｎ尹ｌＯ)，对该序列做Ｎ＝１０的ＤＦＴ，结果如图ｌ所示。其中ＤＦＴ变换后，频率∞＝兰等。因为《ｎ)由两个频率ｆ，０１＿０．４６兀，ｏ垆ｏ．５０兀厶成，１０点ＤＦＴ的频率采样间隔ｄ＝２ｘ．／１０＝０．２，而需要分辨的频率间隔为△∞＝纰?∞，＝０．０４ｎ，所以无法分辨这两个频率分量。由图Ｉ可知，ＤＦＴ变换后频率集中在(ｏ＝』：；等０．６ｎ(ｋ＝３)处，与实际情况不符。

为使频谱更加接近实际情况，在不改变序列长度ＮＯ的情况下，增加采样点数Ｎ，我们在ｘ(ｎ)(Ｎｏ＝１０)的后面补９０个零，构成长度Ｎ－－１００的序列，然后做１００点ＤＦｒ变换，结果如图２所示。

可以看出，图２比图ｌ的频谱更加清晰，频谱密度明显增大，这是由于我们增加了采样点数，减小了“栅栏”间的距离引起的。此时频率采样间隔ｄ＝２，ｃ，ｌｏｏ＝０．０２ｎ，小于需要分辨的频率间隔Ａｅ＝０．０４ｎ，但是ＤＦＴ变换后的频谱频率集中在∞；』笔争０．５，ｃ(ｋ＝２５)处，虽然稍微向实际情况靠近了一些，但仍与实际情况不符。

于是引出如下定义：若频率采样间隔ｄ≤△(Ｉ)仍不能得到分辨频率间隔△(Ｉ)的频谱，则将该频谱定义为高密度频谱。图２就是“高密度频谱”。由以上的分析可知。高密度频谱是改善栅栏效应的有效方法。

２截断效应与“高分辨率频谱”

实际中遇到的信号ｘ(ｎ)可能是无限长的，用ＤＦＴ对其进行谱分析时，必须将其截断，形成有限长度的序列ｙ(ｎ)＝ｘ(ｎ)ｗ(ｎ)，其中ｗ(ｎ)是长度为Ｔｐ窗函数，由卷积定理可知，ｙ(ｎ)频谱是ｘ(ｎ)的频谱Ｘ(ｅ一)与ｗ(ｎ)的频谱ｗｏｏ)的卷积。以“ｎ)取矩形窗为例，取Ｘ(ｎ)＝ＣＯＳ(Ｏ．５０ｎｎ)，其加窗前和加窗后的幅频特性如图３所示。从图中可以看出，加窗前，ｘ(ｎ)是离散谱线，加窗后，原来的离散谱线向附近展宽(通常称之为泄露)，使频谱变得模糊，分辨率降低。同时，主谱线两边形成

很多旁瓣，引起不同频率分量