心率变异性 matlab,心率变异性好的功率谱分析方面的问题

本帖最后由 天路 于 2018-2-25 21:16 编辑

本人正在学习心率变异性方面的内容，但是按照文献上的方法做出来的结果并不是很理想，文献上说的是心率变异性的频率的范围是0.4以内，但是我做的功率谱上显示频率分布在整个频域内，试了很多种方法也没有效果。不知道怎么回事，下面是我的代码，请各位帮个忙吧，这个问题困扰很长时间了.

clc

clear

load data

data1=data(260000:360000);

%% 去除肌电信号干扰

Fs=1000;                        %采样频率

fp=10;fs=12;                    %通带截止频率，阻带截止频率

rp=1.4;rs=5;                    %通带、阻带衰减

wp=2*pi*fp;ws=2*pi*fs;

[n,wn]=buttord(wp,ws,rp,rs,'s');     %'s'是确定巴特沃斯模拟滤波器阶次和3dB截止模拟频率

Fs=500;             %采样频率

fp=20;fs=25; wp=2*pi*fp;ws=2*pi*fs; rp=1.4;rs=1.6; %通带截止频率，阻带截止频率 通带、阻带衰减

[n,wn]=buttord(wp,ws,rp,rs,'s');     %'s'是确定巴特沃斯模拟滤波器阶次和3dB截止模拟频率

[z,P,k]=buttap(n);   %设计归一化巴特沃斯模拟低通滤波器，z为极点，