MATLAB心电图信号分析实例

在MATLAB中进行心电信号（ECG）分析是一个常见的应用，它可以帮助我们理解心脏活动并检测异常情况。以下是简单的心电图信号分析的一个例子：

1. **导入信号**：首先使用`load`或`readmat`函数从存储设备或文件中读取心电信号数据。数据通常包含两个通道：一个表示电压变化（R波），另一个表示时间。

data = load('ecg_data.mat');
signal = data.signal; % R波电压数据
timestamp = data.time; % 时间戳

2. **预处理**：去除噪声，例如50Hz电源干扰，使用高通滤波器或Butterworth滤波器。

fs = 1000; % 根据实际情况设定采样频率
[~, filtered_signal] = butter(2, [0.5 30]/fs, 'bandpass', fs);
filtered_signal = filter(filtered_signal, 1, signal);

3. **R-R间隔分析**：寻找每个R波后的下一个R波的时间间隔（即心跳周期），可以使用`diff`函数。

rr_intervals = diff(find(signal > threshold)); % 利用阈值找到R波位置

4. **心率计算**：根据R-R间隔计算平均心率。

mean_hr = mean(rr_intervals); % 平均心率

5. **图形展示**：使用`plot`绘制原始信号和处理过的信号，以及心率分布图。

figure;
subplot(2,1,1);
plot(timestamp, signal, timestamp, filtered_signal);
legend('Original Signal', 'Filtered Signal');

subplot(2,1,2);
histogram(rr_intervals, 'Normalization', 'probability');
xlabel('R-R Interval (ms)');
ylabel('Probability');

**相关问题--:**
1. MATLAB中如何识别QRS复杂波形？
2. 对于异常心律，如何使用MATLAB进行诊断？
3. 有没有现成的MATLAB函数库可以直接处理心电图数据？