实验报告-噪声中非平稳信号的频谱分析.docVIP

实验报告三 姓名：贾少波 学号：S311050041 导师：张友文 实验名称：噪声中非平稳信号的频谱分析 实验目的： 学习和掌握短时傅里叶变换的基本概念。 仿真实现短时傅里叶变换对噪声中的线性调频信号的处理。 分析短时傅里叶变换的性能，同时分析影响性能的因素及原因。 实验原理 对于非平稳信号的处理，存在一个二维函数，它既能反映该信号的频率内容，也能反映出该频率内容随时间变化的规律。这一问题的信号处理理论称为信号的联合时频分布。其中最重要的是以Cohen类为代表的双线性时频分布，此分布可表示为： ............................ 1.1 式中是一个二维的窗函数，给定不同的窗函数可得到不同的时频分布。在上式中x t 出现了两次，且是相乘的形式，这一特点称为双线性。 若 1，则上式可简化为Wigner-Ville分布，此分布可表示为 .............................................................. 1.2 若，式中w是一个一维的窗函数，则 1.1 式可简化为如下的谱图（spectrogram）： ........................................... 1.3 式中...................................................... 1.4 称为信号x t 的短时傅里叶变换（STFT）。以上四式都是的二维函数，称它们为时频联合分布。前三个反映了信号能量随时间和频率的分布，第四个反映了信号的频谱随时间和频率的分布。 四、主要编程步骤 周期图法步骤： 1．编制一个非平稳正弦信号，选择三个频率不同的信号首尾相连并且和噪声做和。本实验选择的是80Hz，90Hz，210Hz,根据采样定理选择采样频率为1000Hz。 2．采用经典的功率谱估计方法对此信号进行分析，本实验采用的周期图法，20段welch法和现代谱估计处理方法。 3．对于非平稳信号，可以用双线性时频分布描述信号的频率随时间变化规律， 试验中，选择Specgram函数估计出频率幅度，取模求平方，转化为dB格式，选择hanning窗，重叠64点，nfft为128。 4．改变不同重叠点数，观察实验结果变化。 5．改变不同窗函数，观察实验结果。 五、实验结果及分析 1．设计的含有噪声和非平稳信号序列如下图 2. 采用周期图法：Fs选择1000Hz，nfft 1048，N选择1020点, 功率谱分析结果如下图 采用welch法 采用LD递推法估计 采用自相关法 结果分析：周期图法可以估计出三个频率的功率谱，但是不能确定在什么时间段存在什么频率成分 3．使用短时傅里叶变换进行时频分析，如下图 结果分析：短时傅里叶变化可以看出功率谱随时间和频率的变化，分辨出三段时间依次有三个频率成分。 4.改变不同的重叠点数，观察功率谱在时域上投影判定分辨率性质 重叠0点 重叠250点 结果分析：从上图可以看出，随着重叠点数增大，分辨率降低，上表面更平滑，这是因为重叠点数增加相当于数据加窗，所以导致分辨率降低。 六、结论 短时傅里变换可以对非平稳信号进行时频分析，解决了FFT不能同时进行时域和频域分析的缺陷在一定程度上克服了标准傅立叶变换不具有局部分析能力的缺陷，窗口函数的形状一旦确定便无法改变，其时频分辨率也无法改变，只能实现一定程度的时频局部化。这不符合实际问题中高频信号的时域分辨率应比低频信号高的要求。要提高频域分辨率就得加大时窗长度，同时造成时域分辨率下降。