rddata.m -- MIT-BIH ECG 信号的数据读取Matlab程序.m

% This programm reads ECG data which are saved in format 212.
% （e.g. 100.dat from MIT-BIH-DB cu01.dat from CU-DB...）
% The data are displayed in a figure together with the annotations.
% The annotations are saved in the vector ANNOT the corresponding
% times （in seconds） are saved in the vector ATRTIME.
% The annotations are saved as numbers the meaning of the numbers can
% be found in the codetable “ecgcodes.h“ available at www.physionet.org.
%
% ANNOT only contains the most important information which is displayed
% with the program rdann （available on www.physionet.org） in the 3rd row.
% The 4th to 6th row are not saved in ANNOT.
%
%
%      created on Feb. 27 2003 by
%      Robert Tratnig （Vorarlberg University of Applied Sciences）
%      （email: rtratnig@gmx.at）
%
%      algorithm is based on a program written by
%      Klaus Rheinberger （University of Innsbruck）
%      （email: klaus.rheinberger@uibk.ac.at）
%
%-------------------------------------------------------------------------
clc; clear all;

%------ SPECIFY DATA ------------------------------------------------------
%------ 指定数据文件 -------------------------------------------------------
PATH= ‘H:\比赛\数据1‘; % 指定数据的储存路径
HEADERFILE= ‘100.hea‘;      % .hea 格式，头文件，可用记事本打开
ATRFILE= ‘100.atr‘;         % .atr 格式，属性文件，数据格式为二进制数
DATAFILE=‘100.dat‘;         % .dat 格式，ECG 数据
SAMPLES2READ=2048;          % 指定需要读入的样本数
                            % 若.dat文件中存储有两个通道的信号:
                            % 则读入 2*SAMPLES2READ 个数据 

%------ LOAD HEADER DATA --------------------------------------------------
%------ 读入头文件数据 -----------------------------------------------------
%
% 示例：用记事本打开的117.hea 文件的数据
%
%      117 2 360 650000
%      117.dat 212 200 11 1024 839 31170 0 MLII
%      117.dat 212 200 11 1024 930 28083 0 V2
%      # 69 M 950 654 x2
%      # None
%
%-------------------------------------------------------------------------
fprintf（1‘\\n$> WORKING ON %s ...\n‘ HEADERFILE）; % 在Matlab命令行窗口提示当前工作状态
% 
% 【注】函数 fprintf 的功能将格式化的数据写入到指定文件中。
% 表达式：count = fprintf（fidformatA...）
% 在字符串‘format‘的控制下，将矩阵A的实数数据进行格式化，并写入到文件对象fid中。该函数返回所写入数据的字节数 count。
% fid 是通过函数 fopen 获得的整型文件标识符。fid=1，表示标准输出（即输出到屏幕显示）；fid=2，表示标准偏差。
%
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 ，准备获取下一行数据
for k=1:nosig           % 读取每个通道信号的数据信息
    z= fgetl（fid1）;
    A= sscanf（z ‘%*s %d %d %d %d %d‘[15]）;
    dformat（k）= A（1）;           % 信号格式; 这里只允许为 212 格式
    gain（k）= A（2）;              % 每 mV 包含的整数个数
    bitres（k）= A（3）;            % 采样精度（位分辨率）
    zerovalue（k）= A（4）;         % ECG 信号零点相应的整数值
    firstvalue（k）= A（5）;        % 信号的第一个整数值 （用于偏差测试）
end;
fclose（fid1）;
clear A;

%------ LO