About Site Map Submit Contact Us Log in | Create an account
Create an account Log In
Average Rating
User Rating:
Visitors Rating:
My rating:

Write review
See full specifications

scripts default iconParticle Smoothing Expectation Maximization Procedure (Scripts) Publisher's description

Suppose you have a random process x(t)

Suppose you have a random process x(t), that is generated from time indexed densities N(m1(t),sigma1(t)) with probability alpha, and from density N(m2(t),sigma2(t)) with probability 1-alpha.

The process x(t) is directly unobservable but can be observed through y(t)=f(x(t))+noise

Given only the data y(t), we develop a methodology to estimate the parameters [m1(t),sigma1(t),m2(t),sigama2(t),alpha] that govern this specified process.

We illustrate the application of this methodology using the following process

x(t) = a*x(t-1)+Jump(t) where Jump(t)=N(m1,sigma1) with prob. alpha, and Jump(t)=0 with prob. 1-alpha
y(t) = c*exp(x(t)) + noise;

This model was developed because it can be use to generate Random Walks, AR(1) processess, Mean Reversion Processes, Stochastic Volatility, and Stochastics processes with spikes or Jumps. Our particular interest was in developing a model for estimating price spikes in electricity.

The paper attached 'An estimation Technique for Time Indexed Gaussian Mixture Model', shows how to develop the estimation procedure for this model and the same technique can be applied to any model of the form


The code attached illustrates the application of our methodology, if is extremely computionally intensive as it requires the use of particle smoothing. This is the extension of the first parameter estimation code, in which we assumed that both (y(t) and x(t)) were known. Now we only assume the data y(t) as given.

Run simulationexample.m to obtain 'figure 5' and dataexample.m to obtain 'figure 6'

System Requirements:

MATLAB 7.9 (2009b)
Program Release Status: New Release
Program Install Support: Install and Uninstall

Particle Smoothing Expectation Maximization Procedure (Scripts) Tags:

Click on a tag to find related softwares

Is Particle Smoothing Expectation Maximization Procedure (Scripts) your software?

Manage your software

Most Popular

scripts default icon ASK, OOK, FSK, BPSK, QPSK, 8PSK modulation 1.1
ASK, OOK, FSK, BPSK, QPSK, 8PSK modulation contain several functions for digital modulation simulation
scripts default icon Simulink Communication Labs 1.1
Simulink Communication Labs allows you to learn communication systems in greater depth.
scripts default icon M-QAM modulation and demodulation 1.1
M-QAM modulation and demodulation is the QAM modulation and demodulation tech.
scripts default icon LZW Compression/Decompression 1.1
LZW Compression/Decompression - Updated LZW compressor and decompressor with reasonable performance
scripts default icon InSPIRE utility to plot a 2D displacement field (Scripts) 1.0
This program plots the deformation field (displace vectors) contained in vector.txt.