%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;
12导联 ECG 和 VCG 相互变换的方法
于 2023-03-21 11:43:55 首次发布