python如何提取一个复数信号的双谱特征，请给出代码

以下是Python代码示例，用于提取复数信号的双谱特征：

import numpy as np
from scipy.signal import stft

# 生成复数信号
t = np.linspace(0, 1, 1000)
x = np.exp(2j * np.pi * 10 * t) + np.exp(2j * np.pi * 20 * t)

# 计算短时傅里叶变换
f, t, Zxx = stft(x, nperseg=256)

# 计算双谱特征
S2 = np.abs(Zxx)**2
S3 = np.abs(Zxx)**3
S4 = np.abs(Zxx)**4
D = np.sqrt(S4 / (S2**2) - 1)

# 输出结果
print('双谱特征D：', D)

请注意，这只是一个示例代码，实际应用中需要根据具体情况进行调整。