还在按部就班的算自相关?FFT让你体验飞一般的感觉! |
您所在的位置:网站首页 › matlab求互相关 › 还在按部就班的算自相关?FFT让你体验飞一般的感觉! |
本文同步发布在我的个人博客,希望大家可以到我的个人博客玩耍。http://www.weekreport.cn/archives/212 自相关函数自相关(Autocorrelation),也叫序列相关,是一个信号于其自身在不同时间点的互相关。非正式地来说,它就是两次观察之间的相似度对它们之间的时间差的函数。它是找出重复模式(如被噪声掩盖的周期信号),或识别隐含在信号谐波频率中消失的基频的数学工具。它常用于信号处理中,用来分析函数或一系列值,如时域信号。–百度百科 自相关函数在随机信号处理领域非常重要,一般用下式计算一个随机信号的自相关。
R
(
s
,
t
)
=
E
[
X
(
s
)
x
(
t
)
]
R(s,t) = E[X(s)x(t)]\quad
R(s,t)=E[X(s)x(t)] 如果该随机信号是均方遍历的,则可以使用一个足够长样本函数的时间自相关去估计真实的自相关函数。 对时间自相关函数两边求傅立叶变换并整理得下式: 随机信号 u ( n ) u\left ( n \right ) u(n)的观测样本由一段零均值、方差为 σ 2 \sigma^{2} σ2的复高斯白噪声序列 v ( n ) v\left ( n \right ) v(n)叠加3个归一化频率分别为 f 1 = 0.15 , f 2 = 0.17 , f 3 = 0.26 f_{1}=0.15,f_{2}=0.17,f_{3}=0.26 f1=0.15,f2=0.17,f3=0.26,信噪比分别为 30 d B , 30 d B , 27 d B 30dB,30dB,27dB 30dB,30dB,27dB的复正弦信号构成。 仿真结果令信号观测样本长度
N
=
32
N=32
N=32,基于
F
F
T
FFT
FFT的自相关函数快速算法估计出的自相关函数
r
0
^
(
m
)
\hat{r_0}\left( m \right)
r0^(m),并与由信号观测样本直接计算得到的自相关函数进行比较,结果如下,
r
1
r_1
r1代表基于
F
F
T
FFT
FFT的自相关函数快速算法估计出的自相关函数,
r
0
r_0
r0代表由信号观测样本直接计算得到的自相关函数,可以看到,二者结果是一样的。 利用
M
a
t
l
a
b
Matlab
Matlab中的
t
i
c
,
t
o
c
tic,toc
tic,toc函数,得到两种方法的运行时间,可以看到基于
F
F
T
FFT
FFT的自相关函数快速算法要比由信号观测样本直接计算自相关函数快一个数量级左右,这是很大的提升,特别是在样本点数
N
N
N很大时。 |
今日新闻 |
推荐新闻 |
CopyRight 2018-2019 办公设备维修网 版权所有 豫ICP备15022753号-3 |