在R中已經有現成的粒子群演算法套件,可以直接網路下載安裝後使用,
避免重複造輪子之外也可以參考別人的寫法。
避免重複造輪子之外也可以參考別人的寫法。
首先要安裝和載入PSO package
install.packages("pso")
library("pso")
之後就可以開心使用了,首先設定目標函數:
fit <- function(x) 100 * (x[2] - x[1]^2)^2 + (1 - x[1])^2
開始一個簡單的最佳化吧!
install.packages("pso")
library("pso")
fit <- function(x) 100 * (x[2] - x[1]^2)^2 + (1 - x[1])^2