R/stepplot.R

In benubah/control: A Control Systems Toolbox

# Step Response Plot
#' @export
stepplot <- function (sys, t = NULL, input = 1) {

  if(is.null(t)) {
    tf <- seq(0, 5, 0.01)
  }
  if(!is.null(t)) {
    tf <- t

  }

  res <- step(sys, tf, input[1])
  y <- res$y
  t <- res$t
  if (length(input) == 1) {
  # Plot the results for all rows of y against t
     if (nrow(y) == 1) {
       graphics::plot(t, y, type = "l", col = "blue", xlab = "Time, sec", ylab = paste("y"), main = "Step response")
       graphics::grid(5, 5)
      } else {
        graphics::par(mfrow = c(1, nrow(y)))
        for (i in 1:nrow(y)) {
          graphics::plot(t, y[i, ], type = "l", col = "blue", xlab = "Time, sec", ylab = paste("y", i), main = "Step response")
          graphics::grid(5, 5)
        }
        graphics::par(mfrow = c(1,1))
      }
  }
  if (length(input) > 1) {
    for(i in 1:length(input)){
      stepplot(sys, tf, input = i)
      }
  }
}