论文变量分布图的绘制

   在近红外变量选择的论文中，需要对选择的波段进行可视化，以下介绍两种变量分布图的绘制。

MATLAB绘制点图

直接上代码

%数据输入，一行多列
load('cor1_huatu.mat')
load('RSR1_huatu.mat')
load('SSC1_huatu.mat')
load('sopsrl1_huatu.mat')data = {};
data{1} = cor1_huatu';
data{2} = RSR1_huatu';
data{3} = SSC1_huatu';
data{4} = sopsrl1_huatu';linewidth = 1;
linelength = 0.06;
% color = 'rcbrmcmky';color(i+1)
% c = colorcube(32);figure(1);
len = length(data);
for i = 1:lencc = 0.5- i*0.1;%0.2 * i ;%+ 0.2;dx = cell2mat(data(i))+1000;dy = ones(1, length(dx))-0.4;plot([dx;dx], [dy-i/5; dy-i/5-linelength],'color',[cc cc+0.1 cc+0.4], 'linewidth', linewidth);hold on;
end% 光谱图如下
plot(1001:1800,X);

画出来的光谱图如下：

最上面为全谱区的近红外光谱图，底下为四种算法筛选出的变量数在全谱区的变量分布。

MATLAB绘制线图

代码如下：

% 数据输入，一行多列
load('cor1_huatu.mat')
load('RSR1_huatu.mat')
load('SSC1_huatu.mat')
load('sopsrl1_huatu.mat')load('cor2_huatu.mat')
load('RSR2_huatu.mat')
load('SSC2_huatu.mat')
load('sopsrl2_huatu.mat')load('cor3_huatu.mat')
load('RSR3_huatu.mat')
load('SSC3_huatu.mat')
load('sopsrl3_huatu.mat')data = {};
data{1} = cor1_huatu';
data{2} = RSR1_huatu';
data{3} = SSC1_huatu';
data{4} = sopsrl1_huatu';data{5} = cor2_huatu';
data{6} = RSR2_huatu';
data{7} = SSC2_huatu';
data{8} = sopsrl2_huatu;data{9} = cor3_huatu';
data{10} = RSR3_huatu';
data{11} = SSC3_huatu';
data{12} = sopsrl3_huatu;
% data{5} = randi([0,800], 1 , 600);
% data{6} = randi([0,800], 1 , 500);
% data{7} = randi([0,800], 1 , 80);
% data{8} = randi([0,800], 1 , 80);name = {'A''B''C''D''E''F''G''H''I''J''K''L'};
high = max(max(X))+0.1;
low = min(min(X))-0.2;
linewidth = 1;
% color = 'rykbgmcbrgbykcm';
color = [0.5 0.5 0.8 ];
% subplot(ceil(len/4), i*(len<4)+4*(len>3),i);
len = length(data);
num = 4;
order = [1     5     9     2     6    10     3     7    11     4     8    12];
% box offfor i = 1:lensubplot(ceil(len/num), len*(len<num)+(num)*(len>num-1),i);dx = cell2mat(data(i)) + 1000;dy = zeros(1, length(dx));plot([dx;dx], [dy+high; dy+low], 'color', color, 'linewidth', linewidth);set(gca,'YLim',[0,2.6]);title(name(i));hold on;plot(1001:1800,X); xlabel('Wavelength/nm');ylabel('Absorbance');
end% 光谱图如下