Prediction confidence intervals using a state-space model
显示 更早的评论
Hello all,
Suppose one has estimated a state-space model using the following.
id = iddata(data',[],Ts);
sys = ssest(id, 2, 'Ts',id.Ts, 'DisturbanceModel','estimate');
Where data is a (output-only) time-series. Then, forecasting is as follows.
frc = forecast(sys, id, K);
My question is, how can one get the confidence intervals of the forecasted output, given that the system has uncertainties? The uncertainties of the model can be obtained with
[pvec, pvec_sd] = getpvec(sys) % parameter values and std deviations
An overview of the identified system is available with
present(sys)
Thank you in advance!
采纳的回答
Rajiv Singh
2013-6-7
For discrete-time systems, the forecasting model is same as the simulation (original) model. So you can use commands such as SIM and SIMSD for uncertainty analysis.
For generating response uncertainty bounds by Gauss Approximation formula:
[frc,x0e] = forecast(sys, id, K);
zu = iddata([],zeros(K,0), id.ts, 'tstart',frc.tstart);
opt = simOptions('InitialCondition',x0e);
[y,ysd] = sim(sys, zu, opt);
t = y.SamplingInstants;
plot(t, y.y, t, y.y+ysd.y, 'r', t, y.y-ysd.y,'r')
For Monte-Carlo simulation of the estimated model:
opt = simsdOptions('InitialCondition',x0e);
simsd(sys, zu, 20, opt)
更多回答（0 个）
类别
在 帮助中心 和 File Exchange 中查找有关 Uncertainty Analysis 的更多信息
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!
Translated by 
