正余弦信号的谱分析课程设计.doc

正余弦信号的谱分析课程设计

PAGE

PAGE3

目录

TOC\o1-3\h\u14063目录 2

22984摘要 3

7941一．DFT的简介 4

194261.1概述 4

321871.2．DFT的定义 4

297811.3谱分析的原理 5

3017二．用DFT对连续信号进行谱分析 6

13460三．用DFT进行谱分析的误差问题 6

25471.混叠现象 6

7762.栅栏效应 7

130263.截断效应 7

32243四．设计实现 8

187001.设计内容 8

282372.用Matlab软件实现 8

2200五．结果分析 13

641六．结束语 14

17819七致谢 14

摘要

数字信号处理方法的一个重要用途是在离散时间域中确定一个连续时间信号的频谱，通常称为频谱分析，更具体的说它也包括能量谱或功率谱，所谓信号的谱分析就是计算信号的傅里叶变换，而DFT的实质是有限长序列傅里叶变换的有限点离散采样，从而实现了频域离散化，使数字信号处理可以在频域采样数值运算的方法进行，这样就大大提高了数字信号处理的灵活性，从而使信号的实时处理和设备的简化得以实现。利用Matlab软件对正余弦信号进行设计程序分析并画出频谱图，所以说DFT不仅在理论上有重要意义，而且在各种信号的处理中亦起着核心的作用，数字频谱分析可以应用在很广的领域。

关键字：Matlab频谱分析DFT

折叠频率的频率成分，以免发生频谱混叠现象。

2.栅栏效应

N点DFT是在频率区间［0,2］上对时域离散信号的频谱进行N点等间隔采样，而采样点之间的频谱是看不到的。这就好像从N个栅栏缝隙中观看信号的频谱情况，仅得到N个缝隙中看到的频谱函数值，这就是栅栏效应。由于栅栏效应可能漏掉大的频谱分量。故对于有限长序列，可以在原序列尾部补零；对于无限长序列，可以增大截取长度及DFT变换区间长度，从而使频域采样间隔变小，增大频域采样点数和采样点位置，使原来漏掉的某些频谱分量被检测出来。

3.截断效应

实际中遇到的序列x（n）可能是无限长的，用DFT对其进行谱分析师必须将其截断，形成有限长序列y（n）=x(n)w(n),长度为N。w（n）=,称为矩形窗函数。截断后对谱分析的影响主要表现在以下两点：

（1）.泄露：原来序列x（n）的频谱是离散谱线，经截断后，是原来的离散谱线向附近展宽，通常称这种展宽为泄露。泄露可以是频谱变模糊，使谱分辨率降低。

（2）.普间干扰：在主谱线两边形成很多旁瓣，引起不同频率分量间的干扰，特别是强信号谱的旁瓣可能湮没弱信号的主谱线，或者把强信号谱的旁瓣误认为是另一频率的信号的谱线，从而造成假信号，这样就会使谱分析产生较大偏差。

截断效应就是有以上两种影响对信号截断引起的。

四．设计实现

1.设计内容

（1）对一个频率为10Hz，采样频率为64Hz的32点余弦序列进行谱分析，画出其频谱图；若将频率改为11Hz，其他参数不变，重新画出该序列的频谱图，观察频谱泄漏现象，分析原因；

（2）考察DFT的长度对双频率信号频谱分析的影响。设待分析的信号为

令两个长度为16的正余弦序列的数字频率为及。取N为四个不同值16，32，64，128。画出四个DFT幅频图，分析DFT长度对频谱分辨率的影响。

（3）在上题中若把两个正弦波的频率取得较近，令，，试问怎样选择FFT参数才能在频谱分析中分辨出这两个分量？

2.用Matlab软件实现

（1）当频率为10Hz时

F=input(输入信号频率);t=0:0.001:0.2;x1=cos(2*pi*F*t);subplot(3,1,1);plot(t,x1);

title(x1连续余弦信号);n=0:31;

x2=cos(2*pi*F*n*1/64);

subplot(3,1,2),stem(n,x2);xlabel(n),ylabel(x1(n));title(x2采样后的余弦信号);k=0:31;

X=abs(fft(x2,32));subplot(3,1,3);stem(k,X);

xlabel(k),ylabel(X(k));

string=[num2str(32),点FFT幅频曲线];title(strin);

当频率为11Hz时

F=input(输入信号频率);t=0:0.001:0.2;x1=cos(2*pi*F*t);subplot(3,1,1);plot(t,x1);

title(x1连续余弦信号);n=0:31;

x2=cos(2*pi*F*n*1/64);