血压fft算法代码

以下是基于Python实现的血压FFT算法代码：

import numpy as np
import matplotlib.pyplot as plt

# 读取数据
data = np.loadtxt("blood_pressure_data.txt")

# 采样频率
fs = 100  # Hz

# FFT变换
fft_data = np.fft.fft(data)
freq = np.fft.fftfreq(len(data), 1.0/fs)

# 选取主频率范围
low_freq = 0.5  # Hz
high_freq = 4.0  # Hz
idx = np.where((freq>=low_freq) & (freq<=high_freq))

# 计算主频率下的幅值
amp = 2.0*np.abs(fft_data[idx])/len(data)

# 绘制频率谱图
plt.plot(freq[idx], amp)
plt.xlabel('Frequency (Hz)')
plt.ylabel('Amplitude')
plt.title('Blood Pressure FFT')
plt.show()

需要注意的是，需要将血压数据保存在名为"blood_pressure_data.txt"的文本文件中，每行为一个数据点。此外，需要安装NumPy和Matplotlib库。