Matlab 绘制相关系数像素图

1.代码 

%变量
para_all=[input_train_2017,input_train_2018,input_train_2019...
   ;output_train_2017,output_train_2018,output_train_2019 ];
%求相关系数矩阵
corr1=corrcoef(para_all');
%相关系数矩阵是对称的,令下三角矩阵保持不变,上三角置为0
corr2=tril(corr1)
%令上三角置为NaN
corr2(corr2==0)=NaN;

img=imagesc(corr3); 
%NaN默认显示为数值最小的颜色,这里让它显示为白色
set(img,'alphadata',~isnan(corr3));

colormap(jet)
c=colorbar;
caxis([-1,1])
set(c,'YTick',(-1:0.4:1)); 

xlabel('Variables');
ylabel('Variables');
%x轴刻度标签
set(gca,'XTick',(1:9),'xticklabel',[string('Year'),string('Doy'),string('Hour'),...
string('P'),string('T'),string('RH'),string('ZTD'),string('HTG'),string('Rain')],'FontSize',10); 
%Y轴刻度标签
set(gca,'YTick',(1:9),'yticklabel',[string('Year'),string('Doy'),string('Hour'),...
string('P'),string('T'),string('RH'),string('ZTD'),string('HTG'),string('Rain')],'FontSize',10); 
%图片大小
set(gcf,'unit','centimeters','position',[6,4,12,9])
%保存为分辨率300的tiff图片
sfname='corr';
print(gcf,sfname,'-r300','-dtiff');

2.结果

 

 

评论 4
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值