matlab画相关系数图(长时序--分时期分月份)

成图如下:

a图为例子,纵坐标是1-12月,横坐标是不同的因子。如左下第一格表示的就是x1与AO因子在十二月的相关性。

(a/b/c是x1因子,d/e/f是x2因子)

首先要计算出每个相关系数(可以在SPSS中完成),然后将它们按矩阵型输入到excel中,如下:

接着在matlab中读入数据,直接画图:

Matlab 复制代码
corr = readtable('AONAODA.xls','Sheet','Sheet2');
A1 = table2array(corr(1:12,2:4));
A2 = table2array(corr(1:12,7:9));
A3 = table2array(corr(1:12,12:14));
A4 = table2array(corr(15:26,2:4));
A5 = table2array(corr(15:26,7:9));
A6 = table2array(corr(15:26,12:14));
% CMap=slanCM(102);
CMap=slanCM(95);


figure;
sub1 = subplot(2,3,1);
imagesc(A2);
caxis([-1 1]);
xticklabels({'AO','NAO','DA'});
yticks([1:2:12]);
yticklabels({'Jan','Mar','May','Jul','Sep','Nov'});
pos=set(sub1,'Position', [0.06 0.58 0.22 0.38]);  %设置位置
txt = {'(a) 1984-2000'};
text(0.1,1,txt)

% yticklabels({'Jan','Feb','Mar','Apr','May','Jun','Jul','Aug','Sep','Oct','Nov','Dec'});

sub2 = subplot(2,3,2);
imagesc(A3);
caxis([-1 1]);
xticklabels({'AO','NAO','DA'});
yticks([1:2:12]);
yticklabels({'Jan','Mar','May','Jul','Sep','Nov'});
pos=set(sub2,'Position', [0.36 0.58 0.22 0.38]);
txt = {'(b) 2001-2022'};
text(0.1,1,txt)

sub3 = subplot(2,3,3);
imagesc(A1);
caxis([-1 1]);
xticklabels({'AO','NAO','DA'});
yticks([1:2:12]);
yticklabels({'Jan','Mar','May','Jul','Sep','Nov'});
txt = {'(c) 1984-2022'};
text(0.1,1,txt)
pos=set(sub3,'Position', [0.66 0.58 0.22 0.38]);

sub4 = subplot(2,3,4);
imagesc(A5);
caxis([-1 1]);
xticklabels({'AO','NAO','DA'});
yticks([1:2:12]);
yticklabels({'Jan','Mar','May','Jul','Sep','Nov'});
pos=set(sub4,'Position', [0.06 0.1 0.22 0.38]);
txt = {'(d) 1984-2000'};
text(0.1,1,txt)

sub5 = subplot(2,3,5);
imagesc(A6);
caxis([-1 1]);
xticklabels({'AO','NAO','DA'});
yticks([1:2:12]);
yticklabels({'Jan','Mar','May','Jul','Sep','Nov'});
pos=set(sub5,'Position', [0.36 0.1 0.22 0.38]);
txt = {'(e) 2001-2022'};
text(0.1,1,txt)

sub6=subplot(2,3,6);
imagesc(A4);
caxis([-1 1]);
xticklabels({'AO','NAO','DA'});
yticks([1:2:12]);
yticklabels({'Jan','Mar','May','Jul','Sep','Nov'});
pos=set(sub6,'Position', [0.66 0.1 0.22 0.38]);
txt = {'(f) 1984-2022'};
text(0.1,1,txt)

colormap(CMap);


h=colorbar('eastoutside','fontsize',12);
set(get(h,'Title'),'string','Corr','fontsize',9)
set(h,'YTick',-1:1:1,'fontsize',9);%设置colorbar单位
pos3=set(h,'Position', [0.92 0.1 0.022 0.85]);  %设置colorbar位置

上面的代码大部分都是重复的,只要弄懂第一个子图怎么画,后面的改变量就行。

(CMap是我下的颜色库,也可以用matlab自带的)

另外,多子图的colorbar的设置需要在最后单独写几行代码,可根据自己需求更改。

祝大家画出好看的图。

相关推荐
ytttr8737 小时前
混合动力汽车MATLAB建模实现方案
开发语言·matlab·汽车
尤利乌斯.X8 小时前
复杂网络仿真从入门到精通:0 学习路线
网络·学习·matlab·仿真·复杂网络
feifeigo12313 小时前
MATLAB的无线传感器网络(WSN)算法仿真
网络·算法·matlab
CUMT_DJ1 天前
matlab计算算法的运行时间
开发语言·算法·matlab
weixin_514221851 天前
FDTD与matlab、python耦合
python·学习·matlab·fdtd
Evand J3 天前
组合导航的MATLAB例程,二维平面上的CKF滤波,融合IMU和GNSS数据,仿真,观测为X和Y轴的坐标,附代码下载链接
开发语言·matlab·平面·imu·组合导航
CappuccinoRose3 天前
MATLAB学习文档(二十三)
matlab·信息可视化·数据挖掘·数据分析
民乐团扒谱机3 天前
【微实验】激光测径系列(四)关于硬件上的一些实验
计算机视觉·matlab·激光测径
MATLAB代码顾问4 天前
MATLAB计算标准径流指数(Standard Runoff Index,SRI)
数据结构·算法·matlab
listhi5204 天前
基于MATLAB的高斯混合模型(GMM)实现
开发语言·matlab