心电信号采集与分析MATLAB

心电信号处理是国内外近年来迅速发展的一个研究热点，是现代生命科学研 究的重要组成部分，其目的是为了从获得的信号中提取有用信息。心电图反应人 体心脏工作状况，各个波形的不同形式往往体现了某些病变。本实验通过对心电 信号的产生、采集、分析、处理，可以有效的监测人的心脏和血压的健康状况。

clc;
clear;
close all;

%% 参数设置
N=1024;
fs=300;
ff=input（‘input f/Hz\n‘）
Ts=1/fs;
n=0:N-1;
fp=n/N*fs;

%% 时、频波形
[FileNamePathName] = uigetfile（‘USB.dat‘‘Select the USB.dat file‘）;
f = fullfile（PathNamefilesepFileName）;
fid = fopen（f‘r‘）;
data = fscanf（fid‘%x‘）;
fclose（fid）;
data = data（1:2:end）*256 + data（2:2:end）;
figure（1）;
plot（data）;     
data=data-mean（data）;
b=max（data）;
a=min（data）;
data=data/（b-a）;
F=fft（data）;
figure（2）;
plot（abs（F））;

%%  滤波后的时、频波形
fp=0.05;fc=120;As=80;Ap=1;Fs=300;
wc=2*pi*fc/Fs; wp=2*pi*fp/Fs;
wd=wc-wp;
beta=0.1102*（As-8.7）;
N=ceil（（As-7.95）/2.286/wd）;
wn= kaiser（N+1beta）; 
ws=（wp+wc）/2/pi;
b=fir1（Nwswn）;
x=fftfilt（bdata）;
X=fft（x）;
figure（3）;
plot（x）;title（‘低通滤波后的信号‘）;
figure（4）;
plot（abs（X））;

%% 信噪比计

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     目录           0  2019-03-28 15:01  程序+采集数据\
     文件        8192  2017-08-29 15:17  程序+采集数据\0.8hz.dat
     文件        8192  2017-08-28 14:58  程序+采集数据\1.5hz.dat
     文件      448539  2019-03-28 14:59  程序+采集数据\1.docx
     文件        8192  2017-08-29 15:12  程序+采集数据\1hz.dat
     文件        8192  2017-08-29 15:04  程序+采集数据\2hz.dat
     文件        8192  2017-08-30 10:34  程序+采集数据\400.dat
     文件        8192  2017-08-30 11:00  程序+采集数据\500.dat
     文件        8192  2017-08-29 15:48  程序+采集数据\800m.dat
     文件        8192  2017-08-30 11:06  程序+采集数据\USB.dat
     文件        8192  2017-08-29 11:20  程序+采集数据\USB_100.dat
     文件        3232  2017-09-01 13:57  程序+采集数据\a1.m