summaryrefslogtreecommitdiff
path: root/R/iv.R
diff options
context:
space:
mode:
authorSuraj Yerramilli2016-02-13 17:50:34 +0530
committerSuraj Yerramilli2016-02-13 17:50:34 +0530
commit1242ed4d247a9447db66278a72b5187cc034592a (patch)
tree174b24efb1c2f7e3e8f77f3cbc0071c50f0e6d82 /R/iv.R
parentfced802d3097d85d020e629a1ddc3a11e2535072 (diff)
downloadSysID-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.R10
1 files changed, 6 insertions, 4 deletions
diff --git a/R/iv.R b/R/iv.R
index b0f9aff..9257766 100644
--- a/R/iv.R
+++ b/R/iv.R
@@ -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);