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;%分段重拍后的數組元素總個數
%對沒有采樣值但需要出