SystemIdentification.jl

module SystemIdentification
if !isdefined(:DEBUG); DEBUG = false; end
export
Model,LinearModel,NonLinearModel,
Network,
Polynom,PolynomMatrix,
AR,
ARX,
RBFARX,
FitStatistics,
FitResult,
IdData,
# Functions
ar,arx,getARregressor,getARXregressor,find_na,
toeplitz, kalman

## Fit Methods =================
:LS
:LS_reg
:LM

## Types =======================
abstract Model
abstract LinearModel <: Model
abstract NonLinearModel <: Model
abstract Network <: NonLinearModel
typealias Polynom{T<:Real} Array{T,1}
typealias PolynomMatrix Array{Polynom,1}



"""
`a::Vector{Float64}`: The polynomial coeffs A(z) (not including the first 1)\n
`na::Int`\n
`bias::Bool`\n
`λ::Float64`\n
"""
type AR <: LinearModel
    a::Polynom
    na::Int
    bias::Bool
    λ::Float64
end

"""
`a::Vector{Float64}`: The polynomial coeffs A(z) (not including the first 1)\n
`b::VecOrMat{Float64}`: The polynomial coeffs B(z)\n
`na::Int`\n
`nb::Vector{Int}`\n
`bias::Bool`\n
`λ::Float64`\n
"""
type ARX <: LinearModel
    a::Polynom
    b::PolynomMatrix
    na::Int
    nb::Vector{Int}
    bias::Bool
    λ::Float64
end

type RBFARX <: Network
    na::Int
    n_centers::Int
    n_state::Int
    outputnet::Bool
    normalized::Bool
end

type FitStatistics
    RMS::Float64
    FIT::Float64
    AIC::Float64
    FitStatistics(y::Vector{Float64},yh::Vector{Float64},d) = new(rms(y-yh),fit(y,yh),aic(y-yh,d))
end

type FitResult
    error::VecOrMat{Float64}
    fit::FitStatistics
    method::Symbol
    function FitResult(y,yh,d::Int,method::Symbol)
        error = y-yh
        new(error, FitStatistics(y,yh, d),method)
    end
end

type IdData
    y::VecOrMat{Float64}
    u::VecOrMat{Float64}
    Ts::Float64
    IdData(y::VecOrMat{Float64}) = new(y,[],0)
    IdData(y::VecOrMat{Float64} ,Ts::Float64) = new(y,[],Ts)
    IdData(y::VecOrMat{Float64}, u::VecOrMat{Float64}) = new(y,u,0)
end

## Helper functions
rms(x) = sqrt(mean(x.^2))
sse(x) = sum(x.^2)
fit(y,yh) = 100 * (1-rms(y-yh)./rms(y-mean(y)));
aic(x,d) = log(sse(x)) + 2d/length(x)

include("armax.jl")
include("kalman.jl")
include("PCA.jl")
include("kernelPCA.jl")
include("toeplitz.jl")

end