将时间序列转成图像——相对位置矩阵方法 Matlab实现
admin
2024-01-28 20:25:43
0

目录

1 方法

2 Matlab代码实现

3.结果


【若觉文章质量良好且有用,请别忘了点赞收藏加关注,这将是我继续分享的动力,万分感谢!】

其他:

1.时间序列转二维图像方法及其应用研究综述_vm-1215的博客-CSDN博客

2.将时间序列转成图像——格拉姆角场方法 Matlab实现_vm-1215的博客-CSDN博客

3.将时间序列转成图像——马尔可夫转移场方法 Matlab实现_vm-1215的博客-CSDN博客

4.将时间序列转成图像——递归图方法 Matlab实现_vm-1215的博客-CSDN博客

5.将时间序列转成图像——图形差分场方法 Matlab实现_vm-1215的博客-CSDN博客

1 方法

相对位置矩阵(Relative Position Matrix, RPM)包含了原始时间序列的冗余特征,使转换后的图像中,类间和类内的相似度信息更容易被捕捉。 对于一个时间序列X=(x_t,t=1,2,\ldots,N),可以通过以下步骤得到RPM图:

  1. 针对原始时间序列,通过以下z-分值标准化的方法得到一个标准正态分布Z

    z_{t}=\frac{x_{t}-\mu}{\sigma}, t=1,2, \ldots, N

    其中\mu表示X的平均值,\sigma表示X的标准差。
  2. 采用分段聚合近似(PAA)方法,选择一个合适的缩减因子k,生成一个新的平滑时间序列\tilde X,将维度N减少到m

    \begin{aligned} &\tilde x_{i}=\left\{\begin{array}{l} \frac{1}{k} \sum_{j=k *(i-1)+1}^{k * i} z_{j}, i=1,2, \ldots, m,\left\lceil\frac{N}{k}\right\rceil-\left\lfloor\frac{N}{k}\right\rfloor=0 \\ {\left\{\begin{array}{l} \frac{1}{k} \sum_{j=k *(i-1)+1}^{k * i} z_{j}, i=1,2, \ldots, m-1 \\ \frac{1}{N-k *(m-1)} \sum_{j=k *(m-1)+1}^{N} z_{j}, i=m \end{array},\left\lceil\frac{N}{k}\right\rceil-\left\lfloor\frac{N}{k}\right\rfloor>0\right.} \end{array}\right.\\ &m=\left\lceil\frac{N}{k}\right\rceil \end{aligned}

    通过计算分段常数的平均值进行降维,可以保持原始时间序列的近似趋势,最终新的平滑时间序列\tilde X的长度为m
  3. 计算两个时间戳之间的相对位置,将预处理后的时间序列X转换为二维矩阵M

    M=\left[\begin{array}{cccc} \tilde x_{1}-\tilde x_{1} & \tilde x_{2}-\tilde x_{1} & \cdots & \tilde x_{m}-\tilde x_{1} \\ \tilde x_{1}-\tilde x_{2} & \tilde x_{2}-\tilde x_{2} & \cdots & \tilde x_{m}-\tilde x_{2} \\ \vdots & \vdots & \ddots & \vdots \\ \tilde x_{1}-\tilde x_{m} & \tilde x_{2}-\tilde x_{m} & \cdots & \tilde x_{m}-\tilde x_{m} \end{array}\right]

    如上所示,该矩阵表征了时间序列中每两个时间戳之间的相对位置关系。其每一行和每一列都以某一个时间戳为参考,进一步表征整个序列的信息;
  4. 最后利用最小-最大归一化将M转换为灰度值矩阵,最终得到相对位移矩阵F

    F=\frac{M-\min (M)}{\max (M)-\min (M)} \times 255

2 Matlab代码实现

clc
clear
close all% %%
% % 閫熷害
% speed_original = round(rand(1,111)*40+40);
% 
% m = length(speed_original);
% 
% % 浣庨�熷尯
% low_xu = round((rand()+0.1)*(m - 10));
% low_num = round((rand()+0.2)*15);
% 
% % 鏁版嵁
% speed = speed_original;
% speed(low_xu:low_xu+low_num-1) = round(rand(1,low_num)*20);
% 
% % X = speed;
% 
% % 璁烘枃閲岀殑鏁版嵁
% X = [0 1 2 1 2 3 4 3 2 3 2 1];
% 
% N = length(X);%% 鐢熸垚鏁版嵁
speed = xlsread('3_1_link6_28_5_30min.csv');
% speed = xlsread('3_1_link1_1_5_30min.csv');%% z-鍒嗗�兼爣鍑嗗寲
X = speed';
mu = mean(X);
sigma = sqrt(var(X));Z = (X-mu)/sigma;%% PAA
k = 2;
N = length(X);
m = ceil(N/k);if ceil(N/k)-floor(N/k) == 0for i = 1:mX2(i) = 1/k * sum(Z(k*(i-1)+1:k*i));end
elsefor i = 1:m-1X2(i) = 1/k * sum(Z(k*(i-1)+1:k*i));endX2(m) = 1/(N-k*(m-1)) * sum(Z(k*(m-1)+1:N));
end%% 鐢熸垚RPM鍥�
M = repmat(X2,m,1) - repmat(X2',1,m);F = (M - min(M(:))) / (max(M(:))) - min(M(:)) * 255;h = figure(2);
set(gcf,'position',[800 300 500 420])
imagesc(F)
title('鐩稿浣嶇Щ鐭╅樀(RPM)')
saveas(h,'RPM_11.bmp')

3.结果

【若觉文章质量良好且有用,请别忘了点赞收藏加关注,这将是我继续分享的动力,万分感谢!】

相关内容

热门资讯

诺基亚8刷安卓系统,解锁无限可... 你手中的诺基亚8是不是已经有点儿落伍了呢?别急,今天就来给你支个招,让你的老伙计焕发新生,变身安卓小...
安卓系统能不能,可以。 你有没有想过,安卓系统到底能不能?这个问题,就像是在问一个老朋友,他是不是真的懂你。安卓系统,这个陪...
安卓系统恢复误删视频,轻松找回... 手机里的视频突然不见了,是不是你也遇到了这样的尴尬情况?别急,今天就来教你如何用安卓系统恢复误删的视...
华为安卓系统的siri,华为安... 你知道吗?华为最近在安卓系统上搞了个大动作,那就是推出了自己的Siri——华为助手。这可真是让人眼前...
wp模拟安卓系统界面,畅游虚拟... 你有没有想过,在电脑上也能体验到安卓系统的流畅与便捷呢?没错,这就是今天我要跟你分享的神奇小玩意——...
安卓系统的开发团队,谷歌开发团... 你知道吗?在科技的世界里,有一个团队可是默默无闻地创造了无数奇迹,他们就是安卓系统的开发团队。这个团...
俄语流利说安卓系统,轻松掌握俄... 你有没有想过,学习一门新语言竟然可以变得如此轻松有趣?没错,我要给你安利一款神器——俄语流利说安卓系...
安卓P系统原装铃声,唤醒科技之... 你有没有发现,手机里的那些原装铃声,有时候比我们自己的手机铃声还要动听呢?尤其是安卓P系统的原装铃声...
稳定无广告安卓系统,探索稳定无... 你有没有想过,手机系统就像是我们生活的环境,有时候干净整洁,有时候却满是杂乱无章的广告?今天,我要给...
安卓系统隔离运行app,技术革... 你知道吗?在智能手机的世界里,安卓系统可是个超级明星呢!它不仅功能强大,而且兼容性极好,几乎所有的手...
佳博3120安卓系统,引领移动... 你有没有听说过佳博3120安卓系统?这款设备最近可是火得一塌糊涂呢!想象一台集成了安卓系统的打印机,...
安卓系统放音乐全屏,沉浸式听觉... 你有没有发现,用安卓手机放音乐的时候,有时候屏幕会自动全屏显示,这可真是挺有趣的。你知道吗?这个小小...
安卓子系统是win,基于Win... 你知道吗?在科技的世界里,总是充满了惊喜和未知。今天,我要给你揭秘一个你可能没听说过的秘密:安卓子系...
金刚导航升级安卓系统,畅享智能... 你知道吗?最近金刚导航可是来了一次大变身呢!没错,就是那个我们平时出行时离不开的导航神器——金刚导航...
安卓系统有话筒软件,畅享便捷沟... 你有没有发现,手机里的安卓系统里竟然藏着这么一个神奇的小玩意儿——话筒软件!没错,就是那个可以让你随...
安卓备份系统版本更新,版本更新... 你知道吗?最近安卓备份系统又来了一次大更新,这可真是让人兴奋不已呢!想象你的手机备份功能变得更加智能...
htc怎么降级安卓系统,轻松恢... 你有没有发现,有时候手机系统更新后,新功能虽然多了,但速度却慢了下来,甚至有些功能还不太好用?这不,...
索尼电视安卓系统优点,畅享无限 亲爱的电视迷们,你是否在寻找一款既能满足你对画质追求,又能让你畅享智能生活的电视呢?今天,就让我带你...
智能驭领系统安卓,引领未来智能... 你有没有发现,最近手机界又掀起了一股热潮?没错,就是那个让人眼前一亮的智能驭领系统安卓!今天,就让我...
安卓u15.1系统,功能升级与... 你有没有发现,最近你的安卓手机更新到了U15.1系统?别小看了这个小小的升级,它可是带来了不少惊喜呢...