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??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