R/obsv.R

In benubah/control: A Control Systems Toolbox

#' @title Observability Matrix
#'
#' @description
#' This function creates the observability matrix.
#'
#' @param A State-space matrix, A
#' @param C State-space matrix, C
#'
#' @return  \code{obsv(A, C)} returns the observability matrix, \code{obsvm}.
#' where
#'  obsvm =  | C CA CA^2 ... CA^(n-1) |
#'
#'@seealso \code{\link{ctrb}}
#'
#' @examples
#' A <- rbind(c(0,1), c(-25,-4))
#' C <- rbind(c(1,0), c(0,1))
#'  obsv(A, C)
#'
#' @export

obsv <- function(A, C) {
  obsm <- t(ctrb (t(A), t(C)))
  return(obsm)
}