ARMA model for temperature simulation

3 次查看(过去 30 天)
Peter
Peter 2012-6-20
Hello, I really checked all relevant posts and I could not find the right answer. I'm also not an expert in ARMA models so please apologize my (maybe) naive question. I want to model daily temperature using an ARMA or AR model. So I think I know how to fit it now (in the system identification toolbox). So I fitted an ARMA(1,1) model to 31 January temperature values:
model=armax(JanTemp(:,1),[1 1]);
Discrete-time IDPOLY model: A(q)y(t) = C(q)e(t)
A(q) = 1 - 0.9658 q^-1
C(q) = 1 - 0.653 q^-1
But how can I now simulate from this fitted model stochastically? What I want to do is simulating say 31 artifical January days that follow the statistical behaviour of my historical dataset but represent a new realization? I know there is a simulation routine "sim" but I got stuck when trying to understand the documentation.
I unfortunately don't have the econometrics toolbox which seems to have some ARMA algorithms.
Your help would be highly appreciated! - Peter

回答(2 个)

Rajiv Singh
Rajiv Singh 2012-6-27
See SIM and SIMSD functions in System Identification Toolbox.
When using SIM, add the 'noise' qualifier as input argument to perform stochastic simulation, as in:
y = sim(model, zeros(100,0), 'noise')
Post examples of things you tried, then I may be able to give you more targeted advice.
rajiv

grapevine
grapevine 2012-6-27
If you have enough measurements of the inputs and the outputs of the process you might perform a system identification There are severals techniques :
Least means squares identification it's the easy one, I return a linear estimation of the transfer function
Luenberger observer it returns an estimation of the state variables
Kalman filter it's a powerful algorithm which uses a series of measurements observed over time, containing noise, and produces and estimation of the states
<http://personnel.supaero.fr/alazard-daniel/Pdf/cours_Kalman.pdf> (if you speak French)
you could be interested at this Matlab toolbox
and a this function: * tfest*
wish u the best

类别

Help CenterFile Exchange 中查找有关 Linear Model Identification 的更多信息

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by