судна
Метод трапецій:
program
s: real;
h: real;
m: integer;: = (ba)/n;: = (f (a) + f (b))/2; m: = 1 to n-1 do s: = s + f (a + m * h);: = s * h;; ('vvedi A, B'); (a, b);: = 100;: = trap (a, b, n); ('r = ', r: 6:3);;.
Результат програми:
vvedi A, B
.994
r = 49.111
Час розгону до Vст t = 49,111 сек.
5.4 Знаходження шляху розгону судна
Метод Ейлера 1 модифікація:
Program x0, xn, y0, z0, V, D '); (x0, xn, y0, z0, V, D); (' Ti vvel x0 = ', x0: 6:2,' xn = ', xn: 6 : 2, 'D =', D: 6:2, 'y0 =', y0: 6:2, 'z0 =', z0: 6:2, 'V =', V: 6:2);: = x0; y: = y0; z: = z0;: = V/D;: = xn/D; x <= xn do; ('t =', x: 8:4, 'S =', y: 8 : 4, 'V =', z: 8:4);: = x + h/2;: = y + h/2 * z;: = z + N/2;: = x + h;: = y + h * z1;: = z + N;
end;;.
Результат програми:
В
Шлях розгону до Vст S = 441,8517 м
Графік S (t)
В
Графік V (t)
В
5.5 Знаходження енергії розгону судна
Обчислюємо енергію розгону за формулою: E = m/2? v 2 (t) dt
V (t) = 0,36639 t
Метод трапецій:
program to n-1 do s: = s + f (a + m * h);: = s * h;; ('vvedi A, B'); (a, b);: = 100;: = trap (a , b, n); ('r =', r: 6:3);;.
Результат програми:
vvedi A, B
.111
r = 107071869.240
Енергія розгону судна до Vст E = 107,1 МДж
6. Модельна задача 3
6.1 Апроксимація функцій R (V) і T (v) поліномом 3 ступеня
В
6.2 Знаходження стаціонарної швидкості
F (x) = -52,567 x 3 +1687,1 x 2 -17668,6 x +82315
Відділення кореня кроковим методом:
program shag;
var a, b, h '); (a, b, h);: = 0;: = a + h; x <= b do: = xh;: = x;: = y (x1);: = y ( x2); y1 * y2 <0 then: ...