matlab fftshift_Matlab数字信号处理基础三

问题1：对于正弦信号，请写一个Matlab程序验证对该正弦信号抽样时应遵循的几个准则要求：

1. 参数A、、 的数值自行设定；
2. 按准则①-④要求实现对该信号的抽样，比较抽样信号的频谱和x(t)的频谱；
3. 分别比较准则①-④其中一个不满足时抽样信号的频谱和波形。

关注我一起学习

（1）、由x(t)=Asin(2πf_0 t+φ)设A=1，f_0=100Hz，φ=0

得：x(t)=sin(2π*100t)；周期:T=1/f_0=0.01s；

（2）、按准则①-④要求实现对该信号的抽样的结果如下图所示：

由生成图像可知，此时的正弦信号满足四个条件，故做离散时间傅里叶变换的时候，频谱是在-100、100处的两条线谱，与原信号一致不会造成频谱泄露。

（3）、不满足条件一：生成抽样频率fs不满足正弦信号频率f0的整数倍的图像如下图：

当不满足条件一时，即抽样频率不为正弦信号频率的整数倍，fs=512Hz，f0=100Hz，

每个周期采fs/f0=5.12采样点N=64，采样点数为2的整次幂，包含4个完整的周期。该抽样信号的波形出现失真，但是频谱同样为在 f0处的两条线谱，无泄漏现象。

不满足条件二：生成截断后N的长度不包括完整的周期的信号的图像如下图所示：

当不满足条件二时，即截取长度不包含完整周期，该抽样信号的波形正确，但是出现频谱泄漏现象。

不满足条件三：生成不满足每个周期抽样为2的整次幂的信号的图像如图所示：

当不满足条件三时，即每一个周期抽样不满足2的整次幂。fs=300;f0=100;每一个周期采样点数为3个，截取长度为30,包含完整的10个周期。该抽样信号的波形出现严重失真，但是频谱同样为在 f0处的两条线谱，无泄漏现象。

不满足条件四：生成截断后补零的信号的图像如图所示：

当不满足条件四时，即截短后补零。该抽样信号的波形正确，但是出现频谱泄漏现象。

所以对该正弦信号x(t)=Asin(2πf_0 t+φ)做抽样时，只有同时满足四个条件，所得频谱才无泄漏，频谱为在 f0处的两条线谱。（实验所用MATLAB程序见附录一）。

问题2：设x(n)=sin⁡(2πf_1 n/f_s)+sin⁡(2πf_2 n/f_s)+sin⁡(2πf_3 n/f_s)，其中f_1=2Hz,f_2=2.02Hz,f_3=2.07Hz,f_s=10Hz。利用Matlab实现如下信号的分析：

1. 取该信号的长度为N=256、512和1024，分别画出该信号的幅度谱；
2. 对N=256和512两种情况，通过补零作1024点FFT并画出其幅度谱。
3. 在（1）和（2）中默认窗函数为矩形窗，现在改用汉明窗重复（1）和（2）的过程。
4. 根据以上结果说明DFT的物理频率分辨率和计算频率分辨率的区别。比较汉明窗和矩形窗在信号的频谱分析中的区别。

（1）截取信号长度N=256,512,1024的频谱图

当N1=256时，频率分辨率 =10/256≈0.039,因为f2-f1=0.02< ,f3-f1=0.07> ,所以仅能分辨出f1,f3的正弦分量，对应于频谱图的两个谱峰；当N2=512时，频率分辨率 =10/512≈0.0195,f2-f1=0.02> ,f3-f1=0.07> ，故能分辨f1，f2，f3的正弦分量，对应频谱图的三个谱峰。因为f2-f1与频率分辨率很接近，使f2的正弦分量的一部分没被分辨出来，对应的f2的谱峰较低当N3=1024时，频率分辨率 =10/1024≈0.0097,因为f2-f1=0.02> ,f3-f1=0.07> ，所以能分辨出f1，f2，f3的正弦分量，对应频谱图的三个谱峰。

（2）截取信号长度N=256,512补零作1024点FFT的频谱图

N=256个点补零至1024点后的信号频谱仍只有两个谱峰，物理分辨率没有改变，谱线间隔有所改变。

N=512个点补零至1024点后的信号频谱有3个谱峰，物理分辨率没有改变，谱线间隔有所改变。提高了频率的细分度，故f2正弦分量对应的谱峰幅度增大了。

（3）N=256,512,1024点分别加汉明窗后信号的频谱

当N=256点加汉明窗后，因为汉明窗主瓣加宽，旁瓣明显衰减，所以频率分辨率降低，只能分辨出f1的正弦分量对应频谱的谱峰也只有一个，并且频谱泄漏现象和频谱幅度都有所减少。

当N=512点加汉明窗后，有3个谱峰可以分辨出3个正弦分量，旁瓣、频谱泄漏现象和频谱幅度均有所减少。

当N=1024点加汉明窗后，同样有3个谱峰可以分辨出3个正弦分量，旁瓣、频谱泄漏现象和频谱幅度均有所减少。

N=256,512点补零并加汉明窗后信号的频谱

当N=256点补零至1024点后并加汉明窗的频谱仍然只有两个谱峰，正弦分量f1,f3可以分辨，由于汉明窗的旁瓣有显著衰减的特性，该频谱的泄漏现象明显减少，但频谱幅度较矩形窗要低一些。

当N=512点补零至1024点后并加汉明窗的频谱有三个谱峰， f1,f2,f3的正弦分量可以分辨，旁瓣明显减少，频谱的泄漏现象和频谱幅度也减少一些。

（4）由以上结果可知：

信号实际长度越长，物理分辨率越高，也就是物理分辨率与信号的实际长度成正比。但信号补零后不能增加数据的有效长度，所以补零不能提高物理分辨率，但能提高DFT的计算分辨率，减小栅栏效应。

汉明窗主瓣较宽，频谱幅度和频率分辨率较低，旁瓣明显衰减的特性；矩形窗主瓣窄，频谱幅度和频率分辨率较高，同时旁瓣也很大，会造成频谱泄漏。汉明窗适用于分析有较强干扰噪声的窄带信号，矩形窗适用于要求精确读出主瓣频率。

%1、%画出原函数时域与频域的图像