由心跳信号求心率的python代码

以下是由心跳信号求心率的基本Python代码：

1. 加载数据

import scipy.io as sio

data = sio.loadmat('heartbeatsignal.mat') # 加载心跳信号数据
signal = data['signal'].reshape(-1) # 获取信号数据并转换为一维数组
fs = data['fs'][0][0] # 获取采样率

2. 滤波

from scipy.signal import butter, filtfilt

b, a = butter(2, [0.5, 30]/(fs/2), 'bandpass') # 2阶巴特沃斯带通滤波器设计
filtered_signal = filtfilt(b, a, signal) # 应用滤波器

3. 检测R峰

from biosppy.signals import ecg

out = ecg.ecg(signal=filtered_signal, sampling_rate=fs, show=False) # 使用biosppy库的ecg函数进行R峰检测
peaks = out['rpeaks'] # 获取R峰位置索引

4. 计算心率

rrInterval = np.diff(peaks) / fs # 计算相邻R峰的时间间隔
heartRate = 60 / rrInterval # 计算心率

注：此代码使用了`biosppy`库进行R峰检测，需要先通过`pip install biosppy`命令安装。同时，数据源（如示例代码中的`heartbeatsignal.mat`）需要根据实际情况进行替换。