项目作者: compops

项目描述 :
Hierarchical Bayesian approaches for robust inference in ARX models
高级语言: Matlab
项目地址: git://github.com/compops/rjmcmc-sysid2012.git
创建时间: 2013-03-19T07:21:42Z
项目社区:https://github.com/compops/rjmcmc-sysid2012

开源协议:GNU General Public License v3.0

下载


rjmcmc-sysid2012

Reversible-Jump Markov Chain Monte Carlo sampler for inference in ARX models with Student-t innovations

This code was downloaded from < https://github.com/compops/rjmcmc-sysid2012 > and contains the code used to produce the results in

J. Dahlin, F. Lindsten, T. B. Schön and A. Wills, Hierarchical Bayesian ARX models for robust inference. In the Proceedings of the 16th IFAC Symposium on System Identification, Brussels, Belgium, July 2012.

which is available from < http://www.johandahlin.com/files/publications/Dahlin2012-RJMCMC.pdf >.

Included files

RUNME_simulationstudy
The code reproduces the comparision presented in Section 5.1. Note that the
results from that section is based on 25 000 randomly generated systems and
the provided code only simulates 100 systems.

RUNME_EEG
The code reproduces a similar output as presented in Fig 3 of Section 5.3.
Using real EEG data, kindly provided by Eline Borch Petersen and Thomas Lunner
at Eriksholm Research Centre, Oticon A/S, Denmark.

Helpfiles and subroutines are found in the folder “helpers”:

  • buildPhi/buildY: builds the Phi and Y matrix containing the values of u and y.
  • dexprnd: generates a discrete exponentially distributed random number
  • genrndARXmodel: generates random ARX models with Gaussian/Student’s t-noise, outliers and missing data.
  • igamrnd: generates inverse gamma distributed random variables
  • method_createdatasets: divides the output from the random ARX generation function into estimation and validation data.
  • method_validate: validates the model and calculates the performance measures (RMSE and VAF)