12导联 ECG 和 VCG 相互变换的方法

%12导联ECG 计算 VCG






clc;close all; close all;

D = [-0.515 0.157 -0.917
    0.044 0.164 -1.387
    0.882 0.098 -1.277
    1.213 0.127 -0.601
    1.125 0.127 -0.086
    0.831 0.076 0.230
    0.632 -0.235 0.059
    0.235 1.066 -0.132];









T = (D'*D)\(D');
% T =[-0.172 -0.074 0.122 0.231 0.239 0.194 0.156 -0.010
%     0.057 -0.019 -0.106 -0.022 0.041 0.048 -0.227 0.887
%     -0.229 -0.310 -0.246 -0.063 0.055 0.108 0.022 0.102]




% read data of 12-lead ECG
cd 'G:\data\physionet-data\ptb-diagnostic-ecg-database-1.0.0\ptb-diagnostic-ecg-database-1.0.0\patient003'
record = 's0017lre';
[sig, Fs, tm] = rdsamp(record,[],2000);






s = sig(:, [7:12 1 2])';
VCG = T * s;
VCG = VCG';

%figure;