自相关函数的频率估算matlab,科学网—3、自相关函数估计 - 柏世平的博文

对于做科研工作或工程应用的人来说，拿到一列或数列随机信号观测序列，首先想到的就是看看它们的自相关、互相关函数，自功率谱、互功率谱密度函数估计吧。

这一篇对上一篇中各数据矩阵TW4615_12、TW2672_12、GY2672_12、DY2672_12、MB2672_12数据作一下相关函数估计，并给出其图形。

MATLAB程序如下(注:为了以示区别,本博客各博文中运行的程序(包括源代码及注释)全部采用斜体)。

%先将各矩阵转换为零均值(时间平均)的一维时间序列TW461512_0=TW4615_12';4615天体温

TW461512_0=TW461512_0(:);

TW461512_0=TW461512_0-mean(TW461512_0);

subplot(2,2,1)%绘制子图

plot(TW461512_0)

TW267212_0=TW2672_12';%TW4615_12后半部分2672天体温

TW267212_0=TW267212_0(:);

TW267212_0=TW267212_0-mean(TW267212_0);

GY267212_0=GY2672_12';%高压

GY267212_0=GY267212_0(:);

GY267212_0=GY267212_0-mean(GY267212_0);

subplot(2,2,2)

plot(GY267212_0)

DY267212_0=DY2672_12';%低压

DY267212_0=DY267212_0(:);

DY267212_0=DY267212_0-mean(DY267212_0);

subplot(2,2,