我有电力系统频率的原始数据。 86 400个数字。
frequency=a$Ist_Frq
plot.ts(frequency, main="System frequency [Hz]", xlab="Time [s]")
参见示例:
现在,我必须确定四分之一小时的时间间隔。
frequency=ts(a$Ist_Frq, start=1, frequency=900)
[quarter-hour time interval][2]
我的问题是:
有没有办法如何确定每个季度小时的标准偏差?
感谢您的回答。
答案 0 :(得分:0)
这个问题可能有几个解决方案:这里有一个
#some data
x <- rnorm(10000)
#identify quarter hour segments
y <- rep(1:ceiling(length(x)/(15 * 60)), each = 15 * 60)[1:length(x)]
#use tapply to find sd of x for every value of y
tapply(x, y, sd)
nb最后一个值可能基于少于900个值