Алгоритм программы приведен на рисунке 3

6.1 Алгоритм программы приведен на рисунке 3

		ВЫХОД

6.2 Программа (написана на языке TURBO BASIC)

LET h = .00001

Bo = 1

La = .00235

m1 = 100

D = .8

w = 40

C = .00015

ksi = 2 * 10 ^ 7

nu = 4000

E = 10

R = .00105

Ra = .05

Pi = 3.141592654#

z = 1

t1 = .00007

t2 = .00621:

t56 = t1

integral = 0

integral2 = 0

LET schet = 1

INPUT "параметры выводить на экран? n-нет"; q1$

IF q1$ = "n" OR q1$ = "N" THEN q = 0 ELSE q = 1

SCREEN 12

Uc = 0

 LOCATE 1, 45: PRINT "время t(мс)"

 LOCATE 2, 45: PRINT "белая линия I(А)"

 LOCATE 3, 45: PRINT "синяя линия Uc(В)"

 LOCATE 4, 45: PRINT "фиолетовая линия X(мм)"

 LOCATE 5, 45: PRINT "красная линия V(мм/с)"

 LOCATE 6, 45: PRINT "зеленая линия a(m/S)"

 LOCATE 7, 45: PRINT "коричневая линия Ua(В)"

0 LET i1 = I

LET Uc1 = Uc

LET x1 = X

LET V1 = V

LET xc = 0

1 LET k1i = E / La - (R / La) * i1 - Uc / La - (Bo * Pi * D * w) * V / La

LET i1 = I + h * .5 * k1i

LET k2i = E / La - (R / La) * i1 - Uc / La - (Bo * Pi * D * w) * V / La

LET i1 = I + k2i * h * .5

LET k3i = E / La - (R / La) * i1 - Uc / La - (Bo * Pi * D * w) * V / La

LET i1 = I + k3i * h

LET k4i = E / La - (R / La) * i1 - Uc / La - (Bo * Pi * D * w) * V / La

LET di = h * (2 * k2i + k1i + 2 * k3i + k4i) / 6

2 LET k1Uc = I / C

LET Uc1 = Uc + h * .5 * k1Uc

LET k2Uc = I / C

LET Uc1 = Uc + h * .5 * k2Uc

LET k3Uc = I / C

LET Uc1 = Uc + h * k3Uc

LET k4Uc = I / C

LET dUc = h * (2 * k2Uc + k1Uc + 2 * k3Uc + k4Uc) / 6

3 LET k1x = V

LET x1 = X + k1x * h * .5

LET k2x = V

LET x1 = X + k2x * h * .5

LET k3x = V

LET x1 = X + k3x * h

LET k4x = V

LET dx = (k1x + 2 * k2x + 2 * k3x + k4x) * h / 6

4 LET k1V = (Bo * Pi * D * w * I) / m1 - (nu / m1) * V1 - (ksi / m1) * X

LET V1 = V + k1V * h * .5

LET k2V = (Bo * Pi * D * w * I) / m1 - (nu / m1) * V1 - (ksi / m1) * X

LET V1 = V + k2V * h * .5

LET k3V = (Bo * Pi * D * w * I) / m1 - (nu / m1) * V1 - (ksi / m1) * X

LET V1 = V + k3V * h

LET k4V = (Bo * Pi * D * w * I) / m1 - (nu / m1) * V1 - (ksi / m1) * X

LET dv = h * (2 * k2V + k1V + 2 * k3V + k4V) / 6

LET I = I + di

LET Uc = Uc + dUc

LET X = X + dx

LET V = V + dv

LET integral2 = integral2 + ABS(I) * h

LET integral = integral + h * (ABS(I)) ^ 2

LET a = (Bo * Pi * D * w * I - nu * V - ksi * X) / m1

LET Ua = La * di / h + Bo * Pi * D * w * V + I * Ra

LET tall = tall + h

IF q = 1 THEN

LOCATE 1, 1: PRINT "t="; tall * 1000, " "

LOCATE 2, 1: PRINT "I="; I, " "

LOCATE 3, 1: PRINT "Uc="; Uc, " "

LOCATE 4, 1: PRINT "X="; X * 1000, " "

LOCATE 5, 1: PRINT "V="; V * 1000, " "

LOCATE 6, 1: PRINT "a="; a, " "

LOCATE 7, 1: PRINT "Ua="; Ua, " "

END IF

PSET (t * 20000 * .01 / t2, 250 - I * .5)

PSET (t * 20000 * .01 / t2, 250 - Uc * .1), 3

PSET (t * 20000 * .01 / t2, 250)

PSET (t * 20000 * .01 / t2, 250 - V * 100), 4

PSET (t * 20000 * .01 / t2, 250 - X * 100000), 5

PSET (t * 20000 * .01 / t2, 250 - a / 3), 2

PSET (t * 20000 * .01 / t2, 250 - Ua * .1), 6

LET t = t + h

IF t > t1 THEN

IF (((I - ikontr) / I) < .001) AND (param = 0) THEN

LOCATE 14, 45

PRINT "Тпп(мс)="; tall * 1000; : INPUT zxc

LET param = 1

END IF

LET ikontr = I

GOTO 5

END IF

 GOTO 0

5 i2 = I

V2 = V

x2 = X

7 LET k1i = E / La - (R / La) * i2 - (Bo * Pi * D * w * V / La)

LET i2 = I + h * .5 * k1i

LET k2i = E / La - (R / La) * i2 - (Bo * Pi * D * w * V / La)

LET i2 = I + k2i * h * .5

LET k3i = E / La - (R / La) * i2 - (Bo * Pi * D * w * V / La)

LET i2 = I + k3i * h

LET k4i = E / La - (R / La) * i2 - (Bo * Pi * D * w * V / La)

LET di = h * (2 * k2i + k1i + 2 * k3i + k4i) / 6

8 LET k1x = V

LET x2 = X + k1x * h * .5

LET k2x = V

LET x2 = X + k2x * h * .5

LET k3x = V

LET x2 = X + k3x * h

LET k4x = V

LET dx = (k1x + 2 * k2x + 2 * k3x + k4x) * h / 6

9 LET k1V = (Bo * Pi * D * w * I) / m1 - (nu / m1) * V2 - (ksi / m1) * X

LET V2 = V + k1V * h * .5

LET k2V = (Bo * Pi * D * w * I) / m1 - (nu / m1) * V2 - (ksi / m1) * X

LET V2 = V + k2V * h * .5

LET k3V = (Bo * Pi * D * w * I) / m1 - (nu / m1) * V2 - (ksi / m1) * X

LET V2 = V + k3V * h

LET k4V = (Bo * Pi * D * w * I) / m1 - (nu / m1) * V2 - (ksi / m1) * X

LET dv = h * (2 * k2V + k1V + 2 * k3V + k4V) / 6