脉搏信号处理y05.doc

数字信号处理综合训练 说 明 书 题目：脉搏信号间期序列的时域、频域分析 学院：电气工程与信息工程学院 班级：电子信息科学与技术（1）班 姓名：马艳霞 学号 2011年7月15日 目录 1 脉搏信号处理的基本思路…………………………….1 信号预处理……………………………………………..2 3 脉搏信号的时域分析…………………………………..5 4 脉搏信号频域分析……………………………………..8 程序清单………………………………………………..9 6 心得体会………………………………………………..23 参考文献………………………………………………..24 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 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;%通带波纹系数; 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; 脉搏信号通过零相位带通滤波器后，高频信号和基线漂移信号有所衰减，信号变得清晰。 图2 %小波包消噪； %进行一维小波消噪 figure(3); z=wden(y1,heursure,s,mln,3,sym8);%[XD,CXD,LXD]=wden(X,tptr,sorh,scal,n,’wavename’)使用小波系数阈值，返回输入信号X除噪后的信号XD，输出参数[CXD,LXD]表示XD的小波分解结构。 %输入参数中，tptr同thselect（）函数；sorh为‘ s’或‘h’表示软硬阈值；n表示在n层上的小波分解；wavename指定小波名称；scal定义阈值调整比例： %‘one’不设定比例； %‘sln’使用的基于第单层系数噪声估计，设置比例； %‘mln’用噪声层的层相关估计，调整比例。 subplot(2,1,1); plot(z); title(小波消噪后的信号); xlabel(时间/ms); ylabel(幅度); grid on; %小波包消噪