Моделирование переходного режима энергосистемы (Лабораторная работа № 9), страница 7

6.5. Какие виды вероятности характерны для разных способов решения дифференциальных уравнений.

6.6. Что такое одношаговые и многошаговые числовые методы.Чем определяется порядок числового метода.

6.7. Какой геометрический смысл имеют явный и неявный методы Эйлера.

6.8. Что такое числовое интегрирование функции.

6.9. Как распространяются алгоритмы на системы дифференциальных уравнений.

6.10. Объясните структуру, назначения операторов и реализацию алгоритмов использованных вами в  программах.

6.11. Объясните результаты, получившиеся вами при выполнении работы.

6.12. Покажите, как исполняются отношения между параметрами режима по первому и второму законам Кирхгофа на получившихся вами векторных диаграммах и осциллограммах.

6.13. Соблюдите проверку совпадения результатов моделирования установленного послеаварийного режима с параметрами переходного режима в конце переходного процесса.

Дополнение 1

C         ПРОГРАММА 'DIFRK409' 16.09.97;

C         РЕШЕНИЕ  СиСТеМЫ  ДИФ. УРАВНЕНИЙ

С         dY/dx = F(x,Y)  МЕТОДОМ РУНГЕ-КУТТЫА ЧЕТВЕРТОГО ПАРЯДКА  С  ПОСТОЯННЫМ  ШАГОМ  ИЗМЕНЕНИЯ   НЕЗАВИСИМОЙ 

С         ПЕРЕМЕННОЙ.   АЛГОРИТМ  НА   ОДНОМ  ВРЕМЕННОМ ШАГЕ ДЛЯ  КОМПОНЕНТЫ ВЕКТОРНОЙ  фУНКЦИИ Y ИМЕЕТ  ВИД:

C            yi+1=yi+(h/6)(k1+2k2+2k3+k4),

C                k1=f(xi,yi,

C                k2=f(xi+h/2,yi+(h/2)k1),

C                k3=f(xi+h/2,yi+(h/2)k2),

C                k4=f(xi+h,yi+hk3),

C            i=0,1,2, ... ,m-1 (m-КОЛИЧЕСТВО ШАГОВ ПО ОСИ x)

C            ОБОЗНАЧЕНИЕ ПЕРЕМЕННЫХ В ПРОГРАММЕ:

C  N - КОЛИЧЕСТВО   ДИФФЕРЕНЦИАЛЬНЫХ  УРАВНЕНИЙ В СИСТЕМЕ-

С МЕ (ВХОДНАЯ   ВЕЛИЧИНА).  МАКСИМАЛЬНОЕ КОЛИЧЕСТВО  УРАВНЕНИЙ  СИСТЕМЫ  ПРИНЯТО  В    ПРОГРАММЕ  РАВНОЙ

С         16, НО ЭТО ОГРАНИЧЕНИЕ СВЯЗАНО ТОЛЬКО С ПРИНЯТЫС         МИ  В  ПРОГРАММЕ  РАЗМЕРНОСТЯМИ  МАССИВОВ   Y(N), Y0(N)  I

С         РАЗМЕРАМИ  ТАБЛИЦ   РЕЗУЛЬТАТОВ   AAAA.REZ   I   AAAA.GRF

C         INTWR - ИНТЕРВАЛ ПЕЧАТИ - КОЛИЧЕСТВО  ШАГОВ  ПО  НЕЗАС         ВИСИМОЙ  ПЕРЕМЕННОЙ    x,  КОТОРЫЕ   ПРОПУСКАЮТСЯ  ПРИ

С        ВЫВОДЕ    РЕУЛЬТАТОВ  РЕШЕНИЯ НА ЭКРАН  ИЛИ  ПРИНТЕР 

С         (ВХОДНАЯ ВЕЛИЧИНА).  НАПРИМЕР,  ПРИ  X0=0,  XKAN=

С         0.1, H=0.0001 РЕШЕНИЕ СИСТЕМЫ   УРАВНЕНИЙ    БУДУТ   

С        ВЫЧИСЛЯТСЯ  В  (XKAN-XO)/H=(0.1-0.0)/0.0001=1000     УЗЛАХ

С         ОСИ   ПЕРЕМЕННОЙ.   НО   ДЛЯ   АНАЛИЗА    ПЕРЕХОДНОГО

С         ПРОЦЕССА  ИССЛЕДОВАТЕЛЬ СЧИТАЕТ ДОСТАТОЧНЫМ СДЕЛАТЬ ВЫС         ВОД  РЕШЕНИЯ СИСТЕМЫ  В 100 УЗЛАХ. ТОГДА INTWR=1000/

С         100=10 (БУДЕТ   СДЕЛАН   ВЫВОД   РЕЗУЛЬТАТОВ РАСЧЕТА   В 

С         КАЖДОМ  ДЕСЯТОМ УЗЛЕ НЕЗАВИСИМОЙ ПЕРЕМЕННОЙ x);

C  X - НЕЗАВИСИМАЯ   ПЕРЕМЕННАЯ  x   (РАБОЧАЯ   И    ВЫХОДНАЯ

С         ВЕЛИЧИНА)

C  X0 - НАЧАЛЬНОЕ ЗНАЧЕНИЕ   x(0)   НЕЗАВИСИМОЙ    ПЕРЕМЕН-  

C         НОЙ  x (ВХОДНАЯ ВЕЛИЧИНА);

C  XKAN - ОКОНЧАТЕЛЬНОЕ ЗНАЧЕНИЕ НЕЗАВИСИМОЙ ПЕРЕМЕННОЙ

C         (ВХОДНАЯ ВЕЛИЧИНА);

C  H - ШАГ ИЗМЕНЕНИЯ НЕЗАВИСИМОЙ  ПЕРЕМЕННОЙ x  (ВХОДС         НАЯ ВЕЛИЧИНА);

C  Y0(16) -  МАССИВ НАЧАЛЬНЫХ ЗНАЧЕНИЙ ВЕКТОРА Y(x), КОМС         ПОНЕНТАМИ  КОТОРОГО  ЯВЛЯЮТСЯ   ЗНАЧЕНИЯ ЗАВИСИМЫХ

С         ПЕРЕМЕННЫХ  y1(0),y2(0),y3(0), ...ДЛЯ  х=0  (ВХОДНЫЕ ВЕЛИС         ЧИНЫ);

C  Y(16) - МАССИВ ТЕКУЩИХ ЗНАЧЕНИЙ ВЕКТОРОВ Y(x),  КОТОРЫЕ  СЧИТАЮТС         СЯ  В  КОНЦЕ КАЖДОГО ШАГА МЕТОДА РУНГЕ-УГОЛС        СЯ (РАБОЧИЕ ПЕРЕМЕННЫЕ И ВЫХОДНЫЕ ВЕЛИЧИНЫ);

C  Z1,Z2,Z3 - СИСВОЛЬНЫЕ ПЕРЕМЕННЫЕ ИМЕН   ФАЙЛОВ  С  ВС         ХОДНЫМИ И ВЫХОДНЫМИ ДАННЫМИ ПРОГРАММЫ:

C NREZ - СЧЕТЧИК КОЛ-ВА ЗАПИСЕЙ, КОТОРЫЕ  ВЫВОДЯТСЯ  В

С        ВЫХОДНЫЕ ФАЙЛЫ

C PRAV - ИМЯ ФАЙЛА, В КОТОРЫХ  НАХОДИТСЯ ПОДПРОГРАММА ВЫЧИСС         ЛЕНИЯ ПРАВЫХ ЧАСТЕЙ  F(x,Y) СИСТЕМЫ ДИФ.

-С        УРАВНЕНИЙ. ЭТО ИМЯ ПРИОБЩЕННО В СПИСОК  ФАКТИЧ.С        ПАРАМЕТРОВ ОПЕРАТОРА CALL RNGKT4(...,PRAV) ЭТОЙ

C         ПРОГРАММЫ. ПО ПРИЧИНЕ ТОГО, ЧТО  В  СПИСКЕ  ФОРМАЛЬ-

С         НЫХ ПАРАМЕТРОВ ОПЕРАТОРА SUBROUTINE  RNGKT(...,FNC) 

С         ИДЕНТИФИКАТОРОВ  PRAV СООТВЕТСТВУЕТ  ВТОРОЕ ИМЯ (FNC),  ТО 

С         ИМЯ PRAV ОПИСАНО В ОПЕРАТОРЕ EXTERNAL

REAL L3,I1,I2,I5,I6

DIMENSION Y(16),Y0(16)

C         Оператор  общей  области памяти для  передачи в подпрограмму  и  с

C         подпрограммы  PRAV параметров системы (11)  дифференциальных  урав-