2017-02-22 76 views
-1

我模擬每週和每年一個月的威布爾分佈風電。由於風速/功率每週都有差異,我想將它們全部結合起來,然後複製結果以查看特定風力位置在特定時間內可產生多少能量。第一次模擬產生第一週的結果,第二週的第二次模擬等等,目標有一個月的合併結果,然後將結果複製爲期望的月數。合併兩個或多個模擬

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  

現在的問題是,當我複製我得到相同的結果,但結果應該是不同的。有沒有簡單而快速的解決方案。

+0

你能解釋一下多一點?預期的輸出是什麼樣的?我沒有看到你正在通過任何功能來複制,當你說複合時,它的混淆和使用複製......'cbind'應該結合起來 – user5249203

+0

我是模擬每週和每個月的風能與威布爾分佈年。因爲風速/功率每週都有差異,所以我希望將它們全部結合起來,然後複製結果以查看在特定時間內特定風力位置可產生多少能量 – kelamahim

+0

,或者您應該將其傳遞給像'replicate(2,c (rweibull(7,3.28,2.14),rweibull(7,3.33,2.112)))'。而不是儲存在一個變量和後來通過它 – user5249203

回答

1

而不是將變量傳遞到replicate你可以做到這一點。

代碼

fun <- function(){ 
    first <- rweibull(7,3.28,2.14) 
    second <- rweibull(7,3.33,2.112) 
} 

replicate(2, fun())