我正在计算信号的PSD,我希望从0Hz到20Hz的频率范围内供电。这是我尝试使用linspace
df = pd.read_csv(path)
df = pd.DataFrame(df)
x = np.linspace(0, 20, 41)
dt = x[1] - x[0]
fs = 1 / dt
f,P = signal.welch(df, fs=5, nperseg=30, noverlap=5,axis=0)
在这里,我从0Hz到2.5Hz获得6个频率分量,但我想计算频率范围的功率,直到20Hz
任何人都可以帮我在这里做同样的事情。
PSD仅从0到fs / 2,您应该在fs参数中指定正确的采样频率,在您的情况下为1024。
此示例说明如何获取正弦信号的PSD:
import numpy as np
import matplotlib.pyplot as plt
from scipy.signal import welch
Fs = 1024 # Hz
Ts = 1/Fs
time = np.arange(0, 2, Ts) # 2 seconds
freqs = [20, 50, 100] # frequencies in Hz
x = np.zeros(len(time))
for f in freqs:
x += np.sin(2 * np.pi * f * time)
plt.plot(x)
f, P = welch(x, fs=Fs)
plt.figure()
plt.stem(f, P)