好多地方都有求数据包络线的程序这里不再赘述,比如链接:
添加链接描述
图像展示:
但是对于时间序列数据我们求包络线的过程中有时候想保留对应的时间值,而不仅仅是看个热闹。
上面链接中第三步:
第一个问题:
“包络线的生成,需要将原数据分成等宽度的小段,每段取一个最高点所对应的横坐标和纵坐标”,其实它根根没有取每段最高点多对应的横坐标,从最后一步生成的横坐标可以看出来。
第二个问题:
“比如设定每小段的长度为 d = 4 或 8 或 16(或2048的其他约数),”这个每小段的长度设定有限制,导致程序的实用性大大降低。
对此我们做出改进
d = 5; % 可以设定每小段长度范围 2 ~n
Loopn = round(length(Error2)/d);
ClusteringData = zeros(2,Loopn);
for k=1:Loopn-1
[ma,index] = max(Error2((k-1)*d+1:(k)*d));
ClusteringData(1,k) = ma; % 找到每次分段中的最大值
ClusteringData(2,k) = T2(index + d*(k-1)); % 储存对应的时间值
end
[ma,index] = max(Error2((Loopn-1)*d+1:length(Error2)));%最后一段单独处理
ClusteringData(1,Loopn) = ma; % 找到分段中的最大值
ClusteringData(2,Loopn) = T2(index + d*(Loopn-1)); % 储存对应的时间值
plot(ClusteringData(2,:),ClusteringData(1,:),"-*b");%图像展示
说明:原始数据为:
T2(横坐标),Error2(纵坐标)
最终得到的上侧包络线:ClusteringData(2,:),ClusteringData(1,:)
至于下侧包络线和上面一样的思路。