Python,计算数据中心的滞后

时间:2014-11-22 09:34:01

标签: python numpy signals distribution lag

我的两个数据集是:

fnamerp1=([  93,   87,   96,   93,   90,  123,  111,   82,   87,  115,  103,
        101,   93,   92,  111,  107,  114,  106,  116,  106,  128,  115,
        141,  134,  120,  149,  140,  166,  152,  171,  192,  207,  227,
        266,  270,  286,  355,  385,  397,  488,  462,  531,  579,  622,
        711,  720,  801,  858,  906,  915,  915,  956, 1004, 1012, 1045,
        1076, 1063, 1013,  985,  924,  959,  838,  766,  763,  742,  642,
        587,  557,  484,  393,  353,  341,  284,  240,  221,  209,  147,
        109,  113,  102,   71,   63,   63,   50,   29,   39,   36,   25,
        30,   23,   27,   23,   19,   19,   24,   15,   23,   21,   26,
        15])

fnamerp2=([  105,  89, 120, 121, 103, 105, 113,  94, 104, 115, 122, 116, 121,
   129, 118, 126, 138, 146, 161, 163, 178, 192, 194, 222, 268, 272,
   285, 342, 380, 378, 373, 448, 493, 511, 571, 603, 691, 772, 738,
   796, 839, 832, 883, 930, 963, 975, 972, 931, 947, 941, 934, 964,
   871, 869, 826, 793, 733, 708, 606, 610, 515, 483, 409, 352, 358,
   264, 266, 205, 191, 167, 136, 138,  99, 102,  82,  57,  65,  53,
    51,  32,  26,  27,  39,  21,  29,  23,  25,  24,  16,  17,  27,
    33,  19,  13,  24,  26,  18,  22,  18,  20])

我想找到两个峰的中心之间的滞后(不仅仅是它们的最大值)。我的计划是使用np.argmax(signal.correlate(fnamerp1,fnamerp2))

从数学的角度来看,从Python中优雅地做到这一点的正确方法是什么?

0 个答案:

没有答案