Question

我正在使用R（库mco）上的NSGA2包。

我的NSGA2代码需要永远运行，所以我想知道：

1）有没有办法限制解决方案值的精度（比方说，可能最多3位小数）而不是无限？

2）如何设置一个等式约束（在线所有似乎都是关于＆gt; =或＆lt; = than =）？不确定我是否做得对。

我的所有相关代码供参考，以便于跟踪：https://docs.google.com/document/d/1xj7OPng11EzLTTtWLdRWMm8zJ9f7q1wsx2nIHdh3RM4/edit?usp=sharing

此处转载的相关代码部分代码：

VTR = get.hist.quote(instrument = 'VTR',
        start="2010-01-01",  end = "2015-12-31",
        quote = c("AdjClose"),provider = "yahoo",
        compress = "d") 

ObjFun1 <- function (xh){
    f1 <- sum(HSVaR_P(merge(VTR, CMI, SPLS, KSS, DVN, MAT, LOE, KEL, COH, AXP), xh, 0.05, 2))
    tempt = merge(VTR, CMI, SPLS, KSS, DVN, MAT, LOE, KEL, COH, AXP)
    tempt2 = tempt[(nrow(tempt)-(2*N)):nrow(tempt),]
    for (i in 1:nrow(tempt2))
    {
        for (j in 1:ncol(tempt2))
        {
        if (is.na(tempt2[i,j])) 
        {
        tempt2[i,j] = 0
        }
        }
    }
    f2 <- ((-1)*abs(sum((xh*t(tempt2)))))
    c(f1=f1,f2=f2)
}

Constr <- function(xh){ 
    totwt <- (1-sum(-xh))
    totwt2 <- (sum(xh)-1)
    c(totwt,totwt2)
}


Solution1 <- nsga2(ObjFun1, n.projects, 2,
                  lower.bounds=rep(0,n.projects), upper.bounds=rep(1,n.projects),
                  popsize=n.solutions,  constraints = Constr, cdim=1,
                  generations=generations)

函数HSVaR_P返回矩阵（x，2 * 500,1）。

即使我设置代数= 1，代码似乎也没有运行。显然，代码中应该存在一些错误，但我不完全确定NSGA2算法的机制。

感谢。

R中的NSGA2遗传算法

0 个答案: