Матэматычная мадэль сінхроннага генератара (Лабараторная работа № 10), страница 11

WRITE(4,103)T,PSId,PSIq,PSIf,PSIud,PSIuq,PSIa,PSIb,PSIc            285

WRITE(4,104)Id,Iq,If,Iud,Iuq,Ia,Ib,Ic                              286

WRITE(4,105)S,GAM,Mmex,Mem                                         287

103 FORMAT(10F8.4)                                                     288

104 FORMAT(8X,8F9.4)                                                   289

105 FORMAT(8X,4F8.4/)                                                  290

C        ВЫВОД ИНФОРМАЦИИ НА ЭКРАН ДИСПЛЕЯ ДЛЯ ТЕКУЩЕГО КОНТРОЛЯ

WRITE(5,106)T,NREZ                                                 295

106 FORMAT(1X,'T=',F6.4,'    NREZ=',I3)                                296

C        ЗАДАНИЕ ГРАНИЧНЫХ УСЛОВИЙ ДЛЯ ВНЕШНЕЙ ЦЕПИ ГЕНЕРАТОРА, СООТC        ВЕТСТВУЮЩИХ РЕЖИМУ ТРЕХФАЗНОГО КОРОТКОГО ЗАМЫКАНИЯ, ПО ИСТЕC        ЧЕНИИ ДВУХ ПЕРИОДОВ РАСЧЕТА НОРМАЛЬНОГО РЕЖИМА

5 IF(T.LT.0.04)GO TO 6                                               300

Ud=0.                                                              310

UQ=0.                                                              320

IF(T.GT.0.04+H)GO TO 6                                             330

WRITE(3,*)'   АВАРИЙНЫЙ РЕЖИМ'                                     331

C        ОБРАЩЕНИЕ К ПОДПРОГРАММЕ РЕШЕНИЯ ДИФФЕРЕНЦИАЛЬНЫХ УРАВНЕНИЙ

6 CALL RNGKT4(NUR,TAU,H1,Y,Y0,PRAV)                                  340

C        ПЕРИОДИЗАЦИЯ УГЛА GAM

IF(Y(7).LT.6.2831853)GO TO 9                                       345

Y(7)=Y(7)-6.2831853                                                346

Y0(7)=Y(7)                                                         347

C         ПРОДОЛЖАТЬ ЛИ ВЫЧИСЛЕНИЯ?

9 IF(NREZ.LT.502)GO TO 10                                            350

WRITE(5,*)'SINGEN: РАБОТА ЗАКОНЧЕНА ПО ПРИЧИНЕ НАПОЛНЕНИЯ ВЫХОДНОГ

*О МАССИВА'

GO TO 11                                                           360

10 IF ((T-TKON).GT.0.)GO TO 11                                        370

C        НАРАЩИВАНИЕ НЕЗАВИСИМОЙ ПЕРЕМЕННОЙ НА ВЕЛИЧИНУ ШАГА

TAU=TAU+H1                                                         380

T=TAU/OMEGA                                                        390

GO TO 4                                                            400

C        ЗАКРЫТИЕ ВЫХОДНЫХ  ФАЙЛОВ

11 CLOSE(UNIT=3)                                                      410

CLOSE(UNIT=4)                                                      411

WRITE(5,101) Z2,Z3                                                 420

101 FORMAT(1X,'SINGEN: РАБОТА ЗАКОНЧЕНА. Р-ТЫ ПОМЕЩЕНЫ В ФАЙЛЫ ',A8,   421

*' И ',A8)                                                          422

STOP                                                               430

END                                                                431

C

SUBROUTINE RNGKT4(N,X,H,Y,Y0,FNC)                                  600

C ПОДПРОГРАММА РЕШЕНИЯ СИСТЕМЫ ДИФФЕРЕНЦИАЛЬНЫХ УРАВНЕНИЙ dY/dx=F(x,Y)

C С НАЧАЛЬНЫМИ УСЛОВИЯМИ (ЗАДАЧА КОШИ) МЕТОДОМ РУНГЕ-КУТТА ЧЕТВЕРТОГО

C ПОРЯДКА С ЗАДАННЫМ ШАГОМ ПО ВРЕМЕНИ НА ОДНОМ ВРЕМЕННОМ ШАГЕ

C (

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

C    N - КОЛИЧЕСТВО УРАВНЕНИЙ В СИСТЕМЕ (ВХОДНАЯ ВЕЛИЧИНА)

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

C        ВЕЛИЧИНА)

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

C    Y0(N) - МАССИВ ЗНАЧЕНИЙ ИНТЕГРИРУЕМЫХ ПЕРЕМЕННЫХ Y НА НАЧАЛО КАЖДОГО

C        ШАГА. ПРИ ПЕРВОМ ОБРАЩЕНИИ К ПОДПРОГРАММЕ В ЭТОМ МАССИВЕ ДОЛЖC        НЫ НАХОДИТЬСЯ НАЧАЛЬНЫЕ ЗНАЧЕНИЯ ИНТЕГРИРУЕМЫХ ПЕРЕМЕННЫХ

C        (ВХОДНЫЕ И РАБОЧИЕ ВЕЛИЧИНЫ)

C    Y(N) - МАССИВ ТЕКУЩИХ ЗНАЧЕНИЙ ИНТЕГРИРУЕМЫХ ПЕРЕМЕННЫХ Y ВНУТРИ

C        ШАГА (РАБОЧИЙ МАССИВ)

C    YF(N) - МАССИВ ЗНАЧЕНИЙ ИНТЕГРИРУЕМЫХ ПЕРЕМЕННЫХ Y, ИСПОЛЬЗУЕМЫХ

C        В КАЧЕСТВЕ АРГУМЕНТОВ ПРИ ВЫЧИСЛЕНИИ ФУНКЦИЙ F(x,Y) В ПРАВЫХ

C        ЧАСТЯХ УРАВНЕНИЙ (РАБОЧИЙ МАССИВ)

C    F(N) - МАССИВ ЗНАЧЕНИЙ ФУНКЦИЙ F(x,Y) В ПРАВЫХ ЧАСТЯХ УРАВНЕНИЙ

C        (РАБОЧИЙ МАССИВ).ЗНАЧЕНИЯ ФУНКЦИЙ ВЫЧИСЛЯЮТСЯ ВО ВНЕШНЕЙ

C        ПОДПРОГРАММЕ FNC, КОТОРАЯ ОПРЕДЕЛЯЕТСЯ ВИДОМ ПРАВОЙ ЧАСТИ

C        СИСТЕМЫ УРАВНЕНИЙ И ДЛЯ КАЖДОЙ НОВОЙ СИСТЕМЫ СОСТАВЛЯЕТСЯ

C        ПОЛЬЗОВАТЕЛЕМ ЗАНОВО

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

C        ФУНКЦИЙ F(x,Y) В ПРАВОЙ ЧАСТИ СИСТЕМЫ. ИМЯ ЭТОЙ ПОДПРОГРАММЫ