基于Matlab的信号平稳性检验系统

基于Matlab的信号平稳性检验系统

3 数据仿真和分析
　　
　　将实验数据取为调频信号 x(t)=sin(sin(t／8)πt)。t起始为0，采样周期为0.1 s，数据长度为400个点，其时域波形和频域波形如图3所示。

　　
　　由式(6)产生的替代数据s(t)的时域波形和频域波形如图4所示。

　　
　　由图3和图4不难看出，替代数据与原数据的傅里叶变换幅值相同，但替代数据傅里叶变换的相位是随机的。

　　
　　图5显示了由式(7)计算得到的原数据和替代数据的时频分布图。由图5中可见，原数据的时频分布图有明显的结构性。它表明是非平稳的，而替代数据的时频分布图的结构性较原数据有明显减弱，表明替代数据的平稳性增加。
　　
　　由式(11)计算得到的?0的概率密度函数f(?0) 如图6所示，其中一共计算了1 000次替代数据。