如何计算时频域renyi熵？

Question

这就是任义熵在时频域的定义。 我将使用此函数来计算频谱图的 Renyi 熵。 px(t,f) 将是该频谱图的结果。 例如， array([[5.85959171e-01, 7.95627554e+00, 1.22937578e-01, ..., 1.82902033e+00, 5.96424909e-01, 3.20172883e+00], [8.12774889e+02, 8.03338062e+02, 6.11289787e+02, ..., 5.36101542e+02, 5.27223899e+01, 1.55881316e+02], [2.19951771e+01, 4.11390671e+02, 3.04105310e+02, ..., 1.01747040e+02, 2.36618805e+01, 8.56933907e+01], ..., [6.25524855e-04, 4.55459276e-04, 2.70485227e-04, ..., 1.40541362e-03, 1.38247253e-03, 4.60767521e-04], [1.16518569e-04, 1.56417512e-05, 6.71531675e-05, ..., 4.81702798e-04, 2.16510385e-04, 3.43761986e-04], [5.82912053e-05, 2.61144404e-05, 5.96264388e-05, ..., 1.23245212e-04, 3.72841550e-04, 1.46325938e-04]])那么，谁能帮我用这个数组编写这个函数的代码？

Answer 1

你可以简单地尝试以下（在 matlab 符号中）

log(sum(sum((tfd.^3))))其中tfd ：时频二维数组