1.
-
21
: -05-01 . .
: . .
2009
:
;
;
.
:
21
, 0, 2,9
5 ()
11 (-)
10 (-)
7 ()
12 (-)
Y(1) 0,80
Y(2) 0,05
Y(3) 0,03
Y(4) 0,03
Y(5) 0,09
4,03129
1214,645
221,205
- 4,12768
1461,925
215,073
- 4,11810
1474,679
213,686
4,07427
1345,087
219,516
- 4,11103
1454,328
215,411
:
- ;
- .
.
∑(Y(i)/K(i))=1, (1)
Y(i) i- , , ;
K(i) i- ;
N .
K(i)=P(i)/P0, (2)
P(i) , ;
P0 , .
P(i) :
lgP(i)=A(i)-B(i)/(C(i)+t), (3)
(i), B(i),C(i) ;
t , .
(3)
t=B(i)/(A(i)-lgP0)-C(i), (4)
i- 0.
(1) , , , :
Z=(∑(Y(i)/K(i)))-1=0. (5)
:
1. .
1. .
:
|
|
1. 0 (4). [t1;t2]. Z1=f(t1) Z2=f(t2).
2.
t=(t1+t2)/2 (6)
3. P(i) K(i) (3) (2).
4. Z (5), |t-t2|<ε, t .
5. |t-t2|>ε, Z*Z1<0, t2=t, Z2=Z, t1=t, Z1=Z. , 2, .
2. .
2. .
:
1. 0 (4). [t1;t2].
2. Z1=f(t1) Z2=f(t2).
3.
t=t1-(t2-t1)*Z1/(Z2-Z1). (7)
4. P(i) K(i) (3) (2).
5. Z (5), |t-t2|<ε, t .
6. |t-t2|>ε, Z*Z2>0, t1=t, Z1=Z, t2=t, Z2=Z. , 2, .
(7) () Z=0.
3. -
program lab1;
uses crt;
label 1,2,3,4,5;
var a,b,c,ti,y,k1,k2,p1,p2,p,ki:array [1..4] of real;
p0,t1,m,t,t2,z1,z2,z,s,s1,s2,e:real;
i,k,n:integer;
begin
clrscr;
writeln(' ');readln(n);
writeln(' , 0, ');readln(p0);
for i:=1 to n do begin
writeln(' ',i);
writeln(' ,,');readln(a[i],b[i],c[i]);
writeln(' Y');readln(y[i]);
end;
writeln(' ');
for i:=1 to n do begin
ti[i]:=b[i]/(a[i]-(ln(p0))/(ln(10)))-c[i];
writeln('T',i,'=',ti[i]:5:2);
end;
writeln(' ');readln(t1,t2);
writeln(' ');readln(e);
for i:=1 to n do begin
p1[i]:=exp((a[i]-b[i]/(c[i]+t1))*ln(10));
p2[i]:=exp((a[i]-b[i]/(c[i]+t2))*ln(10));
k1[i]:=p1[i]/p0;
k2[i]:=p2[i]/p0;
end;
s1:=0; s2:=0;
for i:=1 to n do begin
s1:=s1+y[i]/k1[i];
s2:=s2+y[i]/k2[i];
end;
z1:=s1-1; z2:=s2-1;
if z1*z2<0 then goto 1 else goto 2;
1: k:=0; m:=t2;
3: k:=k+1;
t:=(t1+t2)/2;
for i:=1 to n do begin
p[i]:=exp((a[i]-b[i]/(c[i]+t))*ln(10));
p2[i]:=exp((a[i]-b[i]/(c[i]+t2))*ln(10));
ki[i]:=p[i]/p0;
k2[i]:=p2[i]/p0;
end;
s:=0; s2:=0;
for i:=1 to n do begin
s:=s+y[i]/ki[i];
s2:=s2+y[i]/k2[i];
end;
z:=s-1; z2:=s2-1;
if abs(t-m)<e then goto 4; m:=t;
if z*z2<0 then t1:=t else t2:=t; goto 3;
2: writeln(' '); goto 5;
4: writeln(' T=',t:5:2);
writeln(' ',k:5);
|
|
5: readln;
end.
.
T1=119,14 ;
2=183,78 ;
3=189,71 ;
4=152,89 ;
T5=183,18 .
e=0,1 .
t=141,137 .
k=220.
1. 0:
0, | ||||||
, | ||||||
2,9 | 119,14 | 183,78 | 189,71 | 152,89 | 183,18 | 141,137 |
3,0 | 120,55 | 185,39 | 191,34 | 154,42 | 184,80 | 142,647 |
3,1 | 121,92 | 186,96 | 192,93 | 155,90 | 186,37 | 144,022 |
3,2 | 123,26 | 188,49 | 194,48 | 157,35 | 187,91 | 145,364 |
3,3 | 124,57 | 189,99 | 195,99 | 158,77 | 189,41 | 146,674 |
3,4 | 125,86 | 191,45 | 197,48 | 160,15 | 190,88 | 147,955 |
3,5 | 127,11 | 192,87 | 198,92 | 161,51 | 192,31 | 149,208 |
3,6 | 128,33 | 194,27 | 200,34 | 162,83 | 193,71 | 150,434 |
4. .
, 0
2. .
ti pi ki yi 1, 2 k1, k2 ,, z t1 t2 t e z1, z2 p0 i n k | i- , i- , i- i- , , i- t1 t2, i- t1 t2 i- ( (5)) , , , z t1 t2 ( (5)) , |
3. .
, | k | , | |
0 200 | 142,05 | ||
100 200 | 141,38 | ||
140-150 | 141,21 | ||
140-145 | 141,20 |
5. - ()
program lab2;
uses crt;
label 1,2,3,4,5;
var a,b,c,ti,y,k1,k2,p1,p2,p,ki:array [1..4] of real;
p0,t1,m,t,t2,z1,z2,z,s,s1,s2,e:real;
i,k,n:integer;
begin
clrscr;
writeln(' ');readln(n);
writeln(' , 0, ');readln(p0);
for i:=1 to n do begin
writeln(' ',i);
writeln(' ,,');readln(a[i],b[i],c[i]);
writeln(' Y');readln(y[i]);
end;
writeln(' ');
for i:=1 to n do begin
ti[i]:=b[i]/(a[i]-(ln(p0))/(ln(10)))-c[i];
writeln('T',i,'=',ti[i]:5:2);
end;
writeln(' ');readln(t1,t2);
writeln(' ');readln(e);
for i:=1 to n do begin
p1[i]:=exp((a[i]-b[i]/(c[i]+t1))*ln(10));
p2[i]:=exp((a[i]-b[i]/(c[i]+t2))*ln(10));
k1[i]:=p1[i]/p0;
k2[i]:=p2[i]/p0;
end;
s1:=0; s2:=0;
for i:=1 to n do begin
s1:=s1+y[i]/k1[i];
s2:=s2+y[i]/k2[i];
end;
z1:=s1-1; z2:=s2-1;
if z1*z2<0 then goto 1 else goto 2;
|
|
1: k:=0; m:=t2;
3: k:=k+1;
t:=(z1*t2-z2*t1)/(z1-z2);
for i:=1 to n do begin
p[i]:=exp((a[i]-b[i]/(c[i]+t))*ln(10));
p2[i]:=exp((a[i]-b[i]/(c[i]+t2))*ln(10));
ki[i]:=p[i]/p0;
k2[i]:=p2[i]/p0;
end;
s:=0; s2:=0;
for i:=1 to n do begin
s:=s+y[i]/ki[i];
s2:=s2+y[i]/k2[i];
end;
z:=s-1; z2:=s2-1;
if abs(t-m)<e then goto 4; m:=t;
if z*z2<0 then t1:=t else t2:=t; goto 3;
2: writeln(' '); goto 5;
4: writeln(' T=',t:5:2);
writeln(' ',k:5);
5: readln;
end.