weltch法画功率谱密度

要绘制功率谱密度图，可以使用Python中的scipy库中的signal模块。下面是一个简单的例子：

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

# 生成一个信号
t = np.linspace(0, 10, 1000, endpoint=False)
x = np.sin(2*np.pi*3*t) + np.sin(2*np.pi*5*t)

# 计算功率谱密度
f, Pxx = signal.welch(x)

# 绘制功率谱密度图
plt.semilogy(f, Pxx)
plt.xlabel('frequency [Hz]')
plt.ylabel('PSD [V**2/Hz]')
plt.show()

这段代码会生成一个包含两个正弦波的信号，并计算该信号的功率谱密度。然后，使用matplotlib库绘制功率谱密度图。其中，semilogy函数将y轴设置为对数刻度，使得图像更易于观察。