tic?%?計(jì)時(shí)器
%%?清空環(huán)境，準(zhǔn)備數(shù)據(jù)
close?all
clear
clc
format?compact
%?載入測(cè)試數(shù)據(jù)wine其中包含的數(shù)據(jù)類(lèi)別數(shù)為3；wine:178*13的矩陣wine_labes:178*1的列向量
load?wine
%?選定訓(xùn)練集和測(cè)試集
%?將第一類(lèi)的1-30第二類(lèi)的60-95第三類(lèi)的131-153做為訓(xùn)練集
train_wine?=?[wine（1:30:）;wine（60:95:）;wine（131:153:）];
%?相應(yīng)的訓(xùn)練集的標(biāo)簽也要分離出來(lái)
train_wine_labels?=?[wine_labels（1:30）;wine_labels（60:95）;wine_labels（131:153）];
%?將第一類(lèi)的31-59第二類(lèi)的96-130第三類(lèi)的154-178做為測(cè)試集
test_wine?=?[wine（31:59:）;wine（96:130:）;wine（154:178:）];
%?相應(yīng)的測(cè)試集的標(biāo)簽也要分離出來(lái)
test_wine_labels?=?[wine_labels（31:59）;wine_labels（96:130）;wine_labels（154:178）];
%?數(shù)據(jù)預(yù)處理
%?數(shù)據(jù)預(yù)處理將訓(xùn)練集和測(cè)試集歸一化到[01]區(qū)間
[mtrainntrain]?=?size（train_wine）;
[mtestntest]?=?size（test_wine）;
dataset?=?[train_wine;test_wine];

[dataset_scaleps]?=?mapminmax（dataset‘01）;
dataset_scale?=?dataset_scale‘;

train_wine?=?dataset_scale（1:mtrain:）;
test_wine?=?dataset_scale（?（mtrain+1）:（mtrain+mtest）:?）;
%%?CS-SVM
time=20;
n=20;?%?n為巢穴數(shù)量
pa=0.25;?%?被宿主發(fā)現(xiàn)的概率
dim?=?2;?%?需要尋優(yōu)的參數(shù)個(gè)數(shù)
Lb=[0.010.01];?%?設(shè)置參數(shù)下界
Ub=[100100];?%?設(shè)置參數(shù)上界

%?隨機(jī)初始化巢穴
nest=zeros（ndim）;
for?i=1:n?%?遍歷每個(gè)巢穴
????nest（i:）=Lb+（Ub-Lb）.*rand（size（Lb））;?%?對(duì)每個(gè)巢穴，隨機(jī)初始化參數(shù)
end

fitness=ones（1n）;?%?目標(biāo)函數(shù)值初始化
[fminbestnestnestfitness]=get_best_nest_svm（nestnestfitnesstrain_wine_labelstrain_winetest_wine_labelstest_wine）;?%?找出當(dāng)前最佳巢穴和參數(shù)

%%?迭代開(kāi)始
for?t=1:time
????new_nest=get_cuckoos（nestbestnestLbUb）;?%?保留當(dāng)前最優(yōu)解，尋找新巢穴
????[~~nestfitness]=get_best_nest_svm（nestnew_nestfitnesstrain_wine_labelstrain_winetest_wine_labelstest_wine）;?%?找出當(dāng)前最佳巢穴和參數(shù)
????new_nest=empty_nests（nestLbUbpa）;?%?發(fā)現(xiàn)并更新劣質(zhì)巢穴
????%?找出當(dāng)前最佳巢穴和參數(shù)
????[fnewbestnestfitness]=get_best_nest_svm（nestnew_nestfitnesstrain_wine_labelstrain_winetest_wine_labelstest_wine）;?
????if?fnew????????fmin=fnew;
????????bestnest=best?;
????end
end
%%?打印參數(shù)選擇結(jié)果
bestobjfun=fmin;
bestc=bestnest（1）;
bestg=bestnest（2）;
disp（‘打印參數(shù)選擇結(jié)果‘）;
str=sprintf（‘Best?c?=?%g，Best?g?=?%g‘bestcbestg）;
disp（str）
%%?利用最佳的參數(shù)進(jìn)行SVM網(wǎng)絡(luò)訓(xùn)練
cmd_gwosvm?=?[‘-c?‘num2str（bestc）‘?-g?‘num2str（bestg）];
model_gwosvm?=?svmtrain（train_wine_labelstrain_winecmd_gwosvm）;
%%?SVM網(wǎng)絡(luò)預(yù)測(cè)
[predict_labelaccuracy]?=?svmpredict（test_wine_labelstest_winemodel_gwosvm）;
%?打印測(cè)試集分類(lèi)準(zhǔn)確率
total?=?length（test_wine_labels）;
right?=?sum（predict_label?==?test_wine_labels）;
disp（‘打印測(cè)試集分類(lèi)準(zhǔn)確率‘）;
str?=?sprintf（?‘Accuracy?=?%g%%?（%d/%d）‘a(chǎn)ccuracy（1）righttotal）;
disp（str）;
%%?結(jié)果分析
%?測(cè)試集的實(shí)際分類(lèi)和預(yù)測(cè)分類(lèi)圖
figure;
hold?on;
plot（test_wine_labels‘o‘）;
plot（predict_label‘r*‘）;
xlabel（‘測(cè)試集樣本‘‘FontSize‘12）;
ylabel（‘類(lèi)別標(biāo)簽‘‘FontSize‘12）;
legend（‘實(shí)際測(cè)試集分類(lèi)‘‘預(yù)測(cè)測(cè)試集分類(lèi)‘）;
title（‘測(cè)試集的實(shí)際分類(lèi)和預(yù)測(cè)分類(lèi)圖‘‘FontSize‘12）;
grid?on
snapnow
%%?顯示程序運(yùn)行時(shí)間
toc

?屬性????????????大小?????日期????時(shí)間???名稱(chēng)
-----------?---------??----------?-----??----
?????目錄???????????0??2016-09-05?17:22??CS-SVM\
?????文件????????3051??2016-09-05?17:22??CS-SVM\CS_SVM_exmp.m
?????文件????????2716??2016-09-05?17:00??CS-SVM\CS_SVR_exmp.m
?????文件????????6483??2016-09-05?16:52??CS-SVM\cuckoo_ori_with_chinese_note.m
?????文件?????????999??2016-09-05?16:49??CS-SVM\empty_nests.m
?????文件?????????426??2016-09-05?16:49??CS-SVM\fobj.m
?????文件?????????676??2016-09-05?17:22??CS-SVM\get_best_nest.m
?????文件?????????698??2016-09-05?17:21??CS-SVM\get_best_nest_svm.m
?????文件????????1629??2016-09-05?16:49??CS-SVM\get_cuckoos.m
?????文件??????????70??2016-09-05?16:49??CS-SVM\mymae.m
?????文件??????????85??2016-09-05?16:49??CS-SVM\mymape.m
?????文件??????????74??2016-09-05?16:49??CS-SVM\mymse.m
?????文件?????????402??2016-09-05?17:13??CS-SVM\objfun.m
?????文件?????????376??2016-09-05?16:49??CS-SVM\simplebounds.m
?????文件???????20168??2016-09-05?17:17??CS-SVM\wine.mat
?????文件?????????447??2016-09-05?16:49??CS-SVM\wndspd.mat
?????目錄???????????0??2016-09-05?16:45??PSO-SVM\
?????文件?????????426??2016-09-05?16:14??PSO-SVM\fobj.m
?????文件?????????108??2016-09-05?15:04??PSO-SVM\fun.m
?????文件??????????70??2016-09-05?16:14??PSO-SVM\mymae.m
?????文件??????????85??2016-09-05?16:14??PSO-SVM\mymape.m
?????文件??????????74??2016-09-05?16:14??PSO-SVM\mymse.m
?????文件?????????402??2016-09-05?17:12??PSO-SVM\objfun.m
?????文件????????4106??2016-09-05?16:43??PSO-SVM\PSO_SVM_exmp.m
?????文件????????3814??2016-09-05?17:00??PSO-SVM\PSO_SVR_exmp.m
?????文件????????2276??2016-09-05?16:11??PSO-SVM\pso_test.m
?????文件???????20168??2016-09-05?16:37??PSO-SVM\wine.mat
?????文件?????????447??2016-09-05?16:14??PSO-SVM\wndspd.mat