91av视频/亚洲h视频/操亚洲美女/外国一级黄色毛片 - 国产三级三级三级三级

  • 大小: 6KB
    文件類型: .m
    金幣: 1
    下載: 0 次
    發(fā)布日期: 2021-06-14
  • 語言: Matlab
  • 標(biāo)簽:

資源簡介

給出了最小二乘支持向量機在matlab中的仿真程序

資源截圖

代碼片段和文件信息

function?[bestCVaccuarcybestcbestgpso_option]?=?psoLSSVMcgForClass(trainsettrainset_labelYcpso_option)
%?psoLSSVMcgForClass?
%by?Tangxiaobiao??QQ?444646122?blog.sina.com.cn/lssvm
%?2010.05.31

%%?參數(shù)初始化
if?nargin?==?3
????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));
end
%?c1:初始為1.5pso參數(shù)局部搜索能力
%?c2:初始為1.7pso參數(shù)全局搜索能力
%?maxgen:初始為200最大進化數(shù)量
%?sizepop:初始為20種群最大數(shù)量
%?k:初始為0.6(k?belongs?to?[0.11.0])速率和x的關(guān)系(V?=?kX)
%?wV:初始為1(wV?best?belongs?to?[0.81.2])速率更新公式中速度前面的彈性系數(shù)
%?wP:初始為1種群更新公式中速度前面的彈性系數(shù)
%?v:初始為3SVM?Cross?Validation參數(shù)
%?popcmax:初始為100SVM?參數(shù)c的變化的最大值.
%?popcmin:初始為0.1SVM?參數(shù)c的變化的最小值.
%?popgmax:初始為1000SVM?參數(shù)g的變化的最大值.
%?popgmin:初始為0.01SVM?參數(shù)c的變化的最小值.
Yc=Yc;
Vcmax?=?pso_option.k*pso_option.popcmax;
Vcmin?=?-Vcmax?;
Vgmax?=?pso_option.k*pso_option.popgmax;
Vgmin?=?-Vgmax?;

eps?=?10^(-3);
[Yccodebookold_codebook]?=?code(trainset_label‘code_MOC‘);
%%?產(chǎn)生初始粒子和速度
for?i=1:pso_option.sizepop
????
????%?隨機產(chǎn)生種群和速度
????pop(i1)?=?(pso_option.popcmax-pso_option.popcmin)*rand+pso_option.popcmin;
????pop(i2)?=?(pso_option.popgmax-pso_option.popgmin)*rand+pso_option.popgmin;
????V(i1)=Vcmax*rands(1);
????V(i2)=Vgmax*rands(1);
????
????%?計算初始適應(yīng)度
?????gam=pop(i1);sig2=pop(i2);
?????model=initlssvm(trainsetYc‘c‘gamsig2‘RBF_kernel‘);
?????model=trainlssvm(model);
?????Yd0=simlssvm(modeltrainset);
?????predict_label?=?code(Yd0old_codebook[]codebook);%解碼分類結(jié)果
?????testnum=size(trainset_label1);
?????right=sum(trainset_label==predict_label);
?????n?=?sum(sum(trainset_label~=predict_label));
?????fitness(i)?=?(1-n/prod(size(trainset_label)))*100;
????fitness(i)?=?-fitness(i);
end

%?找極值和極值點
[global_fitness?bestindex]=min(fitness);?%?全局極值
local_fitness=fitness;???%?個體極值初始化

global_x=pop(bestindex:);???%?全局極值點
local_x=pop;????%?個體極值點初始化

%?每一代種群的平均適應(yīng)度
avgfitness_gen?=?zeros(1pso_option.maxgen);

%%?迭代尋優(yōu)
for?i=1:pso_option.maxgen
????
????for?j=1:pso_option.sizepop
????????
????????%速度更新
????????V(j:)?=?pso_option.wV*V(j:)?+?pso_option.c1*rand*(local_x(j:)?-?pop(j:))?+?pso_option.c2*rand*(global_x?-?pop(j:));
????????if?V(j1)?>?Vcmax
????????????V(j1)?=?Vcmax;
????????end
????????if?V(j1)?????????????V(j1)?=?Vcmin;
????????end
????????if?V(j2)?>?Vgmax
????????????V(j2)?=?Vgmax;
????????end
????????if?V(j2)?????????????V(j2)?=?Vgmin;
????????end
????????
????????%種群更新
????????pop(j:)=pop(j:)?+?pso_option.wP*V(j:);
????????if?pop(j1)?>?pso_option.popcmax
????????????pop(j1)?=?pso_option.popcmax;
????????end
????????if?pop(j1)????????????

評論

共有 條評論