matlab 计算心电信号周期并作图

计算心电信号周期可以通过以下步骤实现：

1. 导入心电信号数据
2. 对信号进行预处理，如滤波、去基线漂移等
3. 对预处理后的信号进行QRS波群检测，提取QRS峰点位置
4. 计算相邻QRS峰点之间的时间差，即为心电信号的周期

以下是一个简单的示例代码：

% 导入心电信号数据
ecg = load('ecg_data.mat').ecg_data;

% 滤波
fs = 1000; % 采样率为1000Hz
fc = 50; % 滤波截止频率为50Hz
[b, a] = butter(2, fc/(fs/2), 'high');
ecg_filt = filtfilt(b, a, ecg);

% QRS波群检测
[qrs_amp_raw, qrs_i_raw] = pan_tompkin(ecg_filt, fs, 0);

% 计算QRS峰点之间的时间差，即为心电信号的周期
period = diff(qrs_i_raw) / fs;

% 作图
plot(period);
xlabel('QRS峰点序号');
ylabel('周期（秒）');
title('心电信号周期图');

其中，`pan_tompkin`函数是一个QRS波群检测算法，这里使用了Matlab自带的`ecg-kit`工具箱中的实现。如果没有安装该工具箱，可以自行寻找其他QRS波群检测算法的实现。