def hjorth_cal(x, axis=-1):
#输入1维或N维数据
#返回activity、Mlbility、complexity
x = np.asarray(x)
# Calculate derivatives
dx = np.diff(x, axis=axis)
ddx = np.diff(dx, axis=axis)
# 方差计算
x_var = np.var(x, axis=axis)
dx_var = np.var(dx, axis=axis)
ddx_var = np.var(ddx, axis=axis)
act =x_var
# Mobility and complexity
mob = np.sqrt(dx_var / x_var)
com = np.sqrt(ddx_var / dx_var) / mob
return act, mob, com