i im每周和每年的每月都会使用Weibull分布来模拟风力。由于风速/功率每周有差异,所以我想将它们全部结合在一起,然后复制结果,以查看特定风能在某些时间内可以产生多少能量。第一个模拟在第一周产生了第二周的第二个模拟结果,依此类推,目标的合并结果为一个月,然后在所需的几个月内复制结果。
first<-rweibull(7,3.28,2.14) #time series for days in first week
second<-rweibull(7,3.33,2.112)# time series for days in second week
combined<-replicate(2,c(first,second))#combined
[,1] [,2]
[1,] 2.3314627 2.3314627
[2,] 2.4219863 2.4219863
[3,] 2.2993156 2.2993156
[4,] 0.9576821 0.9576821
[5,] 0.8065639 0.8065639
[6,] 2.2820711 2.2820711
[7,] 0.6315118 0.6315118
[8,] 1.8152144 1.8152144
[9,] 0.9296449 0.9296449
[10,] 1.6707807 1.6707807
[11,] 1.7307920 1.7307920
[12,] 2.2075004 2.2075004
[13,] 2.0368642 2.0368642
[14,] 2.5320678 2.5320678
现在的问题是,当我复制时,我得到相同的结果,但是结果应该不同。是否有任何简单而快速的解决方案。
,而不是将变量传递给replicate
。您可以做到这一点。
代码
fun <- function( ){
first <- rweibull(7,3.28,2.14)
second <- rweibull(7,3.33,2.112)
}
replicate( 2, fun() )