数据分析-基于R（潘文超）第14至18章 程序代码汇总.docxVIP

目录 TOC \o 1-3 \h \z \u 第十四章 1 14.2 1 例一 1 例二 1 14.3 2 例一 2 例二 2 14.4 2 例一 2 例二 2 第十五章 2 15.1 2 15.2 3 例一 3 例二 4 第十四章 14.2 例一 install. packages(mcga) library(mcga) f-function(x){ return( - 25 + x^2) } m - mcga( popsize= 200, chsize=5, minval=0.0, maxval 9, maxiter=2500, crossprob=1. 0, mutateprob=0.01, evalFunc=f) cat( Best chromosome: \n ) print(m$ population[1, ]) cat(Cost:,m$ costs[1], \n) 例二 install.packages(mcga) library(mcga) V=rnorm(4) H=rnorm(4) W=100 f - function(x) { return(W* ((H[1]-x[1])^2+ (V[1]- x[2])^2)^0.5+W* ((H[2]-x[1])^2+ (V[2] - x[2])^2)^0.5+W* ((H[3] - x[1])~2+ (V[3]- x[2])~2)~0.5+W* ((H[4]-x[1])~2+ (V[4]-x[2])~2)~0.5) } m- mcga( popsize= 200, chsize= 5, minval=0.0, maxval = 9, maxiter = 1000, crossprob= 1.0, mutateprob= 0.01, evalFunc= f) cat( Best chromosome: \n ) print(m $ population[1, ]) cat( Cost:,m$ costs[1],\n) 14.3 例一 install. packages( pso) library( pso) f - function(x) - cos(x[1])￥cos(x[2]) * exp( - ((x[1]-pi)^2 + (x[21-pi)^2)) psopt. im( rep(NA, 2), f, lower= - 100, upper= 100) 例二 install. packages( pso ) l ibrary( pso ) V= rnorm(4) H= rnorm(4) W= 100 o= psopt im( rep(NA, 2)，function(x)W* ((H[1]- x[1])^2+ (V[1] - x[2])^2)^0.5 + W* ((H[2]- x[1])~2 + (V[2]- x[2])~2)^0.5+W* ((H[3]- x[1])~2+ (V[3]- x[2])~2)^0.5+W* ((H[4]-x[1])~2+ (V[4]- x[2])~2)^0.5,lower=c( - 100, - 100)，upper = c( 100, 100)，control = list(s= 50, maxit= 1000, trace= 1, REPORT= 1, trace. stats = TRUE) ) plot(o$ stats$ it,o$ stats$ error,1og= y,xlab=It, ylab = Cost) 14.4 例一 install. packages( ABCoptim) 1 ibrary( ABCoptim) fun - function(x){ - cos(x[1])￥cos(x[2])￥exp(-((x[1]- pi)^2+ (x[2]- pi)~2))} abc_ optim( rep(0, 2), fun,1b= - 20, ub= 20, maxCycle = 1000, criter = 1000) 例二 install. packages( ABCopt inm) 1 ibrary( ABCoptim) V= rnorm(4) H= rnorm(4 ) W= 100 fun - function(x){W* ((H[1]- x[1])~2+ (V[1]- x[2])^2)^0.5+ W* ((H[2]- x[1])~2 + (V[2]- x[2])~2)^0.5+W* ((H[3]-x[1])~2+ (V[3]- x[2])~2)^0.5+W* ((H[4]- x[