求一段时序数据的波峰波谷个数

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

如何求一段序列数字的波峰波谷个数
例如a<-c(2,3,6,5,7,9,2,3,6,5,7,2,3,6,5,7,9,10)

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏2 回帖

a<-c(2,3,6,5,7,9,2,3,6,5,7,2,3,6,5,7,9,10)
peak <- function(x) {
    n <- length(x)
    if(n < 3) {
        return(0)
    }
    count <- 0
    for(i in 2:(n - 1)) {
        if(x[i] > x[i - 1] && x[i] > x[i + 1]) {
            count <- count + 1
        }
    }
    return(count)
}
peak(a)
peak(-a) ##波谷

1. a<-c(2,3,6,5,4,7,9,2,3,6,5,7,2,3,6,5,7,9,10)
   
2. findpeak_boo <- function(x){
   
3.   r <- rle(x)
   
4.   pks <- which(rep(x = diff(sign(diff(c(-Inf, r$values, -Inf)))) == -2,
   
5.                     times = r$lengths))
   
6.   boos <- which(rep(x = diff(sign(diff(c(-Inf, r$values, -Inf)))) == 2,
   
7.                     times = r$lengths))
   
8.   n_pks <- length(pks)
   
9.   n_boos <- length(boos)
   
10.   return(list(
    
11.     "pks"=pks,
    
12.     "boos"=boos,
    
13.     "n_pks"=n_pks,
    
14.     "n_boos"=n_boos
    
15.   ))
    
16. }
    
17. findpeak_boo(a)

复制代码

1. > findpeak_boo(a)
   
2. $pks
   
3. [1]  3  7 10 12 15 19
   
4. 
   
5. $boos
   
6. [1]  5  8 11 13 16
   
7. 
   
8. $n_pks
   
9. [1] 6
   
10. 
    
11. $n_boos
    
12. [1] 5

复制代码

像这种

菊花冰糖水 发表于 2018-3-5 18:07
像这种

Temp<-mydata$齿轮箱油池温度_平均值
findpeak_valley1 <- function(x){
  r <- rle(x)
  pks_position <-which(rep(x = diff(sign(diff(c(-Inf, r$values, -Inf)))) == -2, times = r$lengths))
  valley_position <-which(rep(x = diff(sign(diff(c(-Inf, r$values, -Inf)))) == 2, times = r$lengths))
  pks<-x[pks_position]
  valley<-x[valley_position]
  n_pks <- length(rle(pks)$value)
  n_valley <- length(rle(valley)$value)
  return(list(
    "pks_position"=pks_position,
    "valley_position"=valley_position,
    "pks"=pks,
    "valley"=valley,
    "n_pks"=n_pks,
    "n_valley"=n_valley
  ))
}
result<-findpeak_valley1(Temp)
paste('The pks number is %d',result$n_pks)
paste('The valley number is %d',result$n_valley)

你好，我中间稍微改动了点程序，将波峰波谷中重复的点去掉，然后计算一段油温数据的波峰波谷，最后得出的波峰和波谷的数量是不同，这是为什么呢？理论上讲，一段序列数据，其波峰和波谷的数量不应该是相同或者差1吗？可是我的差很多

yan_huili 发表于 2018-3-6 10:15
Temp

程序没什么大问题。没有数据我也不清楚哪里出的错。 找波峰波谷的逻辑感觉也没啥问题。你可能需要多花点功夫具体问题具体分析了。 有时候能用的代码并不一定在所有的数据上都适用。

yan_huili 发表于 2018-3-6 10:15
Temp

程序没什么大问题。没有数据我也不清楚哪里出的错。 找波峰波谷的逻辑感觉也没啥问题。你可能需要多花点功夫具体问题具体分析了。 有时候能用的代码并不一定在所有的数据上都适用。

yan_huili 发表于 2018-3-6 10:15
Temp

程序没什么大问题。没有数据我也不清楚哪里出的错。 找波峰波谷的逻辑感觉也没啥问题。你可能需要多花点功夫具体问题具体分析了。 有时候能用的代码并不一定在所有的数据上都适用。

finds<-function(x){
  r<-rle(x)
   pks<-which(diff(sign(diff(c(-Inf,r$values,-Inf))))==-2)
  boos<-which(diff(sign(diff(c(-Inf,r$values,-Inf))))==2)
  n_pks<-length(pks)
  n_boos<-length(boos)
  return(list(
    "pks"=pks,
    "boos"=boos,
    "n_pks"=n_pks,
    "n_boos"=n_boos
    ))
}