iPLS算法

clc
clear
close all

data=xlsread("F:\2023.6.3.实验测试代码\ＣＧＬ＿ＮＩＲ\CGL_NIR.xlsx");
load normalize_maxain.mat
normalize_maxain=Columnset;
nir=normalize_maxain(:,1:117);
nir_label=data(:,118);

plot(nir')

[Xcal,Xcal_1,Xtest,Xtest_1,cal_indx,test_indx] =kennard_stone_example(data,nir,162);
train_label = Xcal_1(:,118);
train = Xcal;
test_label = Xtest_1(:,118);
test = Xtest;

xaxis = 900:801/117:1700;

Model=ipls(train,train_label,10,'mean',1,xaxis,'syst123',5); 

plsrmse(Model,0)

iplsplot(Model,'intlabel',4)

%iplsplot(Model,'wavlabel',4)

plspvsm(Model,9,1)