心率变异性分析 （Frequency_HRV）

利用AR模型和Lomb算法的心率变异性分析。

clear all

rr2=load（‘D:\1001.txt‘）;

% [rr_statdetr] = detrend（rr2301000）


N=length（rr2）;
mean2=mean（rr2）/1000;
mhrf=1/mean2;
[ae]=arburg（rr210）;
[Pf]=pburg（rr210[2048][mhrf]）;
% P（1）=P（2）
% plot（f10*log10（P））
figure
plot（f10*log10（P））
% figure
% plot（log（f）log（P）‘-*‘）
% [v lf hf] = lfhf（P）;
% v
% lf
% hf
space=f（2）-f（1）;
ivf=find（（0<=f）&（f<0.04））
ilf = find（ （0.04 <= f） & （f <= 0.15） ）;
ihf = find（ （0.15  < f） & （f <= 0.40） ）;
Pvlf=sum（P（ivf）*space）
Plf = sum（P（ilf））*space
Phf = sum（P（ihf））*space
Ptf=Pvlf+Plf+Phf
v = Plf/Phf;


% m = mean（rr）;
% lambda = 30;
% I = speye（N）;
% D2 = spdiags（ones（N-21）*[1 -2 1][0:2]N-2N）;
% M = I-inv（I+lambda^2*D2‘*D2）;
% rr_stat = （I-inv（I+lambda^2*D2‘*D2））*rr;

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     目录           0  2014-07-27 05:22  Frequency_HRV\
     文件         760  2014-07-27 05:21  Frequency_HRV\AR.m
     文件        2368  2011-10-11 22:13  Frequency_HRV\lomb.m