function[HR_UU_posU_ampP_posR_amp]=RPT_analyse2（FILENAME）
%?RPT_analyse2.m???????脈搏信號(hào)特征提取和心率計(jì)算
%?HR_U????????????????心率
%?U_pos???????????????升支起點(diǎn)出現(xiàn)時(shí)刻
%?U_amp???????????????升支起點(diǎn)幅值
%?P_pos???????????????主波峰值出現(xiàn)時(shí)刻
%?R_amp???????????????主波峰幅值
%?FILENAME????????????輸入脈搏信號(hào)
%?RPT_analyse2???????????用紅色圓圈標(biāo)出全部主波峰值點(diǎn)Ｐ，用藍(lán)色圓圈標(biāo)出全部升支起點(diǎn)Ｕ

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

%%?平滑曲線
RPT_smooth1=smooth（FILENAME）;
RPT_smooth2=smooth（RPT_smooth1）;
%%?找出所有主波峰值，計(jì)算脈搏周期
%?找出所有主波峰值
PRT_extrmax_ind=find（diff（sign（diff（RPT_smooth2）））==-2）+1;%差分法找所有極大值點(diǎn)?該數(shù)組中存放的是極大值對(duì)應(yīng)的元素下標(biāo)?
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;%找出主波峰值點(diǎn)在[存放極值橫坐標(biāo)的數(shù)組]（PRT_extrmax_ind）中對(duì)應(yīng)的下標(biāo)
rpt_max_real_x=PRT_extrmax_ind（rpt_max_real）;%主波峰值點(diǎn)的橫坐標(biāo)
%?計(jì)算脈搏周期
tt_diff=diff（rpt_max_real_x）;%差分法?計(jì)算各主峰間距
period=mean（tt_diff）;%取平均值
%%?通過(guò)周期將RPT所有脈搏值分段??利用reshape對(duì)數(shù)組重排以分段??

%?說(shuō)明：max函數(shù)只能在完整的一組數(shù)據(jù)中找出最大值及其下標(biāo)，不能對(duì)一個(gè)數(shù)組分段操作，若想分段尋找仍需要將該數(shù)組中所有元素重新放到幾個(gè)數(shù)組中，故尚未修改此方法。

l_RPT=length（FILENAME）;

%不在最后一段找最大最小值的程序（可選擇）
n=floor（l_RPT/period）;%分段（周期）個(gè)數(shù)，不滿(mǎn)1的舍去
nn=n*period;%分段重拍后的數(shù)組元素總個(gè)數(shù)
%將需要查找最大最小值的采樣點(diǎn)值復(fù)制到一個(gè)新的數(shù)組
for?jj=1:nn
RPT_1（jj）=FILENAME（jj）;
end

%保留最后一段的程序
%n=ceil（l_RPT/period）;%分段（周期）個(gè)數(shù)，不滿(mǎn)1的加一
%nn=n*period;%分段重拍后的數(shù)組元素總個(gè)數(shù)
%對(duì)沒(méi)有采樣值但需要出