脉搏信号特征提取

function[HR_UU_posU_ampP_posR_amp]=RPT_analyse2（FILENAME）
% RPT_analyse2.m       脉搏信号特征提取和心率计算
% HR_U                心率
% U_pos               升支起点出现时刻
% U_amp               升支起点幅值
% P_pos               主波峰值出现时刻
% R_amp               主波峰幅值
% FILENAME            输入脉搏信号
% RPT_analyse2           用红色圆圈标出全部主波峰值点Ｐ，用蓝色圆圈标出全部升支起点Ｕ

% Zhou Yihui Jiang Wangshu Cui Xi Liu Xue 修改于2018-05-30

%% 平滑曲线
RPT_smooth1=smooth（FILENAME）;
RPT_smooth2=smooth（RPT_smooth1）;
%% 找出所有主波峰值，计算脉搏周期
% 找出所有主波峰值
PRT_extrmax_ind=find（diff（sign（diff（RPT_smooth2）））==-2）+1;%差分法找所有极大值点 该数组中存放的是极大值对应的元素下标 
rpt_max=max（RPT_smooth2（PRT_extrmax_ind））;%找到所有极大值中的最大值
rpt_max_pre=RPT_smooth2（PRT_extrmax_ind）;%存放所有极大值
rpt_max_real= abs（rpt_max_pre-rpt_max）<0.1;%找出主波峰值点在[存放极值横坐标的数组]（PRT_extrmax_ind）中对应的下标
rpt_max_real_x=PRT_extrmax_ind（rpt_max_real）;%主波峰值点的横坐标
% 计算脉搏周期
tt_diff=diff（rpt_max_real_x）;%差分法 计算各主峰间距
period=mean（tt_diff）;%取平均值
%% 通过周期将RPT所有脉搏值分段  利用reshape对数组重排以分段  

% 说明：max函数只能在完整的一组数据中找出最大值及其下标，不能对一个数组分段操作，若想分段寻找仍需要将该数组中所有元素重新放到几个数组中，故尚未修改此方法。

l_RPT=length（FILENAME）;

%不在最后一段找最大最小值的程序（可选择）
n=floor（l_RPT/period）;%分段（周期）个数，不满1的舍去
nn=n*period;%分段重拍后的数组元素总个数
%将需要查找最大最小值的采样点值复制到一个新的数组
for jj=1:nn
RPT_1（jj）=FILENAME（jj）;
end

%保留最后一段的程序
%n=ceil（l_RPT/period）;%分段（周期）个数，不满1的加一
%nn=n*period;%分段重拍后的数组元素总个数
%对没有采样值但需要出