.


:




:

































 

 

 

 





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.





:


: 2016-10-22; !; : 837 |


:

:

, .
==> ...

815 - | 653 -


© 2015-2024 lektsii.org - -

: 0.043 .