%---etalonmudel
% kõik on ühikringi sees
% hea sumbuvusega vt zgrid
td=1
z1=0.5, z2=0.48, am1=-z1-z2, am2=z1*z2
bm1=1;
bm2=-0.1;
Beta=(1+am1+am2)/(bm1+bm2);
bm1=bm1*Beta
bm2=bm2*Beta
Zm=roots([1 am1 am2])
if all(isreal(Zm)), plot(Zm,0,'*'); else plot(Zm,'*');end; zgrid; line([0 -1;0 1],[-1 0;1 0]); 

%---parameetrid
% term obj T1=1, T2=7, K=0.015 deg/W, td=1 [min]
a1=-1.5
a2=0.7
b1=1
b2=0.5
Z=roots([1 a1 a2])
a0=1
r00=1
R=(r00/b1)*[b1 b2]
s1=(am1-a1)/b1
s2=(am2-a2)/b1
%--- lambda on seotud iteratiivse parameetrite hindamisega
lamda=0.5;
aw_dsim2
