使用FIR滤波器实现对心电信号特征提取的程序，使用MATLAB软件

下面是使用MATLAB软件实现对心电信号特征提取的程序，其中使用了FIR滤波器。

% 读取心电信号数据
data = load('ecg_data.mat');
ecg_signal = data.ecg_signal;

% 设计FIR滤波器
fs = 200; % 采样频率
f1 = 5;   % 截止频率1
f2 = 50;  % 截止频率2
n = 100;  % 滤波器阶数
h = fir1(n, [f1/(fs/2), f2/(fs/2)]);

% 对心电信号进行滤波
filtered_ecg_signal = filter(h, 1, ecg_signal);

% 计算心电信号的特征
R_peaks = ecg_r_peak_detection(filtered_ecg_signal, fs);
Q_peaks = ecg_q_peak_detection(filtered_ecg_signal, R_peaks);
S_peaks = ecg_s_peak_detection(filtered_ecg_signal, R_peaks);
T_peaks = ecg_t_peak_detection(filtered_ecg_signal, R_peaks);

% 绘制心电信号及其特征点
t = (1:length(filtered_ecg_signal))/fs;
plot(t, filtered_ecg_signal);
hold on;
plot(t(R_peaks), filtered_ecg_signal(R_peaks), 'ro');
plot(t(Q_peaks), filtered_ecg_signal(Q_peaks), 'gx');
plot(t(S_peaks), filtered_ecg_signal(S_peaks), 'bx');
plot(t(T_peaks), filtered_ecg_signal(T_peaks), 'kx');
legend('心电信号', 'R峰', 'Q峰', 'S峰', 'T峰');
xlabel('时间 (s)');
ylabel('电压 (mV)');

其中，`ecg_r_peak_detection`、`ecg_q_peak_detection`、`ecg_s_peak_detection`和`ecg_t_peak_detection`分别是检测心电信号中R峰、Q峰、S峰和T峰的函数。您可以根据自己的需求来实现这些函数。