它使用
rmoutliers
删除的合成曲线,然后计算的半最大值宽度。你可以尽情发挥
rmoutliers
得到预期的曲线或使用不同的方法如果你找到更合适的。
其他approaces:
- 符合高斯分布的曲线和使用合适的参数来计算应用。
- 使用findpeaks随着半高宽度参考。看到演示。
data =负载(“FWHM_Help_Plot_Data.mat”);
%绘制原始数据
h(1) =(数据。X data.Z“DisplayName的”,“RawData”);
持有在
%去除离群值
Zm评选,idx] = rmoutliers (data.Z“movmedian”,300);
h(2) =情节(data.X (~ idx), Zm评选,“DisplayName的”,“rmOutlier”);
%计算应用
halfMax = max (Zm评选)/ 2;%简单自基线接近0
lrIdx =[找到(Zm评选> = halfMax 1“第一”),发现(Zm评选> = halfMax 1“最后一次”));% L / R指数宽度
应用= diff (data.X (lrIdx));
%显示范围
h (3) = yline (halfMax,“k——”,“DisplayName的”,“HalfHeight”);
h(4) =参照线(data.X (lrIdx (1)),“DisplayName的”,“WidthBounds”);
参照线(data.X (lrIdx (2)));
传奇(h)
标题(sprintf (= % g ' '的半最大值宽度、应用)
