脉搏信号处理y105

数字信号处理综合训练 说 明 书 题目：脉搏信号间期序列的时域、频域分析 学院：电气工程与信息工程学院 班级：电子信息科学与技术（1）班 姓名：*** 学号：******** 2011年7月15日 8 目录 1 脉搏信号处理的基本思路… … … … … … … … … … … .1 2 信号预处理… … … … … … … … … … … … … … … … … ..2 3 脉搏信号的时域分析… … … … … … … … … … … … … ..5 4 脉搏信号频域分析… … … … … … … … … … … … … … ..8 5 程序清单… … … … … … … … … … … … … … … … … … ..9 6 心得体会… … … … … … … … … … … … … … … … … … ..23 7 参考文献… … … … … … … … … … … … … … … … … … ..24 9 5 程序清单 %加载采样数据； y=[y10 y11 y12 y13 y14 y15]; t=1:6000; %时域范围； figure(1); subplot(2,1,1); y=-y; plot(t,y); grid on; title('脉搏信号'); xlabel('时间/ms'); ylabel('相对幅度'); %截取整数倍周期 data1=min(y(1:500)); data2=min(y(5000:5500)); for i=1:500 if y(i)==data1 m1=i; end end for j=5000:5500 if y(j)==data2 m2=j; end end 1 0 mb=m2-m1; y=y(m1:m2); t=1:mb+1; subplot(2,1,2); plot(t,y); title('脉搏信号'); xlabel('时间/ms'); ylabel('相对幅度'); grid on; 原始的脉搏信号和经过整数倍提取后的信号如下图所示 图1 %带通滤波器,虑除工频.基线漂移和肌电干扰； figure(2); fs=1500;%采样频率； wp=[0.9,30]/(fs/2);%设置通带截止频率; ws=[0.1,200]/(fs/2);%设置阻带截止频率; rp=0.5;%通带波纹系数; 1 1 rs=40;%阻带波纹系数; [N,wc]=buttord(wp,ws,rp,rs); [num,den]=butter(N,wc);%滤波器分子分母系数向量； [H,W]=freqz(num,den); subplot(2,1,1); plot(fs*W/(2*pi),abs(H)); title('butterworth带通滤波器幅频响应'); xlabel('频率/Hz'); ylabel('幅度'); grid on; y1=filtfilt(num,den,y);%脉搏信号通过零相位滤波器； subplot(2,1,2); plot(y1); title('脉搏信号通过带通滤波器'); xlabel('时间/ms'); ylabel('幅度'); grid on; 脉搏信号通过零相位带通滤波器后，高频信号和基线漂...