心率变异性 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为极点,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_lvbo=filter(bz,az,data1);
%% RR间期序列提取
data_qjx=qjxpj1(data_lvbo); %去除基线漂移
[data_Rf,t]=RFtiqu(data_qjx); %提取R峰
s=1:length(data1);
data_rr=diff(t)/Fs; %生成RR间期序列
figure
plot (s,data_qjx,t,data_Rf,'r+') %画R峰图
data_rr(find(data_rr>1.3*mean(data_rr)))=[];
data_rr(find(data_rr<0.7*mean(data_rr)))=[];
t(find(data_rr>1.3*mean(data_rr)))=[];
t(find(data_rr<0.7*mean(data_rr)))=[];
%%
%均匀重采样
fs_resample=10;%重采样频率
f_RR_resample=t(1):Fs/fs_resample:t(end-1);%确定插值频率
t(end)=[];
RR_interp=interp1(t,data_rr,f_RR_resample,'spline');
% RR_interp=interp1(t,real_RR,f_RR_resample,'spline');%插值,重采样
figure;
plot(f_RR_resample/fs,RR_interp,'.');
% RR_interp=RR_interp-mean(RR_interp); %该语句的实质作用?去除直流分量
RR_interp=QZL(RR_interp);
%%窗函数法计算功率谱并画图
nfft=2^(nextpow2(length(RR_interp)-1));
window=hamming(length(RR_interp));
[Pxx,f]=pwelch(RR_interp,window,0,nfft,fs_resample,'onesided');
Pxx=Pxx*2;
figure;
plot(f,Pxx);%实信号,频谱有正负 ,乘以2是真实功率谱
xlabel('f/Hz');ylabel('功率');title('瞬时心率的功率谱')
[px3,f3]=pyulear(data_rr,16,nfft,10);
figure
plot(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 上传
汉明窗求功率谱
心率变异性 matlab,心率变异性好的功率谱分析方面的问题相关推荐
- ar谱matlab实验,用MATLAB进行AR模型功率谱分析
用MATLAB 进行AR 模型功率谱分析 随机信号序列x(n)是均值为0方差为1的高斯型白噪声经过AR 模型 ()4 3219606.01697.29403.22137.211 ----+-+-= z ...
- matlab ar谱分析,用MATLAB进行AR模型功率谱分析
用MATLAB进行AR模型功率谱分析 用MATLAB进行AR模型功率谱分析 随机信号序列x(n)是均值为0方差为1的高斯型白噪声经过AR模型 H z 1 1 2.2137z 1 2.9403z 2 2 ...
- 心率变异性 matlab,心率变异性及其相关算法的实现.ppt
心率变异性及其相关算法的实现 附加功能参考文献 附加功能 该算法除了可以计算心率变异性之外,还有以下2种附加功能: 对心率的正常与否进行判断,输出有心率正常.心动过速.心动过缓3种情况: 对心脏的早搏 ...
- matlab实验——信号和噪声产生及其功率谱分析
文章目录 前言 一.高斯白噪声自相关函数及功率谱 1.matlab代码 2.运行结果 二.均匀白噪声自相关函数及功率谱 1.matlab代码 2.运行结果 三.正弦波与高斯白噪声叠加 1.matlab ...
- 脑电功率谱分析matlab,基于matlab的功率谱分析方法研究 毕业论文.doc
基于matlab的功率谱分析方法研究 摘 要 数字信号处理(DSP)重要的应用领域之一,是建立在周期信号和随机信号基础上的功率谱估计.在实际应用中往往不能获得具体信号的表达式,需要根据有限的数据样本来 ...
- Matlab实现经典功率谱分析和估计
Matlab 实现经典功率谱分析和估计 Matlab 实现经典功率谱分析和估计 功率谱 Matlab使用 1.直接法 2.间接法 3.改进直接法:Bartlett法 4.Welch法 附上谋篇论文,分 ...
- 功率谱分析笔记-------脑电相关
1:功率谱分析的方法介绍 功率谱分析的方法大致可以分为两大类:第一类是经典的功率谱计算方法,第二类是现代功率谱计算方法,如图1所示.其中第一类经典功率谱分析方法,又可以分为直接法.间接法和改进的直接法 ...
- matlab验证dtft移位性质,matlab 实验四 信号的谱分析.doc
matlab 实验四 信号的谱分析.doc 实验四信号的谱分析一.实验目的:1.掌握DTFT原理及其程序实现,学习用DTFT对信号进行谱分析.2.掌握DFT原理及其程序实现,学习用DFT对信号进行谱分 ...
- 风电功率预测matlab,MatlAB毕业论文:风电功率预测问题.doc
您所在位置:网站首页 > 海量文档  > 学术论文 > 自然科学论文 MatlAB毕业论文:风电功率预测问题.doc2 ...
最新文章
- Windows10看视频出现蓝屏重启
- 二叉树的深搜(DFS)与广搜(BFS)
- Spring Data(二)查询
- unity实现图片轮播效果_Unity3D实现列表拖拽轮播分页滚动功能
- BugkuCTF-MISC题这是一张单纯的图片
- 一步设置Intellij IDEA 热部署处理方法
- Qone 自动删除说说脚本
- 《WCF技术内幕》翻译35:第2部分_第6章_通道:通道功能
- 使用SAX读取XML文件
- excel常用快捷键汇总
- 无光驱、软驱、USB全NTFS格式硬盘上安装WinXP
- uc看视频显示服务器有点忙,uc浏览器常见问题集锦(一)
- 移动端项目实现手机定位
- 编程入门先要学什么?打怪升级走起来!
- N卡和A卡各型号参数对比收集
- 计算机 蓝牙鼠标卡顿,无线鼠标卡顿不流畅是什么原因_无线鼠标卡顿不流畅的处理方法...
- docker 网络配置
- IDEA插件系列(87):Base64 image encoder插件——以base64编码的形式查看图像
- 《我的眼睛--图灵识别》第十一章:实战演练:图像类识别
- java毕业设计项目源代码javaweb在线商城电商项目|购物
热门文章
- c语言写modbus协议,C语言编写modbus协议
- c51的八位智力抢答器C语言编程,基于AT89C51单片机的八位智力抢答器设计.doc
- 【flutter】 集成HMS需要将证书文件打包到APK中,请直接将assets目录拷贝到应用工程根目录。
- 推行六西格玛管理的成本效益高吗?
- 原生混合云 — 经政企打磨方能赢得政企信任
- Instrumentation框架分析及其使用
- 学生用计算机怎么关屏,电脑自带锁屏使用方法
- Linux下nano编辑器的快捷键使用
- linux 远程解压tar.gz,linux解压tar.gz
- python爬虫抓包工具_「docker实战篇」python的docker爬虫技术-fiddler抓包软件详细配置(七)...