Программирование прямоугольной волны и пилообразной волны в R

Question

Я пытаюсь построить функцию в R, представляющую квадратные и пилообразные волны со 100 Гц в диапазоне [0,1]. Я попытался это:

squarewave <- function (t) { 

    # 0.01 == 100Hz(=2Pi) -> 1 Period of the Squarewave 
    # 0.005 == Pi   -> Half Period of the Squarewave 
    # if t smaller than a half period -> 1 
    # if t greater or equal than half a period -> 0 

    if ((t %% 0.01) < 0.005) 
    return (1) 
    else if ((t %% 0.01) >= 0.005) 
    return (0) 

} 

Когда я пытаюсь построить эту функцию:

plot(squarewave) 

Я получаю следующее сообщение об ошибке:

> plot(squarewave) 
Error in curve(expr = x, from = from, to = to, xlim = xlim, ylab = ylab, : 
    'expr' has not been evaluated to an object of length 'n' 
In addition: Warning message: 
In if ((t%%0.01) < 0.005) return(1) else if ((t%%0.01) >= 0.005) return(0) : 
    the condition has length > 1 and only the first element will be used 

Так почему это не работает?

Answer 1

Вам необходимо передать векторную функцию в plot. Поэтому либо используйте Vectorize, чтобы сделать это автоматически, либо используйте ifelse вместо if.

plot(Vectorize(squarewave)) 

или

squarewave2 <- function (t) { 

    # 0.01 == 100Hz(=2Pi) -> 1 Period of the Squarewave 
    # 0.005 == Pi   -> Half Period of the Squarewave 
    # if t smaller than a half period -> 1 
    # if t greater or equal than half a period -> 0 

    ifelse(((t %% 0.01) < 0.005),1,0) 
} 

plot(squarewave2) 

Чтобы увеличить разрешение plot, используйте параметр n см ?curve для деталей.

Answer 2

Первоначально я был неправ, хотел обновить.

Если не указано, это правильно векторизованное, вы не можете построить необработанную функцию и вместо этого построить график функции. Ниже приведен простой способ сделать это с помощью вашей конкретной функции.

sequence <- seq(from = 0,to = 0.01, by = 0.00001) 
plot(sapply(X = sequence,FUN = squarewave),type = "o")