精华内容
下载资源
问答
  • 2022-02-28 19:38:02

    信号根据 " 能量 " 可以分为 " 能量信号 "" 功率信号 " ;

    • 信号能量定义 : 整个轴上的能量先进行平方 , 然后求积分 ; 如果 能量 小于 无穷 , 则该信号 是 能量信号 ; 有限区间内的信号称为能量信号 ;

    • 信号功率定义 : 在一个信号周期内 , 进行积分求和操作 ; 如果 功率 小于 无穷 , 则该信号 是 功率信号 ; 周期信号 , 随机信号 是功率信号 ;

    本篇博客中的 互相关函数自相关函数 , 都是 " 能量信号 "相关函数 ;





    一、互相关函数



    互相关函数 表示的是 两个不同的信号 之间的相关性 ;

    x ( n ) x(n) x(n) y ( n ) y(n) y(n)" 互相关函数 " 如下 ,

    r x y ( m ) = ∑ n = − ∞ + ∞ x ∗ ( n ) y ( n + m ) r_{xy}(m) = \sum_{n=-\infty}^{+\infty} x^*(n) y(n + m) rxy(m)=n=+x(n)y(n+m)

    其中 y ( n ) y(n) y(n) 进行了移位 , 向左移动了 m m m 单位 ,

    该 " 互相关函数 " 求的是 y ( n ) y(n) y(n) 移位 m m m 后的序列 x ( n ) x(n) x(n) 序列之间的关系 ;

    注意这里的 n n n 表示的是时刻 , m m m 表示的是信号移动的间隔 ;


    " 互相关函数 " 表示的是 x ( n ) x(n) x(n) 信号 , 与 隔了 m m m 时间后的 y ( n ) y(n) y(n) 信号之间的关系 ;

    2 2 2 个信号 ( 序列 ) 之间 " 关系 " 是一个 函数 , 函数的自变量是 m m m 间隔 , 不是 n n n ;





    二、自相关函数



    自相关函数 ( Autocorrelation Function ) :

    r x x ( m ) = ∑ n = − ∞ + ∞ x ∗ ( n ) x ( n + m ) = r x ( m ) r_{xx}(m) = \sum_{n=-\infty}^{+\infty} x^*(n) x(n + m) = r_x(m) rxx(m)=n=+x(n)x(n+m)=rx(m)

    " 自相关函数 " 是 " 自己信号 " 与 " 隔一段时间后的 自己信号 " 之间的 相关性 ;


    如果 m = 0 m = 0 m=0 时 , " 自己信号 " 与 " 隔一段时间 m m m 后的自己信号 " 完全相等 , 该值就是 信号的能量 ;

    r x ( 0 ) = ∑ n = − ∞ + ∞ ∣ x ( n ) ∣ 2 = E r_{x}(0) = \sum_{n=-\infty}^{+\infty} |x(n)|^2= E rx(0)=n=+x(n)2=E



    更多相关内容
  • 自相关函数及互相关函数编程
  • 自相关函数及互相关函数编程
  • 此 m 文件查找并绘制生成的长度为 31 的 Gold 代码的自相关和互相关函数(ACF 和 CCF)。互相关是三个值。
  • 本文提出了一种归一化互相关函数基音检测算法(NCCFPDA——Normalized Cross-Correlation Function Pitch Detection Algorithm),该算法在基音检测主过程的前端和末端加入了有效的预处理和后处理技术,实验分析 表明...
  • 能够产生良好的m序列的自相关以及互相关结果图,适合初学者
  • xcorrTD 的结果已针对 MatLAB 的 xcorr 函数进行了验证。 有关频域中的互相关,请参阅 xcorrFD。 语法:[lags,ck,cc,td] = xcorrTD(x,y) 输入: x = 输入信号 1(必须是 Nx1 或 1xN 向量) y = 输入信号 2(必须是 ...
  • 自相关函数和互相关函数的matlab 计算和作图 1. 首先说说自相关和互相关的概念 这个是信号分析里的概念他们分别表示的是两个时间序列之间和同一个时 间序列在任意两个不同时刻的取值之间的相关程度即互相关函数是...
  • 理想规则图像特征经标准互相关函数匹配后,相似函数C*(x)可用确定的解析式Z(x)表达。但是对于含噪声图像,相似函数C(x)较之C*(x)发生了变化,但存在“零相似不变性”。将此原理应与于显微视觉中图像物面分辨率的在线...
  • 时间序列分析可以定义为在给定先前值的情况下预测随机过程的未来值。 建模的一个重要部分是决定应该使用多少先行值来预测未来。 自相关函数显示两个系列之间的相关系数... 对于两个过程,还添加了互相关和部分互相关
  • 互相关函数

    2014-05-03 00:39:26
    互相关函数
  • matlab互相关函数

    2014-10-13 11:50:25
    程序是matlab写的互相关仿真实例,可以正常运行
  • m序列少值互相关函数一直都是研究者感兴趣的方向之一,但这方面取得的成果并不完善。有限域上多元高次方程的求解成为解决该问题的关键。对于采样因子形式为d=(pl+1)/(pk+1)的m序列互相关函数,目前已有研究大多是针对...
  • 一种基于FPGA的互相关函数加速器设计.pdf
  • 由于MATLAB本身自带的相关函数在扩频通信中并不适合,性能欠佳。本程序是我自己编写的求自相关或互相关的MATLAB函数。可直接调用该函数。已通过验证。
  • 互相关函数的计算

    千次阅读 2021-10-19 13:30:14
    目录互相关函数的定义互相关函数的计算存在具体项目参数时 整理不易,有条件的点个关注、点个赞呗!感恩各位大哥! 互相关函数的定义 描述两个不同的信号在不同时期上的相关性的函数,主要应用: 混有周期成分数据...

    整理不易,有条件的点个关注、点个赞呗!感恩各位大哥!

    互相关函数的定义

    描述两个不同的信号在不同时期上的相关性的函数,主要应用:

    • 混有周期成分数据(信号)的频率(周期)提取,例如两列数据在其中一列数据滞后三期时相关性最高,则该类数据的周期为3。

    互相关函数的计算

    公式:
    C 1 , 2 ( t ) ≈ ∫ 0 t c v 1 ( τ ) v 2 ( t + τ ) d τ (1) C_{1,2}(t)\approx\int_{0}^{t_{c}}v_{1}(\tau)v_{2}(t+\tau)d\tau\tag{1} C1,2(t)0tcv1(τ)v2(t+τ)dτ(1)
    实际计算举例,本实验中主要针对采用傅里叶变换的循环互相关,因此后续以循环互相关为例进行讲解!
    循环互相关的思想如下:
    在这里插入图片描述
    注意到在计算时要用到超出原始数据索引范围的数据,其数据补充方式并不是“补零”而是“周期延拓”。
    循环互相关实际则通过傅里叶变换进行计算,将上述两列数据分别变换到频率域,从而得到频率域同等长度的两列数据 [ 10. + 0. j − 2.5 + 3.4409548 j − 2.5 + 0.81229924 j ] , [ 10. + 0. j − 2.5 + 3.4409548 j − 2.5 + 0.81229924 j ] [10. +0.j -2.5+3.4409548j -2.5+0.81229924j],[10. +0.j -2.5+3.4409548j -2.5+0.81229924j] [10.+0.j2.5+3.4409548j2.5+0.81229924j][10.+0.j2.5+3.4409548j2.5+0.81229924j],则输出互相关矩阵为第一个元素为:
    ( 10 ∗ 10 + 0 j ∗ 0 j ) + ( 10 ∗ 0 j − 0 j ∗ 10 ) = ( 100 + 0 j ) (10*10+0j*0j)+(10*0j-0j*10)=(100+0j) (1010+0j0j)+(100j0j10)=(100+0j),同理第二个元素为 ( − 2.5 ∗ − 2.5 + 3.4409548 j ∗ 3.4409548 j ) + ( − 2.5 ∗ 3.4409548 j + 2.5 ∗ 3.4409548 j ) (-2.5*-2.5+3.4409548j *3.4409548j )+(-2.5*3.4409548j+2.5*3.4409548j) 2.52.5+3.4409548j3.4409548j+2.53.4409548j+2.53.4409548j)

    存在具体项目参数时

    计算方式一:

                for c in range(len(irover)):
                    # ================= Non center Station Index =================
                    Atr=st[irover[c]]
                    # ================= Location =================
                    # ================= Do FFT at non center =================
                    freq_rover, time_rover, Sxx_rover_ = sg.spectrogram(Atr.data  ,
                                                                        nperseg=wd, # 窗口大小,对应fftlen
                                                                        noverlap=0, # noverlap对应overlaprate,但是noverlap是具体值,overlaprate是重叠率
                                                                        nfft=nft, # 实际就是窗口大小,表示是否对数据进行填充。
                                                                        fs=Fs, # 频率,暂时不清楚频率在计算中的作用
                                                                        scaling='spectrum', # 输出设置为功率谱
                                                                        mode='complex')
                    
                    # ======== Find Data in range output frequency ===========
                    Sxx_rover_=Sxx_rover_[idx]
                    
                    # ================= Measure AutoCorrelation Ratio ================
                    rat_autocorr=np.zeros((len(Sxx_pusat_[:,0]),len(Sxx_pusat_[0,:])), dtype=np.complex)
                    for nwd in range(len(Sxx_pusat_[0,:])):
                        Abs_Sxx_pusat=(abs(Sxx_pusat_[:,nwd]))
                        Smooth_Abs_Sxx_pusat=sg.convolve(Abs_Sxx_pusat,smooth,mode='same')
                        
                        Abs_Sxx_rover=(abs(Sxx_rover_[:,nwd]))
                        Smooth_Abs_Sxx_rover=sg.convolve(Abs_Sxx_rover,smooth,mode='same')
                        
                        rat_autocorr[:,nwd]=((Sxx_pusat_[:,nwd]*np.conj(Sxx_rover_[:,nwd]))/
                                    (Smooth_Abs_Sxx_pusat*Smooth_Abs_Sxx_rover))
                        
                    time_autocorr=np.mean(rat_autocorr,axis=1)
                    autocorr[:,c]=time_autocorr   
    

    计算方式二:

    for(long i=0;i<d.nsta;i++){
                if(((k*d.steplen)>=d.startend[i*2])&&((k*d.steplen+d.fftlen)<d.startend[i*2+1]))
                {
                    ifCC[i] = 1;
                }
                else
                {
                    ifCC[i] = 0;   
                }
                
                if (ifCC[i]){
                    for(long j=0;j<d.fftlen;j++){
                        in[i][j] = d.data[j+k*d.steplen+d.npts*i];
                    }
                    fftwf_execute(p[i]);
                    if(d.ifspecwhitenning){
                        for (long j=0;j<d.nf;j++){
                            absv = sqrtf(out[i][j*d.fstride][0]*out[i][j*d.fstride][0]+out[i][j*d.fstride][1]*out[i][j*d.fstride][1]);
                            absv = max(absv,1e-8f);
                            if(absv!=0){
                                datar[j+d.nf*i] = out[i][j*d.fstride][0]/absv;
                                datai[j+d.nf*i] = out[i][j*d.fstride][1]/absv;
                            }
                        }
                    }
                    else{
                        for(long j=0;j<d.nf;j++){
                            datar[j+d.nf*i] = out[i][j*d.fstride][0];
                            datai[j+d.nf*i] = out[i][j*d.fstride][1];
                        }
                    }
                }
            }
    #ifdef useOMP
    #pragma omp parallel
    #pragma omp for
    #endif
            for(long i=0;i<npairs; i++){
                if(ifCC[d.Pairs[i*2]]&ifCC[d.Pairs[i*2+1]]){
                    for(long j=0;j<d.nf;j++){
                        d.ncfsr[j+i*d.nf] += datar[j+d.Pairs[i*2]*d.nf]*datar[j+d.Pairs[i*2+1]*d.nf] + datai[j+d.Pairs[i*2]*d.nf]*datai[j+d.Pairs[i*2+1]*d.nf];
                        d.ncfsi[j+i*d.nf] += datar[j+d.Pairs[i*2]*d.nf]*datai[j+d.Pairs[i*2+1]*d.nf] - datai[j+d.Pairs[i*2]*d.nf]*datar[j+d.Pairs[i*2+1]*d.nf];
                    }
                    CCnumbers[i] = CCnumbers[i] + 1;
                }
            }
        }
    

    相同之处:均采用傅里叶变换求解
    不同之处:方法二窗口的指定是在傅里叶变换函数包的外面,且其将数据集切分为每步再传入傅里叶变换函数。方法一在傅里叶变换中指定数据窗口大小和每次移动步长,暂时不知道这两种方式是否具有差异。

    展开全文
  • 使用matlab验证时域和频域计算互相关函数
  • 自相关函数与互相关函数

    万次阅读 多人点赞 2021-01-16 23:49:31
    1 概念 相关函数是描述信号X(s),Y(t)(这两个信号可以是随机的,也可以是确定的)在任意两个不同时刻s、t... 相关函数分为自相关和互相关。 自相关函数是描述随机信号 x(t) 在任意不同时刻 t1,t2 的取值之间的相...

    目录

    1 概念

    2 自相关函数

    2.1 定义

    2.2 性质

    3 互相关(cross-correlation)函数

    3.1 定义

    3.2 性质

    3.3 线性互相关(linear cross-correlation)

    3.4 循环互相关(Circular Cross-Correlation)的定义和计算

    3.5 用线性互相关处理周期性信号

    3.6 相关问题QA

    3.7 参考资料


    1 概念

          相关函数是描述信号X(s)Y(t)(这两个信号可以是随机的,也可以是确定的)在任意两个不同时刻s、t的取值之间的相关程度。两个信号之间的相似性大小用相关系数来衡量。定义:

            

            称为变量 X 和 Y 的相关系数。若相关系数 = 0,则称 X与Y 不相关。相关系数越大,相关性越大,但肯定小于或者等于1.。

            

            相关函数分为自相关互相关

    • 自相关函数是描述随机信号 x(t) 在任意不同时刻 t1,t2 的取值之间的相关程度。
    • 互相关函数是描述随机信号 x(t)、y(t) 在任意两个不同时刻s,t的取值之间的相关程度

    自相关函数

    2.1 定义

            自相关函数是描述随机信号 x(t) 在任意不同时刻 t1,t2的取值之间的相关程度。自相关函数,是对信号自身的互相关, 表示同一序列不同时刻的相关程度。是用寻找重复模式的数字工具,就如一个存在被覆盖噪声的周期信号,或识别丢失的基频。它经常被用于信号处理中的分析函数或序列,如时域信号 。定义式:

            

    或者:

            

    2.2 性质

    主要性质如下:

    自相关系数:

    3 互相关(cross-correlation)函数

    3.1 定义

            自相关是互相关的一种特殊情况.。互相关函数是描述随机信号 x(t)、y(t) 在任意两个不同时刻 s,t 的取值之间的相关程度,其定义为: 

               

            对于连续函数,有定义:

          

           对于离散的,有定义:

          

            从定义式中可以看到,互相关函数卷积运算类似,也是两个序列滑动相乘,但是区别在于:互相关的两个序列都不翻转,直接滑动相乘,求和;卷积的其中一个序列需要先翻转,然后滑动相乘,求和。所以,f(t)g(t) 做相关等于 f*(-t)g(t) 做卷积。

            在图象处理中,自相关和互相关函数的定义如下:设原函数是 f(t),则自相关函数定义为 R(u)=f(t)*f(-t),其中*表示卷积;设两个函数分别是f(t)和g(t),则互相关函数定义为R(u)=f(t)*g(-t),它反映的是两个函数在不同的相对位置上互相匹配的程度。

    3.2 性质

    互相关函数的性质:

    互相关系数:

            正如卷积有线性卷积(linear convolution)循环卷积(circular convolution)之分;互相关也有线性互相关(linear cross-correlation)循环互相关(circular cross-correlation)。线性互相关和循环互相关的基本公式是一致的,不同之处在于如何处理边界数据。其本质的不同在于它们对原始数据的看法不同。

    3.3 线性互相关(linear cross-correlation)

            假设我们手里有两组数据,分别为个和个,表示为:长,即。序列之间的线性互相关操作表示为,其结果也是一个序列,表示为。具体的操作是用这两个序列进行的一种类似“滑动点积”的操作,如图1和图2所示。

    图1. 线性互相关的计算过程示意

    图2. 线性互相关结果序列中单个值计算示意

            得到的互相关序列总长度为,该序列的前和后个数值是无效的,有效的数据共个。线性互相关的有效数据第个分量的值为:

    注意,线性互相关并不满足交换律,即:

            

            一个简单的应证是,等式两侧操作所得结果的有效数据个数都不一致。

            线性相关的实际意义是,向量中的各个与向量等长的子向量与向量的相似程度。这样,中值最大的索引就是与向量中与最相似的子向量的起始索引。通常,为了获得有效的互相关数据,我们总是用较短的数据去滑动点积较长的数据。

            用一个实际的应用例子来验证一下吧。如图3的第一个子图表示雷达声纳发射了一个探测信号。经过一段时间之后,收到了如图3的第二个子图所示的回波(带有一定的噪声)。此时我们关注的是如何确定回波中从何时开始是对探测信号的响应,以便计算目标距雷达的距离,这就需要用到线性互相关。在第三个子图中的‘Valid’曲线即是有效互相关数据,其中清晰地呈现出两处与探测信号相似的回波的位置。

    图3. 相关计算的一个例子:雷达回波分析

            线性互相关中,还有一些概念值得注意:

    • 一是补零。由线性相关的计算式不难发现,为了计算出个完整的相关系数序列(包含那些“无效数据”在内的所有结果),需要用到一些“不存在”的点。这就需要人为地对这些值进行补充,在线性相关的计算中,对这些超出原始数据储存的区域取值为零。
    • 二是末端效应。由图1可以发现,一头一尾的个互相关数据并没有完全“嵌入”两个原始数组的全部信息,它们或多或少地受到了人为补零的影响。因此一般认为这些数据是不可用的。
    • 三是计算模式的选择。这个问题其实是由问题二衍生而来的,就Python语言中的函数而言,至少有两个可以直接计算线性相关:

    numpy.correlate(a, v, mode)scipy.signal.correlate(a, v, mode)

    它们的调用参数完全相同。在调用时有三种模式可供选择,它们计算的内容是相同的,但是返回值长度各不相同:
    mode = ‘valid’:只返回有效的那一部分相关数据,共$M-N+1$个;
    mode = ‘same’:只返回与 等长的那一部分相关数据,共$N$个;
    mode = ‘full’:返回全部相关数据,共$M+N-1$个。
    图3的第三个子图展示了这三种模式的计算结果,在那个例子中,‘valid’模式是最合适的。

    3.4 循环互相关(Circular Cross-Correlation)的定义和计算

            循环互相关是表征两组等长周期性数据之间相似性的操作,其与线性互相关的区别也正由“等长”“周期性”这个两特点产生。在循环互相关中,被处理的原始数据是等长的,即。序列之间的线性互相关操作表示为,其结果也是一个序列,表示为。其计算式与线性互相关的写法是一致的:

            

            只是得到的互相关序列长度也为。循环互相关的计算的具体过程如图4所示,注意到在计算时要用到超出原始数据索引范围的数据,其数据补充方式并不是“补零”而是“周期延拓”:即。这意味着对于循环互相关,不存在不同的计算模式之分,所有的数据都是有效数据

    图4. 循环互相关的计算过程示意

            注意,循环互相关也不满足交换律

            这里给出了一个关于循环相关的算例。两路原始数据分别由如下函数生成:

            

            

            如果视为某个线性系统的周期输入信号,而视为这个线性系统的输出信号。由于存在外接干扰,因此输出信号不完全由输入信号决定。此时,循环互相关的实际意义是,分辨输出信号中的哪一个部分(频率成分)是由该输入信号产生的。

     

    图5. 时域数据,从上到下:和他们的循环互相关

    图6. 频谱,从上到下:和他们的循环互相关

            从图5和图6可以看出,循环互相关的频谱准确地说明了那些测试信号的相关性。

            遗憾的是,在Python几大数值计算库中,并没有直接可计算循环相关的函数。但是可以采用如下代码构造出一个可用的(经过归一化的)cxcorr(a, v)函数出来:

    def cxcorr(a,v):
    nom = np.linalg.norm(a[:])*np.linalg.norm(v[:])
    return fftpack.irfft(fftpack.rfft(a)*fftpack.rfft(v[::-1]))/nom

            图4中的数据就是通过这个函数计算出来的。其中用到了傅里叶变换反变换来计算循环互相关,这是可行的。它们之间的关系在第四小节的QA中专门讨论。

    3.5 用线性互相关处理周期性信号

            实际上,线性相关也可以处理周期信号,前提是将两组信号采样成长度差异较大的序列。这样,其有效线性互相关也可以完美地反应数据之间的相关性。

            同样采用第二节中的例子。这时为了保证足够的有效线性互相关数据,两组数据的长度故意不一致(但都足够表征其特征),如图7所示。它们的频谱如图8所示,仍然完美地体现了测试数据的相关性。

    图7. 时域数据,从上到下:和他们的线性互相关

    图8. 频谱,从上到下:和他们的线性互相关

            既然线性互相关也能处理周期性数据,为什么还要专门搞一个基于等长序列和周期延拓的循环互相关呢?实际上,正如后文QA中专门讨论的,这是为了利用快速傅利叶变换加速计算。

    3.6 相关问题QA

            至此,两种常用的互相关评价方法及其计算已经总结完毕。然而其中还有一些细节尚待分辨。例如,序列之间的互相关的计算式:

            

            与卷积(convolution)的定义式:

            

    如此类似,如果再联想起傅里叶变换的卷积定理,那么,至少会产生如下的问题:

    Q.1:它们之间有更深意义上的联系吗?
    A.1:文献[1]的答复是坚决的:“不要让求卷积和互相关的数学相似性迷惑你,它们描述了不同的信号处理过程。卷积是系统输入信号、输出信号和冲激响应之间的关系互相关是一种在噪声背景下检测已知信号的方法。二者在数学上的相似仅仅是一种巧合。”实际上,只要注意到卷积操作是满足交换律的,而互相关操作并不满足交换律。仅此一点也许就能说明它们有着本质的不同吧。

    Q.2:可以利用Python中计算卷积的函数来计算互相关吗?
    A.2可以,但是只能用以计算线性互相关。Python中的numpy.convolve()函数就可以计算两个序列之间的卷积。在卷积的计算过程中也会自动进行补零(而不是周期延拓,这就是为什么只能计算线性相关的原因),这种卷积有时被称为线性卷积,同样涉及末端效应、有效数据长度等考虑。具体地,根据相关和卷积的表达式,如果希望计算序列之间的线性互相关序列。等效地,只需要计算序列之间的卷积。表示序列的“反置”,即将序列[1,2,3]反置为[3,2,1]。

    Q.3:可以根据傅立叶变换的性质中有卷积定理,利用傅立叶正/逆变换计算互相关吗?
    A.3可以,但是只能用于计算循环互相关。傅立叶变换的卷积定理中所涉及的卷积是循环卷积。与前述的线性卷积是不同的。实际上不同的并不是卷积本身,它们的计算式是一致的,而是在如何看待参与卷积计算的数据,线性卷积认为参与计算的序列之外都是零,而循环卷积认为参与计算的序列是一个无限循环的数据的一段——这导致了它们对“越界”数据的补齐方式不一样。正如线性互相关和循环互相关的区别!先将循环互相关等效为一个循环卷积,再利用快速傅里叶变换计算卷积即可。实际上本文给出的cxcorr(a, v)函数正是利用这一性质来计算循环相关的。其对计算速度的提升是相当明显的。

    Q.4:怎样进行归一化(normalization),以便于比较互相关数据?
    A.4:根据参考[4],用公式:

    3.7 参考资料

    [1] Steven W. Smith. Digital Signal Processing: A Practical Guide for Engineering and Scientists [M].
    张瑞峰, 詹敏晶 等译. 实用数字信号处理,从原理到应用[M]. 人民邮电出版社, 北京, 2010.
    [2] Mark Owen. Practical Signal Processing [M].
    丘天爽, 李丽, 赵林 译. 实用信号处理 [M]. 电子工业出版社, 北京, 2009.
    [3] 关于MATLAB中的xcorr() 的论述
    http://www.mathworks.cn/cn/help/signal/ref/xcorr.html
    [4] 关于MATLAB中的cxcorr() 的论述
    http://www.mathworks.com/matlabcentral/fileexchange/4810-circular-cross-correlation
    [5] 网络论坛Stackoverflow关于此问题的讨论
    http://stackoverflow.com/questions/6991471/computing-cross-correlation-function
    http://stackoverflow.com/questions/12323959/fast-cross-correlation-method-in-python
    http://stackoverflow.com/questions/9281102/n-fold-fft-convolution-and-circular-overlap
    http://stackoverflow.com/questions/6855169/convolution-computations-in-numpy-scipy
    http://stackoverflow.com/questions/4688715/find-time-shift-between-two-similar-waveforms
    [6] 关于Cross-correlation的定义
    http://mathworld.wolfram.com/Cross-Correlation.html
    http://paulbourke.net/miscellaneous/correlate/
    http://en.wikipedia.org/wiki/Cross-correlation
    [7] 关于 Circular Cross-correlation的定义
    http://en.wikipedia.org/wiki/Circular_convolution
    http://cnx.org/content/m22974/latest/

     

     

     

    自相关函数与互相关函数:http://www.doc88.com/p-5129647069822.html

    互相关(cross-correlation)及其在Python中的实现:https://blog.csdn.net/icameling/article/details/85238412

    第五章:自相关: https://thinkdsp-cn.readthedocs.io/zh_CN/latest/05-autocorrelation.html

    展开全文
  • 一、相关函数最大值、 1、自相关函数最大值、 2、互相关函数最大值、 二、能量有限信号的相关函数在 m 趋近无穷时为 0





    一、相关函数最大值




    1、自相关函数最大值


    自相关函数 在 自变量 m = 0 m = 0 m=0 , 永远大于其它 m ≠ 0 m \not= 0 m=0 的值 ;

    r x ( 0 ) ≥ r x ( m ) r_x(0) \geq r_x(m) rx(0)rx(m)

    也就是说 , 自相关函数 的 最大值 , 就是 m = 0 m = 0 m=0 时的值 ;


    2、互相关函数最大值


    互相关函数最大值是 r x ( 0 ) r y ( 0 ) \sqrt{r_x(0)r_y(0)} rx(0)ry(0) ,

    r x ( 0 ) r_x(0) rx(0) x ( n ) x(n) x(n) 信号的 能量 ;

    r y ( 0 ) r_y(0) ry(0) y ( n ) y(n) y(n) 信号的 能量 ;

    ∣ r x y ( m ) ∣ ≤ r x ( 0 ) r y ( 0 ) = E x E y |r_{xy}(m)| \leq \sqrt{r_x(0)r_y(0)} = \sqrt{E_xE_y} rxy(m)rx(0)ry(0) =ExEy





    二、能量有限信号的相关函数在 m 趋近无穷时为 0



    如果 信号 x ( n ) x(n) x(n)信号 y ( n ) y(n) y(n) 都是 能量信号 ,

    能量信号 指的是 能量有限 的 信号 , 能量是 绝对可和 的 ,

    与之 相对的 是 功率信号 , 功率信号 能量无限 ,

    能量信号 一定 不是 周期信号 ,


    lim ⁡ m → ∞ r x ( m ) = 0 \lim\limits_{m \rightarrow \infty} r_x(m) = 0 mlimrx(m)=0

    lim ⁡ m → ∞ r x y ( m ) = 0 \lim\limits_{m \rightarrow \infty} r_{xy}(m) = 0 mlimrxy(m)=0


    物理意义 : m = 0 m = 0 m=0 , 这两个序列的 相关性最大 , 但是 随着 m m m 增加到 无穷大 ∞ \infty , 则相关性直接变为 0 0 0 ,

    有限序列 , 一旦平移 , 总有 错开的时候 , 一旦错开 , 就任何相关性也没有了 , 相关性为 0 0 0 ;

    展开全文
  • 这个是信号分析里的概念,他们分别表示的是两个时间序列之间和同一个时间序列在任意两个不同时刻的取值之间的相关程度,即互相关函数是描述随机信号x(t),y(t)在任意两个不同时刻t1,t2的取值之间的相关程度,自相关...
  • 信号相关分析原理自相关函数,互相关函数 * 5.1 信号的互能量与互能谱 (一)信号的能量与功率 信号的能量 指信号f(t)的归一化能量即信号的电 压电流加在1?电阻上所消耗的能量 5.11 若f(t)为实数 由公式 当R=1?时即可得...
  • 随机信号的自相关函数和互相关函数Matlab
  • 一、相关函数共轭对称性质、 1、实信号自相关函数偶对称、 2、复信号自相关函数共轭对称、 3、复信号互相关函数共轭对称
  • 自相关函数&互相关函数

    千次阅读 2019-08-15 23:22:45
    这个是信号分析里边的概念。...实函数和,如为能量有限信号,它们之间互相关函数定义为:(注:下角数字,前面的领先) 互相关函数是两信号时间差的函数。 一般 如果和是同一信号,可记为,这...
  • PAGE / NUMPAGES 自相关函数和互相关函数的matlab计算和作图 1. 首先说说自相关和互相关的概念 这个是信号分析里的概念他们分别表示的是两个时间序列之间和同一个时间序列在任意两个不同时刻的取值之间的相关程度即...
  • 互相关函数、互功率谱和卷积之间的关系 1.我们要实现怎样的目标? 如果有两个复信号, 连续信号表示为y1(t)y_1(t)y1​(t)和 y2(t)y_2(t)y2​(t); 离散信号表示为y1(n)y_1(n)y1​(n)和y2(n)y_2(n)y2​(n); 两个信号...
  • 一、功率信号的互相关函数、 二、功率信号的自相关函数
  • 自相关、互相关函数学习笔记

    千次阅读 多人点赞 2019-07-17 15:56:50
    为什么相关性有效 参考文献 【1】自相关的物理意义 【2】自相关和互相关函数计算方法总结及心得体会
  • 互相关函数python实现的三种方法

    千次阅读 2022-01-20 18:53:47
    用三种不同方法计算信号的互相关函数,并得到了几乎一致的结果。为了更好的理解整个过程,可在CSDN免费下载源码

空空如也

空空如也

1 2 3 4 5 ... 20
收藏数 56,859
精华内容 22,743
关键字:

互相关函数