Диференційні рівняння як основа математичного опису енергетичної системиЕкспертна система контролю

“Диференційні рівняння як основа математичного опису енергетичної системи. Експертна система контролю працездатності енергетичної системи.”

Зміст:

1.1Вихідні данні для реалізації системи звичайних диференційних рівнянь.

1.2. Математична основа засобу Рунге – Кутта.

1.3. Блок схема програми.

1.4. Реалізація програми за засобом Рунге – Кутта.

1.5. Результати реалізації системи диференційних рівнянь за засобом Рунге – Кутта.

1.6. Графічне представлення результатів за засобом Рунге – Кутта.

1.7.Математична основа способу Мілна

1.8.Блок-схема реалізації головного модуля, програми та графічної частини зазасобом Мілна

1.9. Реалізація програми за способом Мілна

1.10.Результати реалізації програми за способом Мілна

1.11. Графічне представлення результатів за способом Мілна

1.12.Висновок по першій частині

2. Реалізація контролю працездатності енергетичної системи

2.1. Архітектура експертних систем

2.2. Математичне моделювання енергетичної систем

2.3. Логічне моделювання контролю працездатності

2.4. Початкові дані для реалізації контролю працездатності енергетичної системи

2.5. Алгоритм реалізації контролю працездатності

2.6. Реалізація контролю працездатності енергетичної системи

2.7.Графічне представлення результатів реалізації

Висновок.

1. 1. Вихідні дані для реалізації системи звичайних диференційних рівнянь

Системадиференційних рівнянь :

Початкові умови : А=0, В=1

t(0)=0, x(0)=0, y(0)=0

Задана точність: Е=

Обираємо с=6

1.2 Математична основа засобу Рунге–Кутта

Засіб Рунге -Кутта можливо получити, якщо разкласти у ряд Тейлора значення у(х)

y(x₀ +h)=y(x₀ )+h(x₀ )h3 +hn yn (x₀ )

x_i =x(0)+Ih

y_i+1 =y_i +∙(K_1i +2K_2i +2K_3i +2K_4i )

K_1i =h∙f(x_i ,y_i )

K_2i =h∙f(x_i +∙y_i +)

K_3i =h∙f(x_i +

K_4i =h∙f(x_i +h∙y_i +K₃ )

Блок - схема головного модуля по Рунге - Кутту:

Реалізація програми за засобом Рунге – Кутта. :

DECLARE SUB KUTT (T!, X!, Y!, A%, B%, C%, E!, H!, N%, T(), X(), Y(), K1X!, K1Y!, K2X!, K2Y!, K3X!, K3Y!, K4X!, K4Y!)

DECLARE SUB GRAF (T!, X!, Y!, A%, B%, C%, E!, H!, N%, T(), X(), Y(), K1X!, K1Y!, K2X!, K2Y!, K3X!, K3Y!, K4X!, K4Y!)

INPUT "C"; C%

E! = C% * 10 ^ (-4)

H! = E! ^ (1 / 4)

CONST A% = 0: CONST B% = 1

DIM SHARED T!(2000), X!(2000), Y!(2000), K1X!(2000), K1Y!(2000), K2X!(2000), K2Y!(2000), K3X!(2000), K3Y!(2000), K4X!(2000), K4Y!(2000)

T(0) = 0: X(0) = 0: Y(0) = 0

M1: CALL KUTT(T!, X!, Y!, A%, B%, C%, E!, H!, N%, T(), X(), Y(), K1X!, K1Y!, K2X!, K2Y!, K3X!, K3Y!, K4X!, K4Y!)

FOR I% = 0 TO N%

X1(I%) = X(I%)

Y1(I%) = Y(I%)

NEXT I%

H! = H! / 2

CALL KUTT(T!, X!, Y!, A%, B%, C%, E!, H!, N%, T(), X(), Y(), K1X!, K1Y!, K2X!, K2Y!, K3X!, K3Y!, K4X!, K4Y!)

FOR I% = 0 TO N%

IF ABS(X1(I%) - X(I%)) * (16 / 15) > E! THEN

GOTO M1

ELSE GOTO M2

END IF

IF ABS(Y1(I%) - Y(I%)) * (16 / 15) > E! THEN

GOTO M1

ELSE GOTO M2:

END IF

NEXT I%

M2: FOR I% = 1 TO N%

PRINT T(I%), X(I%), Y(I%)

NEXT I%

PRINT "H"; H!

INPUT K!

CALL GRAF(T!, X!, Y!, A%, B%, C%, E!, H!, N%, T(), X(), Y(), K1X!, K1Y!, K2X!, K2Y!, K3X!, K3Y!, K4X!, K4Y!)

END

SUB GRAF (T!, X!, Y!, A%, B%, C%, E!, H!, N%, T(), X(), X1(), K1X!, K1Y!, K2X!, K2Y!, K3X!, K3Y!, K4X!, K4Y!)

SCREEN 2

VIEW (170, 50)-(470, 150)

WINDOW (-1, 1.5)-(1, -1.5)

FOR I% = 0 TO N% - 1

PSET (T(I%), X(I%))

PSET (T(I%), Y(I%))

LINE (T(I%), X(I%))-(T(I% + 1), X(I% + 1))

LINE (T(I%), Y(I%))-(T(I% + 1), Y(I% + 1))

NEXT I%

LINE (-1, 0)-(1, 0)

LINE (0, -1.5)-(0, 1.5)

END SUB

SUB KUTT (T!, X!, Y!, A%, B%, C%, E!, H!, N%, T(), X(), Y(), K1X!, K1Y!, K2X!, K2Y!, K3X!, K3Y!, K4X!, K4Y!)

N% = (B% - A%) / H!

FOR I% = 0 TO N%

T(I%) = T(0) + I% * H!

K1X(I%) = H! * (-2 * X(I%) + 5 * Y(I%))

K1Y(I%) = H! * ((EXP(.5 * Y(I%) + T(I%)) - EXP(-.5 * Y(I%) + T(I%))) / 3 + .5 * Y(I%))

K2X(I%) = H! * (-2 * (X(I%) + K1X(I%) / 2) + 5 * (Y(I%) + K1Y(I%) / 2))

K2Y(I%) = H! * ((EXP(.5 * (Y(I%) + K1Y(I%) / 2) + (T(I%) + H! / 2) - EXP(-.5 * Y(I%) + K1Y(I%) / 2) - (T(I%) + H! / 2))) / 3 + .5 * (Y(I%) + K1Y(I%) / 2))

K3X(I%) = H! * (-2 * (X(I%) + K2X(I%) / 2) + 5 * (Y(I%) + K2Y(I%) / 2))

K3Y(I%) = H! * ((EXP(.5 * (Y(I%) + K2Y(I%) / 2) + (T(I%) + H! / 2) - EXP(-.5 * Y(I%) + K2Y(I%) / 2) - (T(I%) + H! / 2))) / 3 + .5 * (Y(I%) + K2Y(I%) / 2))

K4X(I%) = H! * (-2 * (X(I%) + K3X(I%)) + 5 * (Y(I%) + K3Y(I%)))

K4Y(I%) = H! * ((EXP(.5 * (Y(I%) + K3Y(I%)) + (T(I%) + H!) - EXP(-.5 * Y(I%) + K3Y(I%)) - (T(I%) + H!))) / 3 +

Внимание, отключите Adblock

Вы посетили наш сайт со включенным блокировщиком рекламы!
Ссылка для скачивания станет доступной сразу после отключения Adblock!

Скачать