信号的相关

news/2024/11/28 17:39:56/

信号的相关

  • 从向量的内积谈起
  • 定义
    • 互相关
    • 互相关的性质
    • 自相关
    • 自相关的性质
    • 相关与卷积
  • 相关的归一化形式
  • 功率和周期信号的相关

  有时候需要将一个信号与多个参考信号进行比较,以便确定每对信号之间的相似性,从而从相似性中提取出额外的信息。

  在雷达应用中,接收到的信号是目标反射回来的信号,该信号是反射信号的延迟形式,通过测量延迟,可以确定目标的位置。

从向量的内积谈起

  按照一般教科书的写法,一般会直接给出一个互相关的公式,比如两个能量信号序列的相关性的度量可以定义为
r x y [ l ] = ∑ n = − ∞ ∞ x [ n ] y [ n − l ] r_{xy}[l]=\sum_{n=-\infty}^{\infty}x[n]y[n-l] rxy[l]=n=x[n]y[nl]先不解释这个公式的组成以及含义,但是我学习这个的时候一直想,为什么这个公式能够给出两个信号之间的相关性的度量。

  所以我打算从向量的内积开始,从向量的内积讲到这个公式,不过由于接下来的都是我的理解,我在数字信号处理的教材中还没有看到相关的介绍(可能是我读书少),所以接下来的数学推导并不是很准确(事实上涉及的数学推导几乎没有),重点是这种思想的理解。

  首先考虑两个向量,这两个向量的模长都是固定的,那么怎么描述两个向量之间的相关程度,在我看来,二者之间的夹角大小可以描述这种相关性,如果两个向量相关性高,说明两向量之间的夹角比较小,根据向量的内积
a ⃗ ⋅ b ⃗ = ∣ a ∣ ∣ b ∣ c o s ϕ \vec{a} \cdot \vec{b}=\vert a \vert \vert b \vert cos\phi a b =abcosϕ其中 ϕ \phi ϕ为向量 a ⃗ \vec{a} a 与向量 b ⃗ \vec{b} b 之间的夹角,如果两个向量相关性越高,那么二者之间的夹角越小,所以 c o s ϕ cos\phi cosϕ越大,由于向量的模长固定,所以相关性越高也就代表着两向量的内积越大。

  如果用坐标形式表示,假设两个向量都是二维向量,那么向量的内积可以表示为
a ⃗ ⋅ b ⃗ = x a x b + y a y b \vec{a} \cdot \vec{b}=x_ax_b+y_ay_b a b =xaxb+yayb

虽然我这里用的二维向量,但是明显可以推广到无限维的情况。

  那么接下来我们在来看之前在文章开头定义的那个公式
r x y [ l ] = ∑ n = − ∞ ∞ x [ n ] y [ n − l ] r_{xy}[l]=\sum_{n=-\infty}^{\infty}x[n]y[n-l] rxy[l]=n=x[n]y[nl]对于序列 x [ n ] x[n] x[n]和序列 y [ n − l ] y[n-l] y[nl]可以看做是一个无限维的向量,而
x [ − ∞ ] , . . . , x [ 0 ] , . . . , x [ ∞ ] x[-\infty], ..., x[0], ..., x[\infty] x[],...,x[0],...,x[]可以看做是"向量" x [ n ] x[n] x[n](这里我把 x [ n ] x[n] x[n]看做是一个向量了)的坐标,相应的
y [ − ∞ ] , . . . , y [ − l ] , . . . y [ ∞ ] y[-\infty], ..., y[-l], ... y[\infty] y[],...,y[l],...y[]可以看做"向量" y [ n − l ] y[n-l] y[nl]的坐标,所以
∑ n = − ∞ ∞ x [ n ] y [ n − l ] \sum_{n=-\infty}^{\infty}x[n]y[n-l] n=x[n]y[nl]可以看做是"向量" x [ n ] x[n] x[n]与"向量" y [ n − l ] y[n-l] y[nl]的内积,而对于取不同的 l l l x [ n ] x[n] x[n]的模与 y [ n − l ] y[n-l] y[nl]的模都是不变的(不论序列 y [ n ] y[n] y[n]是有限长还是无限长),符合我们推导时做的假设,所以在这里我从向量的内积这个方面解释了为什么这个公式可以表示两个信号的相关性程度。

定义


互相关

  上面从向量内积的角度理解了互相关公式,接下来将按照教科书上给出严格的定义。

  互相关序列 r x y [ l ] r_{xy}[l] rxy[l]表示一对能量序列 x [ n ] x[n] x[n] y [ n ] y[n] y[n]之间相似性的度量,其定义为(假设式中的无限求和是收敛的)
r x y [ l ] = ∑ n = − ∞ ∞ x [ n ] y [ n − l ] r_{xy}[l]=\sum_{n=-\infty}^{\infty}x[n]y[n-l] rxy[l]=n=x[n]y[nl]
式中, r x y [ l ] r_{xy}[l] rxy[l]代表了序列 x [ n ] x[n] x[n]与序列 y [ n − l ] y[n-l] y[nl]之间的相关性程度,所以 r x y [ l ] r_{xy}[l] rxy[l]越大,代表序列 x [ n ] x[n] x[n]与序列 y [ n − l ] y[n-l] y[nl]的相似性越高。

互相关的性质

  由互相关公式 r x y [ l ] = ∑ n = − ∞ ∞ x [ n ] y [ n − l ] r_{xy}[l]=\sum_{n=-\infty}^{\infty}x[n]y[n-l] rxy[l]=n=x[n]y[nl]
r y x [ l ] = ∑ n = − ∞ ∞ y [ n ] x [ n − l ] r_{yx}[l]=\sum_{n=-\infty}^{\infty}y[n]x[n-l] ryx[l]=n=y[n]x[nl]
m = n − l m=n-l m=nl,得到
r y x [ l ] = ∑ m = − ∞ ∞ x [ m ] y [ m − ( − l ) ] = r x y [ − l ] \color{red}r_{yx}[l]=\sum_{m=-\infty}^{\infty}x[m]y[m-(-l)]=r_{xy}[-l] ryx[l]=m=x[m]y[m(l)]=rxy[l]

自相关

  同理,我们定义自相关序列,只要将 y [ n ] y[n] y[n]替换为 x [ n ] x[n] x[n]即可
r x x [ l ] = ∑ n = − ∞ ∞ x [ n ] x [ n − l ] r_{xx}[l]=\sum_{n=-\infty}^{\infty}x[n]x[n-l] rxx[l]=n=x[n]x[nl]
由于上面推导得到
r y x [ l ] = r x y [ − l ] r_{yx}[l]=r_{xy}[-l] ryx[l]=rxy[l]
所以得到
r x x [ l ] = r x x [ − l ] r_{xx}[l]=r_{xx}[-l] rxx[l]=rxx[l]
即自相关序列是一个偶序列。

自相关的性质

  另一个需要注意的是,当 l = 0 l=0 l=0时,自相关序列
r x x [ 0 ] = ∑ n = − ∞ ∞ x [ n ] x [ n ] = ∑ n = − ∞ ∞ ∣ x [ n ] ∣ 2 r_{xx}[0]=\sum_{n=-\infty}^{\infty}x[n]x[n]=\sum_{n=-\infty}^{\infty}\vert x[n] \vert^2 rxx[0]=n=x[n]x[n]=n=x[n]2
所以 r x x [ 0 ] r_{xx}[0] rxx[0]的大小就是序列 x [ n ] x[n] x[n]的能量。

  根据我们之前在向量的内积中的定义 ∑ n = − ∞ ∞ x [ n ] x [ n ] \sum_{n=-\infty}^{\infty}x[n]x[n] n=x[n]x[n]可以看做是序列 x [ n ] x[n] x[n]与序列 x [ n ] x[n] x[n]的内积,二者之间的夹角为 0 0 0(为最小夹角,相关性最大),所以
∑ n = − ∞ ∞ x [ n ] x [ n ] ≥ ∑ n = − ∞ ∞ x [ n ] x [ n − l ] \sum_{n=-\infty}^{\infty}x[n]x[n] \geq\sum_{n=-\infty}^{\infty}x[n]x[n-l] n=x[n]x[n]n=x[n]x[nl]
r x x [ l ] ≤ r x x [ 0 ] , − ∞ &lt; l &lt; ∞ r_{xx}[l]\leq r_{xx}[0],\quad -\infty &lt; l &lt; \infty rxx[l]rxx[0],<l<


  在这个的基础上,我们可以简单的谈一下在文章开头提到的雷达检测问题:

在雷达应用中,接收到的信号是目标反射回来的信号,该信号是反射信号的延迟形式,通过测量延迟,可以确定目标的位置。

  这里有个重点是接收到的信号是发射信号的延时,我们记接收到的信号为 h [ n ] h[n] h[n],发射的信号为 x [ n ] x[n] x[n],如果接收到的信号是延迟 2 s 2s 2s的话,那么 h [ n ] = x [ n − 2 ] h[n]=x[n-2] h[n]=x[n2],假设检测的过程是,计算 h [ n ] h[n] h[n] x [ n ] x[n] x[n]及其延迟的互相关值,得到序列 r [ l ] r[l] r[l]
r [ l ] = ∑ n = − ∞ ∞ h [ n ] x [ n − l ] r[l]=\sum_{n=-\infty}^{\infty}h[n]x[n-l] r[l]=n=h[n]x[nl]由于 h [ n ] = x [ n − 2 ] h[n]=x[n-2] h[n]=x[n2],所以序列 r [ 2 ] r[2] r[2]是序列 r [ l ] r[l] r[l]的最大值,所以我们通过观察序列 r [ l ] r[l] r[l]最大值出现的位置,就可以知道延迟的时间,从而确定目标的位置。

  实际上接收到的信号还附有噪声,分析的过程比这个要复杂一些,但是基本的原理还是这个,只是牵涉到了一些数学公式,我觉得把握这种思想就可以了。

相关与卷积

  继续探究互相关序列
r x y [ l ] = ∑ n = − ∞ ∞ x [ n ] y [ n − l ] r_{xy}[l]=\sum_{n=-\infty}^{\infty}x[n]y[n-l] rxy[l]=n=x[n]y[nl]
将其适当的变形可以写成
r x y [ l ] = ∑ n = − ∞ ∞ x [ n ] y [ − ( l − n ) ] = x [ l ] ∗ y [ − l ] r_{xy}[l]=\sum_{n=-\infty}^{\infty}x[n]y[-(l-n)]=x[l]*y[-l] rxy[l]=n=x[n]y[(ln)]=x[l]y[l]
所以互相关序列可以写成序列 x [ l ] x[l] x[l]与序列 y [ − l ] y[-l] y[l]的卷积。

  同理,自相关序列可以写成
r x x [ l ] = x [ l ] ∗ x [ − l ] r_{xx}[l]=x[l]*x[-l] rxx[l]=x[l]x[l]

相关的归一化形式

  在之前我们提到用向量的内积表示相关性的大小,其本质是由其夹角决定的,所以我们可以直接利用夹角的余弦值来定义相关性的大小。在向量理论中:
c o s ϕ = a ⃗ b ⃗ ∣ a ⃗ ∣ ∣ b ⃗ ∣ cos\phi = \frac{\vec{a} \vec{b}}{\vert \vec{a} \vert \vert \vec{b} \vert} cosϕ=a b a b

注意到
∑ n = − ∞ ∞ ∣ x [ n − l ] ∣ 2 = ∑ n = − ∞ ∞ ∣ x [ n ] ∣ 2 = r x x [ 0 ] \sum_{n=-\infty}^{\infty}\vert x[n-l] \vert ^2 = \sum_{n=-\infty}^{\infty}\vert x[n] \vert ^2 = r_{xx}[0] n=x[nl]2=n=x[n]2=rxx[0]
而这个也可以看做事向量模长的平方。

  所以我们定义自相关的归一化形式为
ρ x x [ l ] = r x x [ l ] r x x [ 0 ] r x x [ 0 ] = r x x [ l ] r x x [ 0 ] \rho_{xx}[l]=\frac{r_{xx}[l]}{\sqrt{r_{xx}[0]}\sqrt{r_{xx}[0]}}=\frac{r_{xx}[l]}{r_{xx}[0]} ρxx[l]=rxx[0] rxx[0] rxx[l]=rxx[0]rxx[l]

  同理我们定义互相关的归一化形式为
ρ x y [ l ] = r x y [ l ] r x x [ 0 ] r y y [ 0 ] \rho_{xy}[l]=\frac{r_{xy}[l]}{\sqrt{r_{xx}[0]}\sqrt{r_{yy}[0]}} ρxy[l]=rxx[0] ryy[0] rxy[l]

  所以 ∣ ρ x x [ l ] ∣ ≤ 1 , ∣ ρ x y [ l ] ≤ 1 ∣ \vert \rho_{xx}[l] \vert \leq 1,\vert \rho_{xy}[l]\leq 1 \vert ρxx[l]1,ρxy[l]1

ρ x y [ l ] = 1 \rho_{xy}[l]=1 ρxy[l]=1代表 c o s ϕ = 1 ⇒ ϕ = 0 cos\phi=1 \Rightarrow \phi=0 cosϕ=1ϕ=0,说明这时序列的相关性最大
ρ x y [ l ] = 0 \rho_{xy}[l]=0 ρxy[l]=0代表 c o s ϕ = 0 ⇒ ϕ = π / 2 cos\phi=0 \Rightarrow \phi=\pi /2 cosϕ=0ϕ=π/2,说明这两个序列不相关
ρ x y [ l ] = − 1 \rho_{xy}[l]=-1 ρxy[l]=1代表 c o s ϕ = − 1 ⇒ ϕ = π cos\phi=-1 \Rightarrow \phi=\pi cosϕ=1ϕ=π,说明这时序列线性负相关

功率和周期信号的相关

  在之前的介绍中,介绍的都是能量信号的相关,对于功率信号来说,它的能量为无穷大,所以上述的定义不太适合,类似于功率信号的定义,我们定义功率信号 x [ n ] , y [ n ] x[n],y[n] x[n],y[n]的互相关为
r x y [ l ] = lim ⁡ K → ∞ 1 2 K + 1 ∑ n = − K K x [ n ] y [ n − l ] r_{xy}[l]=\lim\limits_{K \to \infty}\frac{1}{2K+1}\sum_{n=-K}^{K}x[n]y[n-l] rxy[l]=Klim2K+11n=KKx[n]y[nl]功率信号 x [ n ] x[n] x[n]的自相关定义为
r x x [ l ] = lim ⁡ K → ∞ 1 2 K + 1 ∑ n = − K K x [ n ] x [ n − l ] r_{xx}[l]=\lim\limits_{K \to \infty}\frac{1}{2K+1}\sum_{n=-K}^{K}x[n]x[n-l] rxx[l]=Klim2K+11n=KKx[n]x[nl]

  如果 x ~ [ n ] \tilde{x}[n] x~[n] y ~ [ n ] \tilde{y}[n] y~[n]是周期为 N N N的两个周期信号,那么其互相关序列为
r x ~ y ~ [ l ] = 1 N ∑ n = 0 N − 1 x ~ [ n ] y ~ [ n − l ] r_{\tilde{x} \tilde{y}}[l]=\frac{1}{N}\sum_{n=0}^{N-1}\tilde{x}[n]\tilde{y}[n-l] rx~y~[l]=N1n=0N1x~[n]y~[nl]
x ~ [ n ] \tilde{x}[n] x~[n]的自相关序列为
r x ~ x ~ [ l ] = 1 N ∑ n = 0 N − 1 x ~ [ n ] x ~ [ n − l ] r_{\tilde{x} \tilde{x}}[l]=\frac{1}{N}\sum_{n=0}^{N-1}\tilde{x}[n]\tilde{x}[n-l] rx~x~[l]=N1n=0N1x~[n]x~[nl]

  我们很容易的知道,序列 r x ~ y ~ [ l ] r_{\tilde{x} \tilde{y}}[l] rx~y~[l]和序列 r x ~ x ~ [ l ] r_{\tilde{x} \tilde{x}}[l] rx~x~[l]都是周期序列,且周期为 N N N


http://www.ppmy.cn/news/694938.html

相关文章

HDMI上电信号分析

重点是HPD这里检测信号不能错。这样才能读取EDID&#xff0c;TDMS。 要保证TMDS时钟信号完整&#xff0c;关于TDMS信号的阻抗不会影响显示。我因为TMDS0和TMDS2通道反了 HDMI TYPE A接口包含&#xff1a; 3个TMDS数据通道 &#xff0c;1个TMDS时钟通道&#xff0c;CEC控制信号…

小白入门脑电信号

是个脑电信号研究小白&#xff0c;今年3月下旬临时换了研究方向&#xff0c;到现在也看了不少论文&#xff0c;试过一些特征提取的代码还有分类的代码。现在也想总结一下这段时间所学的东西&#xff0c;算是做个中期回顾吧。 Ps:里面的特征提取的代码以及分类的代码均不是我所…

信号

1. 定义 简而言之&#xff0c;信号是一种软件中断&#xff0c;提供了一种处理异步的方法&#xff0c;信号发生是随机的。例如键盘输入中断按键(^C)&#xff0c;它的发生在程序执行过程中是不可预测的。 硬件异常也能产生信号&#xff0c;例如被零除、无效内存引用&am…

心电信号的特征提取、分析与处理

心电信号的特征提取、分析与处理* 数据来源&#xff1a;MIT-BIH数据库&#xff08;可从以下数据中任选两组进行实验&#xff09; 给出4组不同病例的心电信号数据&#xff0c;分别命名为“100-2-3”&#xff0c;“105-2-3”&#xff0c;“109-2-3”&#xff0c;“111-2-3”&…

脑电信号(EEG)简介:原理、应用与分析方法

目录 一、EEG基本原理 二、EEG的应用领域 1.临床应用 2.认知科学 3. 神经反馈治疗 三、EEG信号分析方法 1.预处理 2.时域分析 3.频域分析 4.时频分析 5.机器学习与深度学习 四、挑战与展望 1.挑战 2.展望 总结 参考文献 一、EEG基本原理 脑电信号&#xff0c;也称脑电图&#…

【肌电信号】基于matlab GUI肌电信号处理【含Matlab源码 966期】

⛄一、简介 肌电信号又称为EMG&#xff1b;EMG&#xff08;electromyography&#xff09;是一种生物电流信号&#xff0c;它产生于任意组织器官&#xff0c;一般是一种时间和一系列的振幅、频率和波形函数 。肌电信号是伴随肌肉收缩动作产生的一种生物电信号,采集皮肤表面的肌…

脑电信号情绪识别研究综述

摘要 情绪识别是指通过人的面部表情、行为动作或者生理信号等信息识别人的情绪状态,其成果在医疗辅助、教育、交通安全等方面有很大的应用价值。由于脑电信号的客观真实性等特点,使用脑电信号进行情绪识别研究受到国内外学者们的广泛关注。查阅了大量脑电情绪识别相关文献并进行…

【脑电信号】基于matlab小波变换脑电信号特征提取【含Matlab源码 511期】

⛄一、小波变换简介 小波变换是从20世纪80年代起逐渐发展成熟的一项数学应用技术&#xff0c;具有对时间—频率的双重分析和多分辨率分析能力&#xff0c;目前已经广泛应用于图像处理、模式识别等多个领域。小波变换的窗口大小固定但形状可改变&#xff0c;因此能够满足时域—…