diff options
| author | Suraj Yerramilli | 2016-03-20 13:39:37 +0530 |
|---|---|---|
| committer | Suraj Yerramilli | 2016-03-20 13:39:37 +0530 |
| commit | 24c7e0e8730576f1ddb46af63acdf9cd913aa995 (patch) | |
| tree | 9ee405b1d5140df7ef39c4b77be74b5218aa7bea | |
| parent | 7e376c7936d7d9ea9d4bbdeefbbceadf060d8967 (diff) | |
| download | SysID-R-code-24c7e0e8730576f1ddb46af63acdf9cd913aa995.tar.gz SysID-R-code-24c7e0e8730576f1ddb46af63acdf9cd913aa995.tar.bz2 SysID-R-code-24c7e0e8730576f1ddb46af63acdf9cd913aa995.zip | |
minor corrections
| -rw-r--r-- | R/nonparam.R | 4 |
1 files changed, 2 insertions, 2 deletions
diff --git a/R/nonparam.R b/R/nonparam.R index ad8254d..3096f3e 100644 --- a/R/nonparam.R +++ b/R/nonparam.R @@ -193,7 +193,7 @@ spa <- function(x,winsize=NULL,freq=NULL){ } G <- array(0,c(nout,nin,length(freq))) - spec <- array(0,c(nout,nout,length(freq))) + spec <- array(0,c(nout,1,length(freq))) for(i in 1:nout){ phi_y <- sapply(freq,cov2spec,Ryy[i,i,],M) temp <- phi_y @@ -203,7 +203,7 @@ spa <- function(x,winsize=NULL,freq=NULL){ G[i,j,] <- phi_yu/phi_u temp <- temp - phi_yu*Conj(phi_yu)/phi_u } - spec[i,i,] <- temp + spec[i,1,] <- temp } out <- idfrd(G,matrix(freq),deltat(x),spec) return(out) |