Inlägg
av SvenW »
Här är en något dåligt granskad octaveberäkning:
T_R = [ -30 , 27.9;
-25 , 21.4;
-20 , 16.5 ;
-15 , 12.8 ;
-10 , 10.0 ;
-5 , 7.9 ;
0 , 6.3 ;
5 , 5.0 ;
10 , 4.1 ;
15 , 3.3 ;
20 , 2.7 ]
R= T_R(:,2)*1000
T= T_R(:,1)+273.15
% beta= ( log(R(T1/T2)) ) / ((1/T1) - (1/T2));
beta=(1:8);
beta(1) = ( log(R(1)/R(11 ) )) / (1./T(1) - 1./T(11) );
beta(2) = ( log(R(2)/R(11 ) )) / (1./T(2) - 1./T(11) );
beta(3) = ( log(R(3)/R(11 ) )) / (1./T(3) - 1./T(11) );
beta(4) = ( log(R(4)/R(11 ) )) / (1./T(4) - 1./T(11) );
beta(5) = ( log(R(5)/R(11 ) )) / (1./T(5) - 1./T(11) );
beta(6) = ( log(R(6)/R(11 ) )) / (1./T(6) - 1./T(11) );
beta(7) = ( log(R(7)/R(11 ) )) / (1./T(7) - 1./T(11) );
beta(8) = ( log(R(8)/R(11 ) )) / (1./T(8) - 1./T(11) );
beta
% 3329.3 3346.5 3358.2 3364.8 3366.8 3375.8 3392.3 3349.6