不同采样率的道互相关函数

//* 不同采样率的道互相关函数   
02.void DisDtCorrelativity(float &_af, float *_TracePick, float *_Trace, int _Len,   
03.                        int _Dt, int _Index)   
04.//** _af        :互相关函数   
05.//** _TracePick : trace one 错动道   
06.//** _Trace     : 相关参考道   
07.//** _Index     : 错动样点个数   
08.//** _Len       : 相关段长度   
09.//** _Dt        : 采样率差别倍数   
10.{   
11.    double sum11 = 0.0;   
12.    double sum22 = 0.0;   
13.    double sum1  = 0.0;   
14.    double sum2  = 0.0;   
15.    double sum12 = 0.0;   
16.   
17.    for (int i=0; i<_Len; i++)   
18.    {   
19.        int k = _Dt * i + _Index;   
20.        sum1  = sum1 + _Trace[i];   
21.        sum2  = sum2 + _TracePick[k];   
22.        sum12 = sum12 + _Trace[i] * _TracePick[k];   
23.        sum11 = sum11 + _Trace[i] * _Trace[i];   
24.        sum22 = sum22 + _TracePick[k] * _TracePick[k];   
25.    }   
26.    double aa = (sum11*_Len-sum1*sum1) * (sum22*_Len-sum2*sum2);   
27.    if (aa <= 0)   
28.    {   
29.        _af = 0.0;   
30.        return;   
31.    }   
32.    double a = sqrt(aa);   
33.    double b = (sum12*_Len - sum1*sum2);   
34.    _af = b/a;   
35.} 

转载于:https://my.oschina.net/u/582827/blog/1499074

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值