本文主要是介绍12导联 ECG 和 VCG 相互变换的方法,希望对大家解决编程问题提供一定的参考价值,需要的开发者们随着小编来一起学习吧!
%12导联ECG 计算 VCGclc;close all; close all;D = [-0.515 0.157 -0.9170.044 0.164 -1.3870.882 0.098 -1.2771.213 0.127 -0.6011.125 0.127 -0.0860.831 0.076 0.2300.632 -0.235 0.0590.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;
%plot(tm, sig);figure;
subplot(311);plot(tm, sig(:,13)); hold; plot(tm, VCG(:,1),'r'); title('X');
subplot(312);plot(tm, sig(:,14)); hold; plot(tm, VCG(:,2),'r'); title('Y');
subplot(313);plot(tm, sig(:,15)); hold; plot(tm, VCG(:,3),'r'); title('Z');
legend('XYZ ori', 'XYZ derived');
这篇关于12导联 ECG 和 VCG 相互变换的方法的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!