.


:




:

































 

 

 

 





. ( ) . (, .) , , , . ( .). , .

( ) 15 , ( 2×105 12×105). : W2, W4, W6, W8, W10 W12.

, .

- 1,2 /(.), . . 24 , ( ). . .

0,750,92 /(×Ѱ).

0,00001, , 50 0,5 /. , - .

, , - .

( 80) . , , . , .

 

 

, . : , , , , , , , .

, , ; . , , , , . , . : , , 150 . (, , .).

, , .

, , ( ), , . , , , . , (, , . .), , , . .

. , . , . - : 100, 150, W2, W4.

, :

 

5 510 10-12 12
W4 W6 W8 W12

 

. : F100, F150, F200, F300, F400, F500.

, . . . ( ) , , .

() :

,

; Ts ; r ; t . , .

() . , , .

: - . - .

,

D t . :

.

T (t) t + D t

 

.

, :

1. (t 0, T 0).

2. T (t) t 1 = t 0 + D t: T 1 = T 0 + + D t ( r (T 0 Ts)).

3. , . 2, T (t) t 2 = t 1 + D t . .

, , . :

1. .

2. f (t, T (t)).

3. , , (D t).

4. , .

5. .

6. .

T (0) = 25 MatLab, m -.

m - , , :

 

function F = f(t, T)

% f(t, T) ,

%

F = r(T Ts);

 

function [X, Y]=Euler(T0, t0, t1, N)

% ,

%

%

dt = (t1 t0)/N; %

%

t(1) = t0;

T(1) = T0;

%

for i = 1:N

t(i+1) = t0 + dt*i;

T(i+1) = T(i) + dt*f(t(i),T(i));

end;

%

= t;

Y = T;

f.m, Euler.m. Euler : X, Y , , .

f.m Euler.m, :

>> T0 = 25;

t0 = 0;

t1 = 10;

N = 100;

[X, Y]=Euler(T0, t0, t1, N);

plot(X, Y)

 

.

, (D t)2. N N ×(D t)2. , .

. -, . , , , , . .

, . . .

, , . , , , , .

 

 

 

, , , . , . , ( ) ( ). , , , , . . - , ( ), , , .

 





:


: 2016-11-02; !; : 716 |


:

:

,
==> ...

1865 - | 1706 -


© 2015-2024 lektsii.org - -

: 0.025 .