Реферат Дослідження перехідних процесів в механічній системі

ustify"> .197 -0.003 0.002 0.403 0.145

.204 -0.000 0.003 0.399 0.141

.210 0.003 0.004 0.393 0.136

.217 0.005 0.005 0.384 0.131

.224 0.008 0.006 0.373 0.124

.231 0.010 0.007 0.360 0.117

.238 0.013 0.007 0.345 0.110

.244 0.015 0.008 0.328 0.102

.251 0.017 0.009 0.309 0.093

.258 0.020 0.009 0.289 0.084

.265 0.021 0.010 0.267 0.074

.271 0.023 0.010 0.244 0.064

.278 0.025 0.011 0.220 0.054

.285 0.026 0.011 0.194 0.044

.292 0.028 0.012 0.168 0.033

.298 0.029 0.012 0.141 0.023

.305 0.030 0.012 0.114 0.012

.312 0.031 0.012 0.086 0.002

.319 0.031 0.012 0.058 -0.009

.325 0.032 0.012 0.030 -0.019

.332 0.032 0.012 0.002 -0.029

.339 0.032 0.012 -0.026 -0.039

.346 0.032 0.011 -0.053 -0.049

maxy1 maxy2 maxV1 maxV2

.062 0.025 1.076 0.446

6.4 Програма Модифікований метод Ейлера program ModEiler

typeodnmas = array [1 .. 201] of real;

var then T0: = 3 * T1T0: = 3 * T2;: = T0/n;: = 1; [i]: = 0; [i]: = y10; y2 [i]: = y20; [i]: = V10; V2 [i]: = V20; (f2, 't y1 y2 V1 V2'); [i +1]: = t [i] + ht; [i +1]: = y1 [i] + ht * (V1 [i] + V1 [i] + ht * fv1 (y1 [i], y2 [i], V1 [i], V2 [i], t [i]))/2; [i +1]: = V1 [i] + ht * (fv1 (y1 [i], y2 [i], V1 [i], V2 [i], t [i]) + fv1 (y1 [i] +

{di: = i/10; (frac (di) = 0) or (frac (di) = 0.1) then} (f2, t [i]: 3:3, y1 [i ]: 8:3, y2 [i]: 8:3, V1 [i]: 8:3, V2 [i]: 8:3);: = i +1; i> n +1;

: = y1 [1];: = y2 [1];: = V1 [1];: = V2 [1];

i: = 2 to n doAbs (y1 [i])> maxy1 then: = Abs (y1 [i]); Abs (y2 [i])> maxy2 then: = Abs ( y2 [i]); Abs (V1 [i])> maxV1 then: = Abs (V1 [i]); Abs (V2 [i])> maxV2 then: = Abs (V2 [i]);; ( f2); (f2, 'maxy1', 'maxy2', 'maxV1', 'maxV2'); (f2, maxy1: 8:3, maxy2: 8:3, maxV1: 8:3, maxV2: 8:3) ; (f1); (f2);.

6.5 Файл виведення результату

y1 y2 V1 V2

.000 0.000 0.000 0.000 0.000

.007 0.000 0.000 0.005 0.000

.014 0.000 0.000 0.020 0.004

.020 0.000 0.000 0.043 0.014

.027 0.001 0.000 0.075 0.030

.034 0.001 0.000 0.113 0.052

.041 0.002 0.001 ...