Определить и построить функцию PieceWise в R - PullRequest
0 голосов
/ 11 апреля 2020

Я пытался определить эту функцию с помощью RStudio, но у меня много проблем, и я не могу понять, в чем дело

x <- seq(-10, 10, 0.001)
fx <- function()
if(x < 0 ) {fx<- sin(x)
} 
 else if (x >= 0 && x < 2) 
   {
   fx<- x^2
 } 
  else if(x >= 2 )
    { 
    fx<- 4*exp(x-2)
  }
}
plot(x, fx)

Как я могу это сделать? Я не могу найти ничего полезного на inte rnet

Ответы [ 3 ]

2 голосов
/ 11 апреля 2020

Ваша функция не должна использовать if для векторов. Лучше всего для нас ifelse, или просто используйте векторные возможности R:

fx <- function(x){
    f <- NULL
    f[x<0] <- sin(x[x<0])
    f[x>=0 & x<2] <- x[x>=0 & x<2]^2
    f[x>=2] <- 4*exp(x[x>=2] - 2)
    f
}

x <- seq(-10, 3, 0.01)
plot(x, fx(x), type="l", las=1)

enter image description here

0 голосов
/ 11 апреля 2020

Решение с использованием Vectorize:

x <- seq(-2, 1, 0.001)
f <- function(x) {
  if(x<0) {
    fx <- sin(x)
  } else if (x >= 0 && x < 2) {
    fx <- x^2
  } else if(x >= 2 ) { 
    fx <- 4*exp(x-2)
  }
  fx
}
f <- Vectorize(f)
plot(x, f(x))

enter image description here

0 голосов
/ 11 апреля 2020

Попробуйте это:

x <- seq(-10, 10, 0.001)

fx <- function(x) {
  if (x < 0) {
    return(sin(x))
  } else if (x >= 0 & x < 2) {
    return(x^2)
  } else {
    return(4*exp(x-2))
  }
}

plot(x, sapply(x, fx))

Создано в 2020-04-11 пакетом Представить (v0.3.0)

Добро пожаловать на сайт PullRequest, где вы можете задавать вопросы и получать ответы от других членов сообщества.
...