root/applications/dual/SIDP/bakalarka/SIDPS/matlab/motor/generate_realization.m @ 1351

Revision 1351, 0.6 kB (checked in by zimamiro, 13 years ago)
Line 
1function realization=generate_realization(apriori,H0,H_tilde,C,regulator_parameters,system)
2
3H=add_Hc_tilde0(H0,apriori,regulator_parameters); % nehlidam pozitivni definitnost !
4realization=0;
5system.S=zeros(1,3);
6for o=1:regulator_parameters.horizont_simulace
7    [u system]=simple_control(system,H(:,1));
8    u=u+interpol(H,H_tilde,C);
9    system.x=H(:,1)+sqrt(diag(H(:,2:5))).*randn(4,1);
10    system.x=get_next(system,system.x,u);
11    H=kalman_filter(H,u,system);
12    realization=realization+u'*system.gamma*u+(system.x-system.x_opt)'*system.ksi*(system.x-system.x_opt);
13end
14end
Note: See TracBrowser for help on using the browser.