clear?all
clc

%%?II.?導入數據
load?spectra_data.mat

%%
%?1.?隨機產生訓練集和測試集
n?=?randperm（size（NIR1））;

%%
%?2.?訓練集——50個樣本
p_train?=?NIR（n（1:50）:）;
t_train?=?octane（n（1:50）:）;

%%
%?3.?測試集——10個樣本
p_test?=?NIR（n（51:end）:）;
t_test?=?octane（n（51:end）:）;

%%?III.?數據歸一化
%%
%?1.?訓練集
[pn_traininputps]?=?mapminmax（p_train‘）;
pn_train?=?pn_train‘;
pn_test?=?mapminmax（‘apply‘p_test‘inputps）;
pn_test?=?pn_test‘;

%%
%?2.?測試集
[tn_trainoutputps]?=?mapminmax（t_train‘）;
tn_train?=?tn_train‘;
tn_test?=?mapminmax（‘apply‘t_test‘outputps）;
tn_test?=?tn_test‘;
%%
%[bestCVmsebestcbestgpso_option]?=?psoSVMcgForRegress（t_trainp_trainpso_option）
pso_option?=?struct（‘c1‘1.5‘c2‘1.7‘maxgen‘200‘sizepop‘20?...
????????????????????‘k‘0.6‘wV‘1‘wP‘1‘v‘5?...
?????????????????????‘popcmax‘10^2‘popcmin‘10^（-1）‘popgmax‘10^3‘popgmin‘10^（-2））;
?????????????????%?c1:初始為1.5pso參數局部搜索能力
%?c2:初始為1.7pso參數全局搜索能力
%?maxgen:初始為200最大進化數量
%?sizepop:初始為20種群最大