Python 等角度采样阶次分析代码示例

以下是一个简单的 Python 代码示例，用于执行等角度采样的阶次分析：

import numpy as np
import scipy.signal as signal
import matplotlib.pyplot as plt

# 生成测试数据
t = np.linspace(0, 1, 1000)
x = np.sin(2 * np.pi * 10 * t) + np.sin(2 * np.pi * 20 * t) + np.sin(2 * np.pi * 30 * t)

# 等角度采样
n = len(x)
k = np.arange(n)
theta = 2 * np.pi * k / n
X = np.fft.fft(x) / n

# 计算阶次谱
order_spectrum = np.abs(X) ** 2
order_spectrum = np.fft.fftshift(order_spectrum)

# 绘制阶次谱图像
plt.plot(theta, order_spectrum)
plt.xlabel('Theta (radians)')
plt.ylabel('Order Spectrum')
plt.show()

在这个示例中，我们首先生成了一个包含多个正弦波的测试信号。然后，我们对信号进行等角度采样，计算其阶次谱，并绘制了阶次谱的图像。该图像显示了每个旋转频率的能量分布，其中峰值表示信号中的阶次分量。