. ( ) . (, .) , , , . ( .). , .
( ) 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. , .
. -, . , , , , . .
, . . .
|
|
, , . , , , , .
, , , . , . , ( ) ( ). , , , , . . - , ( ), , , .