如何提取罗斯包中主题的频率(以Hz为单位)？

import numpy as np
import matplotlib.pyplot as plt

# 定义罗斯包方程
def rossler_attractor(x, y, z, a=0.2, b=0.2, c=5.7):
    dx = -y - z
    dy = x + a * y
    dz = b + z * (x - c)
    return dx, dy, dz

# 参数设置
dt = 0.01
num_steps = 100000
a, b, c = 0.2, 0.2, 5.7

# 初始化状态
x, y, z = 0.1, 0.1, 0.1
states = []

# 迭代生成数据
for _ in range(num_steps):
    dx, dy, dz = rossler_attractor(x, y, z, a, b, c)
    x += dx * dt
    y += dy * dt
    z += dz * dt
    states.append((x, y, z))

# 提取x分量进行傅里叶变换
x_values = np.array([state[0] for state in states])
fft_result = np.fft.fft(x_values)
freqs = np.fft.fftfreq(len(x_values), dt)

# 绘制频率谱
plt.figure(figsize=(10, 6))
plt.plot(freqs[:len(freqs)//2], np.abs(fft_result[:len(fft_result)//2]))
plt.xlabel('Frequency (Hz)')
plt.ylabel('Amplitude')
plt.title('Frequency Spectrum of Rossler Attractor')
plt.show()