root/applications/dual/SIDP/transformace(3)/alstr.asv @ 1107

function  [ztrata u y b P]=alstr(regulator_parameters,system,apriori)
horizont=regulator_parameters.horizont;
sigma=system.sigma;

y=zeros(horizont);
b=zeros(horizont);
P=zeros(horizont);
u=zeros(horizont);
ztrata=0;

y(1)=apriori(1);
b(1)=apriori(2);
P(1)=apriori(3);


    for i=1:size(H,1)-1   
        eta=abs(y(i))/sigma;
        beta=abs(b(i))/sqrt(P(i));
        ny=(0.56+beta)/(2.2+0.08*beta+beta^2)*eta+1.9/(1.7+beta^4);
        
        u(i)=-ny*sigma/sqrt(P(i))*sign(y(i)*b(i));
        
        y(i+1)=y(i)+system.b*u(i)+sigma*randn;
        K=u(i)*P(i)/(P(i)*u(i)^2+sigma^2);
        b(i+1)=b(i)+K*(y(i+1)-y(i)-b(i)*u(i));
        P(i+1)=(1-K*u(i))*P(i);   
        
        ztrata=ztrata+(y(i+1)-system.yr(i))^2;
        
    end
 end