心电信号提取

clc; clear all;
PATH= ‘D:/Program Files/work/‘;  % 指定数据的储存路径
HEADERFILE= ‘100.hea‘;      % .hea 格式，头文件，可用记事本打开
ATRFILE= ‘100.atr‘;         % .atr 格式，属性文件，数据格式为二进制数
DATAFILE=‘100.dat‘;         % .dat 格式，ECG 数据
SAMPLES2READ=2048;          % 指定需要读入的样本数
                            % 若.dat文件中存储有两个通道的信号:
                            % 则读入 2*SAMPLES2READ 个数据 
signalh= fullfile（PATH HEADERFILE）;    % 通过函数 fullfile 获得头文件的完整路径
fid1=fopen（signalh‘r‘）;    % 打开头文件，其标识符为 fid1 ，属性为‘r‘--“只读”
z= fgetl（fid1）;             % 读取头文件的第一行数据，字符串格式
A= sscanf（z ‘%*s %d %d %d‘[13]）; % 按照格式 ‘%*s %d %d %d‘ 转换数据并存入矩阵 A 中
nosig= A（1）;    % 信号通道数目
sfreq=A（2）;     % 数据采样频率
clear A;        % 清空矩阵 A ，准备获取下一行数据
         % 读取每个通道信号的数据信息
    z= fgetl（fid1）;
    A= sscanf（z ‘%*s %d %d %d %d %d‘[15]）;
    dformat（1）= A（1）;           % 信号格式; 这里只允许为 212 格式
    gain（1）= A（2）;              % 每 mV 包含的整数个数
    bitres（1）= A（3）;            % 采样精度（位分辨率）
    zerovalue（1）= A（4）;         % ECG 信号零点相应的整数值