功率谱 python

根据提供的引用内容，可以使用Yule-Walker方法来估计平稳随机信号的功率谱。在Python中，可以使用scipy库中的signal模块来实现Yule-Walker方法。具体步骤如下：

1. 导入需要的库和模块：

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

2. 生成一个平稳随机信号：

n = 1024
x = np.random.randn(n)

3. 使用Yule-Walker方法估计功率谱：

order = 30
a, e, k = signal.arburg(x, order)
w, h = signal.freqz(1, a, worN=n, whole=True)
psd = (e / np.abs(1 - np.exp(-1j*w)))**2

4. 绘制功率谱图：

plt.plot(w, psd)
plt.xlabel('Frequency')
plt.ylabel('Power Spectral Density')
plt.show()

以上代码将生成一个平稳随机信号的功率谱图。