矩阵的单元交换算法

Question

感谢this answer 的一些大力帮助，我写了一些应该的代码

1. 在矩阵中随机选择一个单元格matrnum
2. 将单元格从零更改为一，反之亦然
3. 在对角三角形的相应单元格中进行相反的更改
4. 检查条件2，看变化的新矩阵是否低于原矩阵
5. 如果条件较低则返回新矩阵，如果不返回原始矩阵
6. 每次都使用新的rnum 重复此操作，直到满足整体条件1

mat 在哪里：

mat <- matrix(rep(0, 100), nrow = 10)
mat[lower.tri(mat, diag = TRUE)] <- 1

---

while( condition1 > threshold) {

rnum <- sample(1:100, 1)
if(mat[rnum] > 0){
mat2 <- mat
mat2[rnum] <- 0
mat2[length(mat2) - rnum + 1] <- 1
return(mat2)
}else{
mat2 <- mat
mat2[rnum] <- 1
mat2[length(mat2) - rnum + 1] <- 0
return(mat2)
        }   
if(condition2(mat2) < 
   condition2(mat) ) {
mat <- mat2
}else{
mat
    }
}

这里的条件不是关键，为了论证condition1和condition2，它可能是顶行的总和，即sum(mat[,1])，但我收到以下错误：

错误：没有函数可以返回，跳转到顶层

谁能指出我的缺陷？

Answer 1

我发现的错误之一是在您的 else 块中，未定义 mat2。因此，如果随机数在第一次迭代中将我带到 else 块，代码错误：

Error in mat2[rnum] <- 1 : object 'mat2' not found

您提到的错误正在发生，因为 while 循环只是一个语句，您的代码中没有封闭函数，因此 return(mat2) 语句不正确。如果要从函数返回，则只能使用 return。在您的情况下，没有功能，因此出现错误