21
: .. 03-01
..
: . ..
2007
:
;
;
;
.
:
21
K t . K=f(t) K=K0 e -E/RT, K0 ; E ; T ; R , .
, | |||||||||
, -1 | 0,05 | 0,10 | 0,20 | 0,40 | 0,68 | 1,00 | 1,54 | 2,35 | 4,63 |
( , ) :
1. , , : .
2. ( ). ln K = f(t).
3. ln K, K.
4. .
, | |||||||||
ln() | -2,9957 | -2,3026 | -1,6094 | -0,9163 | -0,3857 | 0,0000 | 0,4318 | 0,8544 | 1,5326 |
ln K = f(t).
, , , y=a+bx. Turbo Pascal, .
program lab4;
uses crt;
const n=2;nn=9;tabl_stroki=13; tabl_stolb=12;
fisher: array [0..tabl_stroki,0..tabl_stolb] of real =
((0.0,1.0,2.0,3.0,4.0,5.0,6.0,7.0,8.0,9.0,10,12,16),
(1.0,161,200,216,225,230,234,237,239,241,242,244,246),
(2.0,18.5,19.0,19.1,19.2,19.3,19.3,19.3,19.3,19.3,19.3,19.4,19.4),
(3.0,10.1,9.55,9.28,9.12,9.01,8.94,8.88,8.84,8.81,8.78,8.74,8.69),
(4.0,7.71,6.94,6.59,6.39,6.26,6.16,6.09,6.04,6.00,5.96,5.91,5.84),
(5.0,6.61,5.79,5.41,5.19,5.05,4.95,4.88,4.82,4.78,4.74,4.68,4.60),
(6.0,5.99,5.14,4.76,4.53,4.39,4.28,4.21,4.15,4.10,4.06,4.00,3.92),
(7.0,5.59,4.74,4.35,4.12,3.97,3.87,3.79,3.73,3.68,3.63,3.57,3.49),
(8.0,5.32,4.46,4.07,3.84,3.69,3.58,3.50,3.44,3.39,3.34,3.28,3.20),
(9.0,5.12,4.26,3.86,3.63,3.48,3.37,3.29,3.23,3.18,3.13,3.07,2.98),
(10.0,4.96,4.10,3.71,3.48,3.33,3.22,3.14,3.07,3.02,2.97,2.91,2.82),
|
|
(12.0,4.75,3.88,3.49,3.26,3.11,3.00,2.92,2.85,2.80,2.76,2.69,2.60),
(14.0,4.60,3.74,3.34,3.11,2.96,2.85,2.77,2.70,2.65,2.60,2.53,2.44),
(16.0,4.49,3.63,3.24,3.01,2.85,2.74,2.66,2.59,2.54,2.49,2.42,2.33));
t:array [1..nn] of real =(100,150,200,250,300,350,400,450,500);
var i:integer;
sx,sx2,sy,sxy:real;
y,yy,w:array [1..nn] of real;
eps,sigma_y,sigma_ost,ysr,F:real;
stroka,stolb:integer;
a:array [1..n,1..n] of real;
b,xx:array [1..n] of real;
s,z,zz,ko,e:real;
j,k:integer;
begin
clrscr;
sx:=0;
sx2:=0;
sy:=0;
sxy:=0;
ko:=0;
e:=0;
z:=0;
zz:=0;
{ }
for i:=1 to nn do begin
write ('t= ',t[i]:6:1,'() k = '); readln (w[i]);
y[i]:=ln(w[i]);
end;
{ }
for i:=1 to nn do begin
sx:=sx+t[i];
sx2:=sx2+t[i]*t[i];
sy:=sy+y[i];
sxy:=sxy+t[i]*y[i];
end;
writeln ('sx= ',sx:5:5);
writeln ('sx2= ',sx2:5:5);
writeln ('sy= ',sy:5:5);
writeln ('sxy= ',sxy:5:5);
readln;
b[1]:=sy;
b[2]:=sxy;
a[1,1]:=nn; a[1,2]:=sx;
a[2,1]:=sx; a[2,2]:=sx2;
for k:=1 to n-1 do begin
for i:=(k+1) to n do begin
s:=-a[i,k]/a[k,k];
for j:=k to n do
a[i,j]:=a[i,j]+s*a[k,j];
b[i]:=b[i]+s*b[k];
end; end;
xx[n]:=b[n]/a[n,n];
for i:=(n-1) downto 1 do begin
s:=0;
for j:=(i+1) to n do
s:=s+a[i,j]*xx[j];
xx[i]:=(b[i]-s)/a[i,i];
end;
{ }
writeln (' :');
writeln ('k=a+b*t, ');
writeln ('a= ', xx[1]:5:9);
writeln ('b= ', xx[2]:5:9);
readln;
writeln ('T () ','K () ','K () ',' ');
for i:=1 to nn do begin
yy[i]:=xx[1]+xx[2]*t[i];
eps:=(yy[i]-y[i])*(yy[i]-y[i]);
z:=yy[1];
zz:=yy[9];
writeln (t[i]:7:2,w[i]:10:4,exp(yy[i]):15:4,eps:10:2);
end;
{ }
E:=8.31*373.15*773.15*(zz-z)/(773.15-373.15);
Ko:=exp(zz)/(exp(-E/(8.31*773.15)));
writeln (' Ko=',ko);
writeln (' E=',e:8:5);
readln;
{ }
ysr:=sy/nn;
for i:=1 to nn do
begin
sigma_y:=sigma_y+sqr(y[i]-ysr);
sigma_ost:=sigma_ost+sqr(y[i]-yy[i]);
end;
sigma_y:=sigma_y/(nn-1);
sigma_ost:=sigma_ost/(nn-n);
writeln (' ',sigma_y:10:2);
writeln (' ',sigma_ost:10:2);
F:=sigma_y/sigma_ost;
writeln (' F=',f:10:2);
stroka:=1;stolb:=1;
for i:=1 to tabl_stroki do
if (nn-n)=round (fisher[i,0]) then stroka:=i;
for i:=1 to tabl_stolb do
if (nn-1)=round (fisher[ 0,i]) then stolb:=i;
writeln (' F=',fisher[stroka,stolb]:5:2);
if fisher[stroka,stolb]<f then writeln (' ')
else writeln (' ');
readln;
end.
w | |
y | |
yy | |
sx | |
sx2 | |
sy | |
sxy | ∙ |
eps | |
sigma_y | |
sigma_ost | |
ysr | |
F | |
stroka | |
stolb | |
xx[1] | |
xx[2] | b |
ko | |
e | |
a,b,s,j,k | |
z | 100 |
zz | 500 |
|
|
-
.
ln K = f(t):
ln K = -3.8573 + 0.0109 t.
, | |||||||||
ln() | -2,9957 | -2,3026 | -1,6094 | -0,9163 | -0,3857 | 0,0000 | 0,4318 | 0,8544 | 1,5326 |
ln() (.) | -2,7712 | -2,2281 | -1,6851 | -1,1420 | -0,5990 | -0,0559 | 0,4871 | 1,0301 | 1,5732 |
0,05 | 0,01 | 0,01 | 0,05 | 0,05 | 0,00 | 0,00 | 0,03 | 0,00 |
ln K = f(t).
ln K = f(t). ; .
ln K K. K = f(t).
, | |||||||||
, -1 (.) | 0,05 | 0,10 | 0,20 | 0,40 | 0,68 | 1,00 | 1,54 | 2,35 | 4,63 |
, -1 (.) | 0,0626 | 0,1077 | 0,1854 | 0,3192 | 0,5494 | 0,9456 | 1,6276 | 2,8015 | 4,8220 |
0,05 | 0,01 | 0,01 | 0,05 | 0,05 | 0,00 | 0,00 | 0,03 | 0,00 |
K = f(t). ; .
:
,
.
, . ( ).
: .
: .
, , .
:
2.24
0.03
79.51
3.73
, , .
K = f(t) , , , , K = f(t) .
.
. , :
.
ln K = f(t) , : y=a+bx. . , , , .
K = f(t). , .
|
|
:
,
.
:
.
.