河海大学sp实验3.docVIP

实验三 用FFT对信号进行频谱分析 班级：2012级通信一班 一、实验目的 1、加深对离散信号DFT的理解； 2、掌握FFT算法的流程及其MATLAB实现； 3、利用FFT对典型信号进行频谱分析； 4、结合理论知识，对频谱分析中出现的有关现象进行理论分析。 实验原理 混叠：采样序列的频谱是被采样信号频谱的周期延拓，当采样频率不满足奈奎斯特采样定理的时候，就会发生混叠，使得刺痒后的序列信号的频谱不能真实的反映原采样信号的频谱。 泄露：根据理论分析，一个时间的信号其频带宽度为无限，一个时间无限的信号其频带宽度则为有限。因此对一个时间有限的信号，应用DFT进行分析，频谱混叠难以避免。对一个时间无限的信号虽然频带有限，但在实际运算中，时间总是取有限值，在将信号截断的过程中，出现了分散的扩展谱线的现象，称之为频谱泄露或功率泄露。 栅栏效应：DFT是对单位圆上Z变换的均匀采样，所以它不可能将频谱视为一个连续函数，就在一定意义上看，用DFT来观察频谱就好象通过一个栅栏来观看一个景象一样，只能在离散点上看到真实的频谱，这样就有可能发生一些频谱的峰点和谷点被“尖桩的栅栏”所挡住，不能被我们观察到。 实验中用到的函数：y=exp(x): 以e为底的指数。 conj(x): 取x的共轭，即改变x的虚部符号。 real(x): 取复数x的实部。 rand(1,N): 生成在0和1之间均匀分布的随机序列，长度为N。 randn(1,N): 生成正态分布(高斯分布)的随机序列，长度为N。 X=fft(x,N): 计算序列x的N点FFT。如果x的长度小于N，则在x后面补零；如果x的长度大于N，则对x进行截取；如果不指定参数N，则以x的实际长度作为FFT的点数。 x=ifft(X,N): 计算序列X的N点IFFT。 Y=fftshift(X): 将序列X分成左右两部分并交换位置。 实验内容 1、对高斯序列xa(n)，令p=8， q分别等于2、4和8，观察q值的改变对高斯序列时域特性和幅频特性的影响，并给出理论解释： 2、对衰减正弦序列xb(n)，a=0.1，f=0.0625，观察其时域和幅频特性，检查谱峰位置是否正确；改变f，使f分别等于0.4375和0.5625，比较这两种情况下频谱形状和谱峰位置，并给出理论解释： 3、观察三角波序列xc(n)和反三角波序列xd(n)的时域特性和幅频特性，用8点FFT分析幅频特性，观察两者的序列形状和频谱曲线有什么异同？为什么？在xc(n)和xd(n)的末尾补零，用32点FFT分析幅频特性，观察幅频特性发生了什么变化？这些变化说明了什么？ 4、对连续单频周期信号，按采样频率进行采样，截取长度分别选N=16和N＝20，比较其幅频特性，并给出理论解释；当截取长度选N＝164时，与N＝20比较，频谱发生了哪些变化？为什么？ 5、一个连续信号含两个频率分量，经采样得： 令N=16Δf分别为1/16和1/64，观察其频谱，两个谱峰是否可以区分开？为什么？当N=128时，Δf不变，结果有何不同？为什么？ 四、实验结果与分析 1.实验结果： 理论分析：采样序列的频谱是被采样信号频谱的周期延拓，当采样频率不满足奈奎斯特采样定理的时候，就会发生混叠。xa(n)中的参数p为高斯序列的峰值位置，q则表示高斯序列峰的尖锐度，q值越大，时域图中图象越平缓，序列变化越慢，其幅频特性图中高频分量越少，频谱越窄，在同样的采样频率下越不容易产生混叠。 2.实验结果： 结论：由图可知：a=0.1，f=0.0625时谱峰位置不正确，此时采样频率f太小，导致谱峰出现的位置不正确。频谱主瓣较宽，呈现出主瓣中间较为平缓，两侧较高的现象。 当f分别等于0.4375和0.5625时，函数的幅频特性完全相同，谱峰位置也完全相同，且没有产生明显的混叠和泄露，随采样频率的增加频谱主瓣越来越窄。 实验结果： 理论分析： 由图知：当N=8时，正反三角波的幅频特性相同，这是因为进行FFT运算时，要先对序列进行周期延拓，之后再取主值部分，而两个函数在周期延拓后是相同的，只差一个相位，因此得到的幅频特性图形也是相同的。 当N=32时，两个函数幅频特性不同。这是因为在补零后再，二者的周期延拓就不再相同。栅栏效应，当N=8时，一些谱线被挡住。通过在原序列的末端补零，增加采样的点数，增加了谱线的数量，弱化了栅栏效应。而且反三角波补零后在n=8处会有一个突变，因此它的幅频特性曲线中频谱分量会增多。 4．实验结果： 理论分析：①比较N=16与N=20时FFT的幅频特性，可以发现N=20时频谱出现泄漏，n=16时得到了单一谱线的频谱。这是由于时域中对信号的非整周期截断产生的频谱泄漏。N=20时截取了两个半周期，N