各地残疾人联合会网站建设,四川建设公共查询平台,福鼎网页设计,企业网站建设的目的本帖最后由 天路 于 2018-2-25 21:16 编辑本人正在学习心率变异性方面的内容#xff0c;但是按照文献上的方法做出来的结果并不是很理想#xff0c;文献上说的是心率变异性的频率的范围是0.4以内#xff0c;但是我做的功率谱上显示频率分布在整个频域内#xff0c;试了很多…本帖最后由 天路 于 2018-2-25 21:16 编辑本人正在学习心率变异性方面的内容但是按照文献上的方法做出来的结果并不是很理想文献上说的是心率变异性的频率的范围是0.4以内但是我做的功率谱上显示频率分布在整个频域内试了很多种方法也没有效果。不知道怎么回事下面是我的代码请各位帮个忙吧这个问题困扰很长时间了.clcclearload datadata1data(260000:360000);%% 去除肌电信号干扰Fs1000; %采样频率fp10;fs12; %通带截止频率阻带截止频率rp1.4;rs5; %通带、阻带衰减wp2*pi*fp;ws2*pi*fs;[n,wn]buttord(wp,ws,rp,rs,s); %s是确定巴特沃斯模拟滤波器阶次和3dB截止模拟频率Fs500; %采样频率fp20;fs25; wp2*pi*fp;ws2*pi*fs; rp1.4;rs1.6; %通带截止频率阻带截止频率 通带、阻带衰减[n,wn]buttord(wp,ws,rp,rs,s); %s是确定巴特沃斯模拟滤波器阶次和3dB截止模拟频率[z,P,k]buttap(n); %设计归一化巴特沃斯模拟低通滤波器z为极点p为零点和k为增益[bp,ap]zp2tf(z,P,k); %转换为Ha(p),bp为分子系数ap为分母系数[bs,as]lp2lp(bp,ap,wp) ;%Ha(p)转换为低通Ha(s)并去归一化bs为分子系数as为分母系数[bz,az]bilinear(bs,as,Fs); %对模拟滤波器双线性变换data_lvbofilter(bz,az,data1);%% RR间期序列提取data_qjxqjxpj1(data_lvbo); %去除基线漂移[data_Rf,t]RFtiqu(data_qjx); %提取R峰s1:length(data1);data_rrdiff(t)/Fs; %生成RR间期序列figureplot (s,data_qjx,t,data_Rf,r) %画R峰图data_rr(find(data_rr1.3*mean(data_rr)))[];data_rr(find(data_rr0.7*mean(data_rr)))[];t(find(data_rr1.3*mean(data_rr)))[];t(find(data_rr0.7*mean(data_rr)))[];%%%均匀重采样fs_resample10;%重采样频率f_RR_resamplet(1):Fs/fs_resample:t(end-1);%确定插值频率t(end)[];RR_interpinterp1(t,data_rr,f_RR_resample,spline);% RR_interpinterp1(t,real_RR,f_RR_resample,spline);%插值重采样figure;plot(f_RR_resample/fs,RR_interp,.);% RR_interpRR_interp-mean(RR_interp); %该语句的实质作用去除直流分量RR_interpQZL(RR_interp);%%窗函数法计算功率谱并画图nfft2^(nextpow2(length(RR_interp)-1));windowhamming(length(RR_interp));[Pxx,f]pwelch(RR_interp,window,0,nfft,fs_resample,onesided);PxxPxx*2;figure;plot(f,Pxx);%实信号频谱有正负 乘以2是真实功率谱xlabel(f/Hz);ylabel(功率);title(瞬时心率的功率谱)[px3,f3]pyulear(data_rr,16,nfft,10);figureplot(f3,px3);title(AR MODEL)下面附计算结果的图片AR模型功率谱法.jpg(34.52 KB, 下载次数: 3)2018-2-25 16:56 上传AR模型求功率谱R峰监测.jpg(51.4 KB, 下载次数: 0)2018-2-25 16:56 上传R峰监测窗函数法功率谱.jpg(45.86 KB, 下载次数: 0)2018-2-25 16:59 上传汉明窗求功率谱
