diff options
author | Suraj Yerramilli | 2016-02-13 17:50:34 +0530 |
---|---|---|
committer | Suraj Yerramilli | 2016-02-13 17:50:34 +0530 |
commit | 1242ed4d247a9447db66278a72b5187cc034592a (patch) | |
tree | 174b24efb1c2f7e3e8f77f3cbc0071c50f0e6d82 /R/iv.R | |
parent | fced802d3097d85d020e629a1ddc3a11e2535072 (diff) | |
download | SysID-R-code-1242ed4d247a9447db66278a72b5187cc034592a.tar.gz SysID-R-code-1242ed4d247a9447db66278a72b5187cc034592a.tar.bz2 SysID-R-code-1242ed4d247a9447db66278a72b5187cc034592a.zip |
minor modification
Diffstat (limited to 'R/iv.R')
-rw-r--r-- | R/iv.R | 10 |
1 files changed, 6 insertions, 4 deletions
@@ -1,12 +1,14 @@ #' @export -iv <- function(z,order=c(0,1,0)){ +iv <- function(z,order=c(0,1,0),x=NULL){ y <- outputData(z); u <- inputData(z); N <- dim(y)[1] na <- order[1];nb <- order[2]; nk <- order[3] nb1 <- nb+nk-1 ; n <- max(na,nb1); df <- N-na-nb - # Initial Guess using ARX - mod_arx <- arx(z,order) - x <- matrix(sim(mod_arx$sys,u,sigma=0)) + if(is.null(x)){ + # Initial Guess using ARX + mod_arx <- arx(z,order) + x <- matrix(sim(mod_arx$sys,u,sigma=0)) + } padZeros <- function(x,n) c(rep(0,n),x,rep(0,n)) yout <- apply(y,2,padZeros,n=n); |