matlab 求时间序列数据的包络线

好多地方都有求数据包络线的程序这里不再赘述,比如链接:
添加链接描述
图像展示:
在这里插入图片描述
但是对于时间序列数据我们求包络线的过程中有时候想保留对应的时间值,而不仅仅是看个热闹。
上面链接中第三步:
第一个问题:
“包络线的生成,需要将原数据分成等宽度的小段,每段取一个最高点所对应的横坐标和纵坐标”,其实它根根没有取每段最高点多对应的横坐标,从最后一步生成的横坐标可以看出来。
第二个问题:
“比如设定每小段的长度为 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,:)
至于下侧包络线和上面一样的思路。

评论 2
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值