河海大学字信号处理实验三.docVIP

实验三 用FFT对信号进行频谱分析 班级： 姓名： 学号： 一、实验目的 1、加深对离散信号DFT的理解； 2、掌握FFT算法的流程及其MATLAB实现； 3、利用FFT对典型信号进行频谱分析； 4、结合理论知识，对频谱分析中出现的有关现象进行理论分析。 二、实验原理 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，程序： n=0:1:15; xa1=exp(-((n-8).^2)/2); subplot(3,2,1); plot(n,xa1,-*); xlabel(t/T); ylabel(x(n)); title(时域特性 p=8,q=2); xk1=abs(fft(xa1)); subplot(3,2,2); stem(n,xk1) xlabel(k); ylabel(X(k)); title(频域特性 p=8,q=2); xa2=exp(-((n-8).^2)/4); subplot(3,2,3); plot(n,xa2,-*); xlabel(t/T); ylabel(x(n)); title(时域特性 p=8,q=4); xk2=abs(fft(xa2)); subplot(3,2,4); stem(n,xk2) xlabel(k); ylabel(X(k)); title(频域特性 p=8,q=4); xa3=exp(-((n-8).^2)/8); subplot(3,2,5); plot(n,xa3,-*); xlabel(t/T); ylabel(x(n)); title(时域特性 p=8,q=8); xk3=abs(fft(xa3)); subplot(3,2,6); stem(n,xk3) xlabel(k); ylabel(X(k)); title(频域特性 p=8,q=8); 结果： 理论分析：采样序列的频谱是被采样信号频谱的周期延拓，当采样频率不满足奈奎斯特采样定理的时候，就会发生混叠。xa(n)中的参数p为高斯序列的峰值位置，q则表示高斯序列峰的尖锐度，q值越大，时域图中图象越平缓，序列变化越慢，其幅频特性图中高频分量越少，频谱越窄，在同样的采样频率下越不容易产生混叠。 2， 程序： n=0:1:15; xa1=exp(-0.1*n).*sin(2*pi*0.0625*n); subplot(3,2,1); plot(n,xa1,-*); xlabel(t/T); ylabel(x(n)); title(时域特性 a=-0.1,f=0.0625); xk1=abs(fft(xa1)); subplot(3,2,2); stem(n,xk1) xlabel(k); ylabel(X(k)); title(频域特性 a=-0.1,f=0.0625); xa2=exp(-0.1.*n).*sin(2*pi*0.4375.*n); subplot(3,2