当我尝试为10000次迭代运行以下代码时,我得到以下错误.代表中的错误(G1 [,2],G1 [,3]):无效的'times'参数.所以不知道如何更改代码来修复该错误.基本上只是想创建发电机性能的时间序列,使用时间失败的等式和一年中8736小时的修复时间,这样我就可以在发电机运行的时间内以小时为单位.起始条件是发电机在第一个小时运行.肯定有一个更优雅的解决方案来模拟这个我只是无法找到它.任何评论或帮助将不胜感激.
MTTF<-2940 # MEDIUM TIME TO FAIL(hours)
MTTR<-60 # MEDIUM TIME TO REPAIR (hours)
TTF<--MTTF*log(runif(100))# equation for Time to fail
TTR<--MTTR*log(runif(100))# equation for Time to repair
mix<-rep(0,length(TTF)+length(TTR))
sw<-rep(0,length(TTF)+length(TTR))
for(i in 1:length(TTF)){
mix[2*i-1]<-TTF[i]
sw[2*i-1]<-1
mix[2*i]<-TTR[i]
}
cmix<-cumsum(mix)
ccmix<-cbind(cmix[1:which(cmix>8736)],sw[1:which(cmix>8736)])
ccmix[dim(ccmix)[1],1]<-8736
G1<-round(ccmix)
# transform binary values
G1[G1 == 1] <- 12 # is the capacity of the generator
G1 <- cbind(G1, c(G1[1,1], diff(G1[,1])))
a1 <- rep(G1[,2], G1[,3]) ## GENERATING 8736 Values
因此,当打开时,所需的输出为8736值12,关闭时为0
我有以下问题.我想计算小于或等于零的值的出现次数.在以下数据中的示例我有3次出现1(0,0,0),2(-1,-2),3(0,0).R中是否有任何函数构建来计算连续出现的次数.
a <- c(2,4,5,3,2,4,7,0,0,0,4,3,2,-1,-2,3,2,0,0,4)
我尝试在Anaconda（Python 3.6.8）上设置pyFMI
安装了pyFMI站点上列出的所有必需软件包。fmu加载没有问题，但是当我尝试模拟fmu时出现错误：
Could not find cannot import name 'radau5'
Could not find cannot import name 'dopri5'
Could not find cannot import name 'rodas'
Could not find cannot import name 'odassl'
Could not find ODEPACK functions.
Could not find RADAR5
Could not find GLIMDA.
Traceback (most recent call last):
File "assimulo\solvers\../lib/sundials_callbacks_ida_cvode.pxi", line 240, in assimulo.solvers.sundials.cv_jac
File "C:\Users\d60378\AppData\Local\Continuum\anaconda3\lib\site-packages\pyfmi\simulation\assimulo_interface.py", line 733, in j
A = self._model._get_A(add_diag=True, output_matrix=self._A)
File "src\pyfmi\fmi.pyx", line 6041, in pyfmi.fmi.FMUModelBase2._get_A
File "src\pyfmi\fmi.pyx", line 7592, in pyfmi.fmi.FMUModelME2._get_directional_proxy
File "src\pyfmi\fmi.pyx", line …