Procedura regr_tau.m

function [a1,b1] = regr_tau(i1,i0,ad1,ad0)

zz=zeros(2,1);

aa=[i1 1;i0 1];

zz=inv(aa)*[24./ad1;24./ad0]

a1=zz(1,1);

b1=zz(2,1);

end

Procedura heln.m

function [hel]=heln(II,Tq)

xx =[-78.2625093995407;0.368102370745303;-0.0129841559173611;-0.000292630771293419;

-7.68518789123352e-007; -1.83708621082882e-007; 1.6873894087207e-009];

ad1=24./(0.01896555529634*II-0.09411383207540)

hel= xx(7,1)*Tq*Tq*ad1*ad1+xx(6,1)*Tq*Tq*ad1+xx(5,1)*Tq*ad1*ad1+xx(4,1)*Tq*ad1;

hel=hel +xx(3,1)*Tq+xx(2,1)*ad1 +xx(1,1);

heln=hel;

regr_tau(7.845,8.65,438.99,343.16)

zz =

1.8966e-002

-9.4114e-002

heln(7.845,294)

ad1 =

4.3899e+002

ans =

1.9340e+001

---